|
可以这样& v6 V7 S4 Y+ W$ }' Q6 M; Y" E! x
x = normrnd(0,1,1,1000);. `5 b% `5 S5 {6 r
[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
5 ?, c9 ~' k" J% @( yxx = linspace(min(x),max(x),50);2 d/ w8 Q/ R4 w! {, A f
y = [];8 j/ G, S2 z8 T& _$ S4 B2 S- S& A
for i=1:length(xx)# D4 z4 U" R3 q. {
id = find(xc <= xx(i));
/ T% \4 P( E2 M: z% u if ~isempty(id)
5 j, u2 p# K9 P0 y3 \, u) z4 c& | y(i) = f(id(end));
- |" T8 ~8 r; i z( Z. \ else
7 ]+ b! b0 A4 C: ?, Q' P2 D y(i) = 0;
/ K3 ]3 m7 @' ~$ ~- b: F end/ ?: ?) i) h* L9 l
end
5 V- \6 P) F' o; S" T8 T* z9 X w6 E4 tstairs(xc,f);
3 m, h; E0 z& H! C7 fhold on( R5 G5 F, {1 g6 F
plot(xx,y,'.r'); |
|