s*********n 发帖数: 18 | 1 方程中有一函数依赖于方程的解。
我将这两个联立的二阶常微分方程化为了四个一阶常微分方程,然后用RUNGE-KUTTA 方法
进行
积分,但是解在几个循环后变趋向与无穷。SHOOTING METHOD 好象作用不明星。
谢谢。 |
s*********n 发帖数: 18 | 2 谢谢howell (Beaver)。
我没说清楚,对不起。这是两个联立的ODE。我重新叙述一遍。
方程中有一函数依赖于方程的解。
我将这两个联立的二阶ODE微分方程化为了四个一阶微分方程,然后用RUNGE-KUTTA 方法
进行积分,但是解在几个循环后变趋向与无穷。SHOOTING METHOD 因为需要同时调整两个
斜率,算了很长时间也没结果。
谢谢。
进行 |
s*********n 发帖数: 18 | 3 谢谢muying (心情娃娃∞放飞自由)。
试过MATLAB的ODE SOLVER:ODE45,ODE15S,不行。
进行 |
s*i 发帖数: 5025 | 4 如果不复杂,能帖方程吗?
法
个
法
【在 s*********n 的大作中提到】 : 谢谢howell (Beaver)。 : 我没说清楚,对不起。这是两个联立的ODE。我重新叙述一遍。 : 方程中有一函数依赖于方程的解。 : 我将这两个联立的二阶ODE微分方程化为了四个一阶微分方程,然后用RUNGE-KUTTA 方法 : 进行积分,但是解在几个循环后变趋向与无穷。SHOOTING METHOD 因为需要同时调整两个 : 斜率,算了很长时间也没结果。 : 谢谢。 : : 进行
|
c*******e 发帖数: 8624 | 5 奇怪,我还没见过用stiff solver解不出来的情况
不行就调一调收敛判定什么的
【在 s*********n 的大作中提到】 : 谢谢muying (心情娃娃∞放飞自由)。 : 试过MATLAB的ODE SOLVER:ODE45,ODE15S,不行。 : : 进行
|
m****g 发帖数: 984 | 6 换换初值?
【在 s*********n 的大作中提到】 : 谢谢muying (心情娃娃∞放飞自由)。 : 试过MATLAB的ODE SOLVER:ODE45,ODE15S,不行。 : : 进行
|
m****g 发帖数: 984 | 7 spt
【在 s*i 的大作中提到】 : 如果不复杂,能帖方程吗? : : 法 : 个 : 法
|
s*********n 发帖数: 18 | 8 我在看是否是我当时使用MATLAB SOLVER 的时候的理解不对,再解解看。 请问如何在这
里输入公式啊。
谢谢诸位。
【在 c*******e 的大作中提到】 : 奇怪,我还没见过用stiff solver解不出来的情况 : 不行就调一调收敛判定什么的
|
s*********n 发帖数: 18 | 9 我试过改初值。
刚才我又用MATLAB的ODE SOLVER的ODE15S试了,从时间上来看,它在某一个时间点不停的
循环。从我自己编的程序计算结果发现,方程的解使得内嵌在方程中的函数值(实际是个
压力函数值)变成了复数,方程的解趋于无穷。
【在 m****g 的大作中提到】 : 换换初值?
|
h****l 发帖数: 7290 | 10 还不是跟原来一个意思,不就是迭代解偏微分方程么?
多数时候是会发散的,这个一点都不奇怪,因为显式
的迭代方法收敛条件比较苛刻。
【在 s*********n 的大作中提到】 : 谢谢howell (Beaver)。 : 我没说清楚,对不起。这是两个联立的ODE。我重新叙述一遍。 : 方程中有一函数依赖于方程的解。 : 我将这两个联立的二阶ODE微分方程化为了四个一阶微分方程,然后用RUNGE-KUTTA 方法 : 进行积分,但是解在几个循环后变趋向与无穷。SHOOTING METHOD 因为需要同时调整两个 : 斜率,算了很长时间也没结果。 : 谢谢。 : : 进行
|
h****l 发帖数: 7290 | 11 co-spt
【在 m****g 的大作中提到】 : spt
|
h****l 发帖数: 7290 | 12 可以定义一个符号代表偏微分的符号,就不成问题了。
【在 s*********n 的大作中提到】 : 我在看是否是我当时使用MATLAB SOLVER 的时候的理解不对,再解解看。 请问如何在这 : 里输入公式啊。 : 谢谢诸位。
|
h****l 发帖数: 7290 | 13 就是sqrt domain error是吧,可是如果你用数值方法解的话,
应该有边界条件吧,那边界条件是怎么给的呢?
【在 s*********n 的大作中提到】 : 我试过改初值。 : 刚才我又用MATLAB的ODE SOLVER的ODE15S试了,从时间上来看,它在某一个时间点不停的 : 循环。从我自己编的程序计算结果发现,方程的解使得内嵌在方程中的函数值(实际是个 : 压力函数值)变成了复数,方程的解趋于无穷。
|