文章详情

短信预约-IT技能 免费直播动态提醒

请输入下面的图形验证码

提交验证

短信预约提醒成功

Matlab怎么绘制散点密度图

2023-06-29 08:12

关注

小编给大家分享一下Matlab怎么绘制散点密度图,相信大部分人都还不怎么了解,因此分享这篇文章给大家参考一下,希望大家阅读完这篇文章后大有收获,下面让我们一起去了解一下吧!

效果

Matlab怎么绘制散点密度图

Matlab怎么绘制散点密度图

Matlab怎么绘制散点密度图

Matlab怎么绘制散点密度图

原理也很简单,通过matlab自带的ksdensity获得网格每一点密度,通过密度拟合曲面,再计算每个数据点对应的概率,并将概率映射到颜色即可
为了怕大家找不到函数这次工具函数放到最前面

1工具函数完整代码

function [CData,h,XMesh,YMesh,ZMesh,colorList]=density2C(X,Y,XList,YList,colorList)[XMesh,YMesh]=meshgrid(XList,YList);XYi=[XMesh(:) YMesh(:)];F=ksdensity([X,Y],XYi);ZMesh=zeros(size(XMesh));ZMesh(1:length(F))=F;h=interp2(XMesh,YMesh,ZMesh,X,Y);if nargin<5colorList=[0.2700         0    0.3300    0.2700    0.2300    0.5100    0.1900    0.4100    0.5600    0.1200    0.5600    0.5500    0.2100    0.7200    0.4700    0.5600    0.8400    0.2700    0.9900    0.9100    0.1300];endcolorFunc=colorFuncFactory(colorList);CData=colorFunc((h-min(h))./(max(h)-min(h)));colorList=colorFunc(linspace(0,1,100)');function colorFunc=colorFuncFactory(colorList)x=(0:size(colorList,1)-1)./(size(colorList,1)-1);y1=colorList(:,1);y2=colorList(:,2);y3=colorList(:,3);colorFunc=@(X)[interp1(x,y1,X,'pchip'),interp1(x,y2,X,'pchip'),interp1(x,y3,X,'pchip')];endend

2参数说明

输入:

输出:

3使用方式

假如编写了如下程序:

PntSet1=mvnrnd([2 3],[1 0;0 2],800);PntSet2=mvnrnd([6 7],[1 0;0 2],800);PntSet3=mvnrnd([8 9],[1 0;0 1],800);PntSet=[PntSet1;PntSet2;PntSet3];scatter(PntSet(:,1),PntSet(:,2),'filled');

结果:

Matlab怎么绘制散点密度图

3.1散点赋色

将上面那段代码改写

PntSet1=mvnrnd([2 3],[1 0;0 2],800);PntSet2=mvnrnd([6 7],[1 0;0 2],800);PntSet3=mvnrnd([8 9],[1 0;0 1],800);PntSet=[PntSet1;PntSet2;PntSet3];CData=density2C(PntSet(:,1),PntSet(:,2),-2:0.1:15,-2:0.1:15);scatter(PntSet(:,1),PntSet(:,2),'filled','CData',CData);

Matlab怎么绘制散点密度图

3.2等高线图

PntSet1=mvnrnd([2 3],[1 0;0 2],800);PntSet2=mvnrnd([6 7],[1 0;0 2],800);PntSet3=mvnrnd([8 9],[1 0;0 1],800);PntSet=[PntSet1;PntSet2;PntSet3];[~,~,XMesh,YMesh,ZMesh,colorList]=density2C(PntSet(:,1),PntSet(:,2),-2:0.1:12,-2:0.1:12);colormap(colorList)contourf(XMesh,YMesh,ZMesh,10)

Matlab怎么绘制散点密度图

3.3带直方图的散点图

PntSet1=mvnrnd([2 3],[1 0;0 2],800);PntSet2=mvnrnd([6 7],[1 0;0 2],800);PntSet3=mvnrnd([8 9],[1 0;0 1],800);PntSet=[PntSet1;PntSet2;PntSet3];colorList=[0.9400    0.9700    0.9600    0.8900    0.9300    0.9200    0.8200    0.9100    0.8800    0.6900    0.8500    0.7700    0.5900    0.7800    0.6900    0.5500    0.7500    0.6500    0.4500    0.6500    0.5600    0.4000    0.5800    0.4900    0.3500    0.5100    0.4200    0.2500    0.3600    0.3100    0.1300    0.1700    0.1400];CData=density2C(PntSet(:,1),PntSet(:,2),-2:0.1:15,-2:0.1:15,colorList);set(gcf,'Color',[1 1 1]);% 主分布图ax1=axes('Parent',gcf);hold(ax1,'on')scatter(ax1,PntSet(:,1),PntSet(:,2),'filled','CData',CData);ax1.Position=[0.1,0.1,0.6,0.6];% X轴直方图ax2=axes('Parent',gcf);hold(ax2,'on')histogram(ax2,PntSet(:,1),'FaceColor',[0.78 0.88 0.82],...    'EdgeColor','none','FaceAlpha',0.7)ax2.Position=[0.1,0.75,0.6,0.15];ax2.YColor='none';ax2.XTickLabel='';ax2.TickDir='out';ax2.XLim=ax1.XLim;% Y轴直方图ax3=axes('Parent',gcf);hold(ax3,'on')histogram(ax3,PntSet(:,2),'FaceColor',[0.78 0.88 0.82],...    'EdgeColor','none','FaceAlpha',0.7,'Orientation','horizontal')ax3.Position=[0.75,0.1,0.15,0.6];ax3.XColor='none';ax3.YTickLabel='';ax3.TickDir='out';ax3.YLim=ax1.YLim;

Matlab怎么绘制散点密度图

3.4带直方图的等高线图

PntSet1=mvnrnd([2 3],[1 0;0 2],800);PntSet2=mvnrnd([6 7],[1 0;0 2],800);PntSet3=mvnrnd([8 9],[1 0;0 1],800);PntSet=[PntSet1;PntSet2;PntSet3];colorList=[0.9300    0.9500    0.9700    0.7900    0.8400    0.9100    0.6500    0.7300    0.8500    0.5100    0.6200    0.7900    0.3700    0.5100    0.7300    0.2700    0.4100    0.6300    0.2100    0.3200    0.4900    0.1500    0.2200    0.3500    0.0900    0.1300    0.2100    0.0300    0.0400    0.0700];[~,~,XMesh,YMesh,ZMesh,colorList]=density2C(PntSet(:,1),PntSet(:,2),-2:0.1:13,-2:0.1:13,colorList);set(gcf,'Color',[1 1 1]);% 主分布图ax1=axes('Parent',gcf);hold(ax1,'on')colormap(colorList)contourf(XMesh,YMesh,ZMesh,10,'EdgeColor','none')ax1.Position=[0.1,0.1,0.6,0.6];ax1.TickDir='out';% X轴直方图ax2=axes('Parent',gcf);hold(ax2,'on')[f,xi]=ksdensity(PntSet(:,1));fill([xi,xi(1)],[f,0],[0.34 0.47 0.71],'FaceAlpha',...    0.3,'EdgeColor',[0.34 0.47 0.71],'LineWidth',1.2)ax2.Position=[0.1,0.75,0.6,0.15];ax2.YColor='none';ax2.XTickLabel='';ax2.TickDir='out';ax2.XLim=ax1.XLim;% Y轴直方图ax3=axes('Parent',gcf);hold(ax3,'on')[f,yi]=ksdensity(PntSet(:,2));fill([f,0],[yi,yi(1)],[0.34 0.47 0.71],'FaceAlpha',...    0.3,'EdgeColor',[0.34 0.47 0.71],'LineWidth',1.2)ax3.Position=[0.75,0.1,0.15,0.6];ax3.XColor='none';ax3.YTickLabel='';ax3.TickDir='out';ax3.YLim=ax1.YLim;

Matlab怎么绘制散点密度图

4使用方式扩展&ndash;与ggplot修饰器联动

ggplot风格修饰器:(点击图片跳转链接)

Matlab怎么绘制散点密度图

示例1

PntSet1=mvnrnd([2 3],[1 0;0 2],800);PntSet2=mvnrnd([6 7],[1 0;0 2],800);PntSet3=mvnrnd([8 9],[1 0;0 1],800);PntSet=[PntSet1;PntSet2;PntSet3];ax=gca;ax.XLim=[-1 13];ax.YLim=[-1 13];ax=ggplotAxes2D(ax);CData=density2C(PntSet(:,1),PntSet(:,2),0:0.1:15,0:0.1:15);scatter(PntSet(:,1),PntSet(:,2),'filled','CData',CData);

是不是瞬间有那味了:

Matlab怎么绘制散点密度图

示例2

PntSet1=mvnrnd([2 3],[1 0;0 2],800);PntSet2=mvnrnd([6 7],[1 0;0 2],800);PntSet3=mvnrnd([8 9],[1 0;0 1],800);PntSet=[PntSet1;PntSet2;PntSet3];ax=gca;ax.XLim=[-3 13];ax.YLim=[-3 13];ax=ggplotAxes2D(ax);[~,~,XMesh,YMesh,ZMesh,colorList]=density2C(PntSet(:,1),PntSet(:,2),-2:0.1:12,-2:0.1:12);colormap(colorList)contourf(XMesh,YMesh,ZMesh,10)

Matlab怎么绘制散点密度图

以上是“Matlab怎么绘制散点密度图”这篇文章的所有内容,感谢各位的阅读!相信大家都有了一定的了解,希望分享的内容对大家有所帮助,如果还想学习更多知识,欢迎关注编程网行业资讯频道!

阅读原文内容投诉

免责声明:

① 本站未注明“稿件来源”的信息均来自网络整理。其文字、图片和音视频稿件的所属权归原作者所有。本站收集整理出于非商业性的教育和科研之目的,并不意味着本站赞同其观点或证实其内容的真实性。仅作为临时的测试数据,供内部测试之用。本站并未授权任何人以任何方式主动获取本站任何信息。

② 本站未注明“稿件来源”的临时测试数据将在测试完成后最终做删除处理。有问题或投稿请发送至: 邮箱/279061341@qq.com QQ/279061341

软考中级精品资料免费领

  • 历年真题答案解析
  • 备考技巧名师总结
  • 高频考点精准押题
  • 2024年上半年信息系统项目管理师第二批次真题及答案解析(完整版)

    难度     813人已做
    查看
  • 【考后总结】2024年5月26日信息系统项目管理师第2批次考情分析

    难度     354人已做
    查看
  • 【考后总结】2024年5月25日信息系统项目管理师第1批次考情分析

    难度     318人已做
    查看
  • 2024年上半年软考高项第一、二批次真题考点汇总(完整版)

    难度     435人已做
    查看
  • 2024年上半年系统架构设计师考试综合知识真题

    难度     224人已做
    查看

相关文章

发现更多好内容

猜你喜欢

AI推送时光机
位置:首页-资讯-后端开发
咦!没有更多了?去看看其它编程学习网 内容吧
首页课程
资料下载
问答资讯