| 5 | 1/1 | 返回列表 |
| 查看: 7219 | 回復(fù): 24 | ||
| 本帖產(chǎn)生 1 個 數(shù)學(xué)EPI ,點擊這里進(jìn)行查看 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
zgchen9金蟲 (小有名氣)
|
[求助]
請問如何求解二元一階微分方程組
|
|
|
各位朋友新年好,我需要解一個二元一階微分方程組,但是本人數(shù)學(xué)水平有限,特請交各位朋友。 A,B,C,D,E,F,K為常數(shù),x 和y為t 的函數(shù)。dx/dt和dy/dt為導(dǎo)數(shù),二元一階微分方程組如下: dx/dt=Ax+By+C dy/dt=Dx+Ey+F 邊界條件為t=0時,x=y=K. 請問如何得到x 和y. 謝謝。 方程組 |

鐵桿木蟲 (職業(yè)作家)
數(shù)學(xué)村村長

|
用mathematica軟件可以求解: DSolve[{x'[t] == A x[t] + B y[t] + C, y'[t] == D x[t] + E y[t] + F, x[0] == K, y[0] == K}, {x, y}, t] // Simplify 答案是: x -> Function[{t}, (2 E^(-(1/ 2) (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) (-2 B C D E^( 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + 2 B C D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + 2 B C D E^( Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - A C E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + C E^(2 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + A C E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - C E^(2 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - 2 B C D E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + A C E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - C E^(2 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - A C E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + C E^(2 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + C E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] + C E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[A^2 + 4 B D - 2 A E + E^2] - C E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] - C E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] + A B E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - A B E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - A B E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - B E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + B E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + A B E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + B E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - B E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - B E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F + B E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F - B E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F + B E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F + A B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + 2 B^2 D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - A B D E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - 2 B^2 D E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + A^2 E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + 2 A B E^( 1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + B D E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - A E^(2 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - A^2 E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - 2 A B E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - B D E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + A E^(2 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K + B D E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K - A E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K - A E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K))/(Sqrt[ A^2 + 4 B D - 2 A E + E^2] (-A - E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) (A + E + Sqrt[ A^2 + 4 B D - 2 A E + E^2]))], y -> Function[{t}, -(2 E^(-(1/ 2) (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) (-A C D E^( 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + A C D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + A C D E^( Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + C D E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - C D E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - A C D E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) - C D E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + C D E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) + C D E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[A^2 + 4 B D - 2 A E + E^2] - C D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] + C D E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] - C D E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] + A^2 E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + 2 B D E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - A^2 E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - 2 B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - A^2 E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - 2 B D E^( Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + A E^(1 + Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - A E^(1 + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + A^2 E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + 2 B D E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - A E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F + A E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) F - A E^(1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F + A E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F - A E^(Sqrt[A^2 + 4 B D - 2 A E + E^2] t + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F + A E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] F + A B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - 2 B D^2 E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - A B D E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + 2 B D^2 E^( 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + A^2 E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - 2 A D E^( 1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + B D E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - A E^(2 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - A^2 E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + 2 A D E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - B D E^( 1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K + A E^(2 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) K - B D E^((A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K - B D E^(1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K + A E^(1 + 1/2 (A + E - Sqrt[A^2 + 4 B D - 2 A E + E^2]) t + 1/2 (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K + A E^(1 + (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) t) Sqrt[ A^2 + 4 B D - 2 A E + E^2] K))/(Sqrt[ A^2 + 4 B D - 2 A E + E^2] (-A - E + Sqrt[A^2 + 4 B D - 2 A E + E^2]) (A + E + Sqrt[A^2 + 4 B D - 2 A E + E^2]))]}} |

至尊木蟲 (知名作家)

至尊木蟲 (知名作家)
|
dx/dt=Ax+By+C dy/dt=Dx+Ey+F 基于 crank-Nicloson格式的求解程序編寫思路 第一步 顯示求解出 當(dāng)前的x(n+1)與y(n+1) 以第一個方程為例: x(n+1)=x(n)+dt*( A*x(n)+B*y(n)+C) 第二部校正開始 由于用顯格式求解誤差會越來越大,故此,要用 x(n+1)=x(n)+dt*( A/2*(x(n)+x(n+1))+B*(y(n)+y(n+1))/2. +c) 當(dāng)兩次求解誤差在設(shè)置的誤差范圍內(nèi),結(jié)束迭代過程,否則,重新回到校正. ................. |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +4 | 孅華 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 293求調(diào)劑 +12 | zjl的號 2026-03-16 | 17/850 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|