相平面图 matlab,相平面分析matlab程序

相平面分析matlab程序

1相平面分析 matlab 程序《应用非线性控制》 ，程代展译1、P13：质量-弹簧系统 12、P14：非线性二阶系统 .23、P15：一阶非线性系统 .34、P17：卫星控制系统 .45、P26：课后习题 .7题 2.2 画下列系统的相图 7题 2.4 卫星控制系统 .101、P13：质量-弹簧系统 0 x-2 0 2-2-1012xDx分分分分分clearclcx=1; %修改此值Dx=0;n=1;t=0;Dt=0.001;for i=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( 相平面分析 )hold on2、P14：非线性二阶系统 20.630 xx-10 0 10-20-10010 xDx分分分分分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;for i=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( 相平面分析 )hold on3、P15：一阶非线性系统 34x-2 0 2-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;for i=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( 相平面分析 )hold on4、P17：卫星控制系统u50elsothr-5 0 5-10-50510 xDx分分分分分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;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(1)plot(x_store,Dx_store)xlabel( x )ylabel( Dx )title( 相平面分析 )hold on5、P19 ：范德波尔方程 20.(1)0 xx6-2 0 2-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;for i=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( 相平面分析 )hold on75、P26：课后习题题 2.2 画下列系统的相图(a) 0.5x-10 0 10-10-50510 xDx分分分分分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;for i=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( 相平面分析 )hold on(b) 0.51x-10 0 10-10-50510 xDx分分分分分cl