Info

이 질문은 마감되었습니다. 편집하거나 답변을 올리려면 질문을 다시 여십시오.

Animation of the Lorenz Attractor

조회 수: 1 (최근 30일)
Westin Messer
Westin Messer 2016년 12월 9일
마감: Walter Roberson 2016년 12월 9일
Hi everybody,
I'm supposed to modify this code in MATLAB to do something similar to this https://www.youtube.com/watch?v=0WhkP-VyccY basically it's supposed to be 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);

답변 (0개)

태그

Community Treasure Hunt

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

Start Hunting!

Translated by