svd01.m

来自「利用奇异值法（SVD）对载入数据进行主成分分析」· M 代码 · 共 26 行

% s 环境阵，p 特征向量阵 ,A 特征值阵;并假定方差阵 R 为单位阵,X 遥外测差向量,D 为估值
% S 是6n*m的矩阵 
% S'*S为实对称阵
 %输入S以及向量X;
 load('S.mat');% 载入 S
 load('X.mat');
s=S';
c=s*S;
 [U,A,V]=svd(S);
sum=A(1,1)^2;
 for i=2:27
     sum=sum+A(i,i)^2;
 end
SUM=A(1,1)^2;
for i=2:27
    if (1-SUM/sum)<0.00001
        break
    else SUM=SUM+A(i,i)^2;
        N=i;
    end 
end
sum=zeros(27,1);  
for i=1:N
 sum=sum+(1/A(i,i)^2)*V(:,i)*V(:,i)'*s*X;
end
D=sum;