24小时热门版块排行榜    

查看: 3471  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dr1581

银虫 (小有名气)

[交流] origin8.5 自定义方程的非线性拟合问题 已有2人参与

求助,origin8.5非线性拟合

希望拟合方程q=(q0*b*p)/(1+b*p),其中q和p分别是因变量和自变量,q0和b是想要拟合出来的参数。

从tools-fitting function buildings一路进行设置,对于parameter选的是q0和b,但这时不清楚对于derived parameters怎么选?选这个有什么意义?

之后,我选了derived parameters,用的Q0和B,最后赋值让QO=q0, B=b,但这样最后通不过,说不能用一个derived parameters去计算另一个derived parameters。

貌似8.5和8的拟合方式相差挺大的,也看不到什么编程窗口之类的。

请大神们指点!
回复此楼

» 猜你喜欢

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

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

戴钢盔的猪头

木虫 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
2楼说的对,从数值计算上来说应该先合并取消冗余项。最后计算结果再按照比例自行分配a,b。这样做并没有否定那几项的物理意义。

[ 发自手机版 http://muchong.com/3g ]
6楼2014-01-29 09:42:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

dr1581

银虫 (小有名气)


fegg7502: 金币+1, 鼓励交流 2014-01-29 07:58:59
我试了一下,可以完全不要derived parameter

但有个很大的问题,为什么我拟合一个稍微复杂一点的方程,虽然最后通过设定初值,使得R2很好,有4个9,但是各个参数的标准差非常大?

方程是lnP=ln(N)+(1/T)*(a0+a1*N+a2*N^2+a3*N^3)+(b0+b1*N+b2*N^2),其中自变量是N,因变量是lnP,T是温度设为定值,其余a和b等是需要拟合的参数。
发现拟合结果和初值关系很大,虽然参考文献,但似乎仍然得不到好结果。

拟合结果见上传的图片,具体数据如下。也希望大神们能帮助看看啊!

拟合数据(左边是N,右边是lnP)
0.0168175022        0.1040992373
0.0240559948        0.5300625763
0.042118893        1.0781495772
0.0930130501        1.853428025
0.1733642532        2.5156171944
0.3127851987        3.1841662834
0.553265158        3.8766403338
0.772978586        4.3101619448
0.9284523442        4.5576591297
1.1249832132        4.8253540394
1.279010958        5.0096052024
1.42002795        5.1631805076
1.5513584861        5.2953399159
1.6785145274        5.4154317081
1.7962369417        5.520579487
1.9097339578        5.616860205
2.0181969544        5.7049563053
2.124155131        5.7878900204
2.2229793298        5.8627374436
2.3105427549        5.9265853918
2.4011390267        5.991182067
2.4885272481        6.0516817301
2.5730899466        6.1089160391
2.6547437596        6.1629700566
2.7335961087        6.213997927
2.8101883513        6.2628204124
2.8847415508        6.3093537102
2.9571389754        6.3539343141
3.0275268828        6.396456012
3.0963376808        6.437319232
3.1631638031        6.4765799707
3.2284972448        6.5142962088
3.2923311247        6.5507435387
3.3543795218        6.5858436108
3.4146488384        6.6196405772
3.4735168718        6.6520169999
origin8.5 自定义方程的非线性拟合问题
fit.png

[ Last edited by dr1581 on 2014-1-27 at 22:14 ]
2楼2014-01-27 16:57:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
fegg7502: 金币+2, 应助指数+1 2014-01-29 07:59:08
你的公式有问题,存在过拟合现象:

原公式:lnP=ln(N)+(1/T)*(a0+a1*N+a2*N^2+a3*N^3)+(b0+b1*N+b2*N^2);
完全可以简化等同于:
lnP=ln(N)+(1/T)*(a0+a1*N+a2*N^2+a3*N^3);

原公式有完全多组解,简化后解唯一,拟合效果也完全一样。
3楼2014-01-28 18:59:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dr1581

银虫 (小有名气)

引用回帖:
3楼: Originally posted by dingd at 2014-01-28 18:59:29
你的公式有问题,存在过拟合现象:

原公式:lnP=ln(N)+(1/T)*(a0+a1*N+a2*N^2+a3*N^3)+(b0+b1*N+b2*N^2);
完全可以简化等同于:
lnP=ln(N)+(1/T)*(a0+a1*N+a2*N^2+a3*N^3);

原公式有完全多组解,简化后解唯 ...

理论是这样,但公式分为这三大项是有物理意义的,所以ai与bi意义不同,不能简化
4楼2014-01-29 09:27:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 4/200 2026-09-01 06:38 by o3EMxHUAQDtS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:56 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:45 by o3EMxHUAQDtS
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见