| 查看: 1833 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
priest木虫 (正式写手)
|
[求助]
帮忙：求解一组微分代数方程，奖励100金币 已有1人参与
|
||
| 请各位帮忙求解一组微分代数方程，14个方程，其中7个微分方程，7个代数方程。我的思路是这样的：把代数方程对时间求导，转化成微分方程，这样方程组写成 M(t,x)*dx/dt = F(t,x)，可以用matlab ode15s求解，关键是写出正确的mass matrix M。我遇到的问题就是：（1）matlab会显示 This DAE appears to be of index greater than 1,但是我仔细检查了代数方程的微分形式，例如，dx1-2*dx3 = x1+4*x2,就是所有代数方程的变量对时间的一阶导数＝f(x, dx/dt), 全部都是一阶导数。所以，我认为这组微分代数方程的index为1，不应该有这个错误。（2）如果x1为状态变量， M(1,1) = c1/x1, matlab 代码里应该如何表示，x1作为数组还是标量，是要用c1./x1 还是c1/x1？或者，我可以在方程两边同时乘以x1?（3） 如果x1=c1*sinh(x2-x3), c1=exp(-c2*t), c1只是作为输入函数，这种情况下，x1对时间进行微分，c1也要对对时间求导么，如果这样，c1也变成了一个状态变量。还是把c1当作常熟处理，dx1/dt=c1*(sinh(x2-x3))’。多谢大家的帮助！ |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有152人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有10人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神帮忙看一下MATLAB求解微分方程的问题？
已经有3人回复
- 求高手帮我解一阶常微分方程组，非常感谢！！
已经有5人回复
- 求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~
已经有4人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 求组！求解一个微分方程组的符号解，并画出图像！！！
已经有7人回复
- 帮忙求一个微分方程的数值解或者 真实解
已经有6人回复
- 一阶常微分方程组，matlab进行数值求解
已经有3人回复
- MATLAB求解代数方程组问题，请帮忙
已经有9人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 一个微分方程组，求解
已经有18人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 求解一个含字母系数的微分方程
已经有5人回复
- 谁那里可以运行maple，帮我解一个微分方程组
已经有3人回复
- 【求助】20金币求一个微积分方程
已经有15人回复
- 【求助】求解以二阶微分方程，赠金币！！
已经有13人回复
zhemeizimei
铁虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 67.8
- 红花: 1
- 帖子: 52
- 在线: 21.4小时
- 虫号: 2206356
- 注册: 2012-12-26
- 性别: GG
- 专业: 金属结构材料
8楼2014-10-06 10:55:07
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
3楼2014-09-27 08:16:11
priest
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1743.8
- 散金: 20
- 红花: 2
- 帖子: 658
- 在线: 463.7小时
- 虫号: 442901
- 注册: 2007-10-27
- 性别: GG
- 专业: 催化化学
4楼2014-09-27 17:49:03
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
5楼2014-09-27 17:56:06
40