[an error occurred while processing this directive]
|
Fs=200; % sampling frequency
A=500; % input amplitude, mV
N=50;
fi=2*pi*(1:N)/N; % input sinc-wave angle
F=51; % input frequency
sinc=[0 1 0 -1]';
cosc=[1 0 -1 0]';
tt=(1:4)/Fs; % running time
Vin=[];
Acalc=[];
for nph=1:N,
Vin=A*sin(2*pi*F*tt+fi(nph));
convs=Vin*sinc;
convc=Vin*cosc;
Acalc(nph)=(convs.^2+convc.^2).^0.5/2;
end
plot(fi,Acalc,[0 2*pi],[sum(Acalc)/N sum(Acalc)/N]),
title(sprintf('Phase-separation relation, Fin=%g, Hz',F)),
xlabel('Input phase, rad'),ylabel('Computed value')
E-mail: info@telesys.ru