y*2 发帖数: 81 | 1 有如下常微分方程组
F''' +3F*F''-2(F')^2+R=0
R''+3*P*F*R'=0
其中,F,R为变量t的函数,,F'''表示F的三阶导数,P是变化的常数,变化范围从1到
100,边界条件为:
t=0,F=F'=0,R=1
t->无穷,F'->无穷,R=0
我不懂怎么转换无穷的边界条件到有界,
请高人指点
多谢多谢~ | E*V 发帖数: 17544 | 2 infinity can be approximated by 10, 100, 1000, 10000,,,,,
you got the idea.
anyway, use mathematica or maple. I think they support infinity
【在 y*2 的大作中提到】 : 有如下常微分方程组 : F''' +3F*F''-2(F')^2+R=0 : R''+3*P*F*R'=0 : 其中,F,R为变量t的函数,,F'''表示F的三阶导数,P是变化的常数,变化范围从1到 : 100,边界条件为: : t=0,F=F'=0,R=1 : t->无穷,F'->无穷,R=0 : 我不懂怎么转换无穷的边界条件到有界, : 请高人指点 : 多谢多谢~
| y*2 发帖数: 81 | 3 谢谢指点
但是因为在该方程组中仍然有一个P是变长的
我担心用10,100……这样去逼近会不会有影响?
另外
这个方程组有可能转化成Canonical形式吗?
【在 E*V 的大作中提到】 : infinity can be approximated by 10, 100, 1000, 10000,,,,, : you got the idea. : anyway, use mathematica or maple. I think they support infinity
| E*V 发帖数: 17544 | 4 sorry,我根本没有看你的问题。。。
呵呵
你看看如果boundary变化后,结果变化部大。就差不多了
【在 y*2 的大作中提到】 : 谢谢指点 : 但是因为在该方程组中仍然有一个P是变长的 : 我担心用10,100……这样去逼近会不会有影响? : 另外 : 这个方程组有可能转化成Canonical形式吗?
| y*2 发帖数: 81 | 5 嗯
是个好主意
我明白啦~
谢谢
但是这种混合着二个函数的常微分方程组怎么解啊?
再次感谢~
【在 E*V 的大作中提到】 : sorry,我根本没有看你的问题。。。 : 呵呵 : 你看看如果boundary变化后,结果变化部大。就差不多了
| l******n 发帖数: 9344 | 6 不能这么弄吧
好歹也得写成标准形式,看看性态吧
这和初值也有关系
【在 E*V 的大作中提到】 : sorry,我根本没有看你的问题。。。 : 呵呵 : 你看看如果boundary变化后,结果变化部大。就差不多了
| y*2 发帖数: 81 | 7 能说得稍微具体点吗?
拜谢~~
【在 l******n 的大作中提到】 : 不能这么弄吧 : 好歹也得写成标准形式,看看性态吧 : 这和初值也有关系
|
|