fir.m

来自「利用MATLAB开发软件进行编程」· M 代码 · 共 27 行

Fs1=0.3;Fp1=0.45;Fs2=0.65;Fp2=0.75;
dp=0.10875;ds=0.01;
f=[Fs1 Fp1 Fs2 Fp2];a=[0 1 0];
dev=[ds dp ds];
[M,fo,ao,w]=remezord(f,a,dev);
h=remez(M,fo,ao,w);
w=linspace(0,pi,1000);
mag=freqz(h,[1],w);
subplot(2,2,1)
hd=plot(w/pi,20*log10(abs(mag)));

title('FIR')
grid;

subplot(2,2,2)
plot(w/pi,abs(mag));
title('幅频');
grid

subplot(2,2,3)
plot(w/pi,angle(mag));
title('相频');
grid
subplot(2,2,4)
stem(impz(h,80));
title('冲激响应')
grid