数字滤波器设计步骤
使用MATLAB进行数字滤波器设计的步骤与方法

使用MATLAB进行数字滤波器设计的步骤与方法数字滤波器是用于信号处理的重要工具,它可以对信号进行去噪、频率调整等操作。
而MATLAB作为一种强大的数学计算软件,提供了丰富的数字信号处理工具箱,可以方便地进行数字滤波器的设计与仿真。
本文将介绍使用MATLAB进行数字滤波器设计的步骤与方法。
1. 了解数字滤波器的基本原理在进行数字滤波器设计之前,首先需要了解数字滤波器的基本原理。
数字滤波器根据其频率响应特性可以分为低通、高通、带通和带阻滤波器等。
此外,数字滤波器的设计还需要考虑滤波器的阶数、截止频率以及滤波器类型等因素。
在设计中,我们可以选择滤波器的类型和相应的参考模型,然后利用MATLAB工具箱提供的函数进行设计。
2. 导入MATLAB中的数字信号处理工具箱使用MATLAB进行数字滤波器设计需要先导入数字信号处理工具箱。
通过在MATLAB命令窗口输入`>> toolbox`即可打开工具箱窗口,并可以选择数字信号处理工具箱进行加载。
加载完成后,就可以调用其中的函数进行数字滤波器设计。
3. 设计数字滤波器在MATLAB中,常用的数字滤波器设计函数有`fir1`、`fir2`、`iirnotch`等。
这些函数可以根据系统特性需求设计相应的数字滤波器。
以FIR滤波器为例,可以使用`fir1`函数进行设计。
该函数需要输入滤波器的阶数和截止频率等参数,输出设计好的滤波器系数。
4. 评估滤波器性能设计好数字滤波器后,需要进行性能评估。
可以使用MATLAB提供的`fvtool`函数绘制滤波器的幅频响应、相频响应和群延迟等。
通过观察滤波器在频域的性能表现,可以判断设计的滤波器是否满足要求。
5. 对滤波器进行仿真在对滤波器性能进行评估之后,还可以使用MATLAB进行滤波器的仿真。
通过将需要滤波的信号输入设计好的滤波器中,观察输出信号的变化,可以验证滤波器的去噪效果和频率调整能力。
MATLAB提供了函数`filter`用于对信号进行滤波处理。
数字滤波器的设计方法

数字滤波器的设计方法数字滤波器是一种用于信号处理的重要工具,可以用于信号去噪、频率选择和信号恢复等应用。
本文将介绍数字滤波器的设计方法,包括滤波器的类型、设计步骤和常用的设计工具。
我们需要了解数字滤波器的类型。
数字滤波器可以分为无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器两种。
IIR滤波器的特点是具有无穷长的冲激响应,而FIR滤波器的冲激响应是有限长的。
接下来,我们来看一下数字滤波器的设计步骤。
首先,我们需要确定滤波器的设计要求,包括滤波器的通带和阻带的频率范围,以及在通带和阻带中的衰减要求。
然后,根据这些设计要求选择合适的滤波器类型,比如IIR滤波器或FIR滤波器。
接下来,我们需要进行滤波器的设计和优化,以满足给定的要求。
最后,我们需要对设计的滤波器进行验证和性能评估。
在数字滤波器的设计过程中,我们可以借助一些常用的设计工具来辅助完成。
其中一种常用的工具是Matlab软件,它提供了丰富的信号处理工具箱,可以方便地进行滤波器的设计、分析和仿真。
另外,还有一些开源的信号处理库,如SciPy和Octave,也可以用于数字滤波器的设计。
除了工具之外,还有一些常用的设计方法可以帮助我们实现数字滤波器的设计。
其中一种方法是基于频率响应的设计方法,即通过设定滤波器在不同频率上的增益来满足设计要求。
这种方法可以通过频域分析和优化来实现。
另一种方法是基于时域响应的设计方法,即通过设定滤波器的冲激响应来满足设计要求。
这种方法可以通过时域分析和优化来实现。
在设计数字滤波器时,还需要考虑滤波器的稳定性和实现的复杂度。
稳定性是指滤波器的输出是否有界,即是否会出现无限增长的情况。
实现的复杂度包括滤波器的计算量和存储量等方面的考虑。
通常情况下,FIR滤波器比IIR滤波器更容易设计和实现,但是在一些特定的应用中,IIR滤波器可能更加适用。
总结起来,数字滤波器的设计是一个复杂而关键的过程,需要根据设计要求选择合适的滤波器类型,进行设计和优化,并进行验证和性能评估。
实验7 窗函数法设计FIR数字滤波器

