首页 > 其他分享 >【习题】随机变量与分布函数

【习题】随机变量与分布函数

时间:2024-03-13 11:23:26浏览次数:29  
标签:xi frac 函数 sigma eta rho theta 习题 随机变量

[T0301] 设随机变量 \(\xi\) 取值于 \([0,1]\), 若 \(P\{x\le\xi<y\}\) 只与长度 \(y-x\) 有关 (对一切 \(0\le x\le y\le 1\)). 试证 \(\xi\sim U[0,1]\).

不妨设 \(P\{x\le\xi<y\}=f(y-x)\). 令 \(x=0\), 则有 \(P\{0\le\xi<y\}=f(y)\). 注意到对 \(\forall 0\le y_1<y_2\le1\), 有

\[f(y_1+y_2)=P\{0\le\xi<y_1+y_2\}=P\{0\le\xi<y_1\}+P\{y_1\le\xi<y_1+y_2\}=f(y_1)+f(y_2) \]

于是 \(e^{f(y_1+y_2)}=e^{f(y_1)}\cdot e^{f(y_2)}\), 由 Cauchy 命题知 \(e^{f(x)}=a^x\Rightarrow f(x)=x\ln a\), 其中 \(a\) 待定. 又 \(f(1)=1\), 故 \(a=e\), 从而 \(f(x)=x\), 故

\[P\{0\le\xi<x\}=x, \quad x\in[0,1] \]

即 \(\xi\sim U[0,1]\). #


[T0302] 设 \((\xi,\eta)\) 服从二元正态分布, 参数为 \(\mu_1,\mu_2,\sigma_1^2,\sigma_2^2,\rho\), 以 \(D(\lambda)\) 记下面椭圆的内部:

\[\frac{(x-\mu_1)^2}{\sigma_1^2}-\frac{2\rho(x-\mu_1)(y-\mu_2)}{\sigma_1\sigma_2}+\frac{(y-\mu_2)^2}{\sigma_2^2}=\lambda^2 \]

试求 \(P\{(\xi,\eta)\in D(\lambda)\}\).

即求二重积分

\[\iint\limits_{(x,y)\in D(\lambda)}\frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\exp\left\{-\frac{1}{2(1-\rho^2)}\left[\frac{(x-\mu_1)^2}{\sigma_1^2}-\frac{2\rho(x-\mu_1)(y-\mu_2)}{\sigma_1\sigma_2}+\frac{(y-\mu_2)^2}{\sigma_2^2}\right]\right\}\mathrm{d}x\mathrm dy \]

坐标变换 \(T:\begin{cases}x=\mu_1+\sigma_1r\cos\theta\\ y=\mu_2+\sigma_2r\sin\theta\end{cases}\), 则 \(J=\sigma_1\sigma_2r\). 经过简单的计算, 上面的二重积分可化为

\[\frac{\sqrt{1-\rho^2}}{2\pi}\left(1-e^{-\frac{\lambda^2}{2(1-\rho^2)}}\right)\int_0^{2\pi}\frac{1}{1-2\rho\sin\theta\cos\theta}\mathrm{~d}\theta \]

注意到 \(\lim\limits_{\lambda\to\infty}P\{(\xi,\eta)\in D(\lambda)\}=1\), 因此

\[\frac{\sqrt{1-\rho^2}}{2\pi}\int_0^{2\pi}\frac{1}{1-2\rho\sin\theta\cos\theta}\mathrm{~d}\theta=1\Rightarrow \int_0^{2\pi}\frac{1}{1-2\rho\sin\theta\cos\theta}\mathrm{~d}\theta=\frac{2\pi}{\sqrt{1-\rho^2}} \]

