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

8 次查看（过去 30 天）
Naoki Ishibashi
Naoki Ishibashi 2016-11-7
评论： Takuji Fukumoto 2016-11-18
添付の地表温度データを以下のプログラムで時系列プロットし、それをフーリエ変換しプロットしたいため、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)
宜しくお願い致します。

请先登录，再回答此问题。

采纳的回答

Takuji Fukumoto
Takuji Fukumoto 2016-11-9
编辑：Takuji Fukumoto 2016-11-9
>フーリエ変換してプロットしたい
こちらはパワースペクトル密度のプロットでよろしいでしょうか。
もしそうで、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
  3 个评论
显示 1更早的评论
Naoki Ishibashi
Naoki Ishibashi 2016-11-15
自分の理解が足りていなく申し訳ないのですが、periodogram()でプッロトされたものもまたsin成分とcos成分が重なった値であると思います。 そのためもしご存知でしたら、fft()で出力された値をsinとcos成分に分解しプロットする手法を教えて頂けると幸いです。
Takuji Fukumoto
Takuji Fukumoto 2016-11-18
新しく質問を立てていらっしゃったようなので、そちらにコメントさせていただきました。
https://jp.mathworks.com/matlabcentral/answers/312236-fft

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Filter Analysis 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by