non linear coupled diff equation
3 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
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 Kommentare
Jan
am 20 Jul. 2021
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.
Akzeptierte Antwort
Jan
am 20 Jul. 2021
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.
Do I oversee a point?
0 Kommentare
Weitere Antworten (0)
Siehe auch
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!