EDA365电子论坛网

标题: Matlab如何画出一组随机数的概率分布图？？？ [打印本页]

---

作者: SpencerJ2    时间: 2021-6-28 14:16
标题: Matlab如何画出一组随机数的概率分布图？？？

( j) i, Y; ]* Y! E6 M0 Y9 p用matlab可以画出一组随机数的概率密度图
如
x=normrnd(1,0.5,1,1000);
xmin=min(x);
- b0 S  c- Y+ J& V; H5 Sxmax=max(x);
/ i9 _2 ^) v9 W) [1 ^! Qxp=linspace(xmin,xmax,20);
f=ksdensity(x,xp);
2 M# p! E' r# R2 x( Aplot(xp,f,'*');
3 Z. C0 A' }6 Z( x' \) U如何画出xp所对应的概率分布值呢？
5 s. Q+ G0 Y+ _1 c: }好像cdfplot，只持一个输入参数cdfplot(x,xp)，不支持cdfplot(x,xp)的形式;
- U- l1 w* e# G% Y; \# b谢谢！
+ ]1 W. m$ C% U0 ?7 j

---

作者: NingW    时间: 2021-6-28 15:17
竟然没读懂，难过了

---

作者: piday123    时间: 2021-6-28 15:20
频率直方图（注意不是频数直方图）的画法：

• x=normrnd(1,0.5,1,1000);
• [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
• ecdfhist(f, xc); % 绘制频率直方图
  2 h# U& _5 c% u6 L; J9 O) h

: U# b. O( F# G: S" p5 u/ b/ E

经验分布函数（累积分布）图的画法：

• ecdf(x); %或 cdfplot(x);
  * F- j7 b! R" D5 ?

5 k! b6 [$ Z/ z; ]( o7 A6 R% B( R核密度函数图的画法：

• x = sort(x);
• ksdensity(x,x);
• 
  8 y* S8 S! ~6 ?" L. v& X+ h( W

---

作者: SpencerJ2    时间: 2021-6-28 15:24

piday123 发表于 2021-6-28 15:20
频率直方图（注意不是频数直方图）的画法：

x=normrnd(1,0.5,1,1000);

[f, xc] = ecdf(x); % 调用e ...

( j4 x9 p* J4 s- a! |cdfplot或者ecdf+plot这两个命令无法画出像一楼贴出的那种离散效果，也就是说
; N( A) I, `. H" K5 L我们无法控制x轴的取点的间距
比如
, \  r7 ^; ^) k8 U% [[f,xi]=ecdf(x);
plot(xi,f,'.')

5 W# t% @" g: w% k% G# c3 j
请问能否采用什么方法，控制x轴的取点间距，得的所需要的离散效果图。
. S+ y/ s$ Y: T- \1 _如


( d, }  f2 S8 v* ^
( R" u1 s# h2 L
& F$ H/ n) W1 x

---

作者: piday123    时间: 2021-6-28 15:24

SpencerJ2 发表于 2021-6-28 15:24
cdfplot或者ecdf+plot这两个命令无法画出像一楼贴出的那种离散效果，也就是说
5 `% _# p- P* L/ D7 ]我们无法控制x轴的取点的 ...

" o  t6 A) a$ K+ N" _# A可以这样

• x = normrnd(0,1,1,1000);
• [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
• xx = linspace(min(x),max(x),50);
• y = [];
• for i=1:length(xx)
•     id = find(xc <= xx(i));
•     if ~isempty(id)
•         y(i) = f(id(end));
•     else
•         y(i) = 0;
•     end
• end
• stairs(xc,f);
• hold on
• plot(xx,y,'.r');
  - ]- Z; L+ s4 p" g( `

+ B( y* v1 s( t/ E* N

---

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)

Powered by Discuz! X3.2