相平面分析matlab程序.docx
相平面分析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_store,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=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*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 on4、P17：卫星控制系统9 = u(-5 9>0 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 x>0u=-5;elseu=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=-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;% 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=-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)相平面分析dearclc 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('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+Dx>0 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
