Перейти к содержанию
Форум химиков на XuMuK.ru.

Являются ли ферменты вечными двигателями второго рода?


Рекомендуемые сообщения

🚑 Решение задач, контроши, рефераты, курсовые и другое! Онлайн сервис помощи учащимся. Цены в 2-3 раза ниже! 200 руб. на 1-й заказ по коду vsesdal143982
В 17.12.2025 в 22:42, Мансур Фазлуллин сказал:

Увеличил геометрические размеры модели. Предыдущая модель была размером 12 на 16 клеток. Увеличенная 14 на 20 клеток.

Для увеличения модели пришлось усложнить программный код, что привело к снижению скорости расчётов, что не очень удобно.

Поэтому с помощью ИИ провёл рефакторинг кода, что позволило заметно увеличить скорость вычислений, изменения после строки "while t < 1000:"

Код после рефакторинга.

 

import random
import matplotlib.pyplot as plt

n = 0

n_out = 0

x_n = list(range(21))

y_n = []

for i in x_n:
    y_n += [0]

x_v = list(range(21))

y_v = []

for i in x_v:
    y_v += [14]

y_l = list(range(15))

x_l = []

for i in y_l:
    x_l += [0]
    
y_p = list(range(15))

x_p = []

for i in y_p:
    x_p += [20]


x_b = [6, 6, 6, 6, 7, 8, 9, 10, 11, 12, 13, 14, 14, 14, 14]

y_b = [0, 1, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 1, 0] 

x_b1 = [5, 5, 5, 5, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 15, 15, 15, 15]

y_b1 = [0, 1, 2, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 2, 1, 0]

fnx_l = [9, 8, 7, 6, 5, 6, 7, 8, 9, 10]

fny_l = [4, 5, 6, 7, 8, 8, 7, 6, 5, 4]

fnx_r = [10, 11, 12, 13, 14, 15, 14, 13, 12, 11]

fny_r = [4, 5, 6, 7, 8, 8, 7, 6, 5, 4]

P = 0
k = True

while k:
    P += 1
    t = 0
    x = 1
    y = 1
    xn = x
    yn = y
    x1 = []
    y1 = []


    while t < 1000:  
        t += 1

        z = random.randint(1, 100)

        # Условия для движения вправо при z <= 25

        if z <= 25:
            if x < 20:
                x += 1

            if (x == 6 and y <=3) or (x == 15 and y <=3):
                x -= 1

            if x == 10 and y == 4:
                x -= 1

            if x == 9 and y == 5:
                y -= 1

            if x == 8 and y == 6:
                y -= 1

            if x == 7 and y == 7:
                y -= 1

            if x == 11 and y == 4:
                y += 1

            if x == 12 and y == 5:
                y += 1

            if x == 13 and y == 6:
                y += 1

            if x == 14 and y == 7:
                y += 1

        # Условия для движения влево при 26 <= z <= 50

        elif 26 <= z <= 50:
            if x > 0:
                x -= 1

            if (x == 5 and y <= 3) or (x == 14 and y <= 3):
                x += 1

            if x == 9 and y == 4:
                y += 1

            if x == 8 and y == 5:
                y += 1

            if x == 7 and y == 6:
                y += 1

            if x == 6 and y == 7:
                y += 1

            if x == 10 and y == 4:
                x += 1

            if x == 11 and y == 5:
                y -= 1

            if x == 12 and y == 6:
                y -= 1

            if x == 13 and y == 7:
                y -= 1

        # Условия для движения вверх 51 <= z <= 73

        elif 51 <= z <= 73:
            if y < 14:
                y += 1

            if (y == 4 and 6 <= x < 10) or (y == 4 and 10 < x <= 14):
                y -= 1

            if x == 9 and y == 5:
                x -= 1

            if x == 8 and y == 6:
                x -= 1

            if x == 7 and y == 7:
                x -= 1

            if x == 6 and y == 8:
                x -= 1

            if x == 11 and y == 5:
                x += 1

            if x == 12 and y == 6:
                x += 1

            if x == 13 and y == 7:
                x += 1

            if x == 14 and y == 8:
                x += 1

        # Условия для предотвращения проскока сквозь воронку

        elif x == 9 and y == 4 and 74 <= z <= 100:
            z = 101

        elif x == 11 and y == 4 and 74 <= z <= 100:
            z = 101        

        # Условия для движения вниз при 74 <= z <= 100


        elif 74 <= z <= 100:
            if y > 0:
                y -= 1

            if (y == 3 and 6 <= x < 10) or (y == 3 and 10 < x <= 14):
                y += 1

            if x == 9 and y == 4:
                x += 1

            if x == 8 and y == 5:
                x += 1

            if x == 7 and y == 6:
                x += 1

            if x == 6 and y == 7:
                x += 1

            if x == 11 and y == 4:
                x -= 1

            if x == 12 and y == 5:
                x -= 1

            if x == 13 and y == 6:
                x -= 1

            if x == 14 and y == 7:
                x -= 1


        x1.append(x)
        y1.append(y)


    if 6 <= x1[-1] <= 14 and y1[-1] <= 3:
        n += 1

    if  x1[-1] <= 5 and y1[-1] <=3:
        n_out += 1

    if 15 <= x1[-1] and y1[-1] <=3:
        n_out += 1
     

    if P == 10000:
        k = False


print('n=', n)

print('n/P', n/P)

print('n_out=', n_out)

print('n_out/P', n_out/P)

plt.grid(True)
plt.plot(x_n, y_n, 'k')
plt.plot(x_v, y_v, 'k')
plt.plot(x_l, y_l, 'k')
plt.plot(x_p, y_p, 'k')
plt.plot(x_b, y_b, 'k')
plt.plot(x_b1, y_b1, 'k')
plt.plot(fnx_l, fny_l, 'k')
plt.plot(fnx_r, fny_r, 'k')
plt.plot(x1, y1, 'b')
plt.plot(x1[-1], y1[-1], 'ro')
plt.plot(xn, yn, 'go')
plt.show()    

 

Ссылка на комментарий
  • 2 недели спустя...

Прежде чем утверждать, что результаты моделирования опровергли предположение, что для того чтобы концентрация частиц внутри бутылки превысила концентрацию частиц снаружи бутылки потребуется бесконечное количество времени, следует посмотреть как быстро будет меняться отношение концентраций при увеличении размеров модели.

 

При увеличении размеров модели, процесс роста разницы концентраций частиц замедлился, что ожидаемо. При этом какого-то заметного роста разницы концентраций по сравнению с первоначальной моделью не произошло. Т.е. увеличение размеров воронки не приводит к увеличению разницы концентраций.

При этом увеличение размеров модели, открывает новые возможности для поиска более эффективных форм способных загонять броуновские частицы в бутылку, использую энергию теплового движения атомов и молекул.

