数字信号处理经典例题解析
- 格式:doc
- 大小:354.00 KB
- 文档页数:7
第二章习题解答1、求下列序列的z 变换()X z ,并标明收敛域,绘出()X z 的零极点图。
(1) 1()()2nu n (2) 1()()4nu n - (3) (0.5)(1)nu n --- (4) (1)n δ+(5) 1()[()(10)]2nu n u n -- (6) ,01na a <<解:(1) 00.5()0.50.5nn n n zZ u n z z ∞-=⎡⎤==⎣⎦-∑,收敛域为0.5z >,零极点图如题1解图(1)。
(2) ()()014()1414n nn n z Z u n z z ∞-=⎡⎤-=-=⎣⎦+∑,收敛域为14z >,零极点图如题1解图(2)。
(3) ()1(0.5)(1)0.50.5nnn n zZ u n z z --=-∞-⎡⎤---=-=⎣⎦+∑,收敛域为0.5z <,零极点图如题1解图(3)。
(4) [](1Z n z δ+=,收敛域为z <∞,零极点图如题1解图(4)。
(5) 由题可知,101010910109(0.5)[()(10)](0.5)()(0.5)(10)0.50.50.50.50.50.5(0.5)n n nZ u n u n Z u n Z u n z z z z z z z z z z z --⎡⎤⎡⎤⎡⎤--=--⎣⎦⎣⎦⎣⎦⋅=-----==--收敛域为0z >,零极点图如题1解图(5)。
(6) 由于()(1)nn n a a u n a u n -=+--那么,111()(1)()()()nn n Z a Z a u n Z a u n z z z a z a z a a z a z a ----⎡⎤⎡⎤⎡⎤=---⎣⎦⎣⎦⎣⎦=----=-- 收敛域为1a z a <<,零极点图如题1解图(6)。
(1) (2) (3)(4) (5) (6)题1解图2、求下列)(z X 的反变换。
第一章数字信号处理概述简答题:1.在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器,它们分别起什么作用?答:在A/D变化之前为了限制信号的最高频率,使其满足当采样频率一定时,采样频率应大于等于信号最高频率2倍的条件。
此滤波器亦称为“抗混叠”滤波器。
在D/A变换之后为了滤除高频延拓谱,以便把抽样保持的阶梯形输出波平滑化,故又称之为“平滑”滤波器。
判断说明题:2.模拟信号也可以与数字信号一样在计算机上进行数字信号处理,自己要增加一道采样的工序就可以了。
()答:错。
需要增加采样和量化两道工序。
3.一个模拟信号处理系统总可以转换成功能相同的数字系统,然后基于数字信号处理理论,对信号进行等效的数字处理。
()答:受采样频率、有限字长效应的约束,与模拟信号处理系统完全等效的数字系统未必一定能找到。
因此数字信号处理系统的分析方法是先对抽样信号及系统进行分析,再考虑幅度量化及实现过程中有限字长所造成的影响。
故离散时间信号和系统理论是数字信号处理的理论基础。
第二章 离散时间信号与系统分析基础一、连续时间信号取样与取样定理 计算题:1.过滤限带的模拟数据时,常采用数字滤波器,如图所示,图中T 表示采样周期(假设T 足够小,足以防止混叠效应),把从)()(t y t x 到的整个系统等效为一个模拟滤波器。
(a ) 如果kHz rad n h 101,8)(=π截止于,求整个系统的截止频率。
(b ) 对于kHz T 201=,重复(a )的计算。
解 (a )因为当0)(8=≥ωπωj e H rad 时,在数 — 模变换中)(1)(1)(Tj X Tj X Te Y a a j ωω=Ω=所以)(n h 得截止频率8πω=c 对应于模拟信号的角频率c Ω为8π=ΩT c因此 Hz Tf c c 6251612==Ω=π 由于最后一级的低通滤波器的截止频率为Tπ,因此对T8π没有影响,故整个系统的截止频率由)(ωj e H 决定,是625Hz 。
1设序列x(n)={4,3,2,1} , 另一序列h(n) ={1,1,1,1},n=0,1,2,3 (1)试求线性卷积 y(n)=x(n)*h(n) (2)试求6点圆周卷积。
(3)试求8点圆周卷积。
解:1.y(n)=x(n)*h(n)={4,7,9,10,6,3,1}2.6点圆周卷积={5,7,9,10,6,3}3.8点圆周卷积={4,7,9,10,6,3,1,0} 2二.数字序列 x(n)如图所示. 画出下列每个序列时域序列: (1) x(n-2); (2)x(3-n); (3)x[((n-1))6],(0≤n ≤5);(4)x[((-n-1))6],(0≤n ≤5);n12340.5x(3-n)x[((n-1))]n43210.5n12340.5x[((-n-1))6]3.已知一稳定的LTI 系统的H(z)为)21)(5.01()1(2)(111------=z z z z H试确定该系统H(z)的收敛域和脉冲响应h[n]。
解0.52ReIm系统有两个极点,其收敛域可能有三种形式,|z|<0.5, 0.5<|z|<2, |z|>2 因为稳定,收敛域应包含单位圆,则系统收敛域为:0.5<|z|<211111213/25.013/4)21)(5.01()1(2)(--------=---=z z z z z z H)1(232)()5.0(34)(--+=n u n u n h n n4.设x(n)是一个10点的有限序列x (n )={ 2,3,1,4,-3,-1,1,1,0,6},不计算DFT ,试确定下列表达式的值。
(1) X(0), (2) X(5), (3) ∑=90)(k k X,(4)∑=-95/2)(k k j k X eπ解:(1) (2)(3)(4)5. x(n)和h(n)是如下给定的有限序列 x(n)={5, 2, 4, -1, 2}, h(n)={-3, 2, -1 } (1) 计算x(n)和h(n)的线性卷积y(n)= x(n)* h(n); (2) 计算x(n)和h(n)的6 点循环卷积y 1(n)= x(n)⑥h (n); (3) 计算x(n)和h(n)的8 点循环卷积y 2(n)= x(n)⑧h (n); 比较以上结果,有何结论? 解:(1)5 2 4 -1 2-3 2 15 2 4 -1 210 4 8 -2 4-15 -6 -12 3 -6-15 4 -3 13 -4 3 214][]0[190===∑=n N n x X W 12][][]5[119180510-=-===⎩⎨⎧-=∑∑====奇偶奇数偶数n n n n n n x n x X n n W20]0[*10][][101]0[99===∑∑==x k X k X x k k 0]8[*10][][101]))210[((][]))[((2)10/2(92)10/2(910)/2(===-⇔--=-=-∑∑x k X ek X ex k X e m n x k j k k j k m N k j N πππy(n)= x(n)* h(n)={-15,4,-3,13,-4,3,2} (2)5 2 4 -1 2-3 2 15 2 4 -1 210 4 8 -2 4-15 -6 -12 3 -6-15 4 -3 13 -4 3 22-13 4 -3 13 -4 3 2y 1(n)= x(n)⑥h (n)= {-13,4,-3,13,-4,3}(3)因为8>(5+3-1),所以y 3(n)= x(n)⑧h (n)={-15,4,-3,13,-4,3,2,0} y 3(n)与y(n)非零部分相同。
目录习题一 (3)习题二 (26)习题三 (40)习题四 (61)习题五 (83)习题一1.1序列)(n x 如图T1.1所示,用延迟的单位采样序列加权和表示出这个序列。
图 T1.1 习题1.1图【解答】 任一数字序列都可表达为)()()(k n k x n x k -=∑∞-∞=δ所以图T1-1信号可表达为)3(2)1(3)()3(2)(-+-+-+-=n n n n n x δδδδ1.2 分别绘出以下各序列的图形: (1))(2)(1n u n x n =(2))(21)(2n u n x n⎪⎭⎫⎝⎛=(3)()3()2()nx n u n =-(4))(21)(4n u n x n⎪⎭⎫⎝⎛-=【解答】 用MATLAB 得到的各序列图形如图T1.2所示。
图T1.2习题1.2解答1.3 判断下列每个序列是否是周期性的;若是周期性的,试确定其周期。
(1)⎪⎭⎫ ⎝⎛-=873cos )(ππn A n x(2)⎪⎭⎫⎝⎛=n A n x 313sin )(π(3)⎪⎭⎫⎝⎛-=n j e n x 6)(π(4){}{}/12/18()Re Im jn jn x n e e ππ=+(5)16()cos(/17)jnx n e n ππ=【解答】(1)因为730πω=,而31473220==ππωπ,这是一有理数。
所以)(n x 是周期的,周期为14。
(2)因为3130πω=,而136313220==ππωπ,也为有理数。
所以)(n x 是周期的,周期为6。
(3)注意此序列的10=ω,πωπ220=,是无理数,所以)(n x 是非周期的。
(4)实际上()cos(/12)sin(/18)x n n n ππ=+因此)(n x 有两个频率分量,即1201πω=,1802πω=,而 24122201==ππωπ;02223618πππω==都是有理数,所以)(n x 是两个周期信号之和,第一个周期信号的周期241=N ,第二个周期信号的周期362=N ,因此)(n x 的周期是这两个周期的最小公倍数,即72123624)36,24gcd(3624),gcd(2121=⋅=⋅=⋅=N N N N N(5)()x n 是两个周期序列的乘积,其中132N =,234N =,所以该序列的周期是121232343234544gcd(,)gcd(32,34)2N N N N N ⋅⋅⋅====1.4 已知序列)]6()()[6()(---=n u n u n n x ,画出下面序列的示意图。
已知序)()(5n R n x ,求x(n)的8点DFT 变换。
已知模拟滤波器的传输函数 ,用脉冲响应不变法将其转换为数字滤波器,设T=2。
已知采样周期T=2,用双线性变换法将其转换成数字滤波器,说明双线性变换法的有点和缺点。
已知 ,在Z 平面上画出零极点分布图。
已知FIR滤波器的单位脉冲响应为:N=7,h(n)=[3,-2,1,0,1,-2,3] ,说明其相位特性,求群时延。
利用Z变换法求解差分方程描述系统的系统函数H(z)。
1,0)(),(05.0)1(9.0)(-≤==--nnynunyny写出图中流图的系统函数表达式。
已知序列x(n)如图所示,画x((n-2))5R5(n)的图形。
(选做)y(n)1/2 -83 1/4x(n)2Z-1Z-1Z-1求有限长序列x(n)= 的N点DFT。
用脉冲不变法将转换为H(z),采样周期T。
五、计算题(每题12分,共24分)如图所示的RC低通滤波器(1)用脉冲响应不变法转换成数字滤波器。
并画出相应的网络结构图(2)用双线性变换法转换成数字滤波器。
并画出相应的网络结构图(3)以上两种方法所设计的滤波器各自存在那种失真?已知,求两个序列的N=5的循环卷积。
已知系统的差分方程为)2(31)1(32)2()1(2)()(-+---+-+=n y n y n x n x n x n y , (1)求出系统函数(2)画出直接II 型网络结构(3)画出全部一阶节的级联型结构 (4)画出一阶节的并联结构已知序列}4,3,2,1{)(1=n x ,}1,1,1,1{)(2=n x ,求两个序列的线性卷积,和N=5及N=7点的循环卷积。
一个FIR线性相位滤波器的单位脉冲响应是实数的,且n<0 和n>6 时h(n)=0。
如果H(0)=1且系统函数在z=0.5e jπ/3和z=3 各有一个零点,H(z)的表达式是什么?假如x(n)的z变换代数表示式为:(1)求出系统函数所有的零极点;(2)X(z)可能有多少个不同的收敛域?(3)画出不同情况的收敛域图。
第一章 离散时间信号与系统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、对于一个LTI 系统,设其输入序列为矩形冲激信号x(n)=u(n)-u(n-10),而冲激相应为)(9.0)(n u n h n ,用MATLAB 求解输出信号。
可以直接调用卷积函数来实现。
解:clear allx=[1,1,1,1,1,1,1,1,1,1]; n=[0:9]; y=0.9.^n; z=conv(x,y); N=length(z); stem(0:N-1,z);图像如下:给定冲激信号x (n)设定y 函数绘图对x,y 卷积2、编程求两个序列之间的相关系数。
设序列x (k )={3,11,7,0,-1,4,2},n=[-3,-2,-1,0,1,2,3],将x 进行移位再加上一个白噪声信号,即y(k)=x(k-2)+w(k),其中k 属于n ,需要计算x 序列与y 序列之间的相关系数,可以使用卷积来实现。
解:clear all>> x=[3,11,7,0,-1,4,2]; >> nx=[-3:3];>> [y,ny]=sigshift(x,nx,2); >> w=randn(1,length(y)); >> nw=ny;>> [y,ny]=sigadd(y,ny,w,nw); >> [x,nx]=sigfold(x,nx);>> [rxy,nrxy]=conv_m(y,ny,x,nx); >> subplot(1,1,1); >> stem(nrxy,rxy)>> axis([-5,10,-50,250]); >> xlabel('延迟量1'); >> ylabel('rxy');>> title('噪声序列的互相关')图像如下:给定信号x (n)对x 序列移位设定随机信号w根据x,w 得到y 序列对x,y 卷积绘图3、利用filter函数计算冲激相应和单位阶跃响应。
1:周期序列()()n n x 0cos ~ω=, 0ω6π=,()n x ~是由)(~t x a ()t 0cos Ω=理想抽样而得。
试求(1)()n x ~的周期;(2)()()[]n x F e X j ~=ω (3) ()t x a~=∑∞-∞=n ntj n 0e Ωα;求n α (4)()()[]t x F X a ~=Ω 解:(1) 对于周期性序列()()n n x 0cos ~ω= 因为2ωπ=6/2ππ=112=K N所以序列周期12=N(2):由题意知()n x ~是由()t x a ~理想抽样所得,设抽样间隔为s T ,抽样输出为()t xa ˆ; 易得()()[]t x F X a ~=Ω()[]t F 0cos Ω= ]2[00tj t j e e F Ω-Ω+==π()0Ω+Ωδ+π()0Ω-Ωδ由采样序列()n x ~=()nt xa ˆ,由采样定理知: ()()[]n x F e X j ~=ω=()sTaX /ˆω=ΩΩ =∑∞∞--k ss sT k T X T )2(1πω =∑∞∞--k s s T k X T )2(1πω=)]26()26([1sk s s T k T k T ππωπδππωπδ-++--∑∞∞- =)]26()26([ππωπδππωπδk k k -++--∑∞∞-(3) 由)(~t x a ()t 0cos Ω==200tj tj e e Ω-Ω+=∑∞-∞=n nt j n 0e Ωα得:⎪⎩⎪⎨⎧=±==其他n n n 0121α(4)由(2)得:()ΩX =π()0Ω+Ωδ+π()0Ω-Ωδ2:有限长序列()⎪⎭⎫⎝⎛=n n x 6cos π()n R 12求:(1))]([)(n R F e R n j n =ω(2) ()()[]n x F e X j =ω,用)(ωj N e R 表示; (3)求(2)中()ωj e X 的采样值⎪⎪⎭⎫ ⎝⎛k j e X 122π 110≤≤k ; (4)()()[]n x DFT k X =;(5):求第(3)问中⎪⎪⎭⎫ ⎝⎛k j e X 122π的IDFT 变换; (6):求()()⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛=n R n F eX j 2416cos πω的采样值⎪⎪⎭⎫ ⎝⎛k j e X 2421π230≤≤k ; (7):求第(6)问中的采样序列()n x 1; (8):第(2)问中()ωj e X 的采样值⎪⎪⎭⎫ ⎝⎛k j e X 242π对应的采样序列。
.解:(1))]([)(n R F e R n j n =ω=∑-=1)(N n nj N e n R ω=2/2/2/2/2/2/11ωωωωωωωωj j N j N j j N j j N j ee e e e e e e --------⋅=-- =)2/sin()2/sin(2/2/ωωωωN e e j N j ⋅--容易看出在主值周期内当0=ω时)(ωj n e R N =,当)10(2-≤≤=N k Nk πω时)(ωj n e R =0 (2)根据公式()()[]()()ωωπj j e Y e X n y n x F *=21=θπθωππθd e Y e X j j )()(21)(--⎰又由⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛n F 6cos π=)]26()26([ππωπδππωπδk k k -++--∑∞∞-则 ()()[]()⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛==n R n F n x F eX j 126cos πω=θππθπδππθπδπππθωd e R k k k j ⎰∑-∞-∞=--++--)()]26()26([21)(12 =θπθπδπθπδπππθωd e R j ⎰--++-)()]6()6([21)(12 =+-)(21)6(12πωj e R )(21)6(12πω+j e R(3)易知⎪⎪⎭⎫ ⎝⎛k j e X 122π是对()ωj e X 进行间隔为122π的等间隔采样所得,由(1)知在主值周期内当0=ω时)(ωj n e R N =,当)10(2-≤≤=N k Nk πω时)(ωj n e R =0又因为()ωj e X =+-)(21)6(12πωj eR )(21)6(12πω+j eR 所以在主值周期内πππππ26122,6122=+=-k o k 时即11,1=k 时6212122==⎪⎪⎭⎫ ⎝⎛k j e X π,当k 10,9,8,7,6,5,4,3,2,0=时0122=⎪⎪⎭⎫ ⎝⎛k j e X π; 即:⎪⎪⎭⎫ ⎝⎛k j e X 122π⎩⎨⎧===10,9,8,7,6,5,4,3,2,0011,16k k (4)根据DFT 变换与序列傅里叶变换的关系:()k X =kN πωω2j ) X(e ==⎪⎪⎭⎫ ⎝⎛k Nj e X π2其中N 为()n x 的周期 所以由(3)得: ()()[]n x DFT k X =⎪⎪⎭⎫ ⎝⎛=k j e X 122π⎩⎨⎧===10,9,8,7,6,5,4,3,2,0011,16k k (亦可根据公式直接求解:()()[]n x DFT k X ==kn j n en x 122110)(π-=∑=kn j n n j n jee e122110122122)(21πππ-=-∑+=∑=-+-+110)1(122)1(122)(21n k n j k n j e e ππ ⎩⎨⎧===10,9,8,7,6,5,4,3,2,0011,16k k(5)由(3)(4)易得()()n x DFT e X k j =⎪⎪⎭⎫ ⎝⎛122π; 所以()()n R n n x e X IDFT k j 121226cos ⎪⎭⎫ ⎝⎛==⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛ππ(6)由(1)(2)可得()()⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛=n R n F e X j 2416cos πω=+-)(21)6(24πωj e R )(21)6(24πω+j eR 所以: 在主值周期内πππππ26242,6242=+=-k o k 时即22,2=k 时122242421==⎪⎪⎭⎫ ⎝⎛k j e X π, 当k 23,21,20,19,18,17,16,15,14,13,12,11,10,9,8,7,6,5,4,3,1,0=时0122=⎪⎪⎭⎫ ⎝⎛k j e X π; 即⎪⎪⎭⎫ ⎝⎛k j e X 2421π⎩⎨⎧===23,21,20,19,18,17,16,15,14,13,12,11,10,9,8,7,6,5,4,3,1,0022,212k k (7): 根据DFT 变换与序列傅里叶变换的关系()k X kN πωω2j ) X(e ===⎪⎪⎭⎫ ⎝⎛k Nj e X π2 则⎪⎪⎭⎫ ⎝⎛k j e X 2421π()⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛=n R n DFT 246cos π 所以采样序列()()n R n n x 2416cos ⎪⎭⎫⎝⎛=π(8)易知⎪⎪⎭⎫ ⎝⎛k j e X 242π是对()n x 的24=N 的DFT 变换 所以对应的采样序列()=n x 2⎩⎨⎧≤≤≤≤23n 120110)(n n x3:()⎪⎭⎫⎝⎛=n n x 6cos ~π由)(~t x a ()t 0cos Ω=理想采样所得(1)求()()[]n x DFS k X ~~=,并求出主值周期内的值(2)()()()[]n R n x DFT k X 2424~=,()k X 24能准确的反应)(~t x a ()t 0cos Ω=的频率成分吗?为什么?解:(1)易知()()[]n x DFS k X ~~=为()()()[]n R n x DFT k X 12~=以12=N 周期性延拓得到,因为()()()[]n R n x DFT k X 12~=⎩⎨⎧===10,9,8,7,6,5,4,3,2,0011,16k k所以()()[]n x DFS k X ~~=()⎩⎨⎧±==kl l k 其他为整数0,1126主值周期内的值为:()()()[]n R n x DFT k X 12~=⎩⎨⎧===10,9,8,7,6,5,4,3,2,0011,16k k(2):能,(原因自己分析)。
4:求解系统因果解:()()()1-+=n y n x n y ,其中)()(n n x δ= 解:对方程两边做Z 变换得:)()()(1z Y z z X z Y -+=所以()=z Y 11111)(---=-zz z X 1>z反Z 变换得()()n u z Z n y =⎥⎦⎤⎢⎣⎡-=--1111 5:(1):求线性卷积:()()n R n R 22*(2):求圆周卷积:()n R 2○3()n R 2 解(1)()n R 2的支撑区为[]1,0;则()()n R n R 22*的支撑区为[]2,0 所以()()=*n R n R 22)()(22m n R m R m -∑∞-∞==)()(212m n R m R m -∑=有()()=*0022R R 1)0()(212=-∑=m R m Rm()()=*1122R R 2)1()(212=-∑=m R m R m()()=*2222R R1)2()(212=-∑=m R m R m即:()()=*n R n R 22 ⎪⎩⎪⎨⎧==为其他正整数n n n 0122,01 (2)()n R 2的线性卷积支撑区长度为2=N ,由于圆周卷积的周期313=-+==N N L ,根据线性卷积与圆周卷积的关系可得:()n R 2○3()n R 2=()()()n R n R n R 322]~~[⋅*=()()n R n R 22*=⎪⎩⎪⎨⎧==为其他正整数n n n 0122,01 6:画时域抽样8点FFT 的流程图。
解:见课本图3--5(P 109)7:1(3)中n α与3(1)()k X 满足等式∑∞-∞=+=l lN k N k X )()(α,12=N 。
(自证)。