1

Многомерные распределения

[27/63%]
Показать
LaTeX
Задача 1.1

Плотность точки попадания (X,Y)(X, Y) в круглую мишень задаётся формулой

f(x,y)=c(R−x2+y2)дляx2+y2≤R2. f(x, y) = c(R - \sqrt{x^{2} + y^{2}}) \qquad \text{для} \quad x^{2} + y^{2} \leq R^{2}.

Найдите, (а) константу cc, (б) вероятность pp того, что точка попадания окажется в круге радиуса a<Ra < R с центром в начале координат.

?
(a)
(b)
Задача 1.2

Случайные величины X,YX, Y удовлетворяют линейному соотношению

aX−bY=c, aX - bY = c,

где a,b,ca, b, c — константы. Найдите, (а) корреляцию ρ(X,Y)\rho (X, Y) и (б) отношение σx/σy\sigma_{x}/\sigma_{y} стандартных отклонений.

?
(a)
(b)
Задача 1.3

Матрица рассеяния (дисперсий--ковариаций) трёхмерного нормального распределения равна

(σij)=[52−2263−238]. (\sigma _{ij}) = \begin{bmatrix} 5 & 2 & -2 \\ 2 & 6 & 3 \\ -2 & 3 & 8 \end{bmatrix}.

Если μx=μy=μz=0\mu_{x} = \mu_{y} = \mu_{z} = 0, найдите плотность f(x,y,z)f(x, y, z) и её максимум.

?
Задача 1.4

Совместная плотность XX и YY задаётся формулой

f(x,y)=ye−y(x+1),x>0,y>0. f(x, y) = ye^{-y(x+1)}, \qquad x > 0, \quad y > 0.

Найдите, (а) маргинальные плотности XX и YY, и (б) условную функцию распределения FX(x∣y)F_{X}(x \mid y) величины XX при условии Y=yY = y.

?
(a)
(b)
Задача 1.5

Предположим, что случайная величина XX совпадает со случайной величиной XiX_{i} с вероятностью pip_{i}, то есть,

P(X=Xi)=piиE[Xi]=mi,i=1,2,…. \mathbb {P}\left(X = X_{i}\right) = p_{i} \qquad \text{и} \qquad \mathbb {E}\left[X_{i}\right] = m_{i}, \qquad i = 1, 2, \ldots .

Покажите, что

Var⁡[X]=∑ipiVar⁡[Xi]+Var⁡[M], \operatorname {Var}\left[X\right] = \sum _{i} p_{i} \operatorname {Var}\left[X_{i}\right] + \operatorname {Var}\left[M\right],

где P(M=mi)=pi\mathbb {P}\left(M = m_{i}\right) = p_{i} (i=1,2,…)(i = 1, 2, \ldots ).

?
Задача 1.6

Случайный вектор (X,Y,Z)(X, Y, Z) имеет плотность f(x,y,z)f(x, y, z). X,Y,ZX, Y, Z можно наблюдать только одновременно. Наблюдение дало значения u,v,wu, v, w без знания соответствия между ними и X,Y,ZX, Y, Z. Найдите вероятность того, что X=uX = u, Y=vY = v и Z=wZ = w.

?
Задача 1.7

Совместная функция распределения F(x,y)F(x, y) случайных величин XX и YY задаётся формулой

F(x,y)=sin⁡xsin⁡y,0≤x≤π/2,0≤y≤π/2. F(x, y) = \sin x \sin y, \qquad 0 \leq x \leq \pi /2, \quad 0 \leq y \leq \pi /2.

Найдите, (а) функцию плотности f(x,y)f(x, y), и (б) матрицу рассеяния.

?
(a)
(b)
Задача 1.8

Нормальная плотность точки A=(X1,X2,X3)A = (X_{1}, X_{2}, X_{3}) задаётся формулой

f(x,y,z)=1(2π)3/2σ1σ2σ3exp⁡{−12(x2σ12+y2σ22+z2σ32)}. f(x, y, z) = \frac{1}{(2\pi )^{3/2}\sigma _{1}\sigma _{2}\sigma _{3}} \exp \left\{ -\frac{1}{2}\left(\frac{x^{2}}{\sigma _{1}^{2}} + \frac{y^{2}}{\sigma _{2}^{2}} + \frac{z^{2}}{\sigma _{3}^{2}}\right)\right\} .

