Skip to Content
МетодичкиМатематическая статистикаПриложениеП2. Математическое ожидание и дисперсия

П2. Математическое ожидание и дисперсия

Интеграл Лебега

Определение. Для случайной величины ξ≥0\xi \ge 0

Mξ=∫Ωξ(ω) P(dω)=lim⁡n→∞[∑i=1n2ni−12n P(i−12n≤ξ<i2n)+n P(ξ≥n)].\mathsf M\xi = \int_\Omega\xi(\omega)\,\mathsf P(d\omega) = \lim_{n\to\infty}\Bigl[\sum_{i=1}^{n2^n}\frac{i-1}{2^n}\,\mathsf P\Bigl(\frac{i-1}{2^n} \le \xi < \frac i{2^n}\Bigr) + n\,\mathsf P(\xi \ge n)\Bigr].

Эта формула означает, что если заменить ξ\xi дискретной случайной величиной ξ(n)\xi^{(n)} таким образом, что ξ(n)=(i−1)2−n\xi^{(n)} = (i-1)2^{-n} на множестве {ω ⁣:(i−1)2−n≤ξ<i2−n}\{\omega\colon (i-1)2^{-n} \le \xi < i2^{-n}\} (и ξ(n)=n\xi^{(n)} = n там, где ξ≥n\xi \ge n), то, вычисляя Mξ(n)\mathsf M\xi^{(n)} и устремляя n→∞n \to \infty, приходим к пределу, который и называется Mξ\mathsf M\xi.

По кирпичикам: ξ(n)\xi^{(n)} — это ξ\xi, «округлённая вниз» до ближайшего деления линейки с шагом 2−n2^{-n} и «обрезанная» на высоте nn. У такой величины конечное число значений, и её среднее — обычная сумма «значение × вероятность». С ростом nn линейка мельчает и удлиняется, ξ(n)\xi^{(n)} подползает к ξ\xi снизу, и суммы растут к пределу. Предел может оказаться и бесконечным.

Смысл по-пацански: Риман режет ось аргументов на кусочки и складывает «высота × ширина». Лебег режет ось значений: «сколько народу заработало от 0 до 100 рублей, сколько — от 100 до 200» — и складывает «уровень × сколько народу на этом уровне». Лебегу всё равно, что такое ω\omega — точки отрезка, расклады карт или траектории броуновского движения. Поэтому теория вероятностей и живёт на интеграле Лебега.

Определение. Для произвольной случайной величины ξ\xi

Mξ=Mξ+−Mξ−,где ξ+=max⁡{ξ,0}, ξ−=−min⁡{ξ,0},\mathsf M\xi = \mathsf M\xi^+ - \mathsf M\xi^-, \quad\text{где } \xi^+ = \max\{\xi, 0\},\ \xi^- = -\min\{\xi, 0\},

и хотя бы одно из математических ожиданий конечно. В случае, когда Mξ+=Mξ−=+∞\mathsf M\xi^+ = \mathsf M\xi^- = +\infty, Mξ\mathsf M\xi (интеграл Лебега) не существует.

Это и случай закона Коши (гл. 7): у него и «положительная», и «отрицательная» части имеют бесконечные средние, разность ∞−∞\infty - \infty не определена — и среднего нет.

Математическое ожидание как интеграл Лебега

4,20013n
n = 4: сумма = 0,951; Mξ = 1,000

Сумма из определения Mξ (П2): ось значений режется на ступеньки шириной 2⁻ⁿ до высоты n, ξ округляется вниз до ступеньки, а всё, что выше n, считается равным n. С ростом n сумма растёт и подходит к Mξ (пунктир). У |Коши| математического ожидания нет — сумма растёт без предела. Отличие от интеграла Римана: режем не ось аргумента ω, а ось значений ξ — поэтому ω может быть чем угодно, хоть исходами бросков монеты.

Замена переменных

Определение. Борелевской σ-алгеброй B\mathcal B называется наименьшая σ-алгебра, содержащая всевозможные интервалы на прямой.

