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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stOPt函数字符串有什么问题吗？帮我看看呗下面的代码指出以下问题或者帮我跑跑
111/2返回列表12下一页
查看: 1033  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