中國最活躍的射頻微波天線雷達微信技術(shù)群
電子獵頭：幫助電子工程師實現(xiàn)人生價值！
電子元器件：價格比您現(xiàn)有供應商最少降低10%
射頻微波天線新產(chǎn)品新技術(shù)發(fā)布平臺：讓更多優(yōu)秀的國產(chǎn)射頻微波產(chǎn)品得到最好的宣傳！發(fā)布產(chǎn)品歡迎聯(lián)系管理，專刊發(fā)布！強力曝光！
做OFDM通信少不了頻譜分析，基帶信號DA后的頻譜，以及基帶數(shù)字上變頻后的DA信號都要頻譜分析。我覺得其實做任何工程都是這樣，先規(guī)定實施方案，然后仿真成功，再實際開發(fā)，不過也可以一邊開發(fā)，一邊仿真，開發(fā)結(jié)果要與仿真預期結(jié)果一致。
所以分析與仿真工具MATLAB就很重要了，既可以仿真，又可以通過示波器或其他方法把實際信號采下來分析。
matlab使用FFT函數(shù)分析信號頻譜
一般我使用的FFT分析頻譜流程如下：
其中有3個注意的點：
1.FFT的結(jié)果看的是頻譜，所以怎么把橫坐標的值從原來的FFT點數(shù)0:N-1轉(zhuǎn)換為頻率值呢？
首先要引出頻譜分辨率的概念，即分辨兩個不同頻率信號的最小間隔，F(xiàn)FT結(jié)果相鄰點間的間隔
。因為N點FFT對應采樣率為fs的序列，其頻率分辨率為
,其中Ts為采樣周期，T為整個序列的時間長度。有關(guān)頻率分辨率的就不多說了。所以我們橫坐標轉(zhuǎn)換為：f = (0:length(y)-1)*Fs/length(y);
2.直接FFT的結(jié)果里怎么又多余的信號頻率（鏡像頻率）圖2？
DFT具有對稱性，因為其是周期序列DFS在一個周期內(nèi)的點，時域序列是有限長實序列，DFT的結(jié)果的實部周期偶對稱，虛部周期奇對稱，也就是模值周期偶對稱，相位周期奇對稱。其實從奈奎斯特定律也可以看出，fs>=2f，fs的采樣率最多也就顯示fs/2的真實頻率（感性理解哈哈）。
所以程序處理方式就是周期延拓后取-N/2:N/2-1.用到函數(shù)fftshift()，結(jié)果如圖3.如注釋所述：
%該變換還會生成尖峰的鏡像副本，該副本對應于信號的負頻率。
%為了更好地以可視化方式呈現(xiàn)周期性，可以使用 fftshift 函數(shù)對變換執(zhí)行以零為中心的循環(huán)平移。
其實這和設(shè)計數(shù)字濾波器IIR與FIR也一樣，采樣率為fs的信號，設(shè)計的濾波器的通帶阻代也限制在0-fs/2內(nèi)。
3.程序中的信號幅度值都是1，500點的FFT畫出來的幅度值怎么變成了250，應該是1吧？
是的，應該是1。所以怎么變換為1呢，注意到FFT的結(jié)果是偶對稱的，且其反應的真實頻譜是0-fs/2。所以需要的操作是直接取0-N/2的FFT結(jié)果，乘以2，然后除以N。即2*abs(y(1:N/2+1))/N，在上面的程序下接著寫：
注意到要除以N，也就是FFT的長度，為什么除以N，這個有很多說法，我自己理性理解的也不透徹，所以這里就不解釋了自行百度（我只能結(jié)合本例子感性理解哈哈）。
結(jié)果如圖4，可見橫坐標頻率是0-25，縱坐標信號的幅度為1.
然而，有不少人（包括我）平時分析信號都是直接FFT畫頻譜：圖2，且看幅值都是看相對的大小，或者有沒有頻率分量，就很少做圖3與圖4的變換。但是我不知道其中的緣由，今天終于明白了。
圖4 FFT所反映的真正的頻譜幅值
參考資料
[1]傅里葉變換 - MATLAB & Simulink - MathWorks 中國
[2]快速傅里葉變換- MATLAB fft - MathWorks 中國
來源：CSDN king阿金
歡迎射頻微波雷達通信工程師關(guān)注公眾號
電子萬花筒平臺自營：Xilinx ALTERA ADI TI ST NXP 鎂光 三星 海力士內(nèi)存芯片 等百余品牌的電子元器件，可接受BOM清單，缺料，冷門，停產(chǎn)，以及國外對華禁運器件業(yè)務！
歡迎大家有需求隨時發(fā)型號清單，我們將在第一時間給您呈上最好的報價，微信(QQ同號)：1051197468 也希望您把我們的微信推薦給采購同事，感謝對平臺的支持與信任！
與我們合作，您的器件采購成本將相比原有供應商降低10%以上！！不信？那您就來試試吧！！歡迎來撩！！
