24小时热门版块排行榜    

查看: 307  |  回复: 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的回帖
信息提示
请填处理意见