24小时热门版块排行榜    

查看: 1972  |  回复: 5

wenzimu

新虫 (初入文坛)

[求助] 用1stopt拟合非线性方程结果与其他软件拟合结果差异大

通过1stopt对同一组数据和同一个模型进行非线性方程拟合发现,结果与其他软件给出的结果差别很大,这是为什么?(数据见附件)请有知道的帮个忙,谢谢~
1stopt拟合结果:
迭代数: 44
计算用时(时:分:秒:微秒): 00:00:07:375
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 5.5467429740041
残差平方和(SSE): 799.925298111262
相关系数(R): 0.973134439390573
相关系数之平方(R^2): 0.946990637128004
决定系数(DC): 0.946959669563737
卡方系数(Chi-Square): 23.4967654780386
F统计(F-Statistic): 428.750206750339

参数        最佳估算
----------        -------------
p1                 -0.0351005548573163
p2                 6.55416721757231E-6
p3                 -4.74438332683165E-9
p4                 1.97731430746392E-5
p5                 6.25715387238158E-9
p6                 4.43626088101405E-12
p7                 0.00134857002935478
p8                 -5.8705113024446E-7


其他软件给出的拟合参数为(文献上看到的):
p1                 -0.123
p2                 8.39E-5
p3                 -2.6E-8
p4                 4.378E-5
p5                 -5.85E-5
p6                 2.08E-11
p7                 6.825E-4
p8                 -1.29E-7
回复此楼

» 本帖附件资源列表

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

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
只有数据没见你的拟合公式啊?一起放上来看看。

另外结果的好坏吧参数值返带回去验证比较一下不就明白哪个好哪个差了。
2楼2012-12-10 15:48:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenzimu

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2012-12-10 15:48:21
只有数据没见你的拟合公式啊?一起放上来看看。

另外结果的好坏吧参数值返带回去验证比较一下不就明白哪个好哪个差了。

不好意思,忘了给公式了~
P=p1*s+p2*s^2+p3*s^3+p4*s*d+p5*s^2*d+p6*s^3*d+p7*s*g+p8*s^2*g

两个模型的结果R2都很高,我只是想知道为什么会有差别;而我自己在matlab和origin里面计算,初值用的是文献中给的,最后出来还是和用1stopt的结果一致,是不是文献中给出的结果不可信还是有其他原因呢?
3楼2012-12-10 18:42:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
wenzimu: 金币+3, ★★★★★最佳答案, 感谢 2012-12-10 20:58:04
csgt0: 金币+2, 谢谢 2012-12-11 10:13:19
就你的数据和模型公式,1stOpt得出的结果是最优的,也稳定唯一。把文献参数结果带入验证,与实际差的太远了,能确信你的数据公式与文献已知无误?
4楼2012-12-10 19:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

抱歉,最后一句是“能确信你的数据公式与文献一致无误?”
5楼2012-12-10 19:51:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenzimu

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by dingd at 2012-12-10 19:51:27
抱歉,最后一句是“能确信你的数据公式与文献一致无误?”

按照您的提醒我也验证了一遍    和您的结果一致,后来返回去看文献,从交代数据情况来看,我认为是相同的,但是否其中有没有给出的明确数据处理或其他相关步骤,我也没找到!   十分感谢dingd的回答~!
6楼2012-12-10 20:12:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wenzimu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 h4CP7TrQR8Lg 2026-08-22 4/200 2026-08-23 20:55 by 0p9y1isJp2jM
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 19:55 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 15:33 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见