爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9268|回复: 9

[程序设计] 全最小二乘法的求解程序 基于SVD分解

[复制链接]

新浪微博达人勋

发表于 2011-12-2 23:48:42 | 显示全部楼层 |阅读模式

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

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

x
function [Bhat VEM] = OLSregress(X,y)
% 在观测矩阵X中不必输入常数项1
% 程序会自动给出带有常数项的估计系数Bhat
[n p] = size(X);
dat = [X,y];
mdat = mean([X,y]);  % 均值,用于归零和确定常数项系数
dat = dat - ones(n,1)*mdat;%归零
[U S V] = svd(dat);
A=dat(:,1:p);b=dat(:,p+1);d=diag(S);
B = pinv(A'*A-d(p+1)^2*eye(p))*A'*b;
Bhat = [mdat*[-B;1];B];%正交最小二乘估计系数
VEM = sum(d(1:p))/sum(d);% 线性模型能解释的方差(VEM)

评分

参与人数 2金钱 +18 贡献 +1 收起 理由
topmad + 8 + 1 很给力!
言深深 + 10 matlab功能好强大的说

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-12-3 10:00:26 | 显示全部楼层
飞天也是matlab中人啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-9-5 09:06:33 | 显示全部楼层
全最小二乘是什么意思?和整体最小二乘是一个样子的东西吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-8 05:38:10 | 显示全部楼层
很厉害,很强大!!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-13 06:21:49 | 显示全部楼层
很不错!!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-29 02:37:33 | 显示全部楼层
可以举个例子详细说明一下吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-15 15:29:23 | 显示全部楼层
谢谢楼主
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-20 14:04:10 | 显示全部楼层
不错 用了试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-11 11:15:05 | 显示全部楼层
多谢楼主!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-1-1 19:46:15 | 显示全部楼层
超人 发表于 2012-9-5 09:06
全最小二乘是什么意思?和整体最小二乘是一个样子的东西吗?

同问:全最小二乘和整体最小二乘的区别?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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