概率论与数理统计学习:随机向量(三)——知识总结与C语言实现案例

news/2024/4/29 3:45:09/文章来源:https://blog.csdn.net/weixin_62917800/article/details/127215558

hello,大家好

这里是第八期概率论与数理统计的学习,我将用这篇博客去总结这期的知识点以及实现用C语言去做题的过程。
在这里插入图片描述

本期知识点:

  1. 条件分布
  • 条件分布的概念
  • 离散型随机变量的条件概率分布
  • 连续型随机变量的条件概率密度
  1. 随机变量的独立性

那么首先进入知识总结的环节
在这里插入图片描述

💦 条件分布

☁️ 条件分布的概念

请大家先回忆一下,我们最开始是不是也学过这个啥条件的东西?

对的,在前面的那叫条件概率,是对随机事件而言的,因为那时还没引入随机变量。(在这里呢大家可以看看前面那期有关条件概率的博客→\rightarrow条件概率)

注意了嗷,但这个知识点的标题是条件分布,其实两者概念上是差不多的,只是针对的知识点不同而已。

那么什么是条件分布呢?“分布”一词是与随机变量绑定在一起的,所以,也就是说条件分布一定是与随机变量有关的。

🌱 定义:设有两个随机变量XXXYYY,在给定了YYY取某个值或某些值的条件下,XXX的分布称为XXX的条件分布。同时也要注意❗️,这两个随机变量是有联系的,即它们之间可以构成随机向量的关系。

既然概念已经搞清楚了,又因为它还是涉及到随机变量,而我们所学习的随机变量又分为两种情况:离散型和连续型。所以下面又是对这两种情况分别讲解。
在这里插入图片描述

☁️ 离散型随机变量的条件概率分布

首先,我们需要借鉴一下下条件概率的知识点。对于两个事件AAABBB,我们知道在BBB已发生的条件下AAA发生的概率的求法为:P(A∣B)=P(AB)P(B)P(A|B)=\frac{P(AB)}{P(B)}P(AB)=P(B)P(AB),它表示的是A、BA、BAB同时发生的概率除以BBB发生的概率。

那么根据这个知识点,我们可以推出,对于给定YYY值,求XXX的分布列的求法为:X、YX、YXY两个值都给定的概率(对照上面的说法),除以YYY给定值的概率(就是不管XXX的值是多少),即P{X=xi∣Y=yi}=P{X=xi,Y=yi}P{Y=yi}P\{X=x_{i}|Y=y_{i}\}=\frac{P\{X=x_{i},Y=y_{i}\}}{P\{Y=y_{i}\}}P{X=xiY=yi}=P{Y=yi}P{X=xi,Y=yi}
是不是这样子捏?

又由上一节可以知道啊,P{X=xi,Y=yi}=pijP\{X=x_{i},Y=y_{i}\}=p_{ij}P{X=xi,Y=yi}=pijP{Y=yi}=p⋅jP\{Y=y_{i}\}=p_{·j}P{Y=yi}=pj。于是👇

🌱 定义:对于给定的jjj,在Y=yjY=y_{j}Y=yj条件下随机变量的条件概率分布为:P{X=xi∣Y=yi}=pijp⋅ji=1,2...P\{X=x_{i}|Y=y_{i}\}=\frac{p_{ij}}{p_{·j}}~~~~~~~~i=1,2...P{X=xiY=yi}=pjpij        i=1,2...
同样,对于给定的iii,在X=xiX=x_{i}X=xi条件下随机变量的条件概率分布为:P{Y=yj∣X=xi}=pijpi⋅j=1,2...P\{Y=y_{j}|X=x_{i}\}=\frac{p_{ij}}{p_{i·}}~~~~~~~~j=1,2...P{Y=yjX=xi}=pipij        j=1,2...

那么离散型随机变量的条件概率分布就搞定啦!
在这里插入图片描述

☁️ 连续型随机变量的条件概率密度

注意注意噢!上一个标题是条件概率分布,这个标题是条件概率密度了,这是为啥?

概率分布是离散型随机变量的性质,概率密度是连续型随机变量的性质(把离散型和连续型分清楚就好了)。

对于一个随机变量的给定值,求它的概率在连续型这是行不通的,由连续型随机变量的性质可知:对任意x,yx,yx,yP{X=x}=0,P{Y=y}=0P\{X=x\}=0,P\{Y=y\}=0P{X=x}=0,P{Y=y}=0。所以我们这里要把这个点扩散到一个范围,使它“连续”。也就是使这个点处于一个很小很小很小的一个范围内就行了。总之就是要使它有个范围!

