Matlab 根据状态方程，绘制相轨迹

以二阶状态方程为例
创建函数draw_phase_locus

function draw_phase_locus(a,b,c)% 状态方程 x''+ a*x' + b*x +c=0
global A;global B;global C; %状态方程三个参数
A=a;B=b;C=c;
clf
for j=-20:4:20   %遍历各种初始条件
for i=-20:4:20
[t,y]=ode45(@nolinear,[0,50],[i,j]);%(微分方程，[时间区间]，[初始条件x',x])
hold on;
plot(y(:,1),y(:,2)); %y(:,1),y(:,2)分别代表 x， x'
end
end
end
function dy =nolinear( t,y ) %建立二阶微分方程，详见ode45的帮助文档
global A;global B;global C;
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-A*y(2)-B*y(1)-C;
end

例：求x’’ + x’ + x -1=0
命令行输入
draw_phase_locus(1,1,-1)
得到

一阶、高阶以此类推, 参考微分方程的帮助文档