Raw

Лекція 9. Основи комбінаторики та теорії ймовірностей

Огляд

Усі розділи дискретної математики, які ми проходили досі, описують детермінований світ: множина або містить елемент, або ні (Лекція 1); два об’єкти або перебувають у відношенні, або ні; висловлення або істинне, або хибне. Проте більшість реальних задач фаху звучить інакше: «наскільки ймовірно, що навмання взятий відбиток виявиться бракованим?», «як часто друкарська машина дає збій за зміну?», «якою мірою можна довіряти партії, якщо у вибірці з тиражу трапилося два дефекти?». Тут потрібна математика градуйованої певності — теорія ймовірностей, яка кожному непевному твердженню зіставляє число від 00 до 11, що вимірює, наскільки воно правдоподібне.

Ця лекція будує скінченну (дискретну) ймовірність з трьох цеглин.

  1. Мова. Дослід, простір елементарних наслідків Ω\Omega і події — які виявляться не чим іншим, як множинами з Лекції 1, тож уся алгебра множин переноситься сюди дослівно, лише змінюється прочитання.
  2. Міра. Як зіставити події число P(A)[0,1]P(A) \in [0,1]: класичне (лапласове) означення через відношення m/nm/n і комбінаторні засоби, потрібні, щоб полічити ці mm і nn.
  3. Числення подій. Теореми, що комбінують ймовірності, — правило протилежної події, умовна ймовірність, теореми множення й додавання, формула повної ймовірності та формула Байєса, — кожну з яких ми доводимо, а не лише проголошуємо.

Наприкінці ми вивчаємо випадкові величини — математичне сподівання M(X)M(X), дисперсію D(X)D(X) і середньоквадратичне відхилення σ\sigma — та оглядово перелічуємо типові розподіли. Наскрізна прикладна нитка — контроль якості друку та вибіркові перевірки тиражу: саме ймовірнісні моделі стоять за стандартами приймального контролю, за оцінкою надійності обладнання та за діагностикою джерела дефекту.

Про строгість. Кілька результатів нижче оформлено як теореми й супроводжено доведенням у межах позначок Доведення.\blacksquare. Уся теорія ймовірностей тримається на кількох простих властивостях, а решта — їхні наслідки. Читайте доведення активно: майже кожне з них — це спочатку дія над множинами (з Лекції 1), а потім застосування властивості ймовірності. Коли словесна умова заплутує, перекладіть її мовою множин над Ω\Omega — арифметика піде слідом.


9.1 Випадкові події та простір елементарних наслідків

Дослід, наслідок, простір

Випадковим дослідом (випробуванням) називають процес, окремий результат якого не можна передбачити наперед із певністю, хоча множина всіх можливих результатів відома, а сам дослід можна (принаймні подумки) повторювати за однакових умов. Класичні приклади — підкидання монети, кидання грального кубика, витягання кулі з ящика, витягання лотерейного квитка; так само влаштована й «вибірка одного відбитка з тиражу».

Означення (простір елементарних наслідків). Простором елементарних наслідків Ω\Omega (позначають також UU чи SS) називають множину всіх можливих результатів досліду. Кожен окремий результат ωΩ\omega \in \Omega — це елементарний наслідок (елементарна подія), який уже не розкладається на дрібніші.

Означення (подія). Подія — це будь-яке твердження про дослід, що стає однозначно істинним або хибним, щойно відомий наслідок; формально це підмножина AΩA \subseteq \Omega. Кажуть, що подія AA настала, якщо здійснений наслідок ω\omega належить AA. Подію, що містить більш ніж один елементарний наслідок, називають складеною.

Отже, згідно зі слайдами, «подія — будь-який факт, який може відбутися чи не відбутися», а мовою множин це просто підмножина простору Ω\Omega.

Приклад 9.1 (кубик). Кинемо один симетричний гральний кубик: Ω={1,2,3,4,5,6}\Omega = \{1,2,3,4,5,6\}, Ω=6|\Omega| = 6. Елементарна подія «випало 44» — це синглетон {4}\{4\}. Складена подія A=A = «випало парне число» — це {2,4,6}\{2,4,6\}; подія B=B = «випало більше за 44» — це {5,6}\{5,6\}. Тоді AB={2,4,5,6}A \cup B = \{2,4,5,6\}, AB={6}A \cap B = \{6\}, а A={1,3,5}\overline{A} = \{1,3,5\}.

Приклад 9.2 (дві монети). Підкинемо дві розрізнювані монети: Ω={ГГ,ГР,РГ,РР}\Omega = \{ГГ, ГР, РГ, РР\} (Г — герб, Р — решка), Ω=4|\Omega| = 4. Подія «рівно один герб» — це {ГР,РГ}\{ГР, РГ\}. Зауважте, що цей простір є декартовим добутком {Г,Р}×{Г,Р}\{Г,Р\} \times \{Г,Р\} з Лекції 1: складені простори будують із простих саме операцією добутку.

Простір елементарних наслідків Ω як точки в прямокутнику; подія A ⊆ Ω виділена заштрихованою овальною областю навколо сприятливих наслідків

Події — це множини: алгебра подій

Оскільки подія є підмножиною Ω\Omega, усі операції над множинами з Лекції 1 діють без жодних змін — змінюється лише прочитання. Ця таблиця-словник — найкорисніша в усій лекції; тримайте її перед очима.

Мова множин Мова подій «настає, коли…»
ωA\omega \in A наслідок ω\omega спричиняє подію AA здійснений наслідок лежить у AA
Ω\Omega достовірна подія завжди
\varnothing неможлива подія ніколи
ABA \cup B «AA або BB» (принаймні одна) ωA\omega \in A або ωB\omega \in B
ABA \cap B (також ABAB) «AA і BB» (обидві) ωA\omega \in A і ωB\omega \in B
A=ΩA\overline{A} = \Omega \setminus A «не AA» (протилежна) ωA\omega \notin A
ABA \subseteq B «AA спричиняє BB» коли настає AA, настає й BB
AB=A \cap B = \varnothing A,BA, B несумісні не можуть настати разом

Тому закони де Моргана з Лекції 1 у мові подій читаються так:

AB=AB(«не (A або B=«ані те, ані те»),AB=AB.\overline{A \cup B} = \overline{A} \cap \overline{B} \quad (\text{«не (}A\text{ або }B\text{)»} = \text{«ані те, ані те»}), \qquad \overline{A \cap B} = \overline{A} \cup \overline{B}.

Зауваження (чому «події — це множини» важливо). Це не гарний збіг, який варто відзначити й забути. Кожне перетворення ймовірностей нижче — спершу перетворення множин, і лише потім застосування властивості міри. Тому навичка з Лекції 1 — вправно об’єднувати, перетинати й доповнювати множини — тут працює безпосередньо.


9.2 Типи подій

Те, як події пов’язані одна з одною, вирішує, яке правило застосовувати далі, тож класифікацію варто сформулювати точно (слайд 2).

Достовірна подія — настає в кожному випробуванні. Це весь простір Ω\Omega; згодом доведемо, що її ймовірність дорівнює 11. Приклад: «кубик показав число з {1,,6}\{1,\dots,6\}».

Неможлива подія — не настає ніколи. Це порожня множина \varnothing; її ймовірність дорівнює 00. Приклад: «кубик показав 77».

Сумісні подіїможуть настати в одному випробуванні: ABA \cap B \neq \varnothing. На одному кубику A=A = «парне» і B=B = «більше за 33» мають спільний наслідок 66, тож вони сумісні.

Несумісні подіїне можуть настати разом: AB=A \cap B = \varnothing (диз’юнктні множини з Лекції 1). «Випало 22» і «випало 55» не можуть статися одночасно за одного кидка. Різні елементарні наслідки завжди попарно несумісні.

Означення (повна група подій). Події A1,,AnA_1, \dots, A_n утворюють повну групурозбиттям простору Ω\Omega), якщо вони попарно несумісні (AiAj=A_i \cap A_j = \varnothing при iji \ne j) і в об’єднанні дають увесь простір (A1An=ΩA_1 \cup \dots \cup A_n = \Omega). Тоді в кожному випробуванні настає рівно одна з них.

