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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图
3 _* }/ u- Q5 g) e0 c0 R; x
) y5 P: U! n( w* E9 x6 nx=normrnd(1,0.5,1,1000);
6 i! F3 Y4 r. X2 [6 S. }2 Z" ~xmin=min(x);* p. `3 v+ I' y$ \
xmax=max(x);% `/ Y* G8 E. n
xp=linspace(xmin,xmax,20);; ^7 E; |8 k+ p
f=ksdensity(x,xp);; J; M: I2 G: {: v+ [
plot(xp,f,'*');) A* R5 W% z8 D3 u- D
如何画出xp所对应的概率分布值呢?
' T5 |. K7 z  F8 \$ H+ _" j" W好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
) y* d, m; G1 Y, M" p/ ^谢谢!
' m" l% K( k8 D4 J0 B8 [   G  I' e5 O9 ^1 P
  • 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); % 绘制频率直方图
      * l9 C! J" l" q$ w1 @' C3 F

    ( i4 |* F. }4 w
    8 \9 G7 n: Y5 \5 Z; }
    , `8 s6 p6 P1 I. P
    经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);
      / o* X( d) U. _* N) w' u
    ! I) S+ ]" R5 p  j3 Y
    & \' d- ^9 `3 i* G$ H# J0 h0 N
    ! P2 O& O$ ?7 [4 W
    核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);

    •   M2 V3 t0 k2 u& T
    - Q8 k0 s6 n0 w, D

    + c8 B! q# l2 S' q
    5 W7 V! {8 y9 [# S

    6 G, V( m( q4 Y0 Y5 X: f

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样
      t9 L* a+ _/ o  @7 v# W' cx = normrnd(0,1,1,1000);
    & m; _4 c8 x* c0 f6 R) n) P0 @[f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    8 l" i+ D8 s! i. [2 W# cxx = linspace(min(x),max(x),50);) t' C1 \! l, l# d1 i8 U& R" T
    y = [];9 p2 @, Q3 q# Q) y, Y" M
    for i=1:length(xx)/ ?9 h' K5 L4 D9 W( X/ e
        id = find(xc <= xx(i));8 r, p' _/ q' q
        if ~isempty(id), d5 M# y2 p# [
            y(i) = f(id(end));& z/ K  ^8 T' O' [
        else- x6 K3 K1 {! ~  a
            y(i) = 0;
    $ P8 G4 W; d# T3 z; k: f" A$ ^# J6 l    end& C# M% R0 n7 o. ?" R  |3 B
    end
    & p% z4 ?: |% ~9 K, ~% I: Sstairs(xc,f);
    , H- f- ^* f. O; j: Yhold on2 |) Q* U- F+ M: N: f) g
    plot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-24 03:34 , Processed in 0.156250 second(s), 27 queries , Gzip On.

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

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

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