MATLAB | 如何用MATLAB绘制花里胡哨的山脊图
2024-04-22 10:45:05
本期推送教大家如何绘制各种样式的山脊图，这里做了一个工具函数用来实现好看的山脊图的绘制，编写不易请多多点赞，大体绘制效果如下：
依旧工具函数放在文末。
教程部分
0 数据准备
数据为多个一维向量放在元胞数组中，大概像这样：
X1=normrnd(2,2,1,50);
X2=[normrnd(4,4,1,50),normrnd(5,2,1,50)];
X3=[normrnd(6,2,1,50),normrnd(8,4,1,50)];
X4=[normrnd(12,1,1,50),normrnd(12,4,1,50)];
X5=[normrnd(10,2,1,50),normrnd(10,4,1,50)];
X6=[normrnd(7,2,1,50),normrnd(7,4,1,50)];
X7=[normrnd(4,2,1,50),normrnd(4,4,1,50)];Data={X1,X2,X3,X4,X5,X6,X7};
1 基础绘制:上色类型
通过设置ColorMode属性设置上色方式，可设置为：
- Order
- X
- GlobalX
- Kdensity
- Qt
Order：每个山脊不同颜色
JP=joyPlot(Data,'ColorMode','Order');
JP=JP.draw();
绘制并添加图例(同时可以通过MedLine设置显示中位线)：
JP=joyPlot(Data,'ColorMode','Order','MedLine','on');
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)
可以通过设置ColorList属性或者使用setPatchColor函数设置颜色，推荐前者：
ColorList属性设置颜色：
newColorList=[0.1059 0.6196 0.46670.8510 0.3725 0.00780.4588 0.4392 0.70200.6529 0.4059 0.32940.9020 0.6706 0.00780.6510 0.4627 0.11370.4000 0.4000 0.4000];
JP=joyPlot(Data,'ColorMode','Order','ColorList',newColorList);
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)
setPatchColor函数设置颜色，只能设置面颜色设置不了线条颜色：
JP=joyPlot(Data,'ColorMode','Order');
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)% 设置山脊颜色
newColorList=[0.1059 0.6196 0.46670.8510 0.3725 0.00780.4588 0.4392 0.70200.6529 0.4059 0.32940.9020 0.6706 0.00780.6510 0.4627 0.11370.4000 0.4000 0.4000];
JP=JP.setPatchColor(newColorList)
当然可以全部设置为同一颜色：
JP=joyPlot(Data,'ColorMode','Order','ColorList',[0,0,.6]);
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl(1))
X 每个山脊分别把X映射为颜色
JP=joyPlot(Data,'ColorMode','X','MedLine','on');
JP=JP.draw();
当然可以添加colorbar:
colorbar()
设置为其他颜色：
JP=JP.setPatchColor(winter);
JP=JP.setPatchColor(colorcube(256));
GlobalX 全局性把X坐标映射为颜色
Data={X1,X2,X3,X4,X5,X6,X7};
JP=joyPlot(Data,'ColorMode','GlobalX','MedLine','on');
JP=JP.draw();colorbar
换成别的颜色：
JP=JP.setPatchColor(colorcube(256));
Kdensity 依靠核密度
Data={X1,X2,X3,X4,X5,X6,X7};
JP=joyPlot(Data,'ColorMode','Kdensity','MedLine','on');
JP=JP.draw();colorbar
换成别的颜色：
JP=JP.setPatchColor(turbo);
Qt 分位染色
默认是0.25,0.75分位染色：
JP=joyPlot(Data,'ColorMode','Qt','MedLine','on');
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)
修改分位线(通过设置Quantiles属性)：
JP=joyPlot(Data,'ColorMode','Qt','MedLine','on','Quantiles',[.1,.9]);
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)
绘制分位线(可将QtLine设置为on):
JP=joyPlot(Data,'ColorMode','Qt','MedLine','on','Quantiles',[.1,.9],'QtLine','on');
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)
值得一提的是，所有染色方式都可以绘制QtLine分位线：
JP=joyPlot(Data,'ColorMode','Order','MedLine','on','Quantiles',[.1,.9],'QtLine','on');
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)
更多分位线和更丰富配色：
JP=joyPlot(Data,'ColorMode','Qt','MedLine','on','Quantiles',[.1,.25,.75,.9],'ColorList',turbo(5),'QtLine','on');
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl,{'0~0.1','0.1~0.25','0.25~0.75','0.75~0.9','0.9~1'})
更改配色：
JP=JP.setPatchColor(bone(6));
2 山脊图间距
通过设置Sep属性设置间距，例如小间距：
JP=joyPlot(Data,'ColorMode','Kdensity','Sep',1/30);
JP=JP.draw();colorbar
间距设置为1/60:
设置大间距(1/3)：
3 绘制线状散点
只有Order类型配色不是黑色：
JP=joyPlot(Data,'ColorMode','Order','Scatter','on');
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl
Data={X1,X2,X3,X4,X5,X6,X7};
JP=joyPlot(Data,'ColorMode','GlobalX','MedLine','on','Scatter','on');
JP=JP.draw();colorbar
4 属性单独设置
每个玩意的属性都可以单独设置
例如：
- setRidgePatch 设置山脊面
- setRidgeLine 设置山脊边缘线
- setMedLine 设置中位线
- setQtLine 设置分位线
- setScatter 设置散点
JP=joyPlot(Data,'ColorMode','Order','Scatter','on','QtLine','on','MedLine','on','Sep',1/5);
JP=JP.draw();legendHdl=JP.getLegendHdl();
legend(legendHdl)
设置第3个山脊颜色和透明度：
JP.setRidgePatch(3,'FaceColor',[0,0,0],'FaceAlpha',.5)
循环改变所有的颜色和透明度：
for i=1:length(Data)JP.setRidgePatch(i,'FaceColor',[1,1,1]./length(Data).*i,'FaceAlpha',.5)
end
其他玩意也都类似！
for i=1:length(Data)JP.setRidgePatch(i,'FaceColor',[1,1,1]./length(Data).*i,'FaceAlpha',.5)JP.setRidgeLine(i,'Color',[0,0,.8],'LineWidth',1)JP.setScatter(i,'Color',[0,0,0,.4])JP.setMedLine(i,'Color',[0,0,.8])JP.setQtLine(i,'Color',[0,0,.8])
end
5 多组数据绘图
举个例子：
X1_1=normrnd(-15,2,1,20);
X1_2=[normrnd(4,4,1,10),normrnd(5,2,1,10)];
X1_3=[normrnd(6,2,1,10),normrnd(8,4,1,10)];
X1_4=[normrnd(12,1,1,10),normrnd(12,4,1,10)];
X1_5=[normrnd(-7,2,1,10),normrnd(2,4,1,10)];
X1_6=[normrnd(-7,2,1,10),normrnd(-7,4,1,10)];
Data1={X1_1,X1_2,X1_3,X1_4,X1_5,X1_6};
X2_1=normrnd(-8,2,1,20);
X2_2=[normrnd(2,4,1,10),normrnd(2,2,1,10)];
X2_3=[normrnd(18,2,1,10),normrnd(18,4,1,10)];
X2_4=[normrnd(18,1,1,10),normrnd(18,4,1,10)];
X2_5=[normrnd(5,2,1,10),normrnd(5,4,1,10)];
X2_6=[normrnd(-20,2,1,10),normrnd(-20,4,1,10)];
Data2={X2_1,X2_2,X2_3,X2_4,X2_5,X2_6};JP1=joyPlot(Data1,'ColorMode','Order','ColorList',[12,165,154]./255,'MedLine','on','Scatter','on');
JP1=JP1.draw();JP2=joyPlot(Data2,'ColorMode','Order','ColorList',[151,220,71]./255,'MedLine','on','Scatter','on');
JP2=JP2.draw();% 设置中位线颜色
for i=1:length(Data1)JP1.setMedLine(i,'Color',[12,165,154]./255)
end
for i=1:length(Data2)JP2.setMedLine(i,'Color',[151,220,71]./255)
end% 绘制图例
legendHdl1=JP1.getLegendHdl();
legendHdl2=JP2.getLegendHdl();
legend([legendHdl1(1),legendHdl2(1)],{'AAAAA','BBBBB'})
工具函数完整代码
classdef joyPlot
% @author : slandarer
% gzh : slandarer随笔
% Zhaoxu Liu / slandarer (2023). joyplot
% (https://www.mathworks.com/matlabcentral/fileexchange/125255-joyplot),
% MATLAB Central File Exchange. 检索来源 2023/2/23.
%
%
% =========================================================================
% 1.0 版本创建记录 ——————2023.02.23
% 1.0.0 Update1 设置 'Order'/'X'/'GlobalX'/'Kdensity'/'Qt'5种渲染模式(ColorMode)
% 1.0.0 Update2 设置 中位线(MedLine)、分位线(QtLine)、竖线散点(Scatter)
% 1.0.0 Update3 设置 Seq函数以调整山脊间隙
% 1.0.0 Update4 设置 各个组件修饰函数：
% + setRidgePatch(obj,n,varargin)
% + setRidgeLine(obj,n,varargin)
% + setMedLine(obj,n,varargin)
% + setQtLine(obj,n,varargin)
% + setScatter(obj,n,varargin)
% 1.0.0 Update5 使用 try...Catch结构较为摆烂的解决clim及caxis的版本问题
% -------------------------------------------------------------------------
% 1.5.0 版本更新记录 ——————2023.02.24
% 1.5.0 Update1 使用surf函数代替fill函数解决小区间竖直方向占比过大渲染错误问题
% 1.5.0 Update2 同步更新setPatchColor函数颜色修改代码
% -------------------------------------------------------------------------propertiesax,arginList={'ColorMode','ColorList','Sep','Scatter','MedLine','Quantiles','QtLine'}ColorMode='Order' % 上色模式'Order'/'X'/'GlobalX'/'Kdensity'/'Qt'ColorListdefaultColorList1=[0.3725 0.2745 0.5647; 0.1137 0.4118 0.5882; 0.2196 0.6510 0.6471; 0.0588 0.5216 0.32940.4510 0.6863 0.2824; 0.9294 0.6784 0.0314; 0.8824 0.4863 0.0196; 0.8000 0.3137 0.24310.5804 0.2039 0.4314; 0.4353 0.2510 0.4392];defaultColorList2=[0.0015 0.0005 0.0139; 0.0143 0.0122 0.0705; 0.0415 0.0323 0.1373; 0.0773 0.0535 0.2088; 0.1179 0.0664 0.2854; 0.1661 0.0678 0.3630; 0.2209 0.0609 0.4276; 0.2758 0.0616 0.4677; 0.3279 0.0755 0.4889; 0.3784 0.0954 0.5001; 0.4284 0.1160 0.5058; 0.4788 0.1357 0.5080; 0.5297 0.1541 0.5070; 0.5814 0.1715 0.5028; 0.6338 0.1882 0.4951; 0.6867 0.2051 0.4836; 0.7395 0.2231 0.4679; 0.7914 0.2438 0.4480; 0.8410 0.2692 0.4245; 0.8861 0.3019 0.3992; 0.9240 0.3441 0.3761; 0.9524 0.3958 0.3617; 0.9715 0.4540 0.3610; 0.9837 0.5147 0.3747; 0.9912 0.5758 0.4003; 0.9954 0.6364 0.4350; 0.9972 0.6964 0.4765; 0.9971 0.7559 0.5234; 0.9956 0.8150 0.5747; 0.9930 0.8739 0.6300; 0.9899 0.9327 0.6886; 0.9871 0.9914 0.7495];defaultColorList3=[255,153,154;220,220,220;153,153,253]./255;Sep=1/8; % 两个山脊间距离Scatter='off'; % 是否绘制竖线状散点MedLine='off';QtLine='off';Quantiles=[.25,.75];QtX,QtYridgeNum,Data,minX,maxX,maxY,XiSet,FSetridgePatchHdl,ridgeLineHdlmedLineHdl,scatterHdl;QtLineHdl;QtLegendHdlendmethodsfunction obj=joyPlot(Data,varargin)obj.Data=Data;obj.ridgeNum=length(obj.Data);% 获取其他数据disp(char([64 97 117 116 104 111 114 32 58 32,...115 108 97 110 100 97 114 101 114]))for i=1:2:(length(varargin)-1)tid=ismember(obj.arginList,varargin{i});if any(tid)obj.(obj.arginList{tid})=varargin{i+1};endendif isempty(intersect(obj.ColorMode,{'Order','X','GlobalX','Kdensity','Qt'}))error('The ColorMode should be one of the following: Order \ X \ GlobalX \ Kdensity \ Qt')endswitch obj.ColorModecase 'Order',obj.ColorList=obj.defaultColorList1;case 'X',obj.ColorList=obj.defaultColorList2;case 'GlobalX',obj.ColorList=obj.defaultColorList2;case 'Kdensity',obj.ColorList=obj.defaultColorList2;case 'Qt',obj.ColorList=obj.defaultColorList3;endfor i=1:2:(length(varargin)-1)tid=ismember(obj.arginList,varargin{i});if any(tid)obj.(obj.arginList{tid})=varargin{i+1};endendobj.minX=min(obj.Data{1});obj.maxX=max(obj.Data{1});for i=1:obj.ridgeNumobj.minX=min(obj.minX,min(obj.Data{i}));obj.maxX=max(obj.maxX,max(obj.Data{i}));endendfunction obj=draw(obj)obj.ax=gca;hold on;obj.ax.LineWidth=1;obj.ax.YTick=(1:obj.ridgeNum).*obj.Sep;obj.ax.FontName='Cambria';obj.ax.FontSize=13;obj.ax.YGrid='on';% obj.ax.Box='on';obj.ax.TickDir='out';tYLabel{obj.ridgeNum}='';for i=1:obj.ridgeNumtYLabel{i}=['Class-',num2str(i)];endobj.ax.YTickLabel=tYLabel;view(2)% 调整初始界面大小fig=obj.ax.Parent;fig.Color=[1,1,1];if max(fig.Position(3:4))<690fig.Position(3:4)=1.2.*fig.Position(3:4);fig.Position(1:2)=fig.Position(1:2)./2;end% 绘制patch图像obj.minX=min(obj.Data{1});obj.maxX=max(obj.Data{1});obj.maxY=0;for i=1:obj.ridgeNumtX=obj.Data{i};tX=tX(:)';[F,Xi]=ksdensity(tX);obj.minX=min(obj.minX,min(Xi));obj.maxX=max(obj.maxX,max(Xi));obj.maxY=max(obj.maxY,max(F));endfor i=obj.ridgeNum:-1:1tX=obj.Data{i};tX=tX(:)';[F,Xi]=ksdensity(tX);OXi=Xi;Xi=linspace(min(Xi),max(Xi),1000);F=interp1(OXi,F,Xi);obj.XiSet{i}=Xi;obj.FSet{i}=F;% 绘制竖线散点tXX=[tX;tX;tX.*nan];tYY=[tX.*0+obj.Sep.*i-obj.Sep./10;tX.*0+obj.Sep.*i-obj.Sep./2.5;tX.*nan];if isequal(obj.ColorMode,'Order')obj.scatterHdl(i)=plot(tXX(:),tYY(:),'Color',[obj.ColorList(mod(i-1,size(obj.ColorList,1))+1,:),.5],'LineWidth',.8,'Visible','off');elseobj.scatterHdl(i)=plot(tXX(:),tYY(:),'Color',[0,0,0,.5],'LineWidth',.8,'Visible','off');endif isequal(obj.Scatter,'on'),set(obj.scatterHdl(i),'Visible','on');end% 计算分位线for j=1:length(obj.Quantiles)obj.QtX(i,j+1)=quantile(tX,obj.Quantiles(j));obj.QtY(i,j)=interp1(Xi,F,quantile(tX,obj.Quantiles(j)));endobj.QtX(i,1)=min(Xi)-inf;obj.QtX(i,length(obj.Quantiles)+2)=max(Xi)+inf;switch obj.ColorModecase 'Order'obj.ridgePatchHdl(i)=fill([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...obj.ColorList(mod(i-1,size(obj.ColorList,1))+1,:),'EdgeColor','none','FaceAlpha',.5);obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...'Color',obj.ColorList(mod(i-1,size(obj.ColorList,1))+1,:),'LineWidth',.8);colormap(obj.ColorList);try caxis([1,obj.ridgeNum]),catch,endtry clim([1,obj.ridgeNum]),catch,endcase 'X'tTi=[Xi(1),Xi,Xi(end)]-min(Xi);tTi=tTi./max(tTi);tT=linspace(0,1,size(obj.ColorList,1));tR=interp1(tT,obj.ColorList(:,1),tTi);tG=interp1(tT,obj.ColorList(:,2),tTi);tB=interp1(tT,obj.ColorList(:,3),tTi);tC=cat(3,[tR;tR],[tG;tG],[tB;tB]);obj.ridgePatchHdl(i)=surf([Xi(1),Xi,Xi(end);Xi(1),Xi,Xi(end)],[0,F.*0,0;0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),[tTi;tTi].*0,...'CData',tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...'Color',[0,0,0,.9],'LineWidth',.8);colormap(obj.ColorList);try caxis([-1,1]),catch,endtry clim([-1,1]),catch,end% fill函数老版本代码% tTi=[Xi(1),Xi,Xi(end),Xi(end:-1:1)]-min(Xi);tTi=tTi./max(tTi);% tT=linspace(0,1,size(obj.ColorList,1));% tC=cat(3,interp1(tT,obj.ColorList(:,1),tTi),interp1(tT,obj.ColorList(:,2),tTi),interp1(tT,obj.ColorList(:,3),tTi));% obj.ridgePatchHdl(i)=fill([Xi(1),Xi,Xi(end),Xi(end:-1:1)],[0,F,0,F.*0]+obj.Sep.*(i).*ones(1,length(F)*2+2),...% tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');% obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...% 'Color',[0,0,0,.9],'LineWidth',.8);% colormap(obj.ColorList);% try caxis([-1,1]),catch,end% try clim([-1,1]),catch,endcase 'GlobalX'tTi=[Xi(1),Xi,Xi(end)]-obj.minX;tTi=tTi./(obj.maxX-obj.minX);tT=linspace(0,1,size(obj.ColorList,1));tR=interp1(tT,obj.ColorList(:,1),tTi);tG=interp1(tT,obj.ColorList(:,2),tTi);tB=interp1(tT,obj.ColorList(:,3),tTi);tC=cat(3,[tR;tR],[tG;tG],[tB;tB]);obj.ridgePatchHdl(i)=surf([Xi(1),Xi,Xi(end);Xi(1),Xi,Xi(end)],[0,F.*0,0;0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),[tTi;tTi].*0,...'CData',tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...'Color',[0,0,0,.9],'LineWidth',.8);colormap(obj.ColorList);try caxis([obj.minX,obj.maxX]),catch,endtry clim([obj.minX,obj.maxX]),catch,end% fill函数老版本代码% tTi=[Xi(1),Xi,Xi(end),Xi(end:-1:1)]-obj.minX;% tTi=tTi./(obj.maxX-obj.minX);% tT=linspace(0,1,size(obj.ColorList,1));% tC=cat(3,interp1(tT,obj.ColorList(:,1),tTi),interp1(tT,obj.ColorList(:,2),tTi),interp1(tT,obj.ColorList(:,3),tTi));% obj.ridgePatchHdl(i)=fill([Xi(1),Xi,Xi(end),Xi(end:-1:1)],[0,F,0,F.*0]+obj.Sep.*(i).*ones(1,length(F)*2+2),...% tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');% obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...% 'Color',[0,0,0,.9],'LineWidth',.8);% colormap(obj.ColorList);% try caxis([obj.minX,obj.maxX]),catch,end% try clim([obj.minX,obj.maxX]),catch,endcase 'Kdensity'tTi=[0,F,0];tTi=tTi./obj.maxY;tT=linspace(0,1,size(obj.ColorList,1));tR=interp1(tT,obj.ColorList(:,1),tTi);tG=interp1(tT,obj.ColorList(:,2),tTi);tB=interp1(tT,obj.ColorList(:,3),tTi);tC=cat(3,[tR;tR],[tG;tG],[tB;tB]);obj.ridgePatchHdl(i)=surf([Xi(1),Xi,Xi(end);Xi(1),Xi,Xi(end)],[0,F.*0,0;0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),[tTi;tTi].*0,...'CData',tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...'Color',[0,0,0,.9],'LineWidth',.8);colormap(obj.ColorList);try caxis([0,obj.maxY]),catch,endtry clim([0,obj.maxY]),catch,end% fill函数老版本代码% tTi=[0,F,0,F(end:-1:1)];% tTi=tTi./obj.maxY;% tT=linspace(0,1,size(obj.ColorList,1));% tC=cat(3,interp1(tT,obj.ColorList(:,1),tTi),interp1(tT,obj.ColorList(:,2),tTi),interp1(tT,obj.ColorList(:,3),tTi));% obj.ridgePatchHdl(i)=fill([Xi(1),Xi,Xi(end),Xi(end:-1:1)],[0,F,0,F.*0]+obj.Sep.*(i).*ones(1,length(F)*2+2),...% tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');% obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...% 'Color',[0,0,0,.9],'LineWidth',.8);% colormap(obj.ColorList);% try caxis([0,obj.maxY]),catch,end% try clim([0,obj.maxY]),catch,endcase 'Qt'tTi=[Xi(1),Xi,Xi(end)];tR=tTi.*0;tG=tTi.*0;tB=tTi.*0;for j=1:size(obj.QtX,2)-1tR(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,1);tG(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,2);tB(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,3);endtC=cat(3,[tR;tR],[tG;tG],[tB;tB]);obj.ridgePatchHdl(i)=surf([tTi;tTi],[0,F.*0,0;0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),[tTi;tTi].*0,...'CData',tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...'Color',[0,0,0,.9],'LineWidth',.8);colormap(obj.ColorList);try caxis([-1,1]),catch,endtry clim([-1,1]),catch,end% fill函数老版本代码% tTi=[Xi(1),Xi,Xi(end),Xi(end:-1:1)];% tR=tTi.*0;tG=tTi.*0;tB=tTi.*0;% for j=1:size(obj.QtX,2)-1% tR(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,1);% tG(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,2);% tB(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,3);% end% tC=cat(3,tR,tG,tB);% obj.ridgePatchHdl(i)=fill([Xi(1),Xi,Xi(end),Xi(end:-1:1)],[0,F,0,F.*0]+obj.Sep.*(i).*ones(1,length(F)*2+2),...% tC,'EdgeColor','none','FaceAlpha',.9,'FaceColor','interp');% obj.ridgeLineHdl(i)=plot([Xi(1),Xi,Xi(end)],[0,F,0]+obj.Sep.*(i).*ones(1,length(F)+2),...% 'Color',[0,0,0,.9],'LineWidth',.8);% colormap(obj.ColorList);% try caxis([-1,1]),catch,end% try clim([-1,1]),catch,endend% 绘制中位线tMedX=median(tX);tMedY=interp1(Xi,F,tMedX);obj.medLineHdl(i)=plot([tMedX,tMedX],[0,tMedY]++obj.Sep.*[i,i],'LineStyle','--','LineWidth',1,'Color',[0,0,0],'Visible','off'); if isequal(obj.MedLine,'on'),set(obj.medLineHdl(i),'Visible','on');end% 绘制分位线tQtY=[obj.QtY(i,:);obj.QtY(i,:).*0;obj.QtY(i,:).*nan]+obj.Sep.*i;tQtX=[obj.QtX(i,2:end-1);obj.QtX(i,2:end-1);obj.QtX(i,2:end-1).*nan];obj.QtLineHdl(i)=plot(tQtX(:),tQtY(:),'LineWidth',1,'Color',[0,0,0,.8],'Visible','off');if isequal(obj.QtLine,'on'),set(obj.QtLineHdl(i),'Visible','on');end% % 绘制25，75分位线% tQt25X=quantile(tX,0.25);% tQt75X=quantile(tX,0.75);% tQt25Y=interp1(Xi,F,tQt25X);% tQt75Y=interp1(Xi,F,tQt75X);% obj.qt25LineHdl(i)=plot([tQt25X,tQt25X],[0,tQt25Y]+obj.Sep.*[i,i],'LineWidth',1,'Color',[1,1,1,.8],'Visible','off');% obj.qt75LineHdl(i)=plot([tQt75X,tQt75X],[0,tQt75Y]+obj.Sep.*[i,i],'LineWidth',1,'Color',[1,1,1,.8],'Visible','off');% if isequal(obj.QT25Line,'on'),set(obj.qt25LineHdl(i),'Visible','on');end% if isequal(obj.QT75Line,'on'),set(obj.qt75LineHdl(i),'Visible','on');endendaxis tightobj.ax.YLim(1)=obj.Sep/2;for i=1:size(obj.QtX,2)-1obj.QtLegendHdl(i)=fill(mean(obj.ax.XLim).*[1,1,1,1],mean(obj.ax.YLim).*[1,1,1,1],...obj.ColorList(mod(i-1,size(obj.ColorList,1))+1,:),'EdgeColor','none','FaceAlpha',.9);endend% 获取绘制图例对象function legendHdl=getLegendHdl(obj)if isequal(obj.ColorMode,'Qt')legendHdl=obj.QtLegendHdl;elselegendHdl=obj.ridgePatchHdl;endend% 颜色重设置function obj=setPatchColor(obj,ColorList)obj.ColorList=ColorList;colormap(obj.ColorList);for i=obj.ridgeNum:-1:1Xi=obj.XiSet{i};F=obj.FSet{i};switch obj.ColorModecase 'Order'set(obj.ridgePatchHdl(i),'FaceColor',obj.ColorList(mod(i-1,size(obj.ColorList,1))+1,:));case 'X'tTi=[Xi(1),Xi,Xi(end)]-min(Xi);tTi=tTi./max(tTi);tT=linspace(0,1,size(obj.ColorList,1));tR=interp1(tT,obj.ColorList(:,1),tTi);tG=interp1(tT,obj.ColorList(:,2),tTi);tB=interp1(tT,obj.ColorList(:,3),tTi);tC=cat(3,[tR;tR],[tG;tG],[tB;tB]);set(obj.ridgePatchHdl(i),'CData',tC);case 'GlobalX'tTi=[Xi(1),Xi,Xi(end)]-obj.minX;tTi=tTi./(obj.maxX-obj.minX);tT=linspace(0,1,size(obj.ColorList,1));tR=interp1(tT,obj.ColorList(:,1),tTi);tG=interp1(tT,obj.ColorList(:,2),tTi);tB=interp1(tT,obj.ColorList(:,3),tTi);tC=cat(3,[tR;tR],[tG;tG],[tB;tB]);set(obj.ridgePatchHdl(i),'CData',tC);case 'Kdensity'tTi=[0,F,0];tTi=tTi./obj.maxY;tT=linspace(0,1,size(obj.ColorList,1));tR=interp1(tT,obj.ColorList(:,1),tTi);tG=interp1(tT,obj.ColorList(:,2),tTi);tB=interp1(tT,obj.ColorList(:,3),tTi);tC=cat(3,[tR;tR],[tG;tG],[tB;tB]);set(obj.ridgePatchHdl(i),'CData',tC);case 'Qt'tTi=[Xi(1),Xi,Xi(end)];tR=tTi.*0;tG=tTi.*0;tB=tTi.*0;for j=1:size(obj.QtX,2)-1tR(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,1);tG(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,2);tB(tTi>=obj.QtX(i,j)&tTi<obj.QtX(i,j+1))=obj.ColorList(mod(j-1,size(obj.ColorList,1))+1,3);endtC=cat(3,[tR;tR],[tG;tG],[tB;tB]);set(obj.ridgePatchHdl(i),'CData',tC);endendfor i=1:size(obj.QtX,2)-1set(obj.QtLegendHdl(i),'FaceColor',obj.ColorList(mod(i-1,size(obj.ColorList,1))+1,:));endend% 设置Patch及Line对象其他属性function setRidgePatch(obj,n,varargin)set(obj.ridgePatchHdl(n),varargin{:})endfunction setRidgeLine(obj,n,varargin)set(obj.ridgeLineHdl(n),varargin{:})end% 设置各个分位线属性function setMedLine(obj,n,varargin)set(obj.medLineHdl(n),varargin{:})endfunction setQtLine(obj,n,varargin)set(obj.QtLineHdl(n),varargin{:})end% 设置scatter属性function setScatter(obj,n,varargin)set(obj.scatterHdl(n),varargin{:})endend
% @author : slandarer
% gzh : slandarer随笔
% Zhaoxu Liu / slandarer (2023). joyplot
% (https://www.mathworks.com/matlabcentral/fileexchange/125255-joyplot),
% MATLAB Central File Exchange. 检索来源 2023/2/23.
end
完
编程不易，希望大家该点赞的点赞，该在看的在看！！
未经允许本代码请勿作商业用途，引用的话可以引用我file exchange上的链接，可使用如下格式：
Zhaoxu Liu / slandarer (2023). joyplot (https://www.mathworks.com/matlabcentral/fileexchange/125255-joyplot), MATLAB Central File Exchange. 检索来源 2023/2/23.
若转载请保留以上file exchange链接及本文链接！！！！!
MATLAB | 如何用MATLAB绘制花里胡哨的山脊图相关推荐
- MATLAB三维绘图(四)绘制特殊的三维图
MATLAB三维绘图(四)绘制特殊的三维图 1.使用pie3绘制三维饼状图,示例: %% 绘制三维饼图 clear all; x=[0.2 0.1 0.25 0.15 0.16]; y=[3 5.5 ...
- Matlab批量拟合曲线并绘制在一张图上
文章目录 简单的拟合 批量拟合并呈现在一张图上 最终结果与完整代码 简单的拟合 在命令行输入 cftool 可以启动工具箱 1.在"选择数据"中导入数据 2.选择拟合方法(matl ...
- 解一元二次方程 matlab,如何用Matlab求一元二次方程式解的个数以及解
我们学习数学经常会遇到一元二次方程式,那么如何用Matlab软件求一元二次方程式解的个数以及各个解的值呢?下面小编给大家分享一下. 工具/材料 Matlab 首先打开Matlab软件,点击左上角的Ne ...
- 在线作图|如何绘制一张山脊图
山脊图(Ridgeline chart) 山脊图(Ridgeline chart)作为可视化图形的一种,可以研究不同群组的数值变量的分布情况,展示不同类别数据在同一因素的不同水平下的分布差异,分布可以 ...
- 技术解析：如何用pyecharts绘制时间轮播图
在前天的文章『用python制作动态图表看全球疫情变化趋势』中,由于篇幅原因,在数据处理与数据可视化相关内容上我们只是简单带过,那么我将以python小小白的角度去还原如何处理数据与数据可视化.本文为 ...
- MATLAB | 如何用MATLAB绘制这样有气泡感的网络图
今天给大家带来一款用来绘制有气泡感的网络图的工具函数,绘制效果如下: 花里胡哨的,气泡大小代表流入流出数据量综合,不同颜色的气泡代表属于不同类,两个气泡之间有连线代表有数据流动,连线透明度代表流动数据 ...
- 优化MATLAB中quiver函数绘制箭头图或矢量图（1）-MATLAB开发
目录 利用MATLAB绘制更好看的向量图 MATLAB向量图函数quiver MATLAB绘制3D箭头 MATLAB绘制圆锤 MATLAB绘制圆柱 箭头的方向与颜色控制 光照参数调节 新quiver效 ...
- 中通滤波matlab,如何用MATLAB将低通改成中通滤波器
程序注释:有疑问的我标注出来了???f1=10;f2=200;fs=1000;????????F1=f1/fs;F2=f2/fs;???????%数字频率(此处应该是2*f1/fs,2*f2/fs)? ...
- tbd matlab,如何用matlab实现偏相关分析
matlab偏最小二乘分析(2009-11-20 16:43:07)转载▼ 标签:杂谈 单因变量function y=pls(pz) [row,col]=size(pz); aver=mean(pz) ...
最新文章
- Spring MVC Controller介绍
- 观察多个线程同时运行|| 查看进程线程的方法——未完待续
- ALV输出无法指定STATUS
- Modelsim10.7中文注释乱码问题
- SpringCloud微服务(03)：Hystrix组件，实现服务熔断
- 基于iOS 10、realm封装的下载器
- Nginx+Tomcat无法正确获取的域名和ip地址
- JavaScript内存那点事
- Linux 下java jdk安装
- incrby redis 最大值_Redis 的 8 大数据类型，写得非常好！
- 人力资源SaaS“一哥”，二十年艰难进化
- JavaScript开发者应懂的33个概念js-33-concepts
- 在VM14中建立Win10虚拟机并实现与宿主机的互联
- 基于OCCT7.4的实验性CAD程序
- 奇异值分解的几何原理
- 【赶紧收藏】福利：Python全国计算机二级等级考试题库免费送！！！！
- 点心云pcdn跑量越来越少问题解决方案
- 能ping通但是不能打开网页问题的解决方法
- C++面向对象程序设计教程
- 手机发射机的功能电路（4）---上变频器、功率放大器