Animation of the Lorenz Attractor

1 次查看（过去 30 天）

显示更早的评论

Westin Messer 2016-12-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/316114-animation-of-the-lorenz-attractor

关闭： Walter Roberson 2016-12-9

在 MATLAB Online 中打开

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);