b******i 发帖数: 914 | 1 【 以下文字转载自 Mathematics 讨论区 】
发信人: bufangqi ( 不放弃), 信区: Mathematics
标 题: 请教大牛如何解常微分方程
发信站: BBS 未名空间站 (Wed Mar 16 21:50:55 2011, 美东)
请教各位数学大牛
用什么计算机软件可以模拟解出常微分方程组呢?
比如最简单的形式:
d x(t) / dt = -a * x(t)
但是我现在要解得方程组大概有几十个变量和几十个方程组,边界值如x(0)都已知。请
问有什么方便的
数值或者仿真软件可以计算出来任意t下x(t)的值呢?
谢谢了! |
b******i 发帖数: 914 | 2 PS. 我的这个问题似乎也不是线性的,但最高也是二次,比如存在
d x(t) / dt = -a * x(t) * y(t)
其中x(t), y(t)都未知 |
i**i 发帖数: 79 | 3 简单的我用mathcad, maple做过,matlab应该也可以。
你这个这么复杂,不确定哪个能做。 |
w***n 发帖数: 1084 | 4 离散化,然后搞成线性方程,然后biu一下就解出来了。 |
b******i 发帖数: 914 | 5 Hi,
我试了一下MATLAB,数值解用ode23, ode45等函数,挺快的
谢谢哈
【在 i**i 的大作中提到】 : 简单的我用mathcad, maple做过,matlab应该也可以。 : 你这个这么复杂,不确定哪个能做。
|
b******i 发帖数: 914 | 6 什么叫编程离散化?望指教 用C/C++等编程语言来做么?
赞 biu 一下就解出来呵呵
【在 w***n 的大作中提到】 : 离散化,然后搞成线性方程,然后biu一下就解出来了。
|
w***n 发帖数: 1084 | 7 t是时间吗?方程是time evolution functions吗?边界条件是初始条件吗?如果那样
的话,搞个简单的First order euler都可以。
如果不是,你这个不是线性的。不太好解。 |