应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stOPt函数字符串有什么问题吗？帮我看看呗下面的代码指出以下问题或者帮我跑跑
111/2返回列表12下一页
查看: 1144  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zengyin

新虫 (正式写手)

[求助] 1stOPt函数字符串有什么问题吗？帮我看看呗下面的代码指出以下问题或者帮我跑跑 已有2人参与

Title "PSO-3";
Parameters b2;
Variable x,y;
constant n=1000
Function y=7.925+0.44736771*0.5*Gamma(1.1789)/Gamma(1.5)+0.1118419*(exp)^(300*b2)*((x-300)^(b1))*sum(i=1:n)(((b2*(x-300))^(i))/((Gamma(1+i+b1)));
data;
//x        y
335.00         8.28
336.00         8.33
337.00         8.38
338.00         8.42
339.00         8.43
340.00         8.44
343.00         8.61
344.01         8.66
345.00         8.72
346.00         8.76
347.01         8.81
348.01         8.88
349.00         8.93
350.00         8.99
351.00         9.04
354.00         9.18
355.00         9.26
356.00         9.37
358.00         10.08
359.00         10.14
回复此楼

» 猜你喜欢
1楼 2017-03-10 19:54:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★
zengyin: 金币+2, ★有帮助 2017-03-13 16:39:29
引用回帖:
9楼: Originally posted by zengyin at 2017-03-11 19:17:27
Title "PSO-3";
Parameters b1,b2;
Variable x,y;
Function y=7.752+0.56*(x^0.45)*0.92/0.89-0.14*(2.7^(300*b1))*((x-300)^b2)*(1/(b1*(x-300)))*(1/Gamma(b2));
//x        y
331        8.145
335        8. ...

补充关键词 Data后，5.5试用版即可运行：
Root of Mean Square Error (RMSE): 0.441777266294964
Residual Sum of Squares (SSE): 1.36617007110536
Correlation Coef. (R): 0.907778259223795
R-Square: 0.824061367919383
Determination Coef. (DC): 0.57895112365208
Chi-Square: 0.162743048479003
F-Statistic: 11.8751059101946
Parameters Best Estimate
---------- -------------b1 0.00333701086674267
b2 0.440344141842967
====== Output Results =====
No         Observed        y Calculated y
1            8.145     7.38236085447038
2           8.275     7.96676155539426
3           8.424     8.34531328010477
4           8.661     8.98837858605962
5            8.81      9.26551911167519
6           9.366     9.96920057413036
7           10.24     10.1705468946164
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
10楼2017-03-11 19:57:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zengyin

新虫 (正式写手)
对了 除了参数那里我知道还有个b1，我上传的时候忘了。
Title "PSO-3";
Parameters b1,b2;
Variable x,y;
Constant n=1000
Function y=7.925+0.44736771*(x^(0.5))*Gamma(1.1789)/Gamma(1.5)+0.1118419*(exp)^(300*b2)*((x-300)^(b1))*sum(i=1:n)(((b2*(x-300))^(i))/((Gamma(1+i+b1)));
data;
//x        y
335.00         8.28
336.00         8.33
337.00         8.38
338.00         8.42
339.00         8.43
340.00         8.44
343.00         8.61
344.01         8.66
345.00         8.72
346.00         8.76
347.01         8.81
348.01         8.88
349.00         8.93
350.00         8.99
351.00         9.04
354.00         9.18
355.00         9.26
356.00         9.37
358.00         10.08
赞一下 回复此楼
2楼2017-03-10 19:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
Constant n=1000;

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
一团和气
3楼2017-03-10 23:25:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zengyin

新虫 (正式写手)
引用回帖:
3楼: Originally posted by lipenggg at 2017-03-10 23:25:27
Constant n=1000;

这个我也注意到了，我也已经改了，但是还是不运行啊
赞一下 回复此楼
4楼2017-03-11 10:14:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
1. (exp)^(300*b2)，exp表示以自然常数e为底的指数函数，正确写法是exp(300*b2)
2. (((b2*(x-300))^(i))/((Gamma(1+i+b1)))，最左边多了1个左括号；
CODE:
Parameters b1,b2;
Variable x,y;
Constant n=100;
Function y=7.925+0.44736771*(x^(0.5))*Gamma(1.1789)/Gamma(1.5)+0.1118419*exp(300*b2)*((x-300)^(b1))*Sum(i=1:n)(  ( b2*(x-300))^i/Gamma(1+i+b1)   );
data;
//x        y
335.00         8.28
336.00         8.33
337.00         8.38
338.00         8.42
339.00         8.43
340.00         8.44
343.00         8.61
344.01         8.66
345.00         8.72
346.00         8.76
347.01         8.81
348.01         8.88
349.00         8.93
350.00         8.99
351.00         9.04
354.00         9.18
355.00         9.26
356.00         9.37
358.00         10.08
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2017-03-11 11:43:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

1.5版本这种拟合不运行，可能版本太低。功能不强。可以使用优化模式进行拟合。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
一团和气
6楼2017-03-11 13:19:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

5.0试用版可以运行

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
一团和气
7楼2017-03-11 13:26:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

n值不能如太大值，我取到170，超过就不运行了。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
一团和气
8楼2017-03-11 13:33:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zengyin

新虫 (正式写手)
引用回帖:
5楼: Originally posted by 月只蓝 at 2017-03-11 11:43:27
1. (exp)^(300*b2)，exp表示以自然常数e为底的指数函数，正确写法是exp(300*b2)
2. (((b2*(x-300))^(i))/((Gamma(1+i+b1)))，最左边多了1个左括号；


Parameters b1,b2;
Variable x,y;
Constant n=100;
Fun ...

Title "PSO-3";
Parameters b1,b2;
Variable x,y;
Function y=7.752+0.56*(x^0.45)*0.92/0.89-0.14*(2.7^(300*b1))*((x-300)^b2)*(1/(b1*(x-300)))*(1/Gamma(b2));
//x        y
331        8.145
335        8.275
338        8.424
344        8.661
347        8.81
356        9.366
359        10.24
赞一下 回复此楼
9楼2017-03-11 19:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zengyin 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 通信工程求调剂！！！ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 一志愿A区211，22408 321求调剂 +5 随心所欲☆ 2026-04-15 6/300 2026-04-15 19:58 by PengTangGeoE
[考研] 290调剂生物0860 +37 哇哈哈，。 2026-04-11 43/2150 2026-04-15 18:55 by zhuwenxu
[考研] 297，工科调剂?河南农业大学本科 +11 河南农业大学-能 2026-04-14 11/550 2026-04-15 15:03 by 学员JpLReM
[考研] 22专硕求调剂 +9 haoyun上岸 2026-04-11 11/550 2026-04-15 14:47 by elainzaizai
[考研] 材料相关专业344求调剂双非工科学校或课题组 +21 hualkop 2026-04-12 22/1100 2026-04-14 17:44 by lhj2009
[考研] 恳请有学校收留 +3 柯淮然 2026-04-12 3/150 2026-04-14 16:25 by 逆水乘风
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 一志愿华中农微生物，288分，三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 085410 273求调剂 +10 X1999 2026-04-09 10/500 2026-04-12 09:24 by 逆水乘风
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 本人女孩 +7 吼吼， 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
信息提示
关闭
请填处理意见
关闭