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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求大神用软件工具拟合或者高版本的1stopt跑一下程序，获得参数a，b，c
41/1返回列表
查看: 559  |  回复: 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的回帖

LYP009

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

发自小木虫Android客户端
赞一下 回复此楼
2楼2016-02-23 13:54:05
已阅   回复此楼   关注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的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +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的回帖
相关版块跳转 我要订阅楼主 LYP009 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