2013年6月22日土曜日

Rでアニメーション

Rのanimationライブラリを使ってみた.

plotで出力した画像をくつけてアニメーションにする.

こんな感じ,

library(animation)
f<-4
wave<-function(n){w<-0;t<-seq(0,1,0.01);for(k in 1:n){w<-w+sin((2*k-1)*2*pi*f*t)/(2*k-1)};w<-4/pi*w}
wave2<-function(x){for(n in 1:x){plot(x=seq(0,1,0.01),y=wave(n),type="l",main=paste("Square wave(n=",n,")"),xlab="Time(sec)",ylab="Amplitude")}}
saveMovie(wave2(30),interval=0.5,moviename="wave",movietype="gif",outdir=getwd(),width=640,height=480)

で,できたアニメーションはこちら.

グラフの変化をみせるにはいいツールだね.