Describing the motion of a composite body using system of differential equations

Question

Aleem Andrew 2020-3-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/511643-describing-the-motion-of-a-composite-body-using-system-of-differential-equations

评论： Aleem Andrew 2020-3-19

I am trying to solve a system of differential equations analytically using the following code, but although they correctly describe the motion of a body Matlab returns a message stating it is unable to find an explicit solution. Is there a way to fix this problem or an alternate method I can use to solve the system of differential equations?

syms y(t) z(t)
Dy = diff(y,t); Dz = diff(z,t); 
eqns = [diff(y,t,2) == (29.4*cos(z)+4*Dz*Dz)*sin(z)/(1+3*sin(z)*sin(z)), diff(z,t,2) == -3*(19.62+2*Dz*Dz*cos(z))*sin(z)/(2*(1+3*sin(z)*sin(z)))];
[ySol(t),zSol(t)] = dsolve(eqns)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Stephan 2020-3-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/511643-describing-the-motion-of-a-composite-body-using-system-of-differential-equations#answer_420873

在 MATLAB Online 中打开

Solve numeric:

syms y(t) z(t)
Dy = diff(y,t);
Dz = diff(z,t); 
eqns = [diff(y,t,2) == (29.4*cos(z)+4*Dz*Dz)*sin(z)/(1+3*sin(z)*sin(z)),...
    diff(z,t,2) == -3*(19.62+2*Dz*Dz*cos(z))*sin(z)/(2*(1+3*sin(z)*sin(z)))];
fun = matlabFunction(odeToVectorField(eqns),'Vars',{'t','Y'});
tspan = [0 10]; % Time interval for integration
y0 = [0 1 0 1]; % initial conditions
[t,y] = ode45(fun,tspan,y0);
plot(t,y)