24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 936  |  回复: 2

haohao517112

新虫 (初入文坛)


[交流] 【求助】刚性微分方程的边值问题

我是搞力学的,现在一个动力学方程,碰到问题了。希望和大家讨论。
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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jfili

金虫 (正式写手)


S(t),D(t)都是什么阶数的矩阵?中间的乘法是矩阵的普通乘法吗?
D(t)没什么要求吗?

比如说,我设D(t)=0,在你的边界条件下方程的解也不适定啊
2楼2010-12-24 08:37:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haohao517112

新虫 (初入文坛)


关于方程的具体说明

左端是一个14*1向量的导数,右端是一个14*14的矩阵,乘以这个14*1向量。就是普通的矩阵乘以向量,就是普通的矩阵乘法。
微分矩阵没有什么特点,不是对称的,也不是正规的。只是矩阵中有一个很大的数,数量级远远超出了其他的元素。
3楼2010-12-25 00:40:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 haohao517112 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700,311,求调剂 +14 冬十三 2026-04-04 15/750 2026-04-06 01:57 by BruceLiu320
[考研] 一志愿郑州大学085600求调剂 +17 吃的不少 2026-04-05 20/1000 2026-04-06 00:32 by T可可西里T
[考研] 材料求调剂 +11 一样YWY 2026-04-05 11/550 2026-04-05 23:36 by 来看流星雨10
[考研] 机械专硕274求调剂,不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] 272分求调剂 +4 wangyile2233 2026-04-02 4/200 2026-04-05 22:21 by 286640313
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 385分 生物学(071000)求调剂 +11 qf626 2026-04-01 11/550 2026-04-05 17:35 by Ecowxq666!
[考研] 283求调剂 +10 A child 2026-04-04 10/500 2026-04-05 08:22 by qlm5820
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 求调剂!生物与医药专硕 +4 逆转陆先生 2026-04-01 5/250 2026-04-03 08:33 by Jaylen.
[考研] 求调剂 +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[考研] 环境工程调剂 +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 313求调剂 +6 卖个关子吧 2026-03-31 6/300 2026-03-31 10:58 by Jaylen.
信息提示
请填处理意见