【答案】应助回帖

学之道，贵以专！

11楼2018-05-17 21:16:35

小笨熊999

木虫 (小有名气)

应助: 3 (幼儿园)
金币: 2830.2
散金: 20
红花: 1
帖子: 162
在线: 16.8小时
虫号: 1657623
注册: 2012-03-02
性别: GG
专业: 环境污染化学

送红花一朵

引用回帖:

10楼: Originally posted by 独孤神宇 at 2018-05-17 20:21:37
模型不够简化，x-(1-x)写成 2x-1

此外，模型与数据匹配度很差，检查一下模型是否正确

可以考虑将模型简化，如里面的对数就可以简化。

谢谢，模型是正确的，请问里面的对数怎么简化呢

加油加油

12楼2018-05-18 09:47:00

小笨熊999

木虫 (小有名气)

应助: 3 (幼儿园)
金币: 2830.2
散金: 20
红花: 1
帖子: 162
在线: 16.8小时
虫号: 1657623
注册: 2012-03-02
性别: GG
专业: 环境污染化学

引用回帖:

请问我这样编辑的拟合公式对不对呢

TIM图片20180518095016.png

加油加油

13楼2018-05-18 09:55:19

独孤神宇

版主 (知名作家)

应助: 488 (硕士)
贵宾: 0.008
金币: 29527.8
散金: 802
红花: 120
沙发: 1
帖子: 5336
在线: 850.2小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

方程里面括号不对等，
简化如 ln(ab)=ln(a)+ln(b)

发自小木虫Android客户端

» 本帖已获得的红花（最新10朵）

小笨熊999

数值计算

14楼2018-05-18 11:10:01

小笨熊999

木虫 (小有名气)

应助: 3 (幼儿园)
金币: 2830.2
散金: 20
红花: 1
帖子: 162
在线: 16.8小时
虫号: 1657623
注册: 2012-03-02
性别: GG
专业: 环境污染化学

送红花一朵

引用回帖:

14楼: Originally posted by 独孤神宇 at 2018-05-18 11:10:01
方程里面括号不对等，
简化如 ln(ab)=ln(a)+ln(b)

谢谢，我已经改了，但是还是拟合不出来啊

加油加油

15楼2018-05-18 17:24:18

独孤神宇

版主 (知名作家)

应助: 488 (硕士)
贵宾: 0.008
金币: 29527.8
散金: 802
红花: 120
沙发: 1
帖子: 5336
在线: 850.2小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

引用回帖:

15楼: Originally posted by 小笨熊999 at 2018-05-18 17:24:18
谢谢，我已经改了，但是还是拟合不出来啊...

我用其他软件试了下，拟合效果很差，确定模型和数据都没有问题吗？

发自小木虫Android客户端

数值计算

16楼2018-05-18 17:52:30

独孤神宇

版主 (知名作家)

应助: 488 (硕士)
贵宾: 0.008
金币: 29527.8
散金: 802
红花: 120
沙发: 1
帖子: 5336
在线: 850.2小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小笨熊999: 金币+50, ★★★★★最佳答案, 谢谢大神！ 2018-05-21 14:51:25

引用回帖:

15楼: Originally posted by 小笨熊999 at 2018-05-18 17:24:18
谢谢，我已经改了，但是还是拟合不出来啊...

这个是数据有问题，因为对数函数里面的数值不允许等于零。将第一个和最后一个去掉可以较好的拟合，当然，你也可以将两端的数据（x）稍微修改一下，如用0.01代替0，0.99代替1等。

y= x*(1-x)*a0+x*(1-x)*(2*x-1)*a1+(1-x)*ln(9.82E-13*(1-x))+x*ln(2.00E-12*x)

0.01             -27.649626
0.13             -27.78552162
0.24             -27.75368809
0.34             -27.54950609
0.44             -27.42658152
0.53             -27.32707214
0.62             -27.27793132
0.71             -27.07782936
0.81             -27.03468757
0.9             -27.02121764
0.99             -26.9374842

你可以将下面的参数作为初始值代入origin里面拟合。

Iterations: 34
Elapsed Time (Hr:Min:Sec:Msec): 00:00:00:331
Algorithms: Universal Global Optimation(UGO1)
Stop Reason: Convergence tolerance reached
Root of Mean Square Error (RMSE): 0.034174566384315
Residual Sum of Squares (SSE): 0.0128469108631155
Correlation Coef. (R): 0.993304080074849
R-Square: 0.986652995493342
Determination Coef. (DC): 0.986413503726273
Chi-Square: -0.000469246867385238
F-Statistic: 662.422439065742
Parameters Best Estimate
---------- -------------
a0 2.47759090789504
a1 1.40000733107936
====== Output Results =====
No Observed y Calculated y
1 -27.649626 -27.6871282305372
2 -27.78552162 -27.7800580250749
3 -27.75368809 -27.7104256520963
4 -27.54950609 -27.5929350997321
5 -27.42658152 -27.4530549973461
6 -27.32707214 -27.3254438707016
7 -27.27793132 -27.2093538669554
8 -27.07782936 -27.115100049893
9 -27.03468757 -27.0443588785141
10 -27.02121764 -27.0103043143305
11 -26.9374842 -26.962877560121

TIM截图20180518184723.jpg

赞一下(1人)

» 本帖已获得的红花（最新10朵）

小笨熊999

数值计算

17楼2018-05-18 18:48:19

小笨熊999

木虫 (小有名气)

应助: 3 (幼儿园)
金币: 2830.2
散金: 20
红花: 1
帖子: 162
在线: 16.8小时
虫号: 1657623
注册: 2012-03-02
性别: GG
专业: 环境污染化学

送红花一朵

引用回帖:

17楼: Originally posted by 独孤神宇 at 2018-05-18 18:48:19
这个是数据有问题，因为对数函数里面的数值不允许等于零。将第一个和最后一个去掉可以较好的拟合，当然，你也可以将两端的数据（x）稍微修改一下，如用0.01代替0，0.99代替1等。

y= x*(1-x)*a0+x*(1-x)*(2 ...

问题已解决，真是太感谢了！

加油加油

18楼2018-05-21 14:47:46