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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图
% ^" I7 T6 Z. u6 T! s9 J8 E1 K: \/ k8 Y& \$ Q
x=normrnd(1,0.5,1,1000);
' M( q. C  N- \; D, e7 Uxmin=min(x);
# Q1 B- I6 D: A8 B3 xxmax=max(x);
' y' p; F' Q) }6 j( V$ J0 ]xp=linspace(xmin,xmax,20);1 K: G0 R# H' }0 L- @
f=ksdensity(x,xp);
, A  [4 \1 a. R  splot(xp,f,'*');  V, L3 h2 k6 P. J/ H
如何画出xp所对应的概率分布值呢?( ~" }+ b! |* d2 [
好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
$ h8 d+ Q; [9 i谢谢!
9 I/ ^* R2 y( E
6 H/ c+ r/ B, p$ L
  • 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); % 绘制频率直方图
      % {6 N5 `3 M, c" e& `: q2 O8 T
    3 s" J! W1 Y6 G4 V0 g8 M7 k/ p4 J1 }
    1 i) n9 {4 M$ j, @$ V/ D

    . G0 F" ~8 }; Z经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);6 i% |/ K4 Y1 |( @; K! k
    $ n( d! I5 v4 A) q# y* T" n0 e
    ; W' Y8 `6 P. i6 H, G

    % b1 @. v9 b! D; W核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);

    • : d8 p' I# ?+ m% I9 k

    / e& j" s- V" j* F1 ]5 q) I4 P! k7 s7 S
    8 }, u8 a# N3 \2 J- K" U
    & b9 v& y- I8 ~; t6 w

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样7 _$ x9 ~. `0 Y9 e
    x = normrnd(0,1,1,1000);+ S5 @( X  Q) e2 r  S0 M
    [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    . {* P7 C( `5 B& l/ [5 nxx = linspace(min(x),max(x),50);8 U! Q1 z/ |- J" q% M
    y = [];
    3 T! i" W+ R& g5 B1 w7 F3 Ofor i=1:length(xx)
    ; g! z/ p9 z6 g9 x# m$ f% I& T    id = find(xc <= xx(i));
    $ [) `: h; @* Y( B5 F( l( i    if ~isempty(id)
    ; ?& X9 Z$ l# H6 T7 O        y(i) = f(id(end));
      u! b& Y6 C8 K    else' T! Y2 j3 y$ J( M7 M
            y(i) = 0;
    " R* H) K$ h1 ~1 X7 P    end& |' c7 S! Q& r. w/ x- Z
    end
    / S, I- E3 p1 ]$ Bstairs(xc,f);
    5 |5 e9 L- k4 m0 X1 `hold on6 A1 `' j2 n6 v/ O5 A% X- \. v) {
    plot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

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

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

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

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