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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图
2 O' [6 M: B: W6 X0 u" p/ ~( r6 I1 x$ \4 _- Z
x=normrnd(1,0.5,1,1000);
' n) C/ Z/ I% w* p! C1 Jxmin=min(x);
4 E; d6 x6 ^: x4 N. f: T7 txmax=max(x);, {3 H* {$ O8 y5 k" m, Z
xp=linspace(xmin,xmax,20);
, [" b. _, `! Q/ q3 H' y6 O% ef=ksdensity(x,xp);5 t! r' B5 s) D% ]# ]: B
plot(xp,f,'*');
1 G; g8 @& b) }- v& c如何画出xp所对应的概率分布值呢?. V6 `' ?7 W# b8 g' V
好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
+ }5 F$ U9 l2 _' k3 T+ a谢谢!% j: C) g0 h, D

9 C7 N, _  E2 H/ T, _1 D
  • 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); % 绘制频率直方图
      # c9 e. `4 `: u# d1 {  {

    , S1 o( A# G! C- R6 l$ J& I) X- B0 b6 `$ D
    $ k3 A, e9 c% |5 p  N8 i: @
    经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);3 S$ B8 j4 N& m( Q
    2 f* U4 R1 e  K3 S; o
    * R. a3 y5 U9 Z6 |( a

    * @7 ^4 X' l- i6 Z" d" W$ b  E5 R核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);

    • 2 S  ^3 _: c8 N; x
    3 [4 m/ T1 z! B) v8 ~
    ) k, _4 w5 [  o3 W1 l: Y) w% }/ T

    + E$ w% @! e, F  I  S6 W6 n
    * s$ V/ F4 h; o" J+ y6 I

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样) B- t) }) Q0 s# x
    x = normrnd(0,1,1,1000);
    7 |- t8 m, k% G* Z* j[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    # e: U2 l) }( N7 d  vxx = linspace(min(x),max(x),50);
    4 _. m6 g8 D. b" r3 C" G- p( _y = [];
    ' [& t  I- t2 e  _) ~6 U: Zfor i=1:length(xx)
    ! Z4 L5 a4 V4 {7 y/ ?$ _& m) Q0 a  ?    id = find(xc <= xx(i));
    " W9 V/ Y# V* v6 e3 m    if ~isempty(id)
    + \& t& |, p4 U" P* z7 k5 t        y(i) = f(id(end));
    / a8 \2 W( H# {3 C" `! D8 e    else8 {* F* J4 }1 ^, U: t. X- Q
            y(i) = 0;7 h- g6 H  N. c  Y2 i2 m0 Z0 h
        end+ ]. i7 W1 g. y1 i0 v$ o$ s1 ?
    end
    : Z4 P' ^6 P$ }2 Y+ R# Wstairs(xc,f);! u- z6 N  g4 ~7 o1 R3 k6 o
    hold on$ ?3 U4 z4 F' ^" C9 S5 b0 O, r
    plot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-4 16:45 , Processed in 0.187500 second(s), 26 queries , Gzip On.

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

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

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