|
|
可以这样
9 M, Z3 k: \" Gx = normrnd(0,1,1,1000);
: _1 q0 o" C1 p/ D[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f' I1 e0 R3 J# u% K0 h2 I x# ~
xx = linspace(min(x),max(x),50);: c+ V% i8 I- }1 `# g
y = [];
% V: C3 i/ V' k& d N& bfor i=1:length(xx)" r; |0 H; s) W% l
id = find(xc <= xx(i));( R& S9 B4 C) p9 b" N
if ~isempty(id)
% @& Q; s" P( ~5 d+ J- J( Y y(i) = f(id(end));; z2 V4 K. R$ F5 }* k
else
, c' l- A2 {2 C, }0 i6 g8 D6 ] y(i) = 0;
. F# M# R* f$ I, F end. g1 ?. T, ^( U" R
end8 M5 U; b7 _9 z Y4 \' [
stairs(xc,f);
+ G. F1 s& s- L! y' ahold on, B& _7 z+ n3 @$ @
plot(xx,y,'.r'); |
|