🌱 定义:给定yyy,设对于任意固定的ϵ>0\epsilon >0ϵ>0P{y−ϵ<Y≤y+ϵ}>0P\{y-\epsilon <Y\leq y+\epsilon\}>0P{yϵ<Yy+ϵ}>0,且若对于任意实数xxx,极限lim⁡ϵ→0+P{X≤x∣y−ϵ<Y≤y+ϵ}=lim⁡ϵ→0+P{X≤x,y−ϵ<Y≤y+ϵ}P{y−ϵ<Y≤y+ϵ}\lim_{\epsilon\to0^+}P\{X\leq x|y-\epsilon<Y\leq y+\epsilon\}=\lim_{\epsilon\to0+}\frac{P\{X\leq x,y-\epsilon<Y\leq y+\epsilon\}}{P\{y-\epsilon<Y\leq y+\epsilon\}}ϵ0+limP{Xxyϵ<Yy+ϵ}=ϵ0+limP{yϵ<Yy+ϵ}P{Xx,yϵ<Yy+ϵ}
存在,则称此极限为在条件Y=yY=yY=yXXX的条件分布函数,记为P{X≤x∣Y=y}P\{X\leq x|Y=y\}P{XxY=y}FX∣Y(x∣y)F_{X|Y}(x|y)FXY(xy)

同样注意注意!
在这里插入图片描述
到这里,我们先清楚清除,这期的标题:连续型随机变量的条件概率密度,据我所知,概率密度好像是一个类型f(x)f(x)f(x)这样的“fff”开头的函数吧?

