24小时热门版块排行榜    

查看: 586  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
回复此楼
已阅   回复此楼   关注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 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工学调剂 +7 用户573181 2026-03-20 12/600 2026-03-24 11:17 by 用户573181
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 07化学280分求调剂 +4 722865 2026-03-23 4/200 2026-03-24 00:01 by chixmc
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 材料专硕英一数二306 +8 z1z2z3879 2026-03-18 8/400 2026-03-23 20:49 by baobaoye
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
信息提示
请填处理意见