← 首页

Jacobi 行列式、概率论与卷积

Jacobi 行列式与二重积分的一般换元法

FFGG 是关于变量 xxyy 的可微函数,下面这个行列式称为函数 FFGG 关于 xxyyJacobi 行列式

(F,G)(x,y)=FxFyGxGy\frac{\partial(F, G)}{\partial(x, y)} = \begin{vmatrix} \frac{\partial F}{\partial x} & \frac{\partial F}{\partial y} \\[4pt] \frac{\partial G}{\partial x} & \frac{\partial G}{\partial y} \end{vmatrix} \notag

如果变换

T:  {x=x(u,v)y=y(u,v)T:\;\begin{cases} x = x(u, v) \\ y = y(u, v) \end{cases} \notag

建立了 uOvuOv 平面上的区域 DD' xOyxOy 平面上的区域 DD 之间的双射,且 x(u,v)x(u,v) y(u,v)y(u,v) 具有连续偏导数,Jacobi 行列式 (x,y)(u,v)\dfrac{\partial(x, y)}{\partial(u, v)} DD' 内不为零,则有二重积分的一般换元公式

Df(x,y)dxdy=Df(x(u,v),y(u,v))(x,y)(u,v)dudv\iint_D f(x,y)\,\mathrm{d}x\,\mathrm{d}y = \iint_{D'} f\big(x(u,v), y(u,v)\big) \left| \frac{\partial(x,y)}{\partial(u,v)} \right| \mathrm{d}u\,\mathrm{d}v \notag

二维连续型随机变量的分布

(X,Y)(X,Y) 的概率密度为 f(x,y)f(x,y)Z=g(X,Y)Z = g(X,Y),假设对任意固定的 xxyg(x,y)y \mapsto g(x,y) 关于 yy 严格单调,那么存在映射 hh 满足 y=h(x,z)y = h(x,z)。可通过下面的方法求 fZ(z)f_Z(z)

fZ(z)=dFZ(z)dz=limΔz0FZ(z+Δz)FZ(z)Δz=limΔz0P{zZz+Δz}Δz=limΔz01ΔzDz,Δz(xy)f(x,y)dxdy\begin{aligned} f_Z(z) = \frac{dF_Z(z)}{dz} &= \lim_{\Delta z \to 0} \frac{F_Z(z + \Delta z) - F_Z(z)}{\Delta z} \\[6pt] &= \lim_{\Delta z \to 0} \frac{P\{z \leq Z \leq z + \Delta z\}}{\Delta z}\\[6pt] &= \lim_{\Delta z \to 0} \frac{1}{\Delta z} \iint_{D_{z,\Delta z}^{(xy)}} f(x,y)\,dx\,dy \end{aligned}\notag

其中 Dz,Δz(xy)={(x,y)zg(x,y)z+Δz}D_{z,\Delta z}^{(xy)} = \{ (x,y) \mid z \leq g(x,y) \leq z + \Delta z \},做代换

T:  {x=uy=h(u,v)T:\;\begin{cases} x = u \\ y = h(u, v) \end{cases} \notag fZ(z)=limΔz01ΔzDz,Δz(xy)f(x,y)dxdy=limΔz01ΔzDz,Δz(uv)f(u,h(u,v))(x,y)(u,v)dudv=limΔz01ΔzDz,Δz(uv)f(u,h(u,v))det(xuxvyuyv)dudv=limΔz01ΔzDz,Δz(uv)f(u,h(u,v))h(u,v)vdudv\begin{aligned} f_Z(z) &= \lim_{\Delta z \to 0} \frac{1}{\Delta z} \iint_{D_{z,\Delta z}^{(xy)}} f(x,y)\,dx\,dy \\[8pt] &= \lim_{\Delta z \to 0}\frac{1}{\Delta z}\iint_{D_{z,\Delta z}^{'(uv)}} f(u,h(u,v)) \left|\frac{\partial{(x,y)}}{\partial(u,v)}\right|\,du\,dv\\[8pt] &= \lim_{\Delta z \to 0}\frac{1}{\Delta z}\iint_{D_{z,\Delta z}^{'(uv)}} f(u,h(u,v)) \left|\det \begin{pmatrix} \dfrac{\partial x}{\partial u} & \dfrac{\partial x}{\partial v} \\[8pt] \dfrac{\partial y}{\partial u} & \dfrac{\partial y}{\partial v} \end{pmatrix}\right|\,du\,dv\\[8pt] &= \lim_{\Delta z \to 0}\frac{1}{\Delta z}\iint_{D_{z,\Delta z}^{'(uv)}} f(u,h(u,v)) \left|\frac{\partial h(u,v)}{\partial v}\right|\,du\,dv\\[8pt] \end{aligned}\notag

上面的代换 TT 满足 g(x,y)=vg(x,y)=v,也就是说

Dz,Δz(xy)={(x,y)zg(x,y)z+Δz}D_{z,\Delta z}^{(xy)} = \{ (x,y) \mid z \leq g(x,y) \leq z + \Delta z \}\notag

uvuv 平面中变成

Dz,Δz(uv)={(u,v)zvz+Δz}D_{z,\Delta z}^{'(uv)} = \{ (u,v) \mid z \leq v \leq z + \Delta z \}\notag fZ(z)=limΔz01ΔzDz,Δz(uv)f(u,h(u,v))h(u,v)vdudv=limΔz01Δzzz+Δz[+f(u,h(u,v))h(u,v)vdu]dv=limΔz01Δzzz+ΔzΦ(v)dv=limΔz01ΔzΦ(ξ)Δz(积分中值定理)=Φ(z)=+f(u,h(u,z))h(u,z)zdu\begin{aligned} f_Z(z) &=\lim_{\Delta z\to 0} \frac{1}{\Delta z} \iint_{D_{z,\Delta z}^{\prime(uv)}} f(u,h(u,v)) \left| \frac{\partial h(u,v)}{\partial v} \right| \,du\,dv \\[8pt] &= \lim_{\Delta z\to 0} \frac{1}{\Delta z} \int_z^{z+\Delta z} \left[ \int_{-\infty}^{+\infty} f(u,h(u,v)) \left| \frac{\partial h(u,v)}{\partial v} \right| \,du \right]dv\\[8pt] &= \lim_{\Delta z \to 0} \frac{1}{\Delta z} \int_z^{z+\Delta z} \Phi(v) dv\\[8pt] &= \lim_{\Delta z\to 0} \frac{1}{\Delta z} \Phi(\xi)\Delta z\quad \text{(积分中值定理)}\\[8pt] &= \Phi(z) \\[8pt] &= \int_{-\infty}^{+\infty} f(u,h(u,z)) \left| \frac{\partial h(u,z)}{\partial z} \right| \,du \end{aligned}\notag

上面这个就是二维随机变量函数密度的一般公式

fZ(z)=+f(x,h(x,z))h(x,z)zdx\boxed{f_Z(z) = \int_{-\infty}^{+\infty} f(x, h(x,z)) \left|\frac{\partial h(x,z)}{\partial z}\right| dx}\notag

这个公式并不要求变量 XXYY 相互独立。

跟卷积的关联

如果 Z=X+YZ = X + Y,且 XXYY 相互独立,此时一般公式中的偏导数因子 h(x,z)z=1\left| \dfrac{\partial h(x,z)}{\partial z}\right | = 1 。一般公式可化简为

fZ(z)=+f(x,zx)dx=+fX(x)fY(zx)dx=(fXfY)(z)\begin{aligned} f_Z(z) &= \int_{-\infty}^{+\infty} f(x, z-x) dx \\[8pt] &= \int_{-\infty}^{+\infty} f_X(x)f_Y(z-x)dx \\[8pt] &= (f_X * f_Y) (z) \end{aligned}\notag

上式即为教材上常见的卷积公式。

← 首页