時系列プロットのフーリエ変換について

添付の地表温度データを以下のプログラムで時系列プロットし、それをフーリエ変換しプロットしたいため、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)

宜しくお願い致します。

Answer 1

[Takuji Fukumoto]

>フーリエ変換してプロットしたい

こちらはパワースペクトル密度のプロットでよろしいでしょうか。

もしそうで、Signal Processing Toolboxもっていらっしゃれば、 fft()をperiodogram()に置き換えれば良いかと思います。

https://jp.mathworks.com/help/signal/ref/periodogram.html

fft()の関数から計算を行う場合には下記のページが参考になるかと思います。

https://jp.mathworks.com/help/signal/ug/psd-estimate-using-fft.html

[Naoki Ishibashi]

ありがとうございます。うまくできました

[Naoki Ishibashi]

自分の理解が足りていなく申し訳ないのですが、periodogram()でプッロトされたものもまたsin成分とcos成分が重なった値であると思います。 そのためもしご存知でしたら、fft()で出力された値をsinとcos成分に分解しプロットする手法を教えて頂けると幸いです。

[Takuji Fukumoto]

新しく質問を立てていらっしゃったようなので、そちらにコメントさせていただきました。

https://jp.mathworks.com/matlabcentral/answers/312236-fft