Исследование случайных явлений вероятностно-статистическими методами

/>

Н2
– вызваны 1 троечник, 2 хорошиста;

Н3
- вызваны 2 троечника и 1 хорошист;

Н4
– вызваны 2 троечника и 1 отличник.

Задачу можно решить 2 – мя способами.

1) Способ 1 – классический с использованием формулы сочетания, в данном случае неудобен в связи с большими вычислениями, из-за чего эти вычисления я производил в программной среде MS Excel.

Как видно выше, сумма всех гипотез не равна в точности 1. Это связано с ошибками округления чисел вероятностей гипотез.

1/2 - вероятность хорошиста получить или оценку «хорошо», или оценку «отлично»;

1/3 – вероятность троечника получить или оценку «хорошо», или оценку «отлично», или оценку «неудовлетворительно».

1 – вероятность отличника получить оценку «отлично».

2) Способ 2 – вероятность события вычисляется классическим способом.

1/2 - вероятность хорошиста получить или оценку «хорошо», или оценку «отлично»;

1/3 – вероятность троечника получить или оценку «хорошо», или оценку «отлично», или оценку «неудовлетворительно».

1 – вероятность отличника получить оценку «отлично».

Ответ: Р(А)=0,084.

Задание № 5

Случайная величина – годовой доход наугад взятого налогоплательщика. Плотность распределения вероятностей случайной величины задана в виде:

где a – неизвестный параметр распределения, а величины b и c являются константами, значения которых заданы в таблице вариантов задания.

Требуется :

1) Определить значения параметра «а» и построить график функции f(х).

2) Найти функцию распределения F(х) и построить её график.

3) Определить математическое ожидание Mx
, дисперсию Дx
и среднее квадратическое отклонение годового дохода .

4) Вычислить значения третьего М3
и четвертого М4
центральных моментов, и определить коэффициенты ассиметрии и эксцесса.

5) Определить размер годового дохода Х1
в тыс. HE,., не ниже которого с вероятностью Р окажется годовой доход случайного выбранного налогоплательщика.

Решение

1)

2)

3) Математическое ожидание находим по формуле:

Дисперсию проще всего искать по формуле , где - второй начальный момент, а само

С.К.О. есть положительный квадрат от дисперсии, след.

4) Находим третий и четвёртый центральные моменты по формуле через начальные моменты:

Коэффициенты ассиметрии As
, эксцесса Ex
находим по следующим формулам:

5) Для нахождения вероятности воспользуемся формулой для описания вероятност попадания СВ Х на интервал

Подставляем известные значения, получаем:

Ответ: 1) а=0,2391;

2) ;

3)

4) ;

5)

Задание № 6

Производится «n
» независимых испытаний, в каждом из которых события А
может появиться с вероятностью Р
.

Требуется
: 1)Определить вероятность того, что событие А
появится при n
– испытаниях равно k
- раз.

2)Определить вероятность того, что событие А
появится при n
– испытаниях более m
- раз.

3)Определить вероятность того, что событие А
появится при n
– испытаниях не менее k
1

- раз, но не более k
2

- раз.

4)Вычислить среднее число появления события А
при n
– испытаниях и среднее квадратическое отклонение числа появлений события А
.

5)Определить с какой вероятностью должно появляться события А
в каждом из «n
» - опытов при условии, что вероятность не появления события А
ни в одном из «n
» - опытов равна Р0

.

(Указание.
При поиске ответа только на пятый вопрос принять условие, что число появлений события А
при n
– независимых опытах распределено по закону Пуассона).

Решение

1) Применяем формулу Пуассона, т.к. n – велико, р – не мало.

Подставляем, получаем:

2)

3) Использую интегральную формулу Муавра-Лапласа, потому что .

4) Находим МОЖ по следующей формуле . Ввиду больших вычислений расчёт значений произвёл в программном продукте MS Excel.

Среднее квадратическое отклонение находим через дисперсию по формуле:

