Question

怎么用matlab求解含参数方程的三元二阶微分方程

t为参数，x=x(t),y=y(t),theta=theta(t)，其中x(0)=0,y(0)=2
方程组为
x''=[-8.66×10^-6*(x'²+y'²)sin²(theta- arctan(y'/x'))cos(arctan(y'/x'))+0.00305*(x'²+y'²)cos²(theta- arctan(y'/x'))sin(arctan(y'/x'))]/808;
y''=[-8.66×10^-6*(x'²+y'²)sin²(theta- arctan(y'/x'))sin(arctan(y'/x'))+0.00305*(x'²+y'²)cos²(theta- arctan(y'/x'))cos(arctan(y'/x'))]/808;
theta''=-10[0.00305*(x'²+y'²)cos²(theta- arctan(y'/x'))cos(theta- arctan(y'/x'))+8.66×10^-6*(x'²+y'²)sin²(theta- arctan(y'/x'))sin(theta- arctan(y'/x'))]/808;
其中x'代表dx/dt，x''代表d²x/dt²，以此类推。能给出代码最好，或者讲一下用matlab的具体解法步骤思路也行

Answer 1

你是不是求五一建模的标枪题，我们短学期也在做这个题，能不能加个好友交流一下呜呜呜

Answer 2

同学，非要解这个形式的方程吗？这个方程明显可以简化，首先令u=x'，v=y'实现降阶，然后将含u，v的方程化成极坐标形式。