24小时热门版块排行榜    

查看: 327  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

水中的梦88

木虫 (正式写手)

[求助] 1stopt每次运行结果不一样,求虫友们帮忙看看,谢谢! 已有1人参与

代码如下:
Parameter B,C,De0,kr;
Variables t,y;
OdeFunction y'=3*0.5*(-3.07+35.549*y)/(3.15*30*0.35*1e-6)/((1/(B/y^1.5+C*y^3)+30*1e-6/(De0*ln(1/y))*((1-y)^(-1/3)-1)+1/kr*(1-y)^(-2/3));
Data;

         0    0.3535
    7.0000    0.3550
   14.0000    0.3726
   28.0000    0.3776
   56.0000    0.3905
   90.0000    0.4033
  180.0000    0.4201

没有对代码中的数字计算部分进行化简,因为还有其他方程,个别数字需要调整。
采用5.5试用版,有得时候能运行出结果,运行出的结果每次都不一样,有得时候运行出来的参数是0。
求大神帮忙看看,主要的问题在哪儿,谢谢!
如果需要对参数进行约束的话,约束条件是都大于0。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水中的梦88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2015-10-20 22:25:46
公式左右括号都不匹配,不知道你怎么能计算的

谢谢您的提醒,可能是调试的过程中多的,没有发现,已经修改,如下:
Parameter B,C,De0,kr;
Variables t,y;
OdeFunction y'=3*0.5*(-3.07+35.549*y)/(3.15*30*0.35*1e-6)/(1/(B/y^1.5+C*y^3)+30*1e-6/(De0*ln(1/y))*((1-y)^(-1/3)-1)+1/kr*(1-y)^(-2/3));
Data;

         0    0.3535
    7.0000    0.3550
   14.0000    0.3726
   28.0000    0.3776
   56.0000    0.3905
   90.0000    0.4033
  180.0000    0.4201
3楼2015-10-21 08:49:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
公式左右括号都不匹配,不知道你怎么能计算的
2楼2015-10-20 22:25:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854AI CV方向招收调剂 +3 章小鱼567 2026-03-23 3/150 2026-03-24 20:25 by 汪!?!
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-24 3/150 2026-03-24 19:17 by pswait
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 303求调剂 +4 元夕元 2026-03-20 4/200 2026-03-23 19:00 by macy2011
[考研] 352求调剂 +3 大米饭! 2026-03-22 3/150 2026-03-22 23:28 by king123!
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
信息提示
请填处理意见