控制系统设计和仿真.docVIP

熟悉matlab软件的运行环境，包括命令窗体，workspace等，熟悉绘图命令。 2、采用四阶龙格库塔法求如下二阶系统的在幅值为1脉宽为1刺激下响应的数值解。 ， 程序体： T=5; b=0.5; w=10; h=0.01; A=[0 1;-w^2 -2*b*w]; B=[0;w^2]; Y=[0;0]; YY=Y; for t=0:h:T; if(t=1) u=1; else u=0; end K1=A*YY+B*u; K2=A*[YY+h/2*K1]+B*u; K3=A*[YY+h/2*K2]+B*u; K4=A*[YY+h*K3]+B*u; YY=YY+h/6*(K1+2*K2+2*K3+K4); Y=[Y YY]; end y=Y(1,:); t=0:h:(T+h); plot(t,y) 运行结果： 3、采用四阶龙格库塔法求高阶系统阶单位跃响应曲线的数值解。 ，， 主程序： h=0.01; T=15; A=[0,1,0;0,0,1;-20,-102,-10.2]; B=[0;0;20]; Y=[0;0;0]; YY=Y; for t=0:h:T; K1=A*YY+B; K2=A*[YY+h/2*K1]+B; K3=A*[YY+h/2*