进来看看

牛人们好啊

可以这样
x = normrnd(0,1,1,1000);
! U/ _  p* x1 @; e9 s[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
xx = linspace(min(x),max(x),50);
y = [];
for i=1:length(xx)
8 {0 R$ e! F* l0 G" i    id = find(xc <= xx(i));
    if ~isempty(id)
  S% D% D; [9 Q! i2 m1 b+ G+ e( x        y(i) = f(id(end));
+ K) e( ]9 h3 a5 R- x6 X+ m    else
0 q4 _$ d! i4 b! r$ Z        y(i) = 0;
    end
end
; c; y- w! r  R( L# H: l; `; U3 nstairs(xc,f);
hold on
1 z% Q- c2 y* v: e9 b3 Yplot(xx,y,'.r');

频率直方图（注意不是频数直方图）的画法：

• x=normrnd(1,0.5,1,1000);
• [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
• ecdfhist(f, xc); % 绘制频率直方图
  3 a. Q" [2 [) D/ V- w4 t! J

. W4 S* t! [6 y- R6 j, H

经验分布函数（累积分布）图的画法：

• ecdf(x); %或 cdfplot(x);
  2 v% T7 e8 w: i/ v% ?$ y! ^3 D& \

" ^$ H2 n+ Z/ V2 L6 ^! Y
; N4 u$ f  G. L0 A* }
% C6 `6 z1 n; j: c核密度函数图的画法：

• x = sort(x);
• ksdensity(x,x);
• 
  

& ?8 N. a, C! ~- e$ Z5 s
8 M4 q/ t% o- w