Требуется программист в Зеленограде
- обработка данных с датчиков; ColdFire; 40 тыс. e-mail:jobsmp@pochta.ru |
Fs = 2000;
Ts = 1/Fs;
N = 8192;
t = (0:1:(N-1))*Ts;
xi = sin(2*pi*234*t + 1.23);
%xq = -cos(2*pi*234*t + 1.23);
b1 = firpm(50, [0 800 900 Fs/2]/(Fs/2), [1 1 0 0]);
a1 = 1;
%[H1,F] = freqz(b1,a1,4096,'whole', Fs);
n = (0:1:(length(b1)-1));
b2 = b1 .* exp(1j*2*pi*300*n*Ts);
[H2real,F] = freqz(real(b2),1,4096,'whole', Fs);
[H2imag,F] = freqz(imag(b2),1,4096,'whole', Fs);
subplot(3,1,1)
plot(F, abs(H2real))
subplot(3,1,2)
plot(F, abs(H2imag))
subplot(3,1,3)
plot(F, unwrap(angle(H2real) - angle(H2imag))*180/pi);