clc;clear;
data = randn(1000, 1);%随机产生一个1000行1列的服从标准正态分布的随机数random normal distribution
[counts,centers] = hist(data, 7);%返回指定数据,特定间隔(如7)的各个分割区间的数量和中间值
% 画频率分布直方图
figure
bar(centers, counts / sum(counts))%画直方图,x是各个区间中间值,y是对应的概率(区间数/总数量)
% 分布参数拟合
[mu,sigma]=normfit(data); %计算数据data的正态分布参数特征值,均值μ和标准方差δ
% 画已知分布的概率密度曲线
x1 = -4:0.1:4;%获取密集的x间隔数据
y1 = pdf('Normal', x1, mu,sigma);%probability density functions概率密度函数,
hold on
plot(x1, y1)