24小时热门版块排行榜    

CyRhmU.jpeg
查看: 6879  |  回复: 24
本帖产生 1 个 数学EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zgchen9

金虫 (小有名气)

[求助] 请问如何求解二元一阶微分方程组

各位朋友新年好,我需要解一个二元一阶微分方程组,但是本人数学水平有限,特请交各位朋友。
A,B,C,D,E,F,K为常数,x 和y为t 的函数。dx/dt和dy/dt为导数,二元一阶微分方程组如下:
dx/dt=Ax+By+C
dy/dt=Dx+Ey+F
边界条件为t=0时,x=y=K.


请问如何得到x 和y. 谢谢。

方程组
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

一份耕耘,一份收获
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilac_c

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
zgchen9(金币+2): 有帮助 2012-01-03 19:56:55
我们一般不直接求出解析解,太难了.

 都采用数值格式求解

离散,.......龙格库塔
我生活在一个经常爆发地震的年代
4楼2012-01-03 09:49:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilac_c

至尊木虫 (知名作家)

【答案】应助回帖


soliton923(金币+1): 谢谢参与讨论~~~ 2012-01-04 22:20:52
dx/dt=Ax+By+C
dy/dt=Dx+Ey+F

基于 crank-Nicloson格式的求解程序编写思路

第一步
显示求解出 当前的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)
当两次求解误差在设置的误差范围内,结束迭代过程,否则,重新回到校正.

.................
我生活在一个经常爆发地震的年代
6楼2012-01-03 21:12:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zgchen9 的主题更新
信息提示
请填处理意见