exm080651_1.m

matlab mfiles实例分析

% exm080651_1.m		画普通状态轨迹
disp('按任意键，开始运行。')
disp('(1) 在图形窗 N0:1，画普通状态轨线。')
pause
figure(1)
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 
shg

disp(' ')
disp('按任意键，继续运行。')
disp('(2) 在图形窗 N0:2，画精良状态轨线。')
pause
disp(' ')
disp('    正在计算中，请稍等 ！ ')
h=0.01;						
x1=-2.5:0.25:2.5;x2=x1;		
k=3.5;
[X1,X2,Z]=portraitzzy(x1,x2,h);
figure(2)
quiver(x1,x2,k*X1,k*X2,0)	 
xlabel('x1'),ylabel('x2') 
shg

disp(' ')
disp('按任意键，继续运行。')
disp('(3) 在图形窗 N0:2，画状态向量长度构成的曲面。')
pause
figure(3)
surfc(x1,x2,Z)
view([18,32])
xlabel('x1'),ylabel('x2')
shg

disp(' ')
disp('按任意键，继续运行。')
disp('(4) 求系统平衡点')
pause
xa=trim('exm080634_1',[-1,-2]')
xb=trim('exm080634_1',[1,2]')  

disp(' ')
disp('按任意键，继续运行。')
disp('(5) 运用线性化模型进行系统稳定性分析')
pause
Axa=linmod2('exm080634_1',xa);eig_Axa=(eig(Axa))'
Axb=linmod2('exm080634_1',xb);eig_Axb=(eig(Axb))'