最新MATLAB-IIR数字低通滤波器的设计方法.docx

Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-dateMATLAB-IIR数字低通滤波器的设计方法MATLAB-IIR数字低通滤波器的设计方法MATLAB IIR数字低通滤波器的设计方法1、IIR数字低通滤波器滤波器分两大类：经典滤波器和现代滤波器。经典滤波的概念，是根据傅里叶分析和变换提出的一个工程概念。根据高等数学理论，任何一个满足一定条件的信号，都可以被看成是由无限个正弦波叠加而成。换句话说，就是工程信号是不同频率的正弦波线性叠加而成的，组成信号的不同频率的正弦波叫做信号的频率成分或叫做谐波成分。现代滤波器理论研究的主要内容是从含有噪声的数据记录（又称时间序列）中估计出信号的某些特征或信号本身。经典滤波器分为低通、高通、带通、带阻滤波器。每一种又有模拟滤波器（AF）和数字滤波器（DF）。对数字滤波器，又有IIR滤波器和FIR滤波器。我们要实现的是IIR 数字滤波器。IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。单位响应：IIR滤波器的单位脉冲响应为无限长，网络中有反馈回路。幅频特性：IIR数字滤波器幅频特性精度很高，不是线性相位的，可以应用于对相位信息不敏感的音频信号上IIR DF 的转移函数是：IIR滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。2、滤波器的技术要求低通滤波器：通带截止频率（又称通带上限频率）：阻带下限截止频率：通带允许的最大衰减 ：阻带允许的最小衰减（，的单位dB）：通带上限角频率：阻带下限角频率（，）即 3、IIR 数字滤波器的设计方法IIR数字滤波器的设计方法有两类：间接设计法和直接设计法。间接设计法是借助模拟滤波器设计方法进行设计的，先根据数字滤波器设计指标设计相应的过渡模拟滤波器，再将过渡模拟滤波器转换为数字滤波器。直接设计法师在时域或频域直接设计数字滤波器。IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。设计步骤如下：1）按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。2）根据转换后的技术指标设计模拟低通滤波器；3）再按一定的规则将转换成。4）若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标，然后按上述步骤2）设计出低通，再将转换为所需的。4几种不同类型的滤波器的介绍：因为我们设计的滤波器的冲击响应一般都为实数，所以有这样，如果我们能由，求出，那么就容易得到所需要的。不同类型的的表达式，代表了几种不同类型的滤波器。（1）巴特沃思(Butterworth)滤波器：C为待定常数，N为待定的滤波器阶次。（2）切比雪夫I型(Chebyshev I )滤波器：5巴特沃思模拟低通滤波器的设计由于每一个滤波器的频率范围将直接取决于设计者的应用目的，因此必然是千差万别。为了使设计规范化，我们需要将滤波器的频率参数作归一化处理。设所给的实际频率为（或），归一化后的频率为，对低通模拟滤波器，令显然，。又令归一化复数变量为，显然所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。(1)将实际频率规一化(2)求C和N这样C和N可求。若令3dB，则C1，这样巴特沃思滤波器的设计就只剩一个参数N，这时(3)确定因为，根据上面公式有由解得，1，2，···，2N这样可得求得后，用代替变量，即得实际需要得。6用双线性Z变换法设计IIR数字低通滤波器s平面到z平面的映射关系称为双线性Z变换，由此关系求出及即因为设计滤波器时系数会被约掉，所以又有7、MATLAB设计低通滤波器相关函数 1、Buttord.m本文件用来确定数字低通或模拟低通滤波器的阶次，其调用格式分别是1）N，Wn=buttord(Wp，Ws，Rp，Rs)2）N，Wn=buttord(Wp，Ws，Rp，Rs，s)格式1)对应数字滤波器，式中Wp,Ws分别是通带和阻带的截止频率，实际上它们是归一化频率，其值在01之间，1对应抽样频率的一半。对低通和高通滤波器，Wp,Ws都是标量，对带通和带阻滤波器，Wp,Ws都是1×2的向量。Rp,Rs分别是通带和阻带的衰减，单位为dB。N是求出的相应低通滤波器的阶次，Wn是求出的3dB频率，它和Wp稍有不同。格式2）对应模拟滤波器，式中各个变量的含义和格式1)相同，但Wp,Ws及Wn的单位为rad/s ，因此，它们实际上式频率 。2、Buttap.m本文件用来设计模拟低通原型滤波器 ，其调用格式是z，p，k=buttap(N)N是欲设计的低通原型滤波器的阶次，z，p，k分别是设计出的 的极点、零点及增益。3、Lp2lp.m4、Lp2hp.m5、Lp2bp.m6、Lp2bs.m从文件名可以看出，上述4个文件的功能分别是将模拟低通原型滤波器 转换为实际的低通、高通、带通及带阻滤波器，其调用格式分别为：（1）B，Alp2lp（b，a，Wo）或B，Alp2hp（b，a，Wo）（2）B，Alp2bp（b，a，Wo）或B，Alp2bs（b，a，Wo）式中b，a分别是模拟低通原型滤波器 的分子、分母多项式的系数向量，B，A分别是转换后的 的分子、分母多项式的系数向量；在（1）中，Wo是低通或高通滤波器的截止频率；在（2）中，Wo是带通或带阻滤波器的中心频率，Bw是其带宽。7、Bilinear.m本文件实现双线性变换，即由模拟滤波器 得到数字滤波器 ，而s和z的关系由 给出。其调用格式是：Bz，Az bilinear（B，A，Fs） 式中B，A分别是 的分子、分母多项式的系数向量，Bz，Az分别是 的分子、分母多项式的系数向量， Fs是抽样频率。8、Freqz.m本文件用来在已知 , 的情况下求出系统的频率响应 ，格式是H，wfreqz（b，a，N，whole，Fs）其中N是频率轴的分点数，建议N为2的整次幂；w 是返回频率轴坐标向量，供绘图用；Fs是抽样频率，若Fs1，频率轴给出归一化频率；whole 指定计算的频率范围是从0Fs，缺省时是从0Fs/29、Zp2tf.mb,a=zp2tf(z,p,k) 零极点型传递函数转换到一般传递函数，它用于在极零点已知时求出 和 的系数。8、设计一个IIR 数字滤波器 技术指标：Hz，Hz，dB， dB，抽样频率Hz。 设计步骤：（1） 角频率 （2） 对角频率做预畸变 （3）模拟低通滤波器的阶次，利用函数N，Wn=buttord(Wp，Ws，Rp，Rs，s) 注意：Wp，Ws应该为（2）中的，。（4）设计模拟低通原型滤波器，其调用格式是 z，p，k=buttap(N)。 N是欲设计的低通原型滤波器的阶次，z，p，k分别是设计出的的极点、零点及增益。（5）求模拟低通原型滤波器的分子分母系数，b,a=zp2tf(z,p,k)。（6）求出的分子、分数系数。B，Alp2lp（b，a，Wo）（7）求出的分子、分母系数，利用bilinear函数。（8）求频率响应，利用Freqz函数。（9）画的图 。程序设计：fp=100;fs=300;ap=3;%通带最大衰减 as=20;%阻带最小衰减 Fs1=1000;%抽样频率 wp=2*pi*fp/Fs1ws=2*pi*fs/Fs1 %数字角频率 Fs=Fs1/Fs1;T=1/Fs;Op=2/T*tan(wp/2)Os=2/T*tan(ws/2)%角频率畸变 N,Wn=buttord(Op,Os,ap,as,'s')z,p,k=buttap(N)%G(p)的零点,极点,增益 b,a=zp2tf(z,p,k);%G(p)的分子、分数系数 B,A=lp2lp(b,a,Op);%H(s)的分子、分数系数 Bz,Az=bilinear(B,A,Fs*T/2)%H(z)的分子、分数系数 H,w=freqz(Bz,Az,256,Fs*Fs1)%H plot(w,abs(H),'r');title('低通滤波器');grid on;结果：wp =0.62 ws =1.8850Op =0.6498 Os =2.7528N =2 Wn =0.8727 z = p = -0.7071 + 0.7071i -0.7071 - 0.7071ik =1 Bz =0.1804 0.3607 0.1804 Az =1.0000 -0.4935 0.2150-