应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt每次运行结果不一样，求虫友们帮忙看看，谢谢！
31/1返回列表
查看: 304  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

水中的梦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。
回复此楼

» 猜你喜欢
1楼 2015-10-20 19:49:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
公式左右括号都不匹配，不知道你怎么能计算的
赞一下(2人) 回复此楼
2楼2015-10-20 22:25:46
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 水中的梦88 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