[an error occurred while processing this directive]
|
Fs=200; % sampling frequency
A=1:500; % input amplitude, mV
fi=0.1; % input sinc-wave angle - наплевать (легко проверить)
Fin=49:51; % input frequency
sinc=[0 1 0 -1]';
cosc=[1 0 -1 0]';
tt=(1:4)/Fs; % running time
for F=Fin,
for Ain=A,
Vin(Ain,:)=Ain*sin(2*pi*F*tt+fi);
end
convs=Vin*sinc;
convc=Vin*cosc;
Acalc=(convs.^2+convc.^2).^0.5/2;
plot(A,Acalc),title(sprintf('amlitude-separation relation, Fin=%g, Hz',F)),
xlabel('Input amplitude, mV'),ylabel('Computed value'),axis([0 500 0 500])
pause(2)
end
E-mail: info@telesys.ru