//このページが表示される方は、URLから「action=SOURCE&」を削除してみてください
[[Ｒ備忘録 - 記事一覧]]

!!!ローレンツ方程式のストレンジアトラクタ
*投稿者： みゅ
*カテゴリ： なし
*優先度： 普通
*状態： 完了
*日時： 2009年06月15日 10時00分01秒
//{{bugstate}}
!!内容
*ローレンツ方程式のストレンジアトラクタ
!!ローレンツ方程式のストレンジアトラクタ
 library(odesolve)
 library(scatterplot3d)
 model <- function(t, xx, parms) {
   u <- xx[1]
   v <- xx[2]
   w <- xx[3]
   with(as.list(parms),{
     du <- -10*u + 10*v
     dv <- 28*u - v - u*w
     dw <- -8/3*w + u*v
     list(c(du, dv, dw))
   })
 }
 times <- seq(0, 200, 0.01)
 parms <- c(a=1, b=1, c=10,
 alpha1=0.2, alpha2=1,
 k1=0.05, k2=0.00, wstar=0.006)
 xstart <- c(u=5, v=10, w=0.1)
 
 out <- as.data.frame(lsoda(xstart, times, model, parms))
 
 par(mfrow=c(2,2))
 plot(times, out$u, type="l", col="green")
 lines(times, out$v, type="l", col="blue")
 plot(times, out$w, type="l", col="red")
 plot(out$w[-1], out$w[-length(out$w)], type="l")
 scatterplot3d(out$u, out$v, out$w, type="l")

//{{ref_image Lorenz.png}}
//{{ref_image Lorenz2.png}}
!!コメント
//{{comment}}