MATLAB Answers

How to solve Coupled Differential Equations

조회 수: 2(최근 30일)
Tomasz
Tomasz 2021년 5월 9일
댓글: Walter Roberson 2021년 5월 9일
Hi
how to solve that set of equations
y'-x'=2-x
2x'-y'=3+2y

답변(2개)

Star Strider
Star Strider 2021년 5월 9일
Try this —
syms x(t) y(t) x0 y0
Dx = diff(x);
Dy = diff(y);
ode1 = Dy-Dx == 2 - x
ode1(t) = 
ode2 = 2-Dy - Dy == 3 + 2*y
ode2(t) = 
S = dsolve(ode1, ode2, x(0)==x0, y(0)==y0)
S = struct with fields:
y: [1×1 sym] x: [1×1 sym]
x(t) = simplify(S.x, 500)
x(t) = 
y(t) = simplify(S.y, 500)
y(t) = 
.

Walter Roberson
Walter Roberson 2021년 5월 9일
yprime - xprime == 2-x, 2*xprime - yprime == 3 + 2*y
2*yprime - 2*xprime == 4-2*x, 2*xprime - yprime == 3 + 2*y
2*yprime - 2*xprime + 2*xprime - yprime == 4-2*x + 3 + 2*y
yprime == 7 - 2*x + 2*y
7 - 2*x + 2*y - xprime == 2*x
7 - 2*x + 2*y - 2*x == xprime
xprime = 7 - 4*x + 2*y
So...
function dxy = odefun(t, xy)
dxy = [7 - 4*xy(1) + 2*xy(2); 7 - 2*xy(1) + 2*xy(2)];
end
  댓글 수: 4
Walter Roberson
Walter Roberson 2021년 5월 9일
I did make a mistake, but it was the step after you indicated. I copied as 2*x instead of 2-x . The corrected version is
yprime - xprime == 2-x, 2*xprime - yprime == 3 + 2*y
2*yprime - 2*xprime == 4-2*x, 2*xprime - yprime == 3 + 2*y
2*yprime - 2*xprime + 2*xprime - yprime == 4-2*x + 3 + 2*y
yprime == 7 - 2*x + 2*y
7 - 2*x + 2*y - xprime == 2 - x %corrected
7 - 2*x + 2*y - (2 - x) == xprime
xprime = 5 - x + 2*y
So...
function dxy = odefun(t, xy)
dxy = [5 - xy(1) + 2*xy(2); 7 - 2*xy(1) + 2*xy(2)];
end

댓글을 달려면 로그인하십시오.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by