版块导航: 正在加载中...

客户端APP下载

登录注册

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

24小时热门版块排行榜

haohao517112

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 7.5
帖子: 7
在线: 1.9小时
虫号: 1170510

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

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

» 猜你喜欢

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

1stOpt求解常微分方程边值问题已经有5人回复
matlab数值求解边界条件微分方程组已经有7人回复
微分方程的初值问题已经有7人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】求助一个一阶常微分方程的初值问题已经有16人回复
【求助】常微分方程组求解中系数与某变量值关联的问题已经有12人回复
【求助】解常微分方程的问题已经有3人回复
求张文生《科学计算中的偏微分方程有限差分法》和胡建伟，汤怀民《微分方程数值解法》已经有3人回复
【分享】PDECOL-经典数值求解非线性偏微分方程组源程序已经有142人回复

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

查看全部散金贴

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

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 7.5
帖子: 7
在线: 1.9小时
虫号: 1170510

关于方程的具体说明

左端是一个14*1向量的导数，右端是一个14*14的矩阵，乘以这个14*1向量。就是普通的矩阵乘以向量，就是普通的矩阵乘法。
微分矩阵没有什么特点，不是对称的，也不是正规的。只是矩阵中有一个很大的数，数量级远远超出了其他的元素。

3楼2010-12-25 00:40:31

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

相关版块跳转我要订阅楼主 haohao517112 的主题更新