Протилежні події — окремий випадок повної групи з двох подій: AA і A\overline{A} несумісні й разом вичерпують Ω\Omega. Приклад: «кубик показав 66» і «кубик не показав 66». Далі доведемо, що P(A)+P(A)=1P(A) + P(\overline{A}) = 1.

Незалежні події — настання однієї не змінює ймовірності другої (точне означення — у §9.4). Витягли кулю, повернули її назад, витягли знову — це незалежні випробування.

Залежні події — настання однієї змінює ймовірність другої. Витягання двох куль без повернення робить другий добір залежним від першого.

Три типи взаємного розташування подій: несумісні (кола не перетинаються), сумісні (кола перетинаються, спільна частина заштрихована) та протилежні (подія A і її доповнення в межах Ω)

Типова помилка (незалежність — це не несумісність). Ці поняття майже протилежні. Далі ми доведемо, що дві несумісні події з додатними ймовірностями обов’язково залежні: якщо настала AA, то BB вже неможлива — а це максимальна можлива зміна ймовірності BB. Не ототожнюйте «AB=A \cap B = \varnothing» з «A,BA, B незалежні».


9.3 Класичне означення ймовірності

Коли скінченний дослід має nn рівноможливих елементарних наслідків, з яких рівно mm сприятливі для події AA, класичну (лапласову) ймовірність означують як відношення (слайд 3).

Означення (класична ймовірність). Якщо всі елементарні наслідки рівноможливі, то

P(A)=mn=AΩ=кількість сприятливих наслідківкількість усіх наслідків.P(A) = \frac{m}{n} = \frac{|A|}{|\Omega|} = \frac{\text{кількість сприятливих наслідків}}{\text{кількість усіх наслідків}}.

Гіпотеза «рівноможливості» істотна й ґрунтується на симетрії: у симетричної монети, правильного кубика, добре перетасованої колоди чи ретельно перемішаного ящика немає фізичної причини надавати перевагу якомусь наслідку. Там, де симетрії немає (погнута монета, невідомий технологічний процес), класичне означення не працює, і ми звертаємось до статистичної ймовірності (§9.10).

Приклад 9.3 (кубик). P(парне)=36=12P(\text{парне}) = \tfrac{3}{6} = \tfrac12, бо парними є 33 з 66 рівноможливих граней.

Приклад 9.4 (карти). Зі стандартної колоди на 5252 карти P(туз)=452=113P(\text{туз}) = \tfrac{4}{52} = \tfrac{1}{13}, а P(фігурна карта)=1252=313P(\text{фігурна карта}) = \tfrac{12}{52} = \tfrac{3}{13} (валети, дами, королі).

Основні властивості ймовірності

Оскільки завжди 0mn0 \le m \le n, з означення миттєво випливають три властивості.

Теорема 9.1 (властивості класичної ймовірності). Для будь-якої події AA:

0P(A)1,P(Ω)=1,P()=0.0 \le P(A) \le 1, \qquad P(\Omega) = 1, \qquad P(\varnothing) = 0.

Доведення. Кількість сприятливих наслідків задовольняє 0mn0 \le m \le n, тож 0mn10 \le \tfrac{m}{n} \le 1. Для достовірної події сприятливі всі наслідки, m=nm = n, звідки P(Ω)=nn=1P(\Omega) = \tfrac{n}{n} = 1. Для неможливої події немає жодного сприятливого наслідку, m=0m = 0, звідки P()=0n=0P(\varnothing) = \tfrac{0}{n} = 0. \blacksquare

Отже, «ймовірність ніколи не буває вищою за 11» (слайд 4) — це не окрема аксіома, а простий наслідок означення.

Протилежна подія

Теорема 9.2 (правило протилежної події). Для будь-якої події AA

P(A)=1P(A).P(\overline{A}) = 1 - P(A).

Доведення. Події AA і A\overline{A} несумісні (AA=A \cap \overline{A} = \varnothing) і разом вичерпують простір (AA=ΩA \cup \overline{A} = \Omega). У класичній моделі сприятливі для AA наслідки та сприятливі для A\overline{A} наслідки не перетинаються, а разом дають усі nn наслідків, тож mA+mA=nm_A + m_{\overline{A}} = n. Поділивши на nn, дістаємо P(A)+P(A)=1P(A) + P(\overline{A}) = 1, звідки P(A)=1P(A)P(\overline{A}) = 1 - P(A). \blacksquare

Це — робочий кінь усіх обчислень «принаймні один», бо доповнення «жодного» майже завжди рахувати легше. Ми скористаємося ним у Задачі 4 (§9.8): «виграш хоча б по одному квитку» простіше знайти через «жодного виграшу».

Приклад 9.5. Кинули два кубики. Замість того, щоб рахувати всі способи дістати принаймні одну шістку, знайдемо ймовірність протилежної події «жодної шістки»: кожен кубик уникає шістки з ймовірністю 56\tfrac56, тож P(жодної шістки)=5656=2536P(\text{жодної шістки}) = \tfrac56 \cdot \tfrac56 = \tfrac{25}{36} (незалежність кубиків — §9.5). Тоді P(принаймні одна шістка)=12536=11360,306P(\text{принаймні одна шістка}) = 1 - \tfrac{25}{36} = \tfrac{11}{36} \approx 0{,}306.


9.4 Умовна ймовірність

Часто ми дізнаємось часткову інформацію — відомо, що настала пов’язана подія BB, — і маємо переглянути ймовірність AA з огляду на це.

Означення (умовна ймовірність). Для подій A,BA, B з P(B)>0P(B) > 0 умовною ймовірністю AA за умови BB називають

P(AB)=P(AB)P(B).P(A \mid B) = \frac{P(A \cap B)}{P(B)}.

Інтуїція. Звістка «BB настала» звужує світ з усього Ω\Omega до самого BB: наслідки поза BB стали неможливими. Ми перемасштабовуємо так, щоб BB несла повну ймовірність 11 (ділимо на P(B)P(B)), і питаємо, яка частка ймовірності BB лежить також у AA (чисельник P(AB)P(A \cap B)). У класичній моделі це видно прямо:

P(AB)=AB/ΩB/Ω=ABBP(A \mid B) = \frac{|A \cap B| / |\Omega|}{|B| / |\Omega|} = \frac{|A \cap B|}{|B|}

— буквально «сприятливі-всередині-BB поділити на розмір BB».

Діаграма Венна двох подій A і B, що перетинаються; спільна частина A∩B заштрихована; умовна ймовірність P(A|B) дорівнює частці цієї спільної частини у події B

Приклад 9.6 (кубик). Нехай A=A = «випало 22», B=B = «випало парне». Тоді

P(AB)=P({2})P({2,4,6})=1/63/6=13.P(A \mid B) = \frac{P(\{2\})}{P(\{2,4,6\})} = \frac{1/6}{3/6} = \frac13.

Звістка «число парне» піднімає шанс двійки з 16\tfrac16 до 13\tfrac13, бо лишилося лише три можливі наслідки.

Незалежні події

Якщо звістка про BB нічого не змінює в ймовірності AA, події називають незалежними. Це формалізує §9.2.

Означення (незалежність). Події AA і BB незалежні, якщо

P(AB)=P(A)P(B).P(A \cap B) = P(A)\,P(B).

Рівносильно (коли P(B)>0P(B) > 0) це означає P(AB)=P(A)P(A \mid B) = P(A): знання про BB не додає інформації про AA. Для залежних подій, навпаки, P(AB)P(A)P(A \mid B) \ne P(A), і в обчисленнях слід брати саме умовну ймовірність.

Отже, як зазначено на слайді 5, умовна ймовірність по-різному входить у формули для залежних і незалежних подій: у першому випадку P(AB)P(A \mid B) треба обчислювати з умови задачі, у другому вона просто дорівнює P(A)P(A).

