Python은 Monte Carlo 방법을 구현합니다(코드 예)

import numpy as np

total = 8000000
count = 0

for i in range(total):
    x = np.random.rand()
    y = np.random.rand()
    dis = (x**2+y**2)**0.5
    if dis <p>3 위 프로그램에서는 전달을 위해 8,000,000개의 난수를 사용하므로 프로그램을 실행하는 데 일정 시간이 걸리고 최종 결과가 나옵니다. </p><p><span class="img-wrap"> <img src="https://img.php.cn//upload/image/549/675/172/1547436392667253.png" title="1547436392667253.png" alt="Python은 Monte Carlo 방법을 구현합니다(코드 예)"></span></p><p>4. 아래 그림은 제가 그리기 도구에서 임의로 그린 그림으로, 검은색 부분의 면적을 추정해 볼 수 있습니다. 그림에서. </p><p   style="max-width:90%"><span class="img-wrap"><img src="https://img.php.cn//upload/image/583/174/740/1547436406520588.png" title="1547436406520588.png" alt="Python은 Monte Carlo 방법을 구현합니다(코드 예)"></span></p><p>5. 위 그래프는 불규칙합니다. 난수를 많이 넣었을 때 검은색 부분에 난수가 나타날 확률만 알면 됩니다. 그런 다음 전체 면적을 곱하여 검은색을 추정합니다. 부분. 검은색 RGB 코드가 (0,0,0)이라는 것을 알고 있으므로 RGB 코드가 (0,0,0)일 때 난수가 나올 확률을 계산해야 합니다. </p><p>6. 코드 구현: </p><pre class="brush:php;toolbar:false">from PIL import Image
import numpy as np

im = Image.open("C:/Users/21974/Desktop/handwrite2.PNG")
total = 9000000
count = 0
defin = 0
width = im.size[0]
height = im.size[1]

for i in range(total):    #用蒙特卡罗方法获得估计值
    x = np.random.randint(0, width-1)
    y = np.random.randint(0, height-1)
    k = im.getpixel((x, y))
    if k[0]+k[1]+k[2] == 0:
        count += 1
print(int(width*height*count/total))

for i in range(width):    #用遍历获得准确值
    for j in range(height):
        k = im.getpixel((i, j))
        if k[0] + k[1] + k[2] == 0:
            defin += 1
print(defin)