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

Signal Processing Toolboxフーリエ解析日本語時系列

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

>フーリエ変換してプロットしたい
こちらはパワースペクトル密度のプロットでよろしいでしょうか。
もしそうで、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
Related Question