Математическая модель элиминации
В. А. Грозных
Пусть задано множество элементов – простых
сомножителей (единиц, штук, особей) , которые образуют текст длиной Т
(общий перечень особей 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.