Quiero solicitar bucle. A continuación se muestra mi código. He aplicado el bucle for que está mal. Aquí Longitud (t) son números de muestra. Aquí quiero crear un bucle 'for' para una muestra número 5 en particular. Entonces, ¿cómo cambiar h en el bucle for. Además, 'phi_sync' debe ser el número de muestra de 5 y 'phi_acc' debe ser la suma de 0 a 5 números de muestra.
fs = 200;
ts = 1/fs;
t = 0:ts:1-ts;
fc = 50;
S = cos(2*pi*fc*t); %%% equation 7.1 noisless
sigma = 1.8;
phi_sync = sigma * randn(1,length(t)); %%%% random variable with gaussian distributed
phi_acc = sigma * randn(1,length(t)); %%%% random variable with gaussian distributed
% h = 0;
for i=1:length(S)
h = phi_sync + phi_acc(i);
end
gracias por tus esfuerzos. Este es en realidad un algoritmo para agregar ruido de fase en una señal. después de eso quiero agregar h en la señal S, para obtener el ruido de fase. Sin embargo, cuando haga la nueva señal con S1 e intente analizar la señal en el dominio de la frecuencia. Su matriz muestra que muestra.
S1 = cos(2*pi*fc*t + h);
L = length(S1);
nfft = L*100; %%%%%% Zero padding
res = fft(S1,nfft)/nfft; % resize into nfft nr of element % normalize the amplitude
f = fs/2*linspace(0,1,nfft/2+1);
res = res(1:nfft/2+1);
figure, plot(f,abs(res));
% figure, plot(f,angle(res));
return