Matlab  |  Mathcad  |  Maple  |  Mathematica  |  Statistica  |  Другие пакеты Поиск по сайту
Internet-класс  |  Примеры  |  Методики  |  Банк задач  |  Консультации & Форум  |  Download  |  Ссылки  |  Конкурсы
Научно-практический журнал "Exponenta Pro. Математика в приложениях". Вышел 1/2004 номер журнала
 
Лабораторные работы по курсу "Теория вероятностей"

Вернуться на страницу <Методические разработки>
В начало | Л.р.1 | Л.р.2 | Л.р.3Л.р.5 | Л.р.6 | Л.р.7 | Л.р.8

| Л.р.9 | Л.р.10 | Л.р.11 | Л.р.12 | Л.р.13 | Л.р.14 | Л.р.15

Лабораторная работа 4
Предельные теоремы Муавра-Лапласа

Локальная формула ~ Интегральная формула

Теоремы Муавра-Лапласа применяются для приближенного вычисления вероятностей P[n](m) и P[n](m[1],m[2]) = sum(P[n](m),m = m[1] .. m[2]) в независимых испытаниях Бернулли при больших n, m, m[1], m[2] , когда вероятность p наступления успеха в одном испытании постоянна, величины (m-np)/sqrt(npq) , (m[1]-np)/sqrt(npq) , (m[2]-np)/sqrt(npq) равномерно ограничены по n, m, m[1], m[2] .

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

P[n](m) = phi((m-np)/sqrt(npq))/sqrt(npq) ,

где p - вероятность появления успеха в каждом испытании, q=1-p, phi(x) = exp(-x^2/2)/sqrt(2*Pi) .

Интегральная формула . Вероятность того, что в n независимых испытаниях число успехов находится между m[1] и m[2] , равна