Типова помилка (переставляння умови). P(AB)P(A \mid B) і P(BA)P(B \mid A) — різні числа; їх ототожнення називають помилкою прокурора. «Імовірність такого збігу ДНК за умови невинності — один на мільйон» — це P(збігневинний)P(\text{збіг} \mid \text{невинний}), а не P(невиннийзбіг)P(\text{невинний} \mid \text{збіг}). Формула Байєса (§9.6) — це якраз машина для перетворення однієї з них на іншу.


9.5 Теорема множення ймовірностей

Переписавши означення умовної ймовірності без дробу, дістаємо ймовірність того, що настануть обидві події (слайд 6).

Теорема 9.3 (теорема множення). Для подій A,BA, B з додатними ймовірностями

P(AB)=P(B)P(AB)=P(A)P(BA).P(A \cap B) = P(B)\,P(A \mid B) = P(A)\,P(B \mid A).

Якщо події незалежні, то P(BA)=P(B)P(B \mid A) = P(B), і формула спрощується до

P(AB)=P(A)P(B).P(A \cap B) = P(A)\,P(B).

Доведення. За означенням умовної ймовірності P(AB)=P(AB)P(B)P(A \mid B) = \dfrac{P(A \cap B)}{P(B)}; помноживши обидві частини на P(B)P(B), дістаємо P(AB)=P(B)P(AB)P(A \cap B) = P(B)\,P(A \mid B). Симетрично, міняючи ролями AA і BB, маємо P(AB)=P(A)P(BA)P(A \cap B) = P(A)\,P(B \mid A). Для незалежних подій P(BA)=P(B)P(B \mid A) = P(B) за означенням, звідки P(AB)=P(A)P(B)P(A \cap B) = P(A)\,P(B). \blacksquare

Наочна схема — дерево ймовірностей. Теорему множення зручно читати з дерева: кожна гілка позначена своєю (умовною) ймовірністю, а ймовірність шляху дорівнює добутку ймовірностей уздовж нього. Наприклад, шлях «спершу AA, потім BB» має ймовірність P(A)P(BA)=P(AB)P(A) \cdot P(B \mid A) = P(A \cap B).

Дерево ймовірностей для теореми множення: із кореня дві гілки A і не-A з ймовірностями P(A) та P(не-A); від кожної — гілки B і не-B з умовними ймовірностями; ймовірність шляху дорівнює добутку ймовірностей його ребер

Приклад 9.7 (незалежні кидки). Кинули кубик двічі. Оскільки кидки незалежні, P(перша 6 і друга 6)=1616=136P(\text{перша } 6 \text{ і друга } 6) = \tfrac16 \cdot \tfrac16 = \tfrac{1}{36}, що збігається з прямим переліком 3636 рівноможливих пар.

Загальне правило множення (ланцюгове)

Для більшого числа подій теорема множення узагальнюється на ланцюгове правило:

P(A1A2An)=P(A1)P(A2A1)P(A3A1A2)P(AnA1An1).P(A_1 \cap A_2 \cap \dots \cap A_n) = P(A_1)\,P(A_2 \mid A_1)\,P(A_3 \mid A_1 \cap A_2) \cdots P(A_n \mid A_1 \cap \dots \cap A_{n-1}).

Його доводять індукцією за nn: база n=2n = 2 — це Теорема 9.3, а крок отримують, застосувавши теорему множення до пари (A1An1, An)\big( A_1 \cap \dots \cap A_{n-1},\ A_n \big). Саме ланцюгове правило працює в задачах на послідовне витягання без повернення (Задача 2, §9.8).


9.6 Формула повної ймовірності та формула Байєса

Нехай перша стадія досліду породжує одну з кількох гіпотез (причин, сценаріїв), а пізніша подія AA має відому ймовірність за кожної гіпотези. Формула повної ймовірності складає ці частини в P(A)P(A).

Теорема 9.4 (формула повної ймовірності). Нехай гіпотези H1,,HnH_1, \dots, H_n утворюють повну групу (попарно несумісні, H1Hn=ΩH_1 \cup \dots \cup H_n = \Omega, кожна P(Hi)>0P(H_i) > 0). Тоді для будь-якої події AA

P(A)=i=1nP(Hi)P(AHi).P(A) = \sum_{i=1}^{n} P(H_i)\,P(A \mid H_i).

Доведення. Оскільки гіпотези покривають увесь простір, за дистрибутивністю (Лекція 1)

A=AΩ=A(i=1nHi)=i=1n(AHi).A = A \cap \Omega = A \cap \Big( \bigcup_{i=1}^{n} H_i \Big) = \bigcup_{i=1}^{n} (A \cap H_i).

Події AHiA \cap H_i попарно несумісні, бо (AHi)(AHj)HiHj=(A \cap H_i) \cap (A \cap H_j) \subseteq H_i \cap H_j = \varnothing при iji \ne j. Тому за властивістю адитивності (сума ймовірностей несумісних подій) і теоремою множення

P(A)=i=1nP(AHi)=i=1nP(Hi)P(AHi).P(A) = \sum_{i=1}^{n} P(A \cap H_i) = \sum_{i=1}^{n} P(H_i)\,P(A \mid H_i). \qquad \blacksquare

Читають це як зважене середнє: імовірність AA — це середнє її умовних імовірностей P(AHi)P(A \mid H_i), кожну з яких зважено на правдоподібність відповідної гіпотези P(Hi)P(H_i).

Дерево гіпотез для формул повної ймовірності та Байєса: із кореня три гілки-гіпотези H1, H2, H3 з апріорними ймовірностями, від кожної — гілка до події A з умовною ймовірністю P(A|Hi)

Формула Байєса

Формула повної ймовірності працює «вперед»: від причин HiH_i до наслідку AA. Формула Байєса (слайд 7) обертає її назад — переглядає ймовірність кожної причини після того, як спостережено наслідок. Це — математика навчання на досвіді.

Теорема 9.5 (формула Байєса). Нехай H1,,HnH_1, \dots, H_n — повна група з P(Hi)>0P(H_i) > 0, а AA — подія з P(A)>0P(A) > 0. Тоді для кожного ii

P(HiA)=P(Hi)P(AHi)j=1nP(Hj)P(AHj)=P(Hi)P(AHi)P(A).P(H_i \mid A) = \frac{P(H_i)\,P(A \mid H_i)}{\displaystyle\sum_{j=1}^{n} P(H_j)\,P(A \mid H_j)} = \frac{P(H_i)\,P(A \mid H_i)}{P(A)}.

Доведення. За означенням умовної ймовірності (бо P(A)>0P(A) > 0) і теоремою множення в чисельнику

P(HiA)=P(HiA)P(A)=P(Hi)P(AHi)P(A).P(H_i \mid A) = \frac{P(H_i \cap A)}{P(A)} = \frac{P(H_i)\,P(A \mid H_i)}{P(A)}.

Нарешті підставляємо в знаменник формулу повної ймовірності P(A)=jP(Hj)P(AHj)P(A) = \sum_{j} P(H_j)\,P(A \mid H_j). \blacksquare

Словник. P(Hi)P(H_i) називають апріорною ймовірністю (переконання до спостереження), P(AHi)P(A \mid H_i)правдоподібністю (наскільки гіпотеза ii передбачає спостережене), а P(HiA)P(H_i \mid A)апостеріорною ймовірністю (переглянуте переконання). Формула Байєса — це правило перетворення «апріорі ++ дані» на «апостеріорі».

Приклад 9.8 (яка машина дала брак — повна ймовірність і Байєс). Тираж друкують три машини. Машина I дає 50%50\% накладу, машина II — 30%30\%, машина III — 20%20\%. Частка дефектних відбитків у них становить відповідно 2%2\%, 3%3\% і 5%5\%. Гіпотези H1,H2,H3H_1, H_2, H_3 — «відбиток надрукувала машина I / II / III»; подія AA — «відбиток дефектний».

