1. Učit se
    2. /
  2. Kurzy
    3. /
  3. Statistické simulace v Pythonu
Connected

cvičení

Výpočet hodnoty čísla pí

Teď si vyzkoušíme klasický příklad – odhad hodnoty \(\pi\).

Představ si čtverec se stranou délky \(2\), jehož střed leží v počátku souřadnic \((0, 0)\) a jehož rohové body mají souřadnice \((1, 1), (1, -1), (-1, 1), (-1, -1)\). Obsah tohoto čtverce je \(2\times 2 = 4\). Dovnitř čtverce je dokonale vepsána kružnice s poloměrem \(1\) se středem v počátku. Obsah kruhu je \(\pi \times \text{poloměr}^2 = \pi\).

Pro odhad \(\pi\) náhodně generujeme body uvnitř čtverce a sledujeme, jaký podíl z nich leží uvnitř kružnice (\(x^2 + y^2 <= 1\)). Obsah kruhu pak odpovídá \(4\)-násobku tohoto podílu – a to je náš odhad \(\pi\).

Po dokončení tohoto cvičení budeš vědět, jak simulaci využít k numerickým výpočtům.

Pokyny

100 XP
  • Podívej se na skutečnou hodnotu \(\pi\) pomocí np.pi v konzoli. Inicializuj proměnnou sims na hodnotu 10000 a circle_points na 0.
  • Uvnitř cyklu for vygeneruj bod (souřadnice x a y) pomocí np.random.uniform() v rozsahu od -1 do 1 s parametrem size=2.
  • Ověř, zda bod leží uvnitř jednotkové kružnice podle rovnice \(x^2 + y^2 <= 1\), výsledek ulož do proměnné within_circle a podle toho aktualizuj hodnotu circle_points.
  • Vypiš odhad \(\pi\) jako pi_sim – tedy 4-násobek podílu bodů ležících uvnitř kružnice.