Simulering av antall terningkast for å få samme antall øyne i to kast på rad

Simulering av antall terningkast for å få samme antall øyne i to kast på rad

Ane har en vanlig sekssidet terning. Hun ønsker å finne ut hvor mange ganger hun i gjennomsnitt må kaste terningen for å få det samme antallet øyne i to kast på rad.

Hun har laget tabellen nedenfor.

Tabell 1: Sannsynlighet for at et kast er nødvendig
Kast nummer123456
Sannsynlighet for at kastet er nødvendig1156\frac{5}{6}(56)2\left(\frac{5}{6}\right)^2(56)3\left(\frac{5}{6}\right)^3(56)4\left(\frac{5}{6}\right)^4

Forklar at

1+1+56+(56)2+(56)3+1+1+\frac{5}{6}+\left(\frac{5}{6}\right)^2+\left(\frac{5}{6}\right)^3+\ldots

vil gi det forventede antallet kast Ane må gjøre for å få det samme antallet øyne i to kast på rad. Bestem denne verdien.

Bruk simulering til å bestemme forventningsverdien til summen av antall øyne Ane vil få på terningen i kastene hun bruker for å få det samme antallet øyne i to kast på rad.

Fasit

24,5

Løsningsforslag

I denne oppgaven er jeg veldig usikker på hva som kreves for å forklare at uttrykket i oppgaveteksten er det samme som forventningsverdien. Jeg tror ikke det er meningen at elever skal gjøre det samme som jeg har gjort her – men jeg klarer ikke helt å se en enklere måte å argumentere for at forventningsverdien er eksakt lik summen av «antall kast nødvendig».

Valgtre til oppgave 6 del 2

Vi lar XX være antall kast som trengs før vi har fått 2 like terningkast på rad. Sannsynligheten for å at et terningkast har samme antall øyne som det forrige er 1/61/6, og sannsynligheten for at antall øyne er ulikt er 5/65/6. Vi kan bruke multiplikasjonsprinsippet og sette opp følgende sannsynlighetsfordeling for XX:

xix_{i}P(xi)P(x_{i})
1100
2216\frac{1}{6}
335616\frac{5}{6} \cdot \frac{1}{6}
44(56)216\left( \frac{5}{6} \right)^{2} \cdot \frac{1}{6}
55(56)316\left( \frac{5}{6} \right)^{3} \cdot \frac{1}{6}

Forventningsverdien til XX vil da være

E(X)=limni=1nxiP(xi)=10+216+35616+4(56)216+5(56)316+=16(2(56)0+3(56)1+4(56)2+5(56)3+S)\begin{aligned} \text{E}(X)&=\lim_{ n \to \infty } \sum_{i=1}^n x_{i} \cdot P(x_{i})\\ &=1 \cdot 0 + 2 \cdot \frac{1}{6} + 3 \cdot \frac{5}{6} \cdot \frac{1}{6} +4 \cdot \left( \frac{5}{6} \right)^{2} \cdot \frac{1}{6}+5 \cdot \left( \frac{5}{6} \right)^{3} \cdot \frac{1}{6} + \cdots \\ &= \frac{1}{6} \left(\underbrace{ 2\left( \frac{5}{6} \right)^0 + 3\left( \frac{5}{6} \right)^{1} + 4\left( \frac{5}{6} \right)^{2} + 5\left( \frac{5}{6} \right)^{3} + \cdots }_{ S } \right) \end{aligned}

Vi kaller alt inni parentesen for SS, og omskriver heltallene som står foran 56\frac{5}{6} som en sum av enere:

S=(1+1)2(56)0+(1+1+1)3(56)1(1+1+1+1)4(56)2+\begin{aligned} S = \underbrace{ (\textcolor{seagreen}{1}+\textcolor{steelblue}{1}) }_{ 2 } \cdot \left( \frac{5}{6} \right)^{0} + \underbrace{ ( \textcolor{seagreen}{1} + \textcolor{steelblue}{1}+\textcolor{orange}{1}) }_{ 3 } \cdot \left( \frac{5}{6} \right)^{1} \underbrace{ ( \textcolor{seagreen}{1} + \textcolor{steelblue}{1}+\textcolor{orange}{1} + \textcolor{tomato}{1}) }_{ 4 } \cdot \left( \frac{5}{6} \right)^{2} + \dots \end{aligned}

Vi deler nå opp denne summen i en rekke delsummer slik at S=limnS1+S2++SnS = \lim_{ n \to \infty } S_{1} + S_{2} + \dots + S_{n} hvor

