Файл: Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 216

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

• Исходя из определения математического ожидания и дисперсии случайной величины [см. формулы (2.9), (2.10), (2.16)], естественно в качестве точечных оценок их по результатам п наблюдений при­ нять следующие величины:

эмпирическое, или выборочное, среднее

-1 "

эмпирическая,

или

выборочная,

дисперсия

 

 

(4.6)

 

 

 

 

 

 

 

1

"

(Xi

xf

=

1

"

—2

 

 

 

 

 

а2 = —

2

 

 

V xj —х ,

 

где а

— выборочное среднее квадратическое

отклонение.

 

Очевидно, что все Xi независимы между собой. Каждое отдель­

ное значение Xi можно рассматривать

как самостоятельную случай­

ную

величину,

распределенную

 

по

тому

же

закону,

что и X

с M(xt)

= М(Х),

D(xt)

=

D(X).

Легко

показать, что эмпирическое

среднее

х

несмещенная

и

эффективная

оценка М(Х).

Кроме

того,

поскольку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D(x)=M\x—M(x)f

 

 

={l/n2)

 

2

M[xt

— M(X)]* =

 

 

 

 

 

=

D(X)/n

или

a(x)=-a(X)/i

n,

(4.7)

т . е.

при n

оо дисперсия

D(x) ->• 0,

x — состоятельная

оценка.

Это обстоятельство имеет большое значение для эксперимен­ таторов. Оно показывает, во-первых, что, многократно повторяя даже очень «грубый» опыт, можно получить результат со сколь угодно малой случайной погрешностью, а во-вторых, что средняя квадратическая погрешность среднего х в У п раз меньше погреш­ ности отдельного измерениях^. Правда, это справедливо при усло­ вии отсутствия систематической погрешности в опыте..

Нетрудно показать,

что

 

 

 

 

 

 

 

 

 

M(o*) =

п

^D(X)^D(X),

 

 

 

 

 

 

 

 

 

 

 

 

 

следовательно а2 — смещенная оценка

дисперсии D(X).

Смещение

при и-»- оо исчезает: M(a2 )—v D(X). Оценка

 

 

 

S? =

- 5 - a 2

=

•—

V ( X ; - X ) 2

=

2*?

пх

(4.8)

1

. л — 1

л 1

 

1

 

Я 1

 

 

 

•будет несмещенной. Это следует помнить исследователю, опреде­ ляющему среднее квадратическое отклонение по результатам не-


большого числа наблюдений п. В этих случаях надо пользоваться формулой

a ^ l / ^ l ^ - ^ -

(4-9)

Вычисление о по формуле (4.6) приводит к систематической погреш­ ности

 

e = ( l - | /

i L z i j . 1 0 0 % ,

(4.10)

при п > 50 є < 1 %.

 

 

Можно

показать (см. работу

[29]), что обе оценки дисперсии —

и о2 и о]

— являются состоятельными. Любопытно,

что эффек­

тивной оценкой дисперсии является смещенная оценка о2 , а не о\.

Точечная

оценка коэффициентов

корреляции и регрессии. Для

определения точечной оценки р истинного коэффициента

корреля­

ции р величин X и Y необходимо предварительно на опыте в резуль­

тате наблюдения за X и Y получить п пар значений хи

уг; х.2, у.г\

Тогда

 

 

 

 

 

 

 

 

 

 

 

п

 

 

1

"

 

 

 

 

 

^І(ХІ

— Х){УІ—У)

 

^ХіУі—Ху

 

 

Р=

'= '

-

=

 

 

 

,

(4.11)

 

 

 

пах

a,j

 

ох

o,j

 

 

где х, у, ах

и Сту определяются по

формулам

(4.6).

Величину

cov(XY)——

V (xtХ)(УІу)

называют

эмпирической

или вы-

б о р о ч н о й

к о в а р и а ц и е й [см. формулу

(2.25)].

 

На основе формул (2.25), (4.6) и (4.11) нетрудно получить точеч­

ную оценку r(ylx)

для коэффициента регрессии Y

по X

 

 

 

 

 

 

1

п

 

 

 

 

 

= р ^

=

 

2

ХІУІ—ХУ

 

 

Г(У/Х)

= - ^ 1 = Ц

.

(4.12)

Эту величину

называют

эмпирическим или

в ы б о р о ч н ы м

к о э ф ф и ц и е н т о м р е г р е с с и и Y по X .

 

Параметризация. Обычно интересные для инженерной

практики

законы распределения (см. гл. 3) зависят от одного, максимум двух параметров. Если аналитическая форма закона известна, то для определения его конкретного вида достаточно найти величины этих параметров. Поэтому важной практической задачей является получение оценок параметров закона распределения по результа­ там наблюдений — задача параметризации.


Метод максимального правдоподобия. Возникает естественный вопрос: не существует ли какого-то общего метода нахождения оценок параметров законов распределения, удовлетворяющих свойствам (4.3) — (4.5). В статистике разработано несколько та­ ких методов [7, 20, 30, 31]. Наиболее интересным из них, дающим эффективную и (для практических случаев) состоятельную оценку параметра, является метод максимального правдоподобия. В ос­ нове его лежит понятие функции правдоподобия. Чтобы записать эту функцию, необходимо знать вид (аналитическую форму) закона распределения случайной величины, параметр которой оценивается.

Если рассматриваемая случайная величина дискретна и распре­ делена по закону Р(т, а), то функция правдоподобия записывается в виде

 

 

 

п

 

 

 

 

 

 

L = П Р

и

а),

