IIR数字滤波器的设计
Mrakdown常见数学符号
Markdown基础语法
冲激响应不变法
-
冲激响应不变法:就是用其单位冲激响应序列模仿模拟滤波器的单位冲激响应的抽样值
-
设计的具体步骤及方法
首先要设计一个响应的模拟滤波器。
模拟滤波器的单位冲激响应非周期 ⇒ 频率离散(- ∞-+∞ ),频谱范围可能大于折叠频率,即奈奎斯特频率 ⇒ 取样,频率周期延拓 ⇒ 频谱可能产生混叠。
数字域中极点在单位圆内⇔模拟域中极点在左半平面。
可以通过提高抽样频率来减少混叠,但设计指标若以数字域频率给定时,不能通过提高抽样频率改善混叠现象。因为Ω=Tω ,ω是数字域频率,Ω是模拟域频率,增大抽样频率⇒减小采样周期⇒增大截止频率。
模拟域 :Ha(s) = ∑k=1Ns−skAk 数字域:H(z) = ∑k=1N1−eskTz−1Ak
s平面:s=sk z平面:z=eskT 两者系数相同
当采样周期很小时,数字滤波器增益很大,容易溢出,可以修正,令h(n)=Tha(nT),因为H(ejw)=T1Ha(jTω)
该方法对应的MATLAB程序:[BZ, AZ] = impinvar(B,A,Fs),其中BZ,AZ为数字滤波器单位冲击响应的分子和分母系数;B,A为模拟滤波器的分子和分母的系数,Fs为抽样频率。
-
优点
- 时域逼近良好
- 保持线性关系:ω=ΩT
-
缺点
阶跃响应不变法
- 变换原理:数字滤波器的阶跃响应模仿模拟滤波器的阶跃响应,即g(n)=ga(t)|t=nt=ga(nT).
- g(n)=u(n)*h(n)→ G(z)=z−1zH(z) ,ga(t)=u(t)*ha(t)→Ga(s)=s1Ha(s)
变换步骤:
Ha(s)→Ga(s)→ga(t)→g(n)→G(z)→H(z)
- 优缺点同冲激响应不变法
双线性不变法
在推导过程中有一个中间域,最终的关系式为z=1−s1+s ,为了使模拟滤波器与数字滤波器的某一频率有对应关系,引入一个常数c,则有s=c1+z11−z1 ,即z=c−sc+s。若低频处有确切的对应关系应为c=T2 。
s平面虚轴与z平面单位圆的对应关系Ω =ctan2w Ω 为s域,w为数字域
-
优点:无频率混叠,避免了多值映射。
-
缺点:除零频率附近,严重非线性关系,在临界频率点(通带截止频率和阻带截止频率–这个是技术指标)产生畸变。
-
解决临界频率点的畸变–预畸变
Ω =ctan2w ,在临界频率点根据w去设计Ω ,再去设计模拟滤波器,再用双线性变换法设计数字滤波器
-
该方法对应的MATLAB方法:[BZ, AZ] = bilinear(B,A,Fs)
常用模拟滤波器设计
-
由幅度平方函数确定模拟滤波器的系统函数
∣Ha(jΩ)∣2=Ha(jΩ)Ha∗(jΩ)=Ha(jΩ)Ha(−jΩ)=Ha(s)Ha(−s)∣(s=jΩ) ,Ha(jΩ)可以看做模拟滤波器单位冲激响应的傅里叶变化。
Ha(s)和Ha(−s)的零极点对称,成对出现,且其零极点分别分布在左半平面(只有分布在左半平面系统才是稳定的)和右半平面。
-
Butterworth低通逼近
幅度平方函数为∣Ha(jΩ)∣2=1+(ΩcΩ)2N1 ,其中N为滤波器阶数,Ωc为3dB截止频率,因为Ω=Ωc时,幅度为原来的21,
![IIR数字滤波器的设计 IIR数字滤波器的设计](/default/index/img?u=aHR0cHM6Ly9waWFuc2hlbi5jb20vaW1hZ2VzLzIwMC8xZDZjYTU0MjQ0ZjBlYmU1N2FmMzZlZTRiZTdiMDBlMC5wbmc=)
巴特沃斯幅度特性及其与N的关系
-
∣Ha(jΩ)∣Ω=js2=Ha(s)Ha(−s)=1+(jΩs)2N1 ,令1+(jΩs)2N=0即可解得极点分布。
![IIR数字滤波器的设计 IIR数字滤波器的设计](/default/index/img?u=aHR0cHM6Ly9waWFuc2hlbi5jb20vaW1hZ2VzLzM0NS8xODYzNDIwZWI1MDRkZTZmOWRiOGY0ZDBlZjE3ZTY4MS5wbmc=)
三阶和四阶的极点分布图
极点在s平面呈象限对称,个数为2N,极点间角度间隔为Nπ,极点不落在虚轴上,若N为奇数,实轴上有极点;若为偶数,实轴上无极点。
Ha(s)=∏k=1N(s−sk)ΩcN,sk为左半平面的极点。若令Ωc=1rad/s,得到的各个阶次的系统函数为归一化的系统函数Han(s),其与原系统函数之间的关系为Ha(s)=Han(s)∣s=Ωcs。