Radiologická fyzika pravděpodobnost měření a zpracování dat Měření tlaku – nejjednodušší úkon u lékaře aneb jak zacházet s měřenými hodnotami? Pokus s měřením tlaku Pravděpodobnost aneb matematika náhody Co je to pravděpodobnost aneb proč nesázet Sportku? Víte například, že - pravděpodobnost hlavní výhry ve Sportce je sedm milióntin procenta, tj. ani ne jedna ku deseti miliónům? - pravděpodobnost páté ceny je už skoro dvouprocentní? - pravděpodobnost, že mezi 40 lidmi jsou alespoň dva, kteří mají narozeniny ve stejný den, je skoro devadesátiprocentní? A víte, - co to vůbec je pravděpodobnost? - co je to medián, co průměr, a kolik typů průměrů existuje? - co je správnost a co přesnost měření? - co znamená zápis „hodnota veličiny X je (2,51 ± 0,07) m“ ? Náhodné jevy a jejich pravděpodobnost Definice pravděpodobnosti Pravděpodobnost nastoupení jevu A je podílem počtu případů M, v nichž jev A nastal (čitatel), a počtu N všech možných případů (jmenovatel). Úkol 1: a) Jaká je pravděpodobnost, že při hodu kostkou padne sudé číslo? b) Jaká je pravděpodobnost, že padne číslo dělitelné třemi? c) Jaká je pravděpodobnost, že padne číslo menší než 3? Úkol 2: a) Jaká je pravděpodobnost, že při hodu dvěma kostkami (současně) padne součet sedm? b) Jaká je pravděpodobnost, že součin padnuvších čísel bude lichý? Jakých hodnot může pravděpodobnost v principu nabývat? Řešení předchozích úkolů Řešení úkolu 1b) N = 6, M = 2 (příznivé jsou případy, kdy padne trojka nebo šestka) p = 2/6 = 1/3 Řešení úkolu 2a) N = 36 (každá ze 6 možností, které mohou padnout na první kostce, se nezávisle kombinuje s každou ze 6 možností na druhé kostce), M = 6 (jednička na první a šestka na druhé kostce, nebo naopak, dvojka na první a pětka na druhé kostce, nebo naopak, trojka na první a čtyřka na druhé kostce, nebo naopak), p = 6/36 = 1/6 Řešení úkolu 2b) N = 36, M = 9, p = 9/36 = 1/4 Řešení úkolu 2b sami zdůvodněte. Ještě jeden úkol ryze praktický V zásuvce jsou ponožky tří barev. Červené (Č), zelené (Z) a modré (M). Je jich tam od každé barvy hodně. Student jde na schůzku a chce si vzít čisté ponožky. Náhle zhasne světlo. Student vytáhne potmě dvě ponožky. • Jaká je pravděpodobnost, že budou mít stejnou barvu? • Jaká je pravděpodobnost, že nebudou mít stejnou barvu? • Jaká je pravděpodobnost, že ponožky budou mít stejnou barvu, ale ne červenou? Všimněte si výsledků a) a b). Co pro ně platí? Co myslíte, je to náhoda, nebo to tak být musí? Jev jistý – jev s jednotkovou pravděpodobností . Jev nemožný – jev s nulovou pravděpodobností. Cifry, kostky, karty aneb kombinatorická průprava výběry k-té třídy prvků z n prvků výběry k prvků z množiny obsahující n prvků podle jistých pravidel Jak určit počet možných výběrů Vzorce pro variace Určení vzorce pro variace opakováním - Kolika způsoby lze z n prvků vybrat první? … n způsobů - Kolika způsoby lze z n prvků vybrat druhý? … n způsobů - Kolika způsoby lze vybrat první dva prvky? … n n= n^2 způsobů - Kolika způsoby lze tedy vybrat k prvků? Doplňte sami. Určení vzorce pro variace bez opakování - Kolika způsoby lze z n prvků vybrat první? … n způsobů - Kolika způsoby lze ze zbývajících n – 1 prvků vybrat druhý? … n – 1 způsobů - Kolika způsoby lze tedy vybrat první dva prvky? … n(n – 1) - Kolika způsoby lze nakonec ze zbývajících n – (k – 1) prvků vybrat k-tý? … (n – k + 1) - Kolika způsoby lze tedy vybrat k prvků? … n(n – 1)(n – 2)…(n – k + 1) Vzorce pro kombinace Určení vzorce pro kombinace bez opakování - Kolika způsoby lze uspořádat k neopakujících se prvků? … k! Tyto možnosti, lišící se pouze pořadím, jsou ekvivalentní. - Proto je C[k](n) = V[k](n) / k! = … Výsledný vzorec zapište sami. Určení vzorce pro kombinace s opakováním - Jak vypadají kombinace s opakováním k barevných kuliček při výběru z n možných barev? Zvolme např. n = 7, k = 14 Příklady na pravděpodobnost - I Příklad 1. Jaká je pravděpodobnost hlavní výhry ve Sportce? Tah sportky představuje výběr šesti ze čtyřiceti devíti čísel. Příklad 2. Jaká je pravděpodobnost páté ceny ve Sportce? Ze sedmi tažených je třeba uhodnout tři čísla. Příklady na pravděpodobnost - II Příklad 3. Jaká je pravděpodobnost, že ve hře typu Šance milion uhodnete správně taženou skupinu cifer? (Z každého ze šesti bubnů obsahujících cifry 0, 1, 2, … , 9 se náhodně vybere jedna.) Příklad 4. Jaká je pravděpodobnost v předchozí hře, bude-li k dispozici pouze jeden buben, který obsahuje každou z cifer právě jednou? Jak s pravděpodobnostmi počítat aneb pravděpodobnosti „složených“ jevů Někdy je třeba určit pravděpodobnosti jevů, které jsou nějakým způsobem „složeny“ z jevů jednodušších. Uvažujme o dvou jevech A a B, jejichž pravděpodobnosti jsou známy, p(A), p(B). Definujme nové jevy C a D jako C … jevy A a B nastanou současně D … nastane jev A nebo B (v principu zahrnuje i možnost, že nastanou oba) Za určitých podmínek lze pravděpodobnosti jevů C a D určit Pomocí pravděpodobností p(A) a p(B). Nezávislé jevy Jaká je pravděpodobnost, že při současném hodu dvěma kostkami padne na obou šestka? Uvědomte si: To, co padne na jedné kostce, je nezávislé na výsledku druhé kostky. Jev A: Na první kostce padne šestka … p(A) = 1/6 Jev B: Na druhé kostce padne šestka … p(B) = 1/6 Jev C: Jevy A a B nastanou současně … N = 6 6 = 36 možností, M = 1, p(C) = M / N = 1/36 = p(A) p(B) Pravděpodobnost současného nástupu nezávislých jevů je rovna součinu pravděpodobností těchto jevů. Neslučitelné jevy Jaká je pravděpodobnost, že při hodu kostkou padne některé ze dvou nejvyšších čísel, tj. padne šestka nebo pětka? Uvědomte si: Skutečnosti, že šestka i pětka padnou při stejném hodu, jsou neslučitelné. Jev A: Na kostce padne šestka … p(A) = 1/6 Jev B: Na kostce padne pětka … p(B) = 1/6 Jev C: Nastane buď jev A nebo jev B … N = 6 možností, M = 2, p(C) = M / N = 2/6 = p(A) + p(B) Pravděpodobnost nástupu některého z jevů, z nichž každé dva jsou neslučitelné, je rovna součtu pravděpodobností těchto jevů. Dokážete vysvětlit, proč je součet pravděpodobností jevu A a jevu opačného, tj. že jev A nenastane, je rovna 1 ? Příklady na pravděpodobnost - III Příklad 5. Jaká je pravděpodobnost, že ve skupině k osob mají alespoň dvě narozeniny ve stejný den? Rok má n = 365 dní. Určíme nejprve pravděpodobnost jevu A, že každá z osob má narozeniny v jiný den. Počet případů možných pro tento jev je N = V^/[k](n) (variace s opakováním – v principu může mít každá z osob narozeniny v kterýkoli den). Počet případů příznivých je M = V[k](n) (variace bez opakování – nechceme, aby se narozeninový den zopakoval u více osob). Jev, který nás zajímá, je opačným jevem k jevu A, jeho pravděpodobnost je tedy p = 1 – V[k](n) / V^/[k](n). Vypočtěte si tuto pravděpodobnost pro 40 osob. Důležitý příklad – Bernoulliův pokus - I Provedení pokusu a označení - Nastane-li předem definovaný jev A (například „padne šestka“), nazveme to zdarem, v opačném případě nezdarem. - Pravděpodobnost zdaru označíme p (p = 1/6), pravděpodobnost nezdaru je 1 – p (tedy 5/6). - n-krát nezávisle provedeme pokus (například hod kostkou). Jaký jev nás zajímá - Jev B … Právě při x provedeních z celkového počtu n provedení pokusu nastane zdar. Důležitý příklad – Bernoulliův pokus - II Které další jevy s tím souvisí - Jevy A[j] pro j = 1, …, x ... při j-tém provedení pokusu nastane zdar. - Jevy B[k] pro k = x+1, …, n ... při k-tém provedení pokusu nezdar. - Jev C ... Jevy A[1][ ]až A[x][ ]a B[x+][1] až B[n] [ ]nastoupí současně, tj. právě při prvních x opakováních pokusu nastane zdar, při zbývajících nezdar. - Nezáleží nám ale na tom, při kterých x ze všech n opakování pokusu nastal zdar. Možností, kdy zdar nastal právě při x ze všech opakování pokusu, je C[x](n) = n!/[x!(n – x)!]. - Výsledná pravděpodobnost jevu B - Úkol: Vypočtěte pravděpodobnost Bernoulliova pokusu pro 5 opakování hodu kostkou a dva zdary, a pro 5 opakování hodu mincí a žádný zdar. Pro jaké x je při n hodech mincí pravděpodobnost jevu B největší? Úlohy na pravděpodobnost Úloha 1. Jaká je pravděpodobnost, že při současném hodu 6 kostkami padne • na každé kostce jiné číslo, e) více než tři dvojky, • samé jedničky, f) právě ti dvojky, • alespoň tři dvojky, g) všechna čísla stejná. Úloha 2. Pravděpodobnost, že student A složí úspěšně zkoušku z Radiologické fyziky, je p, pravděpodobnost, že zkoušku složí student B, je q. Jaká je pravděpodobnost, že zkoušku složí • právě jeden ze studentů, c) oba studenti • alespoň jeden ze studentů, d) žádný ze studentů. Jednotlivé situace je třeba „šikovně“ vytvořit pomocí jevů nezávislých, resp. neslučitelných, resp. jevů obojího typu. Měření a zpracování dat aneb jak souvisí pravděpodobnost s měřením Měřené hodnoty veličin jsou „náhodné“ Vzpomeňme si na měření krevního tlaku – co znamená fakt, že při různých opakováních měření naměříme jiné hodnoty? Mění se tlak tak rychle, že jej fakticky nemůžeme určit, nebo lze z různých naměřených hodnot zjistit relevantní informaci? Je správné, že lékař měří pacientovi při dané návštěvě tlak pouze jednou? Řada veličin se řídí náhodnými vlivy, takže i za stejných podmínek mohou nabývat různých hodnot, popř. při jejich opakovaném měření můžeme dostat různé hodnoty. Náhodná veličina s diskrétním rozdělením Náhodná veličina X a její (diskrétní) rozdělení Veličina, která nabývá hodnot (x[1], x[2], …, x[N]) s pravděpodobnostmi (p[1], p[2], …, p[N]), p[1] + p[2] + …+ p[N ] = 1 . Víte proč ?? Uvědomte si neslučitelnost jevů, že veličina nabude dvou různých hodnot současně. Rozdělením náhodné veličiny rozumíme soubor všech dvojic (x[j], p[j]), pro j = 1, …, n.[] Bernoulliovo rozdělení - I Binomické a Poissonovo rozdělení – I Poissonovo rozdělení Limitní případ binomického rozdělení (Bernoulliova pro p = 1/2) pro velký počet pokusů, zajímáme-li se o velmi malý počet zdarů ve srovnání s počtem pokusů. Praktický případ Registrace radioaktivních částic v Geigerově-Millerově trubici. Cs^137 → Ba^137 + elektron + neutrino … asi 8% všech rozpadů Cs^137 → Ba^137^*^ + elektron + neutrino … asi 92% všech rozpadů Cs zdroj s aktivitou 10 μC (1 Curie … za 1 s rozpad 3,7 . 10^10 jader, v našem vzorku to znamená n = 3 700 000 pokusů za 10 sekund). V experimentu při počítání pulsů je nastaveno na cca 1 puls (zdar) za 1 s, počet zdarů v intervalu 10 s je tedy velmi malý proti n. Binomické a Poissonovo rozdělení – II Jak zjistit rozdělení experimentálně? Příklad se střelcem Střelec vystřelí n-krát na terč. Dosažené počty bodů při jednotlivých výstřelech představují hodnoty náhodné veličiny. Jaké jsou pravděpodobnosti jednotlivých hodnot? Pro n = 50 například: Při různých počtech výstřelů n se pravděpodobnosti budou obecně měnit. Pro rostoucí n budeme pozorovat jejich „ustalování“. Která hodnota nejlépe reprezentuje rozdělení? Jak máme zadávat náhodnou veličinu Náhodnou veličinu nejdokonaleji reprezentuje zadání jejího rozdělení. To je ovšem poněkud nepraktické. U střelce jsme viděli, že jeho kvalita je reprezentována hodnotou blízkou devítce. Realizovala se nejčastěji, má největší váhu. Reprezentativní hodnota (aritmetický průměr všech hodnot včetně „násobnosti“) Uvedená hodnota je váženým průměrem hodnot a nazývá se střední hodnotou náhodné veličiny. Úkol: Určete střední hodnotu dosažených bodů v příkladu se střelcem. Další charakteristiky rozdělení Který střelec je lepší? Dva střelci vystřelí n-krát na terč. Pro n = 50 máme jejich tabulky. Rozdělení pro oba střelce Rozptyl rozdělení Odvozené náhodné veličiny Y = f (X) náhodná veličina s rozdělením (y[j], p[j]) = (f(x[j]), p[j]), má-li veličina X rozdělení (x[j], p[j]) Která veličina charakterizuje „odchýlení“ hodnot od střední hodnoty? Určete střední hodnotu veličiny X - . Čekali jste tento výsledek? Rozptylem náhodné veličiny X rozumíme střední hodnotu náhodné veličiny Y = [X - ]^2 . Její odmocnina je tzv. směrodatná odchylka. Úkol: Vypočtěte hodnotu rozptylu a směrodatné odchylky u obou střelců. Jak byste interpretovali výsledek? Směrodatná odchylka vyjde 1,7 pro prvého a 1,2 pro druhého střelce. Medián rozdělení Distribuční funkce Funkce definovaná na R součtem pravděpodobností p[1]+ … + p[s] odpovídajících hodnotám menším než x[s+1]. Medián je hodnota x[s], pro kterou F(x[s]) < 0,5 a F(x[s+1]) ≥ 0,5. Úkol: Určete mediány rozdělení pro oba porovnávané střelce. Je výsledek očekávaný? Distribuční funkce pro střelce Náhodná veličina se spojitým rozdělením Náhodná veličina X a její (spojité) rozdělení Veličina, která nabývá všech reálných hodnot x z intervalu [x[m], x[M]] s elementárními pravděpodobnostmi dp = w (x) dx Rozdělením náhodné veličiny rozumíme funkci w(x) na intervalu [x[m], x[M]]. Též hustota pravděpodobnosti. Střední hodnota, rozptyl, distribuční funkce Úkol: Čemu je roven integrál z hustoty pravděpodobnosti (plocha pod grafem)? Důležitý příklad – normální rozdělení Fyzikální veličina a její chyba aneb co znamená zápis typu X = (2,51 ± 0,07) m Který výsledek je ten pravý? Měření délky ukazovátka - Předpokládejme, že všech n studentů v posluchárně bude měřit délku téhož ukazovátka, nebo ji bude jeden student měřit N-krát. - Budou všechny získané hodnoty stejné? - Proč se budou obecně lišit? - Která z naměřených hodnot je skutečnou délkou ukazovátka? Délka ukazovátka se při měření „chová“ jako náhodná veličina. Chyby, kterých se při měření dopouštíme - hrubé a systematické chyby (předpokládejme, že jsme je eliminovali) - náhodné chyby (jejich vlivem se budeme zabývat) Normální rozložení chyb - Předpokládejme, že existuje nějaká „správná“ hodnota délky ukazovátka x a že student naměřil hodnoty (x[1], x[2], …, x[N]), některé mohou být i stejné. - Odchylky od (zatím neznámé) správné hodnoty označme (ε[1], ε [2], …, ε [N]). Tyto hodnoty jsou hodnotami náhodné veličiny E. - Její hustotu pravděpodobnosti označme w(E). Za jistých podmínek je rozdělením normálním. - Předpokládejme, že odchylky jsou způsobeny m nezávislými vlivy, každý z nich odchýlí měřenou hodnotu od x o stejnou hodnotu α, kladnou nebo zápornou, s pravděpodobností 0,5. - Kladnou odchylku +α nazveme zdarem, zápornou (–α) nezdarem.Výsledná odchylka naměřené hodnoty x[i ]od x leží v intervalu (–mα, mα) a může nabývat pouze celých násobků α. Vliv chybových vlivů – I Vliv chybových vlivů – II Pravděpodobnost odchýlení o j kladných a m – j záporných vlivů (j kladných a m – j nezdarů), tj. pravděpodobnost vzniku odchylky ε = jα + (m – j)(–α) = ε = (2j – m)α je dána binomickým rozdělením (Bernoulliovým pro p =1/2). Pro velká m je lze nahradit rozdělením normálním (Gaussovým). To umožňuje následující zpracování výsledků. Aritmetický průměr a jeho chyba Reprezentativní hodnota měření aritmetický průměr všech naměřených hodnot (střední hodnota veličiny) Směrodatná odchylka příslušná aritmetickému průměru Správná hodnota veličiny nebude určena, ale s pravděpodobností 68,3 % leží v intervalu určeném aritmetickým průměrem a směrodatnou odchylkou takto: Krajní chyba … trojnásobek směrodatné odchylky … odpovídá pravděpodobnostnímu intervalu 97 % Různé typy průměrů aneb jen tak pro zajímavost Aritmetický průměr - I Příklad 1. Student měl ze tří matematických písemek v semestru tři hodnocení B a jedno C. U dvou závěrečných písemek měl A a D, u ústní zkoušky E. Jaká je jeho průměrná známka, jestliže všechny známky mají stejnou váhu? Aritmetický průměr - II Příklad 2. Řešte předchozí příklad za předpokladu, že závěrečná písemka má dvakrát větší váhu než průběžná a ústní zkouška má dvakrát větší váhu než závěrečná písemka. Průměrná rychlost – I Příklad 3. Automobil jel z A do B první úsek rychlostí 130 km/h stejnou dobu druhý úsek průměrnou rychlostí 70 km/h. Jaká byla jeho průměrná rychlost na celé trase? Je odpověď dána aritmetickým průměrem obou hodnot, tj. (130 + 70)/2 = 100 km / h? Je to věc definice. Průměrná je definována jako podíl celkové dráhy a celkové doby jízdy. Dráhu ale neznáme. Víme však, že oba úseky trvaly stejně času. Takže přece jen aritmetický průměr? Zkusme úlohu obměnit. Průměrná rychlost – II Příklad 4. Automobil jel z A do B první úsek rychlostí 130 km/h a druhý úsek rychlostí 70 km/h. Oba úseky byly stejně dlouhé. Jaká byla nyní průměrná rychlost? Jedná se o tzv. harmonický průměr. „Průměrný“ obdélník je čtverec Příklad 5. Určete stranu čtverce, který má stejný obsah jako obdélník o stranách a a b, nebo poloměr koule, která má stejný objem jako elipsoid o poloosách a , b , c. „Průměrný“ elipsoid je koule Výpočet: Jedná se o geometrický průměr. harmonický p. ≤ geometrický p. ≤ aritmetický p. A to ještě zdaleka nejsou všechny typy průměrů. Radiologická fyzika Radioaktivita