24小时热门版块排行榜    

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

忧郁的王二

新虫 (小有名气)

[求助] 求大神帮忙用1stopt计算一下~已有3人参与

求大神帮忙用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算法
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : program.txt
  • 2017-04-19 10:52:14, 253.2 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
忧郁的王二(月只蓝代发): 金币+15, 耗费大量时间,热心应助! 2017-04-20 13:45:17
月只蓝: 金币+10, 感谢应助。本部分来自于本版块金币库奖励。 2017-04-20 13:46:28
月只蓝: 回帖置顶 2017-04-20 13:47:01
使用: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楼2017-04-20 13:37:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
忧郁的王二: 金币+1, ★★★很有帮助, 没有软件啊...在网上下的破解版用不了 2017-04-19 13:31:44
一万多个点,计算比较费时间。自己为什么不动动手?有什么问题吗?
记住把公式:y=a1*exp(a2*t)*(u/10)^(a3*t+a4)*[a6(x/0.4)^a5+(a6-1)(x/0.4)]
改成:y=a1*exp(a2*t)*(u/10)^(a3*t+a4)*(a6*(x/0.4)^a5+(a6-1)*(x/0.4));
2楼2017-04-19 11:29:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

2楼的dingd是1stOpt软件的大神。
你的问题就是普通的多元非线性拟合问题,本身没什么难度,用1.5版本的1stOpt就能自己做。
dingd针对你的问题,把在1stOpt中Function的正确写法都给你写出来了。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2017-04-19 11:56:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忧郁的王二

新虫 (小有名气)

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

1.5版打开没有反应,关都关不了..
4楼2017-04-19 13:33:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见