2017年8月8日火曜日

学習環境

Head First Statistics (Dawn Griffiths (著)、黒川 利明 (翻訳)、木下 哲也 (翻訳)、黒川 洋 (翻訳)、黒川 めぐみ (翻訳)、オライリージャパン)の3章(ばらつきと広がりを測る)、プールパズル(p. 111)を取り組んでみる。

プールパズル(p. 111)

( xμ ) 2 n = ( x 2 2μx+ μ 2 ) n = x 2 n 2μ x n + μ 2 n = x 2 n 2μμ+ n μ 2 n = x 2 n μ 2

コード(Emacs)

Python 3

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

from sympy import pprint, symbols, summation

n = 10
i = symbols('i')
xs = [symbols(f'x_{i}') for i in range(1, n + 1)]
u = sum(xs) / n
eq1 = sum([(x - u) ** 2 for x in xs]) / n
eq2 = sum([x ** 2 for x in xs]) / n - u ** 2

for i,  eq in enumerate([eq1, eq2], 1):
    pprint(eq)

print(eq1.expand() == eq2.expand())

入出力結果(Terminal, IPython)

$ ./sample4.py
                                                      2                       
⎛  x₁   x₁₀   x₂   x₃   x₄   x₅   x₆   x₇   x₈   9⋅x₉⎞    ⎛  x₁   x₁₀   x₂   x
⎜- ── - ─── - ── - ── - ── - ── - ── - ── - ── + ────⎟    ⎜- ── - ─── - ── - ─
⎝  10    10   10   10   10   10   10   10   10    10 ⎠    ⎝  10    10   10   1
─────────────────────────────────────────────────────── + ────────────────────
                           10                                                 

                                  2                                           
₃   x₄   x₅   x₆   x₇   9⋅x₈   x₉⎞    ⎛  x₁   x₁₀   x₂   x₃   x₄   x₅   x₆   9
─ - ── - ── - ── - ── + ──── - ──⎟    ⎜- ── - ─── - ── - ── - ── - ── - ── + ─
0   10   10   10   10    10    10⎠    ⎝  10    10   10   10   10   10   10    
─────────────────────────────────── + ────────────────────────────────────────
       10                                                        10           

              2                                                         2     
⋅x₇   x₈   x₉⎞    ⎛  x₁   x₁₀   x₂   x₃   x₄   x₅   9⋅x₆   x₇   x₈   x₉⎞    ⎛ 
─── - ── - ──⎟    ⎜- ── - ─── - ── - ── - ── - ── + ──── - ── - ── - ──⎟    ⎜-
10    10   10⎠    ⎝  10    10   10   10   10   10    10    10   10   10⎠    ⎝ 
─────────────── + ─────────────────────────────────────────────────────── + ──
                                             10                               

                                                    2                         
 x₁   x₁₀   x₂   x₃   x₄   9⋅x₅   x₆   x₇   x₈   x₉⎞    ⎛  x₁   x₁₀   x₂   x₃ 
 ── - ─── - ── - ── - ── + ──── - ── - ── - ── - ──⎟    ⎜- ── - ─── - ── - ── 
 10    10   10   10   10    10    10   10   10   10⎠    ⎝  10    10   10   10 
───────────────────────────────────────────────────── + ──────────────────────
                         10                                                   

                                2                                             
  9⋅x₄   x₅   x₆   x₇   x₈   x₉⎞    ⎛  x₁   x₁₀   x₂   9⋅x₃   x₄   x₅   x₆   x
+ ──── - ── - ── - ── - ── - ──⎟    ⎜- ── - ─── - ── + ──── - ── - ── - ── - ─
   10    10   10   10   10   10⎠    ⎝  10    10   10    10    10   10   10   1
───────────────────────────────── + ──────────────────────────────────────────
     10                                                        10             

            2                                                         2       
₇   x₈   x₉⎞    ⎛  x₁   x₁₀   9⋅x₂   x₃   x₄   x₅   x₆   x₇   x₈   x₉⎞    ⎛  x
─ - ── - ──⎟    ⎜- ── - ─── + ──── - ── - ── - ── - ── - ── - ── - ──⎟    ⎜- ─
0   10   10⎠    ⎝  10    10    10    10   10   10   10   10   10   10⎠    ⎝  1
───────────── + ─────────────────────────────────────────────────────── + ────
                                           10                                 

                                                  2                           
₁   9⋅x₁₀   x₂   x₃   x₄   x₅   x₆   x₇   x₈   x₉⎞    ⎛9⋅x₁   x₁₀   x₂   x₃   
─ + ───── - ── - ── - ── - ── - ── - ── - ── - ──⎟    ⎜──── - ─── - ── - ── - 
0     10    10   10   10   10   10   10   10   10⎠    ⎝ 10     10   10   10   
─────────────────────────────────────────────────── + ────────────────────────
                       10                                                     

                            2
x₄   x₅   x₆   x₇   x₈   x₉⎞ 
── - ── - ── - ── - ── - ──⎟ 
10   10   10   10   10   10⎠ 
─────────────────────────────
  10                         
  2      2     2     2     2     2     2     2     2     2                    
x₁    x₁₀    x₂    x₃    x₄    x₅    x₆    x₇    x₈    x₉    ⎛x₁   x₁₀   x₂   
─── + ──── + ─── + ─── + ─── + ─── + ─── + ─── + ─── + ─── - ⎜── + ─── + ── + 
 10    10     10    10    10    10    10    10    10    10   ⎝10    10   10   

                                 2
x₃   x₄   x₅   x₆   x₇   x₈   x₉⎞ 
── + ── + ── + ── + ── + ── + ──⎟ 
10   10   10   10   10   10   10⎠ 
True
$

0 コメント:

コメントを投稿