Крок 1 — повна ймовірність дефекту.

P(A)=0,50,02+0,30,03+0,20,05=0,010+0,009+0,010=0,029.P(A) = 0{,}5 \cdot 0{,}02 + 0{,}3 \cdot 0{,}03 + 0{,}2 \cdot 0{,}05 = 0{,}010 + 0{,}009 + 0{,}010 = 0{,}029.

Тобто загалом бракованих близько 2,9%2{,}9\% тиражу.

Крок 2 — Байєс: навмання взятий дефектний відбиток, яка машина найімовірніше винна?

P(H1A)=0,0100,0290,345,P(H2A)=0,0090,0290,310,P(H3A)=0,0100,0290,345.P(H_1 \mid A) = \frac{0{,}010}{0{,}029} \approx 0{,}345, \quad P(H_2 \mid A) = \frac{0{,}009}{0{,}029} \approx 0{,}310, \quad P(H_3 \mid A) = \frac{0{,}010}{0{,}029} \approx 0{,}345.

Хоча машина III дає найбільшу частку браку (5%5\%), вона друкує мало (20%20\%), тож за знайденим дефектом машини I і III однаково ймовірні як джерело. Це типовий сюжет діагностики: апостеріорна ймовірність поєднує «схильність до дефекту» з «обсягом роботи».

Зауваження (базова частота). Формула Байєса пояснює, чому навіть дуже «точний» тест може давати оманливі висновки, коли шукана подія рідкісна. Якщо хвороба (чи дефект) трапляється рідко, то нечисленні справжні спрацювання тесту тонуть у морі хибних спрацювань від численної «здорової» більшості. Нехтування базовою частотою — поширена помилка міркування, яку виправляє саме байєсів підхід.


9.7 Теорема додавання ймовірностей

Теорема множення відповідає на питання «настануть обидві?». Теорема додавання (слайди 8–9) відповідає на питання «настане хоча б одна?».

Несумісні події

Теорема 9.6 (додавання для несумісних подій). Якщо події AA і BB несумісні (AB=A \cap B = \varnothing), то

P(AB)=P(A)+P(B).P(A \cup B) = P(A) + P(B).

Загальніше, для попарно несумісних A1,,AnA_1, \dots, A_n

P(A1An)=P(A1)++P(An).P(A_1 \cup \dots \cup A_n) = P(A_1) + \dots + P(A_n).

Доведення. У класичній моделі сприятливі для AA наслідки й сприятливі для BB наслідки не перетинаються (бо AB=A \cap B = \varnothing), тож кількість сприятливих для ABA \cup B дорівнює mA+mBm_A + m_B. Поділивши на nn, дістаємо P(AB)=mA+mBn=mAn+mBn=P(A)+P(B)P(A \cup B) = \tfrac{m_A + m_B}{n} = \tfrac{m_A}{n} + \tfrac{m_B}{n} = P(A) + P(B). Загальний випадок доводять індукцією за nn. \blacksquare

Приклад 9.9. Кинули кубик. Події «випало 11» і «випало 22» несумісні, тож P({1}{2})=16+16=13P(\{1\} \cup \{2\}) = \tfrac16 + \tfrac16 = \tfrac13. Зокрема, ймовірності всіх шести елементарних наслідків у сумі дають 66=1\tfrac66 = 1 — властивість повної групи.

Сумісні події

Для сумісних подій пряме додавання P(A)+P(B)P(A) + P(B) двічі враховує спільну частину ABA \cap B, тож її треба відняти.

Теорема 9.7 (додавання для сумісних подій). Для будь-яких подій AA і BB

P(AB)=P(A)+P(B)P(AB).P(A \cup B) = P(A) + P(B) - P(A \cap B).

Доведення. Запишемо ABA \cup B як об’єднання двох несумісних частин: AB=A(BA)A \cup B = A \cup (B \setminus A), де A(BA)=A \cap (B \setminus A) = \varnothing. За Теоремою 9.6

P(AB)=P(A)+P(BA).(1)P(A \cup B) = P(A) + P(B \setminus A). \tag{1}

Далі BB теж розпадається на дві несумісні частини: B=(AB)(BA)B = (A \cap B) \cup (B \setminus A), звідки

P(B)=P(AB)+P(BA)  P(BA)=P(B)P(AB).(2)P(B) = P(A \cap B) + P(B \setminus A) \ \Longrightarrow\ P(B \setminus A) = P(B) - P(A \cap B). \tag{2}

Підставивши (2) у (1), дістаємо P(AB)=P(A)+P(B)P(AB)P(A \cup B) = P(A) + P(B) - P(A \cap B). \blacksquare

Діаграма Венна двох сумісних подій A і B; заштриховано об'єднання A∪B; лінза перетину A∩B, яку віднімають, показана темнішим тоном

Зауваження. Це — ймовірнісний двійник формули включень-виключень для потужностей AB=A+BAB|A \cup B| = |A| + |B| - |A \cap B| з Лекції 1: поділивши її на Ω|\Omega| у класичній моделі, дістаємо саме Теорему 9.7. Для несумісних подій AB=A \cap B = \varnothing, доданок P(AB)=0P(A \cap B) = 0 зникає, і Теорема 9.7 переходить у Теорему 9.6.

Приклад 9.10 (карти). Витягли одну карту з 5252. Події «король» і «черва» сумісні (є король черв). Тому

P(король або черва)=P(король)+P(черва)P(король черв)=452+1352152=1652=413.P(\text{король або черва}) = P(\text{король}) + P(\text{черва}) - P(\text{король черв}) = \tfrac{4}{52} + \tfrac{13}{52} - \tfrac{1}{52} = \tfrac{16}{52} = \tfrac{4}{13}.


9.8 Розв’язані задачі (зі слайдів)

Ці чотири задачі зі слайдів охоплюють усі щойно доведені теореми. Кожен крозв’язок наведено повністю.

Задача 1 — дві білі кулі з двох ящиків (незалежні події)

Умова (слайд 10). У першому ящику 22 білі та 1010 червоних куль. У другому — 88 білих і 44 червоні. З кожного ящика витягли по одній кулі. Яка ймовірність, що обидві кулі білі?

Розв’язання. Позначимо AA — «біла куля з 11-го ящика», BB — «біла куля з 22-го ящика». Кожен ящик містить 1212 куль, тож за класичним означенням

P(A)=212=16,P(B)=812=23.P(A) = \frac{2}{12} = \frac16, \qquad P(B) = \frac{8}{12} = \frac23.

Добір із різних ящиків незалежний (вміст другого ящика не залежить від того, що взяли з першого), тому за теоремою множення для незалежних подій

P(AB)=P(A)P(B)=1623=218=190,111.P(A \cap B) = P(A)\,P(B) = \frac16 \cdot \frac23 = \frac{2}{18} = \frac{1}{9} \approx 0{,}111.

Відповідь: P(AB)=190,111P(A \cap B) = \tfrac19 \approx 0{,}111.

Задача 2 — дві білі кулі з одного ящика (залежні події)

Умова (слайд 11). У ящику 1010 білих і 1414 червоних куль (2424 разом). З ящика витягли 22 кулі (без повернення). Яка ймовірність, що обидві білі?

Розв’язання. Позначимо AA — «перша куля біла», BB — «друга куля біла». Тепер добори залежні: перша витягнута куля не повертається, тож умови для другого добору змінюються. Спершу

P(A)=1024.P(A) = \frac{10}{24}.

Якщо перша куля була білою, у ящику лишилося 99 білих із 2323 куль, тож

P(BA)=923.P(B \mid A) = \frac{9}{23}.

За теоремою множення для залежних подій

P(AB)=P(A)P(BA)=1024923=90552=15920,163.P(A \cap B) = P(A)\,P(B \mid A) = \frac{10}{24} \cdot \frac{9}{23} = \frac{90}{552} = \frac{15}{92} \approx 0{,}163.

Відповідь: P(AB)=15920,163P(A \cap B) = \tfrac{15}{92} \approx 0{,}163.

