小木虫

当前位置：首页 > 计算模拟 >求大神帮忙用1stopt计算一下~

求大神帮忙用1stopt计算一下~

作者 忧郁的王二: 来源: 小木虫 850 17 举报帖子

求大神帮忙用1stopt求一下参数，跪谢~~

数据在附件里面，函数是y=a1*exp(a2*t)*(u/10)^(a3*t+a4)*[a6(x/0.4)^a5+(a6-1)(x/0.4)]
有1万多个t,u,x,y的数据，要求的是a1到a6，文献里提到用levenberg marquardt算法

下载提醒：APP中打开可直接下载，点击下载

program.txt

今日热帖

精华评论

smqh5207

拟合参数 [ 1.01182669e+04 -2.72687282e+00 4.48463260e-01 -4.54032463e-01 2.63760305e-01 8.44340423e-01]
拟合相关系数R = 0.985935282255

marquet算法需要较高精度的初值，需联合遗传算法效果较好，祝好运！

parafit.png
smqh5207

引用回帖:
7楼: Originally posted by 月只蓝 at 2017-04-19 22:01:24
已验证过了1.5版本的1stOpt就可以计算，测试的结果见附图。
本身就不是什么复杂复杂的问题。
1万多组数据，很费时间，自己动手算吧。不要用5.0版的，那是4参数限制的试用版，你的问题有6个参数，用1.5的就行了。 ...

另外就本题也想交流下，1stOpt，这款软件在计算本题时，初值是怎么估算的？任意给定初值能得到较准确的可行解么，谢谢！
dingd

1stOpt不需要用户提供初值，结果还是稳定的：
均方差(RMSE): 9.51230534458644
残差平方和(SSE): 910087.598958659
相关系数(R): 0.97882457841853
相关系数之平方(R^2): 0.958097555316214
修正R平方(Adj. R^2): 0.958080882703189
确定系数(DC): 0.957968774836335
卡方系数(Chi-Square): 5681.25892281957
F统计(F-Statistic): 46884.6173514059

参数最佳估算
---------- -------------
a1 10118.2684297031
a2 -2.72687288667116
a3 0.448463251251218
a4 -0.454032452015219
a5 0.263760319916379
a6 0.844340406493767
忧郁的王二

引用回帖:
9楼: Originally posted by smqh5207 at 2017-04-20 09:33:21
拟合参数
拟合相关系数R = 0.985935282255

marquet算法需要较高精度的初值，需联合遗传算法效果较好，祝好运！

parafit.png
...

谢谢~
dingd

数据量大计算比较费时，在1stOpt代码中加一句：
QuickReg = 100;
速度能提高至少10倍以上
dingd

均方差(RMSE): 9.50699846584452
残差平方和(SSE): 909072.413445816
相关系数(R): 0.978820787278873
相关系数之平方(R^2): 0.958090133609233
修正R平方(Adj. R^2): 0.958073458043177
确定系数(DC): 0.95801565987347
卡方系数(Chi-Square): 5662.18450319961
F统计(F-Statistic): 46640.8493485193

参数最佳估算
---------- -------------
a1 10118.269535023
a2 -2.7268730137006
a3 0.448463188991172
a4 -0.454032369072656
a5 0.26376030486613
a6 0.844340422810179
月只蓝

引用回帖:
10楼: Originally posted by smqh5207 at 2017-04-20 09:38:28
另外就本题也想交流下，1stOpt，这款软件在计算本题时，初值是怎么估算的？任意给定初值能得到较准确的可行解么，谢谢！...

不需要软件的操作者提供初值。
dingd

使用：QuickReg = 100;命令：
计算用时(时:分:秒:微秒): 00:00:56:24
优化算法: 通用全局优化算法(UGO2)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 9.50699880953343
残差平方和(SSE): 909072.479173839
相关系数(R): 0.97882082045431
相关系数之平方(R^2): 0.958090198554848
修正R平方(Adj. R^2): 0.958073523014634
确定系数(DC): 0.958015656837906
卡方系数(Chi-Square): 5662.25483854615
F统计(F-Statistic): 46642.8318346108

参数最佳估算
---------- -------------
a1 10118.269511819
a2 -2.72687301190166
a3 0.448463189159152
a4 -0.454032369309838
a6 0.844340422779724
a5 0.263760304897529

不使用：QuickReg = 100;命令：

计算用时(时:分:秒:微秒): 01:46:17:800
优化算法: 通用全局优化算法(UGO2)
计算结束原因: 用户中止
均方差(RMSE): 9.50699846584454
残差平方和(SSE): 909072.413445819
相关系数(R): 0.978820787342682
相关系数之平方(R^2): 0.958090133734147
修正R平方(Adj. R^2): 0.958073458168141
确定系数(DC): 0.95801565987347
卡方系数(Chi-Square): 5662.18449255182
F统计(F-Statistic): 46640.8499754397

参数最佳估算
---------- -------------
a1 10118.2695438513
a2 -2.72687301413388
a3 0.448463189230861
a4 -0.454032369370618
a6 0.844340422723348
a5 0.263760304969923

结果基本一样，但后者计算时间是前者的100倍，还是“用户终止”，否则时间更长，