| 查看: 1549 | 回复: 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))’。多谢大家的帮助！ |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有74人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神帮忙看一下MATLAB求解微分方程的问题？
已经有3人回复
- 求高手帮我解一阶常微分方程组，非常感谢！！
已经有5人回复
- 求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~
已经有4人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 求组！求解一个微分方程组的符号解，并画出图像！！！
已经有7人回复
- 帮忙求一个微分方程的数值解或者 真实解
已经有6人回复
- 一阶常微分方程组，matlab进行数值求解
已经有3人回复
- MATLAB求解代数方程组问题，请帮忙
已经有9人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- 一个微分方程组，求解
已经有18人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 求解一个含字母系数的微分方程
已经有5人回复
- 谁那里可以运行maple，帮我解一个微分方程组
已经有3人回复
- 【求助】20金币求一个微积分方程
已经有15人回复
- 【求助】求解以二阶微分方程，赠金币！！
已经有13人回复
priest
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1743.8
- 散金: 20
- 红花: 2
- 帖子: 658
- 在线: 463.7小时
- 虫号: 442901
- 注册: 2007-10-27
- 性别: GG
- 专业: 催化化学
2楼2014-09-27 06:47:37
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 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
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
5楼2014-09-27 17:56:06
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
7楼2014-10-05 23:21:35
zhemeizimei
铁虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 67.8
- 红花: 1
- 帖子: 52
- 在线: 21.4小时
- 虫号: 2206356
- 注册: 2012-12-26
- 性别: GG
- 专业: 金属结构材料
8楼2014-10-06 10:55:07