P[n](m[1],m[2]) = Phi((m[2]-np)/sqrt(npq))-Phi((m[1... ,

где p- вероятность появления успеха в каждом испытании, q=1-p, Phi(x) = int(exp(-t^2/2),t = -infinity .. x)/sqrt(2... - функция Лапласа.

В частности,

P(abs(mu/n-p) <= epsilon) = 1-Phi(-epsilon*sqrt(n/p... ,

где mu/n - частота появления успеха, p - его вероятность, epsilon - некоторое положительное число.

Приближенные формулы Муавра-Лапласа применяют в случаях, когда p и q не малы, а 9 <= npq .

Задача 1 . Вычислить вероятности того, что при 100-кратном бросании монеты герб выпадет а) ровно 50 раз; б) ровно 60 раз.

Решение. Условия применения локальной теоремы Муавра-Лапласа выполняются, поэтому воспользуемся ею:

> n:=100;p:=0.5;q:=0.5;

n := 100

p := .5

q := .5

> x:=(k-n*p)/sqrt(n*p*q);

x := .2000000000*k-10.00000000

> f(x):=1/sqrt(2*Pi)*exp(-x^2/2);

f(.2000000000*k-10.00000000) := 1/2*sqrt(2)*exp(-1/...

> pn(k):=1/sqrt(n*p*q)*f(x);

pn(k) := .1000000000*sqrt(2)*exp(-1/2*(.2000000000*...

Построим график p[n](k)

> plot(pn(k),k=40..60);

[Maple Plot]

Очевидно, что наивероятнейшее число наступления успеха равняется 50.

Задача 2 . Вероятность успеха в каждом испытании равна 0.25. Какова вероятность, что при 300 испытаниях успех наступит а) ровно 75 раз; б) ровно 85 раз?

Задача 3 . В первые классы должно быть принято 200 детей. Определить вероятность того, что среди них окажется 100 девочек, если вероятность рождения мальчика равна 0.515.

Задача 4 . Построить график функции Лапласа.

Решение. Очистим x:

> x:='x':

Зададим функцию Лапласа

> Phi(x):=1/sqrt(2*Pi)*int(exp(-t^2/2),t=-infinity..x);

Phi(x) := 1/2*sqrt(2)*(1/2*erf(1/2*x*sqrt(2))*sqrt(...

Построим график функции Лапласа

> plot(Phi(x),x);

[Maple Plot]

Нетрудно проверить, что Phi(-x)+Phi(x) = 1 .

Задача 5 . Вычислить вероятность наступления случайного события от 790 до 830 раз в 900 независимых испытаниях, если p=0,9.

Решение. Выразим x1 и x2:

> x2:=(830-900*0.9)/sqrt(900*0.9*0.1);

x2 := 2.222222222

> x1:=(790-900*0.9)/sqrt(900*0.9*0.1);

x1 := -2.222222222

Вычислим искомую вероятность:

> P:=1/sqrt(2*Pi)*int(exp(-t^2/2),t=-infinity..x2)-1/sqrt(2*Pi)*int(exp(-t^2/2),t=-infinity..x1);

P := 1.220391717*sqrt(2)/(sqrt(Pi))

Оценим ее:

> evalf(P);

.9737317088

Задача 6 . Какова вероятность того, что в столбике из 100 наугад отобранных монет число монет, расположенных "гербом" вверх, будет от 45 до 55?

Задача 7 . Производство дает 1% брака. Какова вероятность того, что из взятых на исследование 1100 изделий выбраковано будет не больше 17?

Задача 8 . Сколько нужно произвести опытов с бросанием монеты, чтобы с вероятностью 0.92 можно было ожидать отклонение частоты выпадения "герба" от теоретической вероятности 0.5 на абсолютную величину, меньшую чем 0.01?

Задача 9 . Игральную кость бросают 80 раз. Найти приближенно границы, в которых число выпаданий шестерки будет заключено с вероятностью 0,9973.

Задача 10 . (Чистяков В. П. Курс теории вероятностей: Учеб.-3-е изд., испр.-М.: Наука. Гл. ред. физ.-мат. лит.- 1987.) В поселке А 2500 жителей. Каждый из них примерно 6 раз в месяц ездит на поезде в город В, выбирая дни поездок независимо от остальных жителей. Какой наименьшей вместительностью должен обладать поезд, чтобы он переполнялся в среднем не чаще одного раза в 100 дней? (Поезд идет раз в сутки.)

Решение. Пусть наименьшая вместимость поезда составляет m мест. Вероятность выезда в город В жителя поселка А равна 6/30, вероятность переполнения равна 1/100. Определим функции

> f(t):=1/sqrt(2*Pi)*exp(-t^2/2):F:=x->int(f(t),t=-infinity..x);

F := proc (x) options operator, arrow; int(f(t),t =...

Вероятность переполнения поезда (вероятность попадания количества пассажиров в интервал (m, 2500)) выразится так (интегральная формула Муавра-Лапласа):

> p:=F((2500-2500/5)/sqrt(2500/5*4/5))-F((m-2500/5)/sqrt(2500/5*4/5));

p := 1/2*sqrt(2)*(1/2*sqrt(2)*sqrt(Pi)+1/2*sqrt(2)*...

Найдем искомое значение m

> ceil(solve(p-0.01,m));

547

Команда ceil возвращает наименьшее целое, большее или равное аргументу.

В начало | Л.р.1 | Л.р.2 | Л.р.3 |Л.р.5 | Л.р.6 | Л.р.7 | Л.р.8

| Л.р.9 | Л.р.10 | Л.р.11 | Л.р.12 | Л.р.13 | Л.р.14 | Л.р.15
Вернуться на страницу <Методические разработки>

Карта сайта | На первую страницу | Поиск |О проекте |Сотрудничество |
Exponenta Pro | Matlab.ru

Наши баннеры


Copyright © 2000-2003. Компания SoftLine. Все права защищены.

Дата последнего обновления информации на сайте: 11.05.04
Сайт начал работу 1.09.00

Программное обеспечение Microsoft, Macromedia, VERITAS, Novell, Borland, Symantec, Oracle и др.