Mam za zadanie zaprojektować filtr w dziedzinie czasu i częstotliwości w matlabie czegoś takiego:
t=0:0.001:1;
x=sawtooth(2*pi*50*t);
Filtr w dziedzinie czasu zrobiłem za pomocą polecenia filter, zostało mi w dziedzinie częstotliwości. Domyślam się, że do tego należy użyć funkcji freqz.
Zrobiłem coś takiego:
[N,Wn]=cheb2ord(.34,.27,0.01,100);
[b,a]=cheby2(N,20,Wn,'high');
n=1001;
[h,w]=freqz(b,a,n);
yfft=h.*(x');
yout=ifft(yfft);
plot(real(yout));
i niestety nic nie wychodzi. Chciałem tym powyżej uzyskać wykres funkcji y po przejściu przez filtr w dziedzinie czasu.
Z góry dzięki za wszelkie podpowiedzi.
t=0:0.001:1;
x=sawtooth(2*pi*50*t);
Filtr w dziedzinie czasu zrobiłem za pomocą polecenia filter, zostało mi w dziedzinie częstotliwości. Domyślam się, że do tego należy użyć funkcji freqz.
Zrobiłem coś takiego:
[N,Wn]=cheb2ord(.34,.27,0.01,100);
[b,a]=cheby2(N,20,Wn,'high');
n=1001;
[h,w]=freqz(b,a,n);
yfft=h.*(x');
yout=ifft(yfft);
plot(real(yout));
i niestety nic nie wychodzi. Chciałem tym powyżej uzyskać wykres funkcji y po przejściu przez filtr w dziedzinie czasu.
Z góry dzięki za wszelkie podpowiedzi.
