24小时热门版块排行榜    

查看: 2815  |  回复: 12

670442830

木虫 (正式写手)

[求助] 最小二乘拟合 初值 运算时间

我在用mathcad做参数拟合   用最小二乘法(L-M法),但是得到的结果很差,改变初值对结果的影响挺大,但又不知道怎么确定合适的初值,请教各位,怎么设定合适的初值呢?
PS:我现在做了一组,有198个数据点,42个参数,方程形式也挺复杂,拟合时运算时间在8个小时左右(电脑用的台式机,配置还不错,四核的),这样是正常的吗?
非常感谢各位
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

670442830(金币+10): 2011-12-06 10:49:24
拟合目前最强最好用的当属1stOpt了,不用猜初值。
42个参数,对非线性问题来说很难保证最优解。运算8个小时是长了点。
2楼2011-12-06 10:45:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

670442830

木虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2011-12-06 10:45:24:
拟合目前最强最好用的当属1stOpt了,不用猜初值。
42个参数,对非线性问题来说很难保证最优解。运算8个小时是长了点。

1stOpt不设初值能得到全局最优解吗?
参数个数对于非线性拟合的结果也有影响,越多越难找到最优解?
运算时间怎么减小呢?用的是mathcad自带的Minerr函数解的
3楼2011-12-06 10:49:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fiw

木虫 (正式写手)

【答案】应助回帖

670442830(金币+8): 2011-12-06 15:03:53
你可以先确定你的每个数据的最大测定相对误差,如:5%,然后再以此为基础进行优化。
4楼2011-12-06 12:15:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher

【答案】应助回帖

670442830(金币+8): 额 拟合软件很多,但是mathcad基本学会了,再去学其它的。。。 2011-12-06 21:15:19
多参数可以考虑monte carlo
5楼2011-12-06 20:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

apple0905016

银虫 (正式写手)

【答案】应助回帖

670442830(金币+8): 那还是算了 2011-12-07 08:14:55
1stopt不用设初解就能计算出来,很好用,如果是单方程的模拟免费下载1.5版本的就可以用,如果是方程组的话就需要购买高级版本了
6楼2011-12-06 21:38:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iliketg

金虫 (正式写手)

【答案】应助回帖

670442830(金币+10): 2011-12-07 08:16:06
42个参数!不知是谁说的,只要4个参数就可以得到任何一个假设。不是拟合的问题,是你的模型或者说方程不合适。
7楼2011-12-07 05:56:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhdgzhdg

至尊木虫 (著名写手)

【答案】应助回帖

670442830(金币+15): 可有的就拟合的很好,平均相对误差0.3%,那是不是要换拟合方法?拟牛顿?共轭梯度? 2011-12-07 08:17:54
42个参数?先确定是不是分组拟合,一般8个参数是L-M的极限,初值影响大的可能是无解。
8楼2011-12-07 06:22:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhdgzhdg

至尊木虫 (著名写手)

【答案】应助回帖

670442830(金币+15): 恩 好的 谢谢 2011-12-08 08:03:12
拟牛顿和共轭梯度都比不上L-M,还是根据基本常识先分组再拟合,否则数学家也能做化工的题了。
9楼2011-12-08 07:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

670442830

木虫 (正式写手)

引用回帖:
9楼: Originally posted by zhdgzhdg at 2011-12-08 07:14:43:
拟牛顿和共轭梯度都比不上L-M,还是根据基本常识先分组再拟合,否则数学家也能做化工的题了。

整体拟合是不是要比分组拟合的结果好一些?
10楼2011-12-08 08:28:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 670442830 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 11/550 2026-11-09 01:35 by l48hFd8Nwae2
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 01:31 by l48hFd8Nwae2
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 01:24 by l48hFd8Nwae2
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 2rKfVQg2jQww 2026-11-06 11/550 2026-11-09 00:26 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
请填处理意见