信號系統(tǒng)課程設(shè)計(jì).doc

二九 二一學(xué)年第 二 學(xué)期信息科學(xué)與工程學(xué)院課程設(shè)計(jì)報(bào)告書課程名稱： 信號與系統(tǒng)課程設(shè)計(jì) 班 級： 電信0801班 學(xué) 號： 姓 名： wang 指導(dǎo)教師： 盛玉霞 二一 年 六 月 題目一、FFT對連續(xù)信號進(jìn)行譜分析若信號由20Hz，220 Hz和750Hz三個(gè)正弦頻率成分構(gòu)成，試用FFT計(jì)算其頻譜并畫出其幅度譜（要求橫坐標(biāo)以Hz為單位），驗(yàn)證信號包含的頻率成分。分析：設(shè)由于是周期信號，首先要確定該信號的周期，如果截取信號時(shí)沒有至少一個(gè)完整周期的信號，則會(huì)發(fā)生泄漏現(xiàn)象。信號的周期為因此，截取信號的持續(xù)時(shí)間最好為0.1s的整數(shù)倍。信號的最高頻率為，因此，采樣頻率或采樣間隔需滿足，取采樣頻率和采樣間隔分別為：,DFT的點(diǎn)長N從持續(xù)時(shí)間中算得令t=nT, n=0,1,2,199,代入原信號，得利用FFT計(jì)算200點(diǎn)的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對離散信號進(jìn)行譜分析1. 計(jì)算序列的DTFT和DFT，觀察柵欄效應(yīng)設(shè)，要求用MATLAB實(shí)現(xiàn)： （1）計(jì)算的傅里葉變換，并繪出其幅度譜； （2）分別計(jì)算的4點(diǎn)DFT和8點(diǎn)DFT，并繪出其幅度譜。并說明它們和的關(guān)系。（提示：DFT變換可用MATLAB提供的函數(shù)fft實(shí)現(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. 計(jì)算序列的FFT，觀察頻譜泄漏已知周期為16的信號。（1） 截取一個(gè)周期長度M=16點(diǎn)，計(jì)算其16點(diǎn)FFT其頻譜，并繪出其幅度譜；（2） 截取序列長度M=10點(diǎn)，計(jì)算其16點(diǎn)FFT其頻譜，繪出其幅度譜，并與（1）的結(jié)果進(jìn)行比較，觀察頻譜泄漏現(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