Matlab设计IIR基本示例返回-Read.PDFVIP

中国科学技术大学电子工程与信息科学系 多媒体通信实验室 （Copyright 2003 ） 滤波器设计示例 1、Matlab 设计IIR 基本示例 1 （1）直接设计数字滤波器 1 （2 ）脉冲相应不变法2 （3 ）双线性变换2 （4 ）双线形变换法设计chebyshev 高通滤波器3 （5 ）使用MatLab6 下的Filter Designed Tool 4 2 、Matlab 设计IIR 高级示例4 （1）双线性变换实现ButterWorth 低通4 （2 ）双线性变换法实现Chebyshev 低通 （I 型）5 （3 ）低通变换为高通6 （2 ）低通变换为带通7 （3 ）低通变换为带通（双线性变换）8 3、Matlab 窗函数设计FIR 滤波器基本示例9 （1）低通滤波器实现9 （2 ）带通滤波器实现9 1、Matlab 设计IIR 基本示例 [返回] （1）直接设计数字滤波器 [返回] 低通，采样频率为 1Hz ，通带临界频率fp =0.2Hz ，通带内衰减小于1dB （αp=1 ）；阻 带临界频率fs=0.3Hz ，阻带内衰减大于25dB （αs=25 ）。设计一个数字滤波器满足以上参数。 % Matlab 使用归一化的频率参数（临界频率除以采样频率的1/2）。 %这样临界频率参数的取值范围在0 和1 之间，1 代表Fs/2 （用角频率表示的时候对应π） FS=1 [n,Wn]=buttord(0.2/(FS/2),0.3/( FS /2),1,25); [b,a]=butter(n,Wn); freqz(b,a,512, FS); 如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ustc.edu.cn cxh@ustc.edu.cn 中国科学技术大学电子工程与信息科学系 多媒体通信实验室 （Copyright 2003 ） （2 ）脉冲相应不变法 [返回] 低通，采样频率为 1Hz ，通带临界频率fp =0.2Hz ，通带内衰减小于1dB （αp=1 ）；阻 带临界频率fs=0.3Hz ，阻带内衰减大于25dB （αs=25 ）。设计一个数字滤波器满足以上参数。 FS=1 [n,Wn]=buttord(0.2*2*pi,0.3*2*pi,1,25,s); %临界频率采用角频率表示 [b,a]=butter(n,Wn,s); %freqs(b,a) %设计模拟的 [bz,az]=impinvar(b,a,FS); %映射为数字的 freqz(bz,az,512,FS) （3 ）双线性变换 [返回] 低通，采样频率为 1Hz ，通带临界频率fp =0.2Hz ，通带内衰减小于1dB （αp=1 ）；阻 带临界频率fs=0.3Hz ，阻带内衰减大于25dB （αs=25 ）。设计一个数字滤波器满足以上参数。 FS=1 %通带、阻带截止频率 Fl=0.2;Fh=0.3; %频率预畸 wp=(Fl/FS)*2*pi; %临界频率采用角频率表示 ws=(Fh/FS)*2*pi; %临界频率采用角频率表示 OmegaP=2*FS*tan(wp/2); OmegaS=2*FS*tan(ws/2); [n,Wn]=buttord(OmegaP,OmegaS,1,25,s); [b,a]=butter(n,Wn,s); %freqs(b,a) %设计模拟的 [bz,az]=bilinear(b,a,FS); %映射为数字的 freqz(bz,az,512,FS,whole) 如果您在阅读过程中发现疏漏和错误，请您尽快和编者取得联系 network@ustc.edu.cn c