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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求大神用软件工具拟合或者高版本的1stopt跑一下程序，获得参数a，b，c
51/1返回列表
查看: 661  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

LYP009

新虫 (初入文坛)

[求助] 求大神用软件工具拟合或者高版本的1stopt跑一下程序，获得参数a，b，c 已有2人参与

求大神用软件工具拟合的方法，最终根据数据拟合方程获得参数a，b，c；或者用高版本的1stopt跑一下程序获得参数a，b，c。
方程：y=a*b*x/0.09*[1-exp*(-0.09/b/x)]/(1+0.3/x*c)

以下为拟合数据及1stopt拟合程序：
Title "Type your title here";
Parameters a,b,c;
Variable x,y;
Function y=a*b*x/0.09*[1-exp(-(0.09/b/x))]/(1+0.3/x*c);
Data;
-1.51E-06        4.48E-05
1.00E+00        1.58E-03
2.00E+00        2.55E-03
3.00E+00        3.25E-03
4.00E+00        3.77E-03
5.00E+00        4.17E-03
6.00E+00        4.51E-03
7.00E+00        4.81E-03
8.00E+00        5.08E-03
9.00E+00        5.33E-03
1.00E+01        5.57E-03
1.10E+01        5.79E-03
1.20E+01        6.00E-03
1.30E+01        6.20E-03
1.40E+01        6.39E-03
1.50E+01        6.57E-03
1.60E+01        6.74E-03
1.70E+01        6.90E-03
1.80E+01        7.05E-03
1.90E+01        7.19E-03
2.00E+01        7.32E-03
2.10E+01        7.45E-03
2.20E+01        7.57E-03
2.30E+01        7.68E-03
2.40E+01        7.80E-03
2.50E+01        7.90E-03
2.60E+01        8.00E-03
2.70E+01        8.10E-03
2.80E+01        8.20E-03
2.90E+01        8.29E-03
3.00E+01        8.38E-03
3.10E+01        8.47E-03
3.20E+01        8.55E-03
3.30E+01        8.63E-03
3.40E+01        8.71E-03
3.50E+01        8.78E-03
3.60E+01        8.86E-03
3.70E+01        8.94E-03
3.80E+01        9.00E-03
3.90E+01        9.07E-03
4.00E+01        9.14E-03
4.10E+01        9.21E-03
4.20E+01        9.28E-03
4.30E+01        9.34E-03
4.40E+01        9.41E-03
4.50E+01        9.47E-03
4.60E+01        9.53E-03
4.70E+01        9.59E-03
4.80E+01        9.66E-03
4.90E+01        9.72E-03
5.00E+01        9.78E-03
5.10E+01        9.84E-03
5.20E+01        9.90E-03
5.30E+01        9.96E-03
5.40E+01        1.00E-02
5.50E+01        1.01E-02
5.60E+01        1.01E-02
5.70E+01        1.02E-02
5.80E+01        1.02E-02
5.90E+01        1.03E-02
6.00E+01        1.03E-02
6.10E+01        1.04E-02
6.20E+01        1.05E-02
6.27E+01        1.05E-02
6.40E+01        1.05E-02
6.50E+01        1.06E-02
6.60E+01        1.06E-02
6.70E+01        1.07E-02
6.80E+01        1.07E-02
6.90E+01        1.08E-02
7.00E+01        1.08E-02
7.10E+01        1.09E-02
7.20E+01        1.09E-02
7.30E+01        1.10E-02
7.40E+01        1.10E-02
7.50E+01        1.10E-02
7.60E+01        1.11E-02
7.70E+01        1.11E-02
7.80E+01        1.12E-02
7.90E+01        1.12E-02
8.00E+01        1.12E-02
8.10E+01        1.13E-02
8.20E+01        1.13E-02
8.30E+01        1.14E-02
8.40E+01        1.14E-02
8.50E+01        1.15E-02
8.60E+01        1.15E-02
8.70E+01        1.15E-02
8.80E+01        1.16E-02
8.90E+01        1.16E-02
9.00E+01        1.16E-02
9.10E+01        1.17E-02
9.20E+01        1.17E-02
9.30E+01        1.18E-02
9.40E+01        1.18E-02
9.50E+01        1.18E-02
9.60E+01        1.19E-02
9.70E+01        1.19E-02
9.80E+01        1.20E-02
9.90E+01        1.20E-02
1.00E+02        1.20E-02
1.01E+02        1.21E-02
1.02E+02        1.21E-02
1.03E+02        1.22E-02
1.04E+02        1.22E-02
1.05E+02        1.22E-02
1.06E+02        1.23E-02
1.07E+02        1.23E-02
1.08E+02        1.24E-02
1.09E+02        1.24E-02
1.10E+02        1.24E-02
1.11E+02        1.25E-02
1.12E+02        1.25E-02
1.13E+02        1.26E-02
1.14E+02        1.26E-02
1.15E+02        1.26E-02
1.16E+02        1.27E-02
1.17E+02        1.27E-02
1.18E+02        1.28E-02
1.19E+02        1.28E-02
1.20E+02        1.28E-02
1.21E+02        1.29E-02
1.22E+02        1.29E-02
1.23E+02        1.29E-02
1.24E+02        1.30E-02
1.25E+02        1.30E-02
1.26E+02        1.30E-02
1.27E+02        1.31E-02
1.28E+02        1.31E-02
1.29E+02        1.31E-02
1.30E+02        1.32E-02
1.31E+02        1.32E-02
1.32E+02        1.32E-02
1.33E+02        1.33E-02
1.34E+02        1.33E-02
1.35E+02        1.33E-02
1.36E+02        1.34E-02
1.37E+02        1.34E-02
1.38E+02        1.34E-02
1.39E+02        1.34E-02
1.40E+02        1.35E-02
1.41E+02        1.35E-02
1.42E+02        1.35E-02
1.43E+02        1.35E-02
1.44E+02        1.36E-02
1.45E+02        1.36E-02
1.46E+02        1.36E-02
1.47E+02        1.37E-02
1.48E+02        1.37E-02
1.49E+02        1.37E-02
1.50E+02        1.38E-02
1.51E+02        1.38E-02
1.52E+02        1.38E-02
1.53E+02        1.39E-02
1.54E+02        1.39E-02
1.55E+02        1.39E-02
1.56E+02        1.39E-02
1.57E+02        1.40E-02
1.58E+02        1.40E-02
1.59E+02        1.40E-02
1.60E+02        1.41E-02
1.61E+02        1.41E-02
1.62E+02        1.41E-02
1.63E+02        1.42E-02
1.64E+02        1.42E-02
1.65E+02        1.42E-02
1.66E+02        1.42E-02
1.67E+02        1.43E-02
1.68E+02        1.43E-02
1.69E+02        1.43E-02
1.70E+02        1.44E-02
1.71E+02        1.44E-02
1.72E+02        1.44E-02
1.73E+02        1.44E-02
1.74E+02        1.45E-02
1.75E+02        1.45E-02
1.76E+02        1.45E-02
1.77E+02        1.46E-02
1.78E+02        1.46E-02
1.79E+02        1.46E-02
1.80E+02        1.46E-02
1.81E+02        1.47E-02
1.82E+02        1.47E-02
1.83E+02        1.47E-02
1.84E+02        1.48E-02
1.85E+02        1.48E-02
1.86E+02        1.48E-02
1.87E+02        1.48E-02
1.88E+02        1.49E-02
1.89E+02        1.49E-02
1.90E+02        1.49E-02
1.91E+02        1.50E-02
1.92E+02        1.50E-02
1.93E+02        1.50E-02
1.94E+02        1.50E-02
1.95E+02        1.51E-02
1.96E+02        1.51E-02
1.97E+02        1.51E-02
1.98E+02        1.51E-02
1.99E+02        1.52E-02
2.00E+02        1.52E-02
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-02-23 13:31:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shikang999

新虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
a=0.01669456950455   
b=-0.02877239797410   
c=116.932662072696
赞一下 回复此楼
我还是那个我，过去如此，未来亦如此！
3楼2016-02-23 16:21:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

LYP009

新虫 (初入文坛)
第一行数据有异常，请删除

发自小木虫Android客户端
赞一下 回复此楼
2楼2016-02-23 13:54:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
2楼: Originally posted by LYP009 at 2016-02-23 13:54:05
第一行数据有异常，请删除

不是1stopt版本问题，是代码没写对，已修改如下：
CODE:
Parameters a,b,c;
Variable x,y;
Function y=a*b*x/0.09*(1-exp(-(0.09/b/x)))/(1+0.3/x*c);
Data;

1.00E+00        1.58E-03
2.00E+00        2.55E-03
3.00E+00        3.25E-03
4.00E+00        3.77E-03
5.00E+00        4.17E-03
6.00E+00        4.51E-03
7.00E+00        4.81E-03
8.00E+00        5.08E-03
9.00E+00        5.33E-03
1.00E+01        5.57E-03
1.10E+01        5.79E-03
1.20E+01        6.00E-03
1.30E+01        6.20E-03
1.40E+01        6.39E-03
1.50E+01        6.57E-03
1.60E+01        6.74E-03
1.70E+01        6.90E-03
1.80E+01        7.05E-03
1.90E+01        7.19E-03
2.00E+01        7.32E-03
2.10E+01        7.45E-03
2.20E+01        7.57E-03
2.30E+01        7.68E-03
2.40E+01        7.80E-03
2.50E+01        7.90E-03
2.60E+01        8.00E-03
2.70E+01        8.10E-03
2.80E+01        8.20E-03
2.90E+01        8.29E-03
3.00E+01        8.38E-03
3.10E+01        8.47E-03
3.20E+01        8.55E-03
3.30E+01        8.63E-03
3.40E+01        8.71E-03
3.50E+01        8.78E-03
3.60E+01        8.86E-03
3.70E+01        8.94E-03
3.80E+01        9.00E-03
3.90E+01        9.07E-03
4.00E+01        9.14E-03
4.10E+01        9.21E-03
4.20E+01        9.28E-03
4.30E+01        9.34E-03
4.40E+01        9.41E-03
4.50E+01        9.47E-03
4.60E+01        9.53E-03
4.70E+01        9.59E-03
4.80E+01        9.66E-03
4.90E+01        9.72E-03
5.00E+01        9.78E-03
5.10E+01        9.84E-03
5.20E+01        9.90E-03
5.30E+01        9.96E-03
5.40E+01        1.00E-02
5.50E+01        1.01E-02
5.60E+01        1.01E-02
5.70E+01        1.02E-02
5.80E+01        1.02E-02
5.90E+01        1.03E-02
6.00E+01        1.03E-02
6.10E+01        1.04E-02
6.20E+01        1.05E-02
6.27E+01        1.05E-02
6.40E+01        1.05E-02
6.50E+01        1.06E-02
6.60E+01        1.06E-02
6.70E+01        1.07E-02
6.80E+01        1.07E-02
6.90E+01        1.08E-02
7.00E+01        1.08E-02
7.10E+01        1.09E-02
7.20E+01        1.09E-02
7.30E+01        1.10E-02
7.40E+01        1.10E-02
7.50E+01        1.10E-02
7.60E+01        1.11E-02
7.70E+01        1.11E-02
7.80E+01        1.12E-02
7.90E+01        1.12E-02
8.00E+01        1.12E-02
8.10E+01        1.13E-02
8.20E+01        1.13E-02
8.30E+01        1.14E-02
8.40E+01        1.14E-02
8.50E+01        1.15E-02
8.60E+01        1.15E-02
8.70E+01        1.15E-02
8.80E+01        1.16E-02
8.90E+01        1.16E-02
9.00E+01        1.16E-02
9.10E+01        1.17E-02
9.20E+01        1.17E-02
9.30E+01        1.18E-02
9.40E+01        1.18E-02
9.50E+01        1.18E-02
9.60E+01        1.19E-02
9.70E+01        1.19E-02
9.80E+01        1.20E-02
9.90E+01        1.20E-02
1.00E+02        1.20E-02
1.01E+02        1.21E-02
1.02E+02        1.21E-02
1.03E+02        1.22E-02
1.04E+02        1.22E-02
1.05E+02        1.22E-02
1.06E+02        1.23E-02
1.07E+02        1.23E-02
1.08E+02        1.24E-02
1.09E+02        1.24E-02
1.10E+02        1.24E-02
1.11E+02        1.25E-02
1.12E+02        1.25E-02
1.13E+02        1.26E-02
1.14E+02        1.26E-02
1.15E+02        1.26E-02
1.16E+02        1.27E-02
1.17E+02        1.27E-02
1.18E+02        1.28E-02
1.19E+02        1.28E-02
1.20E+02        1.28E-02
1.21E+02        1.29E-02
1.22E+02        1.29E-02
1.23E+02        1.29E-02
1.24E+02        1.30E-02
1.25E+02        1.30E-02
1.26E+02        1.30E-02
1.27E+02        1.31E-02
1.28E+02        1.31E-02
1.29E+02        1.31E-02
1.30E+02        1.32E-02
1.31E+02        1.32E-02
1.32E+02        1.32E-02
1.33E+02        1.33E-02
1.34E+02        1.33E-02
1.35E+02        1.33E-02
1.36E+02        1.34E-02
1.37E+02        1.34E-02
1.38E+02        1.34E-02
1.39E+02        1.34E-02
1.40E+02        1.35E-02
1.41E+02        1.35E-02
1.42E+02        1.35E-02
1.43E+02        1.35E-02
1.44E+02        1.36E-02
1.45E+02        1.36E-02
1.46E+02        1.36E-02
1.47E+02        1.37E-02
1.48E+02        1.37E-02
1.49E+02        1.37E-02
1.50E+02        1.38E-02
1.51E+02        1.38E-02
1.52E+02        1.38E-02
1.53E+02        1.39E-02
1.54E+02        1.39E-02
1.55E+02        1.39E-02
1.56E+02        1.39E-02
1.57E+02        1.40E-02
1.58E+02        1.40E-02
1.59E+02        1.40E-02
1.60E+02        1.41E-02
1.61E+02        1.41E-02
1.62E+02        1.41E-02
1.63E+02        1.42E-02
1.64E+02        1.42E-02
1.65E+02        1.42E-02
1.66E+02        1.42E-02
1.67E+02        1.43E-02
1.68E+02        1.43E-02
1.69E+02        1.43E-02
1.70E+02        1.44E-02
1.71E+02        1.44E-02
1.72E+02        1.44E-02
1.73E+02        1.44E-02
1.74E+02        1.45E-02
1.75E+02        1.45E-02
1.76E+02        1.45E-02
1.77E+02        1.46E-02
1.78E+02        1.46E-02
1.79E+02        1.46E-02
1.80E+02        1.46E-02
1.81E+02        1.47E-02
1.82E+02        1.47E-02
1.83E+02        1.47E-02
1.84E+02        1.48E-02
1.85E+02        1.48E-02
1.86E+02        1.48E-02
1.87E+02        1.48E-02
1.88E+02        1.49E-02
1.89E+02        1.49E-02
1.90E+02        1.49E-02
1.91E+02        1.50E-02
1.92E+02        1.50E-02
1.93E+02        1.50E-02
1.94E+02        1.50E-02
1.95E+02        1.51E-02
1.96E+02        1.51E-02
1.97E+02        1.51E-02
1.98E+02        1.51E-02
1.99E+02        1.52E-02
2.00E+02        1.52E-02

计算结果：
相关系数(R): 0.983468907853328
相关系数之平方(R^2): 0.967211092714217
决定系数(DC): 0.963748683704118
卡方系数(Chi-Square): 0.00443900431290186
F统计(F-Statistic): 5840.62758445293

参数        最佳估算
----------        -------------
a                  0.0166945694989546
b                  -0.0287723976673137
c                  116.932662127234
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2016-02-23 16:44:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上意见出来了 +6 黄鸟于飞Chao 2026-08-29 12/600 2026-08-29 17:20 by 黄鸟于飞Chao
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 10/500 2026-08-29 16:49 by zICmwzsBXjbN
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:02 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
关闭
请填处理意见
关闭