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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 1stopt软件非线性拟合求助
111/2返回列表12下一页
查看: 1518  |  回复: 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个字符以上)
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