Какова вероятность того, что AA попадёт в эллипсоид с главными полуосями E1,E2,E3E_{1}, E_{2}, E_{3} вдоль осей Ox,Oy,OzOx, Oy, Oz, где E1,E2,E3E_{1}, E_{2}, E_{3} — вероятные ошибки X,Y,ZX, Y, Z соответственно, т.е.,

P(∣X∣<E1)=P(∣Y∣<E2)=P(∣Z∣<E3)=1/2. \mathbb {P}\left(\left|X\right| < E_{1}\right) = \mathbb {P}\left(\left|Y\right| < E_{2}\right) = \mathbb {P}\left(\left|Z\right| < E_{3}\right) = 1/2.
?
Задача 1.9

Вектор (X,Y)(X, Y) имеет равномерное распределение в эллипсе

x2a2+y2b2≤1. \frac{x^{2}}{a^{2}} + \frac{y^{2}}{b^{2}} \leq 1.

Найдите:

?
(a)

маргинальные распределения XX и YY;

(b)

условные распределения XX при условии Y=yY = y, и условные распределения YY при условии X=xX = x;

(c)

Cov⁡[X,Y]=σxy\operatorname {Cov}\left[ X, Y \right] = \sigma_{xy}.

Независимы ли X,YX, Y?

Задача 1.10

Точка (X,Y,Z)(X, Y, Z) равномерно распределена в шаре x2+y2+z2≤R2x^{2} + y^{2} + z^{2} \leq R^{2}. Найдите плотность ZZ и условную плотность f(x,y∣z)f(x, y \mid z). Обобщите на nn-мерный шар.

?
Задача 1.11

Даны fY(y)f_{Y}(y), E[X∣y]\mathbb {E}\left[ X \mid y \right], Var⁡[X∣y]\operatorname {Var}\left[ X \mid y \right]; найдите E[X]\mathbb {E}\left[X\right] и Var⁡[X]\operatorname {Var}\left[X\right].

?
Задача 1.12

Вектор (X,Y)(X, Y) имеет нормальную плотность

f(x,y)=ce−4x2−6xy−9y2. f(x, y) = ce^{-4x^{2} - 6xy - 9y^{2}}.

Найдите плотности fX(x)f_{X}(x), fY(y)f_{Y}(y), fX(x∣y)f_{X}(x \mid y), fY(y∣x)f_{Y}(y \mid x) и константу cc.

?
Задача 1.13

Предположим, что X→\overrightarrow {X} имеет сферически симметричное nn-мерное распределение, т.е. с плотностью f(x12+⋯+xn2)f(x_{1}^{2} + \cdots + x_{n}^{2}). Докажите, что плотность gg «обобщённой χ2\chi^{2}-величины» U=X→′X→U = \overrightarrow {X}'\overrightarrow {X} задаётся формулой

g(u)=12Cnf(u)un/2−1, g(u) = \tfrac {1}{2}C_{n}f(u)u^{n/2-1},

где CnC_{n} — «поверхностная» площадь nn-мерной сферы, т.е.

Cn=2πn/2Γ(n/2). C_{n} = \frac{2\pi ^{n/2}}{\Gamma (n/2)}.
?
Задача 1.14

(Продолжение задачи 1.13.) Для n=2n = 2 покажите, что полярные координаты (R,θ)(R, \theta ), где

X=Rcos⁡θ,Y=Rsin⁡θ, X = R\cos \theta , \qquad Y = R\sin \theta ,

независимы, θ\theta равномерно распределена на [0,2π)[0, 2\pi ), а RR имеет плотность (ср. задачу 221)

h(r)=2πrf(r2). h(r) = 2\pi r f(r^{2}).
?
Задача 1.15

(Продолжение задачи 1.14.) Покажите, что коэффициент корреляции, основанный на единственном наблюдении (X,Y)(X, Y), т.е.

r0=XYX2+Y2 r_{0} = \frac{XY}{X^{2} + Y^{2}}

имеет плотность

