24小时热门版块排行榜    

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

缘牵一线624

木虫 (小有名气)

[求助] 麻烦大神帮看一下 1stopt 拟合结果很不满意呀,帮忙跑一下 已有1人参与

第一个:
Title "Type your title here";
Parameters a,b;
Variable x,y;
Function y=((x*a*((2*pi)^0.5))^(-1))*exp(-(log(x)-b)^2/(2*a^2));
Data;
7.75        0.009555311
9.25        0.013965454
10.75        0.020580669
12.25        0.030871003
13.75        0.035281147
15.25        0.027930908
16.75        0.029400956
18.25        0.023520764
19.75        0.019110621
21.25        0.020580669
22.75        0.020580669
24.25        0.021315693
25.75        0.021315693
27.25        0.009555311
28.75        0.010290334
30.25        0.006615215
31.75        0.005880191
33.25        0.007350239
34.75        0.006615215
36.25        0.009555311
37.75        0.010290334
39.25        0.005880191
40.75        0.007350239
42.25        0.003675119
43.75        0.005880191
45.25        0.002205072
46.75        0.005880191
48.25        0.002205072
49.75        0.000735024
51.25        0.001470048
52.75        0.002940096
54.25        0.007350239
55.75        0.004410143
57.25        0.002940096
58.75        0.005880191
60.25        0.003675119
61.75        0.008820287
63.25        0.002205072
64.75        0.002205072
66.25        0.004410143
67.75        0.002205072
69.25        0.005145167
70.75        0.003675119
72.25        0.004410143
73.75        0.005145167
75.25        0.007350239
76.75        0.006615215
78.25        0.002205072
79.75        0.001470048
81.25        0.005145167
82.75        0.002940096
84.25        0.005145167
85.75        0.005880191
87.25        0.002205072
88.75        0.006615215
90.25        0.003675119
91.75        0.006615215
93.25        0.002205072
94.75        0.004410143
96.25        0.003675119
97.75        0.002940096
99.25        0.002205072
100.75        0.002940096
102.25        0.002205072
103.75        0.001470048
105.25        0.002940096
106.75        0.004410143
108.25        0.005880191
109.75        0.004410143
111.25        0.002205072
112.75        0.005145167
114.25        0.005145167
115.75        0.004410143
117.25        0.005145167
118.75        0.002940096
120.25        0.003675119
121.75        0.002940096
123.25        0.004410143
124.75        0.003675119
126.25        0.002940096
127.75        0.004410143
129.25        0.002940096
130.75        0.005145167
132.25        0.000735024
133.75        0.002205072
135.25        0.003675119
136.75        0.001470048
138.25        0
139.75        0.002940096
141.25        0.000735024
142.75        0.002940096
144.25        0.002205072
145.75        0.001470048
147.25        0.001470048
148.75        0.002940096
150.25        0.002940096
151.75        0.000735024
153.25        0.003675119
154.75        0.002205072
156.25        0.003675119
157.75        0.001470048
159.25        0.001470048
160.75        0.001470048
162.25        0.000735024
163.75        0
165.25        0.001470048
166.75        0.000735024
168.25        0
169.75        0
171.25        0
172.75        0
174.25        0.000735024
175.75        0.000735024
177.25        0
178.75        0.002205072
180.25        0.000735024
181.75        0.001470048
183.25        0
184.75        0.001470048
186.25        0
187.75        0
189.25        0.001470048
190.75        0
192.25        0
193.75        0.000735024
195.25        0.000735024
196.75        0
198.25        0
199.75        0
201.25        0.000735024
202.75        0.000735024
结果:优化算法: 改进差分进化算法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.00403210574738573
残差平方和(SSE): 0.00212978185531124
相关系数(R): 0.795116938932438
相关系数之平方(R^2): 0.63221094657729
决定系数(DC): 0.628563679575216
卡方系数(Chi-Square): 0.123335132084722
F统计(F-Statistic): 221.744533583866

参数        最佳估算
----------        -------------
a                  4.11513127199436
b                  -4.67953427480393

