Скачиваний:
12
Добавлен:
30.09.2023
Размер:
602.98 Кб
Скачать

Министерство образования и науки Российской Федерации Федеральное агентство по образованию

Московский государственный институт электронной техники (технический университет)

______________________________________________________________

В.В.Бардушкин, В.В.Лесин, В.Н.Земсков, Н.Н.Мустафин

Лабораторный практикум по курсу «Теория вероятностей и математическая статистика»

Утверждено редакционно-издательским советом института

Москва 2009

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

УДК 519.2

Рецензент докт. физ.-мат. наук, проф. В.Б. Яковлев

Бардушкин В.В., Лесин В.В., Земсков В.Н., Мустафин Н.Н.

Лабораторный практикум по курсу «Теория вероятностей и математическая статистика». – М.: МИЭТ, 2009. – 116 с.: ил.

Лабораторный практикум состоит из трёх разделов и двух приложений. Первый раздел содержит необходимые теоретические сведения по основным методам обработки статистических данных. Во втором разделе приведены общие сведения о системе программирования Matlab с пакетом прикладных программ по статистике «Statistics Toolbox», элементы программирования на языке Matlab, а также описаны средства пакета для изучения одномерных распределений случайных величин. Третий раздел содержит описание четырёх лабораторных работ. В приложениях приведены варианты индивидуальных заданий.

Предназначен для студентов 2-го курса технических факультетов МИЭТ.

© МИЭТ, 2009

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

Бардушкин Владимир Валентинович Лесин Виктор Васильевич Земсков Владимир Николаевич Мустафин Наиль Нухович

Лабораторный практикум по курсу «Теория вероятностей и математическая статистика»

Редактор Е.Г. Кузнецова. Технический редактор Л.Г. Лосякова. Вёрстка авторов.

Подписано в печать с оригинал-макета 03.03.09. Формат 60х84 1/16. Печать офсетная. Бумага офсетная. Гарнитура

Times New Roman. Усл. печ. л. 6,73.

Уч.-изд. л. 5,8. Тираж 1000 экз. Заказ 25.

Отпечатано в типографии ИПК МИЭТ.

124498, Москва, Зеленоград, проезд 4806, д. 5, МИЭТ.

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

1.Краткие теоретические сведения

1.1.Основные распределения, используемые

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

Равномерное распределение. Случайная величина непрерывного типа Х распределена равномерно на отрезке [a; b] (обозначение X ~ R(a; b) ), если плотность вероятности имеет вид:

ì

0,

 

x Ï[a; b],

ï

1

 

 

f X (x) = í

,

x Î[a; b].

ï

 

 

îb - a

 

 

Графики плотности распределения, построенные в системе Matlab для нескольких случайных величин R(a; b) , приведены на рис.1.1.

Рис.1.1. Графики плотности распределения R(a ; b)

Характеристическая

функция распределения

R(a; b) имеет вид j(t) =

eibt - eiat

, t R ,

it(b - a)

 

 

 

 

математическое ожидание

mX = (a + b) 2 , дисперсия

DX = (b - a)2 12 . Равномерное распределение не

имеет моды, а медиана совпадает с математическим ожиданием.

 

 

С помощью линейного преобразования Y = (b a)X + a от стандартной равномерно распределённой величины X ~ R(0;1) можно перейти к случайной величине Y ~ R(a; b) с произвольными параметрами a и b.

Ошибка, возникающая в результате округления числа с точностью до 10n , удовлетворительно описывается равномерным на промежутке [-0,5×10n ; 0,5×10n ] распределением.

Нормальное распределение. Случайная величина непрерывного типа Х имеет нормальное распределение, или распределение Гаусса, с параметрами m и σ (обозначение X ~ N(m; σ) ), если плотность распределения имеет вид:

1

 

 

(xm)2

 

 

 

 

2σ2

, m R , σ > 0 .

f X (x) =

 

 

 

e

 

 

s

 

 

 

 

2p

 

 

Графики плотности распределения, построенные в системе Matlab для нескольких случайных величин N(m; σ) , приведены на рис.1.2.

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

Рис.1.2. Графики плотности распределения N(m; σ)

 

æ

 

s

2

t

2 ö

 

 

Характеристическая функция распределения N(m; σ)

ç

imt -

 

÷

,

t R ,

имеет вид j(t) = exp

 

 

 

÷

 

ç

 

 

2

 

 

 

è

 

 

ø

 

 

математическое ожидание mX = m , дисперсия DX = s2 .

Мода и медиана X ~ N(m; σ)

 

совпадают с

математическим ожиданием.

Фундаментальная роль нормального распределения объясняется тем, что при некоторых (обычно выполняющихся) условиях закон распределения суммы случайных величин, имеющих произвольные законы распределения, приближается к нормальному закону с ростом числа суммируемых величин. Точная формулировка этого утверждения является содержанием центральной предельной теоремы [1].