f(r0)=1π11−r02,−1<r0<1. f(r_{0}) = \frac{1}{\pi } \frac{1}{\sqrt{1 - r_{0}^{2}}}, \qquad -1 < r_{0} < 1.
?
Задача 1.16

(Продолжение задачи 1.15.) Если X→=(X1,…,Xn)′\overrightarrow {X} = (X_{1}, \ldots , X_{n})' есть N(0,In)\mathscr {N}\left(0, I_{n}\right), и X→\overrightarrow {X} разбит на kk подвекторов X(1),…,X(k)X_{(1)}, \ldots , X_{(k)} с v1,v2,…,vkv_{1}, v_{2}, \ldots , v_{k} компонентами соответственно, покажите, что (X(1)′X(1),…,X(k)′X(k))/∣X∣2(X_{(1)}'X_{(1)}, \ldots , X_{(k)}'X_{(k)})/\left|X\right|^{2} имеет распределение Дирихле (см. (5.11)) с плотностью f(t1,…,tk−1)f(t_{1}, \ldots , t_{k-1}), приведённой ниже.

Формула (5.11), распределение Дирихле: для x1,…,xkx_{1}, \ldots , x_{k} с ∑i=1k+1xi=1\sum_{i=1}^{k+1} x_{i} = 1, xi≥0x_{i} \geq 0,

f(x1,…,xk)=Γ(n1+⋯+nk+1)Γ(n1)+⋯+Γ(nk+1)∏i=1k+1xini−1. f(x_{1}, \ldots , x_{k}) = \frac{\Gamma (n_{1} + \cdots + n_{k+1})}{\Gamma (n_{1}) + \cdots + \Gamma (n_{k+1})} \prod _{i=1}^{k+1} x_{i}^{n_{i}-1}.

Искомая здесь плотность f(t1,…,tk−1)f(t_{1}, \ldots , t_{k-1}) равна

f(t1,…,tk−1)=Γ(n2)∏i=1kΓ(vk2)(1−t1−⋯−tk−1)(vk/2)−1∏i=1k−1ti(vi/2)−1, f(t_{1}, \ldots , t_{k-1}) = \frac{\Gamma \left(\dfrac {n}{2}\right)}{\displaystyle \prod _{i=1}^{k} \Gamma \left(\dfrac {v_{k}}{2}\right)} (1 - t_{1} - \cdots - t_{k-1})^{(v_{k}/2)-1} \prod _{i=1}^{k-1} t_{i}^{(v_{i}/2)-1}, ti≥0,∑ti<1. t_{i} \geq 0, \quad \sum t_{i} < 1.

Отсюда, при v1=v2=⋯=vs=1v_{1} = v_{2} = \cdots = v_{s} = 1, vs+1=n−sv_{s+1} = n - s (k=s+1)(k = s+1),

(∣X1∣,∣X2∣,…,∣Xk∣)∣X∣ \frac{(\left|X_{1}\right|, \left|X_{2}\right|, \ldots , \left|X_{k}\right|)}{\left|X\right|}

имеет плотность

f0(v1,…,vk)=Γ(n2)2kΓ(n−k2)π(1−v12−⋯−vk2)((n−k)/2)−1,vi≥0,∑i=1kvi2<1, f_{0}(v_{1}, \ldots , v_{k}) = \frac{\Gamma \left(\dfrac {n}{2}\right) 2^{k}}{\Gamma \left(\dfrac {n-k}{2}\right)\pi } (1 - v_{1}^{2} - \cdots - v_{k}^{2})^{((n-k)/2)-1}, \qquad v_{i} \geq 0, \quad \sum _{i=1}^{k} v_{i}^{2} < 1,

а плотность (X1,…,Xk)/∣X∣≡(U1,U2,…,Uk)(X_{1}, \ldots , X_{k})/\left|X\right| \equiv (U_{1}, U_{2}, \ldots , U_{k}) равна

f0(u1,…,uk)2kесли∑i=1kui2<1. \frac{f_{0}(u_{1}, \ldots , u_{k})}{2^{k}} \qquad \text{если} \quad \sum _{i=1}^{k} u_{i}^{2} < 1.

