绘制3D相位肖像

时间:2017-01-17 23:20:59

标签: matlab ode

给出方程组:

function test
    A=[9   -21   -10;
       22   -35   -22;
      -22    34    21]; 

    [vals,vector]=eig(A);

    x0=[-0.8,0.8]; %[1,-1, 0.2];%
    y0=[   0,  0]; %[0, 0, 0.4];%
    z0=[-0.9,0.9]; %[1,-1,-0.4];%
    Tmax=2;

    plot3([-10 10],[0 0],[0 0],'k')
    hold on
    grid on
    plot3([0 0],[-10 10],[0 0],'k')
    plot3([0 0],[0 0],[-10 10],'k')

    axis([-10,10,-10,10,-10,10])

    for k=1:length(x0)
        [t,y]=ode45(@sys3d,[0,Tmax],[x0(k),y0(k),z0(k)]);
        plot3(y(:,1),y(:,2),y(:,3))
        [t,y]=ode45(@sys3d,[0,-Tmax],[x0(k),y0(k),z0(k)]);
        plot3(y(:,1),y(:,2),y(:,3))
    end
end

function res=sys3d(t,y)
    A=[-2  1  0;
       -1 -2 -1;
        0  3 -2]; 
    res=A*y;
end

我必须使用以下代码绘制相位图:

x0

我的问题是:

  

y0z0Tmax9x - 21y - 10z + 5来自哪里?如果第一行是   5将{{1}}放在矩阵中哪里?

0 个答案:

没有答案