【B++】在九个点之间跳动

B + +

已于 2023-04-22 11:53:40 修改

阅读量103

点赞数

文章标签：概率论

于 2023-04-21 20:59:11 首次发布

本文链接：https://blog.csdn.net/Barry_Wu_/article/details/130296700

版权

题设

现有九个点O=(0,0)，(0,1)，(0,2)，(1,0)，A=(1,1)，(1,2)，(2,0)，(2,1)，B=(2,2)，设由这九个点的坐标组成的集合为C。动点P从O点出发，每隔1s就会向上或向右跳一格，两种情况的概率都是 $\frac12$ 。当P=B时，它将停止跳动。

1. 原题

如果不允许P跳出集合C，求当P根据原规则从O跳到B的过程中P经过A的概率。

解：
设从某一点开始，P最终抵达点B且（不）经过点A的概率为W（W’），将跳出集合C（且经过点A）的概率为L（L’）。
列出九点概率矩阵：
$\begin{bmatrix}W'&W'&W'\\\frac12W\frac12W'&W&W'\\\frac12W\frac12W'&\frac12W\frac12W'&W'\end{bmatrix}$
因此，所求概率为 $P=\frac12$ 。
验证程序：

from random import randint as r
count = 1000000
imag = 1/2
n = 100
b, B, err = 0, 0, 0
print(f'准确值为{imag}')
for i in range (1,count+1):
    b = 0
    x,y=0,0
    while True:
        if x>2:
            x -= 1
            continue
        elif y>2:
            y -= 1
            continue
        elif (x,y) == (2,2):
            break
        elif (x,y) == (1,1):
            b=1
        if r(0,1):
            x+=1
        else:
            y+=1
    B += b
    if i%(count//n) == 0:
        ans = B/(count//n)
        err += ans - imag
        print(f'实际值为{ans},误差为{ans-imag:.4f}，累计误差为{err:.4f}')
        B = 0

2. 跃出时传送回原点

如果P∉C，则让P跳回原点，之后按原规则跳动。求当P根据原规则从O跳到B的过程中P经过A的概率。

解：
列出九点概率矩阵：
$\begin{bmatrix}\frac14W'\frac34L'&\frac12W'\frac12L'&W'\\\frac14W\frac18W'\frac14L\frac38L'&\frac12W\frac12L&\frac12W'\frac12L'\\\frac14W\frac18W'\frac14L\frac38L'&\frac14W\frac18W'\frac14L\frac38L'&\frac14W'\frac34L'\end{bmatrix}$
根据题意，L等价于W，W’等价于0，所以
$P_O=\frac12W\frac38L'$ ，所求概率 $P=W+L'W+L'^2W+⋯=\frac12∑_{n=0}^{+∞}(\frac38)^n=\frac12·\frac1{1-\frac38}=\frac45$ 。
特别地，如果不将P跳出C的情况计算在内，则所求概率为 $P=\frac{W_O}{W_O+W'_O}=\frac{\frac14}{\frac14+\frac18}=\frac23$ 。
验证程序：

from random import randint as r
count = 1000000
imag = 4/5
n = 100
b, B, err = 0, 0, 0
print(f'准确值为{imag}')
for i in range (1,count+1):
    b = 0
    x,y=0,0
    while True:
        if x>2 or y>2:
            x, y = 0, 0
            continue
        elif (x,y) == (2,2):
            break
        elif (x,y) == (1,1):
            b=1
        if r(0,1):
            x+=1
        else:
            y+=1
    B += b
    if i%(count//n) == 0:
        ans = B/(count//n)
        err += ans - imag
        print(f'实际值为{ans},误差为{ans-imag:.4f}，累计误差为{err:.4f}')
        B = 0

3. 跃出时传送到对面

P跳动时，令x’=(x+1)%3，y’=(y+1)%3。求当P从O跳到B的过程中P经过A的概率。

解：
根据题意，L等价于W，W’等价于0。列出化简后的九点概率矩阵：
$\begin{bmatrix}\frac34L'&\frac12L'&0\\\frac12W\frac38L'&W&\frac12L'\\\frac12W\frac38L'&\frac12W\frac38L'&\frac34L'\end{bmatrix}$
将P跳出C之后传送到(x,y)点的概率记为L_xy，则上述矩阵可写作：
$\begin{bmatrix}\frac12L_{00}\frac14L_{10}&\frac12L_{10}&0\\\frac12W\frac14L_{00}\frac18L_{10}&W&\frac12L_{01}\\\frac12W\frac14L_{00}\frac1{16}L_{01}\frac1{16}L_{10}&\frac12W\frac14L_{00}\frac18L_{01}&\frac12L_{00}\frac14L_{01}\end{bmatrix}$
由对称性可得L₁₀=L₀₁，因此由穿越映射关系（即P∉C时的跳跃规则）可得 $L_{00}=L_{01}=\frac12+\frac14L_{00}+\frac18L_{01}$ 。解之可得所求概率 $P=L_{00}=\frac45$ 。

验证程序：

from random import randint as r
count = 1000000
imag = 4/5n = 100b, B, err = 0, 0, 0
print(f'准确值为{imag}')
for i in range (1,count+1):
    b = 0
    x,y=0,0    
    while True:        
        if (x,y) == (2,2):            
            break       
        elif (x,y) == (1,1):
            b=1        
        if r(0,1):            
            x=(x+1)%3        
        else:            
            y=(y+1)%3    
    B += b    
    if i%(count//n) == 0:        
        ans = B/(count//n)        
        err += ans - imag        
        print(f'实际值为{ans},误差为{ans-imag:.4f}，累计误差为{err:.4f}')        
        B = 0

4. 跃出时传送到对称点

P跳出C时，将P传送到上一步所在点与A成中心对称的点。求当P从O跳到B的过程中P经过A的概率。

解：
仿照上一题，可列出概率矩阵： $\begin{bmatrix}\frac12L_{20}\frac14L_{10}&\frac12L_{10}&0\\\frac12W\frac14L_{20}\frac18L_{10}&W&\frac12L_{01}\\\frac12W\frac18L_{02}\frac18L_{20}\frac1{16}L_{01}\frac1{16}L_{10}&\frac12W\frac14L_{02}\frac18L_{01}&\frac12L_{02}\frac14L_{01}\end{bmatrix}$
按对称性将其化简：
$\begin{bmatrix}\frac12L_{02}\frac14L_{01}&\frac12L_{01}&0\\\frac12W\frac14L_{02}\frac18L_{01}&W&\frac12L_{01}\\\frac12W\frac14L_{02}\frac18L_{01}&\frac12W\frac14L_{02}\frac18L_{01}&\frac12L_{02}\frac14L_{01}\end{bmatrix}$
按穿越映射关系列出方程组：
$\left\{\begin{matrix}L_{01}=\frac12+\frac14L_{02}+\frac18L_{01}\\L_{02}=\frac12L_{02}+\frac14L_{01}\end{matrix}\right.$
解得 $L_{01}=\frac23,L_{02}=\frac13$ 。因此，所求概率 $P=L_{00}=\frac12+\frac14L_{02}+\frac18L_{01}=\frac23$ 。
验证程序：

from random import randint as r
count = 1000000
imag = 2/3
n = 100
b, B, err = 0, 0, 0
print(f'准确值为{imag}')
for i in range (1,count+1):
    b = 0
    x,y=0,0
    while True:
        if (x,y) == (1,1):
            b=1
        if (x,y) == (2,2):
            break
        if x>2:
            x = 0
            y = 2-y
        elif y>2:
            y = 0
            x = 2-x
        elif r(0,1):
            x+=1
        else:
            y+=1
    B += b
    if i%(count//n) == 0:
        ans = B/(count//n)
        err += ans - imag
        print(f'实际值为{ans},误差为{ans-imag:.4f}，累计误差为{err:.4f}')
        B = 0

5. 跃出时传送到随机点

如果P∉C，则让P随机地跳到C中的任意一点，且跳到每个点的概率都是 $\frac19$ ，之后按原规则跳动。求当P根据原规则从O跳到B的过程中P经过A的概率。

解：
列出九点概率矩阵：
$\begin{bmatrix}\frac14W'\frac34L'&\frac12W'\frac12L'&W'\\\frac14W\frac18W'\frac14L\frac38L'&\frac12W\frac12L&\frac12W'\frac12L'\\\frac14W\frac18W'\frac14L\frac38L'&\frac14W\frac18W'\frac14L\frac38L'&\frac14W'\frac34L'\end{bmatrix}$
九个点的总概率式为 $P_C=\frac5{36}W\frac{23}{72}W'\frac5{36}L\frac{29}{72}L'$ 。
根据题意，L等价于W，W’等价于0，所以 $P_O=\frac12W\frac38L',P_C=\frac5{18}W\frac{29}{72}L'$ 。
因此，所求概率 $P=W_O+L'_O(W+L'W+L'^2W+⋯)=\frac12+\frac38∑_{n=0}^{+∞}\frac5{18}(\frac{29}{72})^n=\frac12+\frac38·\frac5{18}·\frac1{1-\frac{29}{72}}=\frac{29}{43}$ 。

验证程序：

from random import randint as r
count = 1000000
imag = 29/43
n = 100
b, B, err = 0, 0, 0
print(f'准确值为{imag}')
for i in range (1,count+1):
    b = 0
    x,y=0,0
    while True:
        if x>2 or y>2:
            x, y = r(0,2), r(0,2)
            continue
        elif (x,y) == (2,2):
            break
        elif (x,y) == (1,1):
            b=1
        if r(0,1):
            x+=1
        else:
            y+=1
    B += b
    if i%(count//n) == 0:
        ans = B/(count//n)
        err += ans - imag
        print(f'实际值为{ans},误差为{ans-imag:.4f}，累计误差为{err:.4f}')
        B = 0

B + +

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
【B++】在九个点之间跳动

现有九个点O=(0,0)，(0,1)，(0,2)，(1,0)，A=(1,1)，(1,2)，(2,0)，(2,1)，B=(2,2)，设由这九个点的坐标组成的集合为C。动点P从O点出发，每隔1s就会向上或向右跳一格，两种情况的概率都是。当P=B时，它将停止跳动。
复制链接

扫一扫