Это распределение любого kk-мерного подвектора nn-мерного вектора

U→=(U1,…,Un) \overrightarrow {U} = (U_{1}, \ldots , U_{n})

равномерно распределённого на (поверхности) единичной сферы

x12+⋯+xn2=1. x_{1}^{2} + \cdots + x_{n}^{2} = 1.
?
Задача 1.17

(Продолжение задачи 1.16.) В общем случае, если X→\overrightarrow {X} имеет сферическое распределение (см. задачу 1.13), поскольку оно инвариантно относительно вращений в nn-мерном пространстве, распределение направления X→/∣X∣\overrightarrow {X}/\left|X\right| совпадает с распределением U→\overrightarrow {U} и не зависит от его длины ∣X∣\left|X\right| (для любого сферического распределения с P(∣X∣=0)=0\mathbb {P}\left(\left|X\right| = 0\right) = 0).

Отсюда заключите, что результаты задачи 1.16 верны, если X→\overrightarrow {X} — произвольный сферически симметричный случайный вектор.

?
Задача 1.18

Пусть X,YX, Y — случайные величины с

E[X]=μ,E[Y]=η,Var⁡[X]=σ12,Var⁡[Y]=σ22,Cov⁡[X,Y]=σ12. \mathbb {E}\left[X\right] = \mu , \qquad \mathbb {E}\left[Y\right] = \eta , \qquad \operatorname {Var}\left[X\right] = \sigma _{1}^{2}, \qquad \operatorname {Var}\left[Y\right] = \sigma _{2}^{2}, \qquad \operatorname {Cov}\left[ X, Y \right] = \sigma _{12}.

Найдите приближённую формулу для

E[g(X,Y)],Var⁡[g(X,Y)], \mathbb {E}\left[g(X, Y)\right], \quad \operatorname {Var}\left[g(X, Y)\right],

используя линейное разложение g(x,y)g(x, y) в окрестности точки (μ,η)(\mu , \eta ) (с точностью до моментов второго порядка). Обобщите результат на nn случайных величин.

?
Задача 1.19

XX и YY — независимые одинаково распределённые случайные величины с плотностью

f(x)=2π1−x2,0≤x≤1. f(x) = \frac{2}{\pi \sqrt{1 - x^{2}}}, \qquad 0 \leq x \leq 1.

Используя задачу 1.18, найдите E[Z]\mathbb {E}\left[Z\right] и Var⁡[Z]\operatorname {Var}\left[Z\right], где

Z=arctan⁡XY. Z = \arctan \frac{X}{Y}.
?
Задача 1.20

Основная частота струны определяется формулой

Ω=12FmL, \Omega = \frac{1}{2}\sqrt{\frac{F}{mL}},

где FF — натяжение, mm — масса струны, а LL — её длина. Если масса считается постоянной, а F,LF, L — случайные величины с E[F]=f\mathbb {E}\left[F\right] = f, E[L]=l\mathbb {E}\left[L\right] = l, Var⁡[F]=σf2\operatorname {Var}\left[F\right] = \sigma_{f}^{2}, Var⁡[L]=σl2\operatorname {Var}\left[L\right] = \sigma_{l}^{2} и Cov⁡[F,L]=σfl\operatorname {Cov}\left[ F, L \right] = \sigma_{fl}, найдите приближение для дисперсии Ω\Omega.

?
Задача 1.21

Пусть вектор (X,Y)(X, Y) имеет плотность f(x,y)f(x, y). Определим комплексные переменные

Z=X+iY,Zt=Zeit, Z = X + iY, \qquad Z_{t} = Ze^{it},

и Zt=ZeitZ_{t} = Ze^{it}. Докажите, что для того чтобы все ZtZ_{t} имели одинаковое распределение, необходимо, чтобы f(x,y)=g(x2+y2)f(x, y) = g(x^{2} + y^{2}) для некоторой плотности gg (см. задачу 1.13).

?
Задача 1.22

Если X→\overrightarrow {X} имеет плотность ff, как в задаче 1.13, покажите, что

