Animation of the Lorenz Attractor

1 回表示 (過去 30 日間)
Westin Messer
Westin Messer 2016 年 12 月 9 日
コメント済み: KSSV 2016 年 12 月 9 日
Hi everybody,
I'm supposed to modify this code so there's a single dot tracing the plot. I think I'm pretty close I just can't completely get it.
function loren3
clear;clf
global A B R
A = 10;
B = 8/3;
R = 28;
u0 = 100*(rand(3,1) - 0.5);
[t,u] = ode45(@lor2,[0,100],u0);
N = find(t>10); v = u(N,:);
x = v(:,1);
y = v(:,2);
plot(x,y);
view(90,90)
M(1)=getframe;
dt=1/100;
for j = 2:1000
t(j)=t(j-1)+dt;
plot(x*t(j),y*t(j),'o','MarkerFaceColor','b','MarkerSize',8)
M(j)=getframe;
if y<=0, break, end
end
pause
movie(M,1)
function uprime = lor2(t,u)
global A B R
uprime = zeros(3,1);
uprime(1) = -A*u(1) + A*u(2);
uprime(2) = R*u(1) - u(2) - u(1)*u(3);
uprime(3) = -B*u(3) + u(1)*u(2);
  1 件のコメント
Jan
Jan 2016 年 12 月 9 日
編集済み: Jan 2016 年 12 月 9 日
The youtube link is not working for me, so I cannot guess,what you want to change. It would be efficient, if you explain this directly instead of letting the readers get this most important detail of your question by using an external web service.

サインインしてコメントする。

採用された回答

KSSV
KSSV 2016 年 12 月 9 日
doc comet
  2 件のコメント
Westin Messer
Westin Messer 2016 年 12 月 9 日
How do I make the comet head go slower?
KSSV
KSSV 2016 年 12 月 9 日
Read the documentation...you can mention time factor as input.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeAnimation についてさらに検索

タグ

タグが未入力です。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by