相平面分析matlab程序

1相平面分析matlab程序《应用非线性控制》，程代展译1、P13：质量-弹簧系统.................................................................................................................12、P14：非线性二阶系统..............................................................................................................23、P15：一阶非线性系统..............................................................................................................34、P17：卫星控制系统..................................................................................................................45、P26：课后习题..........................................................................................................................7题2.2画下列系统的相图.......................................................................................................7题2.4卫星控制系统..............................................................................................................101、P13：质量-弹簧系统0xx-202-2-1012xDx相平面分析clearclcx=1;%修改此值Dx=0;n=1;t=0;Dt=0.001;fori=1:8000DDx=-x;Dx=Dx+DDx*Dt;x=x+Dx*Dt;Dx_store(n)=Dx;x_store(n)=x;2n=n+1;t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel('x')ylabel('Dx')title('相平面分析')holdon2、P14：非线性二阶系统20.630xxxx-10010-20-10010xDx相平面分析clearclc%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:12003Dx=4;n=1;t=0;Dt=0.001;fori=1:1200DDx=-0.6*Dx-3*x-x^2;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')title('相平面分析')holdon3、P15：一阶非线性系统34xxx-202-4-2024xDx相平面分析clearclcx=1.999;%i=1:2000%x=2.001;%i=1:700%x=-1.999;%i=1:30004%x=-2.001;%i=1:700n=1;t=0;Dt=0.001;fori=1:2000Dx=-4*x+x^3;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')title('相平面分析')holdon4、P17：卫星控制系统u505uelseother-505-10-50510xDx相平面分析clearclcx=0.5;%i=1:50005Dx=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;fori=1:5000ifx0u=-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(1)plot(x_store,Dx_store)xlabel('x')ylabel('Dx')title('相平面分析')holdon5、P19：范德波尔方程20.2(1)0xxxx6-202-202xDx相平面分析clearclcx=0;%i=1:20000Dx=10;%x=0;%i=1:7000%Dx=2;%x=0;%i=1:50000%Dx=0.1;n=1;t=0;Dt=0.001;fori=1:20000DDx=-0.2*(x^2-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(1)plot(x_store,Dx_store)xlabel('x')ylabel('Dx')title('相平面分析')holdon75、P26：课后习题题2.2画下列系统的相图(a)0.50xxx-10010-10-50510xDx相平面分析clearclcx=0;%i=1:10000Dx=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=6;%i=1:20000%Dx=-10;n=1;t=0;Dt=0.001;fori=1:200008DDx=-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')title('相平面分析')holdon(b)0.51xxx-10010-10-50510xDx相平面分析clearclcx=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:200009Dx=10;%%x=6;%i=1:20000Dx=-10;n=1;t=0;Dt=0.001;fori=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')title('相平面分析')holdon(c)20.50xxx-6-4-202-20-1001020xDx相平面分析clearclcx=-2;%i=1:4000Dx=20;10%x=0;%i=1:2850%Dx=10;%x=-5;%i=1:10150%Dx=20;%x=-5;%i=1:12000%Dx=7;%x=-6;%i=1:20000%Dx=10;%x=0;%i=1:2600%Dx=20;n=1;t=0;Dt=0.001;fori=1:4000DDx=-Dx^2-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')title('相平面分析')holdon题2.4卫星控制系统u115505uelseother-505-505xDx相平面分析clearclcx=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=0.1;%%%%x=2;%i=1:5000%Dx=2;%%%%x=-2;%i=1:5000%Dx=-2;%%%%%x=3;%i=1:5000%Dx=3;%%%%x=-3;%i=1:5000%Dx=-3;n=1;t=0;Dt=0.001;12fori=1:5000if5*x+Dx0u=-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(1)plot(x_store,Dx_store)xlabel('x')ylabel('Dx')title('相平面分析')holdon