LOADING
730 字
4 分钟
例子2.1.1

解答:求 X1X_1X2X_2 的边缘密度函数

给定二维随机向量 X=(X1,X2)X = (X_1, X_2) 的联合密度函数:

f(x1,x2)=12πe12(x12+x22)[1+x1x2e12(x12+x22)].f(x_1, x_2) = \frac{1}{2\pi} e^{-\frac{1}{2}(x_1^2 + x_2^2)} \left[ 1 + x_1 x_2 e^{-\frac{1}{2}(x_1^2 + x_2^2)} \right].

目标是求 X1X_1X2X_2 的边缘密度函数 fX1(x1)f_{X_1}(x_1)fX2(x2)f_{X_2}(x_2)。由于联合密度函数关于 x1x_1x2x_2 对称(交换 x1x_1x2x_2 表达式不变),只需计算 fX1(x1)f_{X_1}(x_1)fX2(x2)f_{X_2}(x_2) 可由对称性直接得出。


步骤 1:展开联合密度函数

将联合密度拆分为两项:

f(x1,x2)=12πe12(x12+x22)f1(x1,x2)+12πx1x2e(x12+x22)f2(x1,x2).f(x_1, x_2) = \underbrace{\frac{1}{2\pi} e^{-\frac{1}{2}(x_1^2 + x_2^2)}}_{f_1(x_1, x_2)} + \underbrace{\frac{1}{2\pi} x_1 x_2 e^{-(x_1^2 + x_2^2)}}_{f_2(x_1, x_2)}.

其中:

  • f1f_1 对应标准二元正态分布的密度(独立标准正态变量的联合密度);
  • f2f_2 是扰动项,含交叉项 x1x2x_1 x_2

步骤 2:计算 X1X_1 的边缘密度 fX1(x1)f_{X_1}(x_1)

边缘密度定义为:

fX1(x1)=f(x1,x2)dx2=f1(x1,x2)dx2+f2(x1,x2)dx2.f_{X_1}(x_1) = \int_{-\infty}^{\infty} f(x_1, x_2) dx_2 = \int_{-\infty}^{\infty} f_1(x_1, x_2) dx_2 + \int_{-\infty}^{\infty} f_2(x_1, x_2) dx_2.
(a) 计算第一部分 f1(x1,x2)dx2\int_{-\infty}^{\infty} f_1(x_1, x_2) dx_2
12πe12(x12+x22)dx2=12πex122ex222dx2.\int_{-\infty}^{\infty} \frac{1}{2\pi} e^{-\frac{1}{2}(x_1^2 + x_2^2)} dx_2 = \frac{1}{2\pi} e^{-\frac{x_1^2}{2}} \int_{-\infty}^{\infty} e^{-\frac{x_2^2}{2}} dx_2.

已知高斯积分 ex2/2dx=2π\int_{-\infty}^{\infty} e^{-x^2/2} dx = \sqrt{2\pi},代入得:

12πex1222π=12πex122.\frac{1}{2\pi} e^{-\frac{x_1^2}{2}} \cdot \sqrt{2\pi} = \frac{1}{\sqrt{2\pi}} e^{-\frac{x_1^2}{2}}.

此即标准正态分布 N(0,1)\mathcal{N}(0, 1) 的密度函数。

(b) 计算第二部分 f2(x1,x2)dx2\int_{-\infty}^{\infty} f_2(x_1, x_2) dx_2
12πx1x2e(x12+x22)dx2=x12πex12x2ex22dx2.\int_{-\infty}^{\infty} \frac{1}{2\pi} x_1 x_2 e^{-(x_1^2 + x_2^2)} dx_2 = \frac{x_1}{2\pi} e^{-x_1^2} \int_{-\infty}^{\infty} x_2 e^{-x_2^2} dx_2.

分析被积函数 g(x2)=x2ex22g(x_2) = x_2 e^{-x_2^2}

  • g(x2)=x2ex22=g(x2)g(-x_2) = -x_2 e^{-x_2^2} = -g(x_2),故 g(x2)g(x_2)奇函数
  • 积分区间 (,)(-\infty, \infty) 关于原点对称;
  • x2ex22\lvert x_2 e^{-x_2^2} \rvertR\mathbb{R} 上绝对可积(指数衰减主导)。

