多带滤波器matlab设计

gay电的信号处理实验期末试题：

fir2设计任意响应的数字滤波器要求设计一个多带滤波器：其在0到p/8的幅度响应为1，在pi/8 到2p/8幅度响应为1，在2pi8到3p/8幅度响应为09，在30燃到4p/8幅度响应为0.8，在4p/8到5pi/8幅度响应为0.7，在50US到6p/8幅度响应为0.6，在6pi8到7pi8幅度响应为0.5，在7p/8到pi幅度响应为0.1，并且滤波器的阶数为60。画出理想滤波器和设计得到的滤波器的幅度频率响应，在图上加上图标说明、标题（标题内容为“理想滤波器和设计滤波器的幅度频率响应”）、x轴和y轴标注。

%%
%多带滤波器的设计
f=[0  1/8  1/8  2/8  2/8  3/8  3/8  4/8  4/8  5/8  5/8  6/8  6/8  7/8   7/8  1];
m=[1   1   1   1   0.9  0.9  0.8  0.8  0.7  0.7  0.6   0.6   0.5   0.5  0.1  0.1 ];             
%理想滤波器是一个阶梯型滤波器，m中的元素对应f中各频率点的幅度响应，因                                                           
%因为是阶梯型，因此，每个频点上有2个幅度响应值。
b=fir2(60,f,m);                                            
%60为滤波器阶数，f为归一化频率向量，m为对应频率点上的幅度函数fir2                                                          
%用来设计多通带任意响应FIR滤波器，该滤波器的幅频特性由向量对f和m确定，f为归一化                                                                                                                   
[h,w]=freqz(b);                                             
figure(2)
plot(f,m,w/pi,abs(h))                                       
%w/pi,频率归一化
grid on;
legend('理想滤波器','设计滤波器');
title('理想滤波器和设计滤波器的幅度频率响应');