3.源程序
3.1 cwstd.m 总和标准化法标准化矩阵
%cwstd.m,用总和标准化法标准化矩阵
function std=cwstd<vector>
cwsum=sum<vector,1>; %对列求和
[a,b]=size<vector>; %矩阵大小,a 为行数,b 为列数
for i=1:a
for j=1:b
std<i,j>= vector<i,j>/cwsum<j>;
end
end
3.2 cwfac.m 计算相关系数矩阵
%cwfac.m
function result=cwfac<vector>;
fprintf<'相关系数矩阵:\n'>
std=CORRCOEF<vector> %计算相关系数矩阵
fprintf<'特征向量<vec>及特征值<val>:\n'>
[vec,val]=eig<std> %求特征值<val>及特征向量<vec>
newval=diag<val> ;
[y,i]=sort<newval> ; %对特征根进行排序,y 为排序结果,i 为索引