Предположим, что имеется n независимых испытаний с двумя исходами в каждом испытании. Один из исходов будем называть успехом и кодировать цифрой 1, другой исход будем называть неудачей и кодировать цифрой 0. Предполагаем, что вероятность успеха в каждом испытании одна и та же и равна числу p , следовательно, вероятность неудачи равна . Эта схема, очевидно, является обобщением схемы независимого бросания монеты.

Пусть – вероятность того, что общее число успехов равно m. Тогда основная формула схемы Бернулли имеет вид .

Когда числа n и m становятся большими, вычисления по этой формуле становятся затруднительны. Поэтому используются три предельные теоремы: теорема Пуассона, локальная теорема Муавра–Лапласа и интегральная теорема Муавра–Лапласа. Приведем их формулировки.

Теорема Пуассона . (Формулировка приводится в упрощенном виде.) Пусть имеется n независимых испытаний. – вероятность успеха в одном испытании, – вероятность неудачи. Пусть . Тогда для любого фиксированного m справедливо соотношение

при

Комментарий. На практике эта теорема применяется при Это означает, что p должно быть очень малым числом, а n большим.

Локальная теорема Муавра-Лапласа . Пусть имеется n независимых испытаний Бернулли с вероятностью успеха p () в одном испытании и – вероятностью неудачи. Величина не зависит от n . Предположим, что для некоторой постоянной выполнено условие , где Тогда при

.

Комментарий. Эта теорема применяется, когда p отделено от нуля и единицы.

Интегральная теорема Муавра-Лапласа . Пусть имеется n независимых испытаний с вероятностью успеха p () в одном испытании и вероятностью неудачи. Величина не зависит от n . Тогда для любых вещественных чисел при

.

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

Приведем задачи на применение схемы Бернулли и соответствующих предельных теорем.

Задача 30. Случайное блуждание по прямой.Частица движется по целым точкам вещественной прямой, перемещаясь каждую секунду либо на единицу вправо, либо на единицу влево с равными вероятностями. Найти вероятность того, что через n секунд частица вернется в точку 0.

Решение. Очевидно, вернуться в 0 частица может только за четное число секунд. Поэтому считаем, что . Считая успехом движение частицы вправо, заметим, что для возвращения за n секунд должно быть ровно k успехов. Поэтому из формулы Бернулли следует, что вероятность возвращения равна .

Задача 31. Имеется 5 студенческих групп по 25 человек, в каждой из которых по 5 отличников. Из каждой группы выбирается случайным образом по одному студенту. Найти вероятность того, что среди выбранных студентов будет 3 отличника.

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

Задача 32. (Задача Банаха) У рассеянного курильщика в правом и левом карманах пиджака находится по коробку спичек. В каждом коробке по n спичек. Каждый раз, когда ему требуется закурить, курильщик вынимает новую спичку либо из левого, либо из правого кармана с вероятностью 1/2. Найти вероятность того, что в тот момент, когда окажется пустым один из коробков, во втором коробке останется k спичек.

Решение. Пусть A – это событие, сформулированное в вопросе задачи. Будем считать испытанием Бернулли вытаскивание спичек, причем вытаскивание спички из правого кармана будем считать успехом, а из левого – неудачей. Очевидно, вероятность успеха равна 1/2. Поскольку к моменту окончания «эксперимента» из одного коробка вытащили n спичек, а из другого – спичек, то общее число испытаний Бернулли можно считать равным , причем событие A реализуется, если число успехов равно n или k . Поэтому . Здесь использовано свойство биномиальных коэффициентов, согласно которому слагаемые в скобках равны между собой.

Задача 33. Монета бросается 100 раз. Найти приближенно вероятность того, что герб выпадет 40 раз. (Воспользоваться таблицей.)

,

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

