《数字信号处理》作业程佩青(第2版)清华大学出版社课后答案
- 格式:pdf
- 大小:923.75 KB
- 文档页数:55
数字信号处理(程佩青)课后习题解答(1)1. 什么是数字信号处理?数字信号处理(Digital Signal Processing,DSP)是指对数字信号进行滤波、采样、压缩、编码和解码等操作的一种信号处理技术。
数字信号处理通过离散采样将连续时间信号转换为离散时间信号,并利用数学算法对离散时间信号进行处理和分析。
数字信号处理广泛应用于音频处理、图像处理、视频处理、通信系统等领域。
2. 采样定理的原理是什么?采样定理又称为奈奎斯特-香农采样定理(Nyquist-Shannon Sampling Theorem),是指在进行模拟信号的离散化处理时,采样频率必须大于模拟信号中最高频率的两倍。
采样定理的原理是根据信号的频谱特性,将模拟信号转换为离散时间信号时,需要保证采样频率足够高,以避免采样后的信号出现混叠现象,即频域上的重叠造成的信息损失。
根据奈奎斯特-香农采样定理,采样频率必须大于模拟信号中最高频率的2倍,才能完全还原原始信号。
3. 什么是混叠现象?如何避免混叠现象?混叠现象是指在进行模拟信号的采样时,由于采样频率低于模拟信号中的最高频率,导致频域上的重叠,从而造成采样信号中出现与原始信号不一致的频谱。
混叠现象会使得原始信号的高频部分被错误地表示成低频部分,从而损失了原始信号的信息。
为了避免混叠现象,可以采取以下措施:- 提高采样频率:采样频率必须大于模拟信号中最高频率的两倍,以保证信号的频谱不发生重叠。
- 使用低通滤波器:在采样前,先通过低通滤波器将模拟信号中的高频成分滤除,以避免混叠现象。
滤波器的截止频率应该设置为采样频率的一半。
4. 离散时间信号和连续时间信号有哪些区别?离散时间信号和连续时间信号是两种不同的信号表示形式。
离散时间信号是在时间上离散的,通常由序列表示,每个时间点上有对应的取样值。
离散时间信号可以通过采样连续时间信号得到,采样时将连续时间信号在一定时间间隔内进行取样。
连续时间信号是在时间上连续的,可以用数学函数、图像或者波形图来表示,不存在取样点。
数字信号处理答案第二版答案教程答案姓名:网课答案查询年级:V芯恭zhong號椰子答题分数:100.0问:“73855定律”是由谁提出的?答:艾伯特·梅拉比安问:“落后就要挨打”,这句话是谁说的?()答:斯大林问:武装力量动员,通常是指现役部队动员,是武装力量动员的首要对象。
答:√问:在信息时代,核心技术的一级层次是( )。
答:C问:新军事变革是军事发展的必然结果。
答:正确问:创业计划书的读者人群包括:答:员工管理者合作伙伴创业团队问:“水部火灾,金司空大兴土木;南人北相,中书科甚么东西”一联使用了()的技巧。
答:借义问:网络攻击可分为主动攻击和被动攻击,下面不属于主动攻击的是()。
答:截获问:费用效果分析与费用效益分析的差异,下列说法正确的是( )。
答:费用效益分析单位统一,认可度高,结果易于被人们接受;费用效果分析回避了效果定价的难题,最适于效果难于货币化的领域;费用效益分析与费用效果分析使用领域存在差异;费用效益分析与费用效果分析的基本原则是相同的,即最大限度地节约稀缺资源,最大程度地提高经济效果。
问:下列不属于第三周期元素的是?答:铍问:魏良辅究竟对昆山腔做了哪些改造?答:调理腔调和语音的关系完善和提升曲调的音乐性兼容并蓄融合南北曲为一炉伴奏场面和乐队编制的完善问:进行变异数分析时,母体是否须符合常态分配?答:是问:创新按照成果的性质的不同可分为不同类型,电视机的发明属于()。
答:改进创新问:对大学的学习氛围感到不习惯,怎么办?()答:以上都是问:中医认为情志太过会伤及脏腑,一般认为悲、忧伤()。
答:肺问:当经济发展不成熟时,必须由政府进行协调。
()答:正确问:以下哪一个不是我国的诉讼程序法?答:海事诉讼法问:失血量达到20%就会危及病人生命。
()答:错误问:粤菜烩古今中外烹饪技术于一炉,且以()为主答:海味问:葡萄白兰地按酿造原料不同可分为()答:葡萄原汁白兰地葡萄皮渣白兰地葡萄酒泥白兰地问:小明是个动手能力强,做事手脚灵活,动作协调的人。
第一章习题与解答1 .直接计算下面两个序列的卷积和)n (h *)n (x )n (y =,请用公式表示。
⎩⎨⎧<≤=⎩⎨⎧-≤≤=-00 0, 010 , 0n n ,n n ,)n (x n N n a )n (h n n n β其他 2 .已知线性移不变系统的输入为)n (x ,系统的单位抽样响应为)n (h ,试求系统的∑∞-∞=-==m m n h m x n h n x n y )()()(*)()(:解0)()1(0=<n y n n 时当, 1)2(00部分重叠时当-+≤≤N n n n ()∑∑∑==--===-=nn m mnnn n m mn n m nn m m n h m x n y 0)()()(αββααβ()()βαβαβαβααβαβαβ≠--=--=-+-++-,10111n n n n n n n n())(,1)(00βαα=-+=-n n n y n n ()∑∑∑+-=+-=--+===-=-+≥nN n m m nn nN n m mn n m nn m m n h m x n y N n n 111N -00)()()( , 1)3(αββααβ全重叠时当()()()()βααβαβαβαββααβαβαβ==≠--=--=---+++--,)(,100111n n N N n N n n N n n nN n y输出)n (y ,并画图。
)(5.0)(,)1(2 )()4()(5.0)(,)2( )()3()()(,)( )()2()()(,)( )()1(3435n u n h n u n x n R n h n n x n R n h n R n x n R n h n n x n n n =--==-=====δδ3 .已知 10,)1()(<<--=-a n u a n h n ,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
第一章 离散时间信号与系统1 .直接计算下面两个序列的卷积和)n (h *)n (x )n (y =请用公式表示。
分析:①注意卷积和公式中求和式中是哑变量m ( n 看作参量), 结果)(n y 中变量是 n ,; )()()()()(∑∑∞-∞=∞-∞=-=-=m m m n x m h m n h m x n y ②分为四步 (1)翻褶( -m ),(2)移位( n ),(3)相乘,; )( )( 4n y n n y n 值的,如此可求得所有值的)相加,求得一个(③ 围的不同的不同时间段上求和范一定要注意某些题中在 n如此题所示,因而要分段求解。
2 .已知线性移不变系统的输入为)n (x ,系统的单位抽样响应 为)n (h ,试求系统的输出)n (y ,并画图。
)(5.0)(,)1(2 )()4()(5.0)(,)2( )()3()()(,)( )()2()()(,)( )()1(3435n u n h n u n x n R n h n n x n R n h n R n x n R n h n n x n n n =--==-=====δδ分析:①如果是因果序列)(n y 可表示成)(n y ={)0(y ,)1(y ,)2(y ……},例如小题(2)为)(n y ={1,2,3,3,2,1} ;②)()(*)( , )()(*)(m n x n x m n n x n x n -=-=δδ ;0 00 , 01()0 , ,()0,n n n a n N h n n n n x n n n β-⎧≤≤-=⎨⎩⎧≤⎪=⎨<⎪⎩其他()∑∑∑+-=+-=--+===-=-+≥nN n m mn n nN n m mn n m nn m m n h m x n y N n n 111N -00)()()( , 1)3(αββααβ全重叠时当()()()()βααβαβαβαββααβαβαβ==≠--=--=---+++--,)(,100111n n N N n N n n N n n n N n y③卷积和求解时,n 的分段处理。
1.用冲激响应不变法将以下 )(s H a 变换为 )(z H ,抽样周期为T。
为任意正整数 ,)()( )2()()( )1(022n s s As H b a s a s s H na a -=+++=分析:①冲激响应不变法满足)()()(nT h t h n h a nT t a ===,T 为抽样间隔。
这种变换法必须)(s H a 先用部分分式展开。
②第(2)小题要复习拉普拉斯变换公式1!][+=n n S n t L ,n a n t s a S S As H t u n t Ae t h )()()()!1()(010-=⇔-=-,可求出)()()(kT Th t Th k h a kT t a ===,又dz z dX zk kx )()(-⇔,则可递推求解。
解: (1)22111()()2a s a H s s a b s a jb s a jb ⎡⎤+==+⎢⎥+++++-⎣⎦[])( 21)()()(t u e et h t jb a t jb a a --+-+= 由冲激响应不变法可得:[]()()()() ()2a jb nTa jb nT a T h n Th nT ee u n -+--==+ 11011() () 211n aT jbT aT jbT n T H z h n z e e z e e z ∞------=⎡⎤==+⎢⎥--⎣⎦∑ 2211cos 21cos 1 ------+--⋅=ze bT z e bTz e T aT aT aT(2) 先引用拉氏变换的结论[]1!+=n n s n t L可得: na s s As H )()(0-=)()!1()(10t u n t Ae t h n t s a -=-则)()!1()()()(10k u n kT Ae T Tk Th k h n kT s a -⋅==-dzz dX zk kx azk u a ZZk )()( , 11)( 1-−→←-−→←-且按)11()()!1( )()!1( )()(111111000--∞=---∞=----=-==∑∑ze dz d z n AT e z k n T TA z k h z H T s n n k kT s n n k k可得⎪⎪⎩⎪⎪⎨⎧=-=-=•••---,3,2)1(1,1)(111000n z e z e AT n z e AT z H n T s T S n T s ,可以递推求得:2. 已知模拟二阶巴特沃思低通滤波器的归一化系统函数为:2'4142136.111)(ss s H a ++=而3dB 截止频率为50Hz 的模拟滤波器,需将归一化的)('s H a 中的s 变量用502⨯πs来代替424'108696044.928830.444108696044.9)100()(⨯++⨯==s s s H s H a a π 设系统抽样频率为Hz f s 500=,要求从这一低通模拟滤波器设计一个低通数字滤波器,采用阶跃响应不变法。
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列,nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
第四章 快速傅立叶变换运算需要多少时间。
计算需要多少时间,用,问直拉点的,用它来计算每次复加速度为平均每次复乘需如果一台通用计算机的FFT DFT[x(n)]512s 5 s50.1μμ解: 解: ⑴ 直接计算: 复乘所需时间:复加所需时间:⑵用FFT 计算: 复乘所需时间:复加所需时间:运算一次完成。
点试用一个为了提高运算效率值求今需要从值的点实序列是两个已知IFFTN n y n x k Y k X DFT n y n x N k Y k X ,,)(),()(),(,)(),()(),(.2sN T N01152.0512log105log105 2251262261=⨯⨯⨯=⨯⨯=--sT T T sN N T 013824.0 002304.0512log512105.0log105.0 2126262=+=∴=⨯⨯⨯=⨯⨯⨯=--sT T T sN N T 441536.1 130816.0)1512(512105.0)1(105.0 21662=+=∴=-⨯⨯⨯=-⨯⨯⨯=--sN T 31072.1 512105 105 26261=⨯⨯=⨯⨯=--值的过程。
)(),(完成计算点)可用一次()()(综上所述,构造序列)()()()(可得:)()()(再根据都是实序列,)(),(由原题可知:)()()()(()()(性质:又根据可得序列点作对取序列依据题意解 ]Im[ ]Re[ ][][ ][ ).()( )()()( )()();()( ::n y n x IFFT N k jY k X k Z n z n y n z n x n jy n x n z n y n x n jy n x k Y jIDFT k X IDFT k jY k X IDFT DFT n z IFFT N k Z k jY k X k Z k Y n y k X n x +===+=+=+=++=⇔⇔。
输出倒位序顺序频率抽取采用输入自然输出自然数顺序序时间抽取采用输入倒位流图抽取法的按时间抽取法及按频率画出基时), ,,( 2,16.3FFT N -=。
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列,nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2x(m)()h n m -n1 1 1 0 0 0 0 y(n) 0 11 1 1 12 2 1 1 13 3 1 1 1 1 34 0 1 1 1 1 2 511111(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列, nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
第三章 离散傅立叶变换1.如下图,序列x(n)是周期为6的周期性序列,试求其傅立叶级数的系数。
∑∑=-===56265)(~)(~)(X ~:n nkj nkn e n x W n x k π解kj k j k j kj kj e e e e e 562462362262621068101214πππππ-----+++++=计算求得:。
339)5(~; 33)4(~ ; 0)3(~; 33)2(~;339)1(~;60)0(~j X j X X j X j X X +=-==+=-==。
并作图表示试求设)(~),(~)(~ .))(()(~),()(.264k X n x k X n x n x n R n x == ∑∑=-===56265)(~)(~)(~:n nk jnk n en x W n x k X π解k j kj k j e e e πππ---+++=3231。
计算求得: 3)5(~; 1)4(~ ; 0)3(~ ;1)2(~; 3)1(~ ; 4)0(~j X X X X j X X ====-==。
的周期卷积并作图与试求令其它,设 )(~)(~,))(()(~,))(()(~,)2()(,040,1)(.3464n h n x n h n h n x n x n R n h nn n n x ==-=⎩⎨⎧≤≤+=解:在一个周期内的计算值)(~)(~*)(~)(~m n h n h n x n y -==)(~)(~*)(~)(~m n h n h n x n y -==等各序列。
试画出所示如图已知)())((),())3((,))(()())((),())((,))((,13)(.47755633665n R n x n R n x n x n R n x n R n x n x P n x ----)()()5()(x(n)(4)N n 0 ),n -(n )()3()()()2()()(cos )()1()(52000n R n n x n nR n x n R a n x n R n a n x DFT N N N N n N ==<<===δω闭合形式表达式点试求以下有限长序列的])21sin()2sin()21sin()2sin([21])()()()([21)(]1111[)(][)(])([)()(cos )()()(cos )(:0)2(21020)2(2102)2(21)2(21)2(21222)2(21)2(21)2(21222)()(211)(10)(2110211000000000000000000002002002022002ϖπϖϖπϖωωϖπϖϖπϖϖπϖπϖπϖϖϖϖπϖπϖπϖϖϖωωωωωωωωππππππ-++⋅=--+--=--+--=+=+===---+---------+-++-----+---=---=+--=---=-∑∑∑∑k N e N e k N e N e a e ee e e e eeeee e a k R ee ee a k R eea k R e e e a k R en a k X n R n a n x k N j N j k Nj Njk Nj k Nj k Nj NjNjN jk Nj k N j k Nj NjNjNjN k j N j k j N j N N n nj N n nk j N N n nkj n j n j N n N nkj N N N N N N N 解)(111121)(21)()(21)()(cos )( )()(cos )( ) 1 (:)2()2(10)2(10)2(1020010200000k R e ee e a k R e e a k R e e e a k Re n a k X n R n a n x N k N j N j k Nj Nj N N n nN j N n nk N j N N n nk N j n j n j N n N nk N j N ⎥⎥⎦⎤⎢⎢⎣⎡--+--=⎥⎦⎤⎢⎣⎡+=⎥⎦⎤⎢⎣⎡+===--+---=---=+--=---=-∑∑∑∑ωπωωπωωπωππωωπωω解⎥⎥⎥⎦⎤---⎢⎢⎢⎣⎡--=------+-++---)()()()(21)2(21)2(21)2(21222)2(21)2(21)2(212220000000ωπωπωπωωωωπωπωπωωωk Nj k N j k N j N jN jNjk N j k N j k N j Nj N j N j e e e eeee e e e e e a⎥⎥⎥⎥⎦⎤--⎢⎢⎢⎢⎣⎡+⋅=--+--)21sin()2sin()21sin()2sin(210)2(21020)2(21020000ωπωωπωωπωωπωk N e Nek N e Ne a k Nj Njk N j N jk Nj N N n nk NjnN n aea eak X n R a n x ππ210211)()()((2)--=---===∑)( )()( )()()( 0,)()( (3)02102010200k Re k Re n n k Re n x k X N n n n n x Nk n N j NN n nk N j N n NnkN j πππδδ--=--=-=-==<<-=∑∑)(1)( 11)1()())1(()(])1)2( 2[)1( 32()1)(()()()()( )()(411)1(32)1(321)1(110)1(1k R W Nk X N W W N k R W N k R N W N W W W N W W W nW nWW k X k R nW k X W k R nW k X n nR n x N kNkNkN N N n nk N N k N N k N k N k N N kN k N k N N n kn N N n nk Nk NN n N k n N k NN n N nkN N --=∴-=--+--=+--=-+-+++--++++=-=-==∴=∑∑∑∑∑-=---=+-=-=+-=∙∙∙∙∙∙)(kNN N n nkNN W Nk X n nR n x W n k X n R n n x --===∴=∑-=1)()()()4()( )()(5111022,则小题的结论根据第)(221111122)1(232)1(23210)1(2121)1(2)1()2()(12)2()(2)2(2)2()12()1(]1()2(4[)1(94)1)(()(k N kN kNN n nk NN n nkNk N N kN k N k N N k N k N k N N n kn NN n nk NkNN n kn N k NW N W N N k X W NN N k X N N nWN N W n N N W N W W W N W W W W n W nW k X W n k X W ---=∴----=+--=+--=-+--=-+-+++--++++=-=-=∑∑∑∑∑-=-=---=+-=-=+∙∙∙∙∙∙)∙∙∙±±±===∑-=,6,4,2,0)(~)3(?])0([)()2(?)()1(:;)(~1)(~).(~.61)/2(k k x X k X k X ek X Nn x n x N k nkN j 哪些序列列能做到成虚数外除时间原点使所有的哪些序列能够通过选择成为实数时间原点使所有的哪些序列能够通过选择问傅里叶级数这些序列可以表示成列如图画出了几个周期序π条件。
第七章 有限长单位冲激响应(FIR )数字滤波器的设计方法1. 用矩形窗设计一个FIR 线性相位低通数字滤波器。
已知 21,5.0==N c πω。
求出)(n h 并画出)(log 20ωj e H 曲线。
分析:此题给定的是理想线性相位低通滤波器,故⎪⎩⎪⎨⎧<<<<≤≤=-。
-- , , 0- , )(c c c c ωωππωωωωωωαωj j d e eH解:ωπππωωd eeH n h nj j d d ⎰-=)(21)()()](sin[21αωαωπωωπωωωωα--==⎰--n n d eec c c nj j cc⎪⎪⎩⎪⎪⎨⎧≤≤--====-=为其他故:其中n n n n n w n h n h N d c ,0200,)10(]2sin[)()()(5.0 102/)1( πππωα h( 0)= 9.7654073033E-4h( 1)= 3.5358760506E-2 h( 2)= -9.7657600418E-4 h( 3)= -4.5465879142E-2 h( 4)= 9.7651791293E-4 h( 5)= 6.3656955957E-2 h( 6)= -9.7658322193E-4 h( 7)= -1.0610036552E-1 h( 8)= 9.7643269692E-4 h( 9)= 3.1830877066E-1 h( 10)= 4.9902343750E-1 h( 11)= 3.1830900908E-1 h( 12)= 9.7669276875E-4 h( 13)= -1.0610023141E-1 h( 14)= -9.7654142883E-4 h( 15)= 6.3657015562E-2 h( 16)= 9.7660662141E-4 h( 17)= -4.5465819538E-2 h( 18)= -9.7654841375E-4 h( 19)= 3.5358794034E-2 h( 20)= 9.7658403683E-42.用三角形窗设计一个FIR 线性相位低通数字滤波器。
1.用冲激响应不变法将以下 )(s H a 变换为 )(z H ,抽样周期为T。
为任意正整数 ,)()( )2()()( )1(022n s s As H b a s a s s H na a -=+++=分析:①冲激响应不变法满足)()()(nT h t h n h a nT t a ===,T 为抽样间隔。
这种变换法必须)(s H a 先用部分分式展开。
②第(2)小题要复习拉普拉斯变换公式1!][+=n n S n t L ,n a n t s a S S As H t u n t Ae t h )()()()!1()(010-=⇔-=-,可求出)()()(kT Th t Th k h a kT t a ===,又dz z dX zk kx )()(-⇔,则可递推求解。
解: (1)22111()()2a s a H s s a b s a jb s a jb ⎡⎤+==+⎢⎥+++++-⎣⎦[])( 21)()()(t u e e t h tjb a t jb a a --+-+=由冲激响应不变法可得:[]()()()() ()2a jb nTa jb nT a T h n Th nT ee u n -+--==+ 11011() () 211n aT jbT aT jbT n T H z h n z e e z e e z ∞------=⎡⎤==+⎢⎥--⎣⎦∑2211cos 21cos 1 ------+--⋅=ze bT z e bTz e T aT aT aT(2) 先引用拉氏变换的结论[]1!+=n n sn t L可得: na s s As H )()(0-=)()!1()(10t u n t Ae t h n t s a -=-则)()!1()()()(10k u n kT Ae T Tk Th k h n kT s a -⋅==-dzz dX zk kx azk u a ZZk )()( , 11)( 1-−→←-−→←-且按)11()()!1( )()!1( )()(111111000--∞=---∞=----=-==∑∑ze dz d z n AT e z k n T TA z k h z H T s n n k kT s n n k k可得⎪⎪⎩⎪⎪⎨⎧=-=-=•••---,3,2)1(1,1)(111000n z e z e AT n z e AT z H n T s T S n T s ,可以递推求得:2. 已知模拟二阶巴特沃思低通滤波器的归一化系统函数为:2'4142136.111)(ss s H a ++=而3dB 截止频率为50Hz 的模拟滤波器,需将归一化的)('s H a 中的s 变量用502⨯πs来代替424'108696044.928830.444108696044.9)100()(⨯++⨯==s s s H s H a a π 设系统抽样频率为Hz f s 500=,要求从这一低通模拟滤波器设计一个低通数字滤波器,采用阶跃响应不变法。