function?[Fv0Avw]=yinzi1（xnumber）
%此程序是從因子分析直接做的算法
%v0為部分主成分系數，即相關系數矩陣的部分特征向量
%F為因子得分
%A為旋轉后的因子載荷矩陣
%w是對因子評價由大到小的排序
%v是對因子評價排序后元素在原來因子評價向量中的位置。

[x]=bxh（x）;%對數據進行標準化處理
[nm]=size（x）;%n為x的行數m為x的列數，即n為樣本數m為屬性個數。
x1=x;

s=x1（1:）‘*x1（1:）/（n-1）;
for?i=2:n
????s=s+x1（i:）‘*x1（i:）/（n-1）;
end

%s為相關系數矩陣。

[vd]=eig（s‘nobalance‘）;

?for?i=1:m
?????e=1/norm（v（:m-i+1））;
?????%e是向量v（:m-i+1）的模長
?????v1（:i）=e.*v（:m-i+1）;
?end
%將特征向量單位化，同時由特征值從大到小排列.

for?i=m:-1:1?
????d1（m-i+1）=d（ii）;
end
d=d1;

%計算各個主成分的貢獻率
for?i=1:m?
????t（1i）=d（i）/sum（d（:））;
????t（2i）=sum（d（1:i））/sum（d（:））;
end

%因子載荷矩陣v0
?for?i=1:m
????if?（t（2i））????????v0（:i）=sqrt（d（i））.*v1（:i）;
????????t（2i）=z（i）;%向量z是累計貢獻率前幾個
????end???????
?end
?
?A=v0;
?%對A施行方差最大正交變換（旋轉）
?[nx_Any_A]=size（A）;%nx_A為A的行數n