Matlab之使用HeatMap实现画热图

clc; 
clear; 
data = [0,40,0,31,0,0,0,0;    
       30,0,40,0,2,0,1,0;       
       0,0,0,0,0,13,0,60;       
       39,0,31,0,1,0,0,0;        
       0,30,0,39,0,1,0,0;        
      1,0,1,0,56,0,12,0;       
      0,2,0,1,0,56,0,12;       
      0,0,0,0,12,0,59,0];
 xvalues = {'CNC1','CNC2','CNC3','CNC4','CNC5','CNC6','CNC7','CNC8'}; 
yvalues = {'CNC1','CNC2','CNC3','CNC4','CNC5','CNC6','CNC7','CNC8'}; 
yvalues=yvalues';
h=HeatMap(data,'Colormap',redbluecmap,'ColumnLabels',yvalues,'RowLabels',xvalues); 
h.Annotate = true;%在热图中显示数据值
 

X=normrnd(0,10,[100000,1]);
Y=normrnd(0,10,[100000,1]);
numbins = 50;
markersize = 50;
marker = '.';
[values, centers] = hist3([X Y], [numbins numbins]);
centers_X = centers{1,1};
centers_Y = centers{1,2};
binsize_X = abs(centers_X(2) - centers_X(1)) / 2;
binsize_Y = abs(centers_Y(2) - centers_Y(1)) / 2;
bins_X = zeros(numbins, 2);
bins_Y = zeros(numbins, 2);
for i = 1:numbins
	bins_X(i, 1) = centers_X(i) - binsize_X;
	bins_X(i, 2) = centers_X(i) + binsize_X;
	bins_Y(i, 1) = centers_Y(i) - binsize_Y;
	bins_Y(i, 2) = centers_Y(i) + binsize_Y;
end
scatter_COL = zeros(length(X), 1);
onepercent = round(length(X) / 100);       
for i = 1:length(X)
        last_lower_X = NaN;
        last_higher_X = NaN;
        id_X = NaN;
        c_X = X(i);
        last_lower_X = find(c_X >= bins_X(:,1));
        if (~isempty(last_lower_X))
            last_lower_X = last_lower_X(end);
        else
            last_higher_X = find(c_X <= bins_X(:,2));
            if (~isempty(last_higher_X))
                last_higher_X = last_higher_X(1);
            end
        end
        if (~isnan(last_lower_X))
            id_X = last_lower_X;
        else
            if (~isnan(last_higher_X))
                id_X = last_higher_X;
            end
        end
        last_lower_Y = NaN;
        last_higher_Y = NaN;
        id_Y = NaN;
        c_Y = Y(i);
        last_lower_Y = find(c_Y >= bins_Y(:,1));
        if (~isempty(last_lower_Y))
            last_lower_Y = last_lower_Y(end);
        else
            last_higher_Y = find(c_Y <= bins_Y(:,2));
            if (~isempty(last_higher_Y))
                last_higher_Y = last_higher_Y(1);
            end
        end
        if (~isnan(last_lower_Y))
            id_Y = last_lower_Y;
        else
            if (~isnan(last_higher_Y))
                id_Y = last_higher_Y;
            end
        end
        scatter_COL(i) = values(id_X, id_Y);
end   
scatter(X, Y, markersize, scatter_COL, marker);
hold on
colorbar;
box on