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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图
* o7 l5 y" B. m- X+ y2 w. y) K% M1 |4 {
x=normrnd(1,0.5,1,1000);* N0 {) R* P0 }
xmin=min(x);- Y: R# N. \: V  m
xmax=max(x);
$ F7 {$ `' i9 W6 B  V. E) sxp=linspace(xmin,xmax,20);
5 X3 ]& \/ R( Gf=ksdensity(x,xp);/ ]7 z: }4 K8 Q8 m# B
plot(xp,f,'*');
. i1 f0 e. Q' b' [如何画出xp所对应的概率分布值呢?4 }$ J3 M: b3 m3 g# W) @, X- S
好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
; ?6 Y* @' i. ]8 D" A# N8 ~谢谢!
- L; H( r# s' S4 G/ X5 N % o% y( Z$ p4 e
  • 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); % 绘制频率直方图
      ) E0 |8 E- r0 |  `% ~0 e& N
    $ |: F! S8 z- K/ f3 N0 r* x% N

    9 _  T6 O5 u$ Z1 t! t- X9 G& S8 s
    1 p! {9 `1 E3 G8 r1 b: h2 V( N
    经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);/ \+ O, }9 `$ G6 B
    ) t, F8 e# z+ y4 J

    ; a0 y& a5 r7 |; h0 ~" D
    6 I9 s, a0 x+ r5 c; l7 \( g
    核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);
    • ' A' q0 d& R% Z- c
    $ T0 `/ X% D1 S: U7 A4 o* M+ o
    ) ^: \, u8 o- z

    4 c0 r* q" j4 Q: p! m" N2 U: m7 o  S

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样
    7 ]$ E+ p/ u. b( [" wx = normrnd(0,1,1,1000);5 Z9 _6 W( l1 [' [- d
    [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    ' u. Z% N, d/ i  |% xxx = linspace(min(x),max(x),50);; s* P$ E, ]$ v4 X: a2 R
    y = [];
    4 ^7 o% X* H) H3 b( m  x! f) Vfor i=1:length(xx)+ G, h! z* r+ x) F
        id = find(xc <= xx(i));* p; W! B" Q& y+ P' d' z
        if ~isempty(id)
    8 U% i) k' D! v# q) ^! E* _        y(i) = f(id(end));
    : H, N7 t- P$ j% U: N+ q2 q    else
    7 y1 T- o9 F' L4 j3 r& h        y(i) = 0;
    9 q# j; e8 f) i+ w3 @1 C1 K    end6 `+ P" g( |( N% m2 g  w2 y# g
    end0 |! o# W( a  z6 ]* q% Q4 c5 o  A
    stairs(xc,f);
    ! w' r6 T0 h/ Vhold on* s% x+ Y9 d0 J
    plot(xx,y,'.r');

    该用户从未签到

    4#
    发表于 2020-5-6 16:28 | 只看该作者
    牛人们好啊
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-24 20:10 , Processed in 0.171875 second(s), 26 queries , Gzip On.

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

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

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