跳到主要内容

大阪大学 基礎工学研究科 数理科学 (システム創成専攻) 2018年8月実施 数理科学 I [5]

Author

Miyake

Description

题目描述

原文题干缺失。设随机变量 X1,,XnX_1,\ldots,X_n 满足

E[Xi]=μ,V[Xi]=σ2,Cov(Xi,Xj)=θ(ij).E[X_i]=\mu,\qquad V[X_i]=\sigma^2,\qquad \operatorname{Cov}(X_i,X_j)=\theta\quad(i\ne j).

对系数向量 α=(α1,,αn)\alpha=(\alpha_1,\ldots,\alpha_n),定义线性统计量 T(α)=i=1nαiXi.T(\alpha)=\sum_{i=1}^n\alpha_iX_i.

  1. E[T(α)]E[T(\alpha)] 以及 Cov(T(α),T(β))\operatorname{Cov}(T(\alpha),T(\beta))
  2. V[T(α)]V[T(\alpha)]
  3. T(α)T(\alpha)μ\mu 的无偏估计量的约束下,求使方差最小的 α\alpha,并给出最小方差。

Kai

(1)

E[T(α)]=i=1nαiE[Xi]=μi=1nαiCov[T(α),T(β)]=i=1nj=1nαiβjCov[Xi,Xj]=i=1nαiβiV[Xi]+ijαiβjCov[Xi,Xj]=σ2i=1nαiβi+θijαiβj=(σ2θ)i=1nαiβi+θi=1nj=1nαiβj=(σ2θ)i=1nαiβi+θi=1nαij=1nβj \begin{aligned} E [ T(\alpha) ] &= \sum_{i=1}^n \alpha_i E[X_i] \\ &= \mu \sum_{i=1}^n \alpha_i \\ \text{Cov} [ T(\alpha), T(\beta) ] &= \sum_{i=1}^n \sum_{j=1}^n \alpha_i \beta_j \text{Cov} [X_i, X_j] \\ &= \sum_{i=1}^n \alpha_i \beta_i V [X_i] + \sum_{i \neq j} \alpha_i \beta_j \text{Cov} [X_i, X_j] \\ &= \sigma^2 \sum_{i=1}^n \alpha_i \beta_i + \theta \sum_{i \neq j} \alpha_i \beta_j \\ &= \left( \sigma^2 - \theta \right) \sum_{i=1}^n \alpha_i \beta_i + \theta \sum_{i=1}^n \sum_{j=1}^n \alpha_i \beta_j \\ &= \left( \sigma^2 - \theta \right) \sum_{i=1}^n \alpha_i \beta_i + \theta \sum_{i=1}^n \alpha_i \sum_{j=1}^n \beta_j \end{aligned}

(2)

V[T(α)]=Cov[T(α),T(α)]=(σ2θ)i=1nαi2+θi=1nαij=1nαj=(σ2θ)i=1nαi2+θ{i=1nαi}2 \begin{aligned} V[T(\alpha)] &= \text{Cov} [ T(\alpha), T(\alpha) ] \\ &= \left( \sigma^2 - \theta \right) \sum_{i=1}^n \alpha_i^2 + \theta \sum_{i=1}^n \alpha_i \sum_{j=1}^n \alpha_j \\ &= \left( \sigma^2 - \theta \right) \sum_{i=1}^n \alpha_i^2 + \theta \left\{ \sum_{i=1}^n \alpha_i \right\}^2 \end{aligned}

(3)

T(α)T(\alpha)μ\mu の不偏推定量であるための条件は、(1) より、

i=1nαi=1 \begin{aligned} \sum_{i=1}^n \alpha_i = 1 \end{aligned}

であり、このとき、

V[T(α)]=(σ2θ)i=1nαi2+θ \begin{aligned} V[T(\alpha)] = \left( \sigma^2 - \theta \right) \sum_{i=1}^n \alpha_i^2 + \theta \end{aligned}

となる。

これを最小化するために、 ラグランジュの未定乗数 λ\lambda を導入して、

f(α)=i=1nαi2λi=1nαi \begin{aligned} f(\alpha) = \sum_{i=1}^n \alpha_i^2 - \lambda \sum_{i=1}^n \alpha_i \end{aligned}

を考えると、

fαi=2αiλ \begin{aligned} \frac{\partial f}{\partial \alpha_i} = 2 \alpha_i - \lambda \end{aligned}

であるから、

α1=α2==αn=1n,  λ=2n \begin{aligned} \alpha_1 = \alpha_2 = \cdots = \alpha_n = \frac{1}{n} , \ \ \lambda = \frac{2}{n} \end{aligned}

のとき、 f(α)f(\alpha) したがって V[T(α)]V[T(\alpha)] が最小になる。

このとき、

V[T(α)]=(σ2θ)1n+θ=σ2+(n1)θn \begin{aligned} V[T(\alpha)] &= \left( \sigma^2 - \theta \right) \cdot \frac{1}{n} + \theta \\ &= \frac{\sigma^2 + (n-1) \theta}{n} \end{aligned}

である。