OK,从上面定义的结论可知,我们得到了条件分布函数,又由前面的知识可知啊,条件分布函数的导数就是概率密度函数(大概是这样,具体情况具体分析),也就是说有FX∣Y(x∣y)=∫−∞xfX∣Y(u∣y)duF_{X|Y}(x|y)=\int_{-\infty}^{x}f_{X|Y}(u|y)duFXY(xy)=xfXY(uy)du
那么,这个fX∣Y(x∣y)f_{X|Y}(x|y)fXY(xy)就是我们要找的概率密度函数啦!(上面的uuu是为了区分变量xxx和值xxx

好了,真正的定义它来了!!

🌱 定义:设二维连续型随机变量(X,Y)(X,Y)(X,Y)的概率密度函数为f(x,y)f(x,y)f(x,y)YYY的边缘概率密度函数为f(y)f(y)f(y),若f(x,y)f(x,y)f(x,y)在点(x,y)(x,y)(x,y)处连续,f(y)f(y)f(y)yyy处连续,则有:fX∣Y(x∣y)=f(x,y)f(y)f_{X|Y}(x|y)=\frac{f(x,y)}{f(y)}fXY(xy)=f(y)f(x,y)
同样的有:fY∣X(y∣x)=f(x,y)f(x)f_{Y|X}(y|x)=\frac{f(x,y)}{f(x)}fYX(yx)=f(x)f(x,y)
这里就不对它们进行证明咯,有兴趣的大家下去了解吧!
在这里插入图片描述

💦 随机变量的独立性

之前我们学过了两个随机事件的相互独立性,可表示为P(AB)=P(A)P(B)P(AB)=P(A)P(B)P(AB)=P(A)P(B)

那么又来对比,P{X≤x,Y≤y}=P{X≤x}P{Y≤y}P\{X\leq x,Y\leq y\}=P\{X\leq x\}P\{Y\leq y\}P{Xx,Yy}=P{Xx}P{Yy},而这个式子又可以写成为F(x,y)=F(x)F(y)F(x,y)=F(x)F(y)F(x,y)=F(x)F(y)
通常在题目中,有可能要根据F(x,y)F(x,y)F(x,y)F(x)F(x)F(x)F(y)F(y)F(y),这里可以介样子求:F(x)=F(x,+∞)F(x)=F(x,+\infty)F(x)=F(x,+),然后带入+∞+\infty+yyy中即可,同理可求得F(y)F(y)F(y)

(别问为啥要说一下这个,问就是今天做题想了半天)

由前面的知识可知,分布函数是随机变量(随机向量)的特性,也就是说适用于离散型和连续型,于是乎,下面又要分两种情况来介绍~~
在这里插入图片描述
不过这次就比较简单易懂了。

☁️ 离散型随即变量的独立性

依照随机事件的独立性直接给出公式:P{X=xi,Y=yi}=P{X=xi}P{Y=yi}P\{X=x_{i},Y=y_{i}\}=P\{X=x_{i}\}P\{Y=y_{i}\}P{X=xi,Y=yi}=P{X=xi}P{Y=yi}

☁️ 连续型随机变量的独立性

f(x,y)=f(x)f(y)f(x,y)=f(x)f(y)f(x,y)=f(x)f(y)

时刻牢记离散型于连续型的区别和特性!!!!!!
在这里插入图片描述

💦 C语言实现案例

先来一个简单的练练手

  1. 设二维随机向量(X,Y)(X,Y)(X,Y)的概率分布如下表所示,它的条件概率分布
X\Y025
10.150.250.35
30.050.180.02

题目分析:做题的话这个题直接按照公式套就行了,十分简单

#include <stdio.h>
#include <malloc.h>
int main()
{float **XY;				// use a two-dimensional array to restore the probability of (X,Y) ——用二维数组来存储(X,Y)的概率分布int row;				// 行int column;				// 列printf("Please input the values of the row and column:");	// 请输入行和列的值 		scanf("%d,%d",&row,&column);printf("\n");// allocate space for the two-dimensional array ——为二维数组分配空间XY = (float **) malloc (sizeof(float *) * row);// allocate space for the columns of each line ——为每一行分配列的空间for(int i = 0 ; i < row ; i++)XY[i] = (float *) malloc (sizeof(float) * column);printf("Please input the datas according to the line:\n");	// 请按行输入数据for(int i = 0 ; i < row ; i++){for(int j = 0 ; j < column ; j++){scanf("%f",*(XY + i) + j);}}// use two arraies X and Y to restore the serial numbers of X and Y ——用两个一维数组去存储X、Y的序号int *X = (int *) malloc (sizeof(int) * row);int *Y = (int *) malloc (sizeof(int) * column);printf("Please input the serial numbers of X:");	// 请输入X的序号scanf("%d,%d",X,X + 1);printf("Please input the serial numbers of Y:");	// 请输入Y的序号scanf("%d,%d,%d",Y,Y + 1,Y + 2);printf("\n");printf("The distribution of (X,Y) is:\n");    		// (X,Y)的概率分布为for(int i = 0 ; i < row ; i++){for(int j = 0 ; j < column ; j++)printf("%.2f ",XY[i][j]);printf("\n");}printf("\n");// start to calculate according to the distribution of conditional probability ——开始按照条件概率分布计算for(int i = 0 ; i < row ; i++){float sum = 0;for(int j = 0 ; j < column ; j++)sum += XY[i][j];for(int j = 0 ; j < column ; j++){float p = XY[i][j] / sum;printf("P{Y=%d|X=%d}=%.3f\n",Y[j],X[i],p);}}for(int j = 0 ; j < column ; j++){float sum = 0;for(int i = 0 ; i < row ; i++)sum += XY[i][j];for(int i = 0 ; i < row ; i++){float p = XY[i][j] / sum;printf("P{X=%d|Y=%d}=%.3f\n",X[i],Y[j],p);}}return 0;
}

在这里插入图片描述
代码分析:整体来说代码都挺简单的,就是在输出的时候,XXX的值不是0,1而是1,3,YYY的值不是0,1,2而是0,2,5。这里呢我们就需令设两个数组来存储它们的序号。

  1. 求上题中的XXXYYY是否相互独立?

题目分析:将数据带入公式判断即可。

#include <stdio.h>
#include <malloc.h>
#include <stdbool.h>
int main()
{int row;int column;printf("Please input the values of row and column:");scanf("%d,%d",&row,&column);float **XY;XY = (float **) malloc (sizeof(float *) * row);for(int i = 0 ; i < row ; i++)XY[i] = (float *) malloc (sizeof(float) * column);printf("Then please input the datas:\n");for(int i = 0 ; i < row ; i++)for(int j = 0 ; j < column ; j++)scanf("%f",&XY[i][j]);printf("The distrubution of (X,Y) is displayed as below:\n");for(int i = 0 ; i < row ; i++){for(int j = 0 ; j < column ; j++){printf("%.2f ",XY[i][j]);}printf("\n");}printf("\n");bool pos = true;int i = 0,j = 0;while(pos == true && i < row && j < column){float sum1 = 0,sum2 = 0;for(int k = 0 ; k < column ; k++)sum1 += XY[i][k];for(int k = 0 ; k < row ; k++)sum2 += XY[k][j];if(sum1 * sum2 != XY[i][j])pos = false;i++;j++;}if(pos == true)printf("X and Y are relatively independent."); // X和Y相对独立elseprintf("X and Y aren't relatetively independent."); // X和Y不相互独立return 0;
}

在这里插入图片描述
代码分析:代码其实跟上个题差不多的,就是算题的方法不太一样。

这期学习就到这里啦,下期再见
在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.luyixian.cn/news_show_398730.aspx

如若内容造成侵权/违法违规/事实不符,请联系dt猫网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ROS学习笔记三(TF的类)

1.数据类型 数据类型定义在tf/transform_datatypes.h.里 1.1 基本数据类型(Quaternion, Vector, Point, Pose, Transform) TypetfQuaterniontf::QuaternionVectortf::Vector3Pointtf::PointPosetf::PoseTransformtf::Transform 1.2 tf::Stamped tf::Stamped在上面的数据类型…

RocketMQ 5.0:无状态代理模式的探索与实践

本文作者&#xff1a;金吉祥&#xff0c; Apache RocketMQ PMC Member&#xff0c;阿里云智能高级技术专家 背景 首先&#xff0c;让我们来看下是遇到了哪些痛点问题&#xff0c;促使我们去探索一种无状态代理的RocketMQ新架构的&#xff1b; RocketMQ 拥有一套极简的架构&am…

安卓投屏 QtScrcpy

一、电脑安装adb 版本大于1.0.40以上 40不行 adb 1.0.41下载链接 链接&#xff1a;https://pan.baidu.com/s/1WIPI-p7a4ErTLFYHaTC2kw?pwdadbt 提取码&#xff1a;adbt 安装参考 https://blog.csdn.net/M7_xbc/article/details/122957311 二、打开无线调试并且配对 手机打…

驱动开发(10/10-林雪阵)

终端输入1--->LED1点亮 终端输入2--->LED2点亮 终端输入3--->LED3点亮 终端输入0--->LED熄灭 chdev.c (底层驱动代码&#xff09; #include <linux/init.h> #include <linux/module.h> #include <linux/fs.h> #include <linux/uaccess.h>…

【webrtc】rtp 扩展头的ID

ietf 文档地址扩展头的uri 是固定的,因此识别扩展是通过uri地址,而非ID。rtp 扩展头是有个ID的 文档的说法 : 是本地的id,15保留,只可以小于15 本地标识符值 15 保留用于将来的扩展和 不得用作标识符。如果 ID 值 15 是 遇到,它的长度字段应该被忽略,处理 整个扩展应在该…

redis之AOF和RDB持久化

写在前面 因为redis数据是基于内存的&#xff0c;为了避免服务器重启或者是宕机导致数据全部丢失&#xff0c;提供了数据持久化机制&#xff0c;即AOF(Append Only File)日志和RDB快照&#xff0c;接下来我们分别看下。 1&#xff1a;AOF 1.1&#xff1a;AOF日志的实现 首先…

(附源码)计算机毕业设计SSM政府项目管理平台

项目运行 环境配置&#xff1a; Jdk1.8 Tomcat7.0 Mysql HBuilderX&#xff08;Webstorm也行&#xff09; Eclispe&#xff08;IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持&#xff09;。 项目技术&#xff1a; SSM mybatis Maven Vue 等等组成&#xff0c;B/S模式 M…

回顾——PCB绘制

目录 一、原理图库原理图 二、PCB库 三、PCB 一、原理图库原理图 新建工程&#xff1a;文件——New——Project——Name&#xff08;这里复制一下&#xff09;自己选择保存路径 添加文件&#xff1a; 保存工程&#xff1a;(粘贴) 绘制原理图库、原理图 侧边栏消失&…

虚拟社会、区块链和元宇宙

1986年&#xff0c;早期的互联网供应商Quantum Link和娱乐公司Lucasfilm Games发布了第一款MMO游戏名为&#xff1a;《Habitat》基于虚拟角色的社交世界&#xff0c;玩家可以通过300波特的调制解调器(每分钟0.08美元)和用户的Commodore 64(595美元&#xff0c;按今天的价格约为…

Vue 动态换肤

效果如图&#xff1a; 源代码&#xff1a; <template><div :class"[son${temp}]" class"demo3">这是四点零八分的北京<br/>一片手的海浪翻动<br/>这是四点零八分的北京<br/>一声雄伟的汽笛长鸣<br/>北京车站高大的建…

Kafka监控EFAK(Kafka-eagle)部署与踩坑详细记录

环境 阿里云服务器centoskafka 2.4.1 &#xff08;1.1以上版本都能支持&#xff0c;低版本不太清楚了&#xff09;efak 3.0.1 当前时间官网下载的最新版&#xff08;原名叫kafka-eagle&#xff09;efak官网&#xff1a;http://www.kafka-eagle.org/jdk8 部署好的UI 部署流程 …

kafka系列——安装部署,相关命令,配置文件,底层存储结构,log和index文件

点击上方“罗晓胜”&#xff0c;马上关注&#xff0c;您的支持对我帮助很大 / 前言 / Kafka是最初由Linkedin公司开发&#xff0c;用scala语言编写的&#xff0c;是一个分布式、支持分区的&#xff08;partition&#xff09;、多副本的&#xff08;replica&#xff09;&…

时光机特效在哪里?推荐这三个实用软件给你

现如今有一款时光穿梭机的特效软件非常热门&#xff0c;它具有让照片中的人变年轻或者变老的功能&#xff0c;能让我们看到过去以及未来自己的样子&#xff0c;也能给我们枯燥的生活增添一点趣味感。那么大家是不是已经迫不及待地想知道时光穿梭机特效滤镜在哪了呢&#xff1f;…

Linux学习 -- docker的commit命令和本地镜像到云端

commit命令本地镜像到云端 一、commit命令 我们在使用镜像新建容器后&#xff0c;容器只是具有简易的Linux的功能&#xff0c;不具备一些常用的功能&#xff0c;如vim功能&#xff0c;因此我们需要给容器加上一些我们需要的功能。 使用如下命令&#xff1a; docker commit提交容…

不使用第三方库怎么实现【前端引导页】功能?

前言 随着应用功能越来越多&#xff0c;繁多而详细的功能使用和说明文档&#xff0c;已经不能满足时代追求 快速 的需求&#xff0c;而 引导页&#xff08;或分步引导&#xff09; 本质就是 化繁为简&#xff0c;将核心功能以更简单、简短、明了的文字指引用户去使用对应的功能…

Oracle数据库 | SQL语句解析

个人主页&#xff1a;&#x1f497;wei_shuo的个人主页 &#x1f3c0; Hello World &#xff01;&#x1f3c0; 文章目录一.Oracle启动及登录1.1 服务手启动即关闭1.2 SQL* PLUS命令二. 表的创建和维护2.1 创建表2.2 修改表2.3 重命名表2.4 截断表2.5 删除表三. 数据完整性与约…

【C++学习】类与对象(中)

&#x1f431;作者&#xff1a;一只大喵咪1201 &#x1f431;专栏&#xff1a;《C学习》 &#x1f525;格言&#xff1a;你只管努力&#xff0c;剩下的交给时间&#xff01; 类与对象&#xff08;中&#xff09;&#x1f638;构造函数&#x1f639;概念&#x1f639;特性&…

常见数据结构-散列表(上)理论

一&#xff0c;散列表理解 散列表的英文叫“Hash Table”&#xff0c;我们平时也叫它“哈希表”或者“Hash 表”&#xff0c;散列表用的是数组支持按照下标随机访问数据的特性&#xff0c;所以散列表其实就是数组的一种扩展&#xff0c;基于数组演化而来。 散列表是通过散列函…

bp神经网络performance怎么看,BP神经网络用什么软件

1、除了MATLAB能做BP神经网络&#xff0c;还有其他什么软件能做 除了MATLAB能做BP神经网络&#xff0c;还有其他什么软件能做 理论上编程语言都可以&#xff0c;比如VB&#xff0c;C语言&#xff0c;过程也都是建模、量化、运算及结果输出&#xff08;图、表&#xff09;&…

JavaScript设计模式(一):面向对象编程 - 继承

JavaScript设计模式 - 面向对象编程灵活的语言-JavaScript用对象收编变量对象的另一种形式(函数对象)真假对象(闭包和类)一个检测类函数的祖先写的都是看到的-面向对象编程创建一个类&#xff08;三种方式&#xff09;类的属性和方法通过闭包来实现类的静态变量定义&#xff0c…