Animation of the Lorenz Attractor
    5 visualizzazioni (ultimi 30 giorni)
  
       Mostra commenti meno recenti
    
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 Commento
Risposta accettata
Più risposte (0)
Vedere anche
Categorie
				Scopri di più su Animation in Help Center e File Exchange
			
	Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!