С помощью линейного преобразования Y = σX + m от случайной величины со стандартным

нормальным распределением X ~ N(0;1) можно

перейти к случайной величине Y ~ N(m; σ) с

произвольными параметрами m и σ . Функция распределения стандартной нормальной величины

 

 

1

 

 

x

t2

 

F(x) =

 

 

 

òe

2 dt

 

 

 

 

 

2p

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

называется функцией Лапласа, значения которой табулируются (см., например, табл.П1 в [1]).

Нормально распределённая случайная величина с большой вероятностью принимает значения, близкие к своему математическому ожиданию, что выражается правилом k сигм:

ì0,6827 ,

k =1,

ï

,

k = 2 ,

P{| X - m | < ks} = í0,9545

ï

,

k = 3.

î0,9973

Обычно используется правило трёх сигм: значения нормальной случайной величины X ~ N(m; σ) практически достоверно отклоняются от её математического ожидания m менее чем на .

Квантиль u p нормального распределения является корнем уравнения F(u p ) = p , т.е. u p = F−1( p) .

Для обратной функции Лапласа F−1( p) существуют специальные таблицы (см., например, табл.П1 в [1]). В

силу симметричности нормального распределения

u p = -u1− p .

Вероятность попадания случайной величины X ~ N(m; σ)

в интервал (a; b) вычисляется по формуле

æ b - m ö

 

æ a - m ö

P{a < X < b} = Fç

 

÷

- Fç

 

÷ .

s

s

è

ø

 

è

ø

Нормальное распределение композиционно устойчиво, т.е. сумма независимых нормально распределённых случайных величин также имеет нормальное распределение.

Гамма-распределение. Случайная величина непрерывного типа Х имеет гамма-распределение с параметрами a и b (обозначение X ~ Γ(a; b) ), если её плотность распределения имеет вид:

ì

1

 

 

a−1

x

 

 

ï

 

 

b

 

 

f X (x) = í

G(a)ba

 

x

e

 

 

,

x > 0, a > 0, b > 0,

ï

 

0,

 

 

 

x £ 0,

î

 

 

 

 

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

+∞

где G(a) = òtα−1et dt гамма-функция, обладающая свойствами

0

Γ(x +1) = xΓ(x) , Γ(k +1) = k!, где k N ,

Γ(1) = Γ(2) = 1 , G(0,5) = p .

Графики плотности распределения, построенные в системе Matlab для нескольких случайных величин Γ(a; b) , приведены на рис.1.3.

Характеристическая функция гамма-распределения имеет вид j(t) = (1- ibt)a , t R ,

математическое ожидание mX = ab , дисперсия DX = ab2 .

Гамма-распределение композиционно устойчиво при фиксированном параметре b, т.е. сумма независимых случайных величин, имеющих гамма-распределение с параметром b, также имеет гамма- распределение с параметром b.

Рис.1.3. Графики плотности распределения Γ(a; b)

Экспоненциальное (показательное) распределение. Это распределение частный случай гамма-

распределения при a = 1 (обозначение X ~ E(λ) , где l = 1b ). Таким образом, функция плотности гамма- распределения имеет вид:

ì

−λx

, x > 0, l > 0,

ïle

 

f X (x) = í

 

 

ï

0,

x £ 0.

î

Графики плотности распределения, построенные в системе Matlab для нескольких случайных величин E(λ) , приведены на рис.1.4.

 

 

Рис.1.4. Графики плотности распределения E(λ)

 

 

Характеристическая

функция

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

X ~ E(λ)

имеет вид

j(t) = l(l - it)1 ,

t R ,

математическое ожидание mX = 1 l , дисперсия DX

= 1 l2 .

 

 

 

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

χ2 хи-квадрат») с k степенями свободы. Это распределение частный случай

гамма-распределения

при

a = k 2 ,

k N , b = 2

(обозначение

X ~ χ2 (k) ).

Таким образом,

функция

плотности распределения χ2 (k) имеет вид:

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

ì

 

1

 

x k

−1

ï

 

 

 

 

 

 

 

 

e

2 x 2

, x > 0,

2k 2

Г(k 2)

f X (x) = í

 

 

 

 

 

ï

 

0,

 

 

 

 

x £ 0.

î

 

 

 

 

 

Графики плотности распределения, построенные в системе Matlab для нескольких случайных величин χ2 (k) , приведены на рис.1.5.

Рис.1.5. Графики плотности распределения χ2 (k)

Характеристическая функция распределения χ2 (k) имеет вид j(t) = (1- 2it)k2 , t Î R , математическое ожидание mX = k , дисперсия DX = 2k .

