clcclear
load?gj.txt???%把原始數據保存在純文本文件gj.txt中
gj=zscore（gj）;?%數據標準化
r=corrcoef（gj）;??%計算相關系數矩陣
%下面利用相關系數矩陣進行主成分分析，vec1的列為r的特征向量，即主成分的系數
[vec1lamdarate]=pcacov（r）?%lamda為r的特征值，rate為各個主成分的貢獻率
f=repmat（sign（sum（vec1））size（vec11）1）;%構造與vec1同維數的元素為±1的矩陣
vec2=vec1.*f??%修改特征向量的正負號，使得每個特征向量的分量和為正
%?num=4;??%num為選取的主成分的個數
%?df=gj*vec2（:1:num）;??%計算各個主成分的得分
%?tf=df*rate（1:num）/100;?%計算綜合得分
%?[stfind]=sort（tf‘descend‘）;??%把得分按照從高到低的次序排列
%?stf=stf‘?ind=ind‘

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????628??2015-08-06?11:20??PCA\Education.m

?????文件???????1567??2011-04-18?14:36??PCA\gj.txt

?????文件???????1631??2015-08-09?18:13??PCA\Hald.m

?????文件????????224??2011-04-18?14:39??PCA\sn.txt

?????目錄??????????0??2015-08-10?09:48??PCA

-----------?---------??----------?-----??----

?????????????????4050????????????????????5