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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图
0 _8 N6 ^! r1 ?: V: J1 H( u+ V$ @5 e% r4 W$ \! _
x=normrnd(1,0.5,1,1000);
- f7 y- |6 y7 i( C/ b6 gxmin=min(x);0 }( N& }) k1 _) |9 r7 {
xmax=max(x);+ D* G  `$ `8 p$ a* Y" l- |0 l
xp=linspace(xmin,xmax,20);
- I- i% T* t( |/ K$ K8 Of=ksdensity(x,xp);
3 ^$ |6 k: d( @plot(xp,f,'*');4 Q" Q: s9 T1 w5 o4 A
如何画出xp所对应的概率分布值呢?) z7 }* B$ z9 N
好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
. f# h! ]/ g. K4 Z# d谢谢!  r( A7 |$ }( ^! T

5 U% s1 r! M/ J! U+ C& ^3 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); % 绘制频率直方图8 S9 u& O" j( r6 r7 _# W; P: A2 S3 l
    4 n  p9 [; G7 u1 g$ y5 a
    / X- M2 I4 ?) H$ v0 N+ @

    . q5 Z& h8 `1 ~! E+ ~: ?经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);
        F* |* M; f8 G

    * N; S, ~7 Y0 j( k2 p0 k+ M8 I1 W; K! V

    ) g& A2 Y3 {5 F: o5 i核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);

    • # ~# l! ]# m8 M" _$ ~/ e/ S3 u

    : V; j( W4 V+ j+ `1 l
    8 t/ g2 S7 a! M! a, \4 Y6 D/ u# g
    4 K+ S) s' h4 Z% V, c/ G7 d

    5 Y5 O( c: C& Q3 n+ C

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样' e/ x' I7 J4 R
    x = normrnd(0,1,1,1000);6 _  E6 r2 \7 a/ j# L( j
    [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    $ w1 y( j3 D' ]) t, ?  ?' D* a7 Axx = linspace(min(x),max(x),50);, e, L1 g8 u( V
    y = [];
    4 x% `3 x3 b+ {$ Q& q' }) H8 X" F5 afor i=1:length(xx)5 W2 \$ L# I" w( s& @
        id = find(xc <= xx(i));
    + W: P8 C* o5 ]7 B; d8 E" K/ i    if ~isempty(id)# N0 f7 y4 e) h' J; B4 b
            y(i) = f(id(end));* ~7 K! l. P2 r1 V6 A* B' _, A) R' O
        else
    ' z* q" `/ Q- |; A- [        y(i) = 0;7 a8 q* f- v# Y6 g( c6 Z- a7 g
        end* p6 Y+ V5 V0 q" y' f3 [3 j/ R/ |; e
    end% s, o" X! ^, L; h. q, |. p
    stairs(xc,f);2 @0 ~7 E& k* o' |" b
    hold on
    5 [& c7 y! K/ e) l$ k' kplot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

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

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

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

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