应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】刚性微分方程的边值问题
31/1返回列表
查看: 1000  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-12-23 12:36:31
已阅   回复此楼   关注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 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