Article View: fj.sci.math
Article #13588[Q] : about 3D-graphics with Mathematica
From: M.Masahiro
Date: Sun, 28 Jul 1996 00:00
Date: Sun, 28 Jul 1996 00:00
35 lines
1575 bytes
1575 bytes
みかみ@愛媛大学D3といいます。Mathematica のグラフィックに
ついて教えて頂ければ幸いです。
Arnold 方程式
x'(t)=sin(z(t))+0.3cos(y(t))
y'(t)=0.5sin(x(t))+cos(z(t))
z'(t)=0.3sin(y(t))+0.5cos(x(t))
を初期値
x(0)=3.0, y(0)=0.0, z(0)=3.0
に対して Mathematica 2.2 をつかって数値的に解き、その軌道を立
方体 [0,2pi]^3 上で表示させるため、以下のように実行しました。
sol=NDSolve[{x'[t]==Sin[z[t]]+0.3Cos[y[t]],
y'[t]==0.5Sin[x[t]]+Cos[z[t]],
z'[t]==0.3Sin[y[t]]+0.5Cos[x[t]],
x[0]==3,y[0]==0,z[0]==3},{x,y,z},{t,0,90}];
ParametricPlot3D[Evaluate[{x[t]-2Pi Floor[x[t]/(2Pi)]
,y[t]-2Pi Floor[y[t]/(2Pi)]
,z[t]-2Pi Floor[z[t]/(2Pi)]} /.sol]
,{t,0,90},AspectRatio->1
,AxesLabel->{"x","y","z"}];
岩波講座 応用数学 非線形力学(岡本 久、藤井 宏 著)の30頁の
図 3.1 のようなグラフになると思ったのですが、各座標がジャンプする
ときの軌道も表示してしまい思ったグラフが得られませんでした。
どうすればよいか悩んでいます。お知恵をかしてください。
--
M.Masahiro
mikmas@keimac.cs.ehime-u.ac.jp
Technical Details
Message-ID:
<mikmas-2707961336420001@keimac.cs.ehime-u.ac.jp>Path:
arc1-ca.pugleaf.net!unknown!not-for-mailShow Full Headers
Complete Article Headers
X-Google-Language: Japanese,JIS X-Google-Thread: f4956,52fb81ce65dc0cf2,start X-Google-Attributes: gidf4956,public From: mikmas@keimac.cs.ehime-u.ac.jp (M.Masahiro) Subject: [Q] : about 3D-graphics with Mathematica Date: 1996/07/28 Message-ID: <mikmas-2707961336420001@keimac.cs.ehime-u.ac.jp> X-Deja-AN: 170532866 distribution: fj content-type: text/plain; charset=iso-2022-jp organization: Ehime Univ. mime-version: 1.0 newsgroups: fj.sci.math