實驗二實現(xiàn)信號頻譜分析

關(guān)于實驗二實現(xiàn)信號頻譜分析第一頁，共十七頁，編輯于2023年，星期一一、實驗目的進一步加深DFT算法原理和基本性質(zhì)的理解。學習用FFT對連續(xù)信號和時域離散信號進行譜分析的方法，了解可能出現(xiàn)的分析誤差及其原因，以便在實際中正確應用FFT。

應用FFT實現(xiàn)信號頻譜分析第二頁，共十七頁，編輯于2023年，星期一二、實驗原理快速傅立葉變換(FFT)算法長度為N的序列的離散傅立葉變換為：N點的DFT可以分解為兩個N/2點的DFT，每個N/2點的DFT又可以分解為兩個N/4點的DFT。第三頁，共十七頁，編輯于2023年，星期一依此類推，當N為2的整數(shù)次冪時(2N)，由于每分解一次降低一階冪次，所以通過M次的分解，最后全部成為一系列2點DFT運算。以上就是按時間抽取的快速傅立葉變換(FFT)算法。當需要進行變換的序列的長度不是2的整數(shù)次方的時候，為了使用以2為基的FFT，可以用末尾補零的方法，使其長度延長至2的整數(shù)次方。第四頁，共十七頁，編輯于2023年，星期一序列X(k)的離散傅立葉反變換為離散傅立葉反變換與正變換的區(qū)別在于WN變?yōu)閃N-1，并多了一個1/N的運算。因為WN和WN-1對于推導按時間抽取的快速傅立葉變換算法并無實質(zhì)性區(qū)別，因此可將FFT和快速傅立葉反變換（IFFT）算法合并在同一個程序中。第五頁，共十七頁，編輯于2023年，星期一2、MATLAB函數(shù)傅里葉變換函數(shù)常用的FFT及反變換函數(shù)如下表所示。函數(shù)說明fft計算快速離散傅立葉變換fftshift調(diào)整fft函數(shù)的輸出順序，將零頻位置移到頻譜的中心ifft計算離散傅立葉反變換第六頁，共十七頁，編輯于2023年，星期一

fft函數(shù)：調(diào)用方式如下y=fft(x)：計算信號x的快速傅立葉變換y。當x的長度為2的冪時，用基2算法，否則采用較慢的分裂基算法。y=fft(x,n)：計算n點FFT。當length(x)>n時，截斷x，否則補零。第七頁，共十七頁，編輯于2023年，星期一【例2-11】產(chǎn)生一個正弦信號頻率為60Hz，并用fft函數(shù)計算并繪出其幅度譜。第八頁，共十七頁，編輯于2023年，星期一

fftshift函數(shù)：調(diào)用方式如下y=fftshift(x)：如果x為向量，fftshift(x)直接將x的左右兩部分交換；如果x為矩陣（多通道信號），將x的左上、右下和右上、左下四個部分兩兩交換?！纠?-12】產(chǎn)生一個正弦信號頻率為60Hz,采樣率為1000Hz,用fftshift將其零頻位置搬到頻譜中心。例2-1260Hz正弦信號的幅度譜及搬移后的頻譜第九頁，共十七頁，編輯于2023年，星期一

ifft函數(shù)：調(diào)用方式如下y=ifft(x)：計算信號x的傅立葉反變換?！纠?-13】計算方波信號的傅立葉反變換。例2-13矩形信號及其IFFT變換y=ifft(x,n)：計算n點IFFT。如果length(x)>n，以n為長度截短x，否則補零。第十頁，共十七頁，編輯于2023年，星期一三、實驗內(nèi)容考慮長度為5的有限序列，設(shè)采樣周期為0.5s。

要求用FFT來計算其頻譜。

第十一頁，共十七頁，編輯于2023年，星期一三、實驗內(nèi)容用FFT計算下列連續(xù)時間信號的頻譜。

在t<0時，xa(t)=0。已知序列試繪制x(n)及它的離散傅立葉變換|X(k)|圖。第十二頁，共十七頁，編輯于2023年，星期一四、實驗步驟復習DFT的定義、性質(zhì)和用DFT做譜分析的有關(guān)內(nèi)容。復習FFT算法原理與編程思想，并對照DFT-FFT運算流圖和程序框圖，讀懂FFT程序。下圖給出了主程序框圖，供參考第十三頁，共十七頁，編輯于2023年，星期一四、實驗步驟實驗內(nèi)容1的程序運行結(jié)果如下圖所示：第十四頁，共十七頁，編輯于2023年，星期一四、實驗步驟實驗內(nèi)容2的程序運行結(jié)果如下圖所示：第十五頁，共