《第四章—DTFT答案.doc》由会员分享,可在线阅读,更多相关《第四章—DTFT答案.doc(7页珍藏版)》请在金锄头文库上搜索。
1、1a) b)C) 3.产生图形的MATLAB编码如下所示 %a) x = .25 .25 .25 .25 zeros(1,20); X = dft(x); Omega = 0:pi/100:pi*2; Xa = .5*exp(-j*3*Omega/2).*(cos(3*Omega/2)+cos(Omega/2); N = length(x); k = 0:N-1; figure(1) subplot(211),plot(Omega,abs(Xa),2*k*pi/N,abs(X),o) xlabel(Omega) ylabel(part a) %b) x = 1 -2 1 zeros(1,20)
2、; X = dft(x); Xa = 2*exp(-j*Omega).*(cos(Omega)-1); N = length(x); k = 0:N-1; subplot(212),plot(Omega,abs(Xa),2*k*pi/N,abs(X),o) xlabel(Omega) ylabel(part b)%c) n = 0:5; x = 2*(.75).n; X = dft(x); Xa = 2./(1-.75*exp(-j*Omega); N = length(x); k = 0:N-1; figure(2) subplot(311),plot(Omega,abs(Xa),2*k*p
3、i/N,abs(X),o) xlabel(Omega) title(part c) ylabel(n = 0:5) n = 0:10; x = 2*(.75).n; X = dft(x); N = length(x); k = 0:N-1; subplot(312),plot(Omega,abs(Xa),2*k*pi/N,abs(X),o) xlabel(Omega) ylabel(n = 0:10)n = 0:15; x = 2*(.75).n; X = dft(x); N = length(x); k = 0:N-1; subplot(313),plot(Omega,abs(Xa),2*k
4、*pi/N,abs(X),o) xlabel(Omega) ylabel(n = 0:15)a) 和b)的DTF和DTFT,准确匹配因为xn为有限时间序列。 部分c)的DFT和DTFT。这是一个趋于0的无限时间序列。序列的采样点数越多(N越大),DFT越接近于DTFT45.DFT为DTFT在点的离散。当N=6时,这些点位于6.当T越小时图越准确。当NT变大时分辨率变小。注意对于有限但衰减的信号,采样时间应足够长以便x(t)的截断部分是可以忽略的。 % actual X(w) w = 0:.1:20; Xa = 4./(j*w+1); T = 1; N = 10; t = 0:T:T*(N-1)
5、; x = 4*exp(-t); Xi,wi = contfft(x,T); T = 1; N = 20; t = 0:T:T*(N-1); x = 4*exp(-t); Xii,wii = contfft(x,T); T = 0.5; N = 20; t = 0:T:T*(N-1); x = 4*exp(-t); Xiii,wiii = contfft(x,T); T = 0.1; N = 100; t = 0:T:T*(N-1); x = 4*exp(-t); Xiv,wiv = contfft(x,T); subplot(221),plot(w,abs(Xa),wi,abs(Xi),o) title(T=1,N=10) subplot(222),plot(w,abs(Xa),wii,abs(Xii),o) title(T=1,N=20) subplot(223),plot(w,abs(Xa),wiii,abs(Xiii),o) title(T=0.5,N=20) subplot(224),plot(w,abs(Xa),wiv,abs(Xiv),o) axis(0 20 0 5) title(T=0.1,N=100)