1、應(yīng)用傅里葉變換對(duì)信號(hào)進(jìn)行頻譜分析實(shí)驗(yàn)報(bào)告劉運(yùn)鵬5080309226一、 實(shí)驗(yàn)?zāi)康模骸⒓由顚?duì)dft的理解、學(xué)會(huì)應(yīng)用dft對(duì)信號(hào)進(jìn)行頻譜分析的方法二、 實(shí)驗(yàn)內(nèi)容：、復(fù)習(xí)dft和fft有關(guān)內(nèi)容，掌握調(diào)用方法、對(duì)所得結(jié)果加以討論三、 題目一：31為了說(shuō)明高密度頻譜與高分辨率頻譜之間的區(qū)別，考慮序列：求出它基于有限個(gè)樣本的頻譜。a. 當(dāng)時(shí)，確定并畫(huà)出x（n）的離散時(shí)間傅利葉變換。b. 當(dāng)時(shí)，確定并畫(huà)出x（n）的離散時(shí)間傅利葉變換。1. 參數(shù)修改：將0.48改為0.47，0.52改為0.53。2. 程序代碼：n=11;n=0:n-1;x=cos(0.47*pi*n)+cos(0.53*pi*n);k=0:

2、n-1;nk=n*k;wn=exp(j*2*pi*-1/n);wnnk=wn*nk;x=x*wnnk;hx=abs(x);stem(k,hx);n=101;n=0:n-1;x=cos(0.47*pi*n)+cos(0.53*pi*n);k=0:n-1;nk=n*k;wn=exp(j*2*pi*-1/n);wnnk=wn*nk;x=x*wnnk;hx=abs(x);stem(k,hx);3. 實(shí)驗(yàn)結(jié)果和分析：a的dft圖b的dft圖四、 題目二：34設(shè)的基本周期n=50，它的一個(gè)周期由下式給出：的基本周期n=100，它的一個(gè)周期由下式給出：這兩個(gè)周期序列的區(qū)別在于它們的周期，但它們具有相同的非零

3、樣本。a.求出的dfs，用stem函數(shù)畫(huà)出它的幅度和角度樣本。b.求出的dfs，畫(huà)出它的幅度和角度樣本。c.上面兩個(gè)dfs圖有什么不同。提示：本題說(shuō)明時(shí)域補(bǔ)零，則頻域取樣更密1. 參數(shù)修改： 將指數(shù)系數(shù)-0.3改為-0.2。2. 程序代碼：x1=zeros(1,50);for i=1:26 x1(1,i)=(i-1)*exp(-0.2*(i-1);end;x1=fft(x1);x1=x1 x1;k1=0:99;h1=abs(x1);pha1=angle(x1);stem(k1,h1);stem(k1,pha1);x2=zeros(1,100);for i=1:26 x2(1,i)=(i-1)*

4、exp(-0.2*(i-1);end;x2=fft(x2);x2=x2 x2;k2=0:199;h2=abs(x2);pha2=angle(x2);stem(k2,h2);stem(k2,pha2);3. 實(shí)驗(yàn)結(jié)果和分析：后者比前者的采樣更密集x1幅度圖x1相角圖x2幅度圖x2相角圖x幅度圖x相角圖五、 第三題：3612點(diǎn)序列定義為：a求出的dft，畫(huà)出它的幅度和相位曲線。bmatlab畫(huà)出的dtft的幅度和相位曲線。c驗(yàn)證a中的dft是的采樣，采用hold函數(shù)把兩圖放在一幅圖里。d有無(wú)可能從dft重構(gòu)dtft？如果可能，給出重構(gòu)所需的內(nèi)插公式；不可能則說(shuō)明理由。提示：參照題13（c），a調(diào)用fft（），b調(diào)用fplot（）1. 參數(shù)修改：沒(méi)有修改參數(shù)。2. 程序代碼：x=1 2 3 4 5 6 6 5 4 3 2 1;x=fft(x);k=0:11;h=abs(x);pha=angle(x);stem(k*2*pi/12,h);stem(k*2*pi/12,pha);x=1 2 3 4 5 6 6 5 4 3 2 1;n=0:11;w=linspace(0,12,1000000);x=x*exp(-j*n*w);h=abs(x);pha=angle(x);plot(w,h);plot(w,pha);3. 實(shí)驗(yàn)結(jié)果和