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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stOPt函数字符串有什么问题吗？帮我看看呗下面的代码指出以下问题或者帮我跑跑
111/2返回列表12下一页
查看: 1145  |  回复: 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的回帖

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: 金币+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 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试324 中药学 一志愿天中医 求调剂 +3 李的Lucky 2026-04-10 3/150 2026-04-15 23:54 by 陈皮皮
[考研] 药学求调剂 +12 喽哈加油 2026-04-14 14/700 2026-04-15 23:49 by 陈皮皮
[考研] 材料相关专业344求调剂双非工科学校或课题组 +21 hualkop 2026-04-12 23/1150 2026-04-15 22:02 by hualkop
[考研] 297，工科调剂?河南农业大学本科 +11 河南农业大学-能 2026-04-14 11/550 2026-04-15 15:03 by 学员JpLReM
[考研] 310求调剂 +16 666真好 2026-04-11 18/900 2026-04-15 13:28 by 黑科技矿业
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 367求调剂 +11 hffQAQ 2026-04-09 11/550 2026-04-14 17:48 by lhj2009
[考研] 化学070300 求调剂 +23 哈哈哈^_^ 2026-04-12 23/1150 2026-04-14 16:30 by zhouxiaoyu
[考研] 求调剂 +16 雪逢冬 2026-04-10 16/800 2026-04-14 14:27 by 逆水乘风
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 初试分332，一志愿报考西北工业大学， +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 337求调剂 +4 Gky09300550， 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
信息提示
关闭
请填处理意见
关闭