24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1657  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖

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的回帖
普通回帖

月只蓝

主管区长 (职业作家)

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的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
4楼: Originally posted by 忧郁的王二 at 2017-04-19 13:33:23
1.5版打开没有反应,关都关不了.....

右键--以管理员身份运行
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
5楼2017-04-19 14:28:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忧郁的王二

新虫 (小有名气)

引用回帖:
5楼: Originally posted by 月只蓝 at 2017-04-19 14:28:58
右键--以管理员身份运行...

试过了,也不行,用5.0版的按F9提示trial of autofit2.0,大神帮我算一下可以么

发自小木虫IOS客户端
6楼2017-04-19 17:11:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
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的就行了。
求大神帮忙用1stopt计算一下~
附图.png

» 本帖已获得的红花(最新10朵)

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
7楼2017-04-19 22:01:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忧郁的王二

新虫 (小有名气)

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

非常感谢

发自小木虫IOS客户端
8楼2017-04-19 23:59:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smqh5207

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
拟合参数 [  1.01182669e+04  -2.72687282e+00   4.48463260e-01  -4.54032463e-01   2.63760305e-01   8.44340423e-01]
拟合相关系数R =  0.985935282255

marquet算法需要较高精度的初值,需联合遗传算法效果较好,祝好运!
求大神帮忙用1stopt计算一下~-1
parafit.png

» 本帖已获得的红花(最新10朵)

9楼2017-04-20 09:33:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smqh5207

金虫 (小有名气)

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

另外就本题也想交流下,1stOpt,这款软件在计算本题时,初值是怎么估算的?任意给定初值能得到较准确的可行解么,谢谢!
10楼2017-04-20 09:38:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 忧郁的王二 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +10 无情的小羊 2026-04-09 11/550 2026-04-12 12:30 by wp06
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 考研求调剂 +3 ban班小七 2026-04-11 3/150 2026-04-11 20:48 by may_新宇
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 材料085601调剂 +25 何润采123 2026-04-10 27/1350 2026-04-10 23:17 by Ftglcn90
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 一志愿厦大生物学332求调剂 +10 池池池池池池 2026-04-08 10/500 2026-04-09 17:10 by 独醉梦孤城
[考研] 材料专硕322 +14 哈哈哈吼吼吼哈 2026-04-05 14/700 2026-04-09 13:25 by 5268321
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
信息提示
请填处理意见