Математическая модель элиминации

В. А. Грозных

 

Пусть задано множество элементов – простых сомножителей (единиц, штук, особей) , которые образуют текст длиной Т (общий перечень особей U), и пусть каждая особь-штука классифицируется так, что может быть отнесена к тому или иному виду. Тогда общее количество видов S составит словарь объёмом V. Особи одной численности объединяются в популяции (виды, представленные одинаковым количеством особей), образуя группы kєK, называемые кастами [1].

Пусть i=1, 2, 3, ... – возможная численность популяции: количество двоек, троек, пятёрок, семёрок и так далее (из ряда, соответствующего натуральному ряду чисел, последовательно подсчитывается общее количество каждого из встретившихся простых чисел при разложении каждого натурального числа на простые числа-сомножители). Видовое распределение может быть получено из текста Т (факториала N!) непосредственно, если выбрать вначале все виды, состоящие из одной особи (это, собственно все простые числа, большие чем N/2); они образуют тем самым первую (ноеву) касту k=1, общее число видов s в которой составляет , эмпирическая численность особей в касте – . Затем – все виды, представленные двумя особями, тремя и так далее (если все значимые строки нумеровать по порядку, то в этом случае число строк равно числу каст К, где К есть наличествующие популяции). Последовательность  назовём эмпирическим видовым распределением (распределением видов). Будем упрощённо считать однозначными обозначения [2]:

,                                                   (1)

где  – непрерывный аналог мощности (численности) популяций i (i – всегда дискретная величина, i=[x]; α>0 – характеристический показатель); постоянная распределения γ=1+α; , , где  – теоретическое и  – фактическое значение первой точки (ноева каста); А – постоянная распределения, которую находят из условий нормировки (хотя это теоретически и ошибочно из-за отсутствия математического ожидания и бесконечности дисперсии) [5].

Обозначим через  самую мощную (саранчёвую) популяцию (касту), т. е. численность вида, представленного наибольшим количеством особей (в данном случае это будут двойки).

Видовые распределения отличаются характером изменения . Устойчивую зависимость показывают "гипербола" Ω(х); для S(U) характерно относительно более медленное увеличение количества видов при увеличении выборки штук-особей.

Повторяемость d=U/S растёт;  – ноева каста (при увеличении выборки эта величина медленно уменьшается, как того требует теорема Гнеденко–Дёблина для бесконечно делимых негауссовых распределений).

Вторая форма Н-распределения – ранговидовое распределение Λ(r). Оно по определению получается из видового:  – количество особей вида (численность популяции  вида), соответствует рангу r при общем числе особей U (длина текста Т=|U|). Ранг вида s=1, 2, ..., , ..., S – это его порядковый номер (номер строки). Последний номер S, равный количеству рангов, определяет объём словаря V, можно записать V=|S|. Функция  записывается в виде:

, , ,                                 (2)

где В – абсолютная величина и характеристический показатель β>0 – константы ранговидового гиперболического Н-распределения.

Опираясь на формулы (1) и (2), предложим "модель простых чисел" гиперболического Н-распределения для изучения структуры физических, биологических, технических, информационных, социальных ценозов по критериям видового разнообразия и соотношения "крупное-мелкое" [3].

Примем в качестве канонического дискретное распределение простых сомножителей в факториале некоторого числа N. Назовём видом любое простое число , где r – номер простого числа натурального ряда чисел, абстрактно воспринимаемое из ряда: 2, 3, ..., , ..., а особью – появление этого простого числа как сомножителя (единица исключается) в любом из чисел натурального ряда. Тогда каждое натуральное число  представимо следующим образом:

, , (j=0, 1, 2, ..., m),                   (3)

где m – степень (встречаемость) простого числа, r – ранг простого числа. Например,  состоит из сомножителей , где вид  – двойка встретился как особь три раза, вид  – один раз. В целом для факториала, например, =100! двойка (саранчёвый вид) встретилась  встретилась (как особь) , тройка – 48 раз (, ), 10 чисел встретилось по одному разу (ноева каста). Общее число видов в системе S равно максимальному значению ранга – 11, суммарное количество особей равно тексту U – 239. Построение модели элиминации особей производилось при различных заданиях числа уничтожаемых особей, при этом полагалось, что происходит не просто исчезновение особи, а поглощение её более развитыми видами. Таким образом, если вид, состоящий из троек, поглощает вид из двоек, то при полном поглощении для 100! справедливо равенство:

,                                                                          (4)

решением которого является:

.                                                              (5)

Моделирование элиминации производилось различными способами:

1) последовательным увеличением всех видов на одну особь за счёт предыдущих – модель неработоспособна, поскольку в данном случае происходит элиминация высших видов за исключением последнего, число особей которого будет равно двум;

2) увеличением всех видов старше тройки на одну единицу – в данном случае график перемещается вправо, а саранчёвые виды исчезают. Найдена закономерность – оказывается, не только виды подчиняются Н-распределению, но и их потребности в энергии, росте, развитии.

Наиболее эволюционно развитый вид (для 100! это простое число 97) обладает характеристиками, способствующими его выживанию в большей мере, чем все остальные виды. Развитие данного вида, таким образом, осуществляется гораздо быстрее, если происходит поглощение эволюционно более развитых видов, поскольку таковые при выбранных условиях (100!) отсутствуют, и быстрое развитие наблюдается при поглощении себе подобных. Возникает противоречие – развитие отдельно взятой особи наиболее развитого вида при отсутствии ресурсов извне и прочих видов достигается за счёт уменьшения численности вида, т. е. элиминация вида в целом способствует развитию отдельно взятой особи. Стоит вспомнить утверждение Кудрина Б. И. – "наивысшая напряжённость конкурентной борьбы наблюдается именно внутри популяции одного вида". Выполненные азложения 100! (не приводятся), 1000! и 10000! (представлены таблично); элиминации саранчёвых видов с последующим ростом остальных видов на одну особь представлены на рис. 1–3.

