2018年6月10日日曜日

学習環境

数学読本〈6〉線形写像・1次変換/数論へのプレリュード/集合論へのプレリュード/εとδ/落ち穂拾い など(松坂 和夫(著)、岩波書店)の第26章(エピローグ - 落ち穂拾い、など)、26.4(確率分布と平均)、確率変数の平均または期待値、問2.を取り組んでみる。


  1. E a X + b = i = 1 n a x i + b p i = a i = 1 n x i p i + b i = 1 n p i = a E X + b

コード(Emacs)

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, Rational

a, b = symbols('a, b')
xs1 = [i for i in range(1, 7)]
xs2 = [a * x + b for x in xs1]
ps = [Rational(1, 6) for _ in range(6)]


def e(xs, ps):
    return sum([x * p for x, p in zip(xs, ps)])


for t in [a * e(xs1, ps) + b, e(xs2, ps)]:
    pprint(t)
    print()

入出力結果(Terminal, Jupyter(IPython))

$ ./sample2.py
7⋅a    
─── + b
 2     

7⋅a    
─── + b
 2     

$

0 コメント:

コメントを投稿