24小时热门版块排行榜    

查看: 1599  |  回复: 11

瑞aaa

新虫 (正式写手)

[求助] 请高版本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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瑞aaa

新虫 (正式写手)

不好意思着急写错了,没有cf
Variable CA,CB,CC,CD,CE;
2楼2019-04-07 11:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瑞aaa

新虫 (正式写手)

所有参数都大于0,再补充一下

发自小木虫Android客户端
3楼2019-04-07 15:12:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
哪个是因变量?CE?

[ 发自手机版 http://muchong.com/3g ]
4楼2019-04-07 17:49:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瑞aaa

新虫 (正式写手)

就是一个式子,ca,cb,cc,cd,ce是变量,也就是数据,ka,kb,kc,kd,ke,kf是参数,需要求这几个参数

发自小木虫Android客户端
5楼2019-04-07 18:10:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瑞aaa

新虫 (正式写手)

引用回帖:
6楼: Originally posted by 独孤神宇 at 2019-04-07 18:39:55
计算结束原因: 达到收敛判断标准
优化算法: 通用全局优化算法(UGO1)
目标函数值(最小): 8.94776269027934E-5
均方差(RMSE): 0.00472962980459415
残差平方和(SSR): 8.94775923540212E-5
相关系数(R): 0.986304 ...

宇哥,我这个代码怎么写才能使目标与r对比呢?
7楼2019-04-07 19:24:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 瑞aaa at 2019-04-07 19:24:13
宇哥,我这个代码怎么写才能使目标与r对比呢?...

你已经将 r 定义为常数了,除非r也有一组数据

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

数值计算
8楼2019-04-07 20:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

瑞aaa

新虫 (正式写手)

送红花一朵
引用回帖:
8楼: Originally posted by 独孤神宇 at 2019-04-07 20:07:44
你已经将 r 定义为常数了,除非r也有一组数据
...

您看这样行吗,是不是这个代码就可以使我的目标函数变成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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
瑞aaa: 金币+10, ★★★★★最佳答案 2019-04-13 07:07:08
引用回帖:
9楼: Originally posted by 瑞aaa at 2019-04-07 20:53:13
您看这样行吗,是不是这个代码就可以使我的目标函数变成r了吗?
Parameters KA,KB,KC,KD,KE,KF;
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)^ ...

这样是可以计算,但是每一行数据都会计算出一个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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 瑞aaa 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode,4个jtjc了 +7 ziyangfang 2026-08-19 8/400 2026-08-19 11:40 by h初学者f
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +8 启萌科技 2026-08-12 30/1500 2026-08-19 10:22 by 启萌科技
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +52 医学老男孩 2026-08-13 114/5700 2026-08-19 09:38 by 医学老男孩
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见