小木虫 --- 600万学术达人喜爱的学术科研平台

首页 >> 计算模拟 >>查看话题

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

求大神帮忙用1stopt求一下参数,跪谢~~
数据在附件里面,函数是y=a1*exp(a2*t)*(u/10)^(a3*t+a4)*
有1万多个t,u,x,y的数据,要求的是a1到a6,文献里提到用levenberg marquardt算法

用户评论

一万多个点,计算比较费时间。自己为什么不动动手?有什么问题吗?
记住把公式:y=a1*exp(a2*t)*(u/10)^(a3*t+a4)*
改成:y=a1*exp(a2*t)*(u/10)^(a3*t+a4)*(a6*(x/0.4)^a5+(a6-1)*(x/0.4));

不过,想要得到更加专业的回答,我建议还是下载小木虫APP。不要怕麻烦,APP里面有很多非常专业的资料,还有很多各专业的大牛免费解答各种问题,现在安装,下载资源还免费呢,赶快下载【小木虫APP】试试吧!

2楼的dingd是1stOpt软件的大神。
你的问题就是普通的多元非线性拟合问题,本身没什么难度,用1.5版本的1stOpt就能自己做。
dingd针对你的问题,把在1stOpt中Function的正确写法都给你写出来了。

3楼: Originally posted by 月只蓝 at 2017-04-19 11:56:30
2楼的dingd是1stOpt软件的大神。
你的问题就是普通的多元非线性拟合问题,本身没什么难度,用1.5版本的1stOpt就能自己做。
dingd针对你的问题,把在1stOpt中Function的正确写法都给你写出来了。
1.5版打开没有反应,关都关不了..

4楼: Originally posted by 忧郁的王二 at 2017-04-19 13:33:23
1.5版打开没有反应,关都关不了.....
右键--以管理员身份运行

5楼: Originally posted by 月只蓝 at 2017-04-19 14:28:58
右键--以管理员身份运行...
试过了,也不行,用5.0版的按F9提示trial of autofit2.0,大神帮我算一下可以么

6楼: Originally posted by 忧郁的王二 at 2017-04-19 17:11:37
试过了,也不行,用5.0版的按F9提示trial of autofit2.0,大神帮我算一下可以么
...
已验证过了1.5版本的1stOpt就可以计算,测试的结果见附图。
本身就不是什么复杂复杂的问题。
1万多组数据,很费时间,自己动手算吧。不要用5.0版的,那是4参数限制的试用版,你的问题有6个参数,用1.5的就行了。

附图.png

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

拟合参数
拟合相关系数R =  0.985935282255
marquet算法需要较高精度的初值,需联合遗传算法效果较好,祝好运!

parafit.png

7楼: Originally posted by 月只蓝 at 2017-04-19 22:01:24
已验证过了1.5版本的1stOpt就可以计算,测试的结果见附图。
本身就不是什么复杂复杂的问题。
1万多组数据,很费时间,自己动手算吧。不要用5.0版的,那是4参数限制的试用版,你的问题有6个参数,用1.5的就行了。 ...
另外就本题也想交流下,1stOpt,这款软件在计算本题时,初值是怎么估算的?任意给定初值能得到较准确的可行解么,谢谢!

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
...
谢谢~

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

均方差(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,这款软件在计算本题时,初值是怎么估算的?任意给定初值能得到较准确的可行解么,谢谢!...
不需要软件的操作者提供初值。

使用: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倍,还是“用户终止”,否则时间更长。

16楼: Originally posted by dingd at 2017-04-20 13:37:29
使用:QuickReg = 100;命令:
计算用时(时:分:秒:微秒): 00:00:56:24
优化算法: 通用全局优化算法(UGO2)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 9.50699880953343
残差平方和(SSE): 909072.479173839
...
计算时间: 1.30894166448s

初值猜的好计算会快很多。


小木虫APP

学术科研必备,90%的学术科研者都在使用

关于求大神帮忙用1stopt计算一下~的相关话题在小木虫APP已经有20位虫友给出了详细回复。

下载小木虫APP, 立即查看学术达人回复

赶快下载小木虫APP查看回复吧!

学术必备
与600万学术达人在线互动!


扫描下载送金币