По закону χ2 (k) распределена сумма квадратов k независимых стандартных нормальных случайных величин, т.е. если случайные величины X i ~ N(0;1) , i = 1, 2, ..., k , то

k

Y = åX i ~ χ2 (k) .

i=1

Существуют таблицы квантилей распределения χ2 (k) при различных k (см., например, табл.П5 в [1]).

С ростом числа степеней свободы k распределение χ2 (k) , согласно центральной предельной теореме,

приближается к нормальному распределению, т.е. если Y ~ χ2 (k) , то при больших k (на практике при

k ³ 30) можно приближённо считать, что Y ~ N(k ;

 

 

 

 

) .

 

При этом квантили c2p (k) распределения χ2

 

 

2k

могут быть приближённо найдены с помощью квантилей u p

 

 

стандартного нормального закона N(0;1) :

 

c2p (k) » u p

 

+ k .

 

 

 

2k

 

 

Более точные значения квантилей c2p (k)

распределения χ2 можно получить по формулам

c2p (k) » 0,5(

 

 

+ u p )2

 

2k -1

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ö3

 

æ

 

2

 

 

 

 

 

 

 

 

 

 

c2

 

 

 

 

 

 

 

 

2

 

(k) » kç1

-

 

 

 

+ u

p

 

 

÷ .

 

 

 

 

 

p

ç

 

9k

 

 

 

 

 

 

÷

 

è

 

 

 

 

 

 

 

 

9k ø

Один из критериев проверки согласия

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

 

данных с гипотетической функцией

распределения FX (x) основан на изучении c2 -статистики Пирсона

m

(n

j

- np

j

)2

 

c2 = å

 

 

 

 

,

 

 

np

j

 

 

j=1

 

 

 

 

 

 

 

 

 

 

 

где p j = FX (x j ) - FX (x j−1) , x0 = -¥ < x1 < ... < xm = +¥

произвольное разбиение интервала

(−∞; + ∞) , n j число результатов измерений,

попавших

в промежуток [x j−1 ; x j ) . c2 -статистики

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

Пирсона

(в предположении истинности

гипотезы о законе распределения) имеют в пределе при

m

 

 

n = ån j

® ¥ распределение χ2 (m -1)

независимо от вида FX (x) .

j=1

 

 

Распределение Стьюдента. Случайная величина непрерывного типа Х имеет распределение Стьюдента с k степенями свободы (обозначение X ~ St(k) ), если её плотность распределения имеет вид:

 

 

 

æ k +1ö

 

 

 

 

 

 

k+1

 

 

 

 

Гç

 

 

÷

æ

 

x2 ö

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

è

ø

 

 

 

 

f

X

(x) =

 

 

 

 

 

ç1

+

 

÷

 

,

x R , k N .

 

 

 

 

 

 

 

 

 

G(k 2)

pk

ç

 

k

÷

 

 

 

 

 

 

 

è

 

ø

 

 

 

 

Графики плотности распределения, построенные в системе Matlab для нескольких случайных величин St(k) , приведены на рис.1.6.

Рис.1.6. Графики плотности распределения St(k)

Основные числовые характеристики распределения St(k) :

математическое ожидание mX

ì

0,

k >1,

= í

 

k =1;

 

 

 

 

 

îне существует,

 

дисперсия DX

ìk (k - 2),

 

k > 2,

 

 

= í

¥,

k £ 2.

 

 

 

 

î

 

 

Квантили t p (k)

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

 

St(k)

находятся из таблиц (см., например, табл.П6 в [1]). В силу

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

t p (k) = -t1 p (k) .

При больших k распределение Стьюдента приближается к стандартному нормальному распределению.

Поэтому при k ³ 30 можно использовать

приближённое

равенство t p (k) » u p , где u p квантиль

распределения N(0;1) . Более точно квантили t p (k) можно находить по формуле

 

 

 

 

æ

 

 

 

2

 

2 ö

1

 

 

 

 

 

æ

1 ö

 

2

 

t

 

(k) = u

ç

 

 

u p ÷

 

 

p

p ç

ç1-

 

÷

-

 

 

 

.

 

 

2k ÷

 

 

 

è

4k ø

 

 

 

 

 

 

 

 

è

 

 

 

 

ø

 

 

Распределение St(k) широко применяется в математической статистике. Это объясняется тем, что

если случайные величины U и V независимы, причём U ~ N(0;1)

и V ~ χ2

(k) , то отношение

 

U

 

 

Vk

имеет распределение Стьюдента с k степенями свободы:

 

U

 

~ St(k) .

 

 

 

V k

 

 

Распределение Фишера. Случайная величина непрерывного типа Х имеет распределение Фишера с k1 и k2 степенями свободы (обозначение X ~ F(k1 ; k2 ) ), если её плотность распределения имеет вид:

PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

Соседние файлы в папке Ресурсы