[求助]ODE数值解不converge - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Computation - 科学计算

[求助]ODE数值解不converge

[求助]ODE数值解不converge# Computation - 科学计算

h*o2004-10-26 07:10

1 楼

我要解一个ODE, 形式比较复杂点，假设x是自变量，y是函数，y', y", y"'是y对
x的1, 2, 3阶导数，ODE差不多是这样的
F(y",y',y,x)*y" = G(y"',y",y',y,x)*y'
F, G这两个函数里面，y的各阶导数加减乘除纽在一块的(非线性).
为了解它，我把y"'从G里面抽出来，变成
y"' = Y(y",y',y,x)
当然，Y()里面的y",y',y和x都是非线性的纽在一块的，然后用
ODE solver去解，结果总出不converge的问题。
请问，什么样的ODE, 才可以用数值来解呢？

c*e2004-10-26 07:10

2 楼

我一般都写成一阶形式,然后用matlab的solver解一下
注意是不是stiff的.

【在 h***o 的大作中提到】

: 我要解一个ODE, 形式比较复杂点，假设x是自变量，y是函数，y', y", y"'是y对
: x的1, 2, 3阶导数，ODE差不多是这样的
: F(y",y',y,x)*y" = G(y"',y",y',y,x)*y'
: F, G这两个函数里面，y的各阶导数加减乘除纽在一块的(非线性).
: 为了解它，我把y"'从G里面抽出来，变成
: y"' = Y(y",y',y,x)
: 当然，Y()里面的y",y',y和x都是非线性的纽在一块的，然后用
: ODE solver去解，结果总出不converge的问题。
: 请问，什么样的ODE, 才可以用数值来解呢？