Info
This question is closed. Reopen it to edit or answer.
Animation of the Lorenz Attractor
1 view (last 30 days)
Show older comments
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);
1 Comment
Walter Roberson
on 9 Dec 2016
Replaced by better-formatted http://www.mathworks.com/matlabcentral/answers/316118-animation-of-the-lorenz-attractor
Answers (0)
This question is closed.
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!