x(n)rx(n7r)
求X(k),并作图表示x(n),X(k)。(于2 89 2.1)
解:X(k)%(n)exn0N1j2knNen00j2kn7e01
x(n)-707n
3.2求下列序列的傅里叶变换,并分别给出其幅频特性和相频特性。 (1)x1(n)(nn0) (2)x2(n)3() n3 (3)x3(n)au(n)u(nN)
n13n(4)x4(n)anu(n2) a1 解:(1)对x1(n),有
X1(e)j幅频响应为X(e)ejn0jn(nn)e0jnejn0
1,想频响应为()n0
(2)对x2(n),有
3
X2(e)j1njn 3()e3n33其中
n33ejn36cos6cos26cos3
1111()nejn1ej3ejej29ej2ej327ej3
3927n33将上述两项相加即得X2(ej)
(3)对x3(n)an[u(n)u(nN)],有
3X3(e)aejn0njnaenNnjnaen0N1njn1aNejN 1aej幅频特性为
1aNejN1a2N2aNcosNX3(e)
21aej1a2acosj相频特性为
aNsin(N)asin()()arctan()arctan() j1ae1acos()
(4)对x4(n)au(n2),a1,有
nX4(ej)幅频特性为
nau(n2)ejnnn2aa2nejnae2j 1aej2X(e)相频特性为
ja2e2jj21ae
1a2acos()2arctan(asin()1acos())
3.3 已知以下X(k),求IDFT[X(k)]
NjNe k=m 0m 22Nje k=Nm 1. X(k)= 2 0 其他k
NjN e k=m 0m
22Nje k=Nm 2. X(k)= 2 0 其他k
解:
nmn(Nm)1N11Njj2NNjj2Nkn(1)x(n)IDFT(X(k))X(k)WN[eeee]Nk0N221[e22j(nm)Ne2j(nm)N]cos(2nm)N
n0,1,,N12mn)j(mn)1NjmnNj(Nm)n1j(2N(2)x(n)[jeWNjeWN][eeN]N222j
2sin(mn)n0,1,,N1N3.4 证明DFT的对称原理,即假设:(西电 94 4)
X(k)DFT[x(n)]
证明:DFT[x(n)]Nx(Nk)。 证:因为X(k)x(n)Wn0N1knN,所以
DFTX(n)X(n)WN1knNX(n)Wn0N1knNn(mk)x(m)WNm0n0N1N1kmknx(m)WNWNn0m0
N1由于
Nn(mk)WNn00N1mNk
mNk,0mN1k0,1,,N1
所以DFT[x(n)]Nx(Nk)3.5 证明:若x(n)实偶对称,即x(n)x(Nn),则X(k)也实偶对称,若x(n)实奇对称,
即x(n)x(Nn),则X(k)为纯虚函数并奇对称。(西电 94 7) 证:(1)由教材可知,如果将
x(n)表示为x(n)xrn(jx)in,(则
X(k)DF[T(x)en]pX(op,其中
k)XXepk(k)DFT[xr(n)],是X(k)的共轭对称
分量;Xop(k)DFT[jxi(n)],是X(k)的共轭反对称分量。所以,如果x(n)如果为实序列,则Xop(k)DFT[jxi(n)]0,故X(k)DFT[x(n)]Xep(k),即X(k)=X(Nk) (2) 由DFT的共轭对称性可知,如果x(n)xep(n)xop(n),则
*X(k)Re[X(k)]jIm[X(k)],则Re[X(k)]DFT[xep(n)],jIm[X(k)]DFT[xop(n)]所以当x(n)x(Nn)时,等价于上式中 xop(n)0,所以X(k)x(n)中只有xep(n)成分,只有实部,即X(k)为实函数。又由(1)证明结果知道,实序列的DFT必然为共轭函数,即
X(k)=X*(Nk)X(Nk),所以X(k)为实偶对称。同理当x(n)x(Nn)时,等
价于x(n)只有xop(n)成分(即xep(n)0),故X(k)只有纯虚部,且由于x(n)为实序列,即X(k)共轭对称,X(k)=X*(Nk)X(Nk)为纯虚奇函数。 3.6 证明离散帕塞瓦尔定理。若X(k)DFT[x(n)],则:
1N1|x(n)||X(k)|2 Nk0n02N1证:
1N12|X(k)|Nk0N11*X(k)X(k)Nk0*N1N1knX(k)x(n)WNk0k0knNN1*1x(n)Nk0N1k0*N1X(k)Wk0N1n0N1x(n)x(n)|x(n)|23.7 已知f(n)x(n)jy(n),x(n)与y(n)均为N长实序列。设
F(k)DFT[f(n)], 0kN1
1aN1bN 1. F(k) jkk1aWN1bWN 2. F(k)1jN
试求X(k)DFT[x(n)] ,Y(k)DFT[y(n)]以及x(n)和y(n)。 解:由DFT的共轭对称性可知
x(n)X(k)Fep(k)jy(n)jY(k)Fop(k)
1aN1bN1. 令A(k),只要证明,B(k)jkk1aWN1bWNA(k)Fep(k)X(k),B(k)Fop(k)jY(k),因为
1aN1aN*A(k),共轭对称 A(Nk)Nkk1aW1aWNN1bN1bN*jB(k),共轭反对称 B(Nk)jNkk1bWN1bWN**1aN所以,X(k)Fep(k)A(k) k1aWN111bN Y(k)Fopk ()Bk()kjj1bWN2.
1x(n)N1N由于1NN1k0X(k)Wk0N1N1knN1N1aNknWNkko1aWNN1N1mkmknN1m1aWNWNaNkom0m01=0mnmnWk0N1k(mn)N0nN1Wk(mn)N,0m,nN1所以x(n)an同理y(n)bn0nN10nN1x(n)cos(2n)RN(n) N2y(n)sin(n)RN(n)
N
3.8已知两个有限长序列:
用直接卷积法和DFT变换两种方法分别求解f(n)x(n)y(n)。 解:(1)直接卷积法:f(n)x(n)y(n)
221N122sin(m)cos((nm))[sin((mnm))sin((mnm))]
NN2m0NNm0N12(2mn)j(2mn)1N121N12N21N1j2Nsin(n)sin((2mn)sin(n)[eeN] 2m0N2m0N2N4jm0N2sin(n)RN(n) 2NN1(2)DFT变换法:因为
22njnjkn21N1j2NknX(k)cos(n)RN(n)WN[eeN]eNN2n0n0N1j2(1k)n(1k)n1N1j2NN[eeN][(k1)(kN1)] 2n02n022njnjkn21N1j2NknY(k)sin(n)RN(n)WN[eeN]eNN2jn0n0N1N1j2(1k)n(1k)n1N1j2NN[eeN][(k1)(kN1)] 2jn02jn0所以:
N2F(k)X(k)Y(k)[(k1)(kN1)]
4j1N1Nknn(N1)nf(n)IDFT[F(k)]F(k)WN[WNWN]
Nn04j2nj(N1)nNj2NN2[eeN]sin()RN(n) 4j2N3.9已知x(n)RN(n)
(1)求DFT[x(n)]并画出零极点分布图;
jwX(e) 并画出其幅度和相位图; (2)求频谱
(3)求DFT[x(n)]X(k),并画出其幅度和相位图,且与
N1j2knNn0X(ejw)对照。
解:(1) X(k)DFT[x(n)]e1ej2kNN2kN1ej
因为 X(z)zn0N1n1zN,所以得极点为z1。 11zjw(2)
X(k)X(e)w2kNX(e)
jwnejwnej(N1)w/2sin(wN/2)
sin(w/2)幅频响应为:
X(ejw)jwsin(wN/2)sin(w/2)N1warg[sin(wN/2)/sin(w/2)] 相频响应为:arg[X(e)]2(3)X(k)DFT[x(n)]en0N1j2knN1ej2kNN2kN1ej
X(k)与X(ejw)的关系为:X(k)X(ejw)w2kN
X(k)DFT[x(n)],3.10知x(n)是长度为N的有限长序列。现将x(n)的每二点之间补进r1个零值,得到一个长度为rN的有限长序列y(n)
x(n/r) n=ir, i=0,1,2,3,……,N-1 y(n)=
0 其他n
求:DFT[y(n)]与X(k)的关系。(于2 91 2.13)
解:由于 X(k)DFT[x(n)]rN1x(n)Wn0nkrNN1nkN,0kN
可得 Y(k)DFT[y(n)]N1i0y(n)Wn0irk x(ir/r)WrNi0N1ikx(i)WN,0krN1
所以 Y(k)X((k))NRrN(k)
Y(k)是将X(k)(周期为N)延拓r次形成的,即Y(k)周期为rN。
3.11设x(n)是一个8点的有限长序列,y(n)是一个20点的有限长序列。现将每一序列作20
点DFT,然后再乘,再计算IDFT,令r(n)表示它的离散傅立叶反变换,即:
r(n)IDFT[X(k)Y(。指出k)]r(n)中的那些点相当于x(n)与y(n)的线性卷积中的点。(于
2 94 2-23 )
解:利用循环卷积的公式:两个宽度为N的有限长序列x1(n)和x2(n),其离散傅里叶变换为X1(k)和X2(k),可以求得另外一个序列x3(n),使其离散傅里叶变换的系数为
X1(k)X2(k),则x3(n)的表达式为
x3(n)[x1((m))Nx2((nm))N]RN(n)
m0N1所以,我们利用上式可知循环卷积
r(n)[x((m))20y((nm))20]R20(n)
m0201 而x(n)和y(n)的线性卷积为
x(n)y(n)x(m)y(nm)
m0n
其中,0n27(因为20+8-1=27)时有值,其他时为0。
由于线性卷积在n0,1,2,L,27时有值,而循环卷积在n0,1,2,L19时有值,所以
2,我们以n0,1,L1一考虑r(n)和x(n)y(n)异同处,可以得出:对于逐
n0,1,2,3,4,5,6,两者是不同的,而从n=7开始到n=19,两者是相同的。
3.12.设信号x(n){1,2,3,4},通过系统h(n){4,3,2,1},n0,1,2,3, (1)求系统的输出y(n)x(n)h(n); (2)试用循环卷积计算y(n);
(3)简述通过DFT来计算y(n)的思路。 解:(1)由题可知:x(n)(n)2(n1)3(n2)4(n3)
h(n)4(n)3(n1)2(n2)(n)
y(n)x(n)h(n)x(m)h(nm)x(n)[4(n)3(n1)2(n2)(n)]m034x(n)3x(n1)2x(n2)x(n3)
4(n)11(n1)20(n2)30(n3)20(n4)11(n5)4(n6)
(2)按照循环卷积的方法:已知x(n)与h(n)的长度分别为N=4,M=4,则y(n)的长度应为LNM17,不足补零,则x(n){1,2,3,4,0,0,0},h(n){4,3,2,1,0,0,0}, 卷积得:y(n)x(n)h(n)={4,11,20,30,20,11,4} (3)用DFT计算y(n)的思路:
①分别算出x(n)、h(n)的DFT变换X(k)、H(k); ②再计算出Y(k)X(k)H(k);
③由循环卷积定理得X(k)H(k)x(n)h(n); ④对Y(k)求DFT反变换得y(n)。
因篇幅问题不能全部显示,请点此查看更多更全内容