%-----------------------------------------------------------------------
% exa090804_svd.m, for example 9.8.4
% to test svd.m;
% 在MATLAB6.1 和 MATLAB5.3下運行該程序，所得到的 U,V矩陣稍有不同，
% 但奇異值是一樣的
%-----------------------------------------------------------------------
clear all;

% 對給定的矩陣 X，對其作奇異值分解；
X=[1 1 1;2 2 1;3 1 3;1 0 1]
y=[1 2 4 3];
y=y';
[U,S,V]=svd(X)
yp=U'*y
cls=V*(yp(1:3)./diag(S));
els=sum(yp(3+1:4).^2);
cls
els