Определение. Функция φ(x)\varphi(x) — борелевская, если она измерима относительно B\mathcal B (в частности, такова любая непрерывная функция).

Теорема о замене переменных. Для борелевской функции φ\varphi

Mφ(ξ)=∫Ωφ(ξ(ω)) P(dω)=∫−∞+∞φ(x) Fξ(dx),\mathsf M\varphi(\xi) = \int_\Omega\varphi(\xi(\omega))\,\mathsf P(d\omega) = \int_{-\infty}^{+\infty}\varphi(x)\,F_\xi(dx),

где Fξ(x)=P(ξ≤x)F_\xi(x) = \mathsf P(\xi \le x) — функция распределения случайной величины ξ\xi. (Справа стоит интеграл Стилтьеса.)

Следствие. Момент kk-го порядка Mξk=∫−∞+∞xkFξ(dx)\mathsf M\xi^k = \int_{-\infty}^{+\infty}x^kF_\xi(dx).

По кирпичикам: левая часть — интеграл по «невидимому» пространству исходов Ω\Omega, правая — по обычной числовой прямой с весами из функции распределения. Теорема говорит: чтобы посчитать среднее от функции случайной величины, Ω\Omega знать не нужно, достаточно закона распределения. Для величины с плотностью Fξ(dx)=p(x) dxF_\xi(dx) = p(x)\,dx, для дискретной — сумма по значениям с их вероятностями.

Свойства математического ожидания

  1. M(ξ+c)=Mξ+c\mathsf M(\xi + c) = \mathsf M\xi + c,   M(cξ)=c Mξ\;\mathsf M(c\xi) = c\,\mathsf M\xi.
  2. M(ξ+η)=Mξ+Mη\mathsf M(\xi + \eta) = \mathsf M\xi + \mathsf M\eta, если нет неопределённости вида ∞−∞\infty - \infty.
  3. Если ξ≤η\xi \le \eta, то Mξ≤Mη\mathsf M\xi \le \mathsf M\eta.
  4. ∣Mξ∣≤M∣ξ∣|\mathsf M\xi| \le \mathsf M|\xi|.
  5. M(ξη)=Mξ⋅Mη\mathsf M(\xi\eta) = \mathsf M\xi\cdot\mathsf M\eta, когда ξ\xi и η\eta независимы, а Mξ\mathsf M\xi и Mη\mathsf M\eta конечны.

Смысл по-пацански: свойство 2 работает всегда — хоть зависимые величины, хоть нет. Среднее суммы — сумма средних, и точка. А вот свойство 5 про произведение — только для независимых. Путать эти два — классическая ошибка первокурсника.

Дисперсия, ковариация, корреляция

Для случайных величин ξ\xi и η\eta таких, что Mξ2<∞\mathsf M\xi^2 < \infty и Mη2<∞\mathsf M\eta^2 < \infty, определены

  • дисперсия Dξ=M(ξ−Mξ)2=Mξ2−(Mξ)2\mathsf D\xi = \mathsf M(\xi - \mathsf M\xi)^2 = \mathsf M\xi^2 - (\mathsf M\xi)^2,
  • ковариация cov⁡(ξ,η)=M(ξ−Mξ)(η−Mη)=M(ξη)−Mξ Mη\operatorname{cov}(\xi, \eta) = \mathsf M(\xi - \mathsf M\xi)(\eta - \mathsf M\eta) = \mathsf M(\xi\eta) - \mathsf M\xi\,\mathsf M\eta,
  • коэффициент корреляции ρ(ξ,η)=cov⁡(ξ,η)/Dξ Dη\rho(\xi, \eta) = \operatorname{cov}(\xi, \eta)\big/\sqrt{\mathsf D\xi\,\mathsf D\eta}.

По кирпичикам — откуда Mξ2−(Mξ)2\mathsf M\xi^2 - (\mathsf M\xi)^2: раскроем квадрат, M(ξ2−2ξμ+μ2)=Mξ2−2μ Mξ+μ2=Mξ2−μ2\mathsf M(\xi^2 - 2\xi\mu + \mu^2) = \mathsf M\xi^2 - 2\mu\,\mathsf M\xi + \mu^2 = \mathsf M\xi^2 - \mu^2 при μ=Mξ\mu = \mathsf M\xi (свойства 1 и 2). Формула удобна для счёта, но на компьютере опасна: при больших μ\mu вычитаются два почти равных огромных числа и теряется точность.