实验7窗函数法设计FIR数字滤波器一、实验目的掌握窗函数法设计F1R数字滤波器的原理和具体方法二、实验设备与环境计算机、Mat1ab软件环境三、实验基础理论1>基本原理窗函数设计法的基本思想为,首先选择一个适当的理想的滤波器Hd(,3),然后用窗函数截取它的单位脉冲响应%(九),得到线性相位和因果的FIR滤波器,这种方法的重点是选择一个合适的窗函数和理想滤波器,使设计的滤波器的单位脉冲响应逼近理想滤波器的单位脉冲响应。
2、设计步骤(1)给定理想滤波器的频率响应Hd("3),在通带上具有单位增益和线性相位,在阻带上具有零响应。
一个带宽为g(3c<Tr)的低通滤波器由下式给定h(e j^=(eW∣ω∣≤ωc虱)一1Oωc<∣ω∣<π其中α为采样延迟,其作用是为了得到因果的系统。
(2)确定这个滤波器的单位脉冲响应为了得到一个h(n)长度为N的因果的线性相位FIR滤波器,我们令N-Ia=-2-(3)用窗函数截取hd(τι)得到所设计FIR数字滤波器h(n)h(n)=h d(n)w(n)3、窗函数的选择常用的窗函数有矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗等。
Mat1ab提供了一些函数用于产生窗函数,如下表所示:在设计过程中我们需要根据给定的滤波器技术指标,选择滤波器长度N 和窗函数3(n)°表7.2列出了常用的窗函数的一些特性,可供设计时参考。
其中幻是修正的零阶贝塞尔函数,参数B 控制最小阻带衰减,这种窗函数对于相同的N 可以提供不同的过渡带宽。
由于贝塞尔函数比较更杂,这种窗函数的设计方程很难推导,然而幸运的是,有一些经验设计方程可以直接使用。
已知给定的指标叫Msc,Rp 和4,滤波器长度N 和凯瑟窗参数B 可以按如下凯瑟窗设计方程给出过渡带宽:∆ω=ωst -ωp入一7.95 2.285∆ω_(0.1102(4-8.7) ,P=iθ.5842(4-21)04+0.07886(4-21), 四、实验内容1、设计一个数字低通FIR 滤波器,其技术指标如下ωp =0.2τr,RP=0.25dBωst =0.3τr,A s =50dB分别采用矩形窗、汉宁窗、海明窗、布莱克曼窗、凯瑟窗设计该滤波器。
实验五FIR数字滤波器的设计

实验五FIR数字滤波器的设计
FIR数字滤波器的设计可以分为以下几个步骤:
1.确定滤波器的类型和规格:根据实际需求确定滤波器的类型(如低通、高通、带通等)以及滤波器的截止频率、通带衰减以及阻带衰减等规格。
2.选择滤波器的窗函数:根据滤波器的规格,选择合适的窗函数(如矩形窗、汉宁窗、布莱克曼窗等)。
窗函数的选择会影响滤波器的频率响应以及滤波器的过渡带宽度等特性。
3.确定滤波器的阶数:根据滤波器的规格和窗函数的选择,确定滤波器的阶数。
通常来说,滤波器的阶数越高,滤波器的性能越好,但相应的计算和处理也会更加复杂。
4.设计滤波器的频率响应:通过在频率域中设计滤波器的频率响应来满足滤波器的规格要求。
可以使用频率采样法、窗函数法或优化算法等方法。
5. 将频率响应转换为差分方程:通过逆Fourier变换或其他变换方法,将频率响应转换为滤波器的差分方程表示。
6.量化滤波器的系数:将差分方程中的连续系数离散化为滤波器的实际系数。
7.实现滤波器:使用计算机编程、数字信号处理芯片或FPGA等方式实现滤波器的功能。
8.测试滤波器性能:通过输入一组测试信号并观察输出信号,来验证滤波器的性能是否符合设计要求。
需要注意的是,FIR数字滤波器的设计涉及到频率域和时域的转换,以及滤波器系数的选择和调整等过程,需要一定的信号处理和数学背景知识。
数字信号处理实验报告四IIR数字滤波器设计及软件实现

数字信号处理实验报告四IIR数字滤波器设计及软件实现实验目的:本实验的目的是了解IIR数字滤波器的设计原理和实现方法,通过MATLAB软件进行数字滤波器设计和信号处理实验。
一、实验原理IIR数字滤波器是一种使用有限数量的输入样本和前一次输出值的滤波器。
它通常由差分方程和差分方程的系数表示。
IIR滤波器的特点是递归结构,故其频率响应是无限长的,也就是说它的频率响应在整个频率范围内都是存在的,而不像FIR滤波器那样只有在截止频率处才有响应。
根据设计要求选择合适的滤波器类型和滤波器结构,然后通过对滤波器的模型进行参数化,设计出满足滤波要求的IIR滤波器。
常见的IIR滤波器设计方法有模拟滤波器设计方法和数字滤波器设计方法。
在本实验中,我们主要使用数字滤波器设计方法,即离散时间滤波器设计方法。
二、实验内容(一)设计IIR数字滤波器的步骤:1.确定滤波器类型:根据滤波要求选择合适的滤波器类型,如低通滤波器、高通滤波器、带通滤波器、带阻滤波器等。
2.确定滤波器的阶数:根据滤波要求确定滤波器的阶数。
阶数越高,滤波器的频率响应越陡峭,但计算复杂度也越高。
3. 设计滤波器原型:根据滤波要求,设计滤波器的原型。
可以选择Butterworth滤波器、Chebyshev滤波器、Elliptic滤波器等作为原型。
4.选择滤波器结构:根据计算机实现条件和算法复杂度,选择合适的滤波器结构。
常见的滤波器结构有直接形式I、直接形式II、级联形式等。
5.参数化滤波器模型:根据原型滤波器的差分方程,选择合适的参数化方法。
常见的参数化方法有差分方程法、极点/零点法、增益法等。
6.根据参数化的滤波器模型,计算出所有的滤波器系数。
(二)用MATLAB软件实现IIR数字滤波器设计:1.打开MATLAB软件,并创建新的脚本文件。
2. 在脚本文件中,使用MATLAB提供的滤波器设计函数,如butter、cheby1、ellip等,选择合适的滤波器类型进行设计。
实验四FIR数字滤波器的设计

实验四FIR数字滤波器的设计
FIR数字滤波器也称作有限脉冲响应数字滤波器,是一种常见的数字滤波器设计方法。
在设计FIR数字滤波器时,需要确定滤波器的阶数、滤波器的类型(低通、高通、带通、带阻)以及滤波器的参数(截止频率、通带波纹、阻带衰减、过渡带宽等)。
下面是FIR数字滤波器的设计步骤:
1.确定滤波器的阶数。
阶数决定了滤波器的复杂度,一般情况下,阶数越高,滤波器的性能越好,但计算量也越大。
阶数的选择需要根据实际应用来进行权衡。
2.确定滤波器的类型。
根据实际需求,选择低通、高通、带通或带阻滤波器。
低通滤波器用于去除高频噪声,高通滤波器用于去除低频噪声,带通滤波器用于保留一定范围内的频率信号,带阻滤波器用于去除一定范围内的频率信号。
3.确定滤波器的参数。
根据实际需求,确定滤波器的截止频率、通带波纹、阻带衰减和过渡带宽等参数。
这些参数决定了滤波器的性能。
4.设计滤波器的频率响应。
使用窗函数、最小二乘法等方法,根据滤波器的参数来设计滤波器的频率响应。
5.将频率响应转换为滤波器的系数。
根据设计的频率响应,使用逆快速傅里叶变换(IFFT)等方法将频率响应转换为滤波器的系数。
6.实现滤波器。
将滤波器的系数应用到数字信号中,实现滤波操作。
7.优化滤波器性能。
根据需要,可以对滤波器进行进一步优化,如调整滤波器的阶数、参数等,以达到较好的滤波效果。
以上是FIR数字滤波器的设计步骤,根据实际需求进行相应的调整,可以得到理想的滤波器。
如何设计和实现电子电路的数字滤波器

如何设计和实现电子电路的数字滤波器数字滤波器是电子电路设计中常用的一种模块,它可以去除信号中的不需要的频率分量,同时保留所需的信号频率。
本文将介绍数字滤波器的设计和实现方法。
一、数字滤波器的基本原理数字滤波器可以分为两大类:无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。
IIR滤波器的特点是具有无限长的脉冲响应,可以实现更为复杂的滤波功能;而FIR滤波器的脉冲响应是有限长的,适用于对频率响应要求较为严格的应用场景。
数字滤波器的设计思路是将模拟信号进行采样并转换为离散信号,然后利用差分方程实现各种滤波算法,最后将离散信号再次还原为模拟信号。
常见的离散滤波器有低通、高通、带通和带阻四种类型,根据不同的滤波需求选择合适的类型。
二、数字滤波器的设计步骤1. 确定滤波器类型和滤波需求:根据要滤除或保留的频率范围选择滤波器类型,确定截止频率和带宽等参数。
2. 选择合适的滤波器结构:基于具体需求,选择IIR滤波器还是FIR滤波器。
IIR滤波器通常具有较高的性能和更复杂的结构,而FIR滤波器则适用于对相位响应有严格要求的场景。
3. 设计滤波器的差分方程:根据所选滤波器结构,建立差分方程,包括滤波器阶数、系数等参数。
4. 系统状态空间方程:根据差分方程建立系统状态空间方程,包括状态方程和输出方程。
5. 计算滤波器的系数:根据差分方程或系统状态空间方程,计算滤波器的系数。
可以使用Matlab等专业软件进行系数计算。
6. 系统实现和验证:根据计算得到的系数,使用模拟或数字电路实现滤波器。
通过测试和验证,确保滤波器的性能符合设计要求。
三、数字滤波器的实现方法1. IIR滤波器实现方法:IIR滤波器可以通过模拟滤波器转换实现。
首先,将连续系统的模拟滤波器转换为离散滤波器,这一步通常使用差分方程实现。
然后,利用模拟滤波器设计的频响特性和幅频特性,选择合适的数字滤波器结构。
最后,通过转换函数将连续系统的模拟滤波器转换为数字滤波器。
数字滤波器设计方法

数字滤波器设计方法数字滤波器是数字信号处理中重要的一个组成部分,其作用是对数字信号进行滤波处理,消除噪声和干扰,提高信号的质量和可靠性。
数字滤波器的设计是数字信号处理中重要的一个环节,本文将介绍数字滤波器的设计方法及其步骤。
一、数字滤波器的设计方法数字滤波器的设计方法主要分为模拟滤波器设计法和数字滤波器设计法两种。
模拟滤波器设计法是在模拟域内进行滤波器设计,再将其转换为数字域中,而数字滤波器设计法是基于数字信号处理的理论和方法进行设计。
数字滤波器的设计方法可以分为两类,即基于时域设计和基于频域设计。
基于时域设计主要是对数字信号进行时域上的处理,通过调整滤波器传递函数中的系数来实现滤波器设计。
基于频域设计则是对频率响应进行优化设计,通过傅里叶变换将时域信号转换为频率域信号,进而对其进行频率响应设计。
在实际滤波器设计中,两种方法可以相互结合,实现更加灵活有效的数字滤波器设计。
二、数字滤波器设计的步骤数字滤波器设计主要包括以下步骤:1. 滤波器的性能评估首先要明确数字滤波器设计的目的和要求,如要过滤的信号频率范围、所要达到的滤波器性能指标和运算速度等。
在确定这些要素后,可以选择适当的滤波器设计方法和算法。
2. 数字滤波器的类型选择按照数字滤波器传递函数的形式,可将其分为FIR滤波器和IIR滤波器两种类型。
FIR滤波器是有限脉冲响应滤波器,具有线性相位和时域上的线性性质。
其优点在于简单可靠,易于实现,且滤波器响应的改变仅与滤波器系数有关,具有较好的稳定性和可重现性。
而IIR滤波器则是无限脉冲响应滤波器,其传递函数在分母中包含反馈因子,因此具有频域上的非线性性质。
IIR滤波器的优点是设计具有更快的计算速度和更窄的频带滤波器响应,但其稳定性和阶数选择需进行充分考虑。
3. 滤波器的设计在实际滤波器设计中,可以根据所选波形的性质来设计滤波器的系数。
根据所选择的滤波器类型和具体算法,可以采用各种滤波器设计工具进行滤波器系数计算。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
数字信号处理数字滤波器的设计学院计算机与电子信息学院专业电子信息科学与技术班级电子15-2 班学号指导教师利民数字滤波器的设计一、模拟低通滤波器的设计方法1、Butterworth 滤波器设计步骤: ⑴.确定阶次N① 已知Ωc 、Ωs 和As 求Butterworth DF 阶数N② 已知Ωc 、Ωs 和Ω=Ωp()的衰减Ap 求Butterworth DF阶数N③ 已知Ωp 、Ωs 和Ω=Ωp 的衰减Ap 和As 求Butterworth DF 阶数N3dB p Ω≠-/10/1022(/)101,(/)101p s A A N N p c s c ΩΩ=-ΩΩ=-则:⑵.用阶次N 确定 根据公式:在左半平面的极点即为的极点,因而2、切比雪夫低通滤波器设计步骤: ⑴.确定技术指标归一化: ⑵.根据技术指标求出滤波器阶数N 及:⑶.求出归一化系统函数 其中极点由下式求出:()a H s 1,2,2N()()a a H s H s -()a H s 2,,N p Ωp αs Ωs α/1p p p λ=ΩΩ=/s s p λ=ΩΩε0.12101δε=-p δα=或者由和S 直接查表得二、数字低通滤波器的设计步骤:1、确定数字低通滤波器的技术指标:通带截止频率pω、通带最大衰减系数pα、阻带截止频率ω、阻带最小衰减系数s α。
2、将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。
巴特沃斯:切比雪夫:N ()a H p /ss p λ=ΩΩ0.12101δε=-p δα=3、把模拟滤波器变换成数字滤波器,即把模拟滤波器的系数)(S H 映射成数字滤波器的系统函数)(z H 。
实现系统传递函数s 域至z 域映射有脉冲响应不变法和双线性映射两种方法。
(1)脉冲响应不变法。
按照技术要求设计一个模拟低通滤波器,得到模拟低通滤波器的传输函数()s H a 转换成数字低通滤波器的系统函数H(z)。
设模拟滤波器的传输函数为()s H a ,相应的单位冲激响应是()t h a ,()s H a =LT[()t h a ],LT[.]代表拉氏变换,对()t h a 进行等间隔采样,采样间隔为T ,得到()nT h a ,将h(n)= ()nT h a 作为数字滤波器的单位取样响应,那么数字滤波器的系统函数H(z)便是h(n)的Z 变换。
因此脉冲响应不变法是一种时域上的转换方法,它是h(n)在采样点上等于()t h a 。
设模拟滤波器()s H a 只有单阶极点,且分母多项式的阶次高于分子多项式的阶次,将()s H a 用部分分式表示:()∑=-=Ni iia s s A s H 1,式中i s 为()s H a 的单阶极点。
将()s H a 逆拉氏变换得到()t h a :())(1t u e A t h t s Ni i a i ∑== ,式中u(t)是单位阶跃函数。
对()t h a 进行等间隔采样,采样间隔为T,得到:())()(1nT u e A nT h n h nT s Ni i a i ∑===对上式进行Z 变换,得到数字滤波器的系统函数H(z):()∑=--=Ni T s iz eA z H i 111 , 经过一系列变换得到:)(1)(s kae z jk s HTz H sT Ω-=∑=(2)双线性变换法这种变换方法,采用非线性频率压缩方法,将整个频率轴上的频率围压缩到±π/T 之间,再用sTe z =转换到z 平面上。
设Ha (s ),s=j Ω,经过非线性频率压缩后用Ha (s1),1s =j Ω1表示,这里用正切变换实现频率压缩:()T T 15.0tan 2Ω=Ω式中T 仍是采样间隔,当Ω1从-π/T 经过0变化到π/T 时,Ω则由-∞经过0变化到+∞,实现了s 平面上整个虚轴完全压缩到1s 平面上虚轴的±π/T 之间的转换。
这样便有()t s Ts z z T T th T s 111125.021--+-=Ω=再通过sTe z = 转换到z 平面上,得到:11112--+-=z z T ss T s T z -+=22令ωj e z j s =Ω=, ,有ωωj j e e T j --+-=Ω112 ω21tan 2T =Ω 两种方法比较:脉冲响应不变法的优点:1,模拟频率到数字频率的转换时线性的;2,数字滤波器单位脉冲响应的数字表示近似原型的模拟滤波器单位脉冲响应,因此时域特性逼近好。
缺点:会产生频谱混叠现象,只适合带限滤波器 双线性变换法优点:克服多值映射得关系,可以消除频率的混叠 缺点:时域到频域的变换是非线性的,在高频处有较大的失真。
三、数字高通滤波器的设计步骤:① 数字高通滤波器的技术指标为:通带截止频率ωp 阻带截止频率ωs 通带衰减频率αP 阻带衰减频率 αs2/p p sw f F π=2/s s sw f F π=② 预畸变处理,将数字高通指标转换为模拟低通指标2cot()2s s w T Ω=2cot()2p p w T Ω=④ 确定阶数N 由1010lg((101)(101))2lg()ats p N δ--=ΩΩ (可由模拟低通滤波器设计方法可得) ④ 归一化及去归一化 查表令s=s/Ω 归一化模拟低通圆型系统函数G (P )=12+√2P +1⑤ 低通向高通转化 令s1=1/s 由频率变换公式 H a (s )=G (ρ)|ρ=ρρρρρρ 即可得 ⑥ 滤波器数字化 令1-1-z -1z 12s +=T 利用双线性变换化 H (z )=H a (s )|ρ=1−ρ−11+ρ−1带入数据可得数字高通H(z)数字高通不能采用脉冲响应不变法原因是:脉冲响应不变法有频谱周期延拓效应,因此只能用于带限的频响特性,如衰减特性很好的低通或带通。
而高频衰()a H s减越大,频响的混淆效应越小,至于高通滤波器,由于它们在高频部分不衰减,因此将完全混淆在低频响应中。
四、数字带通滤波器的设计:步骤:(1)确定性能指标: 在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标:通带截止频率wc1,wc2、阻带截止频率wr1,wr2、阻带最小衰减αs通带最大衰减αp(2) 对带通数字滤波器H(z)的数字边界频率预畸变Ω=2T∗tan(w2)得到带通模拟滤波器H(s)的边界频率主要是通带截止频率ωp1,ωp2;阻带截止频率ωs1,ωs2的转换。
对双线性变换法一般T=2s通带截止频率wc1=(2/T)*tan(wp1/2) 、wc2=(2/T)*tan(wp2/2) 阻带截止频率wr1=(2/T)*tan(ws1/2)、wr2=(2/T)*tan(ws2/2)阻带最小衰减αs通带最大衰减αp(3)低通到带通频率变换λ=((Ω2)−(Ω02)B∗Ω)将模拟带通滤波器指标转换为模拟低通滤波器指标。
B=wc2-wc1normwr1=(((wr1^2)-(w0^2))/(B*wr1)) normwr2=(((wr2^2)-(w0^2))/(B*wr2)) normwc1=(((wc1^2)-(w0^2))/(B*wc1)) normwc2=(((wc2^2)-(w0^2))/(B*wc2)) 模拟低通滤波器指标: normwc ,normwr ,αp,αs(4)设计模拟低通原型滤波器。
查表得到归一化低通传输函数G(p):G (P )=12+√2P +1用模拟低通滤波器设计方法(由巴特沃斯设计步骤或切比雪夫设计步骤)得到模拟低通滤波器的传输函数Ha(s)(5)模拟低通滤波器转化为模拟带通滤波器。
)(202)()(l u s s p a p G s H Ω-ΩΩ+==(6)利用双线性变换法将模拟带通滤波器Ha(s)转换成数字带通滤波器H(z)由1-1-z -1z 12s +=T 11112202)(--+-=Ω-ΩΩ+=z z s l u s s p 得到 H (z )=H a(s )|ρ=1−ρ−11+ρ−1也可以用脉冲响应不变法:)(1)(s kae z jk s HTz H sT Ω-=∑=两种方法比较比较:脉冲响应不变法数字滤波器单位脉冲响应的数字表示近似原型的模拟滤波器单位脉冲响应,因此时域特性逼近好。
但会产生频谱混叠现象,只适合带限滤波器 双线性变换法可以克服多值映射得关系,可以消除频率的混叠但时域到频域的变换是非线性的,在高频处有较大的失真。
五、数字带阻滤波器的设计:步骤:(1) 确定性能指标:通带截止频率wc1,wc2、阻带截止频率wr1,wr2、阻带最小衰减αs 通带最大衰减αp(2) 对带通数字滤波器H(z)的数字边界频率预畸变Ω=2T ∗tan (w 2) 主要是通带截止频率ωp1,ωp2;阻带截止频率ωs1,ωs2的转换。
对双线性变换法一般T=2s通带截止频率wc1=(2/T)*tan(wp1/2) 、wc2=(2/T)*tan(wp2/2)阻带截止频率wr1=(2/T)*tan(ws1/2)、wr2=(2/T)*tan(ws2/2)阻带最小衰减αs 通带最大衰减αp(由模拟低通滤波器设计方法可得)模拟低通滤波器确定模拟带阻滤波器由模拟低通到模拟带阻的变换 这一模拟低通到带阻的变换关系为()a H ss =Ω02PP 2+Ω02 式中s 为模拟低通原型拉普拉斯变量(s=σ+j Ω),ρ为模拟带阻的拉普拉斯变 量(ρ=σ +j Ω ),Ω0是模拟带阻滤波器的几何中心频率。
令ρ=j Ω可得S =j Ω02ΩΩ02−Ω2故ρ平面的虚轴与s 平面的虚轴相对应,代入s=j Ω,消去j ,可得Ω=n 02ΩΩ02−Ω2低通的阻带映像到带阻的阻带ΩC =Ω02Ω1Ω0−Ω1−ΩC =Ω02Ω2Ω02−Ω22 化简得到:Ω0=Ω1Ω2̅̅̅̅̅̅̅̅̅=Ω2−Ω1=Ω1Ω3ΩC在H ρρ(ρ)中的变换关系,可得到带阻滤波器系统函数ρρρ p =ρρρ ρ |ρ=ρ02ρ(ρ2+ρ02)⁄ 由模拟带阻到数字带阻的变换 仍利用双线性变换P =C 1+z −1̅̅̅̅̅̅̅̅̅1−z −1模拟低通原型滤波器的ρ平面变换成数字带阻滤波器的ρ平面的表达式可得S=D1(1−ρ−1) 1−ρ1ρ−1+ρ−2从模拟低通系统函数Hρρ(ρ),转换数字带阻系统函数HBR p=HLP s|S=D1(1−ρ−1)1−ρ1ρ−1+ρ−2数字带阻滤波器不能用脉冲响应不变法:原因是脉冲响应不变法有频谱周期延拓效应,因此只能用于带限的频响特性,如衰减特性很好的低通或带通。
而高频衰减越大,频响的混淆效应越小,至于带阻滤波器,由于它们在高频部分不衰减,因此将完全混淆在低频响应中。