如何從FFT獲得調幅信號的相位角

Question

我有一個調幅信號的測量。 我用fft（） matlab函數分析了。 在我用“書”計算一切之后，我只有一個問題。 調制信號的相位不合適。 只有當我從計算的相位中減去pi/2 ，才能得到正確的值。 調制信號是第六個組成部分：

X [6] = -8.2257e + 001 -1.6158e + 002i
階段（x [6]）= atan（-8.2257e + 001 / -1.6158e + 002）= 1.0999

真正的階段是： pahse(x[6])-pi/2 = -0.4709

為什么我要減去pi/2 ？

如果我使用atan2(imag(X(6)),real(X(6))) - 第一張圖像

如果我使用atan(imag(X(6))/real(X(6)))-pi/2 - 秒圖像

Answer 1

您遇到象限模糊。 atan()的范圍是[-pi/2 ... +pi/2] ，當超出該范圍時重復。 這意味着，當角度恰好位於圓的“另一側”時，您無法唯一地確定角度的正確象限。

要避免這種情況，請使用angle （或phase ）和/或atan2 （ atan的四象限版本）：

>> X = -8.2257e+001 - 1.6158e+002i;
>> angle(X)
ans =
   -2.041680802478084e+000
>> atan2(imag(X), real(X))
ans =
   -2.041680802478084e+000