爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4561|回复: 7

[程序设计] 关于MATLAB最小二乘法的求解

[复制链接]

新浪微博达人勋

发表于 2015-1-30 22:57:45 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
已知3个行向量,X=[x1 x2 x3 ... xn],Y=[y1 y2 y3 ... yn],Z=[z1 z2 z3 ... zn],并且有Z=f(X,Y)。

具体的函数关系 f 未知,但是 f 中仅有X, Y, XY, X^2, Y^2这些项,如何用最小二乘法求函数 f 的具体形式。

应用MATLAB应该如何求解上述问题?求各路大神不吝赐教。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-31 03:46:15 | 显示全部楼层
x=[1097  1284  1502  1394  1303  1555  1917  2051  2111  2286  2311  2003  2435  2625  2948  3155  3372]';
y=[698  872  988  807  738  1025  1316  1539  1561  1765  1762  1960  1902  2013  2446  2736  2825]';
z=[79.3 200.1 163.1 200.1 146.0 177.7 30.9 291.9 160 339.4 159.6 86.3 237.5 107.2 155 201.4 100.2]';
m=[ones(size(x)) x y x.*y x.*x y.*y];
b = regress(z,m)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-1-31 10:36:28 | 显示全部楼层
tanyueye 发表于 2015-1-31 03:46
x=[1097  1284  1502  1394  1303  1555  1917  2051  2111  2286  2311  2003  2435  2625  2948  3155  3 ...

谢谢您的回复。
在您给的回复中,f=a+b*x+c*y+d*xy+e*x^2+g*y^2(a-g为系数)。
但是f的形式可以多种多样,如f=a+b*x+c*y(其他项系数为0)等等。
如何找到最合适的f的形式,使得拟合的Z与实际的Z的残差平方和最小?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-31 14:39:42 | 显示全部楼层
瑜之姬 发表于 2015-1-31 10:36
谢谢您的回复。
在您给的回复中,f=a+b*x+c*y+d*xy+e*x^2+g*y^2(a-g为系数)。
但是f的形式可以多种多样 ...

这个拟合出来后的值已经是最佳了。如果a-g中有系数为0,b会显示的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-31 23:53:09 | 显示全部楼层
x=[1097  1284  1502  1394  1303  1555  1917  2051  2111  2286  2311  2003  2435  2625  2948  3155  3372]';
y=[698  872  988  807  738  1025  1316  1539  1561  1765  1762  1960  1902  2013  2446  2736  2825]';
z=[79.3 200.1 163.1 200.1 146.0 177.7 30.9 291.9 160 339.4 159.6 86.3 237.5 107.2 155 201.4 100.2]';
m=[ones(size(x)) x y x.*y x.*x y.*y];
b = regress(z,m)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-22 20:50:10 | 显示全部楼层
学习了哈,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-6-23 18:00:12 | 显示全部楼层
hengeuli..............
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-6-23 18:00:25 | 显示全部楼层
hengeiili.............
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表