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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图- U% G+ I. b1 t; A8 N# b8 H) }) |

3 j" q  u1 a, h$ m/ zx=normrnd(1,0.5,1,1000);
: C0 m  q5 ^, C, g$ Yxmin=min(x);9 b7 l" C$ J  n7 e
xmax=max(x);8 q8 m: Y: y+ s
xp=linspace(xmin,xmax,20);5 J; z5 W+ [4 a) a- d
f=ksdensity(x,xp);% T  [9 n, _0 U& K8 F0 @
plot(xp,f,'*');* w) G/ s* T$ L0 O6 x
如何画出xp所对应的概率分布值呢?
  {* b' E$ ]- N$ `9 e好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
: c) @$ B( I; N" ~谢谢!" ?( c; w. r8 o1 w) C

7 |/ A7 g5 ~0 t4 t7 ?) y
  • 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); % 绘制频率直方图
      9 W; v5 H4 b" c0 }# _) r
    ' Q: o; i/ ^# P

    . ^/ o$ Q7 u- [, _. W

    $ u/ A4 e/ q: N8 C2 V* K8 G经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);% Q. w3 s, q! S% Y, r0 Z* L6 w( R& f
    & n/ S" ^* w" m

    $ F/ }* @( o; Y2 v! g0 Y

    ! }8 R* a4 u3 v; ^: X/ \! t7 l核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);
    • 9 f9 [1 V( A' Z* [" D
    0 j3 ~& ^0 w) X# ?# v5 H3 V
    1 {, i( O! B2 v# o6 Y' g' S2 Y
    ( c1 t6 F) r; Z8 d$ b# }

    7 ]0 C$ E  C' V. U  E) Y

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样
    - u' j! |  V* k" r; Y% F: px = normrnd(0,1,1,1000);
    3 f+ }' ]4 m/ A2 J& X' A% i6 n[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
      X8 W" f- n; O- G* d1 ixx = linspace(min(x),max(x),50);
    6 s# f+ J1 o' J1 W7 q' Wy = [];. y' Y* d0 D8 S3 ]5 @5 Q8 S
    for i=1:length(xx)2 q0 s  Q6 G) ?; m" {1 R
        id = find(xc <= xx(i));" K0 Z8 b0 m2 G
        if ~isempty(id)
    8 ~# o, F$ \$ L  K/ b9 ^        y(i) = f(id(end));
    . @4 L3 s& ^# @# ]3 O; O    else
    . z6 m) V2 I0 P2 g1 `- L        y(i) = 0;
    4 T9 K5 u+ y# K$ o( D    end
    " h! d2 u% c! \$ m6 h1 B- Iend3 v5 J9 k5 C) S: k9 d, V. z/ l
    stairs(xc,f);
    ! f2 a0 }/ L1 u- `, a3 Fhold on7 W. k* N  i" @, H
    plot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-10-7 05:36 , Processed in 0.140625 second(s), 26 queries , Gzip On.

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

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

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