Matlab中pwelch的用法總結

剛剛接觸譜分析的內容，用到了matlab中的pwelch函式，一開始並不理解，看了些文獻，對其有了個大概的了解，進行一下總結。

首先給個基本概念吧

隨機訊號在時間上是無限的，在樣本上是無窮多，因此隨機訊號的能量是無限的，不滿足傅利葉變換的絕對可積條件，因此其嚴格意義上的傅利葉變換是不存在的。因此，對隨即訊號的頻域分析，不再是簡單的頻譜，而是功率譜。—— 《基於welch演算法的簡單功率譜估計的matlab分析》。

說實話，上面這段話，我也沒有很深入的理解，下面舉個例子來看看

matlab程式設計

for i=1000:-1:1

a(i,1)=sin(pi/2*i)+sin(pi/20*i)+cos(pi/50*i);

endpwelch(a,125,60,1000,1)

得到了這樣的，對應的頻率就是訊號的頻率，所以隨機訊號譜分析的結果其實跟其他訊號的傅利葉分解效果差不多，都能夠得到主要頻率，但是我不明白的問題就是為什麼這個譜分析的會有那麼大的振盪呢，傅利葉變換得到的影象比這個清楚多了，所以我覺得挺不理解的。