第6章 无限脉冲响应数字滤波器的设计-.pptVIP

③设计巴特沃斯低通滤波器。先计算阶数N及3dB截止频率Ωc。第94页，共138页，星期日，2025年，2月5日取N=6。为求3dB截止频率Ωc,将Ωp和αp代入(6.2.17)式，得到Ωc=0.7032rad/s，显然此值满足通带技术要求，同时给阻带衰减留一定余量，这对防止频率混叠有一定好处。根据阶数N=6，查表6.2.1，得到归一化传输函数为为去归一化，将p=s/Ωc代入Ha(p)中，得到实际的传输函数Ha(s),第95页，共138页，星期日，2025年，2月5日④用脉冲响应不变法将Ha(s)转换成H(z)。首先将Ha(s)进行部分分式，并按照(6.3.11)式、(6.3.12)式，或者(6.3.13)式和(6.3.14)式，得到：第96页，共138页，星期日，2025年，2月5日图6.4.7例6.4.2图——用脉冲响应不变法设计的数字低通滤波器的幅度特性第97页，共138页，星期日，2025年，2月5日(4)带阻滤波器的H(s)为第62页，共138页，星期日，2025年，2月5日6.3用脉冲响应不变法设计IIR

数字低通滤波器为了保证转换后的H(z)稳定且满足技术要求，对转换关系提出两点要求：(1)因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。(2)数字滤波器的频率响应模仿模拟滤波器的频响，s平面的虚轴映射z平面的单位圆，相应的频率之间成线性关系。第63页，共138页，星期日，2025年，2月5日设模拟滤波器的传输函数为Ha(s),相应的单位冲激响应是ha(t)设模拟滤波器Ha(s)只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将Ha(s)用部分分式表示：(6.3.1)式中si为Ha(s)的单阶极点。将Ha(s)进行逆拉氏变换得到ha(t)：(6.3.2)第64页，共138页，星期日，2025年，2月5日式中u(t)是单位阶跃函数。对ha(t)进行等间隔采样，采样间隔为T，得到：(6.3.3)对上式进行Z变换，得到数字滤波器的系统函数H(z)：(6.3.4)设ha(t)的采样信号用ha(t)表示，第65页，共138页，星期日，2025年，2月5日对进行拉氏变换，得到:式中ha(nT)是ha(t)在采样点t=nT时的幅度值，它与序列h(n)的幅度值相等，即h(n)=ha(nT)，因此得到:(6.3.5)第66页，共138页，星期日，2025年，2月5日上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式表示：我们知道模拟信号ha(t)的傅里叶变换Ha(jΩ)和其采样信号的傅里叶变换之间的关系满足(1.5.5)式，重写如下：(6.3.6)第67页，共138页，星期日，2025年，2月5日将s=jΩ代入上式，得由(6.3.5)式和(6.3.8)式得到:(6.3.7)(6.3.8)(6.3.9)第68页，共138页，星期日，2025年，2月5日上式表明将模拟信号ha(t)的拉氏变换在s平面上沿虚轴按照周期Ωs=2π/T延拓后，再按照(6.3.6)式映射关系，映射到z平面上，就得到H(z)。(6.3.6)式可称为标准映射关系。下面进一步分析这种映射关系。设按照(6.3.6)式，得到:因此得到:(6.3.10)第69页，共138页，星期日，2025年，2月5日那么σ=0,r=1σ0,r1σ0,r1另外，注意到z=esT是一个周期函数，可写成为任意整数第70页，共138页，星期日，2025年，2月5日图6.3.1z=esT,s平面与z平面之间的映射关系第71页，共138页，星期日，2025年，2月5日图6.3.2脉冲响应不变法的频率混叠现象第72页，共138页，星期日，2025年，2月5日假设