MATLAB中PCA降维程序

PCA降维程序,可以对训练集和测试集同时降维,也可以对单个矩阵进行降维。P为累计贡献率,可以查看降维后保留了百分之多少的有效信息。

第一个对单独矩阵进行降维
% function [p,trainFeatures] = PCA_DimensionReduct(trainFeatures,dimensionReduct)
% %PCA降维
% [totalTrainNum,~] = size(trainFeatures);
% % [totalTestNum,~] = size(testFeatures);
% [~,score,latent,~] = pca(trainFeatures);
% p=cumsum(latent)/sum(latent);%累积贡献值
% trainFeatures = score(1:totalTrainNum,1:dimensionReduct);
% % testFeatures = score(totalTrainNum+1:totalTrainNum+totalTestNum,1:dimensionReduct);
% end

同时对多个矩阵进行降维
function [p,trainFeatures,testFeatures] = PCA_DimensionReduct(trainFeatures,testFeatures,dimensionReduct)
%PCA降维
[totalTrainNum,~] = size(trainFeatures);
[totalTestNum,~] = size(testFeatures);
[~,score,latent,~] = pca([trainFeatures;testFeatures]);
p=cumsum(latent)./sum(latent);%累积贡献值
trainFeatures = score(1:totalTrainNum,1:dimensionReduct);
testFeatures = score(totalTrainNum+1:totalTrainNum+totalTestNum,1:dimensionReduct);
end
  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,以下是一个MATLAB编写的简单PCA降维程序: 假设我们有一个 $n \times m$ 的数据矩阵 $X$,其 $n$ 是样本数,$m$ 是特征数。我们的目标是使用PCA将 $X$ 降维为 $n \times k$ 的矩阵 $Y$,其 $k < m$。 步骤如下: 1.将数据矩阵 $X$ 的每一列(即每个特征)都减去其均值,使每个特征的均值为 $0$。 ```matlab X_mean = mean(X, 1); % 求每列的均值 X_centered = X - X_mean; % 减去均值 ``` 2.计算 $X$ 的协方差矩阵 $C = \frac{1}{n-1} X_{centered}^T X_{centered}$。 ```matlab C = (1 / (n - 1)) * (X_centered' * X_centered); ``` 3.计算 $C$ 的特征向量和特征值,并按特征值从大到小排序。 ```matlab [V, D] = eig(C); eig_values = diag(D); [eig_values_sorted, eig_values_order] = sort(eig_values, 'descend'); eig_vectors_sorted = V(:, eig_values_order); ``` 4.选择前 $k$ 个特征向量,组成变换矩阵 $P$。 ```matlab P = eig_vectors_sorted(:, 1:k); ``` 5.将数据矩阵 $X$ 乘以变换矩阵 $P$,得到降维后的数据矩阵 $Y$。 ```matlab Y = X_centered * P; ``` 完整的代码如下: ```matlab function Y = pca(X, k) % X: n x m 数据矩阵,n 是样本数,m 是特征数 % k: 降维后的特征数 % Y: n x k 降维后的数据矩阵 % 将每个特征减去其均值,使每个特征的均值为 0 X_mean = mean(X, 1); X_centered = X - X_mean; % 计算协方差矩阵 C = (1 / (n - 1)) * (X_centered' * X_centered); % 计算特征向量和特征值,并按特征值从大到小排序 [V, D] = eig(C); eig_values = diag(D); [eig_values_sorted, eig_values_order] = sort(eig_values, 'descend'); eig_vectors_sorted = V(:, eig_values_order); % 选择前 k 个特征向量,组成变换矩阵 P = eig_vectors_sorted(:, 1:k); % 将数据矩阵 X 乘以变换矩阵 P,得到降维后的数据矩阵 Y Y = X_center

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值