ode45

Question

Nani C 2012-5-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/37626-ode45

Hi, Having solved a second order equation of motion using ode45 function i wonder how could i modify the function to solve a whole system of equations in matrix form [A]{xdoubledot}+[B]{xdot}+[c]{x}={p(t)}, instead of solving individual equations for x vector variables. Would that be possible ?

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

请先登录，再进行评论。

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

Answer 1

Richard Brown 2012-5-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/37626-ode45#answer_46926

在 MATLAB Online 中打开

Yes, certainly possible. The basic approach would be to turn it into a system of first order odes of twice the dimension as follows (where x and xdot are vectors):

y1 = x
y2 = xdot

then, if A is invertible

y1' = y2
y2' = A \ (-B*y2 - c*y1 + p(t))

You'd then write a file to implement this in MATLAB like this

function dy = myode(t, y)
  n = numel(y)/2;
  y1 = y(1:n);
  y2 = y(n+1:2*n);
  dy = [y2; A \ (-B*y2 - c*y1 + p(t))
end

and solve it with ode45 (or whichever solver happens to be most suitable) as usual.

2 个评论
显示无隐藏无

Nani C 2012-5-7

t.y. Richard. it worked. It would be helpful if you can clarify on 'twice the dimension'. Any help links will also be useful.

Jan 2012-5-7

You can convert an ODE of degree 2 for an n-dimensional vector to an ODE of degree 1 for a 2n-dimensional vector. This means "twice the dimensions". In other words: Instead of calculating the 1st and 2nd derivative of the position, you calculate the 1st derivative of the position and the velocity.

请先登录，再进行评论。