Задача 34. Город ежедневно посещают 1000 туристов, которые днем идут обедать. Каждый из них выбирает для обеда один из двух городских ресторанов с равными вероятностями и независимо друг от друга. Владелец одного из ресторанов желает, чтобы с вероятностью близкой к 0,99, все пришедшие в ресторан туристы смогли бы там одновременно пообедать. Сколько мест должно быть для этого в ресторане?

Решение. Обозначим через событие, состоящее в том, что i-й турист пообедал у заинтересованного владельца ресторана i= 1, 2,…, 1000. Наступление события будем называть успехом в i- м испытании. Вероятность успеха . Пусть m – общее число успехов, событие A состоит в переполнении ресторана, k – общее число мест в ресторане. Тогда нам надо подобрать k таким образом, чтобы выполнялось приближенное равенство

Производится n опытов по схеме Бернулли с вероятностью успеха p . Пусть X - число успехов. Случайная величина X имеет область значений {0,1,2,...,n}. Вероятности этих значений можно найти по формуле: , где C m n - число сочетаний из n по m .
Ряд распределения имеет вид:

x	0	1	...	m	n
p	(1-p) n	np(1-p) n-1	...	C m n p m (1-p) n-m	p n

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

Назначение сервиса . Онлайн-калькулятор используется для построения биноминальным ряда распределения и вычисления всех характеристик ряда: математического ожидания, дисперсии и среднеквадратического отклонения. Отчет с решением оформляется в формате Word (пример).

Число испытаний: n = , Вероятность p =
При малой вероятности p и большом количестве n (np формула Пуассона.

Видеоинструкция

Схема испытаний Бернулли

Числовые характеристики случайной величины, распределенной по биноминальному закону

Математическое ожидание случайной величины Х, распределенной по биноминальному закону.
M[X]=np

Дисперсия случайной величины Х, распределенной по биноминальному закону.
D[X]=npq

Пример №1 . Изделие может оказаться дефектным с вероятностью р = 0.3 каждое. Из партии выбирают три изделия. Х – число дефектных деталей среди отобранных. Найти (все ответы вводить в виде десятичных дробей): а) ряд распределения Х; б) функцию распределения F(x) .
Решение . Случайная величина X имеет область значений {0,1,2,3}.
Найдем ряд распределения X.
P 3 (0) = (1-p) n = (1-0.3) 3 = 0.34
P 3 (1) = np(1-p) n-1 = 3(1-0.3) 3-1 = 0.44

P 3 (3) = p n = 0.3 3 = 0.027

x i	0	1	2	3
p i	0.34	0.44	0.19	0.027

Математическое ожидание находим по формуле M[X]= np = 3*0.3 = 0.9
Проверка: m = ∑x i p i .
Математическое ожидание M[X] .
M[x] = 0*0.34 + 1*0.44 + 2*0.19 + 3*0.027 = 0.9
Дисперсию находим по формуле D[X]=npq = 3*0.3*(1-0.3) = 0.63
Проверка: d = ∑x 2 i p i - M[x] 2 .
Дисперсия D[X] .
D[X] = 0 2 *0.34 + 1 2 *0.44 + 2 2 *0.19 + 3 2 *0.027 - 0.9 2 = 0.63
Среднее квадратическое отклонение σ(x) .

