MATLAB求解微分方程组dy(1)/ds=1.2*(0.16*(1-1*cosy(3)+0.5*1*siny(3)))
MATLAB求解微分方程组
dy(1)/ds=1.2*(0.16*(1-1*cosy(3)+0.5*1*siny(3)));
dy(2)/ds=1.2*(0.16*(1-1*cosy(3)+0.5*1*siny(3)))*1*cosy(3)+(y(1)*y(4)*9.8*siny(3))/(300*y(2));
dy(3)/ds=-(1.2*(0.16*(1-1*cosy(3)+0.5*1*siny(3)))*1*siny(3))/y(2)+(y(1)*y(4)*9.8*cosy(3))/(300*y(2));
dy(4)/ds=0;
dy(5)/ds=cosy(3);
dy(6)/ds=siny(3);
初始条件当s=0:y(1)=0.24;y(2)=0.24;y(3)=pi/2;y(4)=144;y(5)=0;y(6)=0