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试了，从时间上来看，它在某一个时间点不停的
: 循环。从我自己编的程序计算结果发现，方程的解使得内嵌在方程中的函数值(实际是个
: 压力函数值)变成了复数，方程的解趋于无穷。
|
