| 查看: 499 | 回复: 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 |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有134人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
yejianling
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1959.8
- 散金: 130
- 帖子: 337
- 在线: 279.6小时
- 虫号: 3156634
- 注册: 2014-04-23
- 性别: GG
- 专业: 自动化
送红花一朵 
3楼2016-11-08 11:35:41
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朵）
2楼2016-11-07 21:01:17
