相平面分析matlab程序.docx

《相平面分析matlab程序.docx》由会员分享，可在线阅读，更多相关《相平面分析matlab程序.docx（11页珍藏版）》请在优知文库上搜索。

1、相平面分析matlab程序应用非线性控制，程代展译1、P13：质量辨簧系统 12、P14：非线性二阶系统 23、P15： 一阶非线性系统 34、P17：卫星控制系统45、P26：课后习题7题2.2画下列系统的相图7题2.4卫星控制系统 101、P13：质量-弹簧系统x + x = 0相平面分析21b 0-1-2.-202Xdearclcx=1;%修改此值Dx=O;n=1；占0；Dt=0.001;for i=1:8000DDx=-x;Dx=Dx+DDx*Dt;x=x+Dx*Dt;Dx_store(n)=Dx;xjstore(n)=x;n=n+1;t=t+Dt;endfigure plot(x_s

2、tore,Dx_store)xlabel(x)ylabel(Dx)title。相平面分析)hold on2、P14：非线性二阶系统x + 0.6x + 3x + X2 = 0相平面分析dearclc% x=-6;%i=1:10000 %Dx=10;%x=-8; %i=10000 %Dx=15;%x=5; %i=1:1700 % Dx=7;%x=4;%i=2000% Dx=8;% x=-7; %i=1:1500% Dx=10;%x=-5; %i=1:1200 Dx=4;n=1;t=0;Dt=0.001 ；for i=1:1200DDx=-0.6*Dx-3*x-x2;Dx=Dx+DDx*Dt;x=

3、x+Dx*Dt;Dx_store(n)=Dx;x_store(n)=x;n=n+1;t=t+Dt;endfigure plot(x_store, Dx_store)xlabel(x)ylabel(Dx)titl(相平面分析)hold on3 P15： 一阶非线性系统x = -4x + X3相平面分析42.IS 0-2 ,，-4-202Xdearclc x=1.999； %i=1:2000%x=2.001;%i=1:700%x=-1.999； %i=1:3000 %x=-2.001;%i=1:700n=1;t=0;Dt=0.001 ；for i=1:2000Dx=-4*x+x3;x=x+Dx*D

4、t;Dx_store(n)=Dx; x_store(n)=x;n=n+1;t=t+Dt;endfigure(1)plot(x_store, Dx_store)xlabel(x,)ylabel(,Dx)Utle(相平面分析)hold on4、P17：卫星控制系统9 = u(-5 90 u =(5 elseother相平面分析dearclcx=0.5;%i=1:5000Dx=0.5;%x=1;%i=1:5000%Dx=1;%x=2; %i=1:5000% Dx=2;%x=3; %i=1:5000% Dx=3;n=1;t=0;Dt=0.001;for i=1:5000if x0u=-5;elseu=

5、5;endDDx=u;Dx=Dx+DDx*Dt;x=x+Dx*Dt;Dx_store(n)=Dx;x_store(n)=x;n=n+1;t=t+Dt;endfigure plot(x_store, Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on5、P19：范德波尔方程x + 0.2(2-1)x + x = 0相平面分析-2，-202Xdearclcx=0;%i=1:20000Dx=10;%x=0;%i=1:7000% Dx=2;% x=O;%i=1:50000% Dx=O. 1;n=1；t=0;Dt=0.001;for i=1:20000DDx=

6、-0.2*(x2-1)*Dx-x;Dx=Dx+DDx*Dt;x=x+Dx*Dt;Dx_store(n)=Dx;x_store(n)=x;n=n+1;t=t+Dt;endfigure plot(x_store, Dx_store)xlabel(,x,)ylabel(Dx)title(相平面分析)hold on5、P26：课后习题题2.2画下列系统的相图(a)x + x + 0.5x = 0相平面分析clearclc x=0; %i=1:10000 Dx=10;% x=0; %i=1:10000%Dx=-10;%x=-5; %i=1:15000% Dx=7;%x=5; %i=1:15000% Dx

7、=-7;% x=-6;%i=1:20000%Dx=10;% x=6j%i=1:20000%Dx=-10;n=1；占0；Dt=0.001;for i=1:20000DDx=-Dx-0.5*x;Dx=Dx+DDx*Dt; x=x+Dx*Dt;Dx_store =Dx; x_store(n)=x;n=n+1;t=t+Dt;endfigure plot(x-store1 Dx_store)xlabel(x)ylabel(Dx)UtIeC相平面分析)hold ong x + x + 0.5x = 1相平面分析clearclc x=0;%i=1:10000Dx=10;x=0; %i=1:10000Dx=-

8、10;%x=-5; %i=1:15000Dx=7;%x=5；%i=1:15000Dx=-7;%x=-6;%i=1:20000Dx=10;%x=6;%i=1:20000Dx=-10;n=1；t=0;Dt=0.001;for i=1:20000DDx=1-Dx-0.5*x;Dx=Dx+DDx*Dt; x=x+Dx*Dt;Dx_store(n)=Dx; x_store(n)=x;n=n+1;t=t+Dt;endfigure(1)plot(x_store, Dx_store) xlabel(x,)ylabel(Dx,)Utle(相平面分析)hold onx X2 + 0.5x = 0(c)相平面分析d

9、earclc x=-2;%i=1:4000 Dx=20;%x=0;%i=1:2850%Dx=10;%x=-5; %i=1J0150% Dx=20;% x=-5; %i=1:12000% Dx=7;% x=-6；%i=1:20000%Dx=10;%x=0;%i=12600% Dx=20;n=1；t=0;Dt=0.001;for i=1:4000DDx=-Dx2-0.5*x;Dx=Dx+DDx*Dt;x=x+Dx*Dt;Dx_store(n)=Dx;x_store(n)=x;n=n+1;t=t+Dt;endfigure plot(x_store, Dx_store)xlabel(x)ylabel(

10、Dx)title(相平面分析)hold on题2.4卫星控制系统9 = u(-5 59 +9 O u =5 elseother相平面分析Xdearclc x=5；%i=1:5000Dx=0.01;%x=-5; %i=1:5000% Dx=-0.01;% %x=2;%i=1:5000 %Dx=-0.1;%x=-2; %i=1:5000% Dx=O. 1;%x=2; %i=1:5000% Dx=2;%x=-2; %=1:5000% Dx=-2;% %x=3; %i=1:5000 % Dx=3;%x=-3; %i=1:5000 % Dx=-3;n=1；t=0;for i=1:5000if 5*x+Dx0 u=-5;elseu=5;endDDx=u;Dx=Dx+DD*Dt;x=x+Dx*Dt;Dx_store =Dx; x_store(n)=x;n=n+1;t=t+Dt;endfigure(1)plot(x_store, Dx_store) xlabel(x,)ylabel(,Dx)Utle(相平面分析) hold on