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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图9 Q- |* H: u2 ?3 d9 k

7 r$ \8 s  ^$ c% J8 {x=normrnd(1,0.5,1,1000);
8 @6 ], ^) Q) ~4 zxmin=min(x);% F7 J* r0 S0 D! M; T4 @& D+ @$ y
xmax=max(x);
7 f4 E9 m4 s9 q" k: J/ P1 Uxp=linspace(xmin,xmax,20);3 o( o2 [4 ^6 h( n: E  r! g' I
f=ksdensity(x,xp);' f1 C1 |# R" t" M$ B' }1 `
plot(xp,f,'*');
$ F, x4 _1 M4 ]% s/ G+ d9 g如何画出xp所对应的概率分布值呢?
. w2 D- T$ l) y好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
1 _. h$ ], y; X) E谢谢!+ Z5 G+ o9 y( {' e5 Y
! d1 |: s8 ~& c; M& w5 Q' Z+ f
  • 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); % 绘制频率直方图
      . B( a% p, Y) I7 `6 n& j

    4 O1 I  r1 L: R4 o4 B1 t0 G' b0 W$ ?6 n" ^; O( l/ Z) A4 N, u
    * w) p: u3 K) w8 T: n, B
    经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);& R% B$ a9 T: R' ^# k" x! V% I( F

    & ~% a% v. D# ^
    # T. ^" p! W& u2 P! i! p* D1 m0 U
    & _' R4 ]2 U. H* W6 Z, {
    核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);
    • 1 Z( }' x3 V1 _. Z& e
    : s% ^7 t6 ]! k& P/ i" \; J
    # o$ N3 L$ {, e: ~; T; j( B
    $ \6 [. z$ s9 i& [+ M3 l2 g
    + j* ~8 Y" w8 N' B# c; y! u/ |

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样
    9 M, Z3 k: \" Gx = normrnd(0,1,1,1000);
    : _1 q0 o" C1 p/ D[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f' I1 e0 R3 J# u% K0 h2 I  x# ~
    xx = linspace(min(x),max(x),50);: c+ V% i8 I- }1 `# g
    y = [];
    % V: C3 i/ V' k& d  N& bfor i=1:length(xx)" r; |0 H; s) W% l
        id = find(xc <= xx(i));( R& S9 B4 C) p9 b" N
        if ~isempty(id)
    % @& Q; s" P( ~5 d+ J- J( Y        y(i) = f(id(end));; z2 V4 K. R$ F5 }* k
        else
    , c' l- A2 {2 C, }0 i6 g8 D6 ]        y(i) = 0;
    . F# M# R* f$ I, F    end. g1 ?. T, ^( U" R
    end8 M5 U; b7 _9 z  Y4 \' [
    stairs(xc,f);
    + G. F1 s& s- L! y' ahold on, B& _7 z+ n3 @$ @
    plot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-24 04:16 , Processed in 0.156250 second(s), 26 queries , Gzip On.

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

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

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