Например, модель с воронкой можно дополнить "зонтиком".

 

Figure_1.png.567b52fd65a27d2d2036523d26ad072f.png

Идея следующая. Зонтик будет сильнее мешать частицам выбираться из воронки чем попасть в неё.

Корректность данного предположения была проверена методом моделирования.

Результаты представлены в виде графика.

 

.thumb.png.ff463f25b2b8f2fac965b51b3e93bf80.png

 

Синий график - результаты моделирования для первоначальной модели 12х16.

Оранжевый график - результаты моделирования для увеличенной модели 14х20.

Жёлтый график - результаты моделирования для увеличенной модели 14х20 дополненной зонтиком.

Расчёты выполнялись для частиц, которые движутся вниз при z от 74 до 100.

Как изменятся результаты при моделировании движения более мелких частиц, которые движутся вниз при z от 75 до 100?

 

Программа для модели с зонтиком и результаты расчётов ниже.

 

import random
import matplotlib.pyplot as plt

n = 0

n_out = 0

x_n = list(range(21))

y_n = []

for i in x_n:
    y_n += [0]

x_v = list(range(21))

y_v = []

for i in x_v:
    y_v += [14]

y_l = list(range(15))

x_l = []

for i in y_l:
    x_l += [0]
    
y_p = list(range(15))

x_p = []

for i in y_p:
    x_p += [20]


x_b = [6, 6, 6, 6, 7, 8, 9, 10, 11, 12, 13, 14, 14, 14, 14]

y_b = [0, 1, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 2, 1, 0] 

x_b1 = [5, 5, 5, 5, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 15, 15, 15, 15]

y_b1 = [0, 1, 2, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 2, 1, 0]

fnx_l = [9, 8, 7, 6, 5, 6, 7, 8, 9, 10]

fny_l = [4, 5, 6, 7, 8, 8, 7, 6, 5, 4]

fnx_r = [10, 11, 12, 13, 14, 15, 14, 13, 12, 11]

fny_r = [4, 5, 6, 7, 8, 8, 7, 6, 5, 4]

x_lbr = [8, 9, 10, 11, 12, 11, 10, 9, 8]

y_lbr = [7, 8, 9, 8, 7, 7, 8, 7, 7]

P = 0
k = True

while k:
    P += 1
    t = 0
    x = 1
    y = 1
    xn = x
    yn = y
    x1 = []
    y1 = []


    while t < 1000:  
        t += 1

        z = random.randint(1, 100)

        # Условия для движения вправо при z <= 25

        if z <= 25:
            if x < 20:
                x += 1

            if (x == 6 and y <=3) or (x == 15 and y <=3):
                x -= 1

            if x == 10 and y == 4:
                x -= 1

            if x == 9 and y == 5:
                y -= 1

            if x == 8 and y == 6:
                y -= 1

            if x == 7 and y == 7:
                y -= 1

            if x == 11 and y == 4:
                y += 1

            if x == 12 and y == 5:
                y += 1

            if x == 13 and y == 6:
                y += 1

            if x == 14 and y == 7:
                y += 1

            # изменения для лабиринта

            if x == 10 and y == 8:
                y += 1

            if x == 11 and y == 8:
                y -= 1

        # Условия для движения влево при 26 <= z <= 50

        elif 26 <= z <= 50:
            if x > 0:
                x -= 1

            if (x == 5 and y <= 3) or (x == 14 and y <= 3):
                x += 1

            if x == 9 and y == 4:
                y += 1

            if x == 8 and y == 5:
                y += 1

            if x == 7 and y == 6:
                y += 1

            if x == 6 and y == 7:
                y += 1

            if x == 10 and y == 4:
                x += 1

            if x == 11 and y == 5:
                y -= 1

            if x == 12 and y == 6:
                y -= 1

            if x == 13 and y == 7:
                y -= 1

            # изменения для лабиринта

            if x == 10 and y == 8:
                y += 1

            if x == 9 and y == 8:
                y -= 1

        # Условия для движения вверх 51 <= z <= 73

        elif 51 <= z <= 73:
            if y < 14:
                y += 1

            if (y == 4 and 6 <= x < 10) or (y == 4 and 10 < x <= 14):
                y -= 1

            if x == 9 and y == 5:
                x -= 1

            if x == 8 and y == 6:
                x -= 1

            if x == 7 and y == 7:
                x -= 1

            if x == 6 and y == 8:
                x -= 1

            if x == 11 and y == 5:
                x += 1

            if x == 12 and y == 6:
                x += 1

            if x == 13 and y == 7:
                x += 1

            if x == 14 and y == 8:
                x += 1

            # изменения для лабиринта            

            if x == 9 and y == 8:
                x += 1

            if x == 10 and y == 9:
                y -= 1

            if x == 11 and y == 8:
                x -= 1

        # Условия для предотвращения проскока сквозь воронку

        elif x == 9 and y == 4 and 74 <= z <= 100:
            z = 101

        elif x == 11 and y == 4 and 74 <= z <= 100:
            z = 101        

        # Условия для движения вниз при 74 <= z <= 100


        elif 74 <= z <= 100:
            if y > 0:
                y -= 1

            if (y == 3 and 6 <= x < 10) or (y == 3 and 10 < x <= 14):
                y += 1

            if x == 9 and y == 4:
                x += 1

            if x == 8 and y == 5:
                x += 1

            if x == 7 and y == 6:
                x += 1

            if x == 6 and y == 7:
                x += 1

            if x == 11 and y == 4:
                x -= 1

            if x == 12 and y == 5:
                x -= 1

            if x == 13 and y == 6:
                x -= 1

            if x == 14 and y == 7:
                x -= 1

            # изменения для лабиринта. 

            if x == 9 and y == 7:
                x -= 1
                
            # При движении вниз из точки(10; 9) частица с равной вероятностью может
            # двигаться влево и вправо. Принял, что частица движется только вправо.      

     

            if x == 10 and y == 8:
                x += 1                

            if x == 11 and y == 7:
                x += 1


        x1.append(x)
        y1.append(y)


    if 6 <= x1[-1] <= 14 and y1[-1] <= 3:
        n += 1

    if  x1[-1] <= 5 and y1[-1] <=3:
        n_out += 1

    if 15 <= x1[-1] and y1[-1] <=3:
        n_out += 1
     

    if P == 1:
        k = False


print('n=', n)

print('n/P', n/P)

print('n_out=', n_out)

print('n_out/P', n_out/P)