Y→=μ→+AX→(A — невырожденная постоянная матрица размера n×n) \overrightarrow {Y} = \overrightarrow {\mu } + A\overrightarrow {X} \qquad \text{(}A\text{ — невырожденная постоянная матрица размера }n \times n\text{)}

имеет плотность вида

c∣Λ∣−1/2f((y→−μ→)′Λ−1(y→−μ→)),Λ=AA′, c\left|\Lambda \right|^{-1/2}f((\overrightarrow {y} - \overrightarrow {\mu })'\Lambda ^{-1}(\overrightarrow {y} - \overrightarrow {\mu })), \qquad \Lambda = AA',

и характеристическую функцию вида

eit′μψ(t′Λt)для некоторой функции ψ. e^{it'\mu }\psi (t'\Lambda t) \qquad \text{для некоторой функции } \psi .
?
Задача 1.23

Пусть ρij=ρ(Xi,Xj)\rho_{ij} = \rho (X_{i}, X_{j}) — коэффициент корреляции между Xi,XjX_{i}, X_{j} для i,j=1,…,ni, j = 1, \ldots , n. Для n=3n = 3 и ρ12=ρ13=ρ23=ρ\rho_{12} = \rho_{13} = \rho_{23} = \rho какие значения может принимать ρ\rho?

?
Задача 1.24

Если pp-мерный вектор X→\overrightarrow {X} имеет распределение N(μ,Λ)\mathscr {N}\left(\mu , \Lambda \right), покажите, что квадратичная форма

χ2=(X→−μ→)′Λ−1(X→−μ→) \chi ^{2} = (\overrightarrow {X} - \overrightarrow {\mu })'\Lambda ^{-1}(\overrightarrow {X} - \overrightarrow {\mu })

имеет распределение χp2\chi_{p}^{2} с pp степенями свободы и плотностью, задаваемой (ср. задача 1.13)

?
Задача 1.25

Энтропия HH многомерного распределения с плотностью f(x→)f(\overrightarrow {x}) определяется формулой

H(f)=−E[log⁡f(X→)]=−∫⋯∫f(x→)log⁡f(x→) dx→. H(f) = -\mathbb {E}\left[\log f(\overrightarrow {X})\right] = -\int \cdots \int f(\overrightarrow {x}) \log f(\overrightarrow {x}) \, d\overrightarrow {x}.

Докажите, что для нормального распределения из задачи 1.24

H=log⁡[(2πe)p/2∣Λ∣−1/2]. H = \log \left[(2\pi e)^{p/2}\left|\Lambda \right|^{-1/2}\right].
?
Задача 1.26

Если X→\overrightarrow {X} имеет распределение N(μ,Λ)\mathscr {N}\left(\mu , \Lambda \right) с плотностью ff, найдите плотность gg случайной величины

Y=f(X→). Y = f(\overrightarrow {X}).
?
Задача 1.27

Если X→=(X1,…,Xn)\overrightarrow {X} = (X_{1}, \ldots , X_{n}) равномерно распределён в симплексе TT:

x1+x2+⋯+xn≤1,xi≥0,i=1,…,n, x_{1} + x_{2} + \cdots + x_{n} \leq 1, \qquad x_{i} \geq 0, \qquad i = 1, \ldots , n,

покажите, что плотность X→\overrightarrow {X} равна

f(x→)=n!дляx→∈T. f(\overrightarrow {x}) = n! \qquad \text{для} \quad \overrightarrow {x} \in T.

Покажите, что (n−1)(n-1)-мерные маргинальные распределения также являются частными случаями распределения Дирихле (5.11).

Формула (5.11), распределение Дирихле: для x1,…,xkx_{1}, \ldots , x_{k} с ∑i=1k+1xi=1\sum_{i=1}^{k+1} x_{i} = 1, xi≥0x_{i} \geq 0,

f(x1,…,xk)=Γ(n1+⋯+nk+1)Γ(n1)+⋯+Γ(nk+1)∏i=1k+1xini−1. f(x_{1}, \ldots , x_{k}) = \frac{\Gamma (n_{1} + \cdots + n_{k+1})}{\Gamma (n_{1}) + \cdots + \Gamma (n_{k+1})} \prod _{i=1}^{k+1} x_{i}^{n_{i}-1}.
?