zremez.m

来自「用MATLAB实现数字信号处理的程序」· M 代码 · 共 13 行

%Remz函数：通过最佳一致逼近法设计FIR滤波器
n=17;
f=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 1];
a=[0 0 0 0 1 1 1 0 0 0 ];
w=[0.9 0.85 0.8 0.2 0.1];
b1=remez(n,f,a);
b2=remez(n,f,a,w);
[H1,W1]=freqz(b1,1,512);
[H2,W2]=freqz(b2,1,512);
plot(f,a,'g-',W1/pi,abs(H1),'b*',W2/pi,abs(H2),'r:');
legend('期望曲线','没有加权时的设计结果','有加权时的设计结果');
grid;