信號系統(tǒng)課程設計.doc

二九 二一學年第 二 學期信息科學與工程學院課程設計報告書課程名稱： 信號與系統(tǒng)課程設計 班 級： 電信0801班 學 號： 姓 名： wang 指導教師： 盛玉霞 二一 年 六 月 題目一、FFT對連續(xù)信號進行譜分析若信號由20Hz，220 Hz和750Hz三個正弦頻率成分構成，試用FFT計算其頻譜并畫出其幅度譜（要求橫坐標以Hz為單位），驗證信號包含的頻率成分。分析：設由于是周期信號，首先要確定該信號的周期，如果截取信號時沒有至少一個完整周期的信號，則會發(fā)生泄漏現(xiàn)象。信號的周期為因此，截取信號的持續(xù)時間最好為0.1s的整數(shù)倍。信號的最高頻率為，因此，采樣頻率或采樣間隔需滿足，取采樣頻率和采樣間隔分別為：,DFT的點長N從持續(xù)時間中算得令t=nT, n=0,1,2,199,代入原信號，得利用FFT計算200點的DFT。fs=2000;T=0.0005;tp=0.1;N=tp/T;n=0:N-1;k=n;w=2*pi/N*k;x=0.25*sin(2*pi*20*n*T)+0.5*sin(2*pi*220*n*T)+sin(2*pi*750*n*T);y=fft(x);stem(w/(2*pi),abs(y),b)hold on;xlabel(hz)title(200dft)圖：題目二、FFT對離散信號進行譜分析1. 計算序列的DTFT和DFT，觀察柵欄效應設，要求用MATLAB實現(xiàn)： （1）計算的傅里葉變換，并繪出其幅度譜； （2）分別計算的4點DFT和8點DFT，并繪出其幅度譜。并說明它們和的關系。（提示：DFT變換可用MATLAB提供的函數(shù)fft實現(xiàn)，也可以自己用C語言或matlab編寫）題一function f=u(t)f=(t=0);T=1;dw=0.1;n=-5:5;w=-7*pi:dw:7*pi;f=u(n)-u(n-4);F=f*exp(-j*n*w)*T;ff=abs(F);plot(w,ff)title(R(N)的傅里葉變換);圖：題二clear;close allN=4;N1=8;n=0:N-1;k=0:N1-1;xn=u(n)-u(n-4);xk1=fft(xn,N);xk2=fft(xn,N1);subplot(121);stem(n,abs(xk1),.);xlabel(n);ylabel(abs(xk1);title(4DFTx1(n);subplot(122);stem(k,abs(xk2),.);xlabel(k);ylabel(abs(xk2);title(8DFTx2(n);2. 計算序列的FFT，觀察頻譜泄漏已知周期為16的信號。（1） 截取一個周期長度M=16點，計算其16點FFT其頻譜，并繪出其幅度譜；（2） 截取序列長度M=10點，計算其16點FFT其頻譜，繪出其幅度譜，并與（1）的結(jié)果進行比較，觀察頻譜泄漏現(xiàn)象，說明產(chǎn)生頻譜泄漏的原因。clear;close allN=16;N1=10;n=0:N-1;k=0:N1-1;xn=cos(10*pi*n/16)+sin(12*pi*n/16);xk1=fft(xn,N);xk2=fft(xn,N1);subplot(121);stem(n,abs(xk1),.);xlabel(n);ylabel(abs(xk1);title(16DFTx1(n);subpl