S1=1(56)0+1(56)1+1(56)2+=1156=116=6S2=1(56)0+1(56)1+1(56)2+=1156=116=6S3=1(56)1+1(56)2+1(56)3+56156=5616=5S4=1(56)2+1(56)3+1(56)4+=(56)2156=(56)216=6(56)2S5=1(56)3+1(56)4+1(56)5+=(56)3156=(56)316=6(56)3\begin{aligned} \color{seagreen} S_{1} & = \textcolor{seagreen}{1} \cdot \left( \frac{5}{6} \right)^{0} + \textcolor{seagreen}{1} \cdot \left( \frac{5}{6} \right)^{1}+ \textcolor{seagreen}{1} \cdot \left( \frac{5}{6} \right)^{2} + \dots = \frac{1}{1-\frac{5}{6}}=\frac{1}{\frac{1}{6}}=\textcolor{seagreen}{6}\\ \color{steelblue} S_{2} & = \textcolor{steelblue}{1} \cdot \left( \frac{5}{6} \right)^{0} + \textcolor{steelblue}{1} \cdot \left( \frac{5}{6} \right)^{1}+ \textcolor{steelblue}{1} \cdot \left( \frac{5}{6} \right)^{2} + \dots = \frac{1}{1-\frac{5}{6}}=\frac{1}{\frac{1}{6}}= \textcolor{steelblue}{6}\\ \color{orange} S_{3} & = \textcolor{orange}{1} \cdot \left( \frac{5}{6} \right)^{1} + \textcolor{orange}{1} \cdot \left( \frac{5}{6} \right)^{2}+ \textcolor{orange}{1} \cdot \left( \frac{5}{6} \right)^{3} + \dots \frac{\frac{5}{6}}{1-\frac{5}{6}}=\frac{\frac{5}{6}}{\frac{1}{6}}=\textcolor{orange}{5}\\ \color{tomato} S_{4} & = \textcolor{tomato}{1} \cdot \left( \frac{5}{6} \right)^{2} + \textcolor{tomato}{1} \cdot \left( \frac{5}{6} \right)^{3}+ \textcolor{tomato}{1} \cdot \left( \frac{5}{6} \right)^{4} + \dots = \frac{\left( \frac{5}{6} \right)^{2}}{1-\frac{5}{6}}=\frac{\left( \frac{5}{6} \right)^{2}}{\frac{1}{6}}=\textcolor{tomato}{6 \cdot \left( \frac{5}{6} \right)^{2}}\\ \color{maroon} S_{5} & = \textcolor{maroon}{1} \cdot \left( \frac{5}{6} \right)^{3} + \textcolor{maroon}{1} \cdot \left( \frac{5}{6} \right)^{4}+ \textcolor{maroon}{1} \cdot \left( \frac{5}{6} \right)^{5} + \dots = \frac{\left( \frac{5}{6} \right)^{3}}{1-\frac{5}{6}}=\frac{\left( \frac{5}{6} \right)^{3}}{\frac{1}{6}}=\textcolor{maroon}{6 \cdot \left( \frac{5}{6} \right)^{3}} \end{aligned}

Forventningsverdien er altså

E(X)=16S=16(S1+S2+S3+S4+S5+)=16(6+6+5+6(56)2+6(56)3+)=1+1+56+(56)2+(56)3+\begin{aligned} \text{E}(X)&=\frac{1}{6}S \\ &=\frac{1}{6}\left( \textcolor{seagreen}{S_{1}}+\textcolor{steelblue}{S_{2}} + \textcolor{orange}{S_{3}} + \textcolor{tomato}{S_{4}} + \textcolor{maroon}{S_{5}} + \cdots \right) \\ &=\frac{1}{6} \left( \textcolor{seagreen}{6} + \textcolor{steelblue}{6} + \textcolor{orange}{5} + \textcolor{tomato}{6 \cdot \left( \frac{5}{6} \right)^{2} } + \textcolor{maroon}{6 \cdot \left( \frac{5}{6} \right)^{3}} + \cdots \right) \\ &= \textcolor{seagreen}{1} + \textcolor{steelblue}{1} + \textcolor{orange}{\frac{5}{6}}+\textcolor{tomato}{\left( \frac{5}{6} \right)^{2}} + \textcolor{maroon}{\left( \frac{5}{6} \right)^{3}} + \dots && \blacksquare \end{aligned}

Hvis vi ser bort fra det aller første leddet (1\textcolor{seagreen}{1}), så er dette en uendelig geometrisk rekke med a1=1a_{1}=1 og k=56k=\frac{5}{6}

s=1+56+(56)2+s= \textcolor{steelblue}{1}+ \textcolor{orange}{\frac{5}{6}}+ \textcolor{tomato}{\left( \frac{5}{6} \right)^{2}} + \cdots

Vi kan finne summen av rekka ss med GeoGebra, eller med formelen for sum av uendelig geometrisk rekke:

s=1156=116=16166=6s=\frac{1}{1-\frac{5}{6}}= \frac{1}{\frac{1}{6}}= \frac{1\cdot 6}{\frac{1}{6}\cdot 6}= 6

Til sammen blir altså E(X)=1+s=1+6=7\text{E}(X)=\textcolor{seagreen}{1}+s=\textcolor{seagreen}{1}+6=7.

Verdien av rekka er 7.

Vi skal simulere forventningsverdien til summen av antall øyne på alle terningene som kastes i jakten på å få to like kast på rad.

from random import randint
N = 100_000                        
sum_øyne = 0                        # totalt antall øyne på terningene

for i in range(N):
    t1 = randint(1,6)               # terningkast 1
    t2 = randint(1,6)               # terningkast 2

    sum_øyne = sum_øyne + t1 + t2   # legger til resultatene til summen
    while t1 != t2:
        t1 = t2                     # flytter t2's verdi til t1
        t2 = randint(1,6)           # ruller t2 på nytt
        sum_øyne = sum_øyne + t2    # legger til nytt resultat til summen

EX = sum_øyne/N                     # forventningsverdi = snitt i det lange løp
print(f"Jeg estimerer forventningsverdien til å være {EX:.3f} etter {N} simuleringer.")

Output: Jeg estimerer forventningsverdien til å være 24.502 etter 100000 simuleringer.

Etter å ha kjørt programmet flere ganger ser det ut til estimatet mitt er stabilt på rundt 24,524{,}5. Det stemmer også godt med at forventningsverdien for en terning er 3,53{,}5 og vi trenger i snitt 77 kast før vi har fått to like på rad.

Jeg estimerer forventningsverdien til summen av antall øyne før Ane får to like terninger på rad til å være 24,5.

Sensorveiledning

1 poeng for å forklare rekken og 1 poeng for å regne ut verdien.

Simulering med en god strategi, men feil svar, kan få 1 poeng. Simulering med riktig strategi, men med små feil som fører til feil svar, kan få 2 poeng. Analytisk løsning med god argumentasjon kan gi 1 poeng.