| 查看: 1314 | 回复: 11 | ||
[求助]
请高版本1stopt代跑 已有2人参与
|
|
这是我自己写的代码,麻烦大佬代跑一下,如果写的有错还请指教!拜托 Parameters KA,KB,KC,KD,KE,KF; Variable CA,CB,CC,CD,CE,CF; Constant r=3.2707; Function r=KF(KA*KB*CA*CB-KC*KD*CC*CD/1.34)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE)^2; Data; 3.914545455 2.146060606 4.249745353 2.737478991 0.588533232 3.828484848 2.232121212 4.308152457 2.706725388 0.56087973 3.767272727 2.293333333 4.340183215 2.703875994 0.531698367 3.711515152 2.349090909 4.37840252 2.683194959 0.514160096 |
» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
物理学I论文润色/翻译怎么收费?
已经有64人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
反铁磁体中的磁性切换:两种不同的机制已成功可视化
已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
新西兰Robinson研究所招收全奖PhD
已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
2楼2019-04-07 11:10:13
3楼2019-04-07 15:12:47
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2019-04-07 17:49:31
5楼2019-04-07 18:10:54
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与,应助指数 +1
|
计算结束原因: 达到收敛判断标准 优化算法: 通用全局优化算法(UGO1) 目标函数值(最小): 8.94776269027934E-5 均方差(RMSE): 0.00472962980459415 残差平方和(SSR): 8.94775923540212E-5 相关系数(R): 0.986304513331934 确定系数(DC): 0.972186335183512 参数 最佳估算 ---------- ------------- ka 11.8201249107019 kb 23.5726553856025 kc 0.843649223024818 kd 0.229359471088254 ke 6.45725464510755 kf 15.6674603816671 ====== 结果输出 ====== 目标 ce 计算 ce 0.588533232 0.585434820888929 0.56087973 0.563011233652204 0.531698367 0.538199947346761 0.514160096 0.508410001555901 ====== 计算结束 ====== |

6楼2019-04-07 18:39:55
7楼2019-04-07 19:24:13
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟

8楼2019-04-07 20:07:44
送红花一朵 |
您看这样行吗,是不是这个代码就可以使我的目标函数变成r了吗? Parameters KA[0,],KB[0,],KC[0,],KD[0,],KE[0,],KF[0,]; Variable CA,CB,CC,CD,CE; Function r=KF(KA*KB*CA*CB-KC*KD*CC*CD/1.34)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE)^2; Data; 3.2707 3.914545455 2.146060606 4.249745353 2.737478991 0.588533232 3.2707 3.828484848 2.232121212 4.308152457 2.706725388 0.56087973 3.2707 3.767272727 2.293333333 4.340183215 2.703875994 0.531698367 3.2707 3.711515152 2.349090909 4.37840252 2.683194959 0.514160096 我再问您几个小问题,怎么选取我的目标函数是哪个呢?这两种不同目标函数的选取对我最后得到的参数有没有影响呢? 拜托了,搞懂这些就去花钱让人代跑,最后一步啦!!!! |
9楼2019-04-07 20:53:13
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
瑞aaa: 金币+10, ★★★★★最佳答案 2019-04-13 07:07:08
瑞aaa: 金币+10, ★★★★★最佳答案 2019-04-13 07:07:08
|
这样是可以计算,但是每一行数据都会计算出一个r值出来。。。 ****************** Parameters KA,KB,KC,KD,KE,KF; ParameterDomain = [0, ]; Variable CA,CB,CC,CD,CE,r; Function r=KF*(KA*KB*CA*CB-KC*KD*CC*CD/1.34)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE)^2; Data; //CA,CB,CC,CD,CE,r; 3.914545455 2.146060606 4.249745353 2.737478991 0.588533232 3.2707 3.828484848 2.232121212 4.308152457 2.706725388 0.56087973 3.2707 3.767272727 2.293333333 4.340183215 2.703875994 0.531698367 3.2707 3.711515152 2.349090909 4.37840252 2.683194959 0.514160096 3.2707 ************************************************* 均方差(RMSE): 0.00180815007608335 残差平方和(SSR): 1.30776267905609E-5 相关系数(R): 0 相关系数之平方(R^2): 0 修正R平方(Adj. R^2): 2 确定系数(DC): 0 卡方系数(Chi-Square): 1.99922523760395E-6 F统计(F-Statistic): -0.4 参数 最佳估算 ---------- ------------- ka 7.26870923240464 kb 76.8033683805592 kc 1.81885884331089E-8 kd 34.1212240739148 ke 80.0470405155196 kf 78.1355455967981 ====== 结果输出 ====== No. 实测值 r 计算值 r 1 3.2707 3.2692790256648 2 3.2707 3.27221542633868 3 3.2707 3.27274334006356 4 3.2707 3.26855834103416 |

10楼2019-04-07 21:33:19













回复此楼