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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
用matlab可以画出一组随机数的概率密度图. W$ l( a' A( T$ c: H% Z

+ S6 A! x# @8 N5 @+ tx=normrnd(1,0.5,1,1000);
/ o5 O2 a3 q: G4 l2 f- ]# sxmin=min(x);
. l$ x+ n  O8 M2 A" E$ x2 ]0 `xmax=max(x);9 `2 X- V6 ?0 ?7 @
xp=linspace(xmin,xmax,20);/ D7 J; H  M+ ^6 u6 ?/ F
f=ksdensity(x,xp);$ b1 e& |) P. M8 m$ V8 P, t# \
plot(xp,f,'*');
3 b: H8 R& u8 J2 e! K6 g  R* d如何画出xp所对应的概率分布值呢?
) t- B7 G! x" @( m, `好像cdfplot,只持一个输入参数cdfplot(x,xp),不支持cdfplot(x,xp)的形式;4 [% ^2 t( M0 ?6 }" `
谢谢!
6 G9 Y& w: a& p: Z
! E& H4 m4 W0 d( }3 s5 Y0 r8 j  m/ ?. b  m
  • 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); % 绘制频率直方图2 k, Y2 e, X% i4 g# J' c& i, l7 P
    1 ^) `. q! R# l( H

    9 ?3 o5 [* w6 }8 Y3 a9 u' Z9 x$ I
    + {9 O1 ~% z; s9 o
    经验分布函数(累积分布)图的画法:
    • ecdf(x); %或 cdfplot(x);7 G6 O2 x# A/ I( g9 |, ~8 s0 {

    7 u# @, A" `, R% j9 c# Z% r# ?) x; R
    . |0 U* _4 p- A$ g& N6 g6 v
    核密度函数图的画法:
    • x = sort(x);
    • ksdensity(x,x);

    • # z9 z2 B0 u. }2 @
    " C' y3 K; h# ?  r! P
    3 U+ D" g4 U' |; X) L

      Q# X" U' U  }8 d
    . V7 c; Z2 ~" k0 L' l& t

    该用户从未签到

    3#
    发表于 2020-5-6 16:19 | 只看该作者
    可以这样! T& r  z7 F' T- U  T; x0 k
    x = normrnd(0,1,1,1000);7 f6 G% {) x% `  }. U4 [
    [f, xc] = ecdf(x); % 调用ecdf函数计算xc处的经验分布函数值f
    3 @( A/ h- d3 ]' N- g! L7 Hxx = linspace(min(x),max(x),50);
    3 H, ~/ f& j! B& s) M- Z$ }y = [];. Q4 a0 h: w5 L7 G
    for i=1:length(xx)& [0 `! R! z- q. b" i( X, o4 k* ?7 u
        id = find(xc <= xx(i));" A4 w6 ]  W4 w! @
        if ~isempty(id)
    0 J" q$ q% p0 o3 @2 l        y(i) = f(id(end));" j% G: Q. {+ r; f! H
        else
    2 r. E) y, f! K' ^  M        y(i) = 0;3 L) S' ^% E, ~* w4 i: o
        end
    1 a9 S6 Z& j( A7 l+ u; Fend
    " w* x: d7 M' F- Dstairs(xc,f);2 @4 Y$ l& Z% ^# m4 w7 O% I, ?5 Y
    hold on
    9 [3 t$ ]7 k8 m$ i3 s8 o1 jplot(xx,y,'.r');

    该用户从未签到

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

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-7-22 13:59 , Processed in 0.125000 second(s), 26 queries , Gzip On.

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

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

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