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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stopt软件非线性拟合求助
111/2返回列表12下一页
查看: 2022  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者天蓝的将赠送您 90 个金币

天蓝的

银虫 (正式写手)

[求助] 1stopt软件非线性拟合求助 已有2人参与

Parameter C[,100], E[5000,], F=88[5,];
Variable x, y;
Function y=C*(1-(165.289*(x/(F*exp(E/6428.74225)))^0.5)*(185.185*(x/(F*exp(E/6428.74225)))^0.5)*(140.252*(x/(F*exp(E/6428.74225)))^0.5));
Data;
//x        y
0        0
10        1.86
22        2.85
34        3.49
46        3.99
58        4.36
71.5        4.75
83        4.96
107        5.35
131        5.79
167        6.32
191        6.63
200        6.73


拟合结果
Iterations: 35
Elapsed Time (Hr:Min:Sec:Msec): 00:00:14:105
Algorithms: Universal Global Optimation(UGO1)
Stop Reason: Convergence tolerance reached
Root of Mean Square Error (RMSE): 1.89309282451297
Residual Sum of Squares (SSE): 46.5894057488923
Correlation Coef. (R): 0
R-Square: 0
Determination Coef. (DC): 1.52511654600143E-16
Chi-Square: 10.6288164821448
F-Statistic: 5
Parameters Best Estimate
---------- -------------
c 4.38331076909242
e 200415.247116993
f 847455.285339582
====== Output Results =====
No Observed y Calculated y
1 0           4.38331076909242
2 1.86      4.38331076909242
3 2.85      4.38331076909242
4 3.49      4.38331076909242
5 3.99      4.38331076909242
6 4.36      4.38331076909242
7 4.75      4.38331076909242
8 4.96      4.38331076909242
9 5.35      4.38331076909242
10 5.79    4.38331076909242
11 6.32    4.38331076909242
12 6.63    4.38331076909242
13 6.73    4.38331076909242
目标根据下面公式拟合出抛物线，但是拟合曲线是一条直线（红色部分）；一开始我没有添加参数范围，可是结果还是一条直线，后来加了大概的参数范围，结果还是直线，希望1stopt软件非线性拟合的大神指导一下，谢谢

目标公式和拟合图 见附件图片
回复此楼

» 猜你喜欢
1楼 2022-03-28 16:02:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天蓝的

银虫 (正式写手)
内容已删除
回复此楼
2楼2022-03-28 16:05:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天蓝的

银虫 (正式写手)
目标公式和拟合图片在下面分享百度网盘里面

链接: https://pan.baidu.com/s/1BlEsbyaUPfhmEyxJb2bsgQ 提取码: my4u 复制这段内容后打开百度网盘手机App，操作更方便哦
赞一下 回复此楼
3楼2022-03-28 16:19:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
这说明 模型和数据不匹配，检查一下模型，可能模型不合适。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

天蓝的
数值计算
4楼2022-03-29 09:27:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天蓝的

银虫 (正式写手)
送红花一朵
引用回帖:
4楼: Originally posted by 独孤神宇 at 2022-03-29 09:27:21
这说明 模型和数据不匹配，检查一下模型，可能模型不合适。

好的，谢谢，模型应该没问题，好多相关文献类似抛物线规律的数据都是很好的拟合这个模型了，感觉可能是参数范围设置的不对，不知道要修改哪些默认参数，重新拟合
赞一下 回复此楼
5楼2022-03-29 22:14:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
天蓝的: 金币+10, ★★★很有帮助 2022-03-31 21:16:52
引用回帖:
5楼: Originally posted by 天蓝的 at 2022-03-29 22:14:50
好的，谢谢，模型应该没问题，好多相关文献类似抛物线规律的数据都是很好的拟合这个模型了，感觉可能是参数范围设置的不对，不知道要修改哪些默认参数，重新拟合...

不约束范围也是一样的结果。
检查一下模型是否输入有误

发自小木虫Android客户端
赞一下 回复此楼
数值计算
6楼2022-03-30 06:54:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
独孤神宇: 金币+10, 鼓励交流 2022-03-31 20:49:47
引用回帖:
5楼: Originally posted by 天蓝的 at 2022-03-29 22:14:50
好的，谢谢，模型应该没问题，好多相关文献类似抛物线规律的数据都是很好的拟合这个模型了，感觉可能是参数范围设置的不对，不知道要修改哪些默认参数，重新拟合...

