如果傅立葉變換提供明顯不準(zhǔn)確的結(jié)果,我能做什么進(jìn)行補(bǔ)救，第4部分

時(shí)間：2024-10-20 22:15:56

關(guān)鍵字：傅立葉變換離散數(shù)據(jù)

手機(jī)看文章

掃描二維碼
隨時(shí)隨地手機(jī)看文章

[導(dǎo)讀]在這一系列中,我們回顧了在微軟的EXECL中實(shí)現(xiàn)的快速傅立葉變換(FFT),并研究了窗口功能。在最后一部分,我們將討論相位測量,但首先,讓我們回顧一下?第三部分 ,我們通過研究窗口不起作用的信號得出結(jié)論。

在這一系列中,我們回顧了在微軟的EXECL中實(shí)現(xiàn)的快速傅立葉變換(FFT),并研究了窗口功能。在最后一部分,我們將討論相位測量,但首先,讓我們回顧一下?第三部分 ,我們通過研究窗口不起作用的信號得出結(jié)論。

為什么窗口不起作用?

圖1a 顯示原始時(shí)域信號和定義它的方程, ω = 2π(39.1 Hz) and φ 表示相位噪聲。計(jì)算 φ ,我用了這個(gè)公式:

φ = -0.5+RAND()

蘭德返回0和1之間的隨機(jī)值,所以表達(dá)式返回-0.5和+0.5之間的隨機(jī)值,表示1拉迪亞尖峰到峰值隨機(jī)相位噪聲。圖1B重復(fù)了第3部分的圖8,但是添加了一個(gè)沒有相位噪聲的信號變換(黑跡),它具有預(yù)期的振幅。這個(gè)問題不是窗口能夠解決的光譜泄漏問題,而是相位噪聲形式的實(shí)際信號損害。

圖1當(dāng)j=0時(shí),樣本信號(a)產(chǎn)生黑色FFT跟蹤(b)。

FFT能測量相位噪聲嗎?

EelelFFT不是這方面的一個(gè)很好的工具,它可以使用前面關(guān)于相位噪聲的系列中討論的技術(shù)來解決(參見?第一部分和?第2部分 )。但FFT能測量出復(fù)FFT中虛部與實(shí)部比例的相角:弧線。

似乎很容易。

當(dāng)心點(diǎn)。考慮藍(lán)色向量圖2 .相角 θ 是0.75/0.5=1.5的弧線,您可以使用elelatan()函數(shù)正確地找到 θ =56.3度。但是嘗試同樣的方法對橙色向量。你的答案是-56.3°(紅色向量),而不是預(yù)期的123.7°。

圖2在計(jì)算Q時(shí)跟蹤象限。

是什么引起的?

一旦你計(jì)算出這個(gè)比率,你就會剝離出 X 和 y 個(gè)體貢獻(xiàn)商的極性。例如,0.75/-0.5(象限2)和-0.75/0.5(象限4)均等于-1.5和ATAN()默認(rèn)值為象限4。作為補(bǔ)償,我們可以添加一些條件編程步驟:"如果x<0和y&tt;0,然后象限2"。幸運(yùn)的是,EXERL包括了ATAN2( X , y 功能,這對我們來說是如此。

我們?nèi)绾涡薷奈覀兊碾娮颖砀駚碛?jì)算相位?

若要在第2部分的圖4電子表格中添加相位計(jì)算,請按紅色的步驟。圖3 .

圖3按照紅色步驟將相位計(jì)算加入FFT。

運(yùn)用圖4A 以波形樣本作為FFT輸入,我們繪制了圖4B的FFT大小和相位--后者在39.1赫茲處約為2.42弧度,即139°。

圖4樣本波形(A)的FFT會產(chǎn)生幅值和相位圖(B)。

啊!我可以買139 在39.1赫茲處°。但余下的階段圖是怎么回事?

讓我們再看一看電子表格,然后向下滾動到39.1赫茲。圖5 )。在任何其他頻率下,復(fù)雜FFT的實(shí)部分和虛部都應(yīng)該是零,而FFT則會返回一個(gè)一點(diǎn)一點(diǎn)的誤差。然而,數(shù)字(藍(lán)色突出顯示)很小,但不是零。我們基本上是用四舍五入誤差和無意義結(jié)果的比率的曲線。