| 查看: 549 | 回复: 2 | ||
yejianling新虫 (正式写手)
|
[求助]
求助高版本的1stOpt代跑个程序 已有1人参与
|
|
代码如下： Parameter p1,p2,p3,p4; Variable t,x1,x2,x3,x4; ODEFunction x2'=(x3/270)*(x1-x2)-p1*exp((-p2+p3*96485*(x4-p4))/(8.314*363))*x2^2; data; t,x1,x2,x3,x4; 0.1 165 69 100.00 361.00 0.3 115 71 110.00 356.00 0.5 103 62 112.50 358.00 0.7 150 69 118.75 380.00 0.9 128 95 112.50 380.00 1.1 140 63 100.00 383.00 1.3 110 36 100.00 383.00 1.5 103 33 118.75 370.00 1.7 96 31 137.50 367.00 1.9 70 25 100.00 356.00 2.1 83 31 125.00 373.00 2.3 90 81 137.50 363.00 2.5 100 88 125.00 353.00 2.7 112 60 137.50 346.00 2.9 119 35 137.50 353.00 3.1 66 30 112.50 353.00 3.3 70 35 125.00 373.00 3.5 160 63 112.50 370.00 3.7 147 74 112.50 370.00 3.9 140 86 112.50 373.00 4.1 108 76 107.50 373.00 4.3 109 79 75.00 363.00 4.5 76 50 125.00 353.00 4.7 79 47 131.00 363.00 4.9 82 43 112.00 362.00 5.1 136 61 135.00 373.00 5.3 145 69 125.00 363.00 5.5 143 52 110.00 366.00 5.7 159 63 125.00 370.00 5.9 178 76 100.00 370.00 6.1 112 56 100.00 360.00 6.3 121 70 81.25 364.00 6.5 152 47 112.50 363.00 6.7 157 60 100.00 373.00 6.9 162 82 93.75 383.00 7.1 170 63 118.75 363.00 7.3 62 49 112.50 291.00 7.5 83 78 112.50 356.00 7.7 227 95 137.50 359.00 7.9 226 75 125.00 360.00 8.1 223 47 125.00 355.00 8.3 116 44 143.75 366.00 8.5 130 70 143.75 370.00 |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有59人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
楼主的公式明显有问题： p1*exp((-p2+p3*96485*(x4-p4))/(8.314*363)) 完全等同于： exp((-p2+p3*96485*(x4))/(8.314*363)) p1和p4多余，去掉后效果是一样的，否则过拟合，参数组值不唯一。 计算用时(时:分:秒:微秒): 00:01:01:312 均方差(RMSE): 17.9229682919401 残差平方和(SSE): 13491.7772805435 相关系数(R): 0.290562575469603 相关系数之平方(R^2): 0.0844266102635286 修正R平方(Adj. R^2): 0.0374741287385814 确定系数(DC): 0.0816922624187681 F统计(F-Statistic): 2.44181044608464 参数 最佳估算 -------------------- ------------- p2 20003.4712902898 p3 0.000151432517221469 |
» 本帖已获得的红花（最新10朵）
yejianling2楼2016-11-07 21:01:17
yejianling
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1959.8
- 散金: 130
- 帖子: 337
- 在线: 279.6小时
- 虫号: 3156634
- 注册: 2014-04-23
- 性别: GG
- 专业: 自动化
3楼2016-11-08 11:35:41
30
