Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Выберем начало координат в точке точного значения математического ожидания случайной величины. В этом случае:
|
|
|
o o |
|
n |
|
|
||
|
X |
|
|||||||
|
|
|
o2 |
|
i=1 M xi2 |
|
|
||
|
M xi2 |
= M xi = Dx, |
|
= nDx; |
(10) |
||||
o |
M [xixj] =oM hxixji = Kxixj = 0, |
|
|
||||||
где xi = xi − Mx, xj |
= xj − Mx |
— центрированные значения |
|||||||
случайной величины X. |
|
|
|
|
|
||||
С учетом полученных соотношений имеем |
|
|
|||||||
|
|
|
M [Dx] = |
n −n |
1 |
Dx. |
|
(11) |
|
Из формулы (11) следует, что оценка дисперсии случайной величины, рассчитываемая по уравнению (5), не будет несмещенной. Чтобы ликвидировать получаемое смещение, достаточно умно-
жить величину Dx в формуле (5) на множитель n −n 1. Таким
образом, окончательный вариант формулы для расчета оценки дисперсии случайной величины X, отвечающей требованию состоятельности и несмещенности, имеет вид
Xn
− ˉ 2
(xi Mx)
ˉ |
i=1 |
|
|
(12) |
Dx = |
|
|
. |
|
|
|
n − 1 |
|
|
Полученное выражение (12) для нормального закона распределения будет асимптотически эффективным. При увеличении числа опытов n отношение дисперсии оценки к минимально возможной неограниченно приближается к единице.
Аналогично может быть записана формула для оценки дисперсии случайной величины Z:
|
n |
|
|
|
|
ˉ |
2 |
|
|
|
(zi − Mz) |
|
|
|
|
i=1 |
|
|
|
Dˉz = |
X |
|
. |
(13) |
n − 1 |
|
|||
|
|
ˉ |
системы |
|
|
|
|
||
Формула для оценки корреляционного момента Kxz |
||||
двух случайных величин X и Z имеет вид
6
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
|
n |
|
|
|
|
ˉ |
ˉ |
|
|
|
(xi − Mx)(zi − Mz) |
|
||
|
i=1 |
|
|
|
Kˉxz = |
X |
|
. |
(14) |
n − 1 |
|
|||
|
|
|
|
|
При обработке результатов стрельбы определяют положение главных осей рассеивания, для которых корреляционный момент
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˉ |
ˉ |
|
равен нулю. В этом случае по оценкам дисперсий Dx, Dz и кор- |
||||||||||||||||||
реляционного момента |
ˉ |
|
|
|
находят угол поворота главных осей |
|||||||||||||
Kxz |
|
|||||||||||||||||
рассеивания α. Получим формулу для угла α. |
|
|
||||||||||||||||
Введем новую систему координат 0ξη. В этой системе случай- |
||||||||||||||||||
|
|
|
|
|
o |
|
o |
|
|
|
|
|
|
|
|
|
|
|
ная точка с координатами x, z будет иметь координаты |
|
|
||||||||||||||||
|
|
ξ |
|
o |
sin α |
|
o |
|
|
|
α, |
|
(15) |
|||||
|
|
= x |
+ z cos |
|
||||||||||||||
|
|
η = |
o |
cos α − |
o |
sin |
α. |
|
(16) |
|||||||||
|
|
x |
z |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Получим формулу для корреляционного момента случайных |
||||||||||||||||||
величин ξ и η, принимая Mξ = M η = 0: |
|
|
|
|||||||||||||||
Kξη = M[(ξ − Mξ)(η − Mη)] = M[ |
o2 |
] sin α cos α+ |
|
|
||||||||||||||
x |
|
|
||||||||||||||||
|
|
|
||||||||||||||||
|
o o |
|
2 |
|
|
|
|
|
2 |
|
|
|
|
o2 |
|
|
||
|
xz |
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|||
+ M[ |
](cos |
|
α − sin α) − M[ |
] sin α cos α = |
|
|||||||||||||
|
|
|
|
|
||||||||||||||
|
|
|
|
|
1 |
|
|
|
α(Dz − Dx) + Kxz cos 2α. |
(17) |
||||||||
|
|
|
= − |
|
sin 2 |
|||||||||||||
|
|
|
2 |
|||||||||||||||
Согласно свойству главных осей рассеивания |
|
|
||||||||||||||||
|
|
|
|
|
|
|
Kξη = 0. |
|
|
|
|
|
(18) |
|||||
Из выражений (17) и (18) получим |
|
|
|
|
|
|
||||||||||||
|
|
|
|
tg 2α = |
|
2Kxz |
. |
|
(19) |
|||||||||
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
Dz − Dx |
|
|
|
|||||
С помощью теоремы о дисперсии линейных функций найдем оценки дисперсий для координат ξ и η. Из формул (15) и (16)
находим |
|
|
|
|
|
|
|
|
Dξ = Dx sin2 α + Dz cos2 α + 2Kxz sin α cos α, |
(20) |
|||||||
Dη = Dx cos2 α + Dz sin2 α − 2Kxz sin α cos α. |
(21) |
|||||||
Переходя к оценкам величин, можем записать: |
|
|
||||||
ˉ |
ˉ |
2 |
ˉ |
2 |
ˉ |
α, |
(22) |
|
Dξ = Dx sin |
|
|
α + Dz cos |
|
α + Kxz sin 2 |
|||
ˉ |
ˉ |
|
2 |
ˉ |
2 |
ˉ |
α. |
(23) |
Dη = Dx cos |
|
α + Dz sin |
|
α − Kxz sin 2 |
||||
|
|
|
|
|
|
|
|
7 |
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Значения оценок получены для ограниченного числа опытов. Для определения точности и надежности результатов рассчитывают доверительные интервалы и доверительные вероятности.
В качестве примера рассмотрим понятие доверительного интервала для некоторой величины a, несмещенная оценка которой равна aˉ. Выберем достаточно большую вероятность β (0,8 или 0,95). Будем считать, что событие с вероятностью β является практически достоверным.
Найдем величину ε, для которой вероятность расхождения между точным значением a и ее оценкой a,ˉ меньшей выбранной величины ε, будет равна β:
p(|aˉ − a| < ε) = β. |
(24) |
В этом случае ошибки в определении a, которые по абсолютной величине больше ε, будут появляться с очень малой вероятностью (1 − β).
Представим выражение (24) |
в ином виде: |
|
p(ˉa − ε < a |
< aˉ + ε) = β. |
(25) |
Как следует из уравнения (25), неизвестное значение величины a с вероятностью β попадает в интервал
Iβ = (ˉa − ε; aˉ + ε). |
(26) |
Другими словами, интервал Iβ с вероятностью β накрывает истинное значение величины a. В этом случае вероятность β называется доверительной вероятностью, а интервал Iβ — доверительным интервалом.
Как правило, нахождение границ доверительного интервала бывает затруднено из-за незнания закона распределения случайной величины, числовую характеристику a которой мы определяем.
Рассмотрим способ определения доверительного интервала для математического ожидания случайной величины Х .
Как следует из центральной предельной теоремы [2], закон распределения оценки математического ожидания независимой слу-
чайной величины |
ˉ |
при достаточно большом числе опытов n |
||
Mx |
||||
близок к нормальному. |
|
|
|
|
Найдем величину εМ , для которой справедливо выражение |
||||
|
p( Mˉx |
− M < εM ) = β. |
(27) |
|
8 |
|
|
|
|
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Выразим вероятность β через функцию Лапласа Φ(x):
r
β = Φ |
εM |
, |
(28) |
σM √2 |
где σ Dx — среднее квадратичное отклонение оценки ˉ .
M = Mx
n
Теоретическое значение дисперсии Dx неизвестно, поэтому при оценке границ доверительного интервала можно воспользо-
ˉ |
|
|||||||
ваться ее оценкой Dx, тогда |
|
|
|
|||||
σM = r |
ˉ |
|
|
|||||
D |
|
|||||||
x |
. |
(29) |
||||||
n |
||||||||
Используя уравнение (28), находим |
|
|||||||
εM = σM √ |
|
Φ−1(β). |
(30) |
|||||
2 |
||||||||
Введем обозначение |
|
|||||||
tβ = √ |
|
Φ−1(β). |
(31) |
|||||
2 |
||||||||
В приложении приведена табл. П1 зависимости параметра tβ от доверительной вероятности β.
Таким образом, окончательно формула оценки границ доверительного интервала для математического ожидания случайной
величины Х имеет вид |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
εMx = tβr |
ˉ |
|
|
|
|
|
|
D |
|
|
|
|
|
||
x |
, |
|
|
(32) |
|||
n |
|
|
|||||
|
ˉ |
− εMx |
; |
ˉ |
(33) |
||
(Iβ)Mx = (Mx |
Mx + εMx ). |
||||||
Рассмотрим схему определения доверительного интервала для дисперсии случайной величины Х с известными оценками мате-
ˉ |
ˉ |
|
|
матического ожидания Mx и дисперсии |
Dx. |
|
|
Как ясно из формулы (12), оценка дисперсии представляет со- |
|||
|
ˉ |
2 |
|
бой сумму n случайных величин типа |
(xi − Mx) |
|
, которые не |
|
n |
|
|
являются независимыми вследствие того, что−в1 каждую из них |
|||
входит оценка математического ожидания. В соответствии с центральной предельной теоремой с увеличением n закон распределения их суммы приближается к нормальному.
9
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Характеристики этого закона будут определяться следующим
образом. Математическое ожидание дисперсии ˉ :
Dx
|
|
M ˉ = Dx. |
|
|
|
||||
|
|
|
|
|
D |
|
|
|
|
ˉ |
: |
|
|
|
|
|
|
|
|
Дисперсия D ˉ |
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
D ˉ = |
1 |
|
μ |
n − 3 |
D2 |
, |
||
|
n |
||||||||
|
D |
|
4 − n(n |
− |
1) |
x |
|
||
|
|
|
|
|
|
|
|
|
|
(34)
(35)
где μ4 — четвертый центральный момент случайной величины Х . Вместо теоретического значения дисперсии Dx в уравнении
(35) можно пользоваться ее оценкой ˉ .
Dx
Серьезные трудности возникают при определении величины μ4. Ее можно было бы определить, используя результаты обработки опытных данных. Однако для ограниченного числа испытаний такая оценка дает слишком большие погрешности. Воспользуемся формулой связи четвертого центрального момента случайной величины μ4 с дисперсией Dx для нормального закона распределения случайной величины Х :
μ4 = 3Dx2. |
(36) |
Подставляя формулу (36) в уравнение (35) и заменяя значение дисперсии ее оценкой, получаем
|
2 |
ˉ2 |
|
|
|
|
|
D ˉ = |
|
D |
, |
D |
n − 1 x |
|
|
r
2
σ ˉ
ˉ = Dx.
D n − 1
(37)
(38)
Оценка границ доверительного интервала для дисперсии случайной величины Х может быть получена с помощью следующих формул:
εDx = tβr |
|
|
|
|
|
|
|
|
2 |
|
Dˉx, |
||
|
n |
− |
1 |
|||
ˉ |
|
|
ˉ |
|||
(Iβ)Dx = (Dx − εDx ; Dx + εDx ).
(39)
(40)
Как известно, наиболее полную информацию о случайной величине дает закон распределения. Поэтому одной из задач при проведении испытаний является установление этого закона, для
10