clc, clear;

% a = load('E:\a-建模\《数学建模算法与应用》课件资源\数学建模算法与应用\程序及数据\10第10章\gj.txt');

a = xlsread('E:\a-建模\第六轮\结果\附件2.汽泵前置泵部分监测数据.xls', 2,'J3:P3602');

% a = xlsread('E:\a-建模\第六轮\结果\附件3.汽轮机本体部分监测数据.xlsx', 2,'A17:X5937');

% a = xlsread('E:\a-建模\第六轮\结果\附件4.一次风机部分监测数据.xls', 2,'L21:P3411');

b = zscore(a); %标准化

r = corrcoef(b); %计算相关系数矩阵

d = pdist(b', 'correlation'); %计算相关系数导出距离方法

z=linkage(d,'average');  %按类平均法聚类

h = dendrogram(z); %画聚类图

set(h,'Color','k','LineWidth',1.3)  %把聚类图线的颜色改成黑色，线宽加粗

xlabel('类别');

% title('气泵前置部分');

% title('汽轮机本体部分');

title('一次风机部分');

ff = 4; %分成类的数量，先用总数运行观察有几类，在来修改这个值

T = cluster(z, 'maxclust', ff);

for i = 1:ff

    tm=find(T==i);  %求第i类的对象

    tm=reshape(tm,1,length(tm)); %变成行向量

    fprintf('第%d类的有%s\n',i,int2str(tm)); %显示分类结果

end