function [mainpca,evals]=pca(data)
Maxtrix=data'*data;
%計(jì)算特征分解
[u,d,v]=svd(Maxtrix);
%計(jì)算特征值對(duì)角線元素
evals=real(diag(d));
%對(duì)角陣元素取實(shí)數(shù)，
mainpca=data*u(:,2:9);
%原數(shù)據(jù)向特征向量的投影，得到主成分
plot(mainpca);
return