fft matlabの位相角の計算

Question

私はMatlabのワークスペースで作成した合成信号のフーリエスペクトルの第2成分の周波数領域（fftを計算した後）で位相角を計算しようとしています。私は位相が0に等しいと確信しています（あなたがコードで見ることができます）が、私が得る結果はpi/2です。コードは次のとおりです。

t = 0:pi/128:(2 * pi - pi/128); 
V = sin(t); 
L = length(V); 
n = 2^nextpow2(L); 
Y = fft(V, n); 
threshold = max(abs(Y))/10000; 
Y(abs(Y)<threshold) = 0; 
mag = abs(Y/n); 
angle = rad2deg(atan2(imag(Y),real(Y))); 

エラーの箇所は表示されません。

Answer 1

実際の周期的な位相の位相は、ビンの中心周波数（および位相オフセットなし）に対応する周波数のゼロ波がゼロであると誤解されています。元のシーケンスの一部であるを表す基底関数は、余弦関数である。 （ - PI/2×）

SIN（X）= COS：余弦波と正弦波のπ/ 2の位相オフセットが減算されなければならないを表現するために

。

したがって、ビン2の位相（元のシーケンスの周波数に相当）は-pi/2です。

（より完全な説明のためには、 this question on DSP.SEを参照してください。）