找回密码
 注册
关于网站域名变更的通知
查看: 445|回复: 4
打印 上一主题 下一主题

请问怎样用Matlab画出一组随机数的概率分布图

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-5-6 13:36 | 只看该作者 |只看大图 回帖奖励 |正序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
用matlab可以画出一组随机数的概率密度图
( |: ^3 P' s+ i  j: I! o1 u1 ^. I8 |6 o
+ b4 Y  N1 P9 d) p  x/ Z0 P6 V' ex=normrnd(1,0.5,1,1000);; I* ~' L; d! ^6 E7 s0 n: m: P
xmin=min(x);
6 o5 b1 Y7 M3 ?& K1 J. Y  sxmax=max(x);
2 e" B+ h" ?- F3 bxp=linspace(xmin,xmax,20);
0 T( ?" [( p* \9 j2 V% I6 Q7 ef=ksdensity(x,xp);
) @$ D' E# j; v5 i" Q- s& r4 w& c1 yplot(xp,f,'*');1 G, k) S4 A1 \* R4 |
如何画出xp所对应的概率分布值呢?
, Y4 T( i/ L7 \* [2 F2 D  `, O好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
8 \; h7 E  p, B: P0 G谢谢!3 `3 o9 J# p/ F6 n6 M- g4 w, I

& b$ l$ n+ }: P7 R3 M0 ?# y" f; r

该用户从未签到

4#
发表于 2020-5-6 16:28 | 只看该作者
牛人们好啊

该用户从未签到

3#
发表于 2020-5-6 16:19 | 只看该作者
可以这样
: Z1 ]. Q$ N  E* [x = normrnd(0,1,1,1000);' r% u6 K4 Z- Y" g2 |. y
[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f, ?: U9 Q/ u2 l) W
xx = linspace(min(x),max(x),50);
2 B) \3 g8 ?' N( H* A& my = [];
3 z/ G- Q) z' c" w$ rfor i=1:length(xx)6 z7 k# w& J/ U( Z. F( _
    id = find(xc <= xx(i));
- ]8 y" q/ `$ y    if ~isempty(id)
: K: m3 U2 c7 N. D4 |0 @        y(i) = f(id(end));; T) L9 Y4 }6 ]: |
    else- J% j* z) \/ H! D1 V
        y(i) = 0;5 H+ L! A" v' J; S
    end) b" h: X" I. b& N4 @* R: g' d
end# w; \9 B; P% M1 N
stairs(xc,f);
! f- ]1 Y" B5 x& khold on
+ q; y2 @$ n9 L2 O  I$ Y( J% Hplot(xx,y,'.r');
  • TA的每日心情

    2019-11-29 15:37
  • 签到天数: 1 天

    [LV.1]初来乍到

    2#
    发表于 2020-5-6 14:05 | 只看该作者
    频率直方图(注意不是频数直方图)的画法:
    • x=normrnd(1,0.5,1,1000);
    • [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    • ecdfhist(f, xc); % 绘制频率直方图( ?; k7 T! |. n/ f

    % ?* x0 L3 r( H$ s6 O/ a7 x- H& B% u5 i5 A% L6 R7 M

    % T9 o  {2 u5 K8 F5 l2 i/ `经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);! ?. B* T" [- x
    . Q' L* n/ O1 t$ ?6 K# K3 N) w" q
    ) i1 {9 Y7 V+ D& O9 h
    : ~8 b8 Z# g. g: P: G- m. o
    核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);

    • 4 x# I9 P0 i  w9 ^7 X# G8 a2 T

    1 O% v  K: r- b/ j- P9 ]1 n9 O0 N& R: O# E) d' k

    4 E4 m2 c) h! _1 b" r, V; s: C
    9 [: O) _, q0 `2 F% j4 i
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-11-24 08:00 , Processed in 0.156250 second(s), 27 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表