第7章 两道滤波器组
- 格式:doc
- 大小:534.50 KB
- 文档页数:11
第6章 有限脉冲响应(FIR)数字滤波器的设计习题1. 已知FIR 滤波器的单位脉冲响应为:(1) h (n )长度N =6h (0)=h (5)=1.5h (1)=h (4)=2h (2)=h (3)=3(2) h (n )长度N =7h (0)=- h (6)=3h (1)=- h (5)=- 2h (2)=-h (4)=1h (3)=0试分别说明它们的幅度特性和相位特性各有什么特点。
2. 已知第一类线性相位FIR 滤波器的单位脉冲响应长度为16, 其16个频域幅度采样值中的前9个为:H g (0)=12, H g (1)=8.34, H g (2)=3.79, H g (3)~H g (8)=0根据第一类线性相位FIR 滤波器幅度特性H g (ω)的特点, 求其余7个频域幅度采样值。
3. 设FIR 滤波器的系统函数为求出该滤波器的单位脉冲响应h (n ), 判断是否具有线性相位, 求出其幅度特性函数和相位特性函数。
4. 用矩形窗设计线性相位低通FIR 滤波器, 要求过渡带宽度不超过π/8 rad 。
希望逼近的理想低通滤波器频率响应函数H d (e j ω)为(1) 求出理想低通滤波器的单位脉冲响应h d (n );(2) 求出加矩形窗设计的低通FIR 滤波器的单位脉冲响应h (n )表达式, 确定)9.01.29.01(101)(4321−−−−++++=z z z z z Hα与N之间的关系;(3)简述N取奇数或偶数对滤波特性的影响。
5.用矩形窗设计一线性相位高通滤波器,要求过渡带宽度不超过π/10 rad。
希望逼近的理想高通滤波器频率响应函数H d(e jω)为(2)用h1(n)和h2(n)分别构成的低通滤波器是否具有线性相位?群延时为多少?题8图9.对下面的每一种滤波器指标,选择满足FIRDF设计要求的窗函数类型和长度。
(1)阻带衰减为20 dB,过渡带宽度为1 kHz,采样频率为12 kHz;(2)阻带衰减为50 dB,过渡带宽度为2 kHz,采样频率为20 kHz;(3)阻带衰减为50 dB,过渡带宽度为500 Hz,采样频率为5 kHz。
第7章 傅里叶变换与滤波器形状 7.1离散时间傅里叶变换基础离散时间傅里叶变换(DTFT )是数字信号分析的一个重要工具。
DTFT 把信号或滤波器从时域变换到频域,主要是为了研究信号或滤波器的频率特性。
该变换主要用于分析信号和滤波器的频谱性质。
对于信号,DTFT 提供的信息称为信号的频谱。
对于滤波系统,DTFT 得到的信息称为滤波器的频率响应(frequency response )。
它由两部分组成:幅度响应(magnitude response )和相位响应(phase response )。
幅度响应给出了滤波器的形状,通过它我们可以深入了解滤波器的工作特性。
信号x[n]的离散时间傅里叶变换定义为:()[]jn n X x n e∞-Ω=-∞Ω=∑,这里Ω为数字频率,单位弧度。
记为(){[]}X x n Ω=F利用欧拉公式,DTFT 变换为()[][](cos()sin())jn n n X x n ex n n j n ∞∞-Ω=-∞=-∞Ω==Ω-Ω∑∑变换()X Ω在每个不同的数字频率上可有不同的值,当信号x[n]与正弦或余弦“共振”时,最大。
也就是说,当x[n]以接近频率Ω变化时,()X Ω较大。
离散时间傅里叶变换反应了信号的频率。
例7.1 求如图信号的离散时间傅里叶变换注意,一般情况,DTFT 是复值。
例7.2 求信号x[n]=4(u[n]-u[n-3])的DTFT 。
离散时间傅里叶变换有两个重要的特性,时延特性和周期性。
00[]()[]()Fjn Fx n X x n n eX -Ω−−→Ω-−−→Ω(2)()X X πΩ+=ΩDTFT 是周期性的,周期为2π。
即离散时间傅里叶变换对所有的数字频率Ω,每2π重复一次,不断重复。
7.2.1 频率响应和差分方程 对差分方程逐项求DTFT∑∑==-=-Mk kN k kk n x b k n y a 0][][][]2[]1[][][]2[]1[][210210M n x b n x b n x b n x b N n y a n y a n y a n y a M N -++-+-+=-++-+-+0101()()()()()()j jN N j jM M a Y a e Y a e Y b X b eX b eX -Ω-Ω-Ω-ΩΩ+Ω++Ω=Ω+Ω++Ω0101()()()()j jN j jM N M a a e a e Y b b e b e X -Ω-Ω-Ω-Ω+++Ω=+++Ω010010()()()Mjk j jM k M k Nj jN jk N k k b eY b b e b eH X a a e a ea e-Ω-Ω-Ω=-Ω-Ω-Ω=Ω+++Ω===Ω+++∑∑例7.3 求差分方程频率响应y[n]=-0.85y[n-1]+0.5x[n].例7.4 求差分方程频率响应y[n]+0.1y[n-1]+0.85y[n-2]=x[n]-0.3x[n-1]7.2.2 频率响应和传输函数10101010()()()MkMkM k NNkN k k b zb b z b zY z H z X z a a z a za z---=---=+++===+++∑∑010010()()()Mjk j jM k M k Nj jN jk N k k b eY b b e b eH X a a e a e a e-Ω-Ω-Ω=-Ω-Ω-Ω=Ω+++Ω===Ω+++∑∑例 7.5 求滤波器的频率响应,它的传输函数21210.2()10.50.9z H z z z ----=++2210.2()10.50.9j j j e H e e-Ω-Ω-Ω-Ω=++频率响应是脉冲响应的DTFT 。
第7章 模拟信号的数字传输7.1 学习指导 7.1.1 要点本章的要点主要有抽样定理;自然抽样和平顶抽样;均匀量化和非均匀量化;PCM 原理,A 律13折线编码,译码;ΔM 原理,不过载条件;PCM ,ΔM 系统的抗噪声性能;PCM 与ΔM 的比较;时分复用和多路数字电话系统原理;1. 概述为了使模拟信号实现数字化传输,首先要通过信源编码使模拟信号转换为数字信号,或称为“模/数转换”即A/D 转换。
模/数转换的方法采用得最早而且应用较广泛的是脉冲编码调制(PCM),PCM 通信系统原理图如图7-1所示。
图7-1 PCM 通信系统原理图抽样量化器编码器模拟信号PCM 信号译码器低通滤波器模拟信号数字通信系统PCM 信号由图7-1可见,PCM 系统由以下三部分组成。
(1) 模/数转换(A/D 转换)模/数转换包括三个步骤:抽样(Sampling)、量化(Quantization)和编码(Coding)。
a. 抽样是把在时间上连续的模拟信号转换成时间上离散的抽样信号,抽样信号在时间上是离散的,但是其取值仍然是连续的,所以是离散模拟信号。
b. 量化。
量化是把幅度上连续的抽样信号转换成幅度离散的量化信号,故量化信号已经是数字信号了,它可以看成是多进制的数字脉冲信号。
c. 是编码。
编码是把时间离散且幅度离散的量化信号用一个二进制码组表示。
(2) 数字方式传输——基带传输或带通传输;(3) 数/模转换(D/A )——将数字信号还原为模拟信号。
包含了译码器和低通滤波器两部分。
2.抽样定理为模拟信号的数字化和时分多路复用(TDM )奠定了理论基础。
根据抽样的脉冲序列是冲激序列还是非冲激序列,抽样可以分为理想抽样和实际抽样。
抽样是按照一定的抽样速率,把时间上连续的模拟信号变成一系列时间上离散的抽样值的过程。
能否由此样值序列重建源信号,取决于抽样速率大小,而描述这一抽样速率条件的定理就是著名的抽样定理。
(1) 低通信号的抽样定理定理:设有一个频带限制在(0,f H )内的连续模拟信号m (t ),若以T s ≤1/(2f H )间隔对它抽样,则m (t )将被这些抽样值所完全确定。
第7章滤波器的设计方法教学目的1.掌握由连续时间滤波器设计离散时间IIR滤波器的方法,包括冲激响应不变法,双线性变换法等;2.了解常用的窗函数,掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法;3.掌握FIR滤波器的逼近原理与设计方法。
教学重点与难点重点:本章是本课程的重中之重,滤波器的设计是核心内容之一。
1.连续时间滤波器设计离散时间IIR滤波器的方法,包括冲激响应不变法,双线性变换法等;2.常用的窗函数,掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法;3.掌握FIR滤波器的逼近原理与设计方法。
难点:1.冲激响应不变法,双线性变换法2.用窗函数法设计FIR滤波器FIR滤波器的逼近原理与设计方法基本概念7.0.1 选频滤波器的分类数字滤波器是数字信号处理的重要基础。
在对信号的过滤、检测与参数的估计等处理中, 数字滤波器是使用最广泛的线性系统。
数字滤波器是对数字信号实现滤波的线性时不变系统。
它将输入的数字序列通过特定运算转变为输出的数字序列。
因此,数字滤波器本质上是一台完成特定运算的数字计算机。
我们已经知道,一个输入序列x(n),通过一个单位脉冲响应为h(n)的线性时不变系统后,其输出响应y(n)为∑∞-)(y))()()(n(nn=m*=xmhnhx将上式两边经过傅里叶变换,可得式中,Y (e j ω)、X (e j ω)分别为输出序列和输入序列的频谱函数, H (ejω)是系统的频率响应函数。
可以看出,输入序列的频谱X (e j ω)经过滤波后,变为X (e j ω)H (e j ω)。
如果|H (e j ω)|的值在某些频率上是比较小的,则输入信号中的这些频率分量在输出信号中将被抑制掉。
因此,只要按照输入信号频谱的特点和处理信号的目的,适当选择H (ej ω),使得滤波后的X (e j ω)H (e j ω)符合人们的要求,这就是数字滤波器的滤波原理。
和模拟滤波器一样,线性数字滤波器按照频率响应的通带特性可划分为低通、高通、带通和带阻几种形式。
第7章信号的运算和处理自测题一、现有电路:A.反相比例运算电路B.同相比例运算电路C.积分运算电路D.微分运算电路E.加法运算电路F.乘方运算电路选择一个合适的答案填入空内。
(1)欲将正弦波电压移相+90o,应选用( C )。
(2)欲将正弦波电压转换成二倍频电压,应选用( F )。
(3)欲将正弦波电压叠加上一个直流量,应选用( E )。
(4)欲实现A u=−100 的放大电路,应选用( A )。
(5)欲将方波电压转换成三角波电压,应选用( C )。
(6)欲将方波电压转换成尖顶波波电压,应选用( D )。
二、填空:(1)为了避免50H z电网电压的干扰进入放大器,应选用( 带阻)滤波电路。
(2)已知输入信号的频率为10kH z~12kH z,为了防止干扰信号的混入,应选用( 带通)滤波电路(3)为了获得输入电压中的低频信号,应选用( 低通)滤波电路。
(4)为了使滤波电路的输出电阻足够小,保证负载电阻变化时滤波特性不变,应选用( 有源)滤波电路。
三、已知图T7.3所示各电路中的集成运放均为理想运放,模拟乘法器的乘积系数k大于零。
试分别求解各电路的运算关系。
(a)(b)图T7.3解:图(a)所示电路为求和运算电路,图(b)所示电路为开方运算电路。
它们的运算表达式分别为:(a) 12413121234()(1)//f I I O f I R u u R u R u R R R R R R =-+++⋅⋅+ 11O O u u dt RC =-⎰(b) '23322144O I O O R R R u u u ku R R R =-⋅=-⋅=-⋅ 2413O I R R u u kR R =⋅习题本章习题中的集成运放均为理想运放。
7.1填空:(1) ( 同相比例 )运算电路可实现A u >1 的放大器。
(2) ( 反相比例 )运算电路可实现A u <0 的放大器。
(3) ( 微分 )运算电路可将三角波电压转换成方波电压。
第7章 两通道滤波器组本本章章首首先先讨讨论论图图77..11..11中中各各信信号号间间的的一一般般关关系系,,然然后后讨讨论论实实现现准准确确重重建建的的途途径径。
也也即即,,如如何何确确定定)(0z H ,,)(1z H ,,)(0z G 和和)(1z G 才才能能去去除除混混叠叠失失真真,,幅幅度度失失真真及及相相位位失失真真,,即即P P R R 。
以下结论设用于F F I I R R //I I I I R R 系统:)()(z H z F -=,则()()()j j F e H e ωωπ±=,()(1)()n f n h n =- (7.2.17a )1()()F z H z -=,则*()()j j F e H e ωω=(实系数),()()f n h n =- (7.2.17b )1()()F z H z -=-,则*()()()j j F e H e ωωπ±=(实系数),()(1)()n f n h n =--(7.2.17c )(1)1()()N F z z H z ---=-,则(1)*()()()j j N j F e e H e ωωωπ--±=(实系数),1()(1)(1)N n f n h N n --=--- (7.2.17d )7.1 两通道滤波器组中各信号的关系图7.1.1两通道滤波器组经推导得到:[]⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡--=)()()()()()()()(21)(ˆ110010z X z X z H z H z H z H z G z G z X(7.1.4))(z X -是是混混叠叠分分量量((由由抽抽取取和和插插值值引引起起的的))。
)由此得:两通道滤波器组的调制性表示:定义)]()()()([21)()]()()()([21)(1100111000z G z H z G z H z F z G z H z G z H z F -+-=+=※(7.1.5)※则)()()()()(ˆ10z X z F z X z F z X -+= (7.1.6)若0)(1=z F ,则可有效的去除混叠失真,这样:)()]()()()([21)()()(ˆ11000z X z G z H z G z H z X z F z X +== ※(7.1.7)※)(0z F 反反应应了了去去除除混混叠叠失失真真后后的的两两通通道道滤滤波波器器组组的的总总的的传传输输特特性性。
系系统统的的幅幅度度失失真真及及相相位位失失真真均均与与)(0z F 有有关关,,因因此此)(0z F 又又称称““失失真真传传递递函函数数((d d i i s s t t o o r r t t i i o o n n t t r r a a n n s s f f e e r r f f u u n n c c t t i i o o n n ))””。
7.2 准确重建的充要条件和)(0z G 和)(1z G 的选择定义: )()()(10z H z H z P -= ……(7.2.3a) 则)()()(10z H z H z P -=-……(7.2.3b)再定义m H z P z P z T det )()()(=--=……(7.2.4))(z P 可以刻划)()()(z P z P z T --=的内部结构。
A) 首先,我们先深入的讨论一下)(0z G 和和)(1z G 的的一一般般选选取取原原则则:由(7.1.4)式可得:(21z X (21z X -⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡-)()()()()()()()()()(21)(ˆ)(ˆ11001010z X z X z H z H z H z H z G z G z G z G z X z X……(7.2.7) 令⎥⎦⎤⎢⎣⎡--=)()()()(1010z G z G z G z G m G , ⎥⎦⎤⎢⎣⎡--=)()()()(1010z H z H z H z H m H[]Tz X z X )(ˆ)(ˆˆ-=X ,[]Tz X z X )()(-=X则X H G XT m m 21ˆ=……(7.2.10) m H 又称“混迭分量(AC )”矩阵。
Xˆ对X 准确重建的【充要条件1】:其中含有)(0z H ,,)(1z H ,,)(0z G 和)(1z G⎥⎦⎤⎢⎣⎡-=--k k Tm m z z )(0021H G (7.2.11)经推导得到【充要条件2】:同样,其中含有)(0z H ,,)(1z H ,,)(0z G 和)(1z G⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-=-=-==--)()(det 2)()()(det 2)(0110z H z T H cz z G z H z T H cz z G m kmk(7.2.14)其中,m H z P z P z T det )()()(=--=,也称为“失真函数”可见,在两通道滤波器组中,)(0z G 、)(1z G 选取是独立于)(0z H 及)(1z H 的,也即,对于给定的0H 和1H ,只要按(7.2.14)式选定)(0z G 和)(1z G ,即可实现PR 。
分析该式,我们发现:1.取)(0z H 和)(1z H 是FIR 的,则)(z T 是FIR 的,但由于)(z T 出现在(7.2.14)式的分母上,那么)(0z G 和)(1z G 将变成IIR 的,这是在应用中所不希望的; 2.保证)(0z G 和)(1z G 是FIR 的,唯一可能而又“方便”(就用c,不用c ’)的做法是令l k cz z T +-=2)(,(7.2.15)则⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧--=-=--)()()()(0110z H z z G z H z z G l l(7.2.16)这也就是【“充要”条件3】,由(7.2.15)((含含有有)(0z H ,,)(1z H ))和(7.2.16)((含含有有)(0z H ,,)(1z H ,,)(0z G 和和)(1z G ))共共同同组组成成。
代入)()]()()()([21)()()(ˆ11000z X z G z H z G z H z X z F z X +==※(7.1.7)※得 )()()(21)(ˆz X cz z X z T z z Xk l --==,正是准确重建。
我们发现,(7.2.16)是去去除除混混叠叠失失真真的的充充分分条条件件。
)(0z G 和)(1z G 最直观、常用的选取方法是l =0:⎭⎬⎫⎩⎨⎧--=-=)()()()(0110z H z G z H z G现在的问题是如何选取)(0z H ,)(1z H ,实现lk cz z T +-=2)((7.2.15)。
B) 现在我们来分析一下)(0z H ,)(1z H ,)(0z G 和)(1z G 在 (7.2.16)式制约下所具有的时域关系(当然前提是都为FIR 滤波器,可复):1) 由(7.2.16)式,有)()1()()()1()(01110l n h n g l n h n g l n l n --=--=+-- (7.2.18)2) 双通道滤波器组中,若按(7.2.16)式选定)(0z G 和)(1z G 。
若l 为奇数,则)(1n g 与)(0n h 的偶次移位正交,)(0n g 和)(1n h 的偶次移位正交;若l 为偶数,则)(1n g 与)(0n h 的奇次移位正交,)(0n g 和)(1n h 的奇次移位正交。
这是两通道滤波器组中的交叉正交关系。
3) 由【“充要”条件3】lk czz P z P z T +-=--=2)()()(,由于()()P z P z --只包含奇次项,所以l k -只能为奇数。
若取0=l ,k 只能是奇数。
这就是说,在0100()()()()()P z H z H z H z G z =-=中,它的奇序号项只能有一项,即在k z -处,于是,我们有li l i c k n h i n g c k n h i n g --=--=--δδ)(),()(),(1100这是反映同一条支路上,分析滤波器和综合滤波器之间的正交关系。
上面这几条时域关系是共性。
后面,共轭正交镜像滤波器组中,将讨论那种情况下的0h 和1h 的关系。
7.4 共轭正交镜像滤波器组)()(10)1(1----=z H z z H N(7.4.1)式中N 为偶数(后面说明原因)。
按【“充要”条件3】,取0=l ,)()()()()(0110)1(10z H z G z H z z H z G N --=-=-=--- (7.4.2)把这四个滤波器时域的关系归纳如下:)1()1()1()1()(01011n N h n N h n h n n N ---=---=+--(7.4.3b))1()(00n N h n g ---= (7.4.3c))()1()(011n h n g n +-=(7.4.3d)—————————————————————————————————————将(7.4.1)式代入(7.2.13)式,有[])()()()(det 100100)1(------+-=z H z H z H z H z N m H(7.4.4)令)()()(100-=z H z H z P C※(7.4.5a)※则)()()()()(111100--=--=-z H z H z H z H z P C ※(7.4.5b)※再令 )()()(z P z P z T C C C -+=※(7.4.6a)※ 则)(det )()1(z T zz T C N m ---==H※(7.4.6b)※请注意此处)(z P C 、)(z T C 的定义与(7.2.3)和(7.2.4)式的区别。
再将CQMFB 的关系代入(7.1.5)式,有)(21)(21)()1(0z T z z T z F C N ---== (7.4.7)为了设计的方便,常取0',')()()(>=-+=c c z P z P z T C C C(7.4.8)则)1(02')(---=N zc z F 为一纯延迟,且倒倒置置。
从而实现)(ˆz X 对)(z X 的准确重建。
在我们利用单带滤波器的方法设计半带滤波器,从而设计出CQMFB 时,实际上取的是1)(=z T C 。
● 另外,这里取c ’=2或1是合理的。
原因有二: (1) 由于)(ωj C eP 恒正,故)(z T C 只能取正的纯延迟形式0',')('>=-c z c z T k C ;(2) 由前节的结论lk cz z P z P z T +-=--=2)()()(,l k -只能为奇数。