给出方程组:
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
我的问题是:
y0
,z0
,Tmax
和9x - 21y - 10z + 5
来自哪里?如果第一行是5
将{{1}}放在矩阵中哪里?