How to solve Coupled Differential Equations
조회 수: 52 (최근 30일)
이전 댓글 표시
Hi
how to solve that set of equations
y'-x'=2-x
2x'-y'=3+2y
댓글 수: 0
답변 (2개)
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
ode2 = 2-Dy - Dy == 3 + 2*y
S = dsolve(ode1, ode2, x(0)==x0, y(0)==y0)
x(t) = simplify(S.x, 500)
y(t) = simplify(S.y, 500)
.
댓글 수: 0
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
2021년 5월 9일
syms xprime yprime x y
eqn = [yprime - xprime == 2-x, 2*xprime - yprime == 3 + 2*y]
sol = solve(eqn, [xprime, yprime])
sol.xprime
sol.yprime
Hmmm, where did I go wrong?
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
참고 항목
카테고리
Help Center 및 File Exchange에서 Numerical Integration and Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!