| 查看: 370 | 回复: 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。 |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有193人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2015-10-20 22:25:46
水中的梦88
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 2378.9
- 散金: 26
- 红花: 5
- 帖子: 323
- 在线: 195.7小时
- 虫号: 1308847
- 注册: 2011-05-28
- 性别: MM
- 专业: 结构工程
|
谢谢您的提醒，可能是调试的过程中多的，没有发现，已经修改，如下： 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