| 查看: 584 | 回复: 2 | ||||
[求助]
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]; |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有292人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
zhang_gang
金虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 2957.8
- 散金: 1948
- 红花: 40
- 帖子: 1296
- 在线: 585.8小时
- 虫号: 2442637
- 注册: 2013-04-30
- 性别: GG
- 专业: 模式识别
2楼2015-07-28 18:00:50
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Carrie_XX: 金币+10, ★★★★★最佳答案, 嗯嗯 明白了 灰常谢谢 2015-07-29 10:17:14
感谢参与，应助指数 +1
Carrie_XX: 金币+10, ★★★★★最佳答案, 嗯嗯 明白了 灰常谢谢 2015-07-29 10:17:14
|
可以求解的。MATLAB有专门用于求解最小二乘问题的函数lsqnonlin。 不过简便起见，用1stopt软件，代码如下： 计算结果： 均方差(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 这两项必然会导致公式过拟合。 |
3楼2015-07-28 19:51:20
投票:
10