- function [x,pdfx,cdfx]=pdf_cdf_gen(x0, M)
- up=max(x0);
- dn=min(x0);
- N=max(size(x0));
- delta=(up-dn)/M;
- pdfx=zeros(M+1,1);
- for n=1:N
- d0=floor((x0(n)-dn)/delta)+1;
- pdfx(d0)=pdfx(d0)+1;
- end
- x=dn:delta:up;
- pdfx=pdfx/(N*delta);
- for m=1:M+1
- cdfx(m)=sum(delta*pdfx(1:m));
- end