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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图# j' T; r. S/ Y" M8 I: o
5 ]" y4 v) W' J- m8 c
x=normrnd(1,0.5,1,1000);6 p/ H# W  S9 v& Q$ T& n
xmin=min(x);3 W7 ?0 E' q1 Z* y; r2 q9 r: D
xmax=max(x);  w& M2 {  [! c$ i3 {! g1 [0 h
xp=linspace(xmin,xmax,20);
' U: t) P- X. c5 o$ n8 M& m9 C" t' Gf=ksdensity(x,xp);
% k3 Y* T* o: o- y+ U7 a$ d* {+ P7 n4 Bplot(xp,f,'*');
  n% P+ g, E; V如何画出xp所对应的概率分布值呢?9 D( D) q+ g) e2 f
好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
0 }! b, E6 {9 }' Z3 i谢谢!
! T" I- q: q" i& y5 [; Z& {
: ]# Q+ F0 H, I2 U5 s; |" O+ \
  • 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); % 绘制频率直方图
      ) \  i. s9 v4 X* P+ _

    1 y4 Y3 \9 X# C# M0 g# V( _# T5 ^% j( l: j0 O+ C2 k9 I

    ( H5 m2 a4 Y& R3 r- F9 ]% l- y经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);( @" g& R- [, _- J. \( m( h
    ) K* `% y) Q' C9 Y0 ~/ {) }  t
    * ^$ y% a* p6 Z8 @

    # b& |0 [2 d7 w2 a/ j* W% Q核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);
    • : m& T5 ^+ {& y- A9 P
    8 Y: C! s! Z. @/ P( e" A* D) L
    * {: J- w# z0 E$ P

    - x2 B, P" \% X/ H
    & H: h  S4 ^- `  V+ }

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样
    ; k( a) S% P# E, ex = normrnd(0,1,1,1000);
      ~5 B! r7 X8 g5 y/ x$ p1 Y' W[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    6 F( J+ {" x" Y' Gxx = linspace(min(x),max(x),50);
    + e2 Y. V# r9 ~  H; l) B" q$ Vy = [];( T( I+ A4 b' r/ n1 @" b4 D
    for i=1:length(xx)7 H. w2 r& K' \
        id = find(xc <= xx(i));
    $ c* [$ O: Y( B' o7 Y    if ~isempty(id)
    $ L. d2 ]0 k! C        y(i) = f(id(end));! z4 ]9 J, A; f3 E2 m" x
        else
    3 Y* n  @0 b" U$ z' K' u        y(i) = 0;
    ; q9 {; \* i5 B    end
    - R: I: L! t0 i" G1 Q1 R9 {4 P& W7 Qend, L5 D' y4 y4 k# k; s, J
    stairs(xc,f);; D. J3 J% q4 v8 m5 Q
    hold on
    ! P; O* W, M+ O8 @5 f2 Y- qplot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-7-23 15:26 , Processed in 0.125000 second(s), 26 queries , Gzip On.

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

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

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