スタイル

R

こちらの記事に書いてあったので同じように設定の変更をさせていただいた画面のスクロールの設定

2011-07-30

分岐

存在から発展へ p120 パラメータBの値をかえたときの挙動定常状態になるまで時間がかかるほど計算させる必要がある領域の大きさの影響もないとは言い切れない以下のソースは図 5.10 の領域内で計算するもの本ではでの計算（）図 5.11 のような回転解…

2011-07-28

時間を閉じる

基礎熱力

こちらのコメントから時間的に閉じるとはどういうかとか空間の閉じた・開いたについては空間内外の間で量のやり取りのないとき閉じた系空間内外の間で量のやり取りがあれば開いた系時間的に閉じるとしたら？ある時刻と別の時刻との間で量の移動がな…

2011-07-27

分岐：ブリュッセレータ

シミュレーション R

こちらから p120の図5.7 中心部にちいさな揺らぎを与えてその後の時間発展をみる空間の大きさに影響を受けるようで、半径の円の中だけでシミュレーションしてみたが、円柱対称な波が非常に小さく、平面の様になってしまうほどの大きさで計算しておいて …

2011-07-26

分岐：ブリュッセレータ

R シミュレーション

存在から発展へこちらのつづき p119-p120 拡散係数を補正する処理を入れる p121 図5.6 拡散係数とシミュレーションする空間の境界の計算も修正中心部に小さな揺らぎをあたえるはじめ全体が振動するがやがて対称性の破れた定常状態となるがともに奇数…

2011-07-25

ブリュッセレータの分岐

R シミュレーション拡散

存在から発展へきのうの記事のつづき p119 2次元でブリュッセレータに拡散があるモデルこちらの反応をこちらの場合のように2次元で考えるおそらく1マスの大きさによって拡散係数は変わってくうと思うので、本に出てくる図と同じになる条件を探索中拡…

2011-07-24

分岐：ブリュッセレーター

R シミュレーション拡散

存在から発展へこちらのつづき p116 ブリュッセレーターに拡散を考慮したものとして p117 の定常状態が得られる下のソースは定常状態にはならずに波をうつムービーを参照 p118ページ相当と思われる A<-2 B<-6 X0<-A Y0<-B/A L<-1 dr<-0.01 Nr<-floor(L…

2011-07-23

誤差の大きさ

R シミュレーション

対称性の破れた構造を作り出すための条件を探す初期値に 0,001 ほどの揺らぎを与えると構造が出現するというような値が出た逆にこれ以下の誤差でははじめにいったような構造がでないことを確かめるしかし、シミュレーション幅 dt を小さくすると、必要な…

2011-07-22

物質の交換

シミュレーション R

存在から発展へこちらの続き p107-p108 物質の移動があるモデル箱を2つ用意してそれぞれの箱の中はブリュッセレーターと同じ挙動を示す 2つの箱の間で物質の移動があるが臨界値以上定常解からずらしてみる p108のようなグラフにならず、リミットサイク…

2011-07-21

ブリュッセレーター

R シミュレーショングラフィクス

こちらのつづきリミットサイクルリミットサイクルになる条件 library(rgl) A<-0.5 B<-2 X0<-A Y0<-B/A dt<-0.1 T<-30 Nt<-T/dt M<-matrix(0,3,Nt+1) N<-c() f1<-function(X){ y<- A+X[1]^2*X[2]-B*X[1]-X[1] return(y) } f2<-function(X){ y<- B*X[1]-X[1…