Hatena::ブログ(Diary)

arupaka-_-arupakaの日記

2017-05-21

Rでスペクトル解析

Rでスペクトル解析。

正弦波以外の周期をもった波が倍音が出力されることに注意

例1)正弦波 30日周期

x<-sin((1:3000)/30*2*pi)

plot(x)

spec.pgram(x,c(10,10))


例2)30日周期のイベント(倍音がでる)

x<-rep(c(rep(0,30)+3*dnorm(1:30,mean=15,sd=2)),100)

plot(x)

spec.pgram(x,c(10,10))

https://jp.mathworks.com/help/signal/examples/practical-introduction-to-frequency-domain-analysis.html

https://jp.mathworks.com/help/signal/examples/analyzing-harmonic-distortion.html