matlab 主成分分析法代码

用lab实现主成分分析的代码有多种方法。以下是两种常用的方法：

方法一：

X = load('shuju.txt');
z = zscore(X); % 数据标准化
M = cov(z); % 协方差
[V, D = eig(M); % 求出协方差矩阵的特征向量、特征值
d = diag(D); % 取出特征值矩阵的对角线元素（即特征值）
eig1 = sort(d, 'descend'); % 将特征值按从大到小排列
v = fliplr(V); % 重新排列特征向量
S = 0;
i = 0;
while S/sum(eig1) < 0.85
    i = i + 1;
    S = S + eig1(i);
end
NEW = z * v(:, 1:i); % 输出产生的新坐标下的数据
W = 100 * eig1 / sum(eig1);
figure(1)
pareto(W); % 画出贡献率的直方图

方法二：

X = load('shuju.txt');
x = zscore(X); % 标准化
 = princomp(x); % 利用princomp处理矩阵
t % 每一组数据在新坐标下到原点的距离
s = 0;
i = 1;
while s/sum(eig) < 0.85
    s = s + eig(i);
    i = i + 1;
end
NEW = x * coef(:, 1:i-1); % 输出新的数据
figure
pareto(eig/sum(eig)); % 输出贡献率直方图

以上两种方法都可以实现主成分分析算法，可以根据自己的需求选择使用。<span class="em">1</span><span class="em">2</span><span class="em">3</span><span class="em">4</span>