Matlab绘制直方图、概率密度函数、累积分布函数

视频学习：https://www.bilibili.com/video/BV1HK411T76d?from=search&seid=16338053070486680597&spm_id_from=333.337.0.0

总结如下：

%%绘制直方图+概率密度函数+累积分布函数
%矩阵A是random生成0-1的随机数，或读入自己的数据
A = rand(1,500)*50;

%求A的最大、最小、均值
maxA = max(A);
minA = min(A);
meanA = mean(A);
%求极差、中位数、偏度、峰度
RangeA = range(A);
medianA = median(A);
skewnessA = skewness(A);
kurtosisA = kurtosis(A);

%生成一个画布
pic = figure; 
%绘制直方图(A中的元素有序划分入 x 轴上介于A的最小值和最大值间的 50 个等间距 bin中),centers表示每个bin中心在x轴
[counts,centers] = hist(A,50);
%绘制条形图，counts / sum(counts)计算频率
pic1 = bar(centers,counts / sum(counts),1);

set(pic1,'FaceColor','y')%设置矩形框颜色，‘b’为蓝色
set(pic1,'BarWidth',1);%设置矩形框宽度
%gcf返回当前axis对象的句柄值（一幅图像的每一个组成部分是一个对象）
%gca返回当前Figure对象的句柄值
set(gcf,'position',[220,200,900,600]);%设置画布的位置（220，200）和大小（宽900高600）
set(gcf,'position',[0.1,0.1,0.8,0.8]);%调整图形的位置和比例

axis([0 50 0 0.035]);%设置坐标轴范围[xmin xmax ymin ymax]

%计算概率密度函数,返回在xi点的概率密度f
[f,xi] = ksdensity(A);
hold on   %保留当前figure,实现图层叠加
pic2 = plot(xi,f,'LineWidth',1,'LineStyle','-','color','K');%绘制概率密度曲线


%计算概率密度积分分布函数
one_A = unique(A); %unique返回不重复的元素，产生的结果按升序排序

n_elements = histc(A,one_A);%统计在给定区间内的值的个数，左闭右开

c_elements = cumsum(n_elements);%计算元素累加的函数  
 %AX是坐标轴的句柄，AX（1）是左边的纵轴，AX（2）是右边的纵轴，H1,H2保存的是图形句柄，即两条线
[AX,H1,H2] = plotyy(0,0,one_A,c_elements);
%设置坐标轴的属性
set(AX(1),'ytick',0:0.005:0.035);%设置左边纵轴的刻度
set(AX(2),'ytick',0:50:500);%设置右边纵轴的刻度
set(AX(1),'ylim',[0,0.035]);%设置左边纵轴的取值范围
set(AX(2),'ylim',[0,500]);%设置右边纵轴的取值范围
%设置累积分布函数的线宽
set(H1,'LineWidth',2,'LineStyle','-','color',[0.9400 0.9400 0.9400]);
set(H2,'LineWidth',2)
%设置概率密度函数的线颜色
set(pic2,'LineWidth',2,'color','b');
%设置轴的颜色
set(AX(1),'ycolor','b');

set(AX(1),'ygrid','on');%打开左边纵轴网格
set(AX(2),'ygrid','on');%打开右边纵轴网格
%set(AX(1),'ygrid','off');%关闭左边纵轴网格

%设置title
title('test PDF')
%获取title的相关属性:get(gca,'title')
%设置字符、字体大小
set(get(gca,'title'),'String','Test PDF','FontSize',14);
set(get(gca,'title'),'Position',[25.001 0.0355 0]);%调整title位置

xlabel('X')%设置横轴标签
%set(get(gca,'xlabel'),'String','FontSize',14);%设置字体大小
set(get(AX(1),'ylabel'),'String','PDF','FontSize',14);%设置左纵轴标签
set(get(AX(2),'ylabel'),'String','CDF','FontSize',14);%设置右纵轴标签
set(get(AX(2),'ylabel'),'String','CDF','rotation',270,'Position',[77 250.0002-1.0000]);%设置旋转角度,标签位置

%添加文本 text(x,y,'')
textA = text(-10,0.025,{['Range:',num2str(RangeA)];['Max:',num2str(maxA)];['Mean:',num2str(meanA)];['Min:',num2str(minA)];['Kurtosis:',num2str(kurtosisA)];['Skewness:',num2str(skewnessA)]},'FontSize',14);
%根据图像显示，调整文本属性
set(textA,'FontSize',13,'Position',[10 0.025]);

%输出保存png
print(pic,'draw_hist.png','-dpng','-r200');%保存为png格式(-djpeg保存为JPEG格式)%‘-r200’是图像大小

运行结果如下：

 注：在编码编写时，可查看已设置的标签属性（必须先运行代码出现图像时，才能显示属性）：

例如：在命令行窗口输入：>>get(AX(2),'ylabel')，可查看对于右侧轴标签的属性