数字信号处理 吴镇扬 第二版 第五章习题答案
- 格式:ppt
- 大小:197.50 KB
- 文档页数:12
习题一 (离散信号与系统)1.1周期序列,最小周期长度为5。
1.2 (1) 周期序列,最小周期长度为14。
(2) 周期序列,最小周期长度为56。
1.5()()()()()()()11s a s s s a n s s a s n X j x t p t X j ΩP j Ω2n τn τj sin j Ωjn e X 2n π2n n τj Sa X j jn e 2T 2πττ∞=-∞∞=-∞Ω==*⎡⎤⎣⎦ΩΩ⎛⎫-=-Ω ⎪⎝⎭ΩΩ⎛⎫-=Ω-Ω ⎪⎝⎭∑∑F 1.6 (1) )(ωj e kX (2) )(0ωωj n j e X e (3) )(21)(2122ωωj j e X e X -+ (4) )(2ωj e X1.7 (1)0n z -(2)5.0||,5.0111>--z z(3)5.0||,5.0111<--z z(4)0||,5.01)5.0(11101>----z zz1.8 (1) 0,)11()(211>--=---z z z z z X N (2) a z az az z X >-=--,)1()(211(3)a z az z a az z X >-+=---,)1()(311211.9 1.10(1))1(2)(1----+n u n u n (2))1(24)()5.0(6--⋅--n u n u n n (3))()sin sin cos 1(cos 000n u n n ωωωω++(4) )()()(1n u a a a n a n ---+-δ1.11 (1) )(1z c X - (2) )(2z X (3) )()1(21z X z -+ (4) -+<<x x R z R z X /1/1),/1(1.12 (1)1,11<-ab ab(2) 1 (3) 00n a n1.13 (1) 该系统不是线性系统;该系统是时不变系统。
第五章习题与上机题5.1 已知序列12()(),0 1 , ()()()nx n a u n a x n u n u n N =<<=--,分别求它们的自相关函数,并证明二者都是偶对称的实序列。
解:111()()()()()nn mx n n r m x n x n m a u n au n m ∞∞-=-∞=-∞=-=-∑∑当0m ≥时,122()1mmnx n ma r m aaa∞-===-∑ 当0m <时,122()1m mnx n a r m aaa -∞-===-∑ 所以,12()1mx ar m a =-2 ()()()()N x n u n u n N R n =--=22210121()()()()()1,0 =1,00, =()(1)x NN n n N mn N n m N r m x n x n m Rn R n m N m N m N m m Nm N m R m N ∞∞=-∞=-∞--=-=-=-=-⎧=--<<⎪⎪⎪⎪=-≤<⎨⎪⎪⎪⎪⎩-+-∑∑∑∑其他从1()x r m 和2()x r m 的表达式可以看出二者都是偶对称的实序列。
5.2 设()e()nTx n u n -=,T 为采样间隔。
求()x n 的自相关函数()x r m 。
解:解:()()()()e()e ()nTn m T x n n r m x n x n m u n u n m ∞∞---=-∞=-∞=-=-∑∑用5.1题计算1()x r m 的相同方法可得2e()1e m Tx Tr m --=-5.3 已知12()sin(2)sin(2)s s x n A f nT B f nT ππ=+,其中12,,,A B f f 均为常数。
求()x n 的自相关函数()x r m 。
解:解:()x n 可表为)()()(n v n u n x +=的形式,其中)2sin()(11s nT f A n u π=,=)(n v 22sin(2)s A f nT π,)(),(n v n u 的周期分别为 s T f N 111=,sT f N 221=,()x n 的周期N 则是21,N N 的最小公倍数。
————第一章———— 时域离散信号与系统理论分析基础本章1.1节“学习要点”和1.2节“例题”部分的内容对应教材第一、二章内容。
为了便于归纳总结,我们将《数字信号处理(第二版)》教材中第一章和第二章的内容合并在一起叙述,这样使读者对时域离散线性时不变系统的描述与分析方法建立一个完整的概念,以便在分析和解决问题时,能全面考虑各种有效的途径,选择最好的解决方案。
1.1 学 习 要 点1.1.1 时域离散信号——序列时域离散信号(以下简称序列)是时域离散系统处理的对象,研究时域离散系统离不开序列。
例如,在时域离散线性时不变系统的时域描述中,系统的单位脉冲响应()n h 就是系统对单位脉冲响应()n δ的响应输出序列。
掌握()n δ的时域和频域特征,对分析讨论系统的时域特性描述函数()n h 和频域特性描述函数()ωj e H 和()z H 是必不可少的。
1. 序列的概念在数字信号处理中,一般用()n x 表示时域离散信号(序列)。
()n x 可看作对模拟信号()t x a 的采样,即()()nT x n x a =,也可以看作一组有序的数据集合。
要点 在数字信号处理中,序列()n x 是一个离散函数,n 为整数,如图1.1所示。
当≠n 整数时,()n x 无定义,但不能理解为零。
当()()nT x n x a =时,这一点容易理解。
当=n 整数时,()()nT x n x a =,为()t x a 在nT t =时刻的采样值,非整数T 时刻未采样,而并非为零。
在学习连续信号的采样与恢复时会看到,()n x 经过低通滤波器后,相邻的()T n nT 1~+之间的()t x a 的值就得到恢复。
例如,()n x 为一序列,取()()2n x n y =,n 为整数是不正确的,因为当=n 奇数时,()n y 无定义(无确切的值)。
2. 常用序列常用序列有六种:①单位脉冲序列()n δ,②矩形序列()n R N ,③指数序列()n u a n,④正弦序列()n ωcos 、()n ωsin ,⑤复指数序列nj eω,⑥周期序列。
【最新整理,下载后即可编辑】习题一 (离散信号与系统)1.1周期序列,最小周期长度为5。
1.2 (1) 周期序列,最小周期长度为14。
(2) 周期序列,最小周期长度为56。
1.5()()()()()()()11s a s s s a n s s a s n X j x t p t X j ΩP j Ω2n τn τj sin j Ωjn e X 2n π2n n τj Sa X j jn e 2T 2πττ∞=-∞∞=-∞Ω==*⎡⎤⎣⎦ΩΩ⎛⎫-=-Ω ⎪⎝⎭ΩΩ⎛⎫-=Ω-Ω ⎪⎝⎭∑∑ 1.6 (1) )(ωj e kX (2) )(0ωωj n j e X e (3) )(21)(2122ωωj j e X e X -+ (4) )(2ωj e X1.7 (1) 0n z -(2) 5.0||,5.0111>--z z (3) 5.0||,5.0111<--z z (4)0||,5.01)5.0(11101>----z z z1.8 (1) 0,)11()(211>--=---z zz z z X N(2) a z az az z X >-=--,)1()(211 (3) a z az z a az z X >-+=---,)1()(311211.91.10 (1))1(2)(1----+n u n u n (2))1(24)()5.0(6--⋅--n u n u n n (3))()sin sin cos 1(cos 000n u n n ωωωω++(4) )()()(1n u a a a n a n ---+-δ 1.11(1))(1z c X - (2) )(2z X (3))()1(21z X z -+ (4)-+<<x x R z R z X /1/1),/1(1.12 (1) 1,11<-ab ab(2) 1 (3)00n a n1.13 (1) 该系统不是线性系统;该系统是时不变系统。
第五章 有限长单位脉冲响应滤波器的设计方法 1.解:(a)⎰⎰⎰-------===πωπαωπαππππαπωωωαπωπωπcn d e e d ed ee Hn h n j j j nj j dd )(221)(21)()()(])[()1(])[(]|)(cos |)([sin )(2)()](sin )([cos )(2c a c n c a c j j j n S n S e n j n n e n d n j n n e c c c ωαπωωαπωαωαωαπααωαωαππαπωππωππαπωππα--=-=----=--+--=---⎰ )(])[()1()()()(n n S n n h n h R c a cnR d ωωαπωω--=⋅= (b) 为了保证线性相位)(])[()1()(2112,12,,21n k n S n h k k k N N N R c a c n ωωπωαα--==-+=+=-=为奇数若h(n)的类型取决于n )1(-,N 为奇数h(n)为偶对称第一类,h(n)必须偶对称于n=α 处,否则不满足N 为奇数的已知条件 若N 为偶数。
即N=2k,则)(])[()1()(21212n n S n h k k R c a c n ωωαπωα--=-=-=h(n)必须奇对称于n=α处,否则不满足N 为偶数的已知条件(c))()]12cos 1)(([2)1()(n N n n S n h R c a c n ωπαωπω----= 2.解: (a)⎰⎰+-----+-=ccd e je d e je n h n j j n j j d ωππωαπωπωπωαπωωπωπ)()(2121)(παωααααωααπααπωωπωππωαπωπωαπαωππωαπωπωαπα)(]2/)[(sin )1()cos(2)1(]||[2][221)()()()(---=----=-+--=+-=++---+---⎰⎰n n n n n e n e e d je d je e c n c n n j n j i n j n j j cc c c)()(]2/)[(sin )1()(*)()(21n n n n n h n h R c n R d ϖπαωααϖ---==+(b)为了保证线性相位21-=N α 若N 为奇数,设N=2k+1则α=k)()(]2/)[(sin 2)1()(21n k n k n n h R c n ϖπω---=+h(n)满足奇对称,即h(n)=-h(N-1-n)属于第III 类FIR 滤波器 若N 为偶数,设N=2k 则α=k-1/2πω)2/1(]2/)2/1[(sin 2)1()(21+-+--=+k n k n n h c nh(n)满足偶对称,即h(n)=h(N-1-n)属于第II 类FIR 滤波器(c))(]12cos 1[)(]2/)sin[()1()(1n N nn n n h R c n ϖππαωα-----=+3.解:⎰⎰+----+--+=ccccd e e d e e n h n j j n j j d ωωωωωωαωωωωωωαωπωπ00002121)(παωαωαπαωωαωωα)()sin()cos(2)()])(sin[()])(sin[(000---=----+-=n n n n n n c c c)()()(n n h n h R d ϖ=(a)N 为奇数时,设N=2k+1,k N =-=21α )()()sin()cos(2)(0n k n k n k n n h R cϖπωω---=h(n)满足于偶对称,属于第I 类FIR 滤波器 (b)N 为偶数时,设N=2k,α=k-1/2)()2/1()2/1sin()2/1cos(2)(n k n k n k n n h R cϖπω+-+-+-=h(n)满足偶对称,属于第II 类FIR 滤波器 (c)N 为奇数时,用升余弦窗设计)(]12cos 46.054.0[)()sin()cos(2)(0n R N nk n k n k n n h N c -----=ππωωN 为偶数,用升余弦窗设计)(]12cos 46.054.0[)()2/1sin()2/1cos(2)(0n R N nk n k n k n n h N c ---+-+-=ππωω4.解:与第三题相比知)(ωϕ由-ωα变为-αω-π/2,所以只需将上题)(n h d 由偶对称变为奇对称即可)()()()sgn()()sin()cos(2)(0n n h n h n n n n n h R d cϖαπαωαωα=----=(a)N 为奇数,α=k)()sgn()()sin()cos(2)(0n k n k n k n k n n h R cϖπωω----=奇对称属于第III 类滤波器 (b)N 为偶数,α=k-1/2)()2/1sgn()2/1()2/1sin()2/1cos(2)(0n k n k n k n k n n h R cϖπωω+-+-+-+-=奇对称属于第IV 类滤波器 (c)用改进升余弦窗设计 N 为奇数)sgn()(]12cos 46.054.0[)()sin()cos(2)(0k n n R N nk n k n k n n h N c ------=ππωωN 为偶数)(]12cos 46.054.0[)2/1sgn()2/1()2/1sin()2/1cos(2)(0n R N nk n k n k n k n n h N c--⨯+-+-+-+-=ππωω5.解:(a)一个带阻滤波器相当于一个全通滤波器减去一个带通滤波器 全通)()(ωϕωj j e e H = 带通)()()(ωϕωωj B j B e H e H =则带阻)()()()](1[)()(ωϕωϕωϕωωωj B j B j j r e H e H e e H -=-= (b)因是线性相位滤波器,不妨设ϕ(ω)=-αω⎰--=ππωωϕωωπd e e H n h n j j B r )()](1[21)()()()sin()(21)(n h n n n h d e E E n j ----=-=⎰--παπαωπππωα6.解: (a))()sgn()()sin()()()()sgn()()sin()(sin 21)(2/)(,1|)(|)()2/(n a n a n a n n n h n h a n a n a n n h n n d e n h a e H e je e H R R d d n j d j d a j a j j d ϖππϖππππωππωωϕππωωπωωω---==---=∴==∴--==∴=-=⎰----(b)N 为奇数时,α=(N-1)/2=k⎩⎨⎧≠==∴a n an n h ,0,1)( N 为偶数时,α=(N-1)/2=k-1/2)()2/1sgn()2/1()2/1sin()(n k n k n k n n h R ϖππ+-+-+-=显然N 为偶数时性能好(c))())]1/(1[1()sgn()()sin()(020ββππI N n I a n a n a n n h ------=7.解: (a))()1()()1(21)(n an n h an d e e j n h R an a n nj a d ϖωωπππωω---=---=-=----⎰ (b)N 为奇数时,a=k)()1()(n kn n h R kn ϖ---=-N 为偶数时,a=k-1/2)(2/1)1()(n k n jn h R k n ϖ+---=-N 为奇数时性能好(c) )())]1/(21[1()1()(020ββI N n I k n n h k n ------=-9.解: (a))1(5.0)14(13,,3,2,0)()(5.0)1(,1)0(),/11()21(2)(*15/)115(14)15/11(d j d d d j d d H e H k k N H k H e H H N k N N kk ====-===--=--=----ππππθ ∑-==12)(1)(N n kn N j d e k H N n h π)5.05.01(1511528151********n j j n j j e e e e ππππ++=- ∑=---=1415/22151)(151)(o k k j j j j e ek H e e H ωωω(b)横截型x(n) 1-z 1-z 1-z 1-z)0(h )1(h )2(h )3(h )14(h频率采样型H(0)ωj e -(c)横截型用的乘法器多,频率采样型用的加法器多 10.解: (a))(2n h 为)(1n h 的圆周移位πθθπk k k k H k H e k H W k H k H jk kmN -==∴===)()(|)(||)(|)()()(2121112(b)如图所示,又5.32/)1(,8=-==N a N 知)(1n h ,)(2n h 均关于n=3.5偶对称,所以属于线性相位滤波器时延为3.5 11.解:⎩⎨⎧≤≤=ωωωω其他,00,1|)(|cj d e H 图略第一类线性相位FIR 滤波器,⎪⎪⎩⎪⎪⎨⎧-≤≤+---≤≤--=-≤≤+-=⎪⎪⎩⎪⎪⎨⎧-≤≤+-=≤≤=--=--==121),21)((2210),21(2)(,121|,)(||)(|,211]2int[,041]2int[0,1|)(|)11(2)21()(,e |H(K)|H(K)(K)j N k N N k N N N k N k NK N k N k N H K H N k N N N k K H Nk k N N K c c ππθπωπωππθθ偶对称⎪⎪⎩⎪⎪⎨⎧≤≤-≤≤-=⎪⎩⎪⎨⎧≤≤≤≤≤≤≤≤=3217),33(3332160,3332)(,3225,12417169,080,1|)(|k k k k K k k k k K H ππθ及设计的过渡带宽3322ππω==∆N 如果边沿设定ϖ(k)为一点,即令ϖ(9)= ϖ(24)=0.39∑+++--+==⎪⎩⎪⎨⎧≤≤=≤≤=-]})332sin(33)]332(33sin[)332sin(33)]332(33sin[[2sin33233sin {)(3342*21610,09,39.080,1|)(|16πωπωπωπωωωππωωk k k k e e H N K K K K H j j 则过渡带宽为12.解: (a)N=33⎪⎪⎩⎪⎪⎨⎧≤≤-≤≤-=⎪⎩⎪⎨⎧≤≤=≤≤≤≤=3217),33(3332160,3332)(,2013,121,12,39.03222,110,0|)(|k k k k k k k k k k H ππθ(b)N=32⎪⎪⎩⎪⎪⎨⎧≤≤-≤≤-=⎪⎩⎪⎨⎧≤≤=≤≤≤≤=3318),34(3433170,3433)(,2113,122,12,39.03323,110,0|)(|k k k k k k k k k k H ππθ13.解:(a)N=33,因为N 为奇数,所以可能是第I ,III 型滤波器⎩⎨⎧≤≤≤≤≤≤≤≤≤≤=2821,125,13229,2013,40,0|)(|k k k k k k H 第I 型⎪⎪⎩⎪⎪⎨⎧≤≤-≤≤-=3217),33(3332160,3332)(k k k k k ππθ第III 型⎪⎪⎩⎪⎪⎨⎧≤≤-+-≤≤-=3217),33(33322160,33322)(k k k k k ππππθ(b)N=34, 可能是第II,IV 型滤波⎩⎨⎧≤≤≤≤≤≤≤≤≤≤=2922,125,13330,2113,40,0|)(|k k k k k k H第II 型⎪⎪⎩⎪⎪⎨⎧≤≤-≤≤-=3318),34(3433170,3433)(k k k k k ππθ 第IV 型⎪⎪⎩⎪⎪⎨⎧≤≤-+≤≤--=3318),34(34332170,34332)(k k k k k ππππθ14.解:(a)N 为偶数,上面正交网络可设计成第IV 型滤波器⎪⎪⎩⎪⎪⎨⎧-+-=--+-=---=1,,1]21[),)(2(212]21[,,0,2212)(N N k k N N N N k N k N k ππππθ(b)N 为奇数,纯虚数幅度响应样本为:⎩⎨⎧-=-==1,,1,0,0)(N k j k k jH r 由于这是一个III 型线性相位滤波器,在ω=π处振幅响应应为零,即0=k H 为了减少波动,在靠近ω=π处(即中点两旁)设过渡点,不妨选值为0.4j⎪⎩⎪⎨⎧-+-=-≤≤+-≤≤--===1,1,2/)1(,2/)3(,4.022/)3(,2/)5(2,2/)1(,0,0N N N k j N k N n k j N k k H k⎪⎪⎩⎪⎪⎨⎧-+-=--+-=---=1,,1]21[),)(2)(21(2]21[,0),2)(21(2)(N N k k N N N N k N k N k ππππθ 15.解:(a)(虚数)幅度样本为:⎪⎪⎩⎪⎪⎨⎧-+-=---=-=1,,1]21[),(2]21[,,0,2)(N N k k N N j N k k N jk jH r ππ N 为奇数时没有突变边沿 N 为偶数时没有突变边沿 (b) N 为偶数时⎪⎪⎩⎪⎪⎨⎧-+-=--+-=---=⎪⎪⎪⎩⎪⎪⎪⎨⎧-+=-=--==1,,1]21[),)(2)(21(2]21[,,0),2)(21(2)(1,,12/),(22/),12(212,,0,2N N k k N N N N k N k N k N N k k N N N k N NN k k H k ππππθπππ。
数字信号处理第五章习题解答————第五章————数字滤波网络5.1 学习要点本章主要介绍数字滤波器的系统函数()z H 与其网络结构流图之间的相互转换方法,二者之间的转换关系用Masson 公式描述。
由于信号流图的基本概念及Masson 公式已在信号与系统分析课程中讲过,所以下面归纳IIR 系统和FIR 系统的各种网络结构及其特点。
5.1.1 IIR 系统的基本网络结构1. 直接型结构如果将系统函数()z H 化为标准形式(5.1)式:()∑∑=-=--=Nk kkMk kkz az bz H 11 (5.1) 则可根据Masson 公式直接画出()z H 的直接II 型网络结构流图如图5.1所示(取N=4,M=3)。
二阶直接II 型网络结构最有用,它是级联型和并联型网络结构的基本网络单元。
优点:可直接由标准形式(5.1)或差分方程()()()∑∑==-+-=Mk kN k kk n x b k n y a n y 01画出网络结构流图,简单直观。
缺点:对于高阶系统:(1)调整零、极点困难;(2)对系数量化效应敏感度高;(3)乘法运算量化误差在系统输出端的噪声功率最大。
2. 级联型结构将(5.1)式描述的系统函数()z H 分解成多个二阶子系统函数的乘积形式()()()()z H z H z H z H m 21?= (5.2) (),1221122110------++=zzzzz H i i i i i i ααβββ m i ,,2,1 = (5.3)画出的级联型方框图如图5.2所示。
图中每一个子系统均为二阶直接型结构,根据()z H 的具体表达式确定()z H i 的系数i i i i 1210,,,αβββ和i 2α后,可画出()z H i 的网络结构流图如图5.3所示。
优点:(1)系统结构组成灵活;(2)调整零、极点容易,因为每一级二阶子系统()z H i 独立地确定一对共轭零点和一对共轭极点;(3)对系数量化效应敏感度低。