| 查看: 321 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 15322002 的 6 个金币 | |||
15322002银虫 (小有名气)
|
[交流]
【求助】MATLAB求解一阶常微分方程
|
||
|
求解一个只有一个变量的常微分方程 方程： dX=zeros(3,1); dX=H/ug*(1-exp(-(3*(1-ep)/r0/(1/Dg+r0*((1-X)^(-1/3)-1)/DH2+1/(RsH2*(1+1/KH2))*(1-X)^(-2/3)))))*b*epb*(PH2-PH2/(1+KH2))/R/Tg/(1-epb)/rB... +H/ug*(1-exp(-(3*(1-ep)/r0/(1/Dg+r0*((1-X)^(-1/3)-1)/DCO+1/(RsCO*(1+1/KCO))*(1-X)^(-2/3)))))*b*epb*(PCO-PCO/(1+KCO))/R/Tg/(1-epb)/rB; 变量是X，其余是常数 t,XB]=ode45(@fun_for_dc,[0:100:3000],[0],options) 总是报错：In an assignment A(I) = B, the number of elements in B and I must be the same. 采用点乘和点除情况依旧。 不知是何原因？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
feichengxyz
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 8535.5
- 散金: 60
- 红花: 4
- 帖子: 1162
- 在线: 908.4小时
- 虫号: 425858
- 注册: 2007-07-28
- 性别: GG
- 专业: 催化化学
2楼2009-09-08 13:09:41
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
★ ★ ★ ★
loovfnd(金币+2,VIP+0):xiexie 9-8 17:05
15322002(金币+2,VIP+0):问题已经解决，全局变量的原因 9-9 11:44
loovfnd(金币+2,VIP+0):xiexie 9-8 17:05
15322002(金币+2,VIP+0):问题已经解决，全局变量的原因 9-9 11:44
|
问题应该出在赋值语句上。检查 dX=H/ug*(1-exp(-(3*(1-ep)/r0/(1/Dg+r0*((1-X)^(-1/3)-1)/DH2+1/(RsH2*(1+1/KH2))*(1-X)^(-2/3)))))*b*epb*(PH2-PH2/(1+KH2))/R/Tg/(1-epb)/rB... +H/ug*(1-exp(-(3*(1-ep)/r0/(1/Dg+r0*((1-X)^(-1/3)-1)/DCO+1/(RsCO*(1+1/KCO))*(1-X)^(-2/3)))))*b*epb*(PCO-PCO/(1+KCO))/R/Tg/(1-epb)/rB; 中各个变量是否都是标量。 |
3楼2009-09-08 15:54:14
30