24小时热门版块排行榜    

查看: 2071  |  回复: 7

wnw1123

金虫 (小有名气)

[求助] 求助用维里方程拟合一组数据...拜谢

如题,求各位高手拟合一组数据...小弟我做好只能做到R2=97%...数据见附件
方程如下x/a0*exp(2*a1*x+3/2*a2*x^2+4/3*a3*x^3)
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 新建文本文档(2).txt
  • 2013-03-07 17:56:03, 3.21 K

» 猜你喜欢

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

咖啡或茶

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-03-08 11:12:44
wnw1123: 金币+2 2013-03-08 12:25:00
刚帮你算了下,不知道能否满足你想要的结果
结果如下:
均方差(RMSE): 4.49646592879876
残差平方和(SSE): 1192.87414508204
相关系数(R): 0.985507642672706
相关系数之平方(R^2): 0.971225313766314
决定系数(DC): 0.971111187247178
卡方系数(Chi-Square): 26.1667442138614
F统计(F-Statistic): 1923.90778599944

参数        最佳估算
----------        -------------
a0                 0.00399345040996115
a1                 -0.908709046790272
a2                 -0.0453715675800201
a3                 0.191983269996326
2楼2013-03-07 22:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wnw1123

金虫 (小有名气)

引用回帖:
2楼: Originally posted by 咖啡或茶 at 2013-03-07 22:56:16
刚帮你算了下,不知道能否满足你想要的结果
结果如下:
均方差(RMSE): 4.49646592879876
残差平方和(SSE): 1192.87414508204
相关系数(R): 0.985507642672706
相关系数之平方(R^2): 0.971225313766314
决定系 ...

感谢虫友,,但是我希望R2有99%以上,,,97%不够啊..能不能在算下
3楼2013-03-08 08:33:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-03-08 11:12:52
这只能从你自身(数据和模型公式)找原因了!

如果模型公式可以随意,下面的公式倒是可以满足你的精度要求:

y = p1+p2/(p3-x)+p4/(p5-x);

均方差(RMSE): 2.03279388911066
残差平方和(SSE): 243.802808740733
相关系数(R): 0.997043442970899
相关系数之平方(R^2): 0.994095627171264
决定系数(DC): 0.994095627171265
卡方系数(Chi-Square): 3.49142210610364
F统计(F-Statistic): 2286.4409805388

参数        最佳估算
----------        -------------
p1        51.5372438009192
p2        4.9353991161098
p3        -0.0579556591057234
p4        6.41878391035308
p5        2.26151919503377
4楼2013-03-08 09:53:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wnw1123

金虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2013-03-08 09:53:23
这只能从你自身(数据和模型公式)找原因了!

如果模型公式可以随意,下面的公式倒是可以满足你的精度要求:

y = p1+p2/(p3-x)+p4/(p5-x);

均方差(RMSE): 2.03279388911066
残差平方和(SSE): 243.802808 ...

感谢感谢 用了你的方程拟合了类似的数据还是可以,,不知道这个是什么方程 有名字么?我们是搞吸附的..平常会拟合一下吸附等温线,不太会拟合  不知道如何去猜初值..不知道有没有什么拟合的文章资料之类的可以推荐..再次拜谢
5楼2013-03-08 12:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

如果拟合用的多,还是推荐用1stOpt吧,不需要初值,简单易用。
6楼2013-03-08 13:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是大帅哥

新虫 (初入文坛)

楼主你的方程x/a0*exp(2*a1*x+3/2*a2*x^2+4/3*a3*x^3)是咋得到的,我也是做吸附的,但是不知道用亨利定律和维里方程咋得到你的方程,求赐教,谢谢
7楼2013-12-10 17:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

y=x/a0*exp(2*a1*x+3/2*a2*x^2+4/3*a3*x^3)
相关系数(R): 0.999903079449834
相关系数之平方(R^2): 0.99980616829326
决定系数(DC): 0.999714665620785
卡方系数(Chi-Square): 1.06510875292582
F统计(F-Statistic): 294012.535675026

参数        最佳估算
----------        -------------
a0                 0.0253286349834672
a1                 2.04679358825356
a2                 -3.44599653754039
a3                 1.16207169402579
8楼2013-12-10 20:47:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wnw1123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 63GmZBFV29mG 2026-08-29 5/250 2026-08-31 03:42 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 zICmwzsBXjbN 2026-08-29 6/300 2026-08-31 03:09 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 6/300 2026-08-31 02:56 by hZiFeudZyoGR
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-31 02:13 by hZiFeudZyoGR
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-31 01:59 by hZiFeudZyoGR
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 6/300 2026-08-31 01:48 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-31 00:21 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 15/750 2026-08-31 00:10 by hZiFeudZyoGR
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:29 by hZiFeudZyoGR
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 19:01 by hZiFeudZyoGR
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:50 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 18:17 by hZiFeudZyoGR
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见