Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Статистические методы в строительной механике

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
18.02 Mб
Скачать

ttс корреляционной матрицей

ПП

К Y.Yk = 2

2 ajaй/1? ^ x* x?*

0.45)

e=l

p=l

 

Необходимо отметить, что формулы (1.44) и (1.45) остаются справедливыми не только для нормального распределения, но и для любых случайных величин, подвергающихся линейному преобразованию .(1-43). Для обоснования формулы (1.44) до­ статочно произвести осреднение по множеству обеих частей со­ отношения (1.43). Чтобы получить формулу (1.45), составим вы­ ражение

*V.y = ( ^ - ^ ) ( П - П ) .

J к

Учитывая (1.44), получим

Я Ъ - ,= 2 « * < * . - * . )

(Хр- Х ,) ~

а-1 P-»

пп

 

= 2

ak? ^ x*xv

 

 

«=i p=i

 

 

что и требовалось

доказать.

 

 

При нелинейном преобразовании случайных функций

 

Y, = ФУ(Х„ Х„

Х„>

(1.46>

 

(/ =

1. 2,

т)

 

соотношения (1.44)

и .(1-45)

не имеют места, а нормальное рас­

пределение не остается таковым. Однако нередко

функции Ф/

в окрестности математических ожиданий случайных величин Xit Хг, ...Хп мало отличаются от линейных функций, а их дисперсия невелика. Тогда функции Фу можно приближенно представить, в виде линейных функций, раскладывая их в степенные ряды в.

окрестности точки Хи Х2,... Хп и отбрасывая нелинейные члены:

 

 

П

_ _ _

 

Y, æ Ф, (Х„ X *

. X.) + V

№ -

X»).

J

 

k=1

дХ к

 

 

 

 

 

Отсюда

 

 

 

 

 

 

Ÿ j~Ф у(Х 1г Х„

.а д ,

(1.47)

Ь, _

Г 1 <ЭФу (Xl, Х21 . . . Х„)

(Xl, ХЯ, . . . Х„)

/4 ,404

 

-----------J z

Щ

(ь48>

0=1 р=1

Если случайные величины Xi, Хц..Хп некоррелированы, то для дисперсий случайных величин Уь Ÿ2...Ym получаем фор­ мулу

D(Y,) »

Г

 

J D(Xj)

(1.49)

Л-1

*■

ft

J

 

Формулы (1.47), (1.48) и (1.49) широко используются в раз­ нообразных практических приложениях.

12. Распределение крайних значений случайных величин

Допустим, что производится п измерений случайной величи­ ны X, имеющей функцию распределения F (х) и плотность веро­ ятности р(х). Результаты измерений дают последовательность чисел хп1« хп2 < ... « хпп. Необходимо найти функцию распре­ деления Fnn(x) и плотность вероятности

=

0 -50)

для максимальных значений случайной величины хпп в сово­ купности п измерений. Функция распределения Fnn(x) есть не что иное, как вероятность обнаружить в каждом из п измерений неравенство Х<х. Если примять схему независимых испытаний, то задача решается чрезвычайно просто. Вероятность обнару­ жить неравенство Х^<. х в результате одного измерения равна, очевидно, F(x). Отсюда по теореме умножения вероятностей

Л .М = P . & < 4

= F(x).

 

Применение формулы

(1.50)

дает

 

 

 

Р„п(*) =

nFn~l (.х)р(х).

 

(1.51)

Аналогичную задачу можно поставить

для

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

минимальных значений

хп1. Соответствующие

характеристики

обозначим через Fnl(x) и рпл{х). Согласно

определению, функ­

ция 1 —Fnl (л:) равна

вероятности

обнаружить неравенство

Х>х в результате каждого из п измерений. Замечая, что по тео­ реме умножения вероятностей

Р„( Х> х ) = [P(X > х)\* = [1 - F (*)]■,

найдем

Fnl( x ) = \ - [ l - F ( x ) ] » .

Эта формула аналогична формуле (1.6). Дифференцируя функ­ цию распределения, получим

Pni (x) = n [ l — F(x)]nл“ 1р(х)- /А

(1.52)

График для плотностей вероятности (1.51) и (1.52) представлен на рис. 5. С увеличением числа испытаний наиболее вероятные значения максимумов случайной величины смещаются, естествен­ но, вправо, а наиболее вероятные значения минимумов — влево. При этом все распределения становятся более компактными.

Задача о распределении крайних значений представляет весьма большой интерес для вопросов прочности. Объясняется это тем, что с точки зрения расчетов на прочность нас больше всего интересуют максимальные значения нагрузок и мини­

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

Если число п велико, то использование формул (1.51) и (1.52) крайне затруднено. Между, тем этот случай как раз типи­ чен для приложений. Например, прочность образца из хрупко­ го материала характеризуется прочностью наиболее дефектного первичного объема, а число этих первичных объемов в образце весьма велико. Поэтому большой интерес представляет изуче­ ние асимптотических свойств распределений при п-*оо, посколь­ ку их можно использовать для приближенной оценки при конеч­ ных, но достаточно больших п. При довольно широких предпо­ ложениях могут быть доказаны две теоремы, формулировку одной из которых, относящуюся к распределению минимальных

значений,

мы приводим ниже.

удовлетворяет

условиям:

Пусть

функция

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

(a) F(x) = 0

при х< х0, но F(x) > 0 при х>х0; (б) при достаточно

малых е> 0

выполняется условие

 

 

 

 

 

 

lim F(*о 4- е) _

 

 

 

 

 

е-*0

 

 

 

где с и а — некоторые положительные числа. Тогда

для функ­

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

минимальных значений

Fnl(x) при больших

п имеет

место асимптотическое соотношение

 

 

 

 

Fnl (•*) = 0 при

х < х0,

 

 

 

F„i (*) ~

1 — ехр [ — сп{х — *„)“]

при х >

х0.

Здесь а и о — математическое ожидание и стандарт для исход­ ного распределения (1.17). Что касается самого распределения, то оно асимптотически описывается «двойным экспоненциаль­ ным законом». Например,

Ftll.(x) ^ 1 — exp [— exp

13. Распределение Пирсона

Среди других теоретических распределений остановимся на распределении Пирсона, которое нам понадобится в дальней­ шем. Английским статистиком К. Пирсоном было предложено

12 типов распределений, среди которых наибольшее значение приобрел, пожалуй, третий тип, или %2-распределение Пирсона;

 

 

 

 

О

 

при и < О,

(1.57)

 

р(и) =

ип-1

и*

 

 

 

 

 

2

при и >• 0.

 

 

 

 

w (п)

 

 

 

 

 

 

 

 

 

 

 

 

Здесь п — некоторая

положительная

величина,

хУ (л)— норми­

ровочная постоянная.

 

(1.9)

найдем,

что

 

Из условия

нормировки

 

 

W (п) -

f ип~1 е

2

du.

 

 

 

 

 

о

 

 

 

 

 

 

Этот интеграл легко выражается через полную гамма-функ­

цию Г(х):

 

 

 

 

п—2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЧГ(л) =

2 2

Г (-у)-

 

(1.58)

Заметим, что Ч;(1) =

^(З) =

J /

-j,

1Р(2) = 1.

При целых по­

ложительных п

 

 

 

 

 

 

 

 

 

 

 

(п — 2)!!

у

 

.если п — нечетное число,

* < » ) -

* 3 , п - г \

 

 

 

 

 

 

(1'59)

 

2 2

 

— J !,если п — четное число,

где (л -2)!! = 1 • 3.. .(л-4) (л-2).

распределения имеет вид

Соответствующая

функция

 

 

 

 

1

1

Л- 1

и*

,

 

 

 

 

 

С

--- о

 

 

 

 

¥

(и) j

«

 

*

da

 

и выражается через

неполную гамма-функцию.

Вместо F(x,n)

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

употребляют

дополнительную

величину

Р(х\п ) , называемую

функцией

х2*распределения Пирсона:

 

 

 

Оо

и»

 

P (x» . n ) = l - F ( x . " ) =

^

J

2 du.

(1.60)

 

 

 

X

 

 

Эта функция протабулирована [43, 101].

Значение х2-распределения Пирсона для приложений опре­ деляется несколькими причинами. Укажем, что распределение

оказывается весьма удобным для описания случайных величин, которые могут принимать лишь положительные значения. При п= 1 оно превращается в распределение, ординаты которого рав­

ны удвоенным

ординатам положительной ветви гауссовского

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

при п = 2 мы получаем распределение Рэлея, при

я = 3 — важное

для статистической физики распределение Мак­

свелла и т. д.

(рис. 8).

Это распределение имеет важное приложение в математиче­ ской статистике. Можно показать, что ему подчиняется случай­

ная величина

 

г =

(1-6 D

вычисленная по п реализациям Хи

случайной величины

X, которая имеет нормальное распределение. На этом основан критерий согласия Пирсона, позволяющий судить о согласован­ ности опытных данных и теоретической кривой распределения (см. ниже п. 15).

Кроме того, как обнаружится в дальнейшем, к интегралам, выражаемым через функцию Р(х2,п), приводят некоторые весь­ ма интересные задачи по исследованию прочности при случай­ ных нагрузках.

14.Предмет математической статистики. Эмпирические распределения и их важнейшие характеристики

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

Частота

Рис. 9

обработки и истолкования характеристик эмпирического распре­ деления. Основная задача статистики состоит в следующем: по известному распределению эмпирических частот найти соответ­ ствующее ему теоретическое распределение вероятностей. В об­ щем виде эта задача весьма трудна. Однако во многих случаях вид функции распределения можно предсказать, исходя из об­ щетеоретических соображений. Тогда задача сводится к при­ ближенному отысканию неизвестных параметров распределения и к оценке надежности найденных значений. Так, если заранее можно сказать, что исследуемое распределение случайной ве­ личины является нормальным, то задача состоит в отыскании математического ожидания и дисперсии, точнее, в установлении доверительных границ, в которых эти параметры лежат с напе­ ред заданной вероятностью.'

'Допустим, что в результате п измерений случайной величины X получена последовательность значений (вариационный ряд)

Хи Х2,...хп. Первичная обработка вариационного ряда состоит в группировке найденных значений по достаточно малым интер­ валам, вычислении средних относительных частот для каждого интервала и графическом представлении результатов в виде ги­ стограммы (рис. 9), полигона (рис. 10) или кумулятивной кри­ вой (рис. 11). Очевидно, что гистограмма и полигон соответству­ ют кривой теоретической плотности вероятности, а кумулятив­

ная кривая — кривой

теоретической

функции

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

Употребительны также термины

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

кривая

плот­

ности вероятности» и «эмпирическая кривая распределения».

Важнейшими характеристиками

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

распределе­

ния будут среднее арифметическое значение

 

 

 

 

 

 

Х =

 

*k

 

 

 

 

 

(1.62)

и среднее квадратическое отклонение

 

 

 

 

 

 

\ X - X f = ^ ~

У < * * -* >

 

 

 

 

(1.63)

 

п

k - i

 

 

 

 

 

 

 

Соответствующие

параметры

теоретического

распределе­

ния— это математическое ожидание М(Х) и дисперсия D(X).

 

 

 

В

отличие

от последних,

 

 

 

характеристики

эмпири­

 

 

 

ческого

 

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

 

 

 

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

 

 

 

личинами. Произведя две

 

 

 

серии

наблюдений

над

 

 

 

случайной

величиной X,

 

 

 

мы получим, вообще гово­

 

 

 

ря,

различные

значения

 

 

 

средних

величин

и сред-

p.ис. 11

 

 

них квадратических

от­

 

 

 

клонений.

Чем

больше

число опытов, тем точнее будут оценены параметры

теоретиче­

ского распределения.

Заметим, что формула (1.63) дает для дисперсии смещенную оценку. Это значит, что при конечных п математическое ожида­

ние определяемой по ней величины оказывается

отличным от

дисперсии D(X). Соответствующая несмещенная

оценка будет

£ > (Л )д а -Ц - 2 (*а- Х )2 = *2.

(1-64)

n ~ l

 

Разница между смещенной и несмещенной оценками становит­ ся ощутимой лишь при малом числе измерений.

Отношение

™ = Y

(1-65)

называется коэффициентом изменчивости; его теоретическим аналогом является, очевидно, величина

w =

VD(X)

__

«

( 1.66)

 

М(Х)

М(Х)

 

зя

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

Допустим, что по результатам п опытов найдены эмпириче­

ские несмещенные оценки X H S для математического

ожидания

а и стандарта а. Спрашивается, какова вероятность

того, что

математическое ожидание а и стандарт а лежат в некотором ин­

тервале, включающем X и 5 соответственно? Обратно, какова ширина «доверительного» интервала, внутри которого могут с заданной вероятностью находиться математическое ожидание и стандарт? Если искомое распределение нормальное, то ответ на эти вопросы дается при помощи формул:

Здесь Sk(t) — табулированная функция распределения Стью- дента-для k степеней свободы

а функция ха‘Распределеиия Пирсона P(t2,k) была введена ра­

нее (п. 13).

15. Оценка близости эмпирического и теоретического распределений

Сопоставление эмпирического и теоретического распределе­ ний представляет собой более трудную задачу, чем это может показаться на первый взгляд. Факт близости двух кривых сам

по себе является недостаточным для того, чтобы делать вывод об удачном представлении эмпирического распределения. По скольку параметры теоретического распределения заранее неиз­ вестны, то приходится принимать, что они совпадают с соответ­ ствующими эмпирическими величинами. Между тем последние являются случайными величинами, обладающими разбросом, тем большим, чем меньше число измерений п. Теоретическая кривая распределения, построенная по одному вариационному ряду и хорошо совпадающая с эмпирической кривой для этого ряда, может существенно отличаться от эмпирических кривых, построенных по другим вариационным рядам той же случайной величины.

Имеется несколько способов оценки близости двух распреде­

лений.

Они называются

критериями

согласия.

Критерий

А. Н. Колмогорова позволяет

судить о близости

теоретической

функции распределения F(x)

и эмпирической функции F(x)

по

наибольшей разности между ними

 

 

 

 

 

 

Dn = m ax|F(x) — ?(*)|.

 

 

 

Если функции F{х) и F(x)

непрерывны, а число п достаточно ве­

лико,

то вероятность того,

что

отклонение DnV n

превысит

зна­

чение К определяется по формуле

 

 

 

 

 

 

P(DnV n > \ ) = \ -

2

( -

1)“ e“ 2e2X'

 

 

 

 

 

 

а - -

00

 

 

 

 

 

 

 

 

 

Функция

 

 

 

 

 

 

 

 

к м = 2

 

 

 

 

 

 

 

 

 

а= ——00

 

 

 

 

 

 

 

протабулирована [43].

 

 

 

 

 

 

Практическое

примене­

 

 

 

 

 

ние критерия

А. Н. Кол­

 

 

 

 

 

могорова сводится к сле­

 

 

 

 

 

дующему.

Вычисляется

 

вис. 12

 

 

 

максимальная

разность

 

 

 

 

Dn, после

чего находится

 

 

 

 

 

 

 

 

 

 

вероятность

 

 

 

P ( D > D n) = l - K ( D ltV n ) .

 

(1.68)

Если эта вероятность мала (например, меньше 0,1), то это зна­ чит, что имеет место маловероятное отклонение, уже не объяс­ нимое случайностью измеренных значений. Расхождение меж­

ду ^(х) и F{x) в этом.случае следует считать существенным. Если же вероятность, вычисляемая по формуле (1.681, окажет­ ся достаточно большой, то расхождение можно признать объ­ яснимым случайным характером измеряемой величины. При этом при относительно небольших числах п допустимы большие отклонения. Напротив, при больших числах п уже сравнительно небольшие отклонения могут указывать на отсутствие согласо­ вания между эмпирическим распределением и выбранной ста­ тистической гипотезой.

Другой способ использования критерия А. Н. Колмогорова

«

состоит в том, что по найденной кумулятивной кривой F(x) строятся доверительные границы, в пределах которых с задан­ ной достоверностью лежит искомая теоретическая кривая рас­ пределения. На рис. 12 показана 99%-ная доверительная об­

Соседние файлы в папке книги