|
|
可以这样
7 ]$ E+ p/ u. b( [" wx = normrnd(0,1,1,1000);5 Z9 _6 W( l1 [' [- d
[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
' u. Z% N, d/ i |% xxx = linspace(min(x),max(x),50);; s* P$ E, ]$ v4 X: a2 R
y = [];
4 ^7 o% X* H) H3 b( m x! f) Vfor i=1:length(xx)+ G, h! z* r+ x) F
id = find(xc <= xx(i));* p; W! B" Q& y+ P' d' z
if ~isempty(id)
8 U% i) k' D! v# q) ^! E* _ y(i) = f(id(end));
: H, N7 t- P$ j% U: N+ q2 q else
7 y1 T- o9 F' L4 j3 r& h y(i) = 0;
9 q# j; e8 f) i+ w3 @1 C1 K end6 `+ P" g( |( N% m2 g w2 y# g
end0 |! o# W( a z6 ]* q% Q4 c5 o A
stairs(xc,f);
! w' r6 T0 h/ Vhold on* s% x+ Y9 d0 J
plot(xx,y,'.r'); |
|