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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图
) Y4 u. C3 X# W9 g9 |: d# f
4 L* W( R. A. D; O7 r1 Z7 lx=normrnd(1,0.5,1,1000);
6 u& a$ V$ C$ Q8 H+ ixmin=min(x);
; X, s9 X  Y- w# H& \% ~  vxmax=max(x);' z/ o5 X' ?" g5 S8 `( n) x: M, L6 L
xp=linspace(xmin,xmax,20);
: u% L0 D2 @+ Ef=ksdensity(x,xp);! n# ]6 B, B  `$ @8 `  S
plot(xp,f,'*');/ [7 S/ F+ [7 W8 W
如何画出xp所对应的概率分布值呢?0 Y% B1 Y. h! u
好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;) |0 I. G: E, G' {% j
谢谢!- P% j/ |, [7 D0 S8 B: }$ l# w# P. g2 F
, E  U; |3 p' b9 P' I  z9 j
  • 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); % 绘制频率直方图! k1 n, A% }4 V! F/ h! y
    : E1 h. M+ }$ L, w0 h
    / Q( r$ H# w: w8 ?- G3 s

    ) A& D( h( u* Z: P8 m# v+ u8 @经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);1 t& F6 K) _' ^# s0 I7 W9 N

    ) d/ }, f, X; ~9 p$ B6 _' y( e$ b' {7 ^( N2 x2 f$ C' c; R) d& S
    " x, C1 ?; T4 G2 m& w
    核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);

    • 2 ~1 N' |! d. Y# U: i
    % ?7 C- D2 |# z/ a6 m/ _
    , f' r4 Z: C0 n& j" R& Z( [

    5 V8 }" L8 \* z" `3 y3 T1 X
    : j- h* _/ h$ _3 L% e

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样4 H, t  T& F2 }1 j! `. L
    x = normrnd(0,1,1,1000);
    ; c6 E1 g: s* v& x' o0 i[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f4 z# B* j# U6 B# V1 F7 B
    xx = linspace(min(x),max(x),50);
    2 K/ c9 t6 b9 V5 l9 wy = [];- c+ i2 a  W$ J0 i- u1 s
    for i=1:length(xx)& r! ]5 Q& z  A, [5 u$ a" V
        id = find(xc <= xx(i));& d" x( F' q. `# }8 ^
        if ~isempty(id)
      d- t( a( C( u) R* v        y(i) = f(id(end));
    ( }- e- l- }3 d" _6 b& B    else
    5 v- W2 j' X3 N2 e9 C4 H0 y        y(i) = 0;
    4 _( M- M: N; Q# m2 E    end% S# v" B6 ?/ }3 B% `& a6 H+ ?
    end
    ' l/ O! q( f3 b8 P2 Y2 gstairs(xc,f);: [5 e3 Q% n, s5 G: A0 f
    hold on$ Z# j" v( W# K+ |3 U  S- D
    plot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-8-23 22:11 , Processed in 0.140625 second(s), 26 queries , Gzip On.

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

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

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