개요

파이톤으로 통계학에 입문한 8.2문제(p173)를 풀다.
8.2<랜덤 워킹> 확률 변수 X1,X2,...,Xn 독립, 확률 분포
P(Xi=1)=p, P(Xi=-1)=q (i=1,2,...,n)
따르다.하지만 p=1-q.
p=0.4시,X1+X2+...+X10,X1+X2+...+도표는 X20의 근사 확률 분포를 나타낸다.

컨디션

Jupter notebook을 사용합니다.

소스 코드

code

import numpy as np
import matplotlib.pyplot as plt
from sympy import *

def ranwal(k):
    ans=factorial(10)/(factorial(k)*factorial(10-k))*(0.4**k)*(0.6**(10-k))
    return ans #関数を設定

def data():
    result=[]
    for k in range(11):
        result.append(ranwal(k))
    return result #0~10の間で回す
plt.plot(data())

code2

import numpy as np
import matplotlib.pyplot as plt
from sympy import *

def ranwal(k):
    ans=factorial(20)/(factorial(k)*factorial(20-k))*(0.4**k)*(0.6**(20-k))
    return ans #関数を設定

def data():
    result=[]
    for k in range(21):
        result.append(ranwal(k))
    return result #0~20の間で回す
plt.plot(data())

결실

다음 도표를 얻었습니다.


잘못이 있으면 알려 주십시오.

참고 문헌

통계학 입문(동경대학 출판회)