| 查看: 1556 | 回复: 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 |
回复此楼» 猜你喜欢
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有221人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有41人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有1人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
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
50