2019-数字信号处理6.4双线性变换法-文档资料
数字信号处理第六章6 双线性变换法
1)线性相位模拟滤波器
非线性相位数字滤波器
2)要求模拟滤波器的幅频响应为分段常数型,不 然会产生畸变 分段常数型模拟滤波器 经变换后仍为分段常数 型数字滤波器,但临界 频率点产生畸变
1 1 / T
2 tg
2012-10-11
1
1 1 c
数字信号处理
s1T
j
1T 2
e e
2 s1T 2
e e
s1T 2
s1T 2
e
s1 T 2
e
1 e 1 e
s1T s1T
1 z 1 z
1 1
z e
s1T
s
1 z 1 z
1 1
z
1 s 1 s
数字信号处理
2012-10-11
为使模拟滤波器某一频率与数字滤波器的任一
预畸变
给定数字滤波器的截止频率 1 ,则
1 c tg
1
2
按 1设计模拟滤 波器,经双线性 变换后,即可得 到 1 为截止频率 的数字滤波器
2012-10-11
数字信号处理
6、模拟滤波器的数字化方法
H (z) H a (s) 1 z Ha c 1 1 z
频率有对应关系,引入系数 c
c tg 1T 2
1 z 1 z
1 1
s c
z
c s cs
2012-10-11
数字信号处理
2、变换常数c的选择
1)低频处有较确切的对应关系:
1T 1T 1 c tg c 2 2
s1T
双线性变换法的原理
双线性变换法的原理
双线性变换法是一种通过将问题转化成一对线性方程组求解的方法,常用于解决二元二次方程或二元二次函数的问题。
其原理可以归纳如下:
1. 假设我们要解决一个二元二次方程或二元二次函数的问题,形式为ax² + bxy + cy² + dx + ey + f = 0(或f(x, y) = 0)。
2. 首先,对于该方程的每一项,我们引入一个新的变量u和v,并将该项表示为一个新的线性方程。
例如,对于ax²,我们将
其表示为au²。
3. 在引入新的变量后,我们得到了一组新的线性方程,形式为Aui + Bvi + Ci + Di + Ei + F = 0,其中i表示第i个线性方程。
4. 接下来,我们要构造一组满足上述线性方程的两个二次式,即f(u, v) = 0。
这里,我们选择f(u, v) = Au² + Buv + Cv² + Du
+ Ev + F。
5. 由于方程组中的每一个线性方程都对应一个二次式,我们可以得到关于u和v的二元二次方程。
我们需要求解这个二元二次方程,从而得到u和v的值。
6. 一旦找到了u和v的值,我们可以将其代入到原方程中,得到x和y的值,从而解决了原始的二元二次方程或二元二次函数问题。
双线性变换法的核心思想是通过引入新的变量,将一个二次式转化为一组线性方程,从而将原问题转化为一对线性方程组,利用线性方程组的解法来求解原问题。
这种方法的优势在于可以利用线性方程组求解的方法解决二次方程或二次函数的问题,而线性方程组求解的方法已经非常成熟和广泛应用。
双线性变换法
双线性变换法双线性变换法(bilinear transofrmation method)是一种通过变换以分析和解决非线性系统的复杂方法。
它最初由Collins,Mitroff和Zinnes提出,其主要特点是将非线性系统转化为线性系统来进行分析。
它把一个非线性系统映射到一个线性系统可以使一些复杂的非线性图像变成简单的线性图像,从而形成简单的表达式来解决复杂的问题。
一、双线性变换法定义双线性变换法是指通过线性常数和相关系数,将一维和多维数据变换为更简单的线性形式,以模拟复杂的非线性系统的运算的一种变换方法。
二、双线性变换法的应用(1)控制论领域。
双线性变换可以将复杂的非线性系统转变为简单的线性系统,使得这些复杂的系统容易控制。
(2)视觉领域。
双线性变换可以解决计算机视觉中的误差传播问题,将非线性的图像识别问题转变为简单的线性问题来处理;另外,在图像处理领域用双线性变换可以实现图像的变换,从而实现复杂的图像变换;(3)机器学习领域。
双线性变换可以将非线性的机器学习问题变换为线性的问题,让算法可以更加简单有效地解决复杂的机器学习问题。
三、双线性变换法的局限性(1)双线性变换法还有一些困难。
例如,当非线性系统出现很多两个变量或多个变量间有联系时,双线性变换也会受到很大影响。
(2)双线性变换法也会遇到数值不稳定的问题,在遇到非线性系统的情况下,很多变量的变化对结果的影响会变得很大,因此会产生数值不稳定的现象。
(3)双线性变换只是一种模拟,它并不能完全模拟出非线性系统的真实行为,因此很多时候双线性变换的结果可能不太准确。
双线性变换法是一种实用性很强的方法,它可以帮助我们更准确地分析和解决非线性系统问题,它也应用于控制论、视觉和机器学习等领域,但由于它有一些限制,如数值不稳定性和无法完全模拟非线性系统,因此我们需要更加谨慎地运用双线性变换法来真正发挥它的优势。
双线性变换法
上节介绍的冲激响应不变法是一种线性映射的 设计方法,其缺点是会产生频率响应的混叠失 真,这主要是因为从s平面到z平面的映射是多值 (多对一)的映射关系。双线性变换法就是从 克服混叠失真的角度出发,寻找从s平面到z平面 的单值映射关系。 其缩射是原的单s理1平值是面的先向,将z可s平平以面面很进进好行行的单压克值缩服映至频射s1率。平响由面应于,的这再混种将叠映压 失真,但值得注意的是,这种映射是非线性的
根据Ω和Ω1, ω的关系,由(5.21)式,及ω= Ω1T 的关系式可得到双线性变换法的模拟角频率Ω和 数字频率ω之间的变换关系为
c tan( )
2
(5.26)
它表示 是单值的一一对应的映射关系
关键频点的对应关系:
Ω
ω
∞
π
0
0
-∞
-π
12
§ 5.4 双线性变换法(频域逼近法)
上述对应关系优点:有效避免了频谱混叠
于数字滤波器的低频特性
7
§ 5.4 双线性变换法(频域逼近法)
2、选择c,使数字滤波器的某一特定频率
(如截止频率ωc=Ω1cT)与模拟滤波器的一 个特定频率Ωc=2fc严格对应,即
c
ctan( 1cT ) 2
ctan(c )
2
此时有
c
cctg
(c
2
)
这种选择的优点在于:在特定的模拟频率和特定 的数字频率处,有严格相等的频率响应,因而可 以较准确地控制截止频率的位置。
s
c
1 1
z z
1 1
c
1 1
e e
j j
j.c
tan(
)
2
双线性变换法公式
双线性变换法(Bilinear Interpolation)是在图像处理中常用的一种插值方法。
公式如下:
f(x,y) = (1-x)(1-y)f(0,0) + (1-x)yf(0,1) + x(1-y)f(1,0) + xyf(1,1)
其中x,y 为目标像素坐标在原图像坐标系中的坐标值,f(0,0),f(0,1),f(1,0),f(1,1) 分别表示目标像素周围4 个像素点的灰度值。
双线性变换法是一种通过线性变换来求解目标像素点灰度值的方法。
它通过对图像进行缩放或旋转时,对于输出图像中缺失的像素点进行插值,来解决图像变形导致的像素点缺失问题。
双线性变换法是一种非常高效的插值方法,其计算量与像素点数量无关。
另外,它还具有较高的精度和较低的计算复杂度。
它在图像处理、图像识别、图像分析、图像压缩等领域有着广泛的应用。
双线性变换法是一种双线性插值法,它基于线性插值法,通过对目标像素周围4个像素点的灰度值进行线性变换来求出目标像素点的灰度值。
其优点是插值效果好,像素质量高,图像变形较小。
双线性变换法在图像缩放、旋转、矫正等操作中都有着广泛的应用。
它在图像处理中常用来解决图像变形导致的像素点缺失问题。
此外还可以用于从低分辨率的图像中重建高
分辨率图像,并且在视频处理中也有着广泛的应用。
双线性变换法原理的解释
s=
2 T
1 - z- 1 1 + z- 1
(5)
则可以直接从 H i (s) 得到对应的 H i (z )。
114 高阶连续时间系统的离散实现
由于一阶环节是高阶连续时间系统的最基本环
节, 若每个一阶环节都用式 (5) 进行了离散化, 则整
个系统都被离散化了。
至此 , 用双线性变换公式 (5) 设计数字滤波器
图 1 连续系统一阶环节
112 积分的数值计算与离散一阶系统
一次积分运算可以用梯形法作数值计算, 即
∫nT
y -∞
∫ ∫ (n- 1) T
nT
x (Σ) dΣ+
x (Σ) dΣ=
-∞
(n- 1) T
∫nT
y [ (n - 1) T ] +
x (Σ) dΣ
(n- 1) T
的原理就比较清晰了。
2 结束语
对于双线性变换法的原理还有其他一些重要的 概念需作介绍和讨论, 如连续时间角频率和离散时 间角频率之间的映射关系等, 一般在现有教材中均 有介绍, 此处不再叙述。
1 双线性变换原理的介绍
111 连续时间系统 H (s) 的最基本环节 连续时间系统 H (s) 的极点有两种情况: 单重节
点和多重节点。 但是一个多重极点环节可以看成由 多个单重极点环节级联构成, 例如对二重极点有
收稿日期: 2003- 11- 24; 修回日期: 2004- 01- 20 作者简介: 丁志中 (1961- ) , 男, 安徽省芜湖人, 硕士, 讲师, 主要从事信号分析与处理、通信等方面的教学与研究。
0 引言
随着M A TLAB 应用的普及, 本科生学习和掌 握数字滤波器设计的原理和方法已经逐渐成为广泛 的需求。 但是离散时间信号和系统的内容理论性较 强, 概念较抽象, 同时初学者对信号与系统中的许多 问题理解还不是很深入, 因此如何将数字滤波器设 计中一些问题的物理概念阐述得简单而清晰一些, 对于该部分内容的教学就显得更为重要。
数字信号处理 第6章
H ( z ) h( n) z n
n 0
N 1
(6.1.2)
(6.1.1)式中的H(z)称为N阶IIR数字滤波器系统函数; (6.1.2) 式中的H(z)称为N-1阶FIR数字滤波器系统函数。这两种 数字滤波器的设计方法有很大区别,因此下面分成两章分 别进行学习。
第6章 无限脉冲响应数字滤波器的设计
s 20 lg
| H (e j0 ) |
j s
dB
(6.1.4b)
p 20 lg | H (e
j p
) | dB
(6.1.5)
s 20 lg | H (e js ) | dB
(6.1.6)
第6章 无限脉冲响应数字滤波器的设计
当幅度下降到 2 / 2 时,标记ω=ωc,此时 p 3dB,称 ωc为3 dB通带截止频率。ωp、ωc和ωs统称为边界频率, 它们是滤波器设计中所涉及到的很重要的参数。对其他 类型的滤波器,(6.1.3b)式和(6.1.4b)式中的H(ej0)应改 成
拟滤波器得到系统函数Ha (s),然后将Ha(s)按某种方法转
换成数字滤波器的系统函数H(z)。这是因为模拟滤波器的 设计方法已经很成熟,不仅有完整的设计公式,还有完善
的图表和曲线供查阅; 另外,还有一些典型的优良滤波
器类型可供我们使用。直接法直接在频域或者时域中设计 数字滤波器,由于要解联立方程,设计时需要计算机辅助 设计。FIR滤波器不能采用间接法,常用的设计方法有窗 函数法、频率采样法和切比雪夫等波纹逼近法。
第6章 无限脉冲响应数字滤波器的设计
图6.1.3所示的单调下降幅频特性,p和s别可以表
示为
p 20 lg
| H (e j0 ) | | H (e
数字信号处理 双线性变换法
河北XX大学课程设计报告学生姓名:学号:专业班级:课程名称:数字信号处理学年学期:学年第学期指导教师:2 0 11 年月课程设计成绩评定表目录一、设计题目二、设计目的三、设计原理四、实现方法五、设计内容及要求六、设计结果及改进建议七、思考题八、设计体会九、参考文献一、设计题目用双线性变换法设计IIR数字滤波器二、设计目的1.熟悉IIR数字滤波器的原理与方法。
2.掌握数字滤波器的计算机仿真方法。
3.通过观察对实际心电图信号的滤波作用,获得数字滤波的感性知识。
三、设计原理1、双线性变换法采用非线性频率压缩方法,将整个频域轴上的频率范围压缩到-π/T~π/T之间,再用z=e sT转换到z平面上。
也就是说,第一步现将整个S平面压缩映射到S1平面的-π/T~π/T一条横带里;第二步再通过标准变换关系1e S T将此横带变换到整个z平面上去。
这样就使S平面与Z平面建立了一一对应的胆汁关系,消除了多只变换性,也就消除了频谱混叠现象,映射关系如图所示。
Z平面为了将S平面的整个虚轴jΩ压缩到S平面轴上的-π/T到π/T段上,1可以通过以上的正切的变换实现Ω=2/Ttan(1ΩT/2)(1-1)式中,T 仍是采样间隔。
当1Ω由-π/T 经过0变化到π/T 时,Ω由-∞经过0变化到+∞,也即映射了整个j Ω轴。
将上式(1-1)写成111122222TTj j TTj j e e j TeeΩΩΩΩ--Ω=∙+ (1-2)将此关系解析延拓到整个S 平面和1S 平面,令j Ω=s, 1j Ω=S1,则得111221tan()21es Ts TS T e TT--==∙+(1-3)再将S1平面通过以下标准变换关系映射到Z 平面:z=S1T e (1-4) 从而得到S 平面和Z 平面的单值映射关系为;11211Z S T Z---=+(1-5) 122122T T S SZ T T SS++==--(1-6)式(1-5)和式(1-6)是S 平面与Z 平面之间的单值映射关系,这种变换都是两个线性函数之比,因此成为双线性变换。
数字信号处理-双线性变换法
(2)2 cosh1(s/c) 5.5338 0.50885
1
1 =
cosh1
1 101.51
cosh1(tan0.15 /tan.2)
=
3.016 = cosh1(1.0196/0.65) 1.0207
cosh1
3.0783
取 N=4,与脉冲不变法相同。
= 1 + 1+ 2 =1.9652+2.205=4.1702
这样模拟滤波器的设计指标为
通带截止频率p =0.65,通带最大衰减p ≤1dB; 阻带边缘频率s =1.019,阻带最小衰减பைடு நூலகம்s ≥15dB; 1 )求 N 、 c 10lg|H (jp)|2 ≥ 1 , lg|H (jp)|2 ≥ 0.1,
|H (j 2 tan(0.1))|2 ≥ 10-0.1 ; 10lg|H (js)|2 ≥ 15 , |H (j 2 tan(0.15))|2 ≥10-1.5。 lg|H (js)| 2 ≥ 1.5,
j j/2 j/2) (e e 2 2 j2 = 1 =j tan 2 T T j /2 j /2 +e ) 2 (e
2 + j = j tan 2 T
比较(6.4-8)等式两边,得到
(6.4-8)
=0
2 = tan 2 T
(6.4-9)
由(6.4-8)看到双线性变换法的映射关系使s平面的虚轴 映射为z平面的单位圆。 而(6.4-9)式频率正切变换关系 实现了频率压缩,使模拟域从~的变化,压缩为 数字域频率 从 ~ 变化。
函数H(z)。
与脉冲不变法一样,设计过程中除了的第一步求数字临 界频率{k}时,要用到取样间隔T或取样频率 fs 以外, 最后的结果与其它各步骤中T 或 f s的取值无关。所以为
数字信号处理双线性变换法
取T=2,利用
= 2 tan( ) 得模拟带阻指标为
T2
p1=6rad, p2=13rad, s1=9rad, s2=11rad,p1dB, s 10dB
第22页/共39页
例:试设计满足下列指标的BW型数字带阻滤波器
p1=2.8113rad/s, p2=2.9880rad/s, p1dB , s1=2.9203rad/s, s2=2.9603rad/s, s 10dB
1 sT 1
2 tan(p / 2)
第10页/共39页
例:用双线性变换法和一阶巴特沃思低通滤波器,设
计一个3dB截频为p的数字滤波器,并与脉冲响
应不变法设计的DF比较。
解:设双线性变换中的参数为T
(3) 用双线性变换法将模拟滤波器转换为数字滤波器
H (s) =
1 sT 1
2 tan(p / 2)
s
计一个3dB截频为p的数字滤波器,并与脉冲响
应不变法设计的DF比较。
解:设双线性变换中的参数为T
(1) 将DF的频率指标转换为AF的频率指标
p
=
2 T
tan(p
2
)
(2) 设计3dB截频为p的一阶BW型模拟低通滤波器,即
N=1, c = p
故
H (s) = 1 = 1 = s / c 1 s / p 1
s 4 198s 2 9801
HBS(s) = Ha (s)
s=
s2
Bs 02
=
s4
4.899s3
210s 2
485s 9801
(5) 由双线性变换模拟带阻滤波器转换成数字带阻滤波器
H (z) = H (s) s= 2 1z1 T 1 z1
双线性变换法
二、性能分析4
• 6.对于分段常数型AF滤波器,经双线性变换后, 仍得到幅频特性为分段常数的DF.但在各个分段 边缘的临界频率点产生畸变,这种频率的畸变, 可通过频率预畸变加以校正。
例1
• 一个线性相位的模拟滤波器经双线性变换后 得到非线性相位的数字滤波器,不再保持原 有的线性相位。如一个模拟微分器将不能通 过双线性变换成为数字微分器。
双线性变换法
• 冲激不变法(和阶跃响应):是使数字滤波器 在时域上模仿模拟滤波器,但它的缺点:产生 频率响应的混叠失真。这是由于从S平面->Z平 面是多值的映射关系所造成的。为了克服这一 缺点,我们采用双线性变换法
2
一、变换原理 1、定义
• 双线性变换法:是从频域出发,使DF的频率响 应与AF的频率响应相似的一种变换法。
(3)变换常数C的选择1
• 调节C,可使AF与DF在不同频率点处有对应的关 系。 • (a)使AF与DF在低频处有较确切的对应关系。
T 1T 1T 1T T tan( ) 2 2 2 2 1T 由 C tan( )C 2 2 2 C T
• 看出在低频处,AF的低频特性近似等于DF的低 频特性。
(3)变换常数C的选择2
(b)利用DF 的某一特定频率(例截止频率 c ) 与AF的某一特定频率 c 严格相对应。 cT 即: c C tan( ) 2 cT 2f c c c ctg ( ) c ctg 2 2 fs f c c cot( ) fs 看出:此方法优点:是在特定AF和特定DF 处, 频率响应是严格相等的, 它可以较准确地控制截止频率的位置。
) |在低频处
ze
s1T
(1)频率压缩
数字信号处理用双线性变换法
数字信号处理用双线性变换法双线性变换双线性变换法双线性变换公式matlab双线性变换线性变换线性变换的矩阵分段线性变换直接线性变换灰度线性变换
一、实验目的
二、实验原理
三、实验步骤
实实验体会
通过这次实验,进一步掌握了双线性变换法设计IIR滤波器的原理和方法,也更熟悉了数字滤波器仿真的方法,对数字滤波器的图形有了一些认识。
chap6-4双线性变换法资料
s1T
e e e e 11 ee j1T
j1T
s1T
s1T
s1T
e 2 e 2 e 2 e 2
(2)数字化
• 将S1平面通过标准变换关系变换到z平面。
z es1T , 此时,s1 z 则可得到s平面 z平面的单值映射关系:
s
1 1
z 1 z 1
或z
1 1
s s
实际中,为了使模拟滤波器的某一频率与数字化滤波器
7所示。
由图6-7看出,在零频率附近,模拟角频率Ω与数字频率ω之 间的变换关系接近于线性关系;但当Ω进一步增加时,ω增长得 越来越慢,最后当Ω→∞时,ω终止在折叠频率ω=π处,因而双线 性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分 去的现象, 从而消除了频率混叠现象。
=
2 T
tan
(2) 按Ω1、Ω2、Ω3和Ω4等指标设计模拟滤波器的系统函数Ha(s)。
(3)将
s
c
1 1
z z
1 1
代入Ha(s),得H(z)为
H
(z)
Ha
(s)
|s
c1 1
z z
1 1
Ha
c
1 1
z 1 z 1
其频率响应为
H (e j )
Ha (
j)
|
c
tan
2
Ha
jc
tan
2
上述这些步骤比用脉冲响应不变法设计滤波器要简便得多。
这样得到的模拟滤波器h经双线性变换后就映射到数字滤波器hz的原转折频率如果给出的是待设计的带通滤波器的模拟域转折频率通阻带截止频率f和采样频率1t则首先利用下式计算数字滤波器的转折频率通阻带截止频率与给定的模拟域转折频率f上述这些步骤比用脉冲响应不变法设计滤波器要简便得多
双线变换法数字低通滤波器word精品文档15页
燕山大学课程设计说明书题目:双线性变换法设计数字低通滤波器学院(系):电气工程学院年级专业: 10级检测2班学号: 100103020192学生姓名:刘培露指导教师:王娜教师职称:讲师课程名称:数字信号处理课程设计基层教学单位:仪器科学与工程系指导教师:王娜学号学生姓名(专业)班级设计题目12、双线性变换法设计数字低通滤波器设计技术参数给定技术指标为:Hzfp100=,Hzfs300=,dBp3=α,dBs20=α,采样频率HzFs1000=。
设计要求设计Butterworth低通滤波器,用双线性变换法转换成数字滤波器。
(buttord,buttap,butter,bilinear)参考资料数字信号处理方面资料MATLAB方面资料周次前半周后半周应完成内容收集消化资料、学习MA TLAB软件,进行相关参数计算编写仿真程序、调试指导教师签字基层教学单位主任签字说明:1、此表一式四份,系、指导教师、学生各一份,报送院教务科一份。
2、学生那份任务书要求装订到课程设计报告前面。
电气工程学院教务科目录一.课题描述 (2)二.设计原理 (2)2.1 IIR数字滤波器设计原理 (2)2.2巴特沃斯低通滤波器的原理 (3)2.3双线性变换法 (4)2.4用双线法设计巴特沃斯数字低通滤波器的步骤 (9)三. MATLAB程序 (10)四.程序中命令介绍 (12)五. 图像结果……………………………………………………………………………. ..14六. 结果分析 (16)七.心得体会: (16)一.课题描述数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。
可以设计系统的频率响应,让它满足一定的要求,从而对通过该系统的信号的某些特定的频率成分进行过滤,这就是滤波器的基本原理。
如果系统是一个连续系统,则滤波器称为模拟滤波器。
如果系统是一个离散系统,则滤波器称为数字滤波器。
数字滤波实质上是一种运算过程,实现对信号的运算处理。
6-4数字信号处理
2 1− z −1 s= ⋅ T 1+ z −1
第7页
2 1 − z −1 = Ha ⋅ T 1 + z −1
jω 频响: 频响 H ( e ) = H a ( jΩ )
ω 2 = H a j ⋅ tan 2 ω 2 Ω = ⋅ tan T
(4)将模拟滤波器 a(s),从s平面转换到 平面,得到数字低 将模拟滤波器H , 平面转换到 平面, 平面转换到z平面 将模拟滤波器 通滤波器系统函数H(z)。 。 通滤波器系统函数
X
第9页
参数采样间隔T的选取 参数采样间隔 的选取 脉冲响应不变法 避免频率混叠,要求所设计的模拟低通带限于± 避免频率混叠,要求所设计的模拟低通带限于±π/T 可选择T 可选择 满足 Ω s < 最高频率在ω=π 按照线性关系 Ω s = 那末一定满足 Ω s < 双线性变换法
通过标准变换关系将此横带变换到整个z由此建立s平面与z平面一一对应的单值关系消64用双线性变换法设计iir数字低通滤波器也因果稳定因果稳定从t经过0变化到t时则由经过0变化到实现了s平面上整个虚轴完全压缩到s平面上虚轴的t之间的转换
第1页
6.4 用双线性变换法设计 用双线性变换法设计IIR数字 数字 低通滤波器
X
当
经过0变化到 经过 变化到π/T时, 时 1 从 -π/T经过 变化到
则由-∞经过 变化到 则由 经过0变化到 经过
第3页
+∞,实现了 平面上整个虚轴完全压缩到 1平面上虚轴 ,实现了s平面上整个虚轴完全压缩到 平面上整个虚轴完全压缩到s 之间的转换。 的±π/T之间的转换。 之间的转换 将这一关系解析扩展至整个S平面, 将这一关系解析扩展至整个 平面,则得到 平面 S平面到 1平面的映射关系: 平面到S 平面的映射关系: 平面到
双线性变换法
H ( z) H (s)
双线性变换法
s
2 1 z 1 T 1 z 1
双线性变换法设计DF的步骤
w tan( ) 2 w ,w Wp,Ws T p s
2
W
设计模拟 滤波器
H(s)
双线性变换
H(z)
利用MATLAB [numd,dend] = bilinear(num,den,Fs)
T
W
H(e )
jW
Wp
双线性变换法
Ws
W
双线性变换法的基本原理
双线性变换法的优缺点 优点:无混叠 缺点:幅度响应不是常数时会产生幅度失真
双线性变换法
双线性变换法设计DF的步骤
1. 将数字滤波器的频率指标{Wk}转换为 模拟滤波器的频率指标{wk} Wk 2 w k tan( ) T 2 2. 由模拟滤波器的指标设计模拟滤波器的H(s)。 3. 利用双线性变换法,将H(s)转换H(z)。
1
例: 利用BW型模拟低通滤波器和双线性变换法设计满足指 标Wp=p/3,Ap=3dB,N=1的数字低通滤波器,并与 脉冲响应不变法设计的DF比较。
解:设双线性变换中的参数为T
(3) 用双线性变换法将模拟滤波器转换为数字滤波器 1 H (s) sT 1 2 tan( W p / 2 ) 2 1 z 1 s T 1 z 1 tan(W p / 2)(1 z 1 ) 0.366 0.366z 1 H ( z) 1 1 0.2697z 1 1 tan(W p / 2) (tan( W p / 2) 1) z 结论:参数T的取值和最终的设计结果无关。 为简单起见,一般取T=2
3dB
Amplitude
脉冲响应不变法 脉冲响应不变法存在频 谱混叠,所设计的DF不满 足给定指标。而双线性变 换法不存在频谱混叠,所 设计的DF满足给定指标。
(完整word版)双线性变换法
脉冲响应不变法的主要弊端是频谱交叠产生的混杂,这是从S 平面到 Z 平面的标准变换z =e sT的多值对应关系致使的 , 为了战胜这一弊端,假想变换分为两步:第一步:将整个 S 平面压缩到 S1 平面的一条横带里;第二步:经过标准变换关系将此横带变换到整个Z 平面上去。
由此成立 S 平面与 Z 平面一一对应的单值关系,除去多值性,也就除去了混杂现象。
图双线性变换的映照关系为了将 s 平面的 j Ω轴压缩到 s1平面 j Ω轴上的 - 一段上,可经过以下的正切变换实现:这里 C 是待定常数,下边会讲到用不一样的方法确立C,可使模拟滤波器的频次特征与数字滤波器的频次特征在不一样频次点有对应关系。
时经过这样的频次变换,当Ω 1由Ω 由即映照了整个 j Ω轴。
将这一关系分析延拓至整个s 平面,则获得 s 平面平面的映照关系:再将 s1平面经过标准变换关系映照到z 平面,即令往常取C=2/T最后得 S 平面与 Z 平面的单值映照关系:此刻我们再来看一看常数 C 的取值方法:双线性换法的主要长处是S 平面与 Z 平面一一单值对应, S 平面的虚轴 ( 整个 j Ω) 对应于Z 平面单位圆的一周, S 平面的Ω=0 处对应于 Z 平面的ω =0 处,对应即数字滤波器的频次响应终止于折迭频次处,所以双线性变换不存在混迭效应。
上边讲到,用不一样的方法确立待定常数C,能够使模拟滤波器的频次特征与数字滤波器的频次特征在不一样频次点有对应关系。
也就是说,常数 C 能够调理频带间的对应关系。
确立 C 的常用方法有两种:①保证模拟滤波器的低频特征迫近数字滤波器的低频特征。
此时二者在低频处有切实的对应关系,即因为Ω和ω都比较小,所以有,所以有Ω =cΩT/2 ,此外,依据归一化数字频次ω与模拟频次Ω的关系,所以, c=2/T②保证数字滤波器的某一特定频次,如截止频次,与模拟滤波器的某一待定频率Ωc 严格对应,即当截止频次较低时,有,所以一般取。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
解:
双线性变换法设计的DF的系统函数为
H 双 (z)=1tan( ta pn /(2 )p /(2 t)a(n 1( z p/1 2 ))1)z1
脉冲响应不变法设计的DF的系统函数为
1ep H脉(z) =1ep z1
取p =0.6p,令z=ej ,可分别获得两者的幅度响应
计高通、带阻等滤波器。
双线性变换法
双线性变换法的基本原理
和 的关系
= 2 tan() =2arctan(T)
T2
2
2 1 z1 s = T 1 z1
2 s
z
=
T 2
s
T
双线性变换法
双线性变换法的基本原理
稳定性分析
s
=
2 T
1 1
z1 z1
z=2Ts 2 T s
令s=+j,则有
z=
(2/T)2 2 (2/T)2 2
双线性变换法
双线性变换法的基本原理
稳定性分析
(2/T)2 2
z=
(2/T)2 2
1) <0, |z|<1 S域左半平面映射到Z域单位圆内 2) =0, |z|=1 S域虚轴映射到Z域单位圆上
3) >0, |z|>1 S域右半平面映射到Z域单位圆外
1
sT
1
2 tan( p / 2)
例:用双线性变换法和一阶巴特沃思低通滤波器,设计一
个3dB截频为p的数字滤波器,并与脉冲响应不变法
设计的DF比较。
解:设双线性变换中的参数为T
(3) 用双线性变换法将模拟滤波器转换为数字滤波器
H(s) =
1
sT 1
2tan(p / 2)
s
=
2 T
1 1
利用MATLAB
H(z) = H(s) s=21z1 T 1z1
[numd,dend] = bilinear(num,den,Fs)
num,den:AF分子、分母多项式的系数向量
Fs:抽样频率
numd,dend:DF分子、分母多项式的系数向量
双线性变换法
例:用双线性变换法和一阶巴特沃思低通滤波器,设计一
z1 z1
H (z)=1tan( ta pn /(2 )p /(2 t)a(n 1( z p/1 2 ))1)z1
结论:参数T的取值和最终的设计结果无关。
为简单起见一般取T=2
例:用双线性变换法和一阶巴特沃思低通滤波器,设计一
个3dB截频为p的数字滤波器,并与脉冲响应不变法
IIR数字滤波器设计的基本思想 模拟低通滤波器设计 模拟域频率变换 脉冲响应不变法 双线性变换法 利用MATLAB设计IIR DF
双线性变换法
问题的提出 双线性变换法的基本原理 双线性变换法设计DF的步骤
双线性变换法
ቤተ መጻሕፍቲ ባይዱ
问题的提出
频率
设计模拟
AF到DF
p,s
变换
p,s 滤波器 H(s)
H(z)
的转换
如何将模拟滤波器转变为数字滤波器? 1. 脉冲响应不变法 2. 双线性变换法
双线性变换法
回顾内容
脉冲响应不变法的优缺点 优点:
1、数字滤波器和模拟滤波器的频率关系为线性
=T
2、数字滤波器的单位脉冲响应完全模仿模拟滤波器 的单位脉冲响应,时域特性逼近好。 缺点:存在频谱混叠,故不能用脉冲响应不变法设
p =0.6p
0.6
1
Normalized frequency
例:用双线性变换法和一阶巴特沃思低通滤波器,设计一
个3dB截频为p的数字滤波器,并与脉冲响应不变法 设计的DF比较。
H双(z)和H脉(z)幅度响应比较的MATLAB实现
Wp=0.6*pi; b=[1-exp(-Wp)];b1=tan(Wp/2)*[1 1]; a=[1 -exp(-Wp)];a1=[1+tan(Wp/2) tan(Wp/2)-1]; w=linspace(0,pi,512); h=freqz(b,a,w);h1=freqz(b1,a1,w); plot(w/pi,(abs(h)),w/pi,(abs(h1)) ); xlabel('Normalized frequency'); ylabel('Amplitude'); set(gca,'ytick',[0 0.7 1]); set(gca,'xtick',[0 Wp/pi 1]); grid;
例:用双线性变换法和一阶巴特沃思低通滤波器,设计一
个3dB截频为p的数字滤波器,并与脉冲响应不变法
设计的DF比较。
1
脉冲响应不变法
0.7
3dB
0 0
Amplitude
脉冲响应不变法存在频 谱混叠,所设计的DF不满 足给定指标。而双线性变 换法不存在频谱混叠,所 设计的DF满足给定指标。
双线性变换法
个3dB截频为p的数字滤波器,并与脉冲响应不变法
设计的DF比较。
解:设双线性变换中的参数为T
(1) 将DF的频率指标转换为AF的频率指标
p
=
2 T
tan(p
2
)
(2) 设计3dB截频为p的一阶BW型模拟低通滤波器,即
N=1, c = p
故
H(s) = 1 = 1 = s / c 1 s / p 1
因果、稳定的AF系统映射为因果、稳定的DF系统
双线性变换法
双线性变换法的基本原理
双线性变换法的优缺点 优点:无混叠 缺点:频率转换非线性(幅度响应不是常数时 会产生幅度失真)
双线性变换法
双线性变换法设计DF的步骤
1. 将数字滤波器的频率指标{k}转换为
模拟滤波器的频率指标{k}
k
=2 T
例:利用AF-BW filter及脉冲响应不变法和双线性变换法设
计一DLF,满足p=0.2p, s=0.6p, p2dB, s15dB
解:
(1) 将数字低通指标转换成模拟低通指标= / T. p=0.2p/T, s=0.6p/T, p2dB, s15dB
(2) 设计模拟低通滤波器 (BW型)
N
log10 (110000..11ps
1) 1
=2
2 log10 (s / p )
tan(k
2
)
2. 由模拟滤波器的指标设计模拟滤波器的H(s)。
3. 利用双线性变换法,将H(s)转换H(z)。
H(z) = H(s) s=21z1 T 1z1 双线性变换法
双线性变换法设计DF的步骤
= 2 p,s T
tan( ) 2p,s
设计模拟 滤波器
双线性变换
H(s)
H(z)