其他-n个互相独立的连续随机变量中第i小的数值期望
提出问题
有\(n\)个互相独立的\(0\)至\(1\)之间等概率生成的随机变量,求从小到大排序后第\(i\)个数的数值期望
一个简化的问题
我们先来求解一个简化的问题:最大值的数值期望是多少?
我们会发现,由于这些变量都是在\(0\)到\(1\)之间等概率生成的,所以一个变量小于等于\(x\)的概率为\(x\)(即\(P(x_0\leq x)=x\)),则这\(n\)个数中最大值为\(x\)的概率为\(x^{n-1}\)(其他\(n-1\)个变量都小于等于\(x\))
再考虑到有\(n\)个数都有可能成为最大值,所以最后答案还要再乘\(\binom n1\)(实际上这个组合数应该放在原式的概率函数\(p(x)\)里的,但为了表达方便,我们将这个组合数提到最外面最后进行计算,后面的运算也是如此)
由于期望的计算公式为
\[E(x)=\sum_{k=1}^{+\infty}x_kp_k\]
套到这题里就是
\[\int_0^1x\cdot x^{n-1}\cdot dx=\int_0^1x^n\cdot dx=\frac 1{n+1}\]
乘上组合数,得到这个简化问题的答案为\(\frac n{n+1}\)
扩展
我们现在求得了最大值(第\(n\)个数)的数值期望为\(\frac n{n+1}\),同理可以计算出最小数(第\(1\)个数)的数值期望为\(\frac 1{n+1}\),大胆猜想第\(i\)个数的数值期望为\(\frac i{n+1}\)
我们下面来证明这个式子
类比上面求最大值的解法,我们可以很容易地列出我们需要的式子
\[\int_0^1x\cdot x^{i-1}\cdot (1-x)^{n-i}\cdot dx=\int_0^1 x^i\cdot (1-x)^{n-i}\cdot dx\]
(第\(i\)个数为\(x\)的概率为前\(i-1\)个数都小于等于\(x\),后\(n-i\)个数都大于等于\(x\),则概率为\(x^{i-1}\cdot (1-x)^{n-i}\))
这个式子在最后还要乘一个\(n\cdot \binom {n-1}{i-1}\)(\(n\)个数都有可能成为第\(i\)个数,还要再选出小于等于\(x\)的\(i-1\)个数)
我们列出了式子,但这个式子并不像\(x^n\)这样好积分;为此,我们需要一些数学工具 妙妙工具
分部积分法
分部积分法由乘法法则推导而得
\[(uv)'=uv'+u'v\]
移项
\[uv'=(uv)'-u'v\]
两边同时积分
\[\int uv'\cdot dx=uv-\int u'v\cdot dx\]
积分
明确目标,我们要求
\[\int_0^1 x^i\cdot (1-x)^{n-i}\cdot dx\]
我们设\(u=(1-x)^{n-i},v=\frac 1{i+1}x^{i+1}\)
则\(u'=-(n-i)(1-x)^{n-i-1},v'=x^i\)
则我们要求的即为
\[\int_0^1uv'\cdot dx=(uv)\big|_0^1-\int_0^1 u'v\cdot dx\]
由于无论\(x\)取\(0\)还是\(1\),\(uv\)都为\(0\),则我们只需要考虑后面的式子即可
\[-\int_0^1 u'v\cdot dx\]
\[=-\int_0^1 -(n-i)(1-x)^{n-i-1}\frac 1{i+1}x^{i+1}\cdot dx\]
\[=\frac {n-i}{i+1}\cdot \int_0^1x^{i+1}(1-x)^{n-i-1}\cdot dx\]
数列
我们发现这个式子和和最初的式子的积分部分很像,可以对比一下:
\[\int_0^1 x^i\cdot (1-x)^{n-i}\cdot dx\]
\[\int_0^1x^{i+1}(1-x)^{n-i-1}\cdot dx\]
发现一个指数上升\(1\),一个下降\(1\),则我们设\(a_i=\int_0^1 x^i\cdot (1-x)^{n-i}\cdot dx\)
则我们可以得到一个有趣的递推式\(a_i=\frac {n-i}{i+1}\cdot a_{i+1}\),而边界条件即为我们一开始证明的式子\(a_n=\frac n{n+1}\)
从而我们可以得到\(a_i\)的通项公式\(a_i=\frac {i!(n-i)!}{(n+1)!}\)
所以我们前面那一长溜的积分式,可以化简为\(\frac {i!(n-i)!}{(n+1)!}\)!(忽略最后一个中文标点)
最后不要忘记我们之前提取出来的\(n\cdot \binom {n-1}{i-1}\)
最终解得的答案为\(\frac {i!(n-i)!}{(n+1)!}\cdot n\cdot \binom {n-1}{i-1}=\frac i{n+1}\)
猜想得证
结论
有\(n\)个互相独立的\(0\)至\(1\)之间等概率生成的随机变量,求从小到大排序后第\(i\)个数的数值期望为\(\frac i{n+1}\)
可以推广,若变量的生成范围为\([l,r]\),则第\(i\)小数的数值期望为\(l+\frac {i\cdot(r-l)}{n+1}\)
最近加了友链的一位同学貌似知道一个比较简单的做法?不过数学之美不恰好体现在各个方法看似不同,实则连通的吗
其他-n个互相独立的连续随机变量中第i小的数值期望的更多相关文章
- n个随机变量中第k小值的期望
Problem 有\(n\)个相互独立的取值在\([0,1]\)的随机变量,问其中第\(k\)小的变量期望是多少? Solution 之前pkuwc的神仙找我换友链,让我受宠若惊.. 我今天突然翻到他 ...
- SQL Server 2008 R2——查找最小nIndex,nIndex存在而nIndex+1不存在 求最小连续数组中的最大值
=================================版权声明================================= 版权声明:原创文章 谢绝转载 请通过右侧公告中的“联系邮 ...
- VS2013编译的exe独立运行在XP中方案
转载知乎 现在,我们深入探讨一下:<如何使用VS 2013发布一个可以在Windows XP中独立运行的可执行文件>. 这个问题是比较常见且容易造成初学者困惑的,作为曾经撞了无数次南墙的初 ...
- “独立”OpenVINO R2019_2 版本中的“super_resolution_demo”例子的,解决由于 R2019_1到R2019_2 升级造成的问题
OpenVINO提供了丰富的例子,为了方便研究和使用,我们需要将这些例子由原始的demo目录中分离出来,也就是“独立”运行,这里我们选择了较为简单的super_resolution_demo来说明问题 ...
- 独立的js文件中不能使用EL表达式取值
在独立的js文件中写了一个EL表达式取值,发现没有取到值,原因在于不能在独立的js文件中使用EL表达式,可以在jsp页面定义全局变量,然后在js文件中引用
- 关于ASIHTTPRequest连续请求,并发连续,间隔时间很小崩溃问题
在不停的刷新ASIHttpRequest的网络请求时,总是在刷新几次之后,整个app崩溃掉.我的app使用的ARC模式,以为可以自动释放到request的请求.经过摸索,还是需要在dealloc函数加 ...
- Javascript如何避免连续调用中取到不存在的属性而导致报TypeError错?
背景: 在最近的 NODEJS 项目中,涉及到数据库的查询,回调函数里返回了查询结果,我这样做处理然后返回给前端: return results.collect_coupon[0].count 但是这 ...
- 581. Shortest Unsorted Continuous Subarray连续数组中的递增异常情况
[抄题]: Given an integer array, you need to find one continuous subarray that if you only sort this su ...
- 关于 n 个 [0,1] 的随机变量第 k 小的期望值
今天做到一道题,感觉里面一个结论有点意思,就到网上扒了篇证明(bushi)下来了. 知乎回答习惯,先抛结论,再给证明(大雾 结论:对于 \(n\) 个取值范围为 \([0,1]\) 的随机变量 \(x ...
随机推荐
- 虚拟机复制问题:如何解决“无法手动启动VMware Tools安装”问题
摘录网址:https://jingyan.baidu.com/article/3c48dd348c4504e10ae35810.html 环境描述: VMware Workstation Ubuntu ...
- entityManager分页
十分操蛋. 需要两步. 第一步,查询一共需要多少条. 第二步 分页得到数据 Query query = this.entityManager.createNativeQuery(sb2.toStr ...
- springMVC的全局拦截器
先说说为什么要使用springMVC的全局拦截器,比如 当我们在访问接口的时候,我们一般都会先判断这个用户是否登陆,我们就要在每个接口的前面都要判断一下,想想是不是很蛋疼,那工作量... 这时候,我们 ...
- Kafka技术内幕 读书笔记之(二) 生产者——新生产者客户端
消息系统通常由生产者(producer ). 消费者( consumer )和消息代理( broker ) 三大部分组成,生产者会将消息写入消息代理,消费者会从消息代理中读取消息 . 对于消息代理而言 ...
- 5句话搞定ES5作用域
JavaScript的作用域一直以来是前端开发中比较难以理解的知识点,对于JavaScript的作用域主要记住几句话,走遍天下都不怕... 一.“JavaScript中无块级作用域” 在Java或C# ...
- 如何计算UDP/TCP检验和checksum
如何计算UDP/TCP检验和checksum 一.下面的图是一个UDP的检验和所需要用到的所有信息,包括三个部分:1.UDP伪首部2.UDP首部3.UDP的数据部分(切记不要遗漏该部分,否则就~吐血了 ...
- [leetcode-108,109] 将有序数组转换为二叉搜索树
109. 有序链表转换二叉搜索树 Given a singly linked list where elements are sorted in ascending order, convert it ...
- hdu 6385
题意是在一个矩形中任给N个点,求这N个点到矩形某边的最短距离和. 一开始想到直接贪心,求出每个点到矩形一边的最短距离,但题中说到线段间不能交叉,这里好像是比较麻烦,但题目中同时说了点与点之间的横纵坐标 ...
- 查看swap占用情况
查看swap被占用的情况 #!/bin/bash # Get current swap usage for all running processes # Erik Ljungstrom 27/05/ ...
- Kettle 和数据建模的几个学习资料
视频课程: 1. 初建军的 [慕课大巴分享]炼数成金——深入BI - Kettle 篇 基础书:1. Kettle 3.0 用户手册, 文件名为: ETL工具Kettle用户手册(上).pdf, ...