matlab繪制繪制頻率分布直方圖和分布擬合

% 2021年7月29日23:51:21
% [e-mail] 
%參考：
% https://jingyan.baidu.com/article/636f38bb6bf4ddd6b9461043.html
% https://blog.csdn.net/qq_31816741/article/details/79360574
clc
clear
close all

%數據 一般至少50個
x = [0.312
0.333
0.334
0.333
0.315
0.345
0.295
0.295
0.256
0.366
0.344
0.344
0.314
0.323
0.275
0.356
0.345
0.314
0.284
0.333
0.333
0.293
0.333
0.295
0.313
0.333
0.313
0.264
0.333
0.334
0.295
0.242
0.344
0.345
0.323
0.322
0.271
0.366
0.324
0.291
0.322
0.413
0.426
0.426
0.419
0.284
0.313
0.315
0.344
0.314
0.345
0.296
0.323
0.322
0.283
0.315
0.304
0.238
0.314
0.21
0.313
0.333
0.261
0.235
0.333
0.417
0.249
0.425
0.224
0.433
0.293
0.426
0.325
0.425
0.305
0.418
0.304
0.425
0.324
0.324
];

% 畫頻率分布直方圖
[counts,centers] = hist(x,5);%柱的個數為5，可根據自己需要調整
x1=centers;%每個柱的中心x坐標
y1=counts / sum(counts);%每個柱的個數（頻數）與數據總個數的比值

% 分布參數擬合
[mu,sigma]=normfit(x);%用正態分布擬合出平均值和標准差

% 畫已知分布的概率密度曲線
x2 = centers(1)*0.5:((centers(end)-centers(1)))/1000:centers(end)*1.5;
% x2 =-centers(end)*1.5:((centers(end)-centers(1)))/1000:centers(end)*1.5;
y2 = pdf('Normal', x2, mu,sigma);%probability density function，求在x2處的pdf值

figure
[hAxes,hBar,hLine]=plotyy(x1,y1,x2,y2,'bar','plot')
set(hLine,'color',[1,0,0],'LineWidth',1,'Marker','o','MarkerSize',2,...
    'MarkerFace','y')
title('my bar and pdf curve')
xlabel('x')
ylabel(hAxes(1),'y1')
ylabel(hAxes(2),'y2')