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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stOPt函数字符串有什么问题吗？帮我看看呗下面的代码指出以下问题或者帮我跑跑
111/2返回列表12下一页
查看: 1128  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 材料求调剂 +4 @taotao 2026-03-21 4/200 2026-03-21 18:52 by 学员8dgXkO
[考研] 280求调剂 +10 咕噜晓晓 2026-03-18 11/550 2026-03-21 17:40 by 学员8dgXkO
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 286求调剂 +8 Faune 2026-03-21 8/400 2026-03-21 17:16 by 学员8dgXkO
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿华南师大 070300（化学）304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
关闭
请填处理意见
关闭