samp7_12.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 14 行

%Samp7_12
clf;n=30;                      %滤波器的阶数
f=[0 0.4 0.5 1];
a=[1 1 0 0];
weit=[1 10];
b=firls(n,f,a,weit);          %采用 firls 设计滤波器
[h,w1]=freqz(b);                %计算滤波器的频率响应
bb=remez(n,f,a,weit);                  %采用 remez 设计滤波器
[hh,w2]=freqz(bb);                   % 计算滤波器的频率响应
figure(1)
plot(w1/pi,abs(h),'b.',w2/pi, abs(hh),'r:',f,a,'m');%绘滤波器幅频响应
xlabel('归一化频率');ylabel('振幅');
legend('firls','remez','理想特性');               %给出图例