




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、實驗3 FFT算法的應(yīng)用實驗?zāi)康模杭由顚﹄x散信號的DFT的理解及其FFT算法的運(yùn)用。實驗原理:N點序列的DFT和IDFT變換定義式如下:N-l1 N-1XI 幻= 叱夕,x=ZX 伙叱丁N &_0- nk利用旋轉(zhuǎn)因子w: =e具有周期性,可以得到快速算法(FFT) 0在MATLAB中,可以用函數(shù)群fft (x, N)和x=ifft (X, N)計算N點序列的DFT正、反 變換。例1對連續(xù)的單一頻率周期信號按采樣頻率采樣,截取長度N分別選N=20和*二16,觀察其DFT結(jié)果的幅度譜。解 此時離散序列,即k=8。用MATLAB計算并作圖,函數(shù)fft用于計算離散傅里葉變換DFT,程序如下:k=8;n
2、l=0:l:19;xal=s in(2*pi*nl/k);subplot (2, 2, 1)plot (nl, xal)xlabel C t/T) ;ylabel (, x(n)J);xkl=fft(xal);xkl=abs(xkl);subplot (2, 2, 2)stem(nl, xkl)xlabel ( k) ; ylabel ( X (k);n2= 0:1:15;xa2=s in(2*pi*n2/k);subplot (2, 2, 3)plot (n2, xa2)xlabel C t/T ) ; ylabel ( x(n);xk2=fft(xa2);xk2=abs(xk2);subp
3、lot (2, 2, 4)stem(n2, xk2)xlabel C k) ; ylabel ( X(k);計算結(jié)果示于圖2.1, (a)和(b)分別是八上20時的截取信號和DFT結(jié)果,由于截取了兩個半周期,頻譜出現(xiàn)泄漏;(c)和(d)分別是P16時的截取信號和DFT結(jié)果, 由于截取了兩個整周期,得到單一譜線的頻譜上述頻譜的誤差主要是由于時域中對信號的 非整周期截斷產(chǎn)生的頻譜泄漏。實驗內(nèi)容:(1) 2N點實數(shù)序列cos(竺 7) + L cos(至 19), = 0,1,2,2N 1N 2 N0,其它N=64o用一個64點的復(fù)數(shù)FFT程序,一次算出X(Z) =??贜卜可,并繪出|X(切(2)已
4、知某序列X()在單位圓上的64等分樣點的Z變換為1-0&,-*包=2-63用N點IFFT程序計算X(n) = IDFTX (%),繪出和X(H)o實驗要求:利用MATLAB編程完成計算,繪出相應(yīng)圖形。并與理論計算相比較,說明實驗結(jié) 果的原因。(1) 用以下代碼實現(xiàn)可得圖6T所示的DFT圖0 N=64; n=0:2*N-l; x=cos(2*pi*7*n/N)+l/2*cos(2*pi*19*n/N); X=fft (x, 128); k=n; stem(k, abs (X) grid xlabel C k*); ylabel C Xk);理論分析如下:由歐拉公式得:xn = cos 有 x 7
5、/0 + geos(三xl9n)對 pn = e,其2N點的DFT變換為:*-j黑盟L黑-Mpm= pne 小= 小“=0=0/ 二 x2M21一6 2當(dāng) 2Z W in 時,“ =0j(2kn)當(dāng)2%=加時,BPp2k = 2N由此可得當(dāng)k=14, 38, 90,1可時有值其余為0(OSk 4 2N-1 )x14 = a-114 = 64, 乂38=凡90 = 32與圖6-1有相同的結(jié)論。(2)用以下代碼可得圖6-2 N=64; k=0:N-l; X=l./(l-0. 8*exp(-j*2*pi*k/N); x = if ft (X, 64); n=k; stem(n, abs (x) gr
6、id圖6-2 xlabelC n );ylabelC xn1);理論分析如下:由M=,,其z變換為X(z) = 一1可得1 一 Z,其 N=64 點的 DFT 為:63 把“ 1-0.8641Xk = y K 乒 再念:1-0.81641-0.8e 得X(q) = X(%) = -j2in,k = J,2,631 UX其IDFT為:m=匚焉乂0.8 5 = 0,1,2 63)比較用IFFT函數(shù)時,得其結(jié)果是復(fù)數(shù)求其模值與理論值進(jìn)行比較得: xl=l/(l-0.8N)*0. 8. delta=xl-abs (x)delta =1.0e-015 *Columns 1 through 900. 11
7、100. 22200. 11100. 22200.16650. 11100. 0278-0. 0555Columns 10 through 18-0. 1110-0. 0416-0. 0833-0. 0833-0. 0625-0. 0278-0. 0347-0. 01390. 0035Columns 19 through 270. 00350. 01730. 0087-0.0347-0. 0755-0. 0668-0. 0572-0. 01000. 0507Columns 28 through 360. 06940. 09820. 12120. 15940.12320. 12310. 13170. 12760. 0911Columns 37 through 450. 03950. 05380. 00980. 0245-0. 0054-0. 0737-0. 0950-0. 1385-0. 1163Columns 46 through 54-0. 0917-0. 0359-0. 0523-0. 0418-0. 0133-0. 0745-0. 0846-0. 0455-0. 0836Columns 55 through 63-0. 0689-0. 0843-0. 1010-0. 0895-0. 0400-0. 0358-0. 01560.
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 開啟2025樂理考試學(xué)習(xí)智庫的有效實踐與運(yùn)用指南試題及答案
- 理清2025年樂理考試中的知識點框架試題及答案
- 如何提升農(nóng)產(chǎn)品電商的社會影響力試題及答案
- 腫瘤早篩技術(shù)臨床應(yīng)用案例研究及2025年市場前景預(yù)測報告
- 系統(tǒng)學(xué)習(xí)樂理的路徑試題及答案
- 新能源汽車電池回收技術(shù)創(chuàng)新與成果鑒定報告2025
- 熱工安全規(guī)程試題及答案
- 家具設(shè)計與社會責(zé)任考核題試題及答案
- 新能源汽車商業(yè)模式的演變與展望試題及答案
- 食品冷鏈物流溫控技術(shù)在生鮮電商中的應(yīng)用前景分析報告
- 六年級語文下冊《(一)字詞積累》期末復(fù)習(xí)課件
- 裝修材料的購銷合同
- 2025年江西金融租賃股份有限公司招聘筆試參考題庫含答案解析
- 臺達(dá)DELTA變頻器VFD-EL系列使用說明書和手冊(完整中文版)VFD007EL23A
- 湖南省長沙市2024-2025學(xué)年高三上學(xué)期新高考適應(yīng)性考試數(shù)學(xué)試題 含答案
- 課題申報書:“四新”建設(shè)背景下教育創(chuàng)新與課程數(shù)字化實踐研究
- 年加工2萬噸再生鋁項目可行性研究報告建議書
- 2025年公司各部門管理規(guī)章制度(4篇)
- 2025年應(yīng)急管理部信息研究院招聘高頻重點提升(共500題)附帶答案詳解
- 2025版《VOCs廢氣處理設(shè)施安全檢查表》(全)
- 普通話水平測試?yán)首x50篇
評論
0/150
提交評論