| 查看: 325 | 回复: 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。 |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有222人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
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
