文档介绍：系统平衡点和普通状态轨线图
【例 -7】求图 -13 所示两个模型的平衡点。模型 (b)输入端比模型 (a)多一个输入口。
1);
DX1=zeros(n,n);DX2=DX1;DP=DX1;
disp('正在逐点计算，请稍等！ ')
for ii=1:n;
for jj=1:n;
opts=simset(opts,'initialstate' ,[x1(ii); x2(jj)]); %<8>
[~,x]=sim( 'exm080707_1',h,opts); %<9>
dx1=x(2,1)-x1(ii);
dx2=x(2,2)-x2(jj);
L=sqrt(dx1^2+dx2^2);DP(jj,ii)=L/h;
if L>-10
DX1(jj,ii)=dx1;DX2(jj,ii)=dx2; %<15>
end
end
end
disp('计算结束 ')
（2）
%
h=;
x1=(-::)';x2=x1;
k=;
set_param('exm080707_1','InitInArrayFormatMsg' , 'None')
xs=trim( 'exm080707_1' ,[-;-]);
xus=trim( 'exm080707_1',[0;1]);
[DX1,DX2,DL]=exm080708_zzy(x1,x2,h);
pcolor(x1,x2,DL)
shading interp
alpha()
colorbar
hold on
quiver(x1,x2,k*DX1,k*DX2,0)
plot(xs(1),xs(2), 'bo',xs(1),xs(2), '+','MarkerSize' ,10)
plot(xus(1),xus(2), 'bo',xus(1),xus(2), '.','MarkerSize' ,10)
grid off
hold off
xlabel('x1'),ylabel('x2')
title( '精良状态轨线斜率图 ')
shg图 -15 精良状态轨线迹斜率图
2 仿真模型和优化指令的协调
【例 -9】题目背景：在迄今的自动控制教材中，凡讨论积分性能指标时，几乎总会提到
所谓的 ITAE 传递函数标准型，并列出相应的分母多项式系数表。但值得指出的是：这些数
据是 20 世纪 50 年代初期，用模拟计算机仿真得到的。因此，这些数据的准确性带有明显
2
的时代缺陷。 与 J ( ISE ) e (t )dt 不同， ITAE 性能函数