path=('G');
filelist=xlsread('G:\AMLmutated.xls');
len =length(filelist);
figure;
% 读取文件
for i = 1:len
filename = [path filelist(i),AMLmutated];
date = filename(15:20);
AML _ data(i,:) = xlsread(filename);%读取
min_val = min(sst_data(i,:));
max_val = max(sst_data(i,:));
pic = reshape(sst_data(i,:),360,180);
pic = rot90(pic);
subplot(4,3,i),imshow(pic,[min_val max_val]),title(date);
colormap;
%???? % 保存数据图像
%???? pic = reshape(sst_data(i,:),360,180);
%???? pic = rot90(pic);???
%???? strI = int2str(i);
%???? save_path = [savepath strI '.tif'];
%???? imwrite(pic,save_path);
end
% 取得样本数据
X = AML_data';
tempX = X;
% 取样本大小:维度m，样本数n
[m n] =size(X);
% 求各样本平均值
meanVal =mean(X);
% 样本矩阵中心化
X = X -repmat(meanVal,64800,1);
% 计算协方差
S = X' * X ./ (m - 1); %".*"表示矩阵元素对应相乘
% 计算特征值eg和特征向量Ev
[Ev eg] =eig(S);
% 由大到小排列对应
Ev = fliplr(Ev);
eg = fliplr(eg);
% % 特征向量转置
% EvT = Ev';
% 得到新的成分
Y = tempX * Ev;
% 将各个成分进行输出
figure;
for i = 1:n
outpic = Y(:,i);
min_num = min(outpic);
max_num = max(outpic);
outpic =reshape(outpic,360,180);
outpic =rot90(outpic);
%outtitle = ['特征值=' num2str(latent(i))];
subplot(4,3,i),imshow(outpic,[min_num max_num]),title(outtitle);%显示，并输出特征值
colorbar;
end