第二个:
Title "Type your title here";
Parameters a,b,c;
Variable x,y;
Function (((x-c)*a*((2*pi)^0.5))^(-1))*exp(-(log(x-c)-b)^2/(2*a^2));
Data;
7.75        0.009555311
9.25        0.013965454
10.75        0.020580669
12.25        0.030871003
13.75        0.035281147
15.25        0.027930908
16.75        0.029400956
18.25        0.023520764
19.75        0.019110621
21.25        0.020580669
22.75        0.020580669
24.25        0.021315693
25.75        0.021315693
27.25        0.009555311
28.75        0.010290334
30.25        0.006615215
31.75        0.005880191
33.25        0.007350239
34.75        0.006615215
36.25        0.009555311
37.75        0.010290334
39.25        0.005880191
40.75        0.007350239
42.25        0.003675119
43.75        0.005880191
45.25        0.002205072
46.75        0.005880191
48.25        0.002205072
49.75        0.000735024
51.25        0.001470048
52.75        0.002940096
54.25        0.007350239
55.75        0.004410143
57.25        0.002940096
58.75        0.005880191
60.25        0.003675119
61.75        0.008820287
63.25        0.002205072
64.75        0.002205072
66.25        0.004410143
67.75        0.002205072
69.25        0.005145167
70.75        0.003675119
72.25        0.004410143
73.75        0.005145167
75.25        0.007350239
76.75        0.006615215
78.25        0.002205072
79.75        0.001470048
81.25        0.005145167
82.75        0.002940096
84.25        0.005145167
85.75        0.005880191
87.25        0.002205072
88.75        0.006615215
90.25        0.003675119
91.75        0.006615215
93.25        0.002205072
94.75        0.004410143
96.25        0.003675119
97.75        0.002940096
99.25        0.002205072
100.75        0.002940096
102.25        0.002205072
103.75        0.001470048
105.25        0.002940096
106.75        0.004410143
108.25        0.005880191
109.75        0.004410143
111.25        0.002205072
112.75        0.005145167
114.25        0.005145167
115.75        0.004410143
117.25        0.005145167
118.75        0.002940096
120.25        0.003675119
121.75        0.002940096
123.25        0.004410143
124.75        0.003675119
126.25        0.002940096
127.75        0.004410143
129.25        0.002940096
130.75        0.005145167
132.25        0.000735024
133.75        0.002205072
135.25        0.003675119
136.75        0.001470048
138.25        0
139.75        0.002940096
141.25        0.000735024
142.75        0.002940096
144.25        0.002205072
145.75        0.001470048
147.25        0.001470048
148.75        0.002940096
150.25        0.002940096
151.75        0.000735024
153.25        0.003675119
154.75        0.002205072
156.25        0.003675119
157.75        0.001470048
159.25        0.001470048
160.75        0.001470048
162.25        0.000735024
163.75        0
165.25        0.001470048
166.75        0.000735024
168.25        0
169.75        0
171.25        0
172.75        0
174.25        0.000735024
175.75        0.000735024
177.25        0
178.75        0.002205072
180.25        0.000735024
181.75        0.001470048
183.25        0
184.75        0.001470048
186.25        0
187.75        0
189.25        0.001470048
190.75        0
192.25        0
193.75        0.000735024
195.25        0.000735024
196.75        0
198.25        0
199.75        0
201.25        0.000735024
202.75        0.000735024
优化算法: 改进差分进化算法
计算结束原因: 达到最大迭代数
均方差(RMSE): 0.00382267882193363
残差平方和(SSE): 0.00191428641221144
相关系数(R): 0.82641411102412
相关系数之平方(R^2): 0.682960282899787
决定系数(DC): 0.666146324132772
卡方系数(Chi-Square): 0.124130918156294
F统计(F-Statistic): 277.889083739702

参数        最佳估算
----------        -------------
a                  16.6356296377931
b                  -37.5941445091556
c                  -5.61028882182676
以下分别是这两种拟合的图片


谁帮着给跑跑哇!1stopt 1.5版本的伤不起呀!结果很不理想额!

麻烦大神帮看一下 1stopt 拟合结果很不满意呀,帮忙跑一下
3.png


麻烦大神帮看一下 1stopt 拟合结果很不满意呀,帮忙跑一下-1
2.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

男儿当自强
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
缘牵一线624: 金币+5 2014-12-08 18:51:21
第二个的结果
"type your title here"

模型公式: y = (((x-c)*a*2.506628275)^(-1))*exp(0-(log(x-c)-b)^2/(2*a^2))
          y = (((x-7.70681005785942)*0.969812630906616*2.506628275)^(-1))*exp(0-(log(x-7.70681005785942)-2.24279570150072)^2/1.881073078)

迭代数: 30
计算用时(时:分:秒:微秒): 00:00:09:411
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE):0.00344030267241807
残差平方和(SSE):0.00155047440459795
相关系数(R): 0.873104529728479
相关系数之平方(R^2): 0.762311519832388
确定系数(DC): 0.729595542228662
卡方系数(Chi-Square): 0.124819633174672
F统计(F-Statistic): 236.682488622729

参数        最佳估算
----------        -------------
a        0.969812630906616
b        2.24279570150072
c        7.70681005785942
3楼2014-12-08 09:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

xxz903

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
缘牵一线624: 金币+5 2014-12-08 18:51:16
第一个的结果
"type your title here"

模型公式: y = ((x*a*2.506628275)^(-1))*exp(0-(log(x)-b)^2/(2*a^2))
          y = ((x*1.02364232986155*2.506628275)^(-1))*exp(0-(log(x)-2.00725248874869)^2/2.095687239)

迭代数: 17
计算用时(时:分:秒:微秒): 00:00:02:577
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE):0.00378409384815714
残差平方和(SSE):0.00187583697896755
相关系数(R): 0.82822321819115
相关系数之平方(R^2): 0.685953699150905
确定系数(DC): 0.672851947983833
卡方系数(Chi-Square): 0.121745751666025
F统计(F-Statistic): 394.31688253985

参数        最佳估算
----------        -------------
a        1.02364232986155
b        2.00725248874869

====== 结果输出 ======
2楼2014-12-08 09:09:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

缘牵一线624

木虫 (小有名气)

引用回帖:
3楼: Originally posted by xxz903 at 2014-12-08 09:10:11
第二个的结果
"type your title here"

模型公式: y = (((x-c)*a*2.506628275)^(-1))*exp(0-(log(x-c)-b)^2/(2*a^2))
          y = (((x-7.70681005785942)*0.969812630906616*2.506628275)^(-1))*e ...

非常感谢xxz903,可以把你跑的结果图给截个吗?  我这边重新绘制了一下,好像拟合的结果不是很好;还有可以问一下你用的是哪个版本的吗?
男儿当自强
4楼2014-12-08 18:50:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见