调用MATLAB工具箱函数frl设计线性相位高通FIR滤波器。要求通带截止频奉为0.6 Πrad,限带截止频率为0.45 Πrad,通带最大衰成为0.2dB.阻带最小衰减为45dB.显示所设计的单位脉冲响应h(n)的数据,并画出损耗函数曲线。
代码:
Fs=80000;
T=1/Fs;
wp=2*pi*4000/Fs;
ws=2*pi*20000/Fs;
rp=0.5;rs=45;
[N,wc]=buttord(wp/pi,ws/pi,rp,rs)
[B,A]=butter(N,wc)
clf;
mpplot(B,A,rs);
functionmpplot(B,A,rs)
ifnargin<3ymin=-80;elseymin=-rs-20;end;
[H,W]=freqz(B,A,1000);
m=abs(H);
p=angle(H);
subplot(2,1,1);
plot(W/pi,20*log10(m/max(m)));gridon;
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,ymin,5]);
title('低通滤波器幅频特性曲线')
subplot(2,1,2);
plot(W/pi,p/pi);
xlabel('\omega/\pi');ylabel('y_w(t)/\pi');gridon;
title('滤除噪声后的信号波形')
截图: