zyy.m

用Matlab实现的简单的二维主成分分析算法

function zyy

x=load('iris.txt');
sig=cov(x);     %得到协方差矩阵
[v d]=eig(sig); %得到特征值对应的特征矢量组成的矩阵v，特征值为对角元素的对角阵d
[m n]=size(v)   %得到矩阵的大小

matr=[v(:,m),v(:,m-1)];
disp('二维主成分为');
y=x*matr
disp('散点图为');
plot(y(1:50,2),y(1:50,1),'ro',y(51:100,2),y(51:100,1),'b+',y(101:150,2),y(101:150,1),'g^');  %绘出散点图