导入%数据
load x_chuli.txt
调用ecdf函数计算xc下的经验分布函数值f_ecdf
[f_ecdf，xc]=ecdf(x_Chuli )；
%创建新的图形窗口并创建频率直方图。 直方图对应于七个单元格之间
图形；
ecdfhist(f_ecdf，xc，97 )；
保持接通；
xlabel(x )； 将%标记为x轴
ylabel(f ) x ) )； %标记y轴
%调用ksdensity函数进行核密度估计
[f_ks1，xi1，U1]=ksdensity(x_Chuli )；
绘制%内核密度估计图，设置线宽为3的黑实线
plot(Xi1，f_ks1，' r '，' linewidth '，2 ) )。
ms=mean(x_Chuli )； %计算x的平均值
ss=STD(x_Chuli )； %计算x的标准偏差
%计算xi1上的正态分布密度函数值，正态分布平均值为ms，标准差为ss
f_norm=normpdf(Xi1，ms，ss )；
创建%正态分布密度函数图表，并设置线宽为3的红色单点划线
plot(Xi1，f_norm，' k-. '，' linewidth '，2 )。
%在绘图中添加尺寸框，尺寸框的位置在坐标系的左上角
legend (“频率直方图”、“核密度估计图”、“正态分布密度图”、“位置”和“通告”)。