| 查看: 1151 | 回复: 2 | |||
[交流]
【求助】刚性微分方程的边值问题
|
|
我是搞力学的，现在一个动力学方程，碰到问题了。希望和大家讨论。 14个一阶微分方程（已无量纲化）及边界条件 用矩阵表示为： d{S(t)}/dt=[D(t)]*{S(t)} (1） 边界条件为：在两个端点处0和L，14个未知函数的前7个分量在为0，后7个分量是任意的。 这是一个变系数的线性微分方程组，微分矩阵[D(t)]中，含有一个很大的元素，[D(t)]的条件数非常大，用常用的计算指数矩阵的方法都不能算。 我的想法是：先把边值问题，用打靶法化为初值问题，然后可以用一些计算刚性方程的如：隐式龙格库塔法等来求解初值问题。这些数值算法，都是从初始点出发，逐一计算的。但是，我的边界条件，在初始点0处，信息不完整，只知道14个未知函数的前7个分量在为0，这样的边界条件如何解决？这样的初始点如何逐一向下计算？ 恳请各位不吝赐教。非常感谢 [ Last edited by haohao517112 on 2010-12-23 at 12:40 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有23人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有22人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stOpt求解常微分方程边值问题
已经有5人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 微分方程的初值问题
已经有7人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】求助一个一阶常微分方程的初值问题
已经有16人回复
- 【求助】常微分方程组求解中系数与某变量值关联的问题
已经有12人回复
- 【求助】解常微分方程的问题
已经有3人回复
- 求张文生《科学计算中的偏微分方程有限差分法》和胡建伟，汤怀民《微分方程数值解法》
已经有3人回复
- 【分享】PDECOL-经典数值求解非线性偏微分方程组源程序
已经有142人回复
2楼2010-12-24 08:37:26
3楼2010-12-25 00:40:31
