h***o 发帖数: 539 | 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*******e 发帖数: 8624 | 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, 才可以用数值来解呢?
| r****y 发帖数: 1437 | 3
why not write one yourself? It would be more reliable.
【在 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, 才可以用数值来解呢?
| h***o 发帖数: 539 | 4 方程可能就是没解的,物理意义上有问题。改了下,现在行了。
【在 r****y 的大作中提到】 : : why not write one yourself? It would be more reliable.
|
|