samp6_7.m

本例程是用MATLAB来实现数字滤波器的设计。

%Samp6_7
wp=0.2*pi;ws=0.3*pi;Rp=1;Rs=15; %数字滤波器截止频率通带波纹和阻带衰减
Fs=100;Ts=1/Fs;Nn=128;  %采样频率
Wp=2/Ts*tan(wp/2.);Ws=2/Ts*tan(ws/2.);   %按频率转换公式进行转换
[N,Wn]=ellipord(Wp,Ws,Rp,Rs,'s');  %计算模拟滤波器的最小阶数
[z,p,k]=ellipap(N,Rp,Rs);  %设计模拟原型滤波器
[Bap,Aap]=zp2tf(z,p,k);  %零点极点增益形式转换为传递函数形式
[b,a]=lp2lp(Bap,Aap,Wn); %低通转换为低通滤波器的频率转换
[bz,az]=bilinear(b,a,Fs);  %运用双线性变换法得到数字滤波器传递函数
[H,f]=freqz(bz,az,Nn,Fs);   %求出频率特性
subplot(2,1,1),plot(f,20*log10(abs(H)));
xlabel('频率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('频率/Hz');ylabel('相位/^o');grid on;