duopindai.m

程序设计了一个多频带的FIR滤波器

%Samp7_9
clf
f=0:0.1:1;                         %归一化频率点数
m=[0 0 1 1 0 0 1 1 1 0 0];          %幅频特性值
Order=10;                            % 滤波器的阶数
b=fir2(Order,f,m,hamming(Order+1));            %设计滤波器
[h,w]=freqz(b,1,128);                  %计算滤波器的频率响应
subplot(2,1,1)
plot(f,m,w/pi,abs(h),'r:')       %绘制理想幅频响应和设计的滤波器幅频响应
legend('理想特性', '实际设计')   %给出图例
title('Order=10');xlabel('归一化频率');ylabel('振幅');
Order=100;
b=fir2(Order,f,m,hamming(Order+1));               %设计阶数为 100 的滤波器
[h,w]=freqz(b,1,128);                 %计算滤波器的频率响应
subplot(2,1,2),plot(f,m,w/pi,abs(h),'r:'); %绘制理想幅频响应和设计的幅频响应
ylim([0 1])
legend('理想特性', '实际设计')   %给出图例
title('Order=100');xlabel('归一化频率');ylabel('振幅');
f1=12;f2=36;                % 输入信号的两种频率成分
t=0:1/100:2;                   % 时间序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);     %输入信号
y=fftfilt(b,x);                  %对输入信号进行滤波
figure(2)
subplot(2,1,1), plot(t,x),title('输入信号')           %绘制输入信号
subplot(2,1,2),plot(t,y)                        %绘制输出信号
hold on;plot(Order/2/100*ones(1,2),ylim, 'r') %绘制延迟到的时刻
title('输出信号'),xlabel('时间/s')