|
|
可以这样
t9 L* a+ _/ o @7 v# W' cx = normrnd(0,1,1,1000);
& m; _4 c8 x* c0 f6 R) n) P0 @[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
8 l" i+ D8 s! i. [2 W# cxx = linspace(min(x),max(x),50);) t' C1 \! l, l# d1 i8 U& R" T
y = [];9 p2 @, Q3 q# Q) y, Y" M
for i=1:length(xx)/ ?9 h' K5 L4 D9 W( X/ e
id = find(xc <= xx(i));8 r, p' _/ q' q
if ~isempty(id), d5 M# y2 p# [
y(i) = f(id(end));& z/ K ^8 T' O' [
else- x6 K3 K1 {! ~ a
y(i) = 0;
$ P8 G4 W; d# T3 z; k: f" A$ ^# J6 l end& C# M% R0 n7 o. ?" R |3 B
end
& p% z4 ?: |% ~9 K, ~% I: Sstairs(xc,f);
, H- f- ^* f. O; j: Yhold on2 |) Q* U- F+ M: N: f) g
plot(xx,y,'.r'); |
|