samp7_13.m

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

%Samp7_13
clf;n=21;              %滤波器的阶数为奇数
f=[0 0.1 0.2 1];a=[0 0 1 1];          %理想滤波器的幅频特性
b=firls(n,f,a,'h');              %采用 firls 设计奇对称系数滤波器
[h,w]=freqz(b,1,512);            %计算滤波器的脉冲响应
bb=remez(n, f,a,'h');           %采用 remez 设计奇对称系数滤波器
[hh,w]=freqz(b,1,512);           %计算滤波器的脉冲响应
figure(1)
plot (w/pi,abs(h),'b.',w/pi, abs(hh),'r--',f,a,'m');     %绘滤波器幅频响应
xlabel('归一化频率');ylabel('振幅');
legend('firls','remez','理想特性');                %给出图例
t=0:1/50:3;x=sin(2*pi*5*t);    %时间序列和输入信号  
figure(2)
subplot(2,1,1),plot(t(1:100),x(1:100))   %绘制输入信号的前100个样本
title('输入信号')
y1=filter(b,1,x);   %运用firls设计的滤波器进行滤波
y2=filter(bb,1,x);  %运用remez设计的滤波器进行滤波
subplot(2,1,2),
plot(t([1:100]+20/2),y1([1:100]+20/2),t([1:100]+20/2),y2([1:100]+20/2),'r:')
%绘制与输入信号对应的输出信号，考虑了延迟效应
legend('firls', 'remez')
xlabel('时间/s')
title('输出信号')