🚀 go-pugleaf

RetroBBS NetNews Server

Inspired by RockSolid Light RIP Retro Guy

Article View: fj.sci.math
Article #13588

[Q] : about 3D-graphics with Mathematica

#13588
From: M.Masahiro
Date: Sun, 28 Jul 1996 00:00
35 lines
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-mail

Show 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