| 查看: 1837 | 回复: 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论文润色/翻译怎么收费?
已经有149人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有11人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神帮忙看一下MATLAB求解微分方程的问题？
已经有3人回复
- 求高手帮我解一阶常微分方程组，非常感谢！！
已经有5人回复
- 求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~
已经有4人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 求组！求解一个微分方程组的符号解，并画出图像！！！
已经有7人回复
- 帮忙求一个微分方程的数值解或者 真实解
已经有6人回复
- 一阶常微分方程组，matlab进行数值求解
已经有3人回复
- MATLAB求解代数方程组问题，请帮忙
已经有9人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 一个微分方程组，求解
已经有18人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 求解一个含字母系数的微分方程
已经有5人回复
- 谁那里可以运行maple，帮我解一个微分方程组
已经有3人回复
- 【求助】20金币求一个微积分方程
已经有15人回复
- 【求助】求解以二阶微分方程，赠金币！！
已经有13人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
5楼2014-09-27 17:56:06
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
priest
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1743.8
- 散金: 20
- 红花: 2
- 帖子: 658
- 在线: 463.7小时
- 虫号: 442901
- 注册: 2007-10-27
- 性别: GG
- 专业: 催化化学
|
(∂(Fx(14)x(1)))/∂t=(-3x(10)-x(8))*ρ_c ∂(Fx(14)x(2))/∂t=(2x(10)-x(9))*ρ_c x(9)=k_K x 〖(12)sinh〗(((x(5)-x(7)))⁄〖2v〗_T ) x(8)=k_Na x(13)x(11)sinh〖(((x(5)-x(6)))⁄〖2v〗_T ) 〗 x(10)=k_NaK tanh(((x(5)+2x(7)-3x(6)-v_ATP ))⁄(2v_T )) (dx(12))/dt=1/τ_K cosh(2 ((x(5)-v_x ))⁄v_T ) {0.5*[1+tanh(2 ((x(5)-v_x ))⁄v_T ) ]-x(12)} (dx(13))/dt=1/τ_Na cosh(2 ((x(5)-h))⁄v_T ) {0.5*[1+tanh(2 ((x(5)-v_h ))⁄v_T ) ]-x(13)} x(11)=0.5*(1+tanh(2(x(5)-v_m )/v_T ) ) x(5)=v_T ln〖(x(3))⁄(x(1))〗 x(6)=v_T ln〖(x(4))⁄(x(2))〗 x(5)=RT/F ln((P_K x(4)+P_Na x(2)+P_Cl 〖Cl〗_i)/(P_K x(3)+P_Na x(1)+P_Cl 〖Cl〗_e )) (∂((1-x(14))x(4)))/∂t=∇D_2 ∇x(4)-∂(x(14)x(2))/∂t (∂((1-x(14))x(3)))/∂t=∇D_1 ∇x(3)-∂(x(14)x(1))/∂t (∂x(14))/∂t=L_p [(x(1)+x(2))-(x(3)+x(4)+Ci)] 全部14个方程如上，x是状态变量，其他都是参数（常数）。尽管x（3）和x（4）随空间变换，但可以通过有限差分进行离散。多谢啦！ |
6楼2014-09-28 06:43:48
20