5) Так как число появлений события А
при n
– независимых опытах распределено по закону Пуассона, то:

Ответ: 1)

2)

3)

4)

5)

Задание № 7

Дискретная двумерная случайная величина (X
,
Y
) описывается законом распределения вероятностей, заданного рядом распределения вероятностей, представленным в таблице:

X Y	X1	X2
Y1	P11	P12
Y2	P21	P22
Y3	P31	P32

Требуется:

1. Определить законы распределения компонент случайного вектора X
и Y
соответственно.

2. Определить условный закон распределения случайной величины X
при условии, что Y
приняла значение yj

.

3. Определить условный закон распределения случайной величины Y
при условии, что X
приняла значение xi

.

4. Определить математические ожидания и дисперсии компонент случайного вектора X
и Y
.

Решение

X Y	X1	X2	Px
Y1	P11=0,11	P12=0,14	0.25
Y2	P21=0,16	P22=0,24	0.4
Y3	P31=0,19	P32=0,16	0.35
Py	0.46	0.54	1

1) Законы распределения компонент случайного вектора X
и Y
вычисляются с помощью складывания вероятностей, находящихся на одной строке – для х, на одном столбце –

для у.

Как видно из таблицы, сумма вероятностей, как и по х, так и по у равна 1.

2) Условный закон распределения случайной величины X
при условии, что Y
приняла значение y
1

определяется как:

3) Условный закон распределения случайной величины Y
при условии, что X
приняла значение x
1

определяется как:

4) Математическое ожидание компонент случайного вектора X
и Y
вычисляем по формуле:

Подставляем значения, получаем:

Дисперсия компонент случайного вектора X
и Y
находиться по формуле, и равна:

Ответ: 1) Законы распределения компонент случайного вектора X
см. в таблице в левом столбце;

законы распределения компонент случайного вектора Y
см. в таблице в нижней строке.

2)

3)

4)

Задание № 8

Случайная величина Y
связана со случайными величинами xi

(i= 1,…,3) функциональной зависимостью вида .

Известны математические ожидания случайных величин и средние квадратические отклонения, Задана также нормированная корреляционная матрица:

Требуется:

1. Вычислить математическое ожидание случайной величины Y
.

2. Вычислить среднее квадратическое отклонение случайной величины Y
.

3. В предположении нормального закона распределения случайных величин xi
записать выражение для плотности распределения случайной величины Y
.

Дано:

mx
1
=2.0, mx2=1.9, mx
3
=1.6;

=0.8, =2.1, =3.2;

Найти:

1)M[Y]-?;

2)-?;

3)f(Y)-?.

Решение

1) В данной задаче применяется 3-й тип постановки задачи. т.к. по известным числовым характеристикам случайных аргументов необходимо найти численные характеристики функций случайных аргументов. Учитывая, что задана линейная функция 3-х случайных аргументов, то для решения задачи применил теорему о математическом ожидании дисперсии линейной функции.

2) Находим корреляционный момент через коэффициент корреляции. Он необходим для расчёта дисперсии, из которой можно будет получить среднее квадратическое отклонение.

3) В предположении нормального закона распределения случайных величин xi
выражение для плотности распределения случайной величины Y
имеет вид:

Ответ: 1) ;

2)

3)

Задание №9

Пусть Х – время задержки момента начала матчей на данном стадионе. Известно среднее значение времени задержки, которое составляет «а» минут.

Требуется:

Оценить вероятность того, что начало матча будет задержано не менее, чем на t минут.Найти минимальное значение времени задержки начала матча t0
, при котором вероятность задержки на время не менее t0
не превышает требуемого значения Ртр
, если дополнительно известно, что среднее квадратичное отклонение времени задержки начала матчей секунд.

Решение

1) Для нахождения вероятности того, что начало матча будет задержано не менее, чем на 2,5 минуты необходимо применить 1-е неравенство П.Л.Чебышева ( А.А. Маркова).

т.е.

