24小时热门版块排行榜    

查看: 1943  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +16 基诺咪客 2026-08-17 20/1000 2026-08-20 12:44 by 201120518022
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 言之凿凿,8月21号(明天)放榜 +12 医学老男孩 2026-08-20 14/700 2026-08-20 10:53 by ran2233
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 125/6250 2026-08-20 09:02 by 医学老男孩
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见