EDA365电子论坛网
标题:
请问怎样用Matlab画出一组随机数的概率分布图
[打印本页]
作者:
smileqq
时间:
2020-5-6 13:36
标题:
请问怎样用Matlab画出一组随机数的概率分布图
用matlab可以画出一组随机数的概率密度图
* a6 g. U7 Y. V. L6 ~
如
8 ?# ]; O" _: a: W8 x
x=normrnd(1,0.5,1,1000);
- j Z$ {& {# C
xmin=min(x);
6 O6 Y* Z$ i: T( I* K$ k g* N- w0 ?2 u
xmax=max(x);
0 K( L' [! q! c* r) ]9 s
xp=linspace(xmin,xmax,20);
2 g, P. t V; M1 _" j
f=ksdensity(x,xp);
6 c, F. Q) b) I9 K) U4 }; F
plot(xp,f,'*');
# c& X" N7 u T4 m/ F: o. j1 H
如何画出xp所对应的概率分布值呢?
! f% h- K3 S/ g& U: r" k
好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
0 U/ v( Q( i7 u8 Y C2 u( {
谢谢!
9 Y$ T8 U5 f1 H6 g- j) z- T* n
14.jpg
(2.87 KB, 下载次数: 8)
下载附件
保存到相册
2020-5-6 13:36 上传
6 f$ E9 t, u& { @* ]# V4 x7 z
作者:
yin123
时间:
2020-5-6 14:05
频率直方图(注意不是频数直方图)的画法:
x=normrnd(1,0.5,1,1000);
[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
ecdfhist(f, xc); % 绘制频率直方图
" Q3 I# d9 J! }4 g
" y" W% _- s. y0 t, b; j1 b( h
9 Q8 w4 V' _& R* J) o
2 Q' }8 P G& Q, f0 i9 Y) @
经验分布函数(累积分布)图的画法:
ecdf(x); %或 cdfplot(x);
, m9 L( |' a; \) p a' Q: J& F: J
5 g9 x7 X- Z9 ^: k; M" p/ {
% A' j& r u' ~# x/ J# B# f2 s! w
+ i# w; O M) P o; j( e* y( v2 [
核密度函数图的画法:
x = sort(x);
ksdensity(x,x);
$ c8 l: I5 W2 G2 W& Q+ r: i' x
2 v6 j. T" p& U2 A- {2 N i6 v7 m
7 L' j- z2 B9 _
1 K0 S- u" F# b& g
% i8 t8 v9 C* M6 ~: L/ P6 R
作者:
SsaaM7
时间:
2020-5-6 16:19
可以这样
+ ~/ l. {$ X* i! s: y# v6 C: D
x = normrnd(0,1,1,1000);
/ U- _ v9 d8 h
[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
! B% e6 B* ~$ X& ^1 b8 L4 j
xx = linspace(min(x),max(x),50);
$ J0 g2 h& N, a& X
y = [];
! t+ ^3 m' }& Q( E% [# T, ~
for i=1:length(xx)
6 I$ Q3 m) v g: O- H+ V& A
id = find(xc <= xx(i));
4 T; L3 M$ K d: B) u9 |
if ~isempty(id)
6 a" [ c& ~* \ G0 Y
y(i) = f(id(end));
J) q4 E8 `# J2 A# \# f
else
1 r4 [# X4 t G. z) [+ K2 \6 b& C
y(i) = 0;
9 K2 Q. S5 P3 N. d h1 @
end
* N- l- F: W# @
end
0 j% v+ h! |, S) y; [; j, Z
stairs(xc,f);
4 ?( Z$ r8 w# U1 @* V3 `
hold on
: P# ?8 H3 U8 m/ u# T
plot(xx,y,'.r');
作者:
NingW
时间:
2020-5-6 16:28
牛人们好啊
作者:
ExxNEN
时间:
2020-5-6 16:48
进来看看
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2