\[P\{(\xi,\eta)\in D(\lambda)\}=1-e^{-\frac{\lambda^2}{2(1-\rho^2)}}.\quad\quad\# \]


[T0303] 若 \(\xi_1,\cdots,\xi_n\) 相互独立同分布于 \(N(0,1)\), 证明 \(\eta=\xi_1^2+\cdots+\xi_n^2\) 服从 \(\chi^2\) 分布.

\[P\{\eta<x\}=P\{\xi_1^2+\cdots+\xi_n^2<x\}=\frac{1}{(2\pi)^{n/2}}{\int\cdots\int}_{x_1^2+\cdots+x_n^2<x}\exp\left\{-\frac{x_1^2+\cdots+x_n^2}{2}\right\}\mathrm dx_1\cdots\mathrm dx_n \]

作变换

\[\begin{cases} x_1=\rho\sin\theta_1\sin\theta_2\cdots\sin\theta_{n-1}\\ x_2=\rho\cos\theta_1\sin\theta_2\cdots\sin\theta_{n-1}\\ \cdots\cdots\\ x_{n-1}=\rho\cos\theta_{n-2}\sin\theta_{n-1}\\ x_n=\rho\cos\theta_{n-1} \end{cases} \]

则 \(J=\frac{\partial(x_1,\cdots,x_n)}{\partial(\rho,\theta_1,\cdots,\theta_{n-1})}=\rho^{n-1}D(\theta_1,\cdots,\theta_{n-1})\), 其中 \(D(\theta_1,\cdots,\theta_{n-1})\) 是 \(\theta_1,\cdots,\theta_{n-1}\) 的函数, 与 \(\rho\) 无关. 于是

\[\begin{aligned} F(x)&=\frac{1}{(2\pi)^{n/2}}\int_{0}^{\sqrt x}e^{-\frac{\rho^2}{2}}\rho^{n-1}\mathrm d\rho \ \int_0^{\pi}\cdots\int_0^{\pi}\int_0^{2\pi}D(\theta_1,\cdots,\theta_{n-1})\mathrm d\theta_1\cdots\mathrm d\theta_{n-1}\\ &=C_n\int_{0}^{\sqrt x}e^{-\frac{\rho^2}{2}}\rho^{n-1}\mathrm d\rho \end{aligned} \]

其中

\[C_n=\frac{1}{(2\pi)^{n/2}} \int_0^{\pi}\cdots\int_0^{\pi}\int_0^{2\pi}D(\theta_1,\cdots,\theta_{n-1})\mathrm d\theta_1\cdots\mathrm d\theta_{n-1} \]

令 \(\rho=\sqrt{t}\), 得

\[F(x)=\frac{C_n}{2}\int_0^xe^{-\frac12t}\cdot t^{\frac n2-1}\mathrm dt \]

再令 \(x\to+\infty\), 得

\[1=\frac{C_n}{2}\int_0^{+\infty}e^{-\frac12t}\cdot t^{\frac n2-1}\mathrm dt\xlongequal{t=2u}2^{\frac n2-1}C_n\int_0^{+\infty}e^{-u}\cdot u^{\frac n2-1}\mathrm du=2^{\frac n2-1}C_n\Gamma\left(\frac n2\right) \]

\[C_n=\frac{1}{2^{\frac n2-1}\Gamma\left(\frac n2\right)} \]

于是 \(\eta\) 的分布函数为

\[F(x)=\frac{1}{2^{\frac n2-1}\Gamma\left(\frac n2\right)}\int_0^xe^{-\frac12t}\cdot t^{\frac n2-1}\mathrm dt, \quad x>0. \]

即 \(\eta\) 服从 \(\chi^2\) 分布. #


[T0304] 设随机变量 \(\xi,\eta\) 独立同分布于 \(N(0,1)\). 试求 \(U=\xi^2+\eta^2\) 与 \(V=\frac{\xi}{\eta}\) 的分布密度函数, 并证明他们是独立的.

易知

\[p_{\xi\eta}(x,y)=\frac{1}{2\pi}e^{-\frac{x^2+y^2}{2}} \]

设 \(u=x^2+y^2, \ v=\frac{x}{y}\), 则 \(J=-\frac{1}{2}\frac{1}{1+v^2}\). 注意到 \((x,y)\) 与 \((-x,-y)\) 对应于同一个 \((u,v)\), 故

\[\begin{aligned} p_{UV}(u,v)&=p(x(u,v),y(u,v))|J|\binom{-\infty<x<\infty}{y\ge0}+p(x(u,v),y(u,v))|J|\binom{-\infty<x<\infty}{y<0}\\ &=\frac{1}{2}e^{-\frac u2}\cdot\frac{1}{\pi(1+v^2)}:=p_U(u)\cdot p_V(v) \end{aligned} \]

故 \(U\) 与 \(V\) 相互独立, 且

\[p_U(u)=\frac{1}{2}e^{-\frac u2},\quad u>0\\ p_V(v)=\frac{1}{\pi(1+v^2)}, -\infty<v<\infty \]

即 \(U\) 服从参数为 \(\frac12\) 的指数分布, \(V\) 服从 Cauchy 分布. #


[T0305] 若 \((\xi,\eta)\) 服从参数为 \(\mu_1,\mu_2,\sigma_1^2,\sigma_2^2,\rho\) 的二元正态分布, 试找出 \(\xi+\eta\) 与 \(\xi-\eta\) 相互独立的充要条件.

显然 \(\xi+\eta\) 与 \(\xi-\eta\) 相互独立 \(\Longleftrightarrow\) \(U=(\xi-\mu_1)+(\eta-\mu_2)\) 与 \(V=(\xi-\mu_1)-(\eta-\mu_2)\) 相互独立. 设

\[\binom uv=\pmatrix{1&1\\1&-1}\binom{x-\mu_1}{y-\mu_2}\Rightarrow\binom{x-\mu_1}{y-\mu_2}=\frac12\pmatrix{1&1\\1&-1}\binom uv \]

则 \(J=\frac{\partial(x,y)}{\partial(u,v)}=\frac12\). 于是 \((U,V)\) 的密度函数为

\[\begin{aligned} q(u,v)&=\frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\cdot\frac12\cdot\exp\left\{-\frac12(u,v)\left[\pmatrix{1&1\\1&-1}\pmatrix{\sigma_1^2&\rho\sigma_1\sigma_2\\\rho\sigma_1\sigma_2&\sigma_2^2}\pmatrix{1&1\\1&-1}\right]^{-1}\binom uv\right\}\\ &=\frac{1}{4\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\exp\left\{-\frac12(u,v)\pmatrix{\sigma_1^2+2\rho\sigma_1\sigma_2+\sigma_2^2&\sigma_1^2-\sigma_2^2\\\sigma_1^2-\sigma_2^2&\sigma_1^2-2\rho\sigma_1\sigma_2+\sigma_2^2}^{-1}\binom uv\right\} \end{aligned} \]

故 \(\xi+\eta\) 与 \(\xi-\eta\) 相互独立 \(\Longleftrightarrow\) \(\sigma_1^2-\sigma_2^2=0 \Longleftrightarrow \sigma_1=\sigma_2\). #


[T0306] 设随机变量 \(\xi\) 与 \(\eta\) 独立同分布, 其密度函数不为 \(0\) 且有二阶导数, 试证: 若 \(\xi+\eta\) 与 \(\xi-\eta\) 相互独立, 则随机变量 \(\xi,\eta,\xi+\eta,\xi-\eta\) 均服从正态分布.

设 \(U=\xi+\eta, \ V=\xi-\eta\), 则

\[p_U(u)p_V(v)=p_{UV}(u,v)=\frac12p_{\xi,\eta}(\frac{u+v}{2},\frac{u-v}{2})=\frac12p(\frac{u+v}{2})p(\frac{u-v}{2}) \]

其中 \(p\) 为 \(\xi\) 和 \(\eta\) 的密度. 两边取对数得

\[\ln p(\frac{u+v}{2})+\ln p(\frac{u-v}{2})-\ln2=\ln p_U(u)+\ln p_V(v) \]

记 \(f(x)=\ln p(x)\), 则 \(f(\frac{u+v}{2})+f(\frac{u-v}{2})-\ln2=\ln p_U(u)+\ln p_V(v)\). 两边对 \(u\) 求偏导, 得

\[\frac12f'(\frac{u+v}{2})+\frac12f'(\frac{u-v}{2})=\frac{p_U'(u)}{p_U(u)} \]

两边再对 \(v\) 求偏导, 得

\[\frac14f''(\frac{u+v}{2})-\frac14f''(\frac{u-v}{2})=0 \]

令 \(u=v=x\) 则 \(f''(x)=f''(0):=-\lambda\), 得 \(f(x)=-\frac{\lambda}{2}x^2+c_1x+c_2\), 于是

\[\ln p(x)=-\frac{\lambda}{2}x^2+c_1x+c_2\Rightarrow p(x)=e^{c_2+\frac{c_1^2}{2\lambda}}e^{-\frac{\lambda}{2}(x-\frac{c_1}{\lambda})^2}:=Ke^{-\frac{\lambda}{2}(x-\frac{c_1}{\lambda})^2} \]

因 \(p(x)\) 是密度函数, 故 \(\lambda>0\). 且 \(c_1=f'(0)=\frac{p'(0)}{p(0)}, \ c_2=f(0)=\ln p(0)\). 故 \(\xi,\eta\) 服从正态分布, 从而 \(\xi+\eta,\xi-\eta\) 也服从正态分布. #


[T0307] 设随机变量 \(\xi\) 和 \(\eta\) 相互独立, 且 \(\xi\sim B(n,p), \ \eta\sim U(0,1)\), 试求 \(\xi+\eta\) 的分布函数和密度函数.

对 \(x\le 0\), 有 \(P\{\xi+\eta<x\}=0\); 对 \(x>n+1\), 有 \(P\{\xi+\eta<x\}=1\); 对 \(k<x\le k+1, \ k=0,1,2,\cdots,n\), 有

\[\begin{aligned} P\{\xi+\eta<x\}&=P\{\xi=k,\xi+\eta<x\}+P\{\xi\le k-1,\xi+\eta<x\}+P\{\xi\ge k+1,\xi+\eta<x\}\\ &=P\{\xi=k,\eta<x-k\}+P\{\xi\le k-1\}\\ &=\binom{n}{k}p^k(1-p)^{n-k}(x-k)+\sum_{\mathscr{l}=0}^{k-1}\binom{n}{\mathscr{l}}p^{\mathscr{l}}(1-p)^{n-\mathscr{l}} \end{aligned} \]

于是

\[F(x)=\begin{cases} 0,& x\le0\\ \binom{n}{k}p^k(1-p)^{n-k}(x-k)+\sum_{\mathscr{l}=0}^{k-1}\binom{n}{\mathscr{l}}p^{\mathscr{l}}(1-p)^{n-\mathscr{l}}, &k<x\le k+1, \ k=0,1,2,\cdots,n\\ 1, & x>n+1 \end{cases} \]

从而密度函数为

\[p(x)=\begin{cases} \binom{n}{k}p^k(1-p)^{n-k},&k<x\le k+1, \ k=0,1,2,\cdots,n\\ 0,& x\le0, \ x>n+1. \end{cases}\quad\quad\# \]


[T0308] 试求顺序统计量 \(\xi_k^*\) 与 \(\xi_l^* \ (k<l)\) 的联合密度函数.

对 \(x<y\), 有

\[\begin{aligned} &P\{x\le\xi_k^*<x+\Delta x, y\le\xi_l^*<y+\Delta y\}\\ =&\frac{n!}{(k-1)!1!(l-k-1)!1!(n-l)!}F(x)^{k-1}\cdot p(x)\Delta x\cdot[F(y)-F(x+\Delta x)]^{l-k-1}\cdot p(y)\Delta y\cdot[1-F(y+\Delta y)]^{n-l} \end{aligned} \]

于是 \(\xi_k^*\) 与 \(\xi_l^* \ (k<l)\) 的联合密度函数为

\[p(x,y)= \begin{cases} \frac{n!}{(k-1)!(l-k-1)!(n-l)!}F(x)^{k-1}\cdot [F(y)-F(x)]^{l-k-1}\cdot[1-F(y)]^{n-l}\cdot p(x)\cdot p(y),&x<y\\ 0,&x\ge y \end{cases} \]

#

标签:xi,frac,函数,sigma,eta,rho,theta,习题,随机变量
From: https://www.cnblogs.com/hznudmh/p/18070212

相关文章

  • SqlServer函数大全三十八:DATEPART函数
    在SQLServer中,DATEPART 函数用于返回日期/时间值的指定部分的整数。与 DATENAME 函数不同,DATEPART 返回的是一个数字,而不是一个字符串。这对于需要进行数学计算或比较的场合特别有用。函数的语法如下:sql复制代码DATEPART(datepart,date)其中:datepart 是你想......
  • SqlServer函数大全三十九:CONVERT函数
    在SQLServer中,CONVERT 函数用于将一种数据类型转换为另一种数据类型。这在处理日期、时间、数字和其他数据类型时非常有用,尤其是当你需要确保数据以特定的格式或类型进行存储或显示时。函数的语法如下:sql复制代码CONVERT(data_type[(length)],expression[,style])......
  • SqlServer函数大全三十五:DATEDIFF(返回日期和时间的边界数)函数
    在SQLServer中,DATEDIFF 函数用于返回两个日期之间的边界数差异。这个函数可以计算两个日期之间的年、月、日、小时、分钟、秒或周数差异。DATEDIFF 函数的语法如下:sql复制代码DATEDIFF(datepart,startdate,enddate)datepart 是指定要返回日期部分的参数,比如......
  • linux Shell 命令行-07-func 函数
    拓展阅读linuxShell命令行-00-intro入门介绍linuxShell命令行-02-var变量linuxShell命令行-03-array数组linuxShell命令行-04-operator操作符linuxShell命令行-05-test验证是否符合条件linuxShell命令行-06-flowcontrol流程控制linuxShell命令行-07-f......
  • Vue3 组合函数 element-plus table数据滚动播放
    Vue滚动播放组合函数import{onMounted,onUnmounted}from"vue";exportfunctioncreateScroll(tableRef){lettimer=null;functionstartScroll(){consttable=tableRef.value.layout.table.refs;consttableWrapper=table.bodyWrapper.f......
  • 为什么defineProps宏函数不需要从vue中import导入?
    前言我们每天写vue代码时都在用defineProps,但是你有没有思考过下面这些问题。为什么defineProps不需要import导入?为什么不能在非setup顶层使用defineProps?defineProps是如何将声明的props自动暴露给模板?举几个例子我们来看几个例子,分别对应上面的几个问题。先来看一个正常的......
  • C语言字符函数和字符串函数
    前言今天这篇博客咱们一起来认识一些特殊的函数,在编程的过程中,我们经常要处理字符和字符串,为了方便字符和字符串,C语言提供了一些库函数,让我们一起看看这些函数都有什么功能吧!!!个人主页:小张同学zkf若有问题评论区见感兴趣就关注一下吧目录 1.字符分类函数2.字符......
  • 函数
    一、字符串函数常用的几个如下:函数功能CONCAT(s1,s2,…,sn)字符串拼接,将s1,s2,…,sn拼接成一个字符串LOWER(str)将字符串全部转为小写UPPER(str)将字符串全部转为大写LPAD(str,n,pad)左填充,用字符串pad对str的左边进行填充,达到n个字符串长度RPAD(......
  • c语言函数传递数组名
    c语言自定义函数中可以在形参中可以使用数组名作为传递代码示例如下#include<stdio.h>floatave(floata[]){ inti; floatb; floatsum=a[0]; for(i=1;i<10;++i) sum=sum+a[i]; b=sum/10; returnb;}intmain(){ floatnum[10],average; inti; for(i=0;i......
  • file_put_content函数下的绕过
    该函数的作用是将文件写入字符串中首先,file_put_content大概有三种情形出现;file_put_contents($filename,"<?phpexit();".$content);file_put_contents($content,"<?phpexit();".$content);file_put_contents($filename,$content."\nxxxxxx"); 思路:将杂糅或者......