应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 求大神帮忙用1stopt计算一下~
181/2返回列表12下一页
查看: 1431  |  回复: 17
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

忧郁的王二

新虫 (小有名气)

[求助] 求大神帮忙用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

» 猜你喜欢
1楼2017-04-19 10:55:54
已阅   回复此楼   关注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));
赞一下(1人) 回复此楼
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倍，还是“用户终止”，否则时间更长。
赞一下(1人) 回复此楼
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的回帖
相关版块跳转 我要订阅楼主 忧郁的王二 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