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.
|
|