Зауваження (порівняйте Задачі 1 і 2). Формально обидві задачі про «дві білі кулі», але в Задачі 1 добори незалежні (P(AB)=P(A)P(B)P(A \cap B) = P(A)P(B)), а в Задачі 2 — залежні (P(AB)=P(A)P(BA)P(A \cap B) = P(A)P(B \mid A)). Різниця в одному множнику: P(B)P(B) проти P(BA)P(B \mid A). Саме тому, читаючи умову, найперше з’ясовуйте: чи повертають вибраний об’єкт назад?

Задача 3 — деталь першого сорту (умовна ймовірність)

Умова (слайд 12). Знайти ймовірність того, що навмання взята деталь є першого сорту, якщо відомо, що 4%4\% деталей браковані, а 75%75\% небракованих деталей — першого сорту.

Розв’язання. Позначимо AA — «деталь небракована», BB — «деталь першого сорту». Бути першого сорту може лише небракована деталь, тож подія «перший сорт» — це насправді ABA \cap B. З умови

P(A)=10,04=0,96,P(BA)=0,75.P(A) = 1 - 0{,}04 = 0{,}96, \qquad P(B \mid A) = 0{,}75.

За теоремою множення

P(AB)=P(A)P(BA)=0,960,75=0,72.P(A \cap B) = P(A)\,P(B \mid A) = 0{,}96 \cdot 0{,}75 = 0{,}72.

Відповідь: ймовірність узяти деталь першого сорту дорівнює 0,72=72%0{,}72 = 72\%.

Задача 4 — виграш хоча б по одному квитку (сумісні залежні події)

Умова (слайд 13). На 100100 лотерейних квитків є 55 виграшних. Куплено 22 квитки. Яка ймовірність виграшу хоча б по одному квитку?

Розв’язання. Позначимо AA — «виграв перший квиток», BB — «виграв другий квиток». Потрібно знайти P(AB)P(A \cup B). Події сумісні (можуть виграти обидва квитки) і залежні (квитки різні, добір без повернення).

Спосіб 1 — через протилежну подію (найшвидший). Протилежна подія «жодного виграшу» означає, що обидва квитки невиграшні. Невиграшних квитків 9595; за теоремою множення для залежних подій

P(AB)=951009499=893099000,902.P(\overline{A} \cap \overline{B}) = \frac{95}{100} \cdot \frac{94}{99} = \frac{8930}{9900} \approx 0{,}902.

Тоді за правилом протилежної події

P(AB)=1P(AB)10,902=0,098.P(A \cup B) = 1 - P(\overline{A} \cap \overline{B}) \approx 1 - 0{,}902 = 0{,}098.

Спосіб 2 — через теорему додавання для сумісних подій (перевірка).

P(A)=P(B)=5100=0,05,P(AB)=P(A)P(BA)=5100499=2099000,00202.P(A) = P(B) = \frac{5}{100} = 0{,}05, \qquad P(A \cap B) = P(A)\,P(B \mid A) = \frac{5}{100} \cdot \frac{4}{99} = \frac{20}{9900} \approx 0{,}00202.

P(AB)=P(A)+P(B)P(AB)=0,05+0,050,002020,098.P(A \cup B) = P(A) + P(B) - P(A \cap B) = 0{,}05 + 0{,}05 - 0{,}00202 \approx 0{,}098.

Обидва способи дають однакову відповідь. Відповідь: P(AB)0,098P(A \cup B) \approx 0{,}098.

Зауваження. Тут ми скористалися законом де Моргана: «жодного виграшу» =AB=AB= \overline{A \cup B} = \overline{A} \cap \overline{B} — обидва квитки невиграшні. Це знову підкреслює: перекласти умову мовою множин — половина розв’язання.


9.9 Основи комбінаторики

Щоб застосувати класичне означення P(A)=m/nP(A) = m/n, треба полічити наслідки. Комбінаторика дає для цього кілька засобів, що спираються на два прості правила.

Правило суми й правило добутку

Правило суми. Якщо об’єкт можна вибрати одним із kk взаємно виключних способів, причому перший спосіб дає m1m_1 варіантів, другий — m2m_2, …, kk-й — mkm_k, то загальна кількість варіантів дорівнює m1+m2++mkm_1 + m_2 + \dots + m_k. («Або / окремі випадки» — додаємо.)

Правило добутку. Якщо процедура складається з kk послідовних кроків, де крок 11 можна виконати n1n_1 способами, крок 22n2n_2 способами незалежно від першого вибору, …, крок kknkn_k способами, то всю процедуру можна виконати

n1n2nk способами.n_1 \cdot n_2 \cdots n_k \ \text{способами.}

І / послідовні кроки» — множимо.)

Обґрунтування правила добутку. Результат процедури — це впорядкований набір (x1,,xk)(x_1, \dots, x_k), де xix_i — один із nin_i варіантів кроку ii. Множина всіх таких наборів — це декартів добуток X1××XkX_1 \times \dots \times X_k з Лекції 1, а потужність скінченного декартового добутку дорівнює добутку потужностей. \blacksquare

Приклад 9.11. Артикул складається з 33 великих латинських літер, за якими йдуть 22 цифри. За правилом добутку різних артикулів 2626261010=263102=175760026 \cdot 26 \cdot 26 \cdot 10 \cdot 10 = 26^3 \cdot 10^2 = 1\,757\,600.

Перестановки

Означення (перестановки). Перестановка nn різних об’єктів — це впорядкований набір усіх nn об’єктів без повторів (слайд 14). Кількість перестановок дорівнює

Pn=n!=123n,за домовленістю 0!=1.P_n = n! = 1 \cdot 2 \cdot 3 \cdots n, \qquad \text{за домовленістю } 0! = 1.

Доведення. Заповнюємо nn позицій зліва направо. Першу позицію можна зайняти nn способами, другу — n1n - 1 (один об’єкт уже використано), …, останню — 11 способом. За правилом добутку кількість способів дорівнює n(n1)1=n!n(n-1) \cdots 1 = n!. \blacksquare

Приклад 9.12. Кількість способів розставити 55 книжок на полиці: P5=5!=120P_5 = 5! = 120.

Розміщення

Означення (розміщення). Розміщення з nn по kk — це впорядкована вибірка kk різних об’єктів з nn (порядок важливий, без повторів; слайд 15). Кількість розміщень

Ank=n!(nk)!=n(n1)(nk+1)k множників.A_n^{k} = \frac{n!}{(n-k)!} = \underbrace{n(n-1) \cdots (n-k+1)}_{k \text{ множників}}.

Доведення. Заповнюємо kk упорядкованих позицій із nn об’єктів: nn варіантів для першої позиції, n1n - 1 для другої, …, nk+1n - k + 1 для kk-ї. Добуток цих kk множників дорівнює n!(nk)!\dfrac{n!}{(n-k)!} після скорочення «хвоста» (nk)(nk1)1(n-k)(n-k-1) \cdots 1. \blacksquare

Приклад 9.13. Скількома способами з 1010 працівників можна призначити керівника проєкту й окремо його заступника (дві різні посади)? Порядок важливий, тож це розміщення: A102=109=90A_{10}^{2} = 10 \cdot 9 = 90.

Сполучення

Означення (сполучення). Сполучення з nn по kk — це вибірка kk об’єктів з nn, у якій порядок не має значення (без повторів; слайд 16), тобто kk-елементна підмножина. Кількість сполучень

Cnk=n!k!(nk)!,0kn.C_n^{k} = \frac{n!}{k!\,(n-k)!}, \qquad 0 \le k \le n.

Доведення. Полічимо впорядковані kk-вибірки двома способами. З одного боку, їх Ank=n!(nk)!A_n^{k} = \dfrac{n!}{(n-k)!}. З іншого боку, щоб побудувати впорядковану вибірку, можна спершу обрати підмножину з kk об’єктів (нехай таких підмножин CnkC_n^{k}), а потім упорядкувати її k!k! способами (перестановки). За правилом добутку Ank=Cnkk!A_n^{k} = C_n^{k} \cdot k!, звідки

