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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stOPt函数字符串有什么问题吗？帮我看看呗下面的代码指出以下问题或者帮我跑跑
111/2返回列表12下一页
查看: 1287  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了，能看出什么问题？ +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息，中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 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
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护，明天很可能放榜！ +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
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
关闭
请填处理意见
关闭