How to solve a system of nonlinear 2nd order differential equations?

Question

Franziska 2013-2-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/64417-how-to-solve-a-system-of-nonlinear-2nd-order-differential-equations

Hi there,

I have a challenge solving a system of differential equations 2nd order. I do not receive an error message but have rather strange results... Could you please have a look on what I did and comment on it. If you think the way it was done is fine, I may have the bug some where else...

 function dHbar = diffH(tspan,H)   
   h1   = H(1);
   dh1  = H(2); 
   ddh1 = H(3);
   h2   = H(4);
   dh2  = H(5);
   ddh2 = H(6);
   h3   = H(7);
   dh3  = H(8);
   ddh3 = H(9);
   % etc
% the equations are very long. In short the look similar to the following where f(h1,dh1) means some expression as a function of ...
% A1 to A3 are known
 ddh1 = ( A1 * f(h1,dh1) + A2 * f(h2,dh2,ddh2) + A3 * f(h3,dh3,ddh3));
 ddh2 = ( A2 * f(h2,dh2) + A1 * f(h1,dh1,ddh1) + A3 * f(h3,dh3,ddh3));
 ddh3 = ( A3 * f(h3,dh3) + A1 * f(h1,dh1,ddh1) + A2 * f(h2,dh2,ddh2));
dHbar = [h1; dh1; ddh1; h2; dh2; ddh2; h3; dh3; ddh3];

Then I call it with

[T,Hbar] = ode45('diffH',tspan,H); % with tspan - time and H - starting values

As mentioned above, I receive some solution, but it seems strange.

Thanks a lot for your help in advance. Cheers, Franziska

7 个评论
显示 5更早的评论隐藏 5更早的评论

Babak 2013-2-25

在 MATLAB Online 中打开

Franziska, now, I'm understanding yoru problem better.

I think you have not formulated your problem in the proper way that ODE45 requires you to.

In the function dHbar = diffH(tspan,H) the return value of this function, should only include the velocity and acceleration. If you see

doc ode45

you will notice that the return of the function that needs to be integrated with ODE45 is the derivative of the input. So you need to have position and velocity as input and velocity and acceleration as output.

You need to write H as (h1, dh1,h2,dh2,h3,dh3) and then write your fucntion like this:

 function dHbar = diffH(tspan,H)      
   h1   = H(1);
   dh1  = H(2); 
   h2 = H(3);
   dh2   = H(4);
   h3  = H(5);
   dh3 = H(6);
% dHbar(1) = dh1;
% dHbar(2) = ddh1;
% dHbar(3) = dh2;
% dHbar(4) = ddh2;
% dHbar(5) = dh3;
% dHbar(6) = ddh3;
% already known, need not be computed:
dHbar(1) = dh1;
dHbar(3) = dh2;
dHbar(5) = dh3;
ddh1 = ( A1 * f(h1,dh1) + A2 * f(h2,dh2,ddh2) + A3 * f(h3,dh3,ddh3));
ddh2 = ( A2 * f(h2,dh2) + A1 * f(h1,dh1,ddh1) + A3 * f(h3,dh3,ddh3));
ddh3 = ( A3 * f(h3,dh3) + A1 * f(h1,dh1,ddh1) + A2 * f(h2,dh2,ddh2));
dHbar = [dh1; ddh1; dh2; ddh2; dh3; ddh3];
 end

Franziska 2013-2-27

Thanks Babak for your help and patience. I have now rewritten the entire system to M(y,t)*y' = f(t,y) and work with ode23t due to singular matrix. Thanks again for your help! Franziska

Babak 2013-3-4

I think after your write your equations in the form I'm mentioning above, you can solve it with any ODE solver like ODE45. It doesn't seem to be a complicated case.

请先登录，再进行评论。

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

Answer 1

Babak 2013-2-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/64417-how-to-solve-a-system-of-nonlinear-2nd-order-differential-equations#answer_76007

You need to use the optimization toolbox and it's fsolve(.) routine. You can embed all your nonlinear functions (can be an integrated ODE or whatever)as the arguments of fsolve, assigning thos functions to be the ones that are needed to be optimized.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Franziska 2013-2-21

Thanks for your answer Babak, according to the matlab-help fsolve is not applicable for systems of differential equations. I need the time dependency.

Cheers, Franziska

请先登录，再进行评论。

Answer 2

Franziska 2013-3-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/64417-how-to-solve-a-system-of-nonlinear-2nd-order-differential-equations#answer_77087

My code is very sensitive to the choice of initial velocity values. I have some questions here:

1) Is there any possibility to search for appropriate initial velocity values similar to decic for ode15i? 2) Does it make sense if I try with ode15i and decic? Are there any limitations due to singular and time dependent mass matrix? I couldn't find a 2nd order problem similar to my challenge but coded in ode15i which might help me.

Thanks.

Franziska

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Alessandro Antonini 2013-6-1

在 MATLAB Online 中打开

I need to solve a system of 3 equations in the variable x1,x2,x3, I do not know how write the ode function that takes into account a term of a second order derivative of x2 in equation 1. I have a system like that:

ddx1=F1(t)-B1*dx1-M3*ddx3-B3*dx3-M2*ddx2-B2*dx2
ddx2=F2(t)-B2*dx2-K2*dx-M1*ddx1-B1*dx1-M2*ddx2-B2*dx2
ddx2=F3(t)-B3*dx3-K3*dx3-M1*ddx1-B1*dx1-K1*dx1-M2*ddx2-B2*dx2

I do not know how write in the ode function for this system. Can you please explain o write an example of the ode function required to solve a non linear system like that? I would be greateful

Best regards Alessandro Antonini

请先登录，再进行评论。