レスラーアトラクタ（方程式）

質問

スレ違いかも知れませんがよろしくお願いします！
function resura(itMax);

h=0.01;a=0.398;b=4;c=4;
x=0.12;y=-4.1;z=0.44;
it=0;
A=plot3(x,y,z,\'.-k\',\'MarkerSize\',2,\'EraseMode\',\'none\');
axis([-20 20 -20 20 -20 50]);view(40,60);

%
while it<itMax;
%
xx=x+h*(-y-z);
yy=y+h*(x+a*y);
zz=z+h*(b-c*z+x*z);
set(A,\'XData\',xx,\'YData\',yy,\'ZData\',zz);drawnow
x=xx;y=yy;z=zz;
%;
it=it+1;
end;

上のようなレスラー方程式のプログラムを用いMATLABを起動させました。このプログラムのa,b,c,x,y,zなどの値を変えて平衛解、周期解、２重周機解、カオス解を求めたいのですがどんなに値を変えてもそれらしい解になりません。どなたか分かる方、a,b,c,x,y,zの初期値を教えてもらえれば幸いです。よろしくお願いします。

回答No.1

rabbit_cat · ベストアンサー

wikipedia参照
http://en.wikipedia.org/wiki/R%C3%B6ssler_map

あと、そのプログラム（1次近似）だと、カオス解はうまくでないかもしれません。
少なくとも、ルンゲ＝クッタぐらいにしたほうがいいと思われます。
matlabなら、ode45とかを使えば勝手にやってくれます。
http://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/ode45.shtml