Cnk=Ankk!=n!k!(nk)!.C_n^{k} = \frac{A_n^{k}}{k!} = \frac{n!}{k!\,(n-k)!}. \qquad \blacksquare

Схема трьох комбінаторних вибірок для n=4: перестановки (усі 4 позиції впорядковані, P4=4!), розміщення (упорядкована вибірка k=2 позицій, A) та сполучення (невпорядкована підмножина з 2 елементів, C)

Приклад 9.14. Комісію з 33 осіб обирають із 1010. Порядок неважливий, тож це сполучення:

C103=10!3!7!=1098321=120.C_{10}^{3} = \frac{10!}{3!\,7!} = \frac{10 \cdot 9 \cdot 8}{3 \cdot 2 \cdot 1} = 120.

Порівняйте з упорядкованим добором A103=720A_{10}^{3} = 720: поділ на 3!=63! = 6 прибирає порядок.

Типова помилка (розміщення чи сполучення?). Єдине питання, що обирає засіб: чи дає перестановка вибраних об’єктів інший результат? Так \Rightarrow розміщення AnkA_n^{k}; ні \Rightarrow сполучення CnkC_n^{k}. «Керівник і заступник» — це впорядкована пара (розміщення); «комісія з трьох» — неупорядкована (сполучення).

Трикутник Паскаля та біном Ньютона

Числа CnkC_n^{k} мають дві властивості, якими постійно користуються.

  • Симетрія: Cnk=CnnkC_n^{k} = C_n^{\,n-k} (обрати kk об’єктів «усередину» — те саме, що обрати nkn - k об’єктів «назовні»).
  • Тотожність Паскаля: Cnk=Cn1k1+Cn1kC_n^{k} = C_{n-1}^{\,k-1} + C_{n-1}^{k}.

Твердження 9.8 (тотожність Паскаля). Cnk=Cn1k1+Cn1kC_n^{k} = C_{n-1}^{\,k-1} + C_{n-1}^{k}.

Доведення (комбінаторне). Зафіксуємо один конкретний об’єкт xx серед nn. Кожна kk-підмножина або містить xx, або ні. Ті, що містять xx, доукомплектовують решту k1k - 1 членів з інших n1n - 1 об’єктів — таких Cn1k1C_{n-1}^{\,k-1}. Ті, що не містять xx, обирають усі kk членів з інших n1n - 1 об’єктів — таких Cn1kC_{n-1}^{k}. Ці випадки не перетинаються й вичерпують усі підмножини, тож їхні кількості додаються. \blacksquare

Ця тотожність породжує трикутник Паскаля: кожне число дорівнює сумі двох чисел над ним. У рядку з номером nn стоять коефіцієнти Cn0,Cn1,,CnnC_n^{0}, C_n^{1}, \dots, C_n^{n}.

Трикутник Паскаля до п'ятого рядка; кожне число дорівнює сумі двох чисел над ним, що ілюструє тотожність Паскаля; підсвічено одне таке додавання

Ті самі числа CnkC_n^{k} є коефіцієнтами розкладу степеня двочлена — бінома Ньютона:

(a+b)n=k=0nCnkakbnk=Cn0bn+Cn1abn1++Cnnan.(a + b)^{n} = \sum_{k=0}^{n} C_n^{k}\, a^{k} b^{\,n-k} = C_n^{0} b^{n} + C_n^{1} a b^{n-1} + \dots + C_n^{n} a^{n}.

Зокрема, при a=b=1a = b = 1 маємо k=0nCnk=2n\sum_{k=0}^{n} C_n^{k} = 2^{n} — це знову кількість усіх підмножин nn-елементної множини (Теорема 1.17 з Лекції 1). А підставивши a=pa = p, b=1pb = 1 - p, дістаємо, що ваги біноміального розподілу Cnkpk(1p)nkC_n^{k} p^{k}(1-p)^{n-k} у сумі дають 11 (це знадобиться в §9.11).


9.10 Статистична (частотна) ймовірність і закон великих чисел

Класичне означення P(A)=m/nP(A) = m/n працює лише за рівноможливих наслідків. Коли симетрії немає (погнута монета, невідома частка браку в новому технологічному процесі), ймовірність оцінюють емпірично (слайд 17).

Означення (відносна частота). Якщо в NN незалежних випробуваннях подія AA настала NAN_A разів, то її відносною частотою називають

W(A)=NAN.W(A) = \frac{N_A}{N}.

Відносна частота — величина дослідна: її обчислюють після серії випробувань, і за різних серій вона трохи різна. Проте спостереження показує, що зі зростанням NN частота стабілізується навколо певного числа.

Означення (статистична ймовірність). Статистичною ймовірністю події AA називають число, навколо якого стабілізується відносна частота W(A)W(A) при необмеженому зростанні кількості випробувань NN.

Закон великих чисел (Я. Бернуллі, спрощене формулювання). При збільшенні числа незалежних випробувань NN відносна частота W(A)=NA/NW(A) = N_A / N як завгодно мало відрізняється (зі скільки завгодно високою певністю) від імовірності P(A)P(A):

W(A)=NAN  N  P(A).W(A) = \frac{N_A}{N} \xrightarrow[\;N \to \infty\;]{} P(A).

Саме закон великих чисел робить осмисленими вислови «ймовірність дощу», «ймовірність відмови жорсткого диска протягом року» чи «частка браку на друкарській лінії» — там, де жодної симетрії немає. Для правильних (симетричних) механізмів статистична й класична ймовірності збігаються: полічена частка парних граней кубика в довгій серії кидків наближається до класичного 12\tfrac12.

Приклад 9.15 (контроль тиражу). Друкарня перевірила N=5000N = 5000 відбитків нового тиражу й виявила NA=65N_A = 65 дефектних. Статистична оцінка ймовірності дефекту —

W(A)=655000=0,013=1,3%.W(A) = \frac{65}{5000} = 0{,}013 = 1{,}3\%.

Оскільки процес новий і симетрії тут немає, класичне означення незастосовне, а частота 0,0130{,}013 — найкраща доступна оцінка ймовірності дефекту, яку далі використовують для планування контролю.


9.11 Випадкові величини: сподівання, дисперсія, розподіли

Часто нас цікавить не який наслідок настав, а число, з ним пов’язане, — кількість дефектних відбитків у вибірці, сума очок на двох кубиках, час роботи механізму до збою.

Означення (випадкова величина). Випадкова величина — це правило XX, що кожному елементарному наслідку ωΩ\omega \in \Omega зіставляє число X(ω)X(\omega). Якщо множина її значень скінченна (чи зліченна), величину називають дискретною; якщо величина набуває всіх значень із проміжку — неперервною.

Дискретну величину задають законом розподілу — таблицею її значень xix_i та відповідних імовірностей pi=P(X=xi)p_i = P(X = x_i), причому pi0p_i \ge 0 і ipi=1\sum_i p_i = 1.

xix_i x1x_1 x2x_2 \dots xkx_k
pip_i p1p_1 p2p_2 \dots pkp_k

Математичне сподівання

Означення (математичне сподівання). Математичним сподіванням дискретної випадкової величини XX називають зважене за ймовірностями середнє її значень (слайд 18):

M(X)=ixipi=x1p1+x2p2++xkpk.M(X) = \sum_{i} x_i\, p_i = x_1 p_1 + x_2 p_2 + \dots + x_k p_k.

(Для неперервної величини суму заступає інтеграл M(X)=xf(x)dxM(X) = \int x\, f(x)\, dx, де ff — щільність; це поза межами дискретного курсу.)

Сподівання — це «довгострокове середнє» значення XX, до якого прямує середнє арифметичне результатів при багаторазовому повторенні досліду (це гарантує закон великих чисел). Основні властивості (наводимо без доведення):

M(C)=C,M(CX)=CM(X),M(X+Y)=M(X)+M(Y),M(C) = C, \qquad M(CX) = C\,M(X), \qquad M(X + Y) = M(X) + M(Y),