2) Минимальное значение времени задержки начала матча t0
, при котором вероятность задержки на время не менее t0
не превышает требуемого значения Ртр
, если дополнительно известно, что среднее квадратичное отклонение времени задержки начала матчей секунд вычисляется по второй теореме П.Л. Чебышева, т.к. она позволяет находить вероятность попадания случайной величины за пределы интервала от .

, т.е.

Ответ: 1) Вероятность того, что мат начало матча будет задержано не менее, чем на t минут составляет не более 44%, т.е. в 1 случае из 100 начало матча будет задержано не более чем на 2,5 минуты.

2) Минимальное значение времени задержки начала матча t0
, при котором вероятность задержки на время не менее t0
не превышает требуемого значения Ртр
=0,4, составляет 2,4703 минуты.

Задание № 10

По заданному в таблице Приложения 4 временному ряду U=¦(T) требуется:

-построить диаграмму рассеивания; найти параметры модели линейного тренда и построить его на диаграмме рассеивания; найти параметры квадратичного тренда и построить его на диаграмме рассеивания;

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

-проверить адекватность построенных моделей по критериям случайности колебаний уровней остаточной последовательности (критерием серий или критерием поворотных точек), критерием соответствия распределения случайной компоненты нормальному закону распределения (критерием ассиметрии и эксцесса или критерием стьюдентизированного размаха), критерием проверки равенства математического ожидания случайной компоненты нулю, критерием проверки независимости значений уровней случайной компоненты (Дарбина-Уотсона);

-выполнить точечный прогноз значений зависимой переменной U по линейному и квадратичному тренду для Τ=8 и T=9;

-выполнить интервальный прогноз с использованием только линейного тренда;

-сформулировать выводы о качестве трендовых моделей.

Дано: Временной ряд:

Время, t	Y
1	6.2
2	4.8
3	7.5
4	9.9
5	14
6	19
7	22

Решение:

Проведём исследование модели временного ряда, которая содержит 3 стадии:

· стадия спецификации модели тренда;

· стадия идентификации модели тренда4

· стадия верификации модели тренда.

1. Проанализируем первую стадию, т.е. спецификацию модели тренда. Суть стадии заключается в том, что подбирается математическая модель, описывающая изменение процессов во времени. Этот процесс начинается с построения диаграммы рассеивания и заканчивается подбором уравнения: линейного, степенного, экспоненциального, полиноминального и прочих. Решается вопрос о линейности (нелинейности) модели тренда и параметров этой модели (коэффициентов).

С помощью пакета программ MS Excel, учитывая исходные данные построил диаграмму рассеивания:

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

2. Стадия идентификации состоит в идентификации коэффициентов модели тренда по экспериментальным данным.

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

Уравнение прямой имеет вид:

При помощи метода наименьших квадратов найдём неизвестные коэффициенты тренда а0
и а1
.

Для уравнения прямой:

1.

2.

3.

4.

5.

6.

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

где

- сумма квадратов отклонений, обусловленная регрессией.

- общая сумма квадратов отклонений зависимой переменной yi
от среднего значения.

- значение зависимой переменной, найденное по уравнению регрессии.

Результаты расчёта занесём в отдельную таблицу

Время, t
1	3,096475
2	6,03575
3	8,975025
4	11,9143
5	14,853575
6	17,79285
7	20,732125

Величина коэффициента детерминации показывает, какая доля зависимой переменной у обусловлена независимой переменной t. Чем ближе R2
к 1, тем лучше регрессия аппроксимирует эмпирические данные. В нашем случае R2
=1,08759≈1, следовательно, полученное уравнение регрессии можно считать вполне достоверным и эффективно оценивающем адекватность регрессионной модели.

Для уравнения параболы:

Уравнение квадратичного тренда имеет вид:

Время, t	y
1	0,3497
2	0,4193
3	0,4890
4	0,5586
5	0,6283
6	0,6980
7	0,7676

Ввиду того, что значения квадратичного тренда крайне малы, я построил его на отдельном графике.