Функция распределения F(X) .
F(xF(0F(1F(2F(x>3) = 1

1. Вероятность появления события в одном испытании равна 0.6 . Производится 5 испытаний. Составить закон распределения случайной величины Х – числа появлений события.
2. Составить закон распределения случайной величины Х числа попаданий при четырех выстрелах, если вероятность попадания в цель при одном выстреле равна 0.8 .
3. Монету подбрасывают 7 раз. Найти математическое ожидание и дисперсию числа появлений герба. Примечание: здесь вероятность появление герба равна p = 1/2 (т.к. у монеты две стороны).

Пример №2 . Вероятность появления события в отдельном испытании равна 0.6 . Применяя теорему Бернулли, определите число независимых испытаний, начиная с которого вероятность отклонения частоты события от его вероятности по абсолютной величине меньше 0.1 , больше 0.97 . (Ответ: 801)

Пример №3 . Студенты выполняют контрольную работу в классе информатики. Работа состоит из трех задач. Для получения хорошей оценки нужно найти правильные ответы не меньше чем на две задачи. К каждой задаче дается 5 ответов из которых только одна правильная. Студент выбирает ответ наугад. Какая вероятность того, что он получит хорошую оценку?
Решение . Вероятность правильно ответить на вопрос: p=1/5=0.2; n=3.
Эти данные необходимо ввести в калькулятор. В ответ см. для P(2)+P(3).

Пример №4 . Вероятность попадания стрелка в мишень при одном выстреле равна (m+n)/(m+n+2) . Производится n+4 выстрела. Найти вероятность того, что он промахнется не более двух раз.

Примечание . Вероятность того, что он промахнется не более двух раз включает в себя следующие события: ни разу не промахнется P(4), промахнется один раз P(3), промахнется два раза P(2).

Пример №5 . Определите распределение вероятностей числа отказавших самолётов, если влетает 4 машины. Вероятность безотказной работы самолета Р=0.99 . Число отказавших в каждом вылете самолётов распределено по биноминальному закону.

Ранее в п. 1.4 введены понятия зависимых и независимых событий. С понятием независимых событий связано и имеет широкое применение понятие независимых опытов или испытаний.

Опыты α 1 , α 2 , … , α n называются независимыми, если любая комбинация их исходов является совокупностью независимых событий. Иначе, если в задаче проводится ряд многократно повторяющихся испытаний α 1 , α 2 , …, α n при неизменном комплексе условий и в каждом испытании некоторые событие А может наступить с некоторой вероятностью p = p (А ) не зависящей от других испытаний, и не наступить с вероятностью p (Ā ), то указанные испытания называются независимыми. Данная схема независимых испытаний носит название схемы Бернулли.

Схема названа в честь Якоба Бернулли – родоначальника семьи выдающихся швейцарских учёных. (Якоб Б., Иоганн Б., Николай Б., Даниил Б. и др.). Якоб Бернулли доказал так называемую теорему Бернулли – важный частный случай закона больших чисел (см. п. 3.11). Указанная теорема относится к рассматриваемой здесь последовательности независимых испытаний.

Примерами независимых испытаний являются: а) многократное (n раз) подбрасывание монеты; б) извлечение (n раз) одинаковых на ощупь шаров из урны с их последующим возвращением; в) любая совокупность независимых испытаний (опытов), в каждом из которых вероятность успешных исходов одинакова, например, серия выстрелов по мишени, выбор n деталей из их совокупности, изучение n анализов горной породы определённого свойства и т.д.

В схеме Бернулли наступление события А с вероятностью p = p (А ) условно называется успехом, а его ненаступление (противоположное событие Ā ) –неудачей. Вероятность неудачи в каждом опыте такого типа равна q = 1 – p .

На практике обычно возникают задачи со сложными событиями, в которых из n опытов, составляющих схему Бернулли, в m опытах (m < n ) событие А наступает (т.е завершается успехом), а в (n – m ) опытах это событие не наступает (завершается неудачей). Пусть P n (k ) – обозначает вероятность того, что при производстве n опытов успех наступает в k опытах (успех реализуется k раз). Ставится следующая задача: пусть в n испыта-ниях, соответствующих схеме Бернулли, k испытаний завершились успехом. Требуется найти вероятность P n (k ) (читается: « P из n испытаний k успешных» ). Данная вероятность рассчитывается по формуле Бернулли, которой соответствует одноименная теорема.

Теорема Бернулли. Если вероятность p наступления события А в каждом из последовательности n испытаний α 1 , α 2 , … , α n постоянна, то вероятность того, что событие А наступит k раз и не наступит n – k раз, вычисляется по формуле Бернулли:

P n (k ) = С n k p k q n-k , (2.1)

где q = 1- p .

Доказательство. Действительно, пусть события A į и Ā į – появление и непоявление соответственно события А в į -ом испытании α i (i = 1, 2, … , n ). Пусть также В k обозначает событие, состоящее в том, что в n независимых испытаниях событие А появилось k раз. При n = 3 и k = 2 событие В 2 выражается через элементарные события А į (į = 1, 2, 3) по формуле:

В 2 = А 1 А 2 Ā 3 + А 1 Ā 2 А 3 + Ā 1 А 2 А 3 .

В общем виде последняя формула будет такой

т.е каждый член суммы (2.2) соответствует появлению события А k раз и (n – k ) раз непоявлений. Число всех комбинаций (слагаемых) в (2.2) равно числу способов выбора из n испытаний k испытаний, в которых событие А произошло, т.е числу сочетаний C n k . Вероятность каждой такой комбинации по теореме умножения вероятностей независимых событий равна p k × q n – k , так как p (А į) = p , p (Ā į) = q , i = 1,2,…,n . Но комбинации в (2.2) являются несовместными событиями. Поэтому по теореме сложения вероятностей получим

Таким образом, имеет место формула Бернулли

P n (k) = C n k p k q n-k .

Что и требовалось доказать.

Замечание 1. Сформулированная выше теорема относится к случаю, когда в каждом испытании вероятность появления события А постоянна. Тогда для расчета вероятности P n (k ) справедлива формула Бернулли (2.1). Если же вероятности наступления события А в испытаниях α 1 , α 2 , … , α n разные, т.е. вероятности составляют значения p 1 , p 2 , … , p n , то тогда вместо (2.1) справедлива формула:

Замечание 6. Вероятность того, что в n опытах, проводящихся по схеме Бернулли, успех наступит от k 1 до k 2 раз , вычисляется по формулеP n (k )) для конкретных значений n и p . Так как аргумент k принимает лишь целые значения, график представляется в виде точек на плоскости (k , P n (k )). Для наглядности точки соединяются ломаной линией, и такой график называется полигоном распределения (рис.2.1). При p = 0,5, n = 6, как показано на рисунке 2.1, полигон симметричен относительно прямой x = np (если p близко к 0,5, то полигон близок к симметричному). При малых p полигон существенно асимметричен, и наивероятнейшими явля-ются частоты, близкие к нулю. На рисунке 2.2 изображен полигон распределения для p = 0,2 при числе испытаний n = 6. При больших p , близких к 1, наиболее вероятны максимальные значения. На рис. 2.3 показан полигон распределения, для p = 0,8 и n = 6.

