samp7_11.m

MATLAB7.x数字信号处理 光盘内容

%Samp7_11
clf;n=50;                               % 滤波器的阶数
f=[0 0.1 0.15 0.25 0.3 0.4 0.45 0.55 0.6 0.7 0.75 0.85 0.9 1];% 频率向量
a=[1 1 0 0 1 1 0 0 1 1 0 0 1 1];           %振幅向量
b=firls(n, f,a);                    % 采用 firls 设计滤波器
[h,w1]=freqz(b);                         % 计算滤波器的频率响应
bb=remez(n, f,a);                     %采用 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','理想特性');               %给出图例
figure(2)
impz(bb,1),title ('脉冲响应')            %给出滤波器的脉冲响应
xlabel('样本数');ylabel('幅度')