首页/文章/ 详情

离散随机变量及其概率分布

3月前浏览1112


 

随机变量的本质

在大多数工程和科学应用中,随机事件的结果是非数值型。为方便进一步研究,将结果用实数轴  上的数值表示是必要的。建立这种映射或转换的函数称为随机变量。形式上讲,随机变量被定义为一个函数,函数的本质就是一种映射,它为随机事件的每个可能结果分配一个实数值。

随机变量通常分为两类:离散型和连续型。

离散型随机变量只能取不同的、通常是整数的值;例如,掷一枚骰 子,出现的结果只能是1到6之间的整数,因此它是离散型随机变量。某一河流某点每年洪水次数的取值也只能是整数,所以它也是离散型随机变量。连续型随机变量的取值则在连续的数值范围内;例如,机器零件的寿命理论上取值是  ,是一个连续型随机变量。由于这两类随机变量的概率计算方式不同,所以要对它们加以区分。

如上所述,一个离散的随机变量是一个函数,它将样本空间的所有结果映射到一个离散的实数集 合,称状态空间(  ),对于结果  ,以及  ,  是随机变量的一个状态。如图1所示,投掷一枚均匀硬币两次,考察正面向上的数量。由四个可能的结果组成样本空间:  。

▲图1

令  为“正面出现的次数”。这是一个可能取值为  的随机变量。作为函数,  给结果  赋值为2,给  和  赋值为1,给  赋值为0。即  ,  ,  。

令  为“反面出现的次数”。用  表示的话,  。换句话说,  和  是同一个随机变量:  对所有  成立。  因此,随机变量  为实验的每个可能结果  赋予一个数值  。

离散型随机变量的分布

给定一个随机变量,我们希望能用概率的语言来描述它的行为。例如,我们可能想回答关于该随机变量落入给定范围的概率问题:如果    是随机选择的市民的年收入,那么    超过十万元的概率是多少?

离散型随机变量的概率由概率质量函数(probability mass function, PMF)给出,它指定了离散型随机变量  等于某个状态值  的概率,记为:

对于所有的状态变量  ,  ,且  

在写    时,我们用    来表示一个事件(event),它由所有被    赋予数值    的结果    组成。这个事件也写作   ;形式上,   定义为   ,但写作    更简短直观。回到上面抛硬币的问题,如果    是正面朝上的次数,那么    由样本结果    和    组成,这是    赋予数值1的两个结果。由于    是样本空间的子集,它是一个事件。因此讨论    是有意义的,更一般地讲,谈论    也是有意义的。如果    不是事件,那么计算其概率就没有意义!写 "  " 是没有意义的;我们只能取一个事件的概率,而不能取随机变量的概率。

继续抛硬币的问题。在这个例子中,我们接下来将求出PMF。以下是我们定义的随机变量及其PMF:

定义  为正面朝上的次数。由于    发生时    等于0,   或    发生时等于1,   发生时等于2,   的PMF是函数   ,由下式给出

显然,  。且对所有其他    值,  。

定义  为反面朝上的次数。如上推理,或利用

  的PMF是

且对所有其他    值,  。

累积质量函数(cumulative mass function)用于表示随机变量  小于或等于  的概率。

根据定义,对于所有小于状态空间中最小值  的值,累积概率定义为零;对于所有大于状态空间中最大值  的值,累积概率定义为1。

投掷一枚均匀硬币三次,记随机变量  表示记录的正面(H)数量。X的可能取值为  。该实验的样本空间是三次抛硬币所有可能结果的集 合。由于每次抛掷的结果要么是正面(H),要么是反面(T),因此通过直接枚举,样本空间是由8个元素组成的集 合    。  的概率质量函数(PMF)和累积分布函数(CDF)如图2所示。

                             
                             
                             
                             

状态空间中的最大值为3,那么  .

▲图2



来源:数值分析与有限元编程
UM理论
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-06-16
最近编辑:3月前
译码当先
本科 慢慢来
获赞 14粉丝 32文章 400课程 0
点赞
收藏
作者推荐

用高斯积分法计算结构位移

高斯积分是常见的数值积分方法,是有限元计算中常用的数学工具。对于简化积分运算大有裨益。现在将其推广到结构力学,在计算结构位移时,用高斯积分来替代图乘法,同时用Python的sympy模块编程计算。一、核心思路在结构力学中,计算梁的挠度时,其莫尔积分公式为:其中:• --实际荷载引起的弯矩方程• -- 虚拟状态施加单位力引起的弯矩方程• --截面抗弯刚度• --梁的长度当 和 是复杂函数时,就需要用数值积分,而高斯积分是其中最精确、最高效的方法之一。二、高斯积分法高斯积分的基本思想是:在积分区间内,精心选择一些非等距的积分点(高斯点),用这些点上的函数值乘以相应的权重系数后求和,来近似积分值。对于一个在标准区间 上的积分,其公式为:其中:• --高斯积分点的个数• --第 个高斯点的坐标(在区间 内)• --第 个高斯点对应的权重• --被积函数重要特性:对于 个高斯点,该方法可以精确积分最高 次的多项式。根据被积函数 的多项式次数,选择足够的高斯点数 以保证精度。如果 和 是多项式,其乘积也是多项式,可以选择能精确积分的 值。常用 或 已能满足大部分工程精度。常用高斯点与权重表( 和 ): 高斯点坐标 权重 2 3 用高斯积分法计算位移的步骤:1. 坐标变换将积分区间从物理坐标 映射到自然坐标 。 于是位移积分变为:2. 定义被积函数令: 这里 和 是实际弯矩和单位弯矩方程用自然坐标 表示的函数。3. 计算积分近似值利用公式: 三、实例应用一根简支梁,长 ,全跨受均布荷载 , 。用高斯积分法求跨中C点的竖向位移。1. 建立弯矩方程 (以左端A为原点)• 实际荷载弯矩 : • 虚拟状态的单位荷载弯矩 (在跨中C点加竖直向下的单位力1): 由于对称,我们可积分一半长度(0到2)再乘以2。2. 坐标变换 (对左半段 )• 区间长度 • 变换: , • 在左半段, 3. 定义被积函数 (对左半段) 4. 选择 高斯积分被积函数最高次为3次,故选择2个高斯点。python代码 import sympy as spx, EI, xi, L1, w = sp.symbols('x EI xi L1 w')Mx = 20*x - 5*x**2mx = 0.5*xf_x = (L1/(2*EI)) * Mx * mx * w # 被积函数f_xi = f_x.subs( x, L1/2 * (1 + xi) ) # 将 x 用xi替换# Gauss-Legendre 积分点(-0.57735 ≈ -1/√3, 1/0.57735 ≈ 1/√3)points = [-0.57735, 0.57735]weight = [1.0, 1.0]delta = 0for i in range(2): # 同时代入当前积分点的 xi 和权重 w delta += f_xi.subs({ xi: points[i], w: weight[i]} )r = delta.subs({EI: 1.6e4, L1: 2})2*r #对称性,左右半跨结果一样,故乘以2 来源:数值分析与有限元编程

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