实验四IIR数字滤波器的设计-代码.docVIP

实验四 用冲激响应不变法将模拟滤波器变换成数字滤波器H(z)， 采样周期。 b=[0,3,2]; a=[2,3,1]; ts=0.1; [bz1,az1]=impinvar(b,a,1/ts); freqz(bz1,az1,128,4000); ，分别用冲激响应不变法及双线性变换法设计一个三阶巴特沃思数字低通滤波器，给出所设计的滤波器的频率特性曲线。 [B,A]=butter(3,2*pi*1000,s); %设计一个的3阶模拟滤波器 [num1,den1]=impinvar(B,A,4000); %冲击响应不变法，转换为数字滤波器 [num2,den2]=bilinear(B,A,4000); %双线性变换，转换为数字滤波器 [h1,f1]=freqz(num1,den1,128,4000); %计算频率响应 [h2,f2]=freqz(num2,den2,128,4000); plot(f1, abs(h1),-.,f2, abs(h2),-); 用冲激响应不变法及双线性变换法设计一个巴特沃思数字低通滤波器，给定指标：，，通带纹波小于3dB, 阻带幅度衰减大于15dB，，采样频率，给出所设计的滤波器的频率特性曲线。给出所设计的滤波器的频率特性曲线。 方法一： 冲激响应不变法 omegap=2000*2*pi;omegast=3000*2*pi; rp=3;rs=15; fs=10000; ts=1/fs;Nn=128; [N,Wn]=buttord(omegap,omegast,rp,rs,s); [B,A]=butter(N,Wn,s); [bz,az]=impinvar(B,A,fs); freqz(bz,az,Nn,fs); % Digital filter frequency response. 双线性变换法 omegap=2000*2*pi; omegast=3000*2*pi; rp=3;rs=15; fs=10000; ts=1/fs;Nn=128; wp= omegap*ts;ws= omegast*ts; %求数字频率 wp1=2*fs*tan(wp/2); ws1=2*fs*tan(ws/2); %利用公式进行预畸变 [N,Wn]=buttord(wp1,ws1,rp,rs,s); [z,p,k]=buttap(N); %设计模拟滤波器 [Bap,Aap]=zp2tf(z,p,k); %将零极点转换为传递函数 [b,a]=lp2lp(Bap,Aap,Wn); %从模拟低通到实际低通 [bz,az]=bilinear(b,a,fs); %双线性不变法 freqz(bz,az,Nn,fs); %计算频率响应 方法二： omegap=2000*2*pi;omegast=3000*2*pi; rp=3;rs=15; fs=10000; ts=1/fs;Nn=128; [N,Wn]=buttord(omegap,omegast,rp,rs,s); [z,p,k]=buttap(N); %设计模拟滤波器 [Bap,Aap]=zp2tf(z,p,k); %进行系数转换 [b,a]=lp2lp(Bap,Aap,Wn); %从模拟低通到实际低通 [bz,az]=impinvar(b,a,fs); %冲激响应不变法 freqz(bz,az,Nn,fs); %计算频率响应 . 注意：方法1和方法2的区别： 方法1：通过阶数N，直接得到实际模拟LPF，转换为数字LPF 方法2：通过阶数N，设计LPF原型，参数转换，lptolp得到实际模拟LPF，转换为数字LPF 用冲激响应不变法设计一个巴特沃思数字低通滤波器，给定指标：通带频率，通带纹波小于1dB, 阻带在内，幅度衰减大于15dB, 采样周期，给出所设计的滤波器的频率特性曲线。 方法一： wp=0.4*pi;wst=0.6*pi; rp=1;rs=15; Ts=0.0001;Nn=128; omegap=wp/Ts;omegast=wst/Ts; [N,Wn]=buttord(omegap,omegast,rp,rs,s); [B,A]=butter(N,Wn,s); [bz,az]=impinvar(B,A,1/Ts); freqz(bz,az,Nn,1/Ts); 方法二： wp=0.4*pi;wst=0.6*pi; rp=1;rs=15; Ts=0.0001;Nn=128; omegap=wp/Ts;omegast=wst/Ts; [N,Wn]=buttord(omegap,omegast,rp,rs,s); [z,p,k]=buttap(N);