3. Этап верификации модели тренда состоит в проверке адекватности модели реальным экспериментальным данным. На этой стадии производится:

- критерий Фишера для проверки гипотезы (Н0
) о значимости коэффициента детерминации R2
,

- критерий поворотных точек,

- критерий ассиметрии и эксцесса по Стьюденту.

1. Значимость коэффициента детерминации определяется критерием Фишера.

,

Для линейного уравнения:

.

Расчётное значение критерия Фишера сравнивается с табличным. При заданном уровне значимости гипотезы Н=0,05 и степенях свободы f1
=n-1 (f1
=6) и f2
=n-m-1 (f2
=5), где n – объём выборки, m – число включённых в модель факторов табличное значение равно 4б95. Сравнивая расчётное и полученное значение, получаем. что 57,0848>4,95, значит модель адекватна исходным данным.

2. Критерий поворотных точек. С помощью этого критерия можно проверить случайность остатков модели. В соответствии с этим критерием каждый уровень ряда сравнивается с двумя соединёнными с ними. Если он больше или меньше их. то эта точка считается поворотной. Затем подсчитывается сумма поворотных точек m. В случайном ряду чисел должно выполняться строгое неравенство:

,

где N – объём выборочной совокупности.

Данное условие выполняется, следовательно, можно утверждать, что ряд остатков моделей не является случайным.

3. Эксцесс характеризует относительную остроконечность или сглаженность распределения по сравнению с нормальным распределением. Положительный эксцесс означает относительно остроконечное распределение. Отрицательный эксцесс обозначает относительно сглаженное распределение.

Эксцесс я высчитал с помощью пакета программ MS Excel, он равен -1,295349862, это говорит об относительно сглаженном распределении.

Асимметрия характеризует степень нессиметричности распределения относительно его среднего. положительная асимметрия указывает на отклонение распределения в сторону положительных значений. Отрицательная асимметрия указывает на отклонение распределения в сторону отрицательных значений.

Асимметрию я рассчитал помощью пакета программ MS Excel, она равна 0,602642408, что говорит об отклонении в сторону отрицательных значений.

Учитывая вышенаписанное, можно утверждать о том, что данное распределение случайной компоненты не соответствует нормальному распределению. т.к. эти величины не малы.

Математическое ожидание равно 12.

Выполним точечный прогноз значений зависимой переменной Y по линейному тренду для t=8 и t=9.

Подставив эти значения в уравнения прямой, мы получаем, что при t=8 y= 23,6714 и для t=9 y= 26,6107.

Вывод о качестве трендовых моделей:

Все модели качественны и адекватны (пригодны для практического применения). С течением времени происходит тенденция возрастания линий трендов.

Ответ
: R2
=, критерий Фишера равен 57,0848, критерий поворотных точек равно 0,09, он равен -1,295349862, это говорит об относительно сглаженном распределении, асимметрия равна 0,602642408, что говорит об отклонении в сторону отрицательных значений, МОЖ равно 12, точечный прогноз значений зависимой переменной Y по линейному и квадратичному тренду при t=8 y= 23,6714 и для t=9 y= 26,6107.

Вывод

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

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

Эта курсовая работа имела большое практическое значение, т.к. решались типичные повседневные задачи, с которыми может столкнуться менеджер на работе.

Так же эта курсовая работа позволила ещё раз отработать навыки и закрепить знания, полученные при изучении дисциплины математический анализ.

Список использованной литературы

1. Вентцель Е.С. Теория вероятностей. Учебник для вузов – М.: «Высшая школа», 2000г.

2. Гмурман В.Е. Руководство к решению задач по теории вероятностей и математической статистике. – М.: «Высшая школа». – 2003 г.

3. Вентцель Е.С., Овчаров Л.А. Прикладные задачи теории вероятностей. – М.: «Радио и связь». – 1983г.

4. Виноградов С.А. и др. Теория вероятностей и математическая статистика. 5. Учебно-методическое пособие для практических занятий. МО, 1998 г.