搜了很多资料,这是其中一个版本。
直接替换Test.mat,运行即可。
% 读取数据
load Test.mat;
% 初始化绘图参数
x = dataSamp(:,1);
y = dataSamp(:,2);
XY = [x,y];
datamin = min(min(x), min(y));
datamax = max(max(x), max(y));
datamin = floor(datamin);
datamax = ceil(datamax);
axismin = datamin;
axismax = datamax;
% 绘制密度图
sizeXYZ = size(XY);
searchR = 1.0;
for i = 1 : sizeXYZ(1)
index_i = find(XY(:, 1) > XY(i, 1) - searchR & XY(:, 1) < XY(i, 1) + searchR ...
& XY(:, 2) > XY(i, 2) - searchR & XY(:, 2) < XY(i, 2) + searchR);
sizeIndexI = size(index_i);
end
c = linspace(1,2,length(x));
scatter(XY(:, 1), XY(:, 2),[], c,'.');
colorbar
colormap('jet')