Рис. 1. Элиминация саранчёвых видов (двоек и троек) с последующим ростом остальных видов

а – на одну особь для 100!

б – на одну особь для 1000!

в – на одну особь для 10000!

 

1. Разложение числа 1000!

Простое

число

Количество раз встречено

Простое

число

Количество раз встречено

Простое

число

Количество раз встречено

2

994

269

3

617

1

3

498

271

3

619

1

5

249

277

3

631

1

7

164

281

3

641

1

11

98

283

3

643

1

13

81

293

3

647

1

17

61

307

3

653

1

19

54

311

3

659

1

23

44

313

3

661

1

29

35

317

3

673

1

31

33

331

3

677

1

37

27

337

2

683

1

41

24

347

2

691

1

43

23

349

2

701

1

47

21

353

2

709

1

53

18

359

2

719

1

59

16

367

2

727

1

61

16

373

2

733

1

67

14

379

2

739

1

71

14

383

2

743

1

73

13

389

2

751

1

79

12

397

2

757

1

83

12

401

2

761

1

89

11

409

2

769

1

97

10

419

2

773

1

101

9

421

2

787

1

103

9

431

2

797

1

107

9

433

2

809

1

109

9

439

2

811

1

113

8

443

2

821

1

127

7

449

2

823

1

131

7

457

2

827

1

137

7

461

2

829

1

139

7

463

2

839

1

149

6

467

2

853

1

151

6

479

2

857

1

157

6

487

2

859

1

163

6

491

2

863

1

167

5

499

2

877

1

173

5

503

1

881

1

179

5

509

1

883

1

181

5

521

1

887

1

191

5

523

1

907

1

193

5

541

1

911

1

197

5

547

1

919

1

199

5

557

1

929

1

211

4

563

1

937

1

223

4

569

1

941

1

227

4

571

1

947

1

229

4

577

1

953

1

233

4

587

1

967

1

239

4

593

1

971

1

241

4

599

1

977

1

251

3

601

1

983

1

257

3

607

1

991

1

263

3

613

1

997

1

Итого простых чисел:

168

 

2. Базовые параметры разложения 1000!

k, ранг

a(i), количество особей, входящих в вид

w(i), число видов, имеющих одинаковое число особей

a(i)w(i), число особей касты

k, ранг

a(i), количество особей, входящих в вид

w(i), число видов, имеющих одинаковое число особей

a(i)w(i), число особей касты

1

1

73

73

17

21

1

21

2

2

28

56

18

23

1

23

3

3

14

42

19

24

1

24

4

4

7

28

20

27

1

27

5

5

8

40

21

33

1

33

6

6

4

24

22

35

1

35

7

7

4

28

23

44

1

44

8

8

1

8

24

54

1

54

9

9

4

36

25

61

1

61

10

10

1

10

26

81

1

81

11

11

1

11

27

98

1

98

12

12

2

24

28

164

1

164

13

13

1

13

29

249

1

249

14

14

2

28

30

498

1

498

15

16

2

32

31

994

1

994

16

18

1

18

 

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

,                                                                    (6)

допустим, произошла полная элиминация вида, представленного тройкой:

.                                                                    (7)

Будем так же полагать, что прямых связей между видами 5 и 2 не существует. В данном случае исчезает питательная база для вида 5, происходит его исчезновение, равно как и вида 7, вид 2 не ограничен другими видами и продолжает своё развитие далее до тех пор, пока не ограничится его ресурсная база или не появится в ходе эволюции вид, способный ограничить его численность.

Природой был найден выход из данной ситуации – не секрет, многие виды потребляют набор из других видов, что предотвращает сильные колебания численностей видов, а также сильные перекосы в экосистеме. Естественно, для большого количества видов смоделировать всёе многообразие межвидовых взаимоотношений невозможно, в общем случае и не нужно. Однако развитие отдельных видов становится вполне предсказуемым. Сравнение смоделированных кривых (изменения численностей видов) с Н-распределением стабильной экосистемы даёт возможность оценить  возникающие по тем или иным причинам патологии.

Введём понятие – высший эквивалент экосистемы, предположим:

,                                                                (8)

Тогда                  ; ;   ,                                                              (9)

.                                                                    (10)

Полученное значение округляем до наименьшего целого. Найденное число упаковывает всю последовательность видов в один наиболее конкурентоспособный и развитый вид и показывает наибольшую теоретически возможную численность особей этого вида. После достижения максимума теоретически возможны два случая: 1) быстрая элиминация; 2) медленная элиминация, которая сопровождается поглощением себе подобных.

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

 

            Литература

1. Кудрин Б. И. Философия технетики: основания постнеклассической философии техники. Вып. 36. "Ценологические исследования". М.: Технетика, 2007.

2. Кудрин Б. И. Классика технических ценозов. Общая и прикладная ценология. Вып. 31. "Ценологические исследования". Томск: Томский государственный университет – Центр системных исследований, 2006.

3. Кудрин Б. И., Ошурков М. Г. Электрика: Объект. Математика. Словарь. – Томск: Изд-во Томск. гос. ун-та, 2004.

4. Электрификация металлургических предприятий Сибири. Вып. 13 / под ред. проф. Б. И. Кудрина. М.: Издательский дом МЭИ, 2007.

5. Кудрин Б. И., Гамазин С. И., Цырук С. А. Электрика: Классика. Вероятность. Ценология. Монографическое издание. Вып. 34. "Ценологические исследования". М.: Технетика, 2007.