相平面分析matlab程序.doc

《相平面分析matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、相平面分析matlab程序《应用非线性控制》，程代展译1、P13：质量-弹簧系统12、P14：非线性二阶系统23、P15：一阶非线性系统34、P17：卫星控制系统45、P26：课后习题7题2.2画下列系统的相图7题2.4卫星控制系统101、P13：质量-弹簧系统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;12n=n+1;t=t+Dt;endfigur

2、e(1)plot(x_store,Dx_store)xlabel('x')ylabel('Dx')title('相平面分析')holdon2、P14：非线性二阶系统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:120012Dx=4;n=1;t=0;Dt=0.001;fori=1:1200DDx=-0.6*Dx-3*x-x^

3、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：一阶非线性系统clearclcx=1.999;%i=1:2000%x=2.001;%i=1:700%x=-1.999;%i=1:300012%x=-2.001;%i=1:700n=1;t=0;Dt=0.001;fori=1:2

4、000Dx=-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：卫星控制系统clearclcx=0.5;%i=1:500012Dx=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

5、.001;fori=1:5000ifx>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(1)plot(x_store,Dx_store)xlabel('x')ylabel('Dx')title('相平面分析')holdon5、P19：范德波尔方程12clearclcx=0;%i=1:20000Dx=10;%x=0;%i=1:7000%Dx=2;%x=0;%i=1:50000%

6、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('相平面分析')holdon125、P26：课后习题题2.2画下列系统的相图(a)clearclcx=0;%i=1:10000Dx=10;%x=0;%i=

7、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:2000012DDx=-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)xlabe

8、l('x')ylabel('Dx')title('相平面分析')holdon(b)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:2000012Dx=10;%%x=6;