西南交大信号与系统实验报告.docxVIP

Southwest Jiaotong University 信号与系统 实验一序列的傅立叶变换 实验目的 研究补零问题； 2 ?快速傅立叶变换（FFT）的应用。 实验步骤 复习DFS和DFT的定义，性质和应用； 熟悉MATLAB语言的命令窗口、编程窗口和图形窗口的使用; 利用提供的程序例了编写实验用程序； 按实验内容上机实验，并进行实验结果分析； 写出完整的实验报告，并将程序附在后面。 实验内容 周期方波序列的频谱 ?、 fl, mN nmN-\-L-\ 刀)=\ [0, zzz/V + Z /? (zz? + l)Ar - 1 m = 0,± 1,±2,… = DFS (1)A = 5, /V = 20;⑵L = 5, W = 40; (3)Z = 5, N = 60; (4)0 = 1, N = 60 试画出下面四种情况下的的幅度频谱，并分析补零后，对信号频谱的影响。 Matlab仿真代码： L=5;N=20;%300 % %L = 5;N = 60;%2000 %L = 5 ;N = 40;%1000 %L = 7;N = 60; %2000 % n=l:N; xn=[ones(1,L),zeros(1, N-L)]; figure (1) subplot (2,1,1); stem(n,xn);xlabel(1n1);ylabel(1xtide(n)1); title ( DFS of SQ.wave 1); % Xk = dfs(n,N); len = length (Xk); magXk=abs([Xk(N/2+1:N) Xk(l:N/2+l)]); k=[一N/2:N/2]; subplot (2,2); stem(k,magXk); axis([-N/2,N/2,0,300]); xlabel();ylabeH(Xtide(k) *); % (1)截图 (2)截图 (3)截图 (4) s Figure 4 Fi-e Edit View Insert Too-S Desktop window He-P 0I3E紇邸ze3@^K?o口spra tn DFS of SQ?wave」LH7?NH60 Xtide(k) 2?血彥木逊 0= xasDFr(l)?x(n)(nHs?10)寻凰圧 XS3逋ax(k)ssia; x(3H cos(0?48mz) + cos(0?52/m) 損(1) ￡3x(n)记李悄§卅爭FHFx(n)岂木也=no21oo)寻II壬 x(n)3遽 ?x(k)3?岡； ?x(n)(no21ooaJ 画圧 xss遍ex(k)3ffi^。淫岳 FFrffih 确步审 mat-ab^甜“ (1) MHle NHle nul:M; xnucos (0?48*pi*n) +COS (0?52*pi*n) j :1:N—l〕- yl 丄 xn(l:l:Mrzeros (1、NIM)〕； figure (1) subplot (3、2* 1)； stem (nryl) xlabel (? n ?)； title (- signal x (n) 、OAHnAHlo - ) j Yl = fft (yl); magYl=abs(Yl(1:1:N/2+1)); kl=O:1:N/2; wl= (2*pi/N)*kl; subplot (3,2,2); title(1 Samples of DTFT Magnitude 1); stem(wl/piAmagYl); axis([0,1,0,10]); xlabel(1 frequency in pi units 1 ); % (2) Ml=10; Nl=100; n=l:Ml; xn=cos(0?48*pi*n)+cos(0?52*pi*n); nl=[0:l:Nl-l]; yl=[xn(1:1:M1),zeros(1,Nl-Ml)]; figure (1) subplot (3,2,3); stem(nl,yl);xlabel ( 1n1); title ( 1 signal x(n) r 0=n=1001); Yl=fft(yl); magYl=abs(Yl(1:1:Nl/2+1)); kl=0:l:Nl/2; wl=2*pi/Nl*kl; subplot (3, 2, 4); title(1 Samples of DTFT Magnitude 1); stem(wl/pizmagYl); axis ([0,1,0,10]); xlabel(frequency in pi uni % (3) M2=100; N2=100; n=l:M2; xn=cos(0.48*pi*n)+cos(0.52*pi*n); nl=[0:l:N2-l]; yl=[