причому остання (адитивність) правдива незалежно від того, залежні XX і YY чи ні.

Приклад 9.16 (кубик). Для числа очок на одному кубику M(X)=k=16k16=1+2+3+4+5+66=216=3,5M(X) = \sum_{k=1}^{6} k \cdot \tfrac16 = \tfrac{1+2+3+4+5+6}{6} = \tfrac{21}{6} = 3{,}5. Саме число 3,53{,}5 жодного разу не випадає — сподівання є середнім, а не прогнозом окремого кидка.

Дисперсія та середньоквадратичне відхилення

Саме лише сподівання не описує розкид значень. Його вимірюють дисперсією (слайд 19).

Означення (дисперсія, середньоквадратичне відхилення). Дисперсія випадкової величини XX — це сподіване значення квадрата відхилення від середнього:

D(X)=M((XM(X))2)=i(xiM(X))2pi.D(X) = M\big( (X - M(X))^2 \big) = \sum_i \big( x_i - M(X) \big)^2 p_i.

На практиці зручніша рівносильна формула

D(X)=M(X2)(M(X))2.D(X) = M(X^2) - \big( M(X) \big)^2.

Середньоквадратичним відхиленням називають σ(X)=D(X)\sigma(X) = \sqrt{D(X)}; воно має ту саму розмірність, що й сама величина XX.

Корисні властивості: D(C)=0D(C) = 0, D(CX)=C2D(X)D(CX) = C^2 D(X), а для незалежних XX і YYD(X+Y)=D(X)+D(Y)D(X + Y) = D(X) + D(Y).

Приклад 9.17 (число дефектів — повний розрахунок MM, DD, σ\sigma). Нехай XX — число дефектних відбитків у невеликій пачці, із законом розподілу

xix_i 00 11 22 33
pip_i 0,500{,}50 0,300{,}30 0,150{,}15 0,050{,}05

Спершу сподівання:

M(X)=00,5+10,3+20,15+30,05=0,3+0,3+0,15=0,75.M(X) = 0 \cdot 0{,}5 + 1 \cdot 0{,}3 + 2 \cdot 0{,}15 + 3 \cdot 0{,}05 = 0{,}3 + 0{,}3 + 0{,}15 = 0{,}75.

Далі M(X2)=020,5+120,3+220,15+320,05=0,3+0,6+0,45=1,35M(X^2) = 0^2 \cdot 0{,}5 + 1^2 \cdot 0{,}3 + 2^2 \cdot 0{,}15 + 3^2 \cdot 0{,}05 = 0{,}3 + 0{,}6 + 0{,}45 = 1{,}35, звідки

D(X)=M(X2)(M(X))2=1,350,752=1,350,5625=0,7875,σ=0,78750,887.D(X) = M(X^2) - (M(X))^2 = 1{,}35 - 0{,}75^2 = 1{,}35 - 0{,}5625 = 0{,}7875, \qquad \sigma = \sqrt{0{,}7875} \approx 0{,}887.

Отже, у пачці в середньому 0,750{,}75 дефектного відбитка з розкидом близько ±0,89\pm 0{,}89.

Види розподілів (оглядово)

Кілька законів розподілу трапляються так часто, що мають власні назви (слайди 20–21).

  • Рівномірний дискретний. Усі nn значень рівноможливі, pi=1np_i = \tfrac1n; це модель правильного кубика чи монети. Тут M(X)=1nixiM(X) = \tfrac1n \sum_i x_i.
  • Біноміальний. Кількість «успіхів» XX у nn незалежних однакових випробуваннях, де успіх у кожному має ймовірність pp:

    P(X=k)=Cnkpk(1p)nk,M(X)=np,D(X)=np(1p).P(X = k) = C_n^{k}\, p^{k} (1-p)^{\,n-k}, \qquad M(X) = np, \qquad D(X) = np(1-p).

    Це основна модель вибіркового контролю: XX — число дефектних у вибірці обсягу nn за частки браку pp.
  • Розподіл Пуассона. Наближення біноміального при великому nn і малому pp (рідкісні події): P(X=k)=λkeλk!P(X = k) = \dfrac{\lambda^{k} e^{-\lambda}}{k!}, де λ=np\lambda = np; тут M(X)=D(X)=λM(X) = D(X) = \lambda. Моделює, наприклад, число рідкісних дефектів на великому аркуші.
  • Нормальний (Гаусса). Неперервний розподіл із симетричною «дзвоноподібною» кривою; до нього прямує сума багатьох дрібних незалежних впливів (центральна гранична теорема). Ним описують сумарні колірні відхилення, коливання щільності фарби тощо. Детально — у подальших курсах.

Стовпчикова діаграма біноміального розподілу числа дефектних копій у вибірці з 8 за частки браку 0,25; стовпці бірюзові, осі та підписи темно-сині; позначено математичне сподівання M(X)=2

Приклад 9.18 (біноміальна модель вибірки). Частка браку на лінії p=0,25p = 0{,}25; беруть вибірку з n=8n = 8 копій. Число дефектних XX має біноміальний розподіл (діаграма вище). Сподіване число дефектних M(X)=np=80,25=2M(X) = np = 8 \cdot 0{,}25 = 2, дисперсія D(X)=np(1p)=80,250,75=1,5D(X) = np(1-p) = 8 \cdot 0{,}25 \cdot 0{,}75 = 1{,}5, а σ=1,51,22\sigma = \sqrt{1{,}5} \approx 1{,}22. Наприклад, ймовірність, що у вибірці не буде дефектних, P(X=0)=C80(0,25)0(0,75)8=0,7580,100P(X = 0) = C_8^0 (0{,}25)^0 (0{,}75)^8 = 0{,}75^8 \approx 0{,}100.


Застосування у видавництві та поліграфії

Ймовірнісні моделі — це щоденний інструмент контролю якості в поліграфії.

  • Вибірковий приймальний контроль тиражу. Перевірити весь тираж дорого, тому з нього беруть вибірку обсягу nn і рахують дефектні відбитки. Число дефектних у вибірці — біноміальна (при доборі без повернення — гіпергеометрична) випадкова величина; за нею ухвалюють рішення «прийняти / відхилити партію». Саме на цьому ґрунтуються стандарти приймального контролю за рівнем якості AQL (напр., ISO 2859-1).
  • Ймовірність бездефектного відбитка та надійність обладнання. Якщо ймовірність дефекту одного відбитка дорівнює qq, то ймовірність, що вся сторінка з rr незалежних елементів вийде бездефектною, — (1q)r (1-q)^{r} ; ймовірність «хоча б одного дефекту» зручно рахувати через протилежну подію (§9.3). Так само оцінюють імовірність збою друкарської машини за зміну.
  • Діагностика джерела дефекту (Байєс). Виявивши брак, за формулою Байєса визначають, яка машина, секція чи фарбова група найімовірніше його спричинила (Приклад 9.8), — ймовірнісний аналог «пошуку винуватця».
  • Колірна відповідність. Відхилення кольору ΔE\Delta E від еталона — випадкова величина; контроль полягає в оцінці ймовірності P(ΔE>допуск)P(\Delta E > \text{допуск}) та сподіваного відхилення M(ΔE)M(\Delta E) і його розкиду σ\sigma.
  • Комбінаторика макета. Кількість варіантів спуску полос, розкладки кольорів чи компонування обкладинки — це задачі на перестановки, розміщення й сполучення (§9.9); вони визначають, скільки макетів доведеться перебрати чи згенерувати.
  • Розпізнавання й фільтрування тексту. Байєсів класифікатор (наївний Байєс) — пряме застосування формули Байєса — лежить в основі OCR і фільтрів спаму, потрібних у редакційно-видавничих системах.

