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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图! w7 Y1 [" U, k0 u" z

0 U0 Q: B2 x- W( Bx=normrnd(1,0.5,1,1000);0 l6 @9 I1 M% _( n" d
xmin=min(x);
3 f* t; w: a& c: fxmax=max(x);( t" [$ P+ y4 w: J
xp=linspace(xmin,xmax,20);, |9 {  E$ s/ w9 u! x% k0 ~3 I
f=ksdensity(x,xp);' Z1 i4 b( ?# ?+ E; A( T+ [
plot(xp,f,'*');
& R" [. R. ~: E6 j+ R如何画出xp所对应的概率分布值呢?
, Z% x6 q( C1 o/ ~好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;
* c, J- n+ j& x5 O6 l! P谢谢!$ i2 r: p; }1 h& i) ~, L
0 t/ R# p; v  d- M& q
  • 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); % 绘制频率直方图
      3 G1 S9 h8 U2 w3 b2 R$ Q- G& ^

    # y' s# }5 E& Q
    2 Q4 i$ U7 [! p2 U/ J
    ; a+ D( j4 b+ N" h! W3 E, W
    经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);
      ' [& @: B9 d7 ]4 F8 u

    3 u( I- v! A5 K, H7 p% c& C
    / O4 h7 }0 U( Y4 X- i" m: o, e8 N' O
      n! `+ i7 b7 `
    核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);
    • : G+ T  P! P/ x

    8 v) b, [" R  b1 E0 W) R; [& z
    9 D) n) M! b5 ^4 v
    ) z+ k' ?9 K# K, j

    & N; B/ M' F% p# E4 E

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样& v6 V7 S4 Y+ W$ }' Q6 M; Y" E! x
    x = normrnd(0,1,1,1000);. `5 b% `5 S5 {6 r
    [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    5 ?, c9 ~' k" J% @( yxx = linspace(min(x),max(x),50);2 d/ w8 Q/ R4 w! {, A  f
    y = [];8 j/ G, S2 z8 T& _$ S4 B2 S- S& A
    for i=1:length(xx)# D4 z4 U" R3 q. {
        id = find(xc <= xx(i));
    / T% \4 P( E2 M: z% u    if ~isempty(id)
    5 j, u2 p# K9 P0 y3 \, u) z4 c& |        y(i) = f(id(end));
    - |" T8 ~8 r; i  z( Z. \    else
    7 ]+ b! b0 A4 C: ?, Q' P2 D        y(i) = 0;
    / K3 ]3 m7 @' ~$ ~- b: F    end/ ?: ?) i) h* L9 l
    end
    5 V- \6 P) F' o; S" T8 T* z9 X  w6 E4 tstairs(xc,f);
    3 m, h; E0 z& H! C7 fhold on( R5 G5 F, {1 g6 F
    plot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-6-22 20:58 , Processed in 0.078125 second(s), 26 queries , Gzip On.

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

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

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