Rでスペクトル解析。
正弦波以外の周期をもった波が倍音が出力されることに注意
例１）正弦波 ３０日周期
x<-sin*1
例２）３０日周期のイベント（倍音がでる）
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
*1:1:3000)/30*2*pi) plot(x) spec.pgram(x,c(10,10
