内容简介作者简介 - 科学与工程计算系

2.0000 -1.0000
(5)
eA=eig(A)',eA1=eig(A1)'
eA =
-1.0000 0
eA1 =
2.5616 -1.5616
8.6.4 系统平衡点的求取
8.6.5 综合算例
8.6.5.1 “一步仿真”和精良状态轨迹斜率图
2 2
⎧x
1 = x1
+ x2
− 4
【例 8.6.5.1-1】求非线性系统 ⎨
的相平面轨迹、平衡点,并进行稳定性分
⎩x
2 = 2x1
− x2
析。本例综合演示:(A)SIMULINK 模型和 MATLAB 指令的配合使用。(B)sim , simset ,
trim 指令的应用。(C)“一步仿真”计算方法。(D)二阶系统相轨迹的精良图形。
(1)
(2)
[exm080651_1.m]
% exm080651_1.m
clf;hold on
xx=[-2,1;-1,1;0,1;1,1;1,0;1,-1;1,-2];
nxx=size(xx,1);
for k=1:nxx
opts=simset('initialstate',[xx(k,1),xx(k,2)]);
[t,x,y]=sim('exm080634_1',10,opts);
plot(x(:,1),x(:,2));
end
xlabel('x1');ylabel('x2'),grid,hold off
(2)
exm080651_1
(4)
[portraitzzy.m]
x2
1.5
1
0.5
0
-0.5
-1
-1.5
传统状态轨迹图
-2
-2 -1.5 -1 -0.5 0 0.5 1 1.5
x1
图 8.6.5.1-1-1
26