厦门网站建设 首选猴子网络,设计公司官方网站,网上商城平台开发,免费空间自带域名本期部分实验效果#xff1a; 这期讲一下如果数据重合严重该咋办(overlap)#xff0c;事先说明#xff0c;本文中的绘图均使用一个几行的简单小代码进行了修饰#xff1a;
function defualtAxes
axgca;hold on;box on
ax.XGridon;
ax.YGridon;
ax.XMinorTickon;
ax.YMinor…本期部分实验效果 这期讲一下如果数据重合严重该咋办(overlap)事先说明本文中的绘图均使用一个几行的简单小代码进行了修饰
function defualtAxes
axgca;hold on;box on
ax.XGridon;
ax.YGridon;
ax.XMinorTickon;
ax.YMinorTickon;
ax.LineWidth.8;
ax.GridLineStyle-.;
ax.FontNameCambria;
ax.FontSize12;
end0 数据说明及基础绘图
假设我们随机构建两列数据
% 随便生成散点
PntSet1mvnrnd([2 3],[1 0;0 2],800);
PntSet2mvnrnd([6 7],[1 0;0 2],800);
PntSet3mvnrnd([8 9],[1 0;0 1],800);
PntSet[PntSet1;PntSet2;PntSet3];
XPntSet(:,1);
YPntSet(:,2);scatter函数绘图效果
% 使用scatter绘图
scatter(X,Y,filled,CData,[36,59,66]./255);% 简单修饰可省略
defualtAxes()可以看到数据重叠严重根本看不出哪里最密集以下讲几个简单的解决数据重叠的方法。 1 设置透明度
这个就很简单设置MarkerFaceAlpha即可
% 使用scatter绘图
scatter(X,Y,filled,CData,[36,59,66]./255,MarkerFaceAlpha,.5);% 简单修饰可省略
defualtAxes()2 根据密度设置CData
分割网格计算核密度可以自行调整网格划分格子数然后通过插值计算每个点所处位置核密度并将其映射为颜色
% 横竖分割一百格计算核密度
n100;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 使用scatter绘图
scatter(X,Y,filled,CData,H);% 简单修饰可省略
defualtAxes()当然可以通过colormap设置其他配色
colormap(summer)自己随便弄个配色
CM[0.2700 0 0.33000.2700 0.2300 0.51000.1900 0.4100 0.56000.1200 0.5600 0.55000.2100 0.7200 0.47000.5600 0.8400 0.27000.9900 0.9100 0.1300];
colormap(CM)颜色有点不连续了插一下值
CM[0.2700 0 0.33000.2700 0.2300 0.51000.1900 0.4100 0.56000.1200 0.5600 0.55000.2100 0.7200 0.47000.5600 0.8400 0.27000.9900 0.9100 0.1300];
CMXlinspace(0,1,size(CM,1));
CMXXlinspace(0,1,256);
CM[interp1(CMX,CM(:,1),CMXX,pchip),...interp1(CMX,CM(:,2),CMXX,pchip),...interp1(CMX,CM(:,3),CMXX,pchip)];
colormap(CM)3 等高线
还是上面的核密度计算方法然后直接画为等高线
% 横竖分割一百格计算核密度
n100;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 绘制等高线图及散点图
hold on
scatter(X,Y,1,filled,CData,[36,59,66]./255);
contour(XMesh,YMesh,ZMesh,20,LineWidth,.8)% 简单修饰可省略
defualtAxes()4 等高线填充
还是和上面几乎一样
% 横竖分割一百格计算核密度
n100;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 绘制等高线填充图
hold on
contourf(XMesh,YMesh,ZMesh,15,EdgeColor,none)我们设置小于一定值就不画了这里设置为1e-3可自行调整
% 横竖分割一百格计算核密度
n100;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 绘制等高线填充图
hold on
levelslinspace(1e-3,max(max(H)),15);
contourf(XMesh,YMesh,ZMesh,levels,EdgeColor,none)依旧改一下配色
colormap(turbo)也可以用自己配色怕乱直接把这部分完整代码放在一起
% 横竖分割一百格计算核密度
n100;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 绘制等高线填充图
hold on
levelslinspace(1e-3,max(max(H)),15);
contourf(XMesh,YMesh,ZMesh,levels,EdgeColor,none)
scatter(X,Y,1,filled,CData,[36,59,66]./255);% 设置colomap
% colormap(turbo)
CM[0.2700 0 0.33000.2700 0.2300 0.51000.1900 0.4100 0.56000.1200 0.5600 0.55000.2100 0.7200 0.47000.5600 0.8400 0.27000.9900 0.9100 0.1300];
CMXlinspace(0,1,size(CM,1));
CMXXlinspace(0,1,256);
CM[interp1(CMX,CM(:,1),CMXX,pchip),...interp1(CMX,CM(:,2),CMXX,pchip),...interp1(CMX,CM(:,3),CMXX,pchip)];
colormap(CM)% 简单修饰可省略
defualtAxes()5 rug图
就是边缘加一些竖线状散点就有点像地毯的边缘
% 使用scatter绘图
hold on
scatter(X,Y,10,filled,CData,[36,59,66]./255);% 绘制边际线条状散点(rug图)
axgca;
XLimax.XLim;YLimax.YLim;
XX(:);YY(:);
LXX[X;X;X.*nan];
LXY[Y.*0YLim(1);Y.*0YLim(1)(diff(YLim))/20;Y.*nan];
plot(LXX(:),LXY(:),Color,[[36,59,66]./255,.3]);LYY[Y;Y;Y.*nan];
LYX[X.*0XLim(1);X.*0XLim(1)(diff(XLim))/20;X.*nan];
plot(LYX(:),LYY(:),Color,[[36,59,66]./255,.3]);% 简单修饰可省略
defualtAxes()6 分bin图
这里横竖都分为30块
% 分X,Y30块的分bin图
binscatter(X,Y,[30 30])
colorbar% 简单修饰可省略
defualtAxes()7 柱状图
在分bin图格子里画一些柱状图
% 绘制散点图及柱状图
hold on
bcHdlbinscatter(X,Y,[20,20],Visible,off);
scatter(X,Y,1,filled,CData,[36,59,66]./255);
XMean(bcHdl.XBinEdges(1:end-1)bcHdl.XBinEdges(2:end))./2;
YMean(bcHdl.YBinEdges(1:end-1)bcHdl.YBinEdges(2:end))./2;
XSepdiff(bcHdl.XBinEdges(1:2));
YSepdiff(bcHdl.YBinEdges(1:2));
for i1:size(bcHdl.Values,1)for j1:size(bcHdl.Values,2)fill([-1,-1,1,1].*XSep./3XMean(i),...[1,0,0,1].*YSep.*bcHdl.Values(i,j)./max(max(bcHdl.Values)).*.95YMean(j),...[36,59,66]./255,FaceAlpha,.9,EdgeColor,none)end
end% 简单修饰可省略
defualtAxes()8 surf曲面
类似分bin图不过这里不是数量统计而是核密度:
% 横竖分割计算核密度
n30;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 绘制surf曲面
hold on
ZMesh(ZMesh1e-3)nan;
surf(XMesh,YMesh,ZMesh,EdgeColor,none);% 加一行[1,1,1]把小数值设置为白色
colormap(parula)% 简单修饰可省略
defualtAxes()9 气泡图
还是类似的不过换成了bubble气泡图
% 横竖分割计算核密度
n30;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 绘制气泡图
ZMesh(ZMesh1e-3)nan;
bubblechart(XMesh(:),YMesh(:),ZMesh(:),ZMesh(:),MarkerEdgeColor,none)
bubblesize([1,12])% 简单修饰可省略
defualtAxes()10 花里胡哨没用的三角剖分
% 横竖分割一百格计算核密度
n100;
XListlinspace(min(X),max(X),n);
YListlinspace(min(Y),max(Y),n);
[XMesh,YMesh]meshgrid(XList,YList);
Fksdensity([X,Y],[XMesh(:),YMesh(:)]);
ZMeshreshape(F,size(XMesh));
Hinterp2(XMesh,YMesh,ZMesh,X,Y);% 没啥用的三角化插值绘图
hold on
DTdelaunay(X,Y);
Z(H(DT(:,1),:)H(DT(:,2),:)H(DT(:,3),:))./3;
trisurf(DT,X,Y,X.*0,CData,Z,EdgeColor,none)% 简单修饰可省略
defualtAxes()
