2-R型聚类

将相似的属性聚合在一起

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