こちらでarrayを使って、３次元格子の時間経過を３次元プロットしている。slice.index()関数を使って３空間軸＋１時間軸の４次元データから、３空間軸の座標を取り出している なるほど、arrayは行列の多次元化したものであって、格子座標が取り出せる これを…

資料 資料のソースをRに書き換えるとfft()の結果と同様の結果になる x<-seq(from=0,to=2*pi,length.out=100) z<-cos(2*x)+2*cos(5*x) par(mfcol=c(2,2)) plot(z,type="l") y<-fft(z) plot(Re(y),type="l") plot(Re(y[1:10]),type="l") a<-fft(y, inverse = T…