|
可以这样! T& r z7 F' T- U T; x0 k
x = normrnd(0,1,1,1000);7 f6 G% {) x% ` }. U4 [
[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
3 @( A/ h- d3 ]' N- g! L7 Hxx = linspace(min(x),max(x),50);
3 H, ~/ f& j! B& s) M- Z$ }y = [];. Q4 a0 h: w5 L7 G
for i=1:length(xx)& [0 `! R! z- q. b" i( X, o4 k* ?7 u
id = find(xc <= xx(i));" A4 w6 ] W4 w! @
if ~isempty(id)
0 J" q$ q% p0 o3 @2 l y(i) = f(id(end));" j% G: Q. {+ r; f! H
else
2 r. E) y, f! K' ^ M y(i) = 0;3 L) S' ^% E, ~* w4 i: o
end
1 a9 S6 Z& j( A7 l+ u; Fend
" w* x: d7 M' F- Dstairs(xc,f);2 @4 Y$ l& Z% ^# m4 w7 O% I, ?5 Y
hold on
9 [3 t$ ]7 k8 m$ i3 s8 o1 jplot(xx,y,'.r'); |
|