因此:

x2ex22dx2=0.\int_{-\infty}^{\infty} x_2 e^{-x_2^2} dx_2 = 0.

第二部分积分为 0

(c) 合并结果
fX1(x1)=12πex122+0=12πex122.f_{X_1}(x_1) = \frac{1}{\sqrt{2\pi}} e^{-\frac{x_1^2}{2}} + 0 = \frac{1}{\sqrt{2\pi}} e^{-\frac{x_1^2}{2}}.

步骤 3:由对称性得 X2X_2 的边缘密度

联合密度 f(x1,x2)f(x_1, x_2) 满足 f(x1,x2)=f(x2,x1)f(x_1, x_2) = f(x_2, x_1)(交换 x1x_1x2x_2 表达式不变),故:

fX2(x2)=f(x1,x2)dx1=12πex222.f_{X_2}(x_2) = \int_{-\infty}^{\infty} f(x_1, x_2) dx_1 = \frac{1}{\sqrt{2\pi}} e^{-\frac{x_2^2}{2}}.

步骤 4:验证与讨论

  1. 非负性验证(确保联合密度有效):
    括号内项 1+x1x2e12(x12+x22)1 + x_1 x_2 e^{-\frac{1}{2}(x_1^2 + x_2^2)} 的最小值分析:

    • u=x1x2e12(x12+x22)u = \lvert x_1 x_2 \rvert e^{-\frac{1}{2}(x_1^2 + x_2^2)},由 AM-GM 不等式 x12+x222x1x2x_1^2 + x_2^2 \geq 2\lvert x_1 x_2 \rvert,得 ux1x2ex1x2u \leq \lvert x_1 x_2 \rvert e^{-\lvert x_1 x_2 \rvert}
    • 函数 h(t)=teth(t) = t e^{-t}t0t \geq 0)的最大值为 h(1)=1/e0.3679<1h(1) = 1/e \approx 0.3679 < 1,故 1u11/e>01 - u \geq 1 - 1/e > 0
    • 因此 f(x1,x2)0f(x_1, x_2) \geq 0 恒成立,联合密度有效。
  2. 联合密度积分为 1 的验证

    R2f(x1,x2)dx1dx2=fX1(x1)dx1=12πex12/2dx1=1,\iint_{\mathbb{R}^2} f(x_1, x_2) dx_1 dx_2 = \int_{-\infty}^{\infty} f_{X_1}(x_1) dx_1 = \int_{-\infty}^{\infty} \frac{1}{\sqrt{2\pi}} e^{-x_1^2/2} dx_1 = 1,

    符合概率密度函数要求(Fubini 定理适用,因被积函数非负且边缘积分收敛)。

  3. 关键观察

    • 尽管 X1X_1X2X_2 不独立(联合密度含交叉项 x1x2x_1 x_2),但扰动项 f2f_2 在边缘积分中因奇函数性质消失。
    • 此例说明:边缘分布为标准正态,不能推出联合分布为二元正态或变量独立。本题是构造非独立但边缘为标准正态的经典例子。

最终答案

X1X_1X2X_2 的边缘密度函数均为标准正态分布:

fX1(x1)=12πex122,fX2(x2)=12πex222\boxed{f_{X_1}(x_1) = \frac{1}{\sqrt{2\pi}} e^{-\frac{x_1^2}{2}}, \quad f_{X_2}(x_2) = \frac{1}{\sqrt{2\pi}} e^{-\frac{x_2^2}{2}}}

X1N(0,1)X_1 \sim \mathcal{N}(0, 1)X2N(0,1)X_2 \sim \mathcal{N}(0, 1)

例子2.1.1
/posts/study/multivariate-statistical-analysis/例子211/
作者
Xs
发布于
2026-04-29
许可协议
CC BY-NC-SA 4.0

部分信息可能已经过时