Свойства дисперсии и ковариации

  1. D(ξ+c)=Dξ\mathsf D(\xi + c) = \mathsf D\xi,   D(cξ)=c2 Dξ\;\mathsf D(c\xi) = c^2\,\mathsf D\xi.
  2. Если ξ\xi и η\eta независимы, то D(ξ+η)=Dξ+Dη\mathsf D(\xi + \eta) = \mathsf D\xi + \mathsf D\eta.
  3. Для случайных величин ξ1,…,ξn\xi_1, \dots, \xi_n с конечными Mξi2\mathsf M\xi_i^2
D(∑i=1nξi)=∑i=1nDξi+2∑i<jcov⁡(ξi,ξj).\mathsf D\Bigl(\sum_{i=1}^n\xi_i\Bigr) = \sum_{i=1}^n\mathsf D\xi_i + 2\sum_{i<j}\operatorname{cov}(\xi_i, \xi_j).
  1. cov⁡(aξ+b,cη+d)=accov⁡(ξ,η)\operatorname{cov}(a\xi + b, c\eta + d) = ac\operatorname{cov}(\xi, \eta) и ρ(aξ+b,cη+d)=ρ(ξ,η)\rho(a\xi + b, c\eta + d) = \rho(\xi, \eta) при ac>0ac > 0.
  2. Пусть случайный вектор ξ=(ξ1,…,ξk)\boldsymbol\xi = (\xi_1, \dots, \xi_k) имеет ковариационную матрицу Σ=∥cov⁡(ξi,ξj)∥k×k\boldsymbol\Sigma = \|\operatorname{cov}(\xi_i, \xi_j)\|_{k\times k}, B\boldsymbol B — произвольная числовая (m×k)(m \times k)-матрица (см. П10). Тогда η=Bξ\boldsymbol\eta = \boldsymbol B\boldsymbol\xi имеет ковариационную матрицу ∥cov⁡(ηi,ηj)∥m×m=BΣBT\|\operatorname{cov}(\eta_i, \eta_j)\|_{m\times m} = \boldsymbol B\boldsymbol\Sigma\boldsymbol B^T.

По кирпичикам — свойство 5: ηi=∑kbikξk\eta_i = \sum_kb_{ik}\xi_k, поэтому cov⁡(ηi,ηj)=∑k∑lbikbjlcov⁡(ξk,ξl)\operatorname{cov}(\eta_i, \eta_j) = \sum_k\sum_lb_{ik}b_{jl}\operatorname{cov}(\xi_k, \xi_l) — ровно элемент (i,j)(i, j) произведения BΣBT\boldsymbol B\boldsymbol\Sigma\boldsymbol B^T. Свойство 3 — частный случай с B=(1,1,…,1)\boldsymbol B = (1, 1, \dots, 1): сумма всех элементов Σ\boldsymbol\Sigma.

Смысл по-пацански: разброс суммы — это не только разбросы слагаемых, но и то, как они «тянут в одну сторону». Если два пацана всегда опаздывают вместе (ковариация плюс), общий разброс опозданий бригады больше суммы разбросов. Если один опаздывает, когда другой приходит рано, — разбросы гасят друг друга.

Проверка

Пусть Dξ=3\mathsf D\xi = 3, Dη=6\mathsf D\eta = 6, cov⁡(ξ,η)=1\operatorname{cov}(\xi, \eta) = 1. Найди D(2ξ−η)\mathsf D(2\xi - \eta).

Загрузка редактора…

Проверь себя

Тест

Чем интеграл Лебега отличается от интеграла Римана по идее построения?

Тест

Какое из утверждений верно для ЛЮБЫХ (возможно, зависимых) случайных величин с конечными средними?

Тест

Почему у закона Коши нет математического ожидания?
Обновлено