MATLAB: Fft関数の出力データの成分分解

fftフーリエ変換日本語

添付の地表温度データを以下のプログラムで時系列プロットし、それをフーリエ変換しプロットしたいため、fft関数を使おうと考えたのですがうまくいきません。例題ページを見たところ、データをsin成分とcos成分に分解してからプロットしなくては行けないのかなと考えたのですが、その方法もわからず、もしご存知の方いらしましたら教えて頂きたいです。 以下時系列プロットのプログラムです
filename =('1.txt');
x = load(filename);
x = timeseries(x);
x.Name = 'Daily Count';
x.TimeInfo.Units = 'days';
x.TimeInfo.StartDate = '04-06-01 00:00';     % Set start date.
%x.TimeInfo.End = '04-08-31 21:00';
x.TimeInfo.Units = 'hours';
x.TimeInfo.Increment = 3;
x.TimeInfo.Format = 'yy-mm-dd HH:MM';       % Set format for display on x-axis.
x.Time = x.Time - x.Time(1);        % Express time relative to the start date.
plot(x)
以下フーリエ変換です
filename =('1.txt');
x = load(filename);
y = fft(x);
plot(y)

Best Answer

画像のプロットの理解については推測が入ってしまいますが、
fft()関数での離散フーリエ変換では
Wn=e^(−2πi)/nの項があり、複素平面で取り扱われています。
e^iθ　＝cosθ+ i sinθであらわせるので、
画像のプロットでは実軸の数値をcosin、虚軸の数値をsinと表示しているのかと推測します。
FFT演算後のｙは複素数になっているので、 実数部と虚数部を取り出してそれぞれをプロットすればよいのではないでしょうか。
実数部を取り出す　→　real（）
虚数部を取り出す　→　imag（）
でできます。
Related Question