non linear coupled diff equation

Question

Sarat 2021-7-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/882303-non-linear-coupled-diff-equation

评论： Sarat 2021-7-22

采纳的回答： Jan

dx/dt=-(r+r1)/L1*x - dy/dt*(L12/L1)

dy/dt= -y*r2 - dx/dt*(L12/L2

r=300 r1= 8.1 L1=13.8 L2 =0.012 L12=0.0013

Kindly guide in Matal core for above equation using ode23 as r2 is function of time

thanks in advance

Sarat Kumar Dash

4 个评论
显示 2更早的评论隐藏 2更早的评论

Jan 2021-7-20

在 MATLAB Online 中打开

As written in my answer:

dx/dt=   -A*x - B*dy/dt
dy/dt= -y*C - dx/dt*D

is equivalent to:

dx/dt= -A*x - B * (-y*C + A*x*D) / (1 -  B*D)
dy/dt = (-y*C + A*x*D) / (1 -  B*D)

This should be easy to be implemented, such that an ODE integrator can handle it. Please try it and ask again, if you have a specfic problem.

Sarat 2021-7-22

thanks for your reply

请先登录，再进行评论。

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

Answer 1

Jan 2021-7-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/882303-non-linear-coupled-diff-equation#answer_750383

在 MATLAB Online 中打开

What about simplifying the equations?

1. dx/dt = -(r+r1)/L1*x - dy/dt * (L12/L1)
2. dy/dt = -y*r2 - dx/dt * (L12/L2)
1. into 2.:
dy/dt = -y*r2 - (-(r+r1) / L1 * x - dy/dt * (L12/L1)) * (L12/L2)
      = -y*r2 + ((r+r1) / L1 * x + dy/dt * (L12/L1)) * (L12/L2)
      = -y*r2 + (r+r1) / L1 * x * (L12 / L2) + dy/dt * L12^2 / (L1*L2)
dy/dt - dy/dt * L12^2 / (L1*L2) = -y*r2 + (r+r1) / L1 * x * (L12 / L2)
dy/dt (1 - L12^2 / (L1*L2))     = -y*r2 + (r+r1) / L1 * x * (L12 / L2)
dy/dt = (-y*r2 + (r+r1) / L1 * x * (L12 / L2)) / (1 - L12^2 / (L1*L2));
Insert this in 1 again.