(4.13)

где

а — параметр

распределения

[см.,

например, формулу

(3.66)];

тх,

т 2 )

т,-,

тп — вариационный ряд значений случайной

величины т, полученных в результате наблюдений над ней. Посколь­

ку mt — независимые случайные величины,

а P(mt, а) — вероят­

ности их появления, то в соответствии

с теоремой умножения (1.8)

величина L есть вероятность, что в результате опыта получим имен­

но заданную совокупность значений mt,

і =

1, 2, .... п.

По аналогии для непрерывных случайных величин функция прав­

доподобия записывается

в виде

 

 

L

= ГЇ f(xh

а),

(4.14)

где а — параметр закона (для иллюстрации выбран однопараметрический закон); хх, х2, xt, хп — вариационный ряд (4.1) для случайной величины X, распределенной по закону f(x, а).

Сущность метода наибольшего правдоподобия заключается в том, что в качестве оценки параметра а принимается такое его зна­ чение а, при котором функции L (или ее логарифм, что одно и то же) максимальна. Иными словами, а определяется из уравнения

dL/da =

0 или д In LI да = 0.

(4.15)

В основе этого метода лежит следующая идея: поскольку в ре­

зультате опыта получена

конкретная совокупность хи

х2, .., хп,

а никакая другая, то вероятность такого события, равная L , долж­ на быть максимальна (ибо это событие произошло, а другое нет, т. е. другое менее вероятно).

Метод моментов. Есть такие законы распределения, для которых вычисление точечных оценок параметров по методу наибольшего правдоподобия приводит к большим вычислительным трудностям



пли вообще невозможно. В этом случае приходится прибегать к дру­

гим методам. Наиболее простым из них является

метод моментов.

Т е о р е т и ч е с к и м

н а ч а л ь н ы м м о м е н т о м пг-го

п о р я д к а случайной

величины

X называется

математическое

ожидание ее т-й степени

 

 

 

v m =

М

(X"') =

' J л-»' / (л-) dx,

. (4.16)

 

 

— со

 

а т е о р е т и ч е с к и м ц е н т р а л ь н ы м м о м е н т о м т-

го п о р я д к а

— математическое

ожидание

 

 

 

 

со

 

 

рт=М[Х-М(Х)]т=

 

^[x-M(X)]"!f(x)dx.

(4.17)

 

 

— со

 

 

Очевидно, что М(Х) = \\,

a D(X)

=

ft.,.

 

Метод моментов заключается в том, что теоретические

моменты

приравниваются

к эмпирическим:

 

 

 

Приравнивается

столько

моментов,

сколько параметров

ищется.

В результате получается система уравнений относительно неизвест­ ных параметров. Эмпирические пли выборочные моменты выража­

ются

соответственно

в

виде

 

 

 

 

 

vm -=

-

V

х?;

f l m = - 2

І-ХГ.

(4.19>

§ 4.2. Нахождение интервальных оценок

 

для характеристик случайных

величин

 

 

Доверительные

интервалы

и вероятности. Если

а — истинное

значение неизвестного параметра, а — его точечная

оценка — слу­

чайная величина

(4.2),

то д о в е р и т е л ь н ы м

и н т е р в а-

л о м

называют отрезок

( а н ш к ш

а верхн)>

построенный около точ­

ки а (рис.- 15), который с заданной вероятностью а накрывает не­

известный

параметр а; причем, н и ж н я я

и в е р х н я я

д о ­

в е р и т е л ь н ы е

г р а н и ц ы

а и и ж 1 1

и а в е Р х П

— это

функции

только результатов

наблюдений хг,

х.2, ..,

хп

и известных

постоян­

ных величин, независящие от оцениваемого параметра

а.

Ины­

ми

словами,

а н ш к н

и

а п е

р х н

являются случайными величинами.

Вероятность

а называется

д о в е р и т е л ь н о й

в е р о я т н о ­

с т ь ю :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

{йнижн <

а <

а в е р х н } =

«•

 

 

(4.20)

Она

характеризует

степень доверия

к утверждению. а п н т н

^

а <

<

я в е р х н .

Очевидно,

что с вероятностью

(1 — а)

истинное

значе­

ние

параметра лежит

вне доверительного

интервала. Эту вероят­

но