気象学研究室 - フォーラム

ts の関数 spectrum

投稿するにはまず登録を

投稿者	スレッド
森	投稿日時: 2005-4-25 21:01
登録日: 2004-7-29 居住地: 地球投稿: 303	ts の関数 spectrum 履歴 `-rw------- 1 mori mori 1047 Dec 6 2001 MEMO_FOR_TS_SPECTRUM` ts の関数 spectrum() についての覚書・taper = 0, fast =F, detrend = F について調べてみる。 sin(x[1:1000]/1002pi) の spectrum をとってみると、 `result <- spectrum(sin(seq(1,1000)/1002pi),fast=F,taper=0,detrend=F)` すると、 `max(result$spec) sum(result$spec)` ともに 250 になる。これは波数を変えても同じで、 `result <- spectrum(sin(seq(1,1000)/100032pi),fast=F,taper=0,detrend=F) result <- spectrum(sin(seq(1,1000)/1000302pi),fast=F,taper=0,detrend=F)` でもそうである。一方、データ数を変えると話は変わってくる。 `result <- spectrum(sin(seq(1,100)/100302pi),fast=F,taper=0,detrend=F)` なら、 `max(result$spec)` は 25 になる。これをデータ数で割ると、振幅の2乗の1/4 になる。・ちなみにちなみに、fft との関係を調べてみると、 `result <- abs(fft(sin(seq(1,100)/10032pi))) result <- abs(fft(sin(seq(1,1000)/100032*pi)))` の結果、 `max(result)` は 50, 500 になる。こちらは2乗していないことに注意。データ数で割ると、振幅の 1/2 倍になる。

投稿するにはまず登録を