MATLAB課程設(shè)計(jì)報(bào)告.doc
一、 設(shè)計(jì)目的通過該設(shè)計(jì),理解傅里葉變換的定義及含義,掌握對(duì)信號(hào)進(jìn)行頻域分析 的方法。 二、 設(shè)計(jì)內(nèi)容1. 信號(hào)的離散傅里葉變換 從連續(xù)到離散: 連續(xù)時(shí)間信號(hào)以及對(duì)應(yīng)的連續(xù)傅里葉變換都是連續(xù)函數(shù)。由于數(shù)字系統(tǒng)只能處理有限長(zhǎng)的離散信號(hào),因此必須將x和都離散化,并且建立對(duì)應(yīng)的傅里葉變換。假設(shè)x(t)時(shí)限于0, L,再通過時(shí)域采樣將離散化,就可以得到有限長(zhǎng)離散信號(hào),記為。設(shè)采樣周期為T,則時(shí)域采樣點(diǎn)數(shù)N=L/T。它的傅里葉變換為這就是在時(shí)域采樣后的連續(xù)傅里葉變換,也就是離散時(shí)間傅里葉變換,它在頻域依然是連續(xù)的。下面將頻域信號(hào)轉(zhuǎn)化為有限長(zhǎng)離散信號(hào)。與對(duì)時(shí)域信號(hào)的處理類似,假設(shè)頻域信號(hào)是帶限的,再經(jīng)過離散化,即可得到有限長(zhǎng)離散信號(hào)。依據(jù)采樣定理,時(shí)域采樣若要能完全重建原信號(hào),頻域信號(hào)應(yīng)當(dāng)帶限于(0,1/T)。由于時(shí)域信號(hào)時(shí)限于0, L,由采樣定理以及時(shí)頻對(duì)偶的關(guān)系,頻域的采樣間隔應(yīng)為1/L。故,頻域采樣點(diǎn)數(shù)為:即頻域采樣的點(diǎn)數(shù)和時(shí)域采樣同為N,頻域采樣點(diǎn)為在DTFT頻域上采樣:令T=1,將其歸一化,就得到離散傅里葉變換。因此,DFT就是先將信號(hào)在時(shí)域離散化,求其連續(xù)傅里葉變換后,再在頻域離散化的結(jié)果。離散傅里葉變換: 離散傅里葉變換(Discrete Fourier Transform,縮寫為DFT),是傅里葉變換在時(shí)域和頻域上都呈離散的形式,將信號(hào)的時(shí)域采樣變換為其DTFT的頻域采樣。在形式上,變換兩端(時(shí)域和頻域上)的序列是有限長(zhǎng)的,而實(shí)際上這兩組序列都應(yīng)當(dāng)被認(rèn)為是離散周期信號(hào)的主值序列。即使對(duì)有限長(zhǎng)的離散信號(hào)作DFT,也應(yīng)當(dāng)將其看作其周期延拓的變換。在實(shí)際應(yīng)用中通常采用快速傅里葉變換計(jì)算DFT。有限長(zhǎng)序列的離散傅里葉變換公式MATLAB函數(shù):fft功能是實(shí)現(xiàn)快速傅里葉變換,fft函數(shù)的格式為: y=fft(x),返回向量x的不連續(xù)fourier變化。ifft功能是實(shí)現(xiàn)快速反傅里葉變換,ifft函數(shù)的格式為: y=ifft(x),返回向量x的不連續(xù)inverse fourier變化。若是一個(gè)N=12的有限序列,利用MATLAB計(jì)算它的傅里葉變換并畫出圖形,然后再對(duì)進(jìn)行離散傅里葉反變換,并求出畫出其波形。2. 頻率分辨率與DFT參數(shù)的選擇在DFT問題中,頻率分辨率是指在頻率軸上所能得到的最小頻率間隔,即最小頻率間隔反比于數(shù)據(jù)的長(zhǎng)度N。若在中的兩個(gè)頻率分別為和的信號(hào),對(duì)用矩形窗截?cái)?,要分辨出這兩個(gè)頻率,N必須滿足 通過下面實(shí)驗(yàn),驗(yàn)證上面的結(jié)論:設(shè)一序列中含有兩種頻率成分,采樣頻率取為,表示為根據(jù)上面結(jié)論,要區(qū)分出著兩種頻率成分,必須滿足N400。 1)取時(shí),計(jì)算的DFT;2)取時(shí),計(jì)算的DFT。三、 總體方案設(shè)計(jì)1. 信號(hào)的離散傅里葉變換 有限長(zhǎng)序列的離散傅里葉變換公式MATLAB函數(shù):fft功能是實(shí)現(xiàn)快速傅里葉變換,fft函數(shù)的格式為: y=fft(x),返回向量x的不連續(xù)fourier變化。ifft功能是實(shí)現(xiàn)快速反傅里葉變換,ifft函數(shù)的格式為: y=ifft(x),返回向量x的不連續(xù)inverse fourier變化。若是一個(gè)N=12的有限序列,利用MATLAB計(jì)算它的傅里葉變換并畫出圖形(見仿真結(jié)果中圖1和圖2),然后再對(duì)進(jìn)行離散傅里葉反變換,并求出畫出其波形(見仿真結(jié)果中 圖3和圖4)。2. 頻率分辨率與DFT參數(shù)的選擇 在DFT問題中,頻率分辨率是指在頻率軸上所能得到的最小頻率間隔,即最小頻率間隔反比于數(shù)據(jù)的長(zhǎng)度N。若在中的兩個(gè)頻率分別為和的信號(hào),對(duì)用矩形窗截?cái)?,要分辨出這兩個(gè)頻率,N必須滿足 通過下面實(shí)驗(yàn),驗(yàn)證上面的結(jié)論:設(shè)一序列中含有兩種頻率成分,采樣頻率取為,表示為根據(jù)上面結(jié)論,要區(qū)分出著兩種頻率成分,必須滿足N400。 1)取時(shí),計(jì)算的DFT,并畫出和的DFT的圖形(見仿真結(jié)果中圖5、圖6和圖7);2)取時(shí),計(jì)算的DFT并畫出和的DFT的圖形(見仿真結(jié)果中圖8、圖9和圖10);試比較兩次實(shí)驗(yàn)是否能區(qū)分出兩種頻率成分。注:圖7和圖8為放大后離散傅里葉變換X(k)四、 主要參數(shù) 1. 是一個(gè)N=12的有限序列。 2. 其中,采樣頻率取為; 1)?。?2)取。五、 源程序1. 離散傅里葉變換及反變換的驗(yàn)證:N=12;n=0:N-1;xn=cos(pi*n/6);k=0:N-1;Wn=exp(-j*2*pi*n*k/N);xk=xn*Wn;figure,subplot(1,2,1),stem(n,xn)subplot(1,2,2),stem(k,abs(xk);y=fft(xn);figure,stem(k,abs(y)Pn=exp(j*2*pi*n*k/N);xn=xk*Pn/N;figure,subplot(1,2,1),stem(n,xk)subplot(1,2,2),stem(k,abs(xn)y1=ifft(xk);figure,stem(k,abs(y1)2. 頻率分辨率與DFT參數(shù)的選擇: f1=2;f2=2.05;fs=10;N=128;n=0:N;xn=sin(2*pi*f1*n/fs)+sin(2*pi*f2*n/fs) figure plot(n,xn) figure y=fft(xn); plot(n,abs(y) N=512;n=0:N;xn=sin(2*pi*f1*n/fs)+sin(2*pi*f2*n/fs); figure plot(n,xn) figure y=fft(xn);plot(n,abs(y)六、 仿真結(jié)果圖1圖2圖3圖4圖5圖6圖7圖8圖9圖10七、 設(shè)計(jì)總結(jié) 通過本次課程設(shè)計(jì),讓我對(duì)信號(hào)的離散傅里葉變換、離散傅里葉反變換、頻率分辨率及DFT參數(shù)的選擇有了更深的認(rèn)識(shí),明白了這一軟件的強(qiáng)大,可以將專業(yè)課信號(hào)與系統(tǒng)里的問題簡(jiǎn)單化、形象化,是一個(gè)好的輔助工具。 另外,對(duì)MATLAB這一軟件的許多專業(yè)領(lǐng)域的功能有了認(rèn)識(shí),強(qiáng)化了對(duì)這一軟件的操作及其想要實(shí)現(xiàn)的功能的編程。與其不斷的翻看課本,不如進(jìn)行實(shí)際的演練,增加了學(xué)習(xí)的趣味性。10