Підсумок

  • Випадковий дослід має простір елементарних наслідків Ω\Omega; подія — це підмножина AΩA \subseteq \Omega. Події є множинами, тож алгебра об’єднання / перетину / доповнення й закони де Моргана з Лекції 1 переносяться дослівно.
  • Типи подій: достовірна (Ω\Omega), неможлива (\varnothing), сумісні проти несумісних (AB=A \cap B = \varnothing), протилежні A,AA, \overline{A}, повна група (розбиття), незалежні проти залежних.
  • Класична ймовірність P(A)=mn=AΩP(A) = \tfrac{m}{n} = \tfrac{|A|}{|\Omega|} для рівноможливих наслідків зводить задачу до підрахунку; звідси 0P(A)10 \le P(A) \le 1, P(Ω)=1P(\Omega) = 1, P()=0P(\varnothing) = 0 та правило протилежної події P(A)=1P(A)P(\overline{A}) = 1 - P(A).
  • Умовна ймовірність P(AB)=P(AB)P(B)P(A \mid B) = \dfrac{P(A \cap B)}{P(B)} звужує простір до BB; для незалежних подій P(AB)=P(A)P(A \mid B) = P(A).
  • Теорема множення: P(AB)=P(A)P(BA)P(A \cap B) = P(A)\,P(B \mid A) (залежні) =P(A)P(B)= P(A)\,P(B) (незалежні); узагальнення — ланцюгове правило.
  • Формула повної ймовірності P(A)=iP(Hi)P(AHi)P(A) = \sum_i P(H_i)\,P(A \mid H_i) (доведена через розбиття) та формула Байєса P(HiA)=P(Hi)P(AHi)jP(Hj)P(AHj)P(H_i \mid A) = \dfrac{P(H_i)\,P(A \mid H_i)}{\sum_j P(H_j)\,P(A \mid H_j)} перетворюють апріорі й правдоподібності на апостеріорі.
  • Теорема додавання: P(AB)=P(A)+P(B)P(A \cup B) = P(A) + P(B) (несумісні) і P(AB)=P(A)+P(B)P(AB)P(A \cup B) = P(A) + P(B) - P(A \cap B) (сумісні) — ймовірнісні включення-виключення.
  • Комбінаторика: правила суми й добутку \Rightarrow перестановки Pn=n!P_n = n!, розміщення Ank=n!(nk)!A_n^{k} = \dfrac{n!}{(n-k)!}, сполучення Cnk=n!k!(nk)!C_n^{k} = \dfrac{n!}{k!(n-k)!} (звідки Ank=Cnkk!A_n^{k} = C_n^{k}\, k!); тотожність Паскаля Cnk=Cn1k1+Cn1kC_n^{k} = C_{n-1}^{k-1} + C_{n-1}^{k} і біном Ньютона (a+b)n=kCnkakbnk(a+b)^n = \sum_k C_n^{k} a^k b^{n-k}.
  • Статистична ймовірність W(A)=NANW(A) = \tfrac{N_A}{N} оцінює P(A)P(A) там, де симетрії немає; закон великих чисел гарантує W(A)P(A)W(A) \to P(A).
  • Випадкові величини: закон розподілу; математичне сподівання M(X)=ixipiM(X) = \sum_i x_i p_i, дисперсія D(X)=M(X2)(M(X))2D(X) = M(X^2) - (M(X))^2 і середньоквадратичне відхилення σ=D(X)\sigma = \sqrt{D(X)}; типові розподіли — рівномірний, біноміальний, Пуассона, нормальний.

Історична довідка. Азартні ігри давні, але їхня математика — ні. Зачинателями теорії ймовірностей вважають Блеза Паскаля і П’єра Ферма, які в листуванні 16541654 року розв’язали кілька гральних задач. Якоб Бернуллі довів перший закон великих чисел (опубл. 17131713). Есе Томаса Баєса про обернену ймовірність вийшло посмертно (17631763), а П’єр-Симон Лаплас систематизував класичне означення (близько 18121812). Сучасну аксіоматичну основу теорії заклав А. М. Колмогоров (19331933).

Вправи

Для розігріву

  1. Кинули один кубик. Випишіть Ω\Omega та як підмножини — події A=A = «непарне» і B=B = «не менше за 44». Обчисліть P(A)P(A), P(B)P(B), P(AB)P(A \cap B), P(AB)P(A \cup B) і перевірте теорему додавання.
  2. Поясніть одним реченням різницю між: а) несумісними і незалежними подіями; б) P(AB)P(A \mid B) і P(AB)P(A \cap B).
  3. Зі стандартної колоди на 5252 карти знайдіть P(червона)P(\text{червона}), P(дама)P(\text{дама}) та P(червона або дама)P(\text{червона або дама}).
  4. Комісію з 44 осіб обирають із 1212. Скільки можливих комісій? А якщо ці четверо обіймають різні посади (голова, секретар, скарбник, вебмайстер) — скільки призначень?
  5. Обчисліть P5P_5, A62A_6^{2} та C62C_6^{2}. Перевірте рівність A62=C622!A_6^{2} = C_6^{2} \cdot 2!.

Стандартні

  1. (Виведення правила.) Спираючись лише на класичне означення й теорему додавання для несумісних подій, доведіть, що для сумісних подій P(AB)=P(A)+P(B)P(AB)P(A \cup B) = P(A) + P(B) - P(A \cap B).
  2. У ящику 55 червоних, 33 зелені та 22 сині кулі; витягли 33 без повернення. Знайдіть ймовірність того, що: а) усі три червоні; б) по одній кулі кожного кольору; в) хоча б одна синя.
  3. У коробці 88 придатних і 22 браковані мікросхеми; витягли 22 без повернення. Знайдіть ймовірність, що обидві придатні, — спершу теоремою множення, а потім через сполучення; переконайтеся, що відповіді збігаються.
  4. (Повна ймовірність.) Дільниці I, II, III дають 50%50\%, 30%30\%, 20%20\% продукції з частками браку 2%2\%, 3%3\%, 5%5\%. Знайдіть ймовірність того, що навмання взятий виріб бракований.
  5. (Байєс.) Для умов вправи 9 навмання взятий виріб виявився бракованим. Яка дільниця найімовірніше його виготовила? Наведіть усі три апостеріорні ймовірності.
  6. Випадкова величина XX задана законом розподілу: значення 1,0,2-1, 0, 2 з імовірностями 0,20{,}2, 0,50{,}5, 0,30{,}3. Обчисліть M(X)M(X), D(X)D(X) і σ(X)\sigma(X).

Підвищеної складності

  1. Кинули два кубики. Знайдіть P(сума=9)P(\text{сума} = 9), P(сума парна)P(\text{сума парна}) та P(сума=9сума непарна)P(\text{сума} = 9 \mid \text{сума непарна}).
  2. (Схема Бернуллі.) Частка браку на лінії p=0,1p = 0{,}1; узяли вибірку з n=6n = 6 відбитків. Використовуючи біноміальний розподіл, знайдіть ймовірність того, що бракованих буде рівно 11, та ймовірність того, що бракованих не буде жодного. Обчисліть M(X)M(X) і D(X)D(X).
  3. (Протилежна подія.) За скількох незалежних кидків одного кубика ймовірність дістати хоча б одну шістку вперше перевищить 12\tfrac12? Обґрунтуйте через протилежну подію.
  4. Доведіть тотожність Cn0+Cn1++Cnn=2nC_n^{0} + C_n^{1} + \dots + C_n^{n} = 2^{n} двома способами: а) як біном Ньютона при a=b=1a = b = 1; б) полічивши всі підмножини nn-елементної множини двома способами (зв’яжіть із Теоремою 1.17 Лекції 1).
  5. (Помилка прокурора.) Тест на певну ознаку, наявну в 2%2\% виробів, має чутливість 95%95\% (виявляє 95%95\% справжніх випадків) і частку хибних спрацювань 10%10\%. Виріб дав позитивний результат. Знайдіть за формулою Байєса ймовірність того, що ознака справді наявна, і поясніть, чому вона значно менша за 95%95\%.

Lectures/ODM-L09.md · 73.3 KB · updated 2026-08-04 14:34