| 查看: 1588 | 回复: 2 | ||
[求助]
能帮忙看一下吧，为什么用1stopt5.5试用版，点运行没反应啊 已有1人参与
|
|
能帮忙看一下吧，为什么用1stopt5.5试用版，点运行没反应啊 ConstStr k=t^b*a*0.473/(2.081-1.608*t); Parameters a,b; Variable t,r; ODEFunction r''=-r'*(1+10*k'/k)/10; Data; 0.228398899 8.32 0.242900416 13.47 0.255672483 15.64 0.267082228 17.69 0.277390305 18.64 0.286788667 18.96 0.295422884 19.35 0.303406043 20.17 0.310827786 21.69 |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有111人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有111人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有68人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
deepseawdd: 金币+5 2017-07-02 20:17:07
感谢参与，应助指数 +1
deepseawdd: 金币+5 2017-07-02 20:17:07
|
先把k的一阶导数求出来代替k': k'=((((b*(t^(b-1)))*a)*0.473)*(2.081-1.608*t)-(((t^b)*a)*0.473)*(-1.608))/((2.081-1.608*t)^2); k'由dk代替表示： a参数似乎无用了，dk/k把a约掉了，也即a可以为任意值。 均方差(RMSE): 0.474551966680738 残差平方和(SSE): 1.80159655264445 相关系数(R): 0.981564892013995 相关系数之平方(R^2): 0.963469637234446 修正R平方(Adj. R^2): 0.948857492128224 确定系数(DC): 0.961944496067107 F统计(F-Statistic): 68.2782010981496 参数 最佳估算 -------------------- ------------- a -67.2023509470802 b 7.18754485098689 r' 初值 406.49264251293 |
2楼2017-07-01 09:40:05
3楼2017-07-01 11:23:31
