BugTrack-Ｒ備忘録/37

Ｒ備忘録 - 記事一覧

ローレンツ方程式のストレンジアトラクタ

• 投稿者： みゅ
• カテゴリ： なし
• 優先度： 普通
• 状態： 完了
• 日時： 2009年06月15日 10時00分01秒

内容

• ローレンツ方程式のストレンジアトラクタ

ローレンツ方程式のストレンジアトラクタ

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")

コメント