24小时热门版块排行榜    

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

Carrie_XX

铁虫 (初入文坛)

[求助] matlab最小二乘法 已有2人参与

求助:请问在matlab中,可以用最小二乘法求解y=k1*C1^a+k2*C2^b+k3*C1^c中的参数a,b,c,k1,k2,k3吗?
其中,C1=[4.04;4.01;3.88;3.73;3.69;3.33;3.60;2.75;2.89;2.53;1.97;1.79;1.68;1.79;1.22;0.93];
C2=[2.31;2.00;2.19;2.05;2.05;1.95;2.25;1.93;2.17;2.08;1.88;1.64;1.52;1.51;1.28;1.03];
y=[-2.47631;-2.45481;-2.4341;-2.41419;-2.39507;-2.37677;-2.32679;-2.24996;-2.20617;-2.19706;-2.20566;-2.21719;-2.28449;-2.41093;-2.68782;-2.95603];
回复此楼

» 猜你喜欢

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

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
Carrie_XX: 金币+10, ★★★★★最佳答案, 嗯嗯 明白了 灰常谢谢 2015-07-29 10:17:14
可以求解的。MATLAB有专门用于求解最小二乘问题的函数lsqnonlin。
不过简便起见,用1stopt软件,代码如下:
CODE:
Variable C1,C2,y;
Function y=k1*C1^a+k2*C2^b+k3*C1^c;
Data;
4.04000000000000        2.31000000000000        -2.47631000000000
4.01000000000000        2        -2.45481000000000
3.88000000000000        2.19000000000000        -2.43410000000000
3.73000000000000        2.05000000000000        -2.41419000000000
3.69000000000000        2.05000000000000        -2.39507000000000
3.33000000000000        1.95000000000000        -2.37677000000000
3.60000000000000        2.25000000000000        -2.32679000000000
2.75000000000000        1.93000000000000        -2.24996000000000
2.89000000000000        2.17000000000000        -2.20617000000000
2.53000000000000        2.08000000000000        -2.19706000000000
1.97000000000000        1.88000000000000        -2.20566000000000
1.79000000000000        1.64000000000000        -2.21719000000000
1.68000000000000        1.52000000000000        -2.28449000000000
1.79000000000000        1.51000000000000        -2.41093000000000
1.22000000000000        1.28000000000000        -2.68782000000000
0.930000000000000        1.03000000000000        -2.95603000000000

计算结果:
均方差(RMSE): 0.0411939983362445
残差平方和(SSE): 0.0271511279828242
相关系数(R): 0.976851818624353
相关系数之平方(R^2): 0.954239475549706
决定系数(DC): 0.95423943455436
卡方系数(Chi-Square): -0.00565378501945285
F统计(F-Statistic): 291.940549593289

原方程过拟合,参数结果不稳定、不唯一:
参数        最佳估算
----------        -------------
k1                 -3.9810461734405
a                  0.463874856031145
k2                 24.5236491353617
b                  0.0187026831926242
k3                 -23.4576701772098
c                  -0.123433176390609

原公式中 k1*C1^a 和k3*C1^c 这两项必然会导致公式过拟合。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2015-07-28 19:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

zhang_gang

金虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
这个应该不能求吧,多项式用最小二乘方法拟合,多项式的阶数要已知呀。

[ 发自小木虫客户端 ]
2楼2015-07-28 18:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +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 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +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 冼亮淀粉酶
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+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
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见