| 查看: 1489 | 回复: 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 |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有210人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 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