exm09651_0.m

用神经网络来跟踪博弈过程

%exm09651_0.m		例 9.6.5.1_1 的主运行程序

clc
%绘制传统状态轨迹图
figure(1),clf
exm09651_1				

disp('按任意键，画精良状态轨迹图。请稍等！'),pause

%绘制精良的状态轨迹斜率图（Quiver plot）
h=0.01;								%设置仿真步长
x1=-2.5:0.25:2.5;x2=x1;				%轨线起始点
k=0.15;
[X1,X2,Z]=portraitzzy(x1,x2,h);
figure(2),clf
quiver(x1,x2,k*X1,k*X2,0)			%调用quiver指令绘制平面上各点处的变化率图 
xlabel('x1'),ylabel('x2')  
title('精良的状态轨迹斜率图')

disp('按任意键，状态导数向量的长度分布曲面。'),pause

%绘制状态导数向量的长度分布曲面
figure(3),clf
surfc(x1,x2,Z),view([18,32])
xlabel('x1'),ylabel('x2')
title('状态导数向量的长度分布曲面')

%利用trim指令求系统平衡点
xa=trim('exm09634_1',[-1,-2]');
xb=trim('exm09634_1',[1,2]');  
disp('该系统有两个平衡点'),disp('       xa       xb'),disp([xa xb])

%运用线性化模型进行系统稳定性分析
Axa=linmod('exm09634_1',xa);eig_Axa=(eig(Axa))';
Axb=linmod('exm09634_1',xb);eig_Axb=(eig(Axb))';
disp('xa 处线性化模型的特征根为')
disp(eig_Axa)
disp('表明该点为“稳定平衡点”')
disp(' '),disp(' ')
disp('xb 处线性化模型的特征根为')
disp(eig_Axb)
disp('表明该点为“不稳定平衡点”')