本系列题集将在本学期长期更新,收集老师课上布置的部分课堂习题、“黑脸”(难度较大的)作业题及部分小测题。
预备知识(概率论回顾与拓展)
设是独立同分布的随机变量,都服从正态分布,令.
(1)求(用标准正态分布的分布函数表示)
点击查看答案
将用表示,即求。又因为均为正态分布且独立同分布,故。
所以设,则,故
(2)验证随机变量和是独立的。
点击查看答案
证明:因为
所以,与不相关。又因为和均服从一维正态分布,所以由不相关可以推出独立。
甲、乙、丙三人在玩下面的游戏。每个人开始时都有且仅有元钱。
每秒钟就会响铃,铃响时每个有钱的玩家随机选择另外两名玩家中的一名并给予该玩家元
(例如:甲和乙可能会同时决定给丙,丙可能会决定把她的钱给甲,这样第一轮游戏结束时丙元,甲元,乙元。
在第二轮比赛中,乙没有钱给,但甲和丙可能会选择对方,那么在第二轮游戏结束时仍然是丙元,甲元,乙元)。
那么,次铃响后,
(1)每个玩家都有元钱的概率是多少?
点击查看答案
设第轮的状态:
:每个玩家都有元钱(种可能);
:一个玩家有元,一个玩家有元,剩下一个玩家没有钱(种可能)。
满足(没有其他可能)
若上一轮状态为,则下一轮有种可能,其中种可能结果为(甲给乙、乙给丙、丙给甲;甲给丙、丙给乙、乙给甲)
若上一轮状态为,则下一轮有种可能(因为没钱的人无法给钱),其中只有种可能结果为(元的玩家给元,元的玩家给元)
综上,,即每轮每个玩家有元钱的概率均为。
(2)玩家甲恰有元的概率?
点击查看答案
设第轮每个玩家都有元钱的状态为(种可能),玩家甲有元(种可能)的状态为,有元(非情况)的状态为(种可能),有元的状态为(种可能)。则
当然也可以用对称性理解:出现的概率是均等的(都有两种可能),且概率和为,故每个概率分别为。
设为取值在上的一个随机变量,证明:。
(此证明思想可用于证明马尔科夫/切比雪夫不等式)
点击查看答案
证明:
\begin{aligned}
\mathrm{Var}X&=\mathbb{E}[(X-\mathbb{E}X)^2]\\
&=\mathbb{E}[(X-\mathbb{E}X)^2\cdot 1_{X>\mathbb{E}X}]+\mathbb{E}[(X-\mathbb{E}X)^2\cdot 1_{X\leq \mathbb{E}X}]\\
&\leq \mathbb{E}[(b-\mathbb{E}X)^2\cdot 1_{X>\mathbb{E}X}]+\mathbb{E}[(a-\mathbb{E}X)^2\cdot 1_{X\leq \mathbb{E}X}]\\
&=(b-\mathbb{E}X)^2\cdot\underbrace{P(X>\mathbb{E}X)}_ y+(a-\mathbb{E}X)^2\cdot\underbrace{P(X\leq \mathbb{E}X)}_{1-y}\\
&\leq\frac{(b-a)^2}{4}
\end{aligned}
当且仅当时取等号。(最后一个不等号可以用二元函数求极值的方法求解)
这里我们直接证明一个一般的不等式:若随机变量的阶矩存在,则对任意,有
P(|X|>x)\leq\frac{\mathbb{E}(|X|^p)}{x^p}
马尔科夫不等式与切比雪夫不等式分别对应与的情形。
证明:令,则
\mathbb{E}Y=\mathbb{E}[Y\cdot 1_{Y>z}]+\mathbb{E}[Y\cdot 1_{Y\leq z}]
又因为,所以
P(Y>z)=P(X^p>x^p)\leq\frac{\mathbb{E}Y}{z}=\frac{|X|^p}{x^p}
又因为等价于,所以不等式成立。
设是独立同分布的随机变量序列,证明当时几乎处处收敛到当且仅当存在。(由此可构造依概率收敛无法推出几乎处处收敛的反例)
点击查看答案
\begin{aligned}
\sum_{k=1}^\infty P(A_k)&=\sum_{k=1}^\infty P(|X_k|>k\epsilon)\\
&=\sum_{k=1}^\infty kP(k\epsilon< |X_k|< (k+1)\epsilon)\\
&=\frac{1}{\epsilon}\sum_{k=1}^nk\epsilon P(k\epsilon< |X_k|< (k+1)\epsilon)\\
&=\frac{1}{\epsilon}\sum_{k=1}^\infty\int_{k\epsilon}^{(k+1)\epsilon}k\epsilon p_{X_k}(x)dx+\frac{1}{\epsilon}\sum_{k=1}^\infty\int_{-(k+1)\epsilon}^{-k\epsilon}k\epsilon p_{X_k}(x)dx\\
&\leq\frac{2}{\epsilon}\sum_{k=1}^\infty\int_{k\epsilon}^{(k+1)\epsilon}|x| p_{|X_k|}(x)dx\\
&=\frac{2}{\epsilon}\int_{\epsilon}^\infty |x|p_{|X_k|}(x)dx\\
&\leq\frac{2}{\epsilon}\mathbb{E}|X_k|< +\infty
\end{aligned}
最后一个小于是因为且存在(积分绝对收敛)。于是由Borel–Cantelli引理,可以无限次发生,
实际上,这里用了一个思想:
\begin{aligned}
\mathbb{E}|X|&=\int_{0}^\infty xdF(x)\\
&=-\int_{0}^\infty xd(1-F(x))\\
&=-x(1-F(x))|_{0}^{+\infty}+\int_{0}^\infty (1-F(x))dx\\
&=\int_{0}^\infty P(|X|>x)dx\approx \sum_{k=1}^{\infty}\epsilon P(|X|>k\epsilon)
\end{aligned}
无论服从何种分布,,所以
只要不存在(如柯西分布),即为依概率收敛无法推出几乎处处收敛的反例。
(秘书问题)某岗位有个应聘者,假设应聘者的面试顺序是完全随机的。每次只面试一人,只有通过和不通过两种结果,且面试结束后立即知晓,没有通过的应聘者,面试官就失去再次招聘他的机会,一旦某位应聘者面试通过,后面的面试就自然取消。假定面试官是公正的,面试通过的人在面试官面试过的应聘者中最优。为了使得面试官能以最大的概率录用到应聘者中的最佳人选,问面试最少要面试几个人?
点击查看答案
最优策略是:
先拒绝前个人,然后选择第一个比前个人都优秀的人。对的计算如下:设最优秀的人顺序为第个(概率为),则若则选到最优的概率为(被拒绝),若,则要保证最优秀的人被选中需要前个候选人中最优秀的人处于前个候选人中,概率为。所以,总概率为
P(r) = \sum_{k=1}^r\frac{1}{n}\cdot 0+\sum_{k=r+1}^n \frac{1}{n} \cdot \frac{r}{k-1}
当时,
\sum_{k=r+1}^n \frac{1}{n} \cdot \frac{r}{k-1}\approx\frac{r}{n}(\ln n-\ln r)=\frac{r}{n}\ln\frac{n}{r}
对关于的函数求导,得极大值为,当且仅当,即先拒绝前的候选人。
关于为什么这个策略是最优的,这里给出两个链接:算法趣谈(1)——秘书问题与37%法则,秘书问题——知乎,涉及到动态规划的知识,在此不赘述。实则我也没大理解ww
一个算法的伪代码如下:
Initialize T in {1,2,...,10} uniformly, S=0
While T < 6
For k in {1,2,...,T}
S=S+k
EndFor
Take T in {1,2,...,10} uniformly
EndWhile
Print S问:
(1)若第6行被注释掉,最后输出的S值平均是多少?
点击查看答案
在这一条件下,算法只对赋了一次值,那么While内语句就最多被执行一次。(而且只有时执行)
利用重期望公式,
(2)若第6行没有被注释掉,最后输出的S值平均是多少?
点击查看答案
将分为两个部分(设第一次):第一次For循环S的增量,与第一次For循环结束到整个While结束S的增量。
那么有
所以。
此题已初见随机过程雏形
设仪器寿命服从参数为的指数分布,假定仪器损坏后能及时更换新仪器,而且仪器之间的寿命相互没有影响。
那么从某个仪器装上后开始算,长为的时间内更换的仪器数量的平均数是多少?
点击查看答案
设为时间内仪器的更换数量(给定后属于随机变量),为仪器的寿命(服从指数分布,设其概率密度函数为),那么
又因为
其中表示时间内仪器的更换数量,且
所以
令得
再令,得
两边求导:
利用一阶线性非齐次微分方程的知识,解得,即.
事实上,这里表示的就是随机过程里面的泊松过程。
设随机变量满足均存在。
(1)试确定参数,使得取到最小值。(由此可得到一般的线性逼近)
点击查看答案
设
分别对和求偏导,得
解得
由此可得
此为最佳线性逼近(会在线性模型中出现)
(2)证明:设为关于的函数,则渠道最小值当且仅当。(即最佳逼近)
点击查看答案
证明:
其中
所以
去等号当且仅当。
补充:在随机向量典型题目第10题中,我们得到二元正态分布对应的条件概率
,这和上面第一问的线性逼近得到的结果几乎相同(除了顺序不一样)
实际上,这说明二元正态分布下最优逼近就是线性逼近。
(Sanov不等式):假设是取值在区间上的独立同分布随机变量,记.
证明:对任意,
其中.
点击查看答案
我们利用概率母函数(设为)证明:
下面对进行估计:
因为,,由函数凸性,
两边取期望得
代入上式:对求导,
解得
代回原不等式,即得
证明完成。
注:这里的也称为KL散度,而这个不等式说明样本均值的偏离程度会随样本增加呈指数级下降。
(Chernoff bound) 若随机变量X满足:对任意,.
证明:对任意,
点击查看答案
证明:因为
又因为上式对任意均成立,所以。(如果连续则可改为最小值)
(可用于估计一些非典型分布的概率)
简单随机模型
当时, 证明Bernoulli过程不是独立增量过程.
点击查看答案
独立增量过程的定义:相互独立
证明:
因为,而
所以当时,,,与不独立,得证。
已知,对任意,,
称为正弦波随机过程,其中为振幅,为角频率都是常数,是随机相位. 证明是严平稳过程。
点击查看答案
证明:
任意时刻,与任意非负整数,考查联合分布与:
因为
而,所以与同分布,
那么
与同分布,进而得到
与同分布,的严平稳性得证。
定义随机过程,其中, 其中,是两个不相关的随机变量,分布不同但.
(1)证明是平稳过程.
(2)若进一步假设,独立且都服从正态分布, 证明是严平稳过程.
(3)求的二维分布概率密度函数族.
点击查看答案
(1)证明:
因为
与无关;又因为(不相关),所以
与无关,所以为宽平稳过程。
(2)证明:
因为相互独立且均服从正态分布,所以服从二维标准正态分布
又
记,则到是单位正交变换,所以同样满足二维标准正态分布。
所以与同分布,进而任意与,
与同分布,为严平稳过程。
(3)即求的概率密度函数。由前两问可知,服从二维正态分布(具体概率密度函数略)
设是-简单随机游动, 其中,,
(1)求的分布(默认初始位置为)。
点击查看答案
因为是简单随机游动且右连续,所以
(2)求.
点击查看答案
其中与分别为从出发首次到的时间与首次到后,从出发到的时间.
然后
而,所以
即.(这也是为什么要求)
(3)证明.(即从任意位置出发,一定能经过有限次游动回到位置)。
点击查看答案
证明思路1:设表示从位置出发,最终经过有限次到达的概率。则
解得
其中为常数。又因为
所以
所以.
注:上面用到的差分方程在处理简单随机游动问题时非常有用,务必牢记。
证明思路2(使用条件期望):
因此
又因为,且,所以只能
当时,.
设是初值为的-简单随机游动,用表示在首次到达之前回到的次数。
已知, 对任意, 证明.
点击查看答案
证明:
设事件表示第次从出发,并在到达前回到。
那么(因为第一步到达后,,此为上题第3问结论),所以
对初值为的简单对称随机游动,求。
这一过程等价于从出发,经过步后第一次回到。
将情况分为种(另一种对称):第一步移动到,之后位置一直大于,到步回到。下面讨论这种情况的路径数。
首先,所有从开始经过步后回到的路径数为(即向上步+向下步)。
然后中间回到过的路径:根据反射原理,
因此路径数为(条上+条下)。
所以最终的路径数量为
而总路径数为,故概率为。
设是右连续随机游动,步长的分布为。记,令,求.
点击查看答案
设,则由条件期望的平滑性及全期望公式:
其中,(即从到经过个首回时),因此
因此。又因为
所以
设分别是强度为的泊松过程且相互独立. 证明是强度为的泊松过程.
点击查看答案
设 和 是两个独立的泊松过程,强度分别为 和 。
根据泊松过程的定义,任意,
且相互独立。因此
又因为与均为独立平稳增量过程,所以也是独立平稳增量过程。因此是强度为的泊松过程。
以表示的第次随机事件发生的时间.对任意,求:
(1);
点击查看答案
当时,;
下面讨论情况:
因为给定时,与均匀分布的最大次序统计量同分布,故的概率密度函数为,因此
综上,.
(2)与;
点击查看答案
可以发现这两个概率相等,这是由指数分布的无记忆性以及泊松过程的平稳增量性决定的。
(3)以表示的第次事件的间隔时间,对任意,求。
点击查看答案
当时,
当时,
又因为与相互独立,,所以
因此
综上,
离散时间随机游动假定了每一次游动都在整数时刻点发生,如果游动可以在任意时刻点发生,我们可以怎样建模?试着用泊松过程和简单-随机游动构造一个随机过程模型,并在时计算从位置出发首次到达位置的平均时间。
提示:连续时间随机游动
点击查看答案
模型建立:设-随机游动的游动时刻服从强度为的泊松过程,即在时间内发生次游动,相邻游动时间间隔服从指数分布,则可以定义连续时间随机游动:
下面求:
因为离散时间下(因为),即首次到达位置游动步数的期望为;又因为游动时间间隔的期望为,所以
某个零件在运行过程中收到撞击,用表示时间段内该零件受到的撞击次数。经验表明为参数的泊松过程。
每次撞击会给零件带来一定的磨损,假定每次磨损量为均服从参数为的指数分布且相互独立。若磨损总量达到或超过就要更换零件,问该零件的平均寿命。
点击查看答案
设零件的寿命为,则
则
又均服从,所以可以将看作强度为的泊松过程的到达时间。于是
因此
离散时间马氏链
(存储模型)设为独立同分布的非负整数值随机变量,令
其中为给定的两个正整数。
(1)证明为离散时间马氏链。
点击查看答案
(1)证明:显然的状态空间。对任意以及,
对任意, 由的定义可知为与的函数,如此递归代入后可知是随机变量的函数。这表明事件
由随机变量确定。再由的独立性,
所以为离散时间马氏链。
(2) 设。当时,写出的转移概率矩阵。
点击查看答案
(2)当时,,故矩阵第一行为
类似可得到剩下三行转移概率(请读者自行推导),最终得到转移概率矩阵为
(3) 在(2)的条件下求的分布。
点击查看答案
(3)初始分布为
则分布为
设是不可约周期马氏链,状态空间为。对任意,令
证明:
(1),互不相交而且。
点击查看答案
(1)先证明:显然,而对任意,因为不可约,所以存在正整数使.
那么取,则,故,所以;
再证明互不相交:若存在使,则存在,于是存在,使
又因为不可约,所以,存在使,那么
所以,于是。
但,故,矛盾。故互不相交。
(2)令,那么是以为一步转移概率矩阵的非周期马氏链。特别地,若,那么是状态空间为的非周期马氏链。
点击查看答案
(2)证明分三步:
- 证明是马氏链,且转移概率矩阵为
由的马氏性,
所以是马氏链,且一步转移概率矩阵为。
2. 证明是非周期马氏链
对任意状态,记中的周期为,而在中如果存在周期,则一定存在,使,所以,是非周期马氏链。
3. 证明若,那么是状态空间为的非周期马氏链
若 ,且 ,,则 。
现 的一步转移对应原链走 步:,余数 。
若初态 ,则
即转移矩阵 仅在 内部转移,不会跑到其余 。
又原链不可约,故对任意 ,存在 使得 ,因此 在 上不可约。
所以是状态空间为的非周期马氏链。
设是马氏链,。证明与独立且同分布。
点击查看答案
对任意正整数,
其中表示从出发首次返回的时间为的概率,恰等于,所以与独立同分布。
在马氏链中,如果常返,且,则常返。
证明:若,则存在正整数,使。又因为常返,所以,且(否则从出发无法到达,说明非常返)。
所以存在正整数,使。于是
所以常返。
在马氏链中,如果非常返,且,则非常返。
一个简单的反例:状态空间,转移概率矩阵,则非常返且,但是常返状态(吸收态)。
假设马氏链的状态空间为,转移概率矩阵为
令,求。
点击查看答案
注意到状态是吸收态,所以
此题关键在于注意到吸收态,由此简化首访条件。
记马氏链的状态空间为,转移概率矩阵为
记,求。
点击查看答案
进行前进一步分解:设分别表示出发状态与第一步到达状态,则
于是我们需要求。同理,
联立这两个方程,解得
设不可约正常返马氏链的转移概率矩阵为为平稳分布。
令,证明是状态空间为的不可约正常返马氏链,平稳分布为。
点击查看答案
- 是马氏链
对任意满足,
因此转移概率与之前的状态无关,故是马氏链。
2. 是不可约的
因为是不可约的,所以任意两个状态,有,即存在一条路径使。
因此对任意,,所以,是不可约的。
3. 存在平稳分布
定义,则
所以是概率分布。又因为对任意,
由于是的平稳分布,有
于是
故是的平稳分布。
4. 是正常返的
因为是不可约马氏链,且有平稳分布,所以一定是不可约正常返马氏链。
综上,命题得证。
证明:不可约马尔可夫链正常返的充要条件是存在非负数列及一个状态使得对任意,而且。
点击查看答案
证明:
- 必要性:若 不可约正常返,对任意 ,令 ,那么任意 ,
且
- 充分性:记,则。任取,对任意,迭代得,
这表明
两边除以并令得
因此为正常返状态,从而是正常返的。
实际上,这道题就是Foster-Lyapunov定理的一个特殊形式,其在稳定性理论中起着重要作用。(Lyapunov这个名字似乎在概率论中也出现过……)
(1)设不可约非常返或零常返马氏链的状态空间,证明:若数列满足
那么对任意,
点击查看答案
证明:任意,由 X 非常返或零常返可知,。注意到
先令,可得,再令,可知此时结论成立。
(2)设不可约链的状态空间,若存在正数及某个状态使得对任意,
而且,证明常返。
点击查看答案
证明:构造一个新马氏链使其转移概率为,
由于是不可约的,因此也是不可约的。对任意,,而且。注意到
若是非常返的,那么令,由上一问结论知,因此,这表明的状态是常返的,从而是常返的,矛盾。
因此一定是常返的,故方程组
的最小非负解恒为,由此可知是常返的。
当非负整数值马氏链的转移概率满足:
(1)
(2)
分别讨论的常返性(正常返/零常返/非常返)。
点击查看第1问答案
(1)由条件可知是不可约的。转移概率矩阵形式为
解法1:注意到这个马氏链(看作随机游动)是左连续的,设,即从状态出发最终到达状态的平均时间。那么可得到关系式
又因为对每个状态,其向下一步的结构相同,故所有相同(记为),代入得
然后再计算,即从出发最终回到的时间期望:
所以是正常返状态。又不可约,故是正常返的。
解法2:对任意,构造。那么对任意,
由上上题结论(Foster-Lyapunov定理)可知是正常返的。
点击查看第2问答案
(2)转移概率矩阵形式为
令,则由可知在内存在唯一解。对任意,令,那么
且对任意,
综上可知存在平稳分布,又由条件可知是不可约的,因此是正常返的。
那么这一问里面平稳分布是如何想到的呢?根据条件,平稳分布需满足
为了将变量与分离,考虑,则,因此只需要即可。【又因为,所以一定有解】
假设是一个不可约马氏链,状态空间为,转移概率满足条件,令为从出发到达之前先到达的概率,求。
点击查看答案
容易验证满足方程组
而且是该方程组的最小非负解。此外由条件可知也是该方程组的非负解。
若且存在某个使得,那么由不可约性,存在使得。由此可得
矛盾表明对所有都成立。
关于上述推导中第二个等号和最后一个等号的解释:
马氏链简单应用
设 G-W 分支过程的后代均值为,后代方差为,假定,求存活过粒子数的方差。
点击查看答案
设存活过粒子数为,则因为,所以。而根据条件方差公式:
又因为,故
综上可得,解得。
进一步可得到二阶矩。
布朗运动
对任意,参考以下结论,求的联合密度函数,其中。
对任意,的联合概率密度函数
点击查看答案
因为与等价,所以的分布函数
又已知的联合概率密度函数,令得
令得
因此的联合密度函数
不过本题还有一个相对简单的证明方法(不使用上面的结论):
的密度函数为。由布朗运动的马氏性和独立增量性可知,在给定的条件下,的分布
因此的条件密度为,其中时退化为。
的联合密度函数。
对任意,记为标准布朗运动首次访问的时间,以及。
(1)证明:。
(2)计算。
点击查看答案
(1)当时,,因此。下面先证明的情形:
利用布朗运动的独立增量性容易证明对任意正整数,可以看成是个独立同分布的之和,从而;类似,对任意正整数,可以看作是个独立同分布的之和,因此,从而。因此对任意有理数.
由布朗运动轨道连续性可知随增加单调不降,因此随单调下降。对任意无理数,存在有理数列使得且。从而
由此可得对任意成立。
又因为也是标准布朗运动(定义为),则其首次访问的时间(与同分布)。所以
对于的情形,因为与同分布,所以。
综上,。
(2)因为的密度函数为,因此
当然也可以用鞅的停时定理计算:可以验证为鞅,且代入后满足停时定理条件,故
因此
设是指数布朗运动,即,求。
点击查看答案
因为,下面进行分类讨论:
当时,,因此
进而。
当时,,因此
进而。
点击查看答案
方法一:由于非负,非负。对任意,
令两式相等,可得。
方法二:对任意,
由此可知。
离散时间鞅
若是鞅,则对任意,。
推导如下:
设关于是鞅,为一停时,证明:
(1)关于也是鞅;
(2),且单调不降。
(3)当时,。
点击查看答案
(1)因为,故可积;
又因为
因此关于也是鞅。
(2)因为是有界停时,故满足有界停时定理,因此。
注:不一定成立(需满足一般停时定理条件)
又因为关于是鞅,所以关于是下鞅,故单调不降。
(3)证明:记,,。那么
由及重期望公式可知对任意,
注意到
上式右端非负,且随的增加单调不增,而且由条件可知
这表明
由的任意性可知
因此
已知独立且都服从参数为的泊松分布。令。
(1)时求。
(2)时,令,对任意,求。
点击查看答案
令。由于是右连续的随机游动,因此若,则。
注意到,对任意,,记。
(1)即求。当时,存在使得,且对任意,。容易验证关于是鞅。由有界停时定理
注意到时,。令得。
再令,得。由此可得。
(2)当时,重复上面讨论。注意到,而且对任意,。即此时上面的,因此。由此可得。取使得
那么。
