| 查看: 416 | 回复: 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。 |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有105人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
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