24小时热门版块排行榜    

查看: 588  |  回复: 1
【悬赏金币】回答本帖问题,作者xbaw3884将赠送您 200 个金币

xbaw3884

新虫 (小有名气)

[求助] 对二元吸附所应用到的统计物理模型方程组拟合 已有1人参与

用Matlab的lsqcurvefit拟合下列方程组:
y1=p1*p2*(x1/p4)^p1/(1+(x1/p4)^p1+(x2/p5)^p3)
y2=p3*p6*(x2/p5)^p3/(1+(x1/p4)^p1)+(x2/p5)^p3)

其中p1-p6是需要拟合的参数,x1,x2是自变量,y1,y2是因变量
问题:1、拟合出来的R2值很低,甚至出现负值,点都不在线上
          2、参数p4和p5一开始拟合出来的数值居然为0,后来调整上下界值,拟合出来的数值结果一直与上界值一致。
          3、lsqcurvefit stopped because the final change in the sum of squares relative to its initial value is less than the default value of the function tolerance.     拟合结果还出现上述这一段文字。
附上数据   x1         y1        x2        y2
               1.30     21.75    1.41     21.48
               3.49    41.28     9.23     26.92
               20.19  49.53     22.48   43.80
               29.58  51.05     32.34   44.16
               39.07  52.32     39.27   51.82
               59.16  52.10     59.21   51.98
               67.68  55.81     69.01   52.48
               77.21  56.98     78.89   52.78
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
第二个方程 多了一个 括号 ?

可能是初始值不合理吧。


均方差(RMSE): 3.59598484490288
残差平方和(SSR): 206.89771207634
相关系数(R): 0.945312412116254
相关系数之平方(R^2): 0.89361555650105
修正R平方(Adj. R^2): 0.752465133149895
确定系数(DC): 0.89356839016106
F统计(F-Statistic): 2.33765227718764

参数                  最佳估算
--------------------        -------------
p1        0.166580875951588
p2        1329780191.67408
p4        2.29622451944733E41
p5        437439243.443519
p3        0.251291588079788
p6        11141.0056564802
数值计算
2楼2021-03-09 17:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xbaw3884 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +46 医学老男孩 2026-08-13 103/5150 2026-08-18 15:33 by dcqxinyang
[基金申请] 明天放榜? +3 Shxjjxjkx 2026-08-18 3/150 2026-08-18 15:32 by a441914426
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 5/250 2026-08-18 14:55 by LNP@mRNA
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见