plt.grid(True)
plt.plot(x_n, y_n, 'k')
plt.plot(x_v, y_v, 'k')
plt.plot(x_l, y_l, 'k')
plt.plot(x_p, y_p, 'k')
plt.plot(x_b, y_b, 'k')
plt.plot(x_b1, y_b1, 'k')
plt.plot(fnx_l, fny_l, 'k')
plt.plot(fnx_r, fny_r, 'k')
plt.plot(x_lbr, y_lbr, 'k')
plt.plot(x1, y1, 'b')
plt.plot(x1[-1], y1[-1], 'ro')
plt.plot(xn, yn, 'go')
plt.show()    

 

Результаты расчётов для увеличенной модели 14х20:

74-100
(14)
         
С воронкой          
в бутылке          
№ серии 20 50 75 100 125
1 0 2 18 54 85
2 0 4 16 28 74
3 0 2 15 43 74
4 0 2 8 35 75
5 0 0 17 46 63
           
среднее 0,0000 2 14,8 41,2 74,2
среднее на
столбик
0 0,22 1,64 4,58 8,24
           
           
С воронкой          
вне бутылки          
№ серии 20 50 75 100 125
1 8161 6649 5886 5419 5091
2 8126 6631 5892 5360 5106
3 8122 6533 5879 5362 5079
4 8122 6613 5905 5466 5012
5 8160 6688 6053 5540 5159
           
среднее 8138,2 6622,8 5923 5429,4 5089,4
среднее на
столбик
678,18 551,9 493,58 452,45 424,12
           
  20 50 75 100 125
Отношение
концентр.
0,0000 0,0004 0,0033 0,0101 0,0194

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 150 175 200 250 300
1 147 186 233 373 508
2 132 182 228 389 497
3 144 183 245 392 533
4 127 189 228 382 480
5 139 189 226 380 517
           
среднее 137,8 185,8 232 383,2 507
среднее на
столбик
15,31 20,64 25,78 42,58 56,33
           
           
С воронкой,          
вне бутылки          
№ серии 150 175 200 250 300
1 4794 4606 4360 4117 3919
2 4813 4610 4438 4149 3877
3 4881 4583 4371 4134 3904
4 4714 4544 4465 4093 3943
5 4882 4542 4405 4222 3910
           
среднее 4816,8 4577 4407,8 4143 3910,6
среднее на
столбик
401,4 381,42 367,32 345,25 325,88
           
  150 175 200 250 300
Отношение
концентр.
0,0381 0,0541 0,0702 0,1233 0,1729

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 350 400 450 500 550
1 659 825 933 1016 1199
2 691 829 906 1046 1171
3 639 804 922 1017 1174
4 672 844 923 1081 1143
5 723 773 975 1074 1120
           
среднее 676,8 815 931,8 1046,8 1161,4
среднее на
столбик
75,2 90,56 103,53 116,31 129,04
           
           
С воронкой,          
вне бутылки          
№ серии 350 400 450 500 550
1 3892 3631 3595 3467 3384
2 3791 3717 3540 3434 3360
3 3736 3557 3586 3472 3366
4 3856 3583 3457 3469 3408
5 3743 3657 3552 3431 3500
           
среднее 3803,6 3629 3546 3454,6 3403,6
среднее на
столбик
316,97 302,42 295,5 287,88 283,63
           
  350 400 450 500 550
Отношение
концентр.
0,2372 0,2994 0,3504 0,4040 0,4550

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 600 650 700 750 800
1 1242 1359 1462 1530 1627
2 1267 1374 1502 1568 1599
3 1193 1358 1398 1461 1512
4 1276 1422 1515 1552 1602
5 1264 1375 1451 1533 1598
           
среднее 1248,4 1377,6 1465,6 1528,8 1587,6
среднее на
столбик
138,71 153,07 162,84 169,87 176,4
           
           
С воронкой,          
вне бутылки          
№ серии 600 650 700 750 800
1 3291 3283 3241 3216 3161
2 3361 3244 3166 3151 3189
3 3384 3174 3259 3209 3060
4 3371 3355 3220 3187 3100
5 3347 3171 3206 3099 3151
           
среднее 3350,8 3245,4 3218,4 3172,4 3132,2
среднее на
столбик
279,23 270,45 268,2 264,37 261,02
           
  600 650 700 750 800
Отношение
концентр.
0,4968 0,5660 0,6072 0,6425 0,6758

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 850 900 950 1000 1050
1 1633 1719 1776 1784 1881
2 1618 1813 1819 1893 1876
3 1661 1713 1823 1907 1883
4 1716 1757 1800 1770 1910
5 1721 1763 1803 1769 1887
           
среднее 1669,8 1753 1804,2 1824,6 1887,4
среднее на
столбик
185,53 194,78 200,47 202,73 209,71
           
           
С воронкой,          
вне бутылки          
№ серии 850 900 950 1000 1050
1 3089 3076 3100 3080 2944
2 3141 3048 2968 2973 2924
3 3117 3111 3009 2993 2969
4 3130 3051 3088 3033 2951
5 3119 3048 2978 2988 2958
           
среднее 3119,2 3066,8 3028,6 3013,4 2949,2
среднее на
столбик
259,93 255,57 252,38 251,12 245,77
           
  850 900 950 1000 1050
Отношение
концентр.
0,7138 0,7621 0,7943 0,8073 0,8533

 

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 1100 1150 1200 1250 1300
1 1926 1951 2028 2012 2052
2 1963 1931 2009 2081 2065
3 1945 1937 2056 1994 2054
4 1936 1991 2019 2075 2042
5 1957 1938 1980 2006 2064
           
среднее 1945,4 1949,6 2018,4 2033,6 2055,4
среднее на
столбик
216,16 216,62 224,27 225,96 228,38
           
           
С воронкой,          
вне бутылки          
№ серии 1100 1150 1200 1250 1300
1 2959 2940 2899 2896 2858
2 2942 2934 2868 2878 2917
3 2922 2962 2891 2872 2852
4 2955 2964 2921 2926 2859
5 2989 2952 2926 2927 2943
           
среднее 2953,4 2950,4 2901 2899,8 2885,8
среднее на
столбик
246,12 245,87 241,75 241,65 240,48
           
  1100 1150 1200 1250 1300
Отношение
концентр.
0,8783 0,8811 0,9277 0,9351 0,9497

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 1350 1400 1450 1500 1600
1 2055 2136 2207 2138 2166
2 2051 2106 2140 2263 2235
3 2092 2097 2097 2117 2214
4 2149 2110 2136 2158 2172
5 2051 2081 2180 2169 2214
           
среднее 2079,6 2106 2152 2169 2200,2
среднее на
столбик
231,07 234 239,11 241 244,47
           
           
С воронкой          
вне, бутылки          
№ серии 1350 1400 1450 1500 1600
1 2802 2787 2789 2855 2855
2 2956 2842 2926 2828 2804
3 2885 2910 2825 2870 2823
4 2826 2828 2809 2867 2759
5 2892 2842 2896 2839 2754
           