楼主公式应该不对吧：
原公式：y=C*(1-(165.289*(x/(F*exp(E/6428.74225)))^0.5)*(185.185*(x/(F*exp(E/6428.74225)))^0.5)*(140.252*(x/(F*exp(E/6428.74225)))^0.5));
等价于：y=C*(1-(165.289*185.185*140.252*(x/(F*exp(E1)))^1.5));
“F*exp(E1)”可视为一个参数（否则无穷多组解），因此该公式又可以简化为：y=C-a*x^1.5;
效果能好吗？还不如一个简单的如：y=a*x^b;
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

天蓝的
7楼2022-03-31 19:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天蓝的

银虫 (正式写手)
引用回帖:
6楼: Originally posted by 独孤神宇 at 2022-03-30 06:54:47
不约束范围也是一样的结果。
检查一下模型是否输入有误
...

谢谢，仔细检查发现，公式输入错误，正确应该是y=C*(1-(1-165.289*Sqrt(x/F)*exp(-E/12857.4845))*(1-185.185*Sqrt(x/F)*exp(-E/12857.4845))*(1-140.252*Sqrt(x/F)*exp(-E/12857.4845)))。还有一个问题，现在拟合是抛物线了，但是单击1次拟合的结果，和单击2次到的拟合结果，为啥参数差别很大呀？相关系数多是0.999，这是说明有好几组最优解吗？热执行和执行2个菜单有啥区别呀？
赞一下 回复此楼
8楼2022-03-31 21:20:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天蓝的

银虫 (正式写手)
送红花一朵
引用回帖:
7楼: Originally posted by dingd at 2022-03-31 19:01:06
楼主公式应该不对吧：
原公式：y=C*(1-(165.289*(x/(F*exp(E/6428.74225)))^0.5)*(185.185*(x/(F*exp(E/6428.74225)))^0.5)*(140.252*(x/(F*exp(E/6428.74225)))^0.5));
等价于：y=C*(1-(165.289*185.185*140.25 ...

谢谢，我的模型参数确实输入不正确，正常的话，应该是y=C*(1-(1-165.289*Sqrt(x/F)*exp(-E/12857.4845))*(1-185.185*Sqrt(x/F)*exp(-E/12857.4845))*(1-140.252*Sqrt(x/F)*exp(-E/12857.4845)))，这样就不是你说的过拟合了。还有一个问题，现在拟合是抛物线了，但是单击1次拟合的结果，和单击2次到的拟合结果，为啥参数差别很大呀？相关系数多是0.999，这是说明有好几组最优解吗？热执行和执行2个菜单有啥区别呀？
赞一下 回复此楼
9楼2022-03-31 21:25:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
引用回帖:
9楼: Originally posted by 天蓝的 at 2022-03-31 21:25:12
谢谢，我的模型参数确实输入不正确，正常的话，应该是y=C*(1-(1-165.289*Sqrt(x/F)*exp(-E/12857.4845))*(1-185.185*Sqrt(x/F)*exp(-E/12857.4845))*(1-140.252*Sqrt(x/F)*exp(-E/12857.4845)))，这样就不是你说的 ...

模型公式中三处共同部分“Sqrt(x/F)*exp(-E/12857.4845)”完全等价于“Sqrt(x/F)*E1”，其中E1=exp(-E/12857.4845)，
“Sqrt(x/F)*E1”中参数F和E1又是过拟合，又完全多组合，可以简化为“H1*sqrt(x)”，一个参数即可，结果就会应该稳定唯一了，效果还一样
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

天蓝的
10楼2022-03-31 22:43:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 天蓝的 的主题更新
111/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +46 医学老男孩 2026-08-13 102/5100 2026-08-18 14:41 by qinwangshen
[基金申请] 重要消息，中午系统在维护 +5 yuleib84 2026-08-18 6/300 2026-08-18 14:34 by 阿荣喝酒
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 KZyU1THERRVo 2026-08-18 4/200 2026-08-18 13:28 by vHWVV9D5Deb5
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 时间戳变了，能看出什么问题？ +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 感觉是下周放榜了 +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
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +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
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
关闭
请填处理意见
关闭