Рис. 2.3.

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ОБРАЗОВАНИЮ

Государственное образовательное учреждение

высшего профессионального образования

«МАТИ»  РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. К.Э. ЦИОЛКОВСКОГО

Кафедра «Моделирование систем и информационные технологии»

Повторение испытаний. Схема бернулли

Методические указания к практическим занятиям

по дисциплине «Высшая математика»

Составители: Егорова Ю.Б.

Мамонов И.М.

Москва 2006 введение

Методические указания предназначены для студентов дневного и вечернего отделения факультета №14 специальностей 150601, 160301, 230102. Указания выделяют основные понятия темы, определяют последовательность изучения материала. Большое количество рассмотренных примеров помогает в практическом освоении темы. Методические указания служат методической основой для практических занятий и выполнения индивидуальных заданий.

СХЕМА БЕРНУЛЛИ. ФОРМУЛА БЕРНУЛЛИ

Схема Бернулли - схема повторных независимых испытаний, при которой какое-то событие А может многократно повторяться с постоянной вероятностью Р (А )= р .

Примеры испытаний, проводимых по схеме Бернулли: многократное подбрасывание монеты или игральной кости, изготовление партии деталей, стрельба по мишени и т.п.

Теорема. Если вероятность наступления события А в каждом испытании постоянна и равна р , то вероятность того, что событие А наступит m раз в n испытаниях (безразлично в какой последовательности), можно определить по формуле Бернулли:

где q = 1 – p .

ПРИМЕР 1. Вероятность того, что расход электроэнергии на протяжении одних суток не превысит установленной нормы, равна р= 0,75. Найти вероятность того, что в ближайшие 6 суток расход электроэнергии в течение 4 суток не превысит нормы.

РЕШЕНИЕ. Вероятность нормального расхода элек­троэнергии на протяжении каждых из 6 суток постоянна и равна р = 0,75. Следовательно, вероятность перерасхода электроэнергии в каждые сутки также постоянна и равна q = 1р = 1  0,75 = 0,25.

Искомая вероятность по формуле Бернулли равна:

ПРИМЕР 2. Стрелок производит по мишени три выстрела. Вероятность попадания в мишень при каждом выстреле равна р= 0,3. Найти вероятность того, что поражена: а) одна мишень; б) все три мишени; в) ни одной мишени; г) хотя бы одна мишень; д) менее двух мишеней.

РЕШЕНИЕ. Вероятность попадания в мишень при каждом выстреле постоянна и равна р =0,75. Следовательно, вероятность промаха равна q = 1 р = 1  0,3= 0,7. Общее число проведенных опытов n =3.

а) Вероятность поражения одной мишени при трех выстрелах равна:

б) Вероятность поражения всех трех мишеней при трех выстрелах равна:

в) Вероятность трех промахов при трех выстрелах равна:

г) Вероятность поражения хотя бы одной мишени при трех выстрелах равна:

д) Вероятность поражения менее двух мишеней, то есть или одной мишени, или ни одной:

1. Локальная и интегральная теоремы муавра-лапласа

Если произведено большое число испытаний, то вычисление вероятностей по формуле Бернулли становится технически сложным, так как формула требует действий над огромными числами. Поэтому существуют более простые приближенные формулы для вычисления вероятностей при больших n . Эти формулы называются асимптотическими и определяются теоремой Пуассона, локальной и интегральной теоремой Лапласа.

Локальная теорема Муавра-Лапласа. А А произойдет m раз в n n (n →∞ ), приближенно равна:

где функция
а аргумент

Чем больше n , тем точнее вычисление вероятностей. Поэтому теорему Муавра-Лапласа целесообразно применять при npq  20.

f ( x ) составлены специальные таблицы (см. приложение 1). При использовании таблицы необходимо иметь в виду свойства функции f(x) :

Функция f(x) является четной f(  x)= f(x) .

При х  ∞ функция f(x)  0. Практически можно считать, что уже при х >4 функция f(x) ≈0.

ПРИМЕР 3. Найти вероятность того, что событие А наступит 80 раз в 400 испытаниях, если вероятность появления события А в каждом испытании равна р= 0,2.

РЕШЕНИЕ. По условию n =400, m =80, p =0,2, q =0,8. Следовательно:

По таблице определим значение функции f (0)=0,3989.

Интегральная теорема Муавра-Лапласа. Если вероятность наступления события А в каждом испытании постоянна и отлична от 0 и 1, то вероятность того, что событие А произойдет от m 1 до m 2 раз в n испытаниях при достаточно большом числе n (n →∞ ), приближенно равна:

где
 интеграл или функция Лапласа,

Для нахождения значений функции Ф( x ) составлены специальные таблицы (например, см. приложение 2). При использовании таблицы необходимо иметь в виду свойства функции Лапласа Ф(x) :

Функция Ф(x) является нечетной Ф(  x)=  Ф(x) .

При х  ∞ функция Ф(x)  0,5. Практически можно считать, что уже при х >5 функция Ф(x) ≈0,5.

Ф (0)=0.

ПРИМЕР 4. Вероятность того, что деталь не прошла проверку ОТК, равна 0,2. Найти вероятность того, что среди 400 деталей окажется непроверенных от 70 до 100 деталей.

РЕШЕНИЕ. По условию n =400, m 1 =70, m 2 =100, p =0,2, q =0,8. Следовательно:

По таблице, в которой приведены значения функции Лапласа, определяем:

Ф(x 1 ) = Ф(  1,25 )=  Ф( 1,25 )=  0,3944; Ф(x 2 ) = Ф( 2,5 )= 0,4938.

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

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

Пусть теперь производится n независимых опытов, в каждом из которых с одной и той же вероятностью p может наступить некоторое событие А . Требуется найти вероятность Р n (к) того, что в n опытах событие А наступит ровно к раз (событие В ).

Описанная схема называется схемой независимых испытаний, или схемой Бернулли, по имени швейцарского математика конца XVII и начала XVIII века Якоба Бернулли, изучавшего её.

Найдем вероятность Р n (к) . Событие В можно представить в виде суммы ряда элементарных событий – вариантов события А . Каждый вариант события А можно записать в виде строки длиной n (число опытов), в которой к компонент соответствуют событию А , а остальные n-к компонент событию . Например, один из возможных вариантов есть

(успех и 1,2,…,k -м опытах и неудача в остальных).

Число всех вариантов равно (числу сочетаний из n элементов по к ), а вероятность каждого варианта в виду независимости опытов равна р к q n -к (где q =1-р ). Отсюда вероятность события В будет равна

Формула (1) носит название формулы Бернулли .

Отсюда следует, что вероятность, хотя бы одного появления события А при n независимых испытаниях (опытах) в одинаковых условиях равна

Пример 1 . Монета бросается 5 раз. Какова вероятность того, что герб выпадет при этом 3 раза?

Решение . В данном случае событием А считается выпадение герба, вероятность p этого события в каждом опыте равна . Отсюда

P = .

Для наглядности условимся каждое наступление события А рассматривать как успех. Если зафиксировать n , то, Р n (к) . есть функция аргумента к , принимающая значения . Выясним, при каком значении к функция Р n (к) принимает наибольшее значение, т.е., какое число успехов к 0 является наиболее вероятным при данном числе опытов n . Оказывается что число к=к 0 можно определить из двойного неравенства.

(3)

Разность граничных значений в этом двойном неравенстве равна 1. Если np+p не является целым числом, то двойное неравенство определяет лишь одно наивероятнейшее значение к 0 .Если же np+p – целое число, то имеются два наивероятнейших значения: и .

Пример 2 . Игральную кость бросают 20 раз. Каково наиболее вероятное число выпадений грани «6» ?

Решение. В данном случае n = 20, откуда . Поскольку nр + р не целое число, то наибольшим среди чисел Р 20 (0), Р 20 (1),…, Р 20 (20) будет число Р 20 (3). Следовательно, наиболее вероятное число выпадений грани «6» будет 3. Найдём, чему равна вероятность такого числа выпадений. По формуле Бернулли имеет:

.

Из формулы (3) видно, что одно из двух ближайших к nр целых чисел является наиболее вероятным числом успехов.

Оказывается, число nр допускает и другую интерпретацию. А именно: nр можно рассматривать, в определенном смысле, как среднее число успехов в n опытах . Будем исходить из частотного истолкования вероятности. Назовем (для краткости) n - кратное повторение данного опыта серией. Пусть мы произвели N серий. Пусть в первой серии было получено к 1 успехов, во второй – к 2 , ….., в N -ой –к N . Составим среднее арифметическое этих чисел

. (4)

Равенство (4) - есть среднее число успехов в N сериях. Оказывается, что с увеличением N указанное среднее арифметическое приближается к некоторому постоянному значению, а именно к числу np .

Действительно запишем (4) в виде:

. (5)

Поскольку каждая серия состоит из n опытов, то производя N серий мы осуществляем данный опыт раз.

Написанная дробь (5) со знаменателем Nn есть нечто иное как отношение общего числа успехов в этих опытах к числу всех опытов. С увеличением N (а значит, и Nn ) эта дробь будет приближаться к числу р - вероятности успеха. Следовательно, число (4) будет приближаться к рn , что и требовалось получить.

Пример 3 . Станок штампует изделия. Вероятность р брака одного изделия равна 0,05. Чему равно среднее число бракованных изделий на сотню?

Решение . Искомое число бракованных изделий равно: .

Замечание 1. Можно рассмотреть более общую схему независимых испытаний. Рассмотрим n независимых испытаний (в различных условиях), причём вероятность события А («успеха») в i -ом опыте равна p i , a q i =1-p i – вероятность неуспеха в i -м испытании (i =1,2,…,n ). Тогда можно показать, что вероятность P n (к) того, что событие А появится в этих n опытах ровно к раз, равна коэффициенту при z k в разложении по степеням z функции

Такую схему независимых испытаний называют схемой Пуассона . Схема Пуассона при p i =p превращается в схему Бернулли. Вероятности P n (к) в схеме Пуассона не записываются в компактном виде аналогичной формуле(1). Из (6) , например, следует:

Замечание 2. Схемы Бернулли и Пуассона допускают обобщение на тот случай, когда в результате каждого опыта возможные не два исхода (А или ), а несколько исходов.

Если производится n независимых опытов (схема Бернулли) причём каждый опыт может иметь к исключающих друг друга исходов , с вероятностями , то вероятность того, что в m 1 опытах появится событие А 1 , в m 2 опытах событие А 2 и т.д., в m k опытах событие А к выражается формулой

Если условия опыта различны (схема Пуассона), т.е.

в i- омопыте событие A j имеет вероятность p ji (i =1,2,…,n ; j =1,2,…,k ), то вероятность вычисляется как коэффициент при члене в разложении по степеням функции:

Пример 4. Завод изготавливает изделия, каждое из которых подвергается четырём видам испытаний. Первое испытание изделия проходит благополучно с вероятностью 0,9; второе с вероятностью 0,95; третье-0,8 и четвертое-0,85. Найти вероятность того, что изделие пройдет благополучно:

A- все четыре испытания

B- ровно два испытания (из четырех)

C- не менее двух испытании (из четырех)

Решение. В условиях задачи проводятся четыре независимых опыта (испытания) в различных условиях. Вероятность события. А – испытание прошло благополучно, в каждом опыте разное. Искомые вероятности находим из формулы (6)

Отсюда получаем:

§12. Вероятности P n (к) при больших значениях n . Приближённые формулы Лапласа и Пуассона.

В приложениях часто возникает необходимость в вычислении вероятностей Р n (к) для весьма больших значений n и k . Рассмотрим, например, такую задачу.

Задача. На некотором предприятии вероятность брака, равна 0,02. Обследуются 500 изделий готовой продукции. Найти вероятность того, что среди них окажется ровно 10 бракованных.

Рассматривая обследование каждого изделия как отдельный опыт, можно сказать, что производиться 500 независимых опытов, причем в каждом их них событие А (изделие оказалось бракованным) наступает с вероятностью 0,02, тогда по формуле Бернулли получаем

Непосредственный подсчет этого выражения представляется сложным. Ещё большую трудность пришлось бы испытать, если бы мы искали вероятность того, что число бракованных изделий среди 500 окажется в пределах, скажем, от 10 до 20. В этом случае потребовалось бы вычислить сумму , что является более сложным делом.

Задачи подобного рода встречаются в приложениях весьма часто. Поэтому возникает необходимость в отыскании приближённых формул для вероятностей Р n (к) , а также для сумм вида

(1)

при больших n .

1. Приближённые формулы Лапласа. Их используют при больших n (порядка сотен или тысяч), вероятностей p или q не слишком близким к 0 или 1 (порядка сотых долей). Обычно условием применения этих приближений является условие npq >9.

а) Локальная приближённая формула Лапласа . При больших n справедливо равенство.

, (2)

где , а φ (х ) обозначает следующую функцию: .

Заметим, что функция φ(х) табулирована, т.е. для нее составлена таблица её значений.

Вторая приближённая формула Лапласа даёт приближённые значения для величины -вероятности того, что число наступлений события А в n опытах (число «успехов») окажется заключенным между заданными границами к 1 и к 2 .

б) Интегральная приближённая формула Лапласа . При больших n справедливо приближённое равенство

, (3)

где Φ(х) обозначает следующую функцию

. (4)

Функция Φ(х) обладает следующими полезными для вычисления свойствами:

1. Φ(х) – нечётная функция: ,

2. при возрастании х от 0 до ∞ функция Φ(х) растет от 0

до 0,5, причем уже при х = 5 значение функции Φ(х)

отличается от 0,5 меньше чем на (т.е. при функция Ф(x) практически равна 0,5).

Пример 1. Монету бросают 100 раз. Какова вероятность того, что герб выпадет ровно 50 раз?

Решение . Имеем: npq = 100· · = 25>9. Воспользовавшись приближённой формулой (2), получим. . Из таблицы для функции φ(x) найдем, что φ(0) = 0,3989…. Отсюда получаем .

Пример 2 . Доведём до конца решение задачи, приведённой в начале этого параграфа. В ней требовалось найти , а также вероятность P 500 (10≤ к ≤20).

Решение. В данном случае npq = 500·0,02·0,98=9,8. Воспользовавшись приближёнными формулами (2) и (3), получим: ,

Замечание. Если мы осуществляем опыт n раз и k - число наступлений события А при этом, то, вообще говоря, дробь -относительная частота наступления события

А – будет близка к р (вероятности события А ). Однако сколь тесной окажется эта близость, предугадать невозможно.

Интегральная теорема Лапласа позволяет оценить вероятность неравенства при достаточно больших n и значениях р не слишком близких к 0 или 1, т.е. определить вероятность того, что отклонение частоты случайного события от его вероятности р по абсолютной величине не превосходит некоторого . Имеем

Таким образом, получаем

(5)

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

На практике надёжность оценки задаётся заранее. Тогда по заданной надёжности можно найти соответствующее значение из уравнения с помощью таблиц функции Лапласа. В этом случае доверительная оценка с заданной надёжностью примет вид р или q к нулю, поэтому, в этом случае используют приближённые формулы Пуассона. При больших n (порядка тысяч, десятков тысяч и больше) и малых р (порядка тысячных долей и меньше) справедливы приближённые равенства. Обычно условием применения этих приближений является условие npq <9.

, (7).

, (8)

где λ =np .

Особенностью формул (7) и (8) является то, что для того, чтобы найти вероятность того или иного числа успехов, вовсе не требуется знать n и р . Всё определяется числом λ=np , которое является (см. §11) средним числом успехов .

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

Пример 5 . Прядильщица обслуживает 1000 веретён. Вероятность обрыва нити на одном веретене в течение одной минуты равна 0,004. Найти вероятность того, что в течение одной минуты обрыв произойдет в пяти веретенах.

Решение. Формула Бернулли приведёт к громоздким вычислениям, поэтому воспользуемся формулой Пуассона (7). Здесь к = 5, р =0.004, n = 1000, тогда λ = np = 4.

Отсюда: .

Пример 6 . Книга в 1000 страниц имеет 100 опечаток. Какова вероятность того, что на случайно выбранной странице будет не менее четырёх опечаток (событие В ).

Решение: Среднее количество опечаток на одну страницу есть . В данном случае следует применить формулу Пуассона. Тогда вероятность p к иметь к опечаток на одной странице будет равна .

Сумма р = p 0 +p 1 +p 2 +p 3 есть вероятность того, что на странице окажется не более трёх опечаток. Пользуясь таблицами (или калькулятором) получаем р = 0,999996 (в данном случае мы пользовались калькулятором, таблицы дадут р =0,9048+0,0905+0,0045+0,0002=1). Вероятность того, что на случайно выбранной странице будет не менее четырёх опечаток, равна 1-р =1-0,999996=0,0000004 (таблицы дадут 1-р =1-1=0). Отсюда можно сделать вывод, что событие В практически невозможно.