среднее 2872,2 2841,8 2849 2851,8 2799
среднее на
столбик
239,35 236,82 237,42 237,65 233,25
           
  1350 1400 1450 1500 1600
Отношение
концентр.
0,9654 0,9881 1,0071 1,0141 1,0481

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 1700 1800 1900 2000 2100
1 2235 2236 2250 2219 2322
2 2236 2195 2164 2280 2278
3 2215 2264 2211 2300 2286
4 2178 2311 2245 2367 2248
5 2173 2237 2251 2196 2312
           
среднее 2207,4 2248,6 2224,2 2272,4 2289,2
среднее на
столбик
245,27 249,84 247,13 252,49 254,36
           
           
С воронкой,          
вне бутылки          
№ серии 1700 1800 1900 2000 2100
1 2777 2808 2815 2856 2765
2 2733 2803 2905 2742 2758
3 2831 2763 2825 2710 2757
4 2915 2759 2784 2722 2845
5 2809 2784 2753 2764 2749
           
среднее 2813 2783,4 2816,4 2758,8 2774,8
среднее на
столбик
234,42 231,95 234,7 229,9 231,23
           
  1700 1800 1900 2000 2100
Отношение
концентр.
1,0463 1,0771 1,0530 1,0983 1,1000

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 2200 2300 2400 2500 2600
1 2242 2328 2347 2323 2409
2 2299 2296 2300 2342 2361
3 2346 2353 2335 2295 2323
4 2272 2291 2394 2297 2402
5 2260 2261 2381 2319 2294
           
среднее 2283,8 2305,8 2351,4 2315,2 2357,8
среднее на
столбик
253,76 256,2 261,27 257,24 261,98
           
           
С воронкой,          
вне бутылки          
№ серии 2200 2300 2400 2500 2600
1 2786 2818 2800 2750 2648
2 2763 2757 2748 2788 2698
3 2702 2733 2794 2746 2856
4 2766 2804 2747 2799 2679
5 2766 2812 2705 2744 2776
           
среднее 2756,6 2784,8 2758,8 2765,4 2731,4
среднее на
столбик
229,72 232,07 229,9 230,45 227,62
           
  2200 2300 2400 2500 2600
Отношение
концентр.
1,1046 1,1040 1,1364 1,1163 1,1510

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 2700 2800 2900 3000 3200
1 2362 2351 2227 2349 2419
2 2280 2226 2390 2312 2303
3 2313 2370 2353 2303 2280
4 2220 2339 2358 2413 2313
5 2288 2362 2411 2359 2294
           
среднее 2292,6 2329,6 2347,8 2347,2 2321,8
среднее на
столбик
254,73 258,84 260,87 260,8 257,98
           
           
С воронкой,          
вне бутылки          
№ серии 2700 2800 2900 3000 3200
1 2808 2727 2738 2747 2717
2 2731 2816 2725 2697 2737
3 2674 2650 2699 2698 2810
4 2767 2779 2843 2741 2811
5 2837 2749 2684 2786 2764
           
среднее 2763,4 2744,2 2737,8 2733,8 2767,8
среднее на
столбик
230,28 228,68 228,15 227,82 230,65
           
  2700 2800 2900 3000 3200
Отношение
концентр.
1,1062 1,1319 1,1434 1,1448 1,1185

 

74-100
(14)
         
С воронкой,          
в бутылке          
№ серии 3400 3600 3800 4000 4500
1 2328 2284 2359 2383 2340
2 2295 2312 2378 2339 2384
3 2332 2316 2339 2347 2307
4 2351 2226 2327 2352 2338
5 2336 2282 2323 2344 2289
           
среднее 2328,4 2284 2345,2 2353 2331,6
среднее на
столбик
258,71 253,78 260,58 261,44 259,07
           
           
С воронкой,          
вне бутылки          
№ серии 3400 3600 3800 4000 4500
1 2780 2809 2713 2665 2716
2 2812 2776 2733 2713 2720
3 2742 2689 2710 2727 2707
4 2681 2760 2643 2663 2753
5 2671 2739 2710 2708 2692
           
среднее 2737,2 2754,6 2701,8 2695,2 2717,6
среднее на
столбик
228,1 229,55 225,15 224,6 226,47
           
  3400 3600 3800 4000 4500
Отношение
концентр.
1,1342 1,1055 1,1574 1,1640 1,1440

 

74-100
(14)
     
С воронкой,      
в бутылке      
№ серии 5000 6000 10000
1 2354 2377 2372
2 2418 2359 2350
3 2371 2380 2382
4 2379 2328 2297
5 2385 2348 2396
       
среднее 2381,4 2358,4 2359,4
среднее на
столбик
264,6 262,04 262,16
       
       
С воронкой,      
вне бутылки      
№ серии 5000 6000 10000
1 2736 2748 2773
2 2719 2669 2754
3 2703 2693 2775
4 2710 2684 2732
5 2766 2765 2725
       
среднее 2726,8 2711,8 2751,8
среднее на
столбик
227,23 225,98 229,32
       
  5000 6000 10000
Отношение
концентр.
1,1644 1,1596 1,1432

 

 

Результаты расчётов для увеличенной модели 14х20 с зонтиком:

 

74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 20 50 75 100 125
1 0 7 21 36 65
2 0 3 13 43 87
3 0 3 22 47 88
4 0 3 20 53 82
5 0 5 10 46 79
           
среднее 0,0000 4,2 17,2 45 80,2
среднее на
столбик
0 0,47 1,91 5 8,91
           
           
С воронкой,          
вне бутылки          
№ серии 20 50 75 100 125
1 8148 6562 5931 5398 5133
2 8070 6698 5882 5490 5174
3 8063 6626 5846 5450 5070
4 8064 6705 5936 5393 5047
5 8114 6732 5894 5445 5102
           
среднее 8091,8 6664,6 5897,8 5435,2 5105,2
среднее на
столбик
674,32 555,38 491,48 452,93 425,43
           
  20 50 75 100 125
Отношение
концентр.
0,0000 0,0008 0,0039 0,0110 0,0209

 

 

74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 150 175 200 250 300
1 135 183 246 414 574
2 137 179 301 430 571
3 124 203 272 424 577
4 128 185 274 404 562
5 122 215 271 429 568
           
среднее 129,2 193 272,8 420,2 570,4
среднее на
столбик
14,36 21,44 30,31 46,69 63,38
           
           
С воронкой,          
вне бутылки          
№ серии 150 175 200 250 300
1 4729 4570 4409 4102 3889
2 4806 4551 4382 4168 3948
3 4767 4619 4407 4136 3860
4 4786 4625 4477 4180 3898
5 4804 4514 4448 4156 4051
           
среднее 4778,4 4575,8 4424,6 4148,4 3929,2
среднее на
столбик
398,2 381,32 368,72 345,7 327,43
           
  150 175 200 250 300
Отношение
концентр.
0,0361 0,0562 0,0822 0,1351 0,1936
           
           
74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 350 400 450 500 550
1 748 902 1010 1114 1277
2 762 902 999 1157 1324
3 735 890 1061 1155 1314
4 680 924 1057 1142 1307
5 773 926 992 1155 1295
           
среднее 739,6 908,8 1023,8 1144,6 1303,4
среднее на
столбик
82,18 100,98 113,76 127,18 144,82
           
           
С воронкой,          
вне бутылки          
№ серии 350 400 450 500 550
1 3731 3682 3529 3447 3295
2 3736 3650 3550 3450 3357
3 3787 3676 3520 3421 3394
4 3746 3668 3487 3562 3369
5 3761 3624 3589 3453 3336
           
среднее 3752,2 3660 3535 3466,6 3350,2
среднее на
столбик
312,68 305 294,58 288,88 279,18
           
  350 400 450 500 550
Отношение
концентр.
0,2628 0,3311 0,3862 0,4402 0,5187
           
74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 600 650 700 750 800
1 1461 1523 1621 1753 1740
2 1365 1534 1587 1762 1725
3 1424 1511 1542 1739 1800
4 1472 1447 1576 1735 1779
5 1392 1506 1522 1699 1821
           
среднее 1422,8 1504,2 1569,6 1737,6 1773
среднее на
столбик
158,09 167,13 174,4 193,07 197
           
           
С воронкой,          
вне бутылки          
№ серии 600 650 700 750 800
1 3255 3278 3271 3109 3085
2 3356 3260 3232 3119 3113
3 3293 3237 3199 3159 3123
4 3276 3266 3273 3114 3087
5 3307 3213 3176 3146 3126
           
среднее 3297,4 3250,8 3230,2 3129,4 3106,8
среднее на
столбик
274,78 270,9 269,18 260,78 258,9
           
  600 650 700 750 800
Отношение
концентр.
0,5753 0,6170 0,6479 0,7403 0,7609

 

74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 850 900 950 1000 1050
1 1857 1892 1974 2052 2078
2 1877 1936 2018 2067 2070
3 1924 1846 2060 2081 2067
4 1866 1989 2023 2030 2162
5 1872 1937 2002 2029 2056
           
среднее 1879,2 1920 2015,4 2051,8 2086,6
среднее на
столбик
208,8 213,33 223,93 227,98 231,84
           
           
С воронкой,          
вне бутылки          
№ серии 850 900 950 1000 1050
1 3064 3029 3026 2950 3038
2 3049 3071 3005 2996 2894
3 3044 3027 2993 2935 2853
4 3123 2957 3025 2947 2941
5 3098 3036 2944 2918 2806
           
среднее 3075,6 3024 2998,6 2949,2 2906,4
среднее на
столбик
256,3 252 249,88 245,77 242,2
           
  850 900 950 1000 1050
Отношение
концентр.
0,8147 0,8466 0,8962 0,9276 0,9572
           
           

 

74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 1100 1150 1200 1250 1300
1 2130 2223 2133 2216 2317
2 2146 2125 2258 2219 2304
3 2219 2242 2245 2269 2306
4 2183 2222 2201 2237 2344
5 2143 2123 2264 2223 2318
           
среднее 2164,2 2187 2220,2 2232,8 2317,8
среднее на
столбик
240,47 243 246,69 248,09 257,53
           
           
С воронкой,          
вне бутылки          
№ серии 1100 1150 1200 1250 1300
1 2908 2847 2876 2858 2880
2 2922 2958 2847 2879 2795
3 2889 2952 2833 2757 2787
4 2839 2825 2846 2852 2759
5 2841 2888 2810 2873 2836
           
среднее 2879,8 2894 2842,4 2843,8 2811,4
среднее на
столбик
239,98 241,17 236,87 236,98 234,28
           
  1100 1150 1200 1250 1300
Отношение
концентр.
1,0020 1,0076 1,0415 1,0469 1,0992
           
74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 1350 1400 1450 1500 1600
1 2314 2314 2456 2467 2387
2 2296 2403 2372 2412 2497
3 2338 2352 2367 2468 2463
4 2289 2365 2482 2464 2413
5 2366 2380 2412 2388 2487
           
среднее 2320,6 2362,8 2417,8 2439,8 2449,4
среднее на
столбик
257,84 262,53 268,64 271,09 272,16
           
           
С воронкой,          
вне бутылки          
№ серии 1350 1400 1450 1500 1600
1 2748 2840 2746 2671 2743
2 2768 2847 2793 2752 2747
3 2871 2738 2802 2695 2713
4 2758 2817 2683 2760 2720
5 2817 2847 2831 2759 2707
           
среднее 2792,4 2817,8 2771 2727,4 2726
среднее на
столбик
232,7 234,82 230,92 227,28 227,17
           
  1350 1400 1450 1500 1600
Отношение
концентр.
1,1081 1,1180 1,1634 1,1927 1,1980
           
74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 1700 1800 1900 2000 2100
1 2489 2479 2529 2594 2594
2 2502 2510 2614 2528 2547
3 2502 2563 2504 2549 2617
4 2534 2489 2538 2569 2595
5 2515 2496 2517 2514 2542
           
среднее 2508,4 2507,4 2540,4 2550,8 2579
среднее на
столбик
278,71 278,6 282,27 283,42 286,56
           
           
С воронкой,          
вне бутылки          
№ серии 1700 1800 1900 2000 2100
1 2774 2722 2661 2667 2640
2 2707 2679 2660 2638 2684
3 2685 2733 2713 2699 2648
4 2680 2684 2639 2677 2628
5 2720 2690 2718 2609 2669
           
среднее 2713,2 2701,6 2678,2 2658 2653,8
среднее на
столбик
226,1 225,13 223,18 221,5 221,15
           
  1700 1800 1900 2000 2100
Отношение
концентр.
1,2327 1,2375 1,2647 1,2796 1,2958
           
74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 2200 2300 2400 2500 2600
1 2636 2647 2616 2666 2675
2 2610 2630 2577 2694 2731
3 2601 2717 2587 2704 2614
4 2530 2638 2669 2644 2662
5 2583 2654 2620 2589 2719
           
среднее 2592 2657,2 2613,8 2659,4 2680,2
среднее на
столбик
288 295,24 290,42 295,49 297,8
           
           
С воронкой,          
вне бутылки          
№ серии 2200 2300 2400 2500 2600
1 2630 2653 2599 2613 2639
2 2633 2614 2716 2618 2653
3 2673 2615 2628 2655 2672
4 2724 2634 2669 2638 2612
5 2670 2557 2620 2664 2604
           
среднее 2666 2614,6 2646,4 2637,6 2636
среднее на
столбик
222,17 217,88 220,53 219,8 219,67
           
  2200 2300 2400 2500 2600
Отношение
концентр.
1,2963 1,3551 1,3169 1,3444 1,3557

 

74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 2700 2800 2900 3000 3200
1 2675 2622 2582 2607 2596
2 2686 2618 2598 2688 2740
3 2612 2627 2634 2682 2677
4 2552 2588 2646 2672 2682
5 2601 2696 2630 2676 2715
           
среднее 2625,2 2630,2 2618 2665 2682
среднее на
столбик
291,69 292,24 290,89 296,11 298
           
           
С воронкой,          
вне бутылки          
№ серии 2700 2800 2900 3000 3200
1 2629 2616 2584 2612 2649
2 2615 2594 2616 2589 2604
3 2576 2651 2655 2579 2638
4 2666 2692 2644 2586 2625
5 2606 2635 2632 2624 2582
           
среднее 2618,4 2637,6 2626,2 2598 2619,6
среднее на
столбик
218,2 219,8 218,85 216,5 218,3
           
  2700 2800 2900 3000 3200
Отношение
концентр.
1,3368 1,3296 1,3292 1,3677 1,3651

 

74-100
(14 + зонт)
         
С воронкой,          
в бутылке          
№ серии 3400 3600 3800 4000 4500
1 2608 2671 2726 2679 2663
2 2615 2543 2706 2680 2635
3 2632 2699 2649 2656 2661
4 2642 2607 2653 2717 2597
5 2680 2641 2640 2574 2680
           
среднее 2635,4 2632,2 2674,8 2661,2 2647,2
среднее на
столбик
292,82 292,47 297,2 295,69 294,13
           
           
С воронкой,          
вне бутылки          
№ серии 3400 3600 3800 4000 4500
1 2673 2607 2596 2599 2681
2 2645 2636 2583 2571 2600
3 2660 2638 2588 2620 2572
4 2635 2554 2601 2672 2677
5 2533 2564 2650 2659 2590
           
среднее 2629,2 2599,8 2603,6 2624,2 2624
среднее на
столбик
219,1 216,65 216,97 218,68 218,67
           
  3400 3600 3800 4000 4500
Отношение
концентр.
1,3365 1,3499 1,3698 1,3521 1,3451

 

74-100
(14 + зонт)
       
С воронкой,        
в бутылке        
№ серии 5000 6000 8000 10000
1 2620 2761 2693 2648
2 2683 2764 2651 2665
3 2617 2587 2695 2697
4 2640 2657 2760 2667
5 2629 2712 2680 2699
         
среднее 2637,8 2696,2 2695,8 2675,2
среднее на
столбик
293,09 299,58 299,53 297,24
         
         
С воронкой,        
вне бутылки        
№ серии 5000 6000 8000 10000
1 2622 2615 2621 2590
2 2559 2576 2646 2665
3 2676 2614 2570 2545
4 2593 2549 2563 2569
5 2615 2644 2633 2566
         
среднее 2613 2599,6 2606,6 2587
среднее на
столбик
217,75 216,63 217,22 215,58
         
  5000 6000 8000 10000
Отношение
концентр.
1,3460 1,3829 1,3790 1,3788

 

 

 

 

 

 

 

 

Ссылка на комментарий
В 31.12.2025 в 04:53, dmr сказал:

Вы прямо на полном серьёзе, искренне считаете, что вот это всё будет кто-то читать? 

Или вы чисто для себя пишите? 

„Все, что мы теперь достоверно знаем из физики, было прежде облечено в догадки, и если б никогда не допускались догадки, даже ошибочные, то мы бы не добыли ни одной истины.“ (Л. Эйлер)

 

200 лет назад Закон всемирного тяготения Ньютона воспринимался практически всеми, как фундаментальный закон природы многократно проверенный экспериментально, как наблюдением за движением небесных тел, так и в лабораторных условиях (опыты Кавендиша 1797 - 1798 годы).

В 1859 году  выяснилось, что с помощью Закона всемирного тяготения нельзя точно рассчитать смещение перигелия Меркурия.

Более 50 лет данную аномалию безуспешно пытались объяснить в рамках Закона всемирного тяготения, а потом появилась теория гравитации Эйнштейна, где сила гравитации - фиктивная сила, и формула всемирного тяготения — просто приближение, которое хорошо работает в окрестности Земли. В общем случае закон всемирного тяготения неправильно описывает действительность.    

Так что один прецедент, когда теорию считающуюся общепризнанным фундаментальным законом природы опровергали уже был.

Изменено пользователем Мансур Фазлуллин
Увидел опечатку
Ссылка на комментарий
В 31.12.2025 в 12:58, Мансур Фазлуллин сказал:

истины

Вы вообще прочитали то что процитировали мои слова? 

 

В 31.12.2025 в 06:53, dmr сказал:

Вы прямо на полном серьёзе, искренне считаете, что вот это всё будет кто-то читать? 

Или вы чисто для себя пишите? 

Такое ощущение что вы просто нажали кнопку цитата дальше развалились в кресле в позе Мавроди руки за голову и глядя в потолок начали рассуждать о своём 

Я вас всего лишь спросил Вы действительно считаете что все ваши длинные посты кто-то читает? 

Ссылка на комментарий
В 31.12.2025 в 12:43, dmr сказал:

Вы вообще прочитали то что процитировали мои слова? 

 

Такое ощущение что вы просто нажали кнопку цитата дальше развалились в кресле в позе Мавроди руки за голову и глядя в потолок начали рассуждать о своём 

Я вас всего лишь спросил Вы действительно считаете что все ваши длинные посты кто-то читает? 

Самая важная информация в моём посте от 29 декабря это то, что идёт в начале - 2 картинки и 13 строк текста.

Далее идёт код программы и результаты расчётов. - Эту информацию я предоставил на случай если кто-то засомневается в полученных мною результатах и захочет их проверить. Может я эти цифры получил не из расчётов, а придумал чтобы подтвердить свою точку зрения. Или может в предложенном мною коде есть какие-то ошибки, которые приводят к росту концентрации частиц в бутылке. 

Ссылка на комментарий
В 29.12.2025 в 20:33, Мансур Фазлуллин сказал:

Как изменятся результаты при моделировании движения более мелких частиц, которые движутся вниз при z от 75 до 100?

Результаты расчётов для более мелких частиц представлены в виде графика ниже.

(75).thumb.png.0366724dd9149059839c4fc49cd82bc6.png

Из графика видно, что при увеличении модели и размеров воронки (ярко-оранжевый график). Скорость изменения отношений концентраций замедлилась, по сравнению с предыдущей, меньшей моделью (синий график), при этом максимальные значения отношений концентраций не изменились. 

При добавлении в модель "зонтика" (жёлтый график) скорость изменения концентраций и их максимальные значения также увеличились, как и для более крупных частиц.

При этом, если сравнивать график, для более крупных частиц (график в посте от 29 декабря) и график в данном посте. Видно, что с уменьшением размеров частиц, скорость изменения отношений концентраций и их максимальные значения увеличились.

Для более крупных частиц я принял, что они движутся вверх при z от 51 до 73 и движутся вниз при z от 74 до 100. Т.е. для крупных частиц вероятность движения вниз выше на 17%, чем вверх.

Для мелких частиц я принял, что они движутся вверх при z от 51 до 74 и движутся вниз при z от 75 до 100. Т.е. для мелких частиц вероятность движения вниз выше на 8%, чем вверх.

Данные допущения я взял "с потолка" и вполне возможно, что для реальных броуновских частиц, разница в вероятности движения вверх или вниз значительно меньше.

Вопрос о разнице в вероятности движения броуновских частиц вверх или вниз я решил прояснить с помощью ИИ.

Для этого я ввёл запрос:

На сколько уменьшается концентрация броуновских частиц при подъёме на 1 мм, 1 см, 1 м, 10 м, 50 м.

На который я получил следующий ответ:

Уменьшение концентрации броуновских частиц с высотой в гравитационном поле можно оценить с помощью некоторых упрощенных моделей, однако точные значения будут зависеть от многих факторов, таких как размер и масса частиц, вязкость среды и температура. Для решения данной задачи можно использовать концепцию термодинамического равновесия и формулу барометрической функции.

Общие соображения

1. Сила тяжести: С увеличением высоты под действием силы тяжести концентрация частиц будет уменьшаться.

2. Диффузия: В нейтральной среде наличие диффузии может также влиять на равновесие частиц.

Эмпирический подход

Для оценки, можно использовать приблизительную зависимость, основанную на экспоненциальном затухании концентрации с высотой:

C(h) = C0∙e-h/H

где:

C(h) – концентрация на высоте h,

C0 начальная концентрация на уровне земли,

 Н – характерная высота (может варьироваться в зависимости от среды, например, для воздуха это может быть порядка 8-10 км).

Пример расчета

Допустим, H составляет 8 км, тогда можно оценить изменение концентрации при разных высотах:

1 мм (0,001 м): C(0,001) C0∙e-0,001/8000 ≈ C0 изменение незначительное.

1 см (0,01 м): C(0,01) C0∙e-0,01/8000 ≈ C0 изменение всё еще незначительное.

1 м: C(1) C0∙e-1/8000 ≈ C0 ∙0,999875 примерно 0,0125% снижение.

10 м: C(10) C0∙e-10/8000 ≈ C0 ∙0,99875 примерно 0,125% снижение.

50 м: C(50) C0∙e-50/8000 ≈ C0 ∙0,99375 примерно 0,625% снижение.

Подводя итоги

Как видно, при малых высотах (несколько метров) концентрация броуновских частиц практически не изменяется. Значительное снижение начнет проявляться только на высотах порядка десятков метров и больше, но для малых высот (до 10 м) оно остается незначительным.

Эти результаты основаны на упрощенных моделях и могут варьироваться в зависимости от реальных условий и характеристик среды.

 

Из ответа ИИ понятно, что влияние гравитации на движение броуновских частиц достаточно мало. Так как рассматриваемая мною модель имеет микроскопические размеры, то концентрации броуновских частиц в нижней и в верхней части будут равны, а значит вполне допустимо принять, что вероятности движения частиц вверх и вниз равны.  

Т.е. скорости увеличения разницы отношений концентраций и их максимальные значения, скорее всего ещё вырастут.

Тут возникает более общий вопрос. 

Если в моих рассуждениях есть какое-то здравое зерно. И при правильном выборе материала бутылки с воронкой и броуновских частиц можно изготовить микроскопические бутылку с воронкой и добиться повышенной концентрации броуновских частиц в бутылке, по сравнению с внешней средой, как из этого извлечь практическую пользу?  

Если воронка позволяет повысить концентрацию броуновских частиц в бутылке, то вполне вероятно, что если развернуть воронку на 180 градусов и соединить с бутылкой широким основанием, то она начнёт работать наоборот и начнёт уменьшать концентрацию броуновских частиц в бутылке, по сравнению с окружающей средой. По крайней мере данное предположение можно проверить моделированием.

Тогда можно собрать следующую конструкцию:

.thumb.PNG.fbd5f077232fc9024099516110bb73d6.PNG

 

На карусели в положении на 3 часа закреплена "бутылка" в верху бутылки воронка, днище бутылки выполнено в виде конуса (воронки развёрнутой на 180 градусов, если перевёрнутую воронку называть днищем, то может возникнуть путаница, т.к потом это днище превращается в крышку, поэтому решил называть её конусом). Также в бутылке есть две заслонки закреплённые на шарнирах, на рисунке заслонки показаны толстыми линиями. Заслонки сопоставимы по размерам с броуновскими частицами, но крупнее их и не способны совершать броуновские движения.

Если посмотреть на бутылку, которая расположена на 3 часа, то заслонка воронки, под воздействием силы тяжести висит вертикально и не мешает броуновским частицам проходить через верхнюю воронку. Заслонка для конуса лежит на нём горизонтально и закрывает отверстие для броуновских частиц в нём.

Если посмотреть на бутылку на 9 часов. То заслонка для конуса висит вертикально и открывает в нём отверстие. Заслонка для воронки лежит на дне и закрывает отверстие воронки.

Допустим в начальный момент времени карусель находится в некоторой замкнутой ёмкости в которой отсутствуют броуновские частицы. Карусель сбалансирована и находится в состоянии покоя.

Затем в замкнутую ёмкость извне вносят броуновские частицы.

Броуновские частицы начинают заполнять ёмкость и попадать в бутылки. Попадать в бутылку на 3 часа им помогает воронка. С ростом концентрации броуновских частиц в  воронке на 3 часа её масса начинает расти.

В бутылку на 9 часов частицы то же попадают, но конус способствует их выходу из бутылки. В итоге концентрация броуновских частиц в бутылке на 9 часов в начальный момент тоже растёт, но медленнее чем в бутылке на 3 часа.

Процесс описанный для бутылки на три часа происходит во всех правых бутылках.

Процесс описанный для бутылки на девять часов происходит во всех левых бутылках. 

В итоге бутылки справа перевешивают левые бутылки и карусель начинается вращаться по часовой стрелке.

Проходя через  нижнюю точку, заслонка воронки закрывает в ней отверстие, заслонка конуса, наоборот, открывает в нём отверстие. Количество частиц в бутылках после прохождения нижнего положения начинает уменьшаться. 

После прохождения верхней точки всё происходит наоборот. Количество частиц в бутылках после верхнего положения начинает расти.

 

Я думаю подобные "вечные крутилки" все помнят из школьной программы. 

Вопрос где ошибки в моих рассуждениях?

 

 

 

 

 

 

 

 

 

Ссылка на комментарий
В 03.01.2026 в 22:41, Мансур Фазлуллин сказал:

Из ответа ИИ понятно, что влияние гравитации на движение броуновских частиц достаточно мало. Так как рассматриваемая мною модель имеет микроскопические размеры, то концентрации броуновских частиц в нижней и в верхней части будут равны, а значит вполне допустимо принять, что вероятности движения частиц вверх и вниз равны.

Результаты расчётов при условии, что вероятности движения вправо, влево, вверх и вниз одинаковы. Представлены на графике. 

 

(1-4).thumb.png.b0be1605703d7b2bc92d572ba0fbded9.png

Какое будет отношение концентраций внутри и вне бутылки если воронку развернуть на 180 градусов, как показано на модели карусели в левой части?

 

Ниже программа для модели 12 х 16 в которой вероятности движения вправо, влево, вверх и вниз равны. В программы для модели 14 х 20 и для модели 14 х 20 + зонт внесены аналогичные изменения.

 

import random
import matplotlib.pyplot as plt

n = 0

n_out = 0

x_n = list(range(17))

y_n = []

for i in x_n:
    y_n += [0]


x_v = list(range(17))

y_v = []

for i in x_v:
    y_v += [12]

    
y_l = list(range(13))

x_l = []

for i in y_l:
    x_l += [0]
    
y_p = list(range(13))

x_p = []

for i in y_p:
    x_p += [16]


x_b = [5, 5, 5, 5, 6, 7, 8, 9, 10, 11, 11, 11, 11]

y_b = [0, 1, 2, 3, 3, 3, 3, 3, 3, 3, 2, 1, 0] 

x_b1 = [4, 4, 4, 4, 4, 5, 6, 7, 8, 9, 10, 11, 12, 12, 12, 12, 12]

y_b1 = [0, 1, 2, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 2, 1, 0]

fnx_l = [7, 6, 5, 4, 5, 6, 7, 8]

fny_l = [4, 5, 6, 7, 7, 6, 5, 4]

fnx_r = [8, 9, 10, 11, 12, 11, 10, 9]

fny_r = [4, 5, 6, 7, 7, 6, 5, 4]

P = 0
k = True

while k:
    P += 1
    t = 0
    x = 1
    y = 1
    xn = x
    yn = y
    x1 = []
    y1 = []


    while t < 10000:
        t += 1

        z = random.randint(1, 4)

                
##        # Условия для движения вправо при z == 1

        if z == 1:
            if x < 16:
                x += 1

            if (x == 5 and y <= 3) or (x == 12 and y <= 3):
                x -= 1

            if x == 8 and y == 4:
                x -= 1

            if x == 7 and y == 5:
                y -= 1

            if x == 6 and y == 6:
                y -= 1

            if x == 9 and y == 4:
                y += 1

            if x == 10 and y == 5:
                y += 1

            if x == 11 and y == 6:
                y += 1


        # Условия для движения влево при z == 2


        elif z == 2:
            if x > 0:
                x -= 1

            if (x == 4 and y <= 3) or (x == 11 and y <= 3):
                x += 1

            if x == 7 and y == 4:
                y += 1

            if x == 6 and y == 5:
                y += 1

            if x == 5 and y == 6:
                y += 1

            if x == 8 and y == 4:
                x += 1

            if x == 9 and y == 5:
                y -= 1

            if x == 10 and y == 6:
                y -= 1

        # Условия для движения вверх z == 3

        elif z == 3:
            if y < 12:
                y += 1

            if (y == 4 and 5 <= x < 8) or (y == 4 and 8 < x <= 11):
                y -= 1

            if x == 7 and y == 5:
                x -= 1

            if x == 6 and y == 6:
                x -= 1

            if x == 5 and y == 7:
                x -= 1

            if x == 9 and y == 5:
                x += 1

            if x == 10 and y == 6:
                x += 1

            if x == 11 and y == 7:
                x += 1


        # Условия для предотвращения проскока


        elif x == 7 and y == 4 and z == 4:
            z = 0

        elif x == 9 and y == 4 and z == 4:
            z = 0        


        # Условия для движения вниз при z == 4


        elif z == 4:
            if y > 0:
                y -= 1

            if (y == 3 and 5 <= x < 8) or (y == 3 and 8 < x <= 11):
                y += 1

            if x == 7 and y == 4:
                x += 1

            if x == 6 and y == 5:
                x += 1

            if x == 5 and y == 6:
                x += 1

            if x == 9 and y == 4:
                x -= 1

            if x == 10 and y == 5:
                x -= 1

            if x == 11 and y == 6:
                x -= 1


        x1.append(x)
        y1.append(y)


    if 5 <= x1[-1] <= 11 and y1[-1] <= 3:
        n += 1

    if  x1[-1] <= 4 and y1[-1] <=3:
        n_out += 1

    if 12 <= x1[-1] and y1[-1] <=3:
        n_out += 1

        

    if P == 10000:
        k = False


print('n=', n)

print('n/P', n/P)

print('n_out=', n_out)

print('n_out/P', n_out/P)


plt.grid(True)
plt.plot(x_n, y_n, 'k')
plt.plot(x_v, y_v, 'k')
plt.plot(x_l, y_l, 'k')
plt.plot(x_p, y_p, 'k')
plt.plot(x_b, y_b, 'k')
plt.plot(x_b1, y_b1, 'k')
plt.plot(fnx_l, fny_l, 'k')
plt.plot(fnx_r, fny_r, 'k')
plt.plot(x1, y1, 'b')
plt.plot(x1[-1], y1[-1], 'ro')
plt.plot(xn, yn, 'go')
plt.show()    

 

Изменено пользователем Мансур Фазлуллин
В программе "8)" превратилась в смайлик
Ссылка на комментарий

Для публикации сообщений создайте учётную запись или авторизуйтесь

Вы должны быть пользователем, чтобы оставить комментарий

Создать аккаунт

Зарегистрируйте новый аккаунт в нашем сообществе. Это очень просто!

Регистрация нового пользователя

Войти

Уже есть аккаунт? Войти в систему.

Войти
  • Последние посетители   0 пользователей онлайн

    • Ни одного зарегистрированного пользователя не просматривает данную страницу
×
×
  • Создать...

Важная информация

Мы разместили cookie-файлы на ваше устройство, чтобы помочь сделать этот сайт лучше. Вы можете изменить свои настройки cookie-файлов, или продолжить без изменения настроек.