Моделирование внутриэтнических эволюционных процессов Текст научной статьи по специальности «Математика»

УДК 517-9

МОДЕЛИРОВАНИЕ ВНУТРИЭТНИЧЕСКИХ ЭВОЛЮЦИОННЫХ ПРОЦЕССОВ MODELING OF INTRAETHNIC EVOLUTIONARY PROCESSES

Р.О. Кенетова R.O. Kenetova

Институт прикладной математики и автоматизации, Россия, 360000, г. Нальчик, ул. Шортанова, д. 89 а Institute of Applied Mathematics and Automation, 89 а, Shortanov St, Nalchik, 360000, Russia

E-mail: raisa.kenetova@mail.ru

Аннотация. В работе впервые с помощью логистических моделей с учетом предыстории описываются внутриэтнические эволюционные процессы на стадии подъема. Чтобы учесть адекватно предысторию или стохастический характер внутриэтнической эволюции приводится математическая модель, основанная на дифференциальном уравнении дробного порядка.

Resume. For the first time using logistic models based on prehistory intraevolutionary processes are described at the stage of recovery. To take into account adequately the background or the stochastic nature in ethnic evolution is a mathematical model based on differential equations of fractional order.

Ключевые слова: синергетика, пассионарность, субпассионарность, внутриэтническая эволюция, теория популяции, логистический закон, закон Мальтуса.

Key words: synergetics, passionarity, subpassionarity, intraethnic evolution, the theory of population, logistics law, Malthusian Law.

Введение

Математическое моделирование на современном этапе представляет собой самостоятельное направление в науке, играет синтезирующую роль, объединяя разные методы и подходы математики, физики, биологии, социологии и вступает в новый этап своего развития, обращаясь к изучению общественных процессов. Изучение развития общества началось с моделирования биосферных процессов Земли. Биосфера получила статус научного предмета благодаря целостному, комплексному подходу к ней, разработанному в трудах В.И. Вернадского [1]. Согласно его теории живое вещество планеты связано с косным веществом посредством биохимической энергии, которую он также называл биогенной миграцией атомов. Эта энергия обеспечивает существование всех живых организмов, запасы которой пополняются за счет солнечной энергии. Модели биосферы строились на основе описания процесса «перетекания» биохимической энергии из одной системы в другую. Жизнь людей также обеспечивается наличием этой энергии. Избыток этой энергии нарушает равновесие в природе, поскольку высвобождение избытка энергии приводит к совершению работы. Вспышки биохимической энергии среди людей называют пассионарными импульсами, т.е. поведенческие импульсы, направленные против инстинкта самосохранения [2].

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

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

Реализуемые модели этнических процессов можно получить, если исходить из математической теории борьбы за существование [3].

Этнические системы и процессы, как правило, являются синергетическими, поэтому и методы нового научного направления, названного Г. Хакеном [4] синергетикой, найдут применение при их моделировании.

Рассмотрим процесс внутриэтнической эволюции. Пусть в этносе из N = N (t) человек, живущих в некоторой географической области с определенной системой связи, в момент времени t содержатся U = U(t) пассионариев и V = V(t) субпассионариев. Допустим, что фаза подъема внутриэтнической эволюции начинается в момент времени t = 0 и заканчивается при t = tm . Очевидно, что области значения R(N), R(U), R(V) функций N,U,V принадлежат множеству целых чисел. Поэтому, вместо разрывных целочисленных функций N,U ,V введем непрерывно дифференцируемые функции n = n(t), u = u(t),v = v(t), имеющие в каждый момент времени те же целые части, что и N,U,V соответственно, т.е. [n(t)] = N(t)„[u(t)] = U(t\[v(t)] = V(t) для всех t е[0, tm ]. Будем предполагать, что

n(t) = u(t) + v(t), 0 < t < tm. (1)

В теории популяции, в том что числе человеческой, наиболее апробированным законом развития является логистический закон роста [5, с. 103]

n\t) = ц n(t) —£ n 2(t), (2)

K

где К - емкость среды обитания, ^ - коэффициент автоприроста. Уравнение (2) можно переписать в виде

d [ 1 ) £ £

л ^ п(/)) к п(/)

Уравнение (3) является линейным относительно функции 1/ п(I) и его единственное решение, удовлетворяющее начальному условию

п(0) = п0, (4)

задается хорошо известной в биологии и медицине формулой

п(/) =-^-. (5)

по + (К - по)ехр(-ц/)

Рис. 1. График функции n(t) при K = 60, n0 = 6, ц = 2, tm = 100.

Fig. 1. The graph of the function n(t) at K = 60,n0 = 6, ц = 2,tm = 100.

График функции (5) (рис. l) можно использовать в качестве приближения к реальному закону роста численности этноса на фазе подъема.

Синергетика пассионериев и субпассионариев носит характер сходный со взаимодействием «хищник-жертва». На фазе подъема в качестве жертвы выступают субпассинарии, а в роли хищника - пассионарии. В соответствии с этим уравнения движения могут записаны в виде

м' (t) = ^ u — y1 uv,

v' (t) = —ц2 v + у 2 uv,

(6)

где 0 < : < :т, м, ц2 - неотрицательные параметры автоприроста пассионариев и субпассионариев, у1, у2 - коэффициенты их взаимодействия.

Из системы (6) следует, что если вектор w = (и, V) - решение (6), то

у 2 [м ' (t) — ^u(t )] = —Yl [v (t) — ^2v(t)].

Формулы (i) и (5) дают основания записать

v(t) = n(t) — u(t),

v'(t) = n ' (t) — u'(t).

(7)

(8)

Подставляя формулы (8) в соотношение (7) получим дифференциальное уравнение роста числен-

ности пассионариев

au' (t) — bu(t) = —у [n(t ) — ц2 n(t)],

(9)

а = У2 -Yl, Ъ = М1У2 +^27Х- (10)

Если суммарная численность п(:) пассионариев и субпассионариев меняется по закону Мальтуса [5, с. 103]

п(:) = по ехр^О, (11)

то из дифференциального уравнения (9) следует, что пассионарии на фазе подъема растут по экспоненциальному закону

u(t) = u(0)exp| bt |.

(12)

Здесь в формуле (12) естественно предположим, что у2 > у, т.е. а > 0.

Когда численность этноса меняется по логистическому закону (5) динамику численности пассионариев в силу уравнения (9) определяется формулой

и(/) = и(0)ехр( —/I—1 [[п'(х)п(х)]ехр

/7 П 3

--(/-х) а

dx.

(13)

Как видно из логистического закона роста (2)

d , и „ ч

—1п п(/) = ц-и п(/),

К

поэтому формула (13) перепишется в виде

и(1) = и(0)ехр| 1п(т)[ц-и! п(х)]ехр

а } а

К

— (I-X) а

—X.

(14)

*

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

и(/) = и(0)ехр| — /1 + !^ а } а

1 ехр

- — (I-х)

а

—х.

(15)

Так как при и < / < /

то из равенства (15) имеем

1 г

-I ехр

а О

— (/-X) а

—X =1 Ь

ехр|--г |-1

и(/) = и(0)ехр( 1 + 11^* I а } —

ехр|--г |-1

(16)

Другими словами, при п(/) = К, численность пассионариев меняется по закону (1), который близок к экспоненциальному.

Логистический закон (2) не учитывает адекватно предысторию или стохастический характер внутриэтнической эволюции. Чтобы учесть эти факторы можно воспользоваться вместо дифференциального уравнения (3) следующей математической моделью

Б0

1

п(п)

п(/)

(17)

где а, и, 8, р - постоянные параметры модели, которые характеризуют тот или иной этнос, причем

1 — г>

0 <а< 1, Р>-1; Б0а,и(х) =

и(х)(/ -X) а - оператор дробного дифференцирования по-

Г(1 -а) ->°

рядка а с началом в начальный момент времени / = 0 [5, с. 28], Г (а) - гамма функция Эйлера.

Относительно функции 2(/) -1/ п(/) уравнение (17) является линейным дифференциальным уравнением дробного порядка а :

Д 2(х) + ^(0 = 8/р, 0 </< 1т. (18)

Решение уравнения дробного порядка (18), удовлетворяющее начальному условию

Нш/1а 2(0 = 0

/^0

имеет вид [6, с. 92]

2(/) = 8Д?хр -ц8Г(а)Б07Е1/2[-ц(/-х)а;аД-У}

0

0

0

™ 2к

где Е (г, м) = У--функция типа Миттаг-Леффлера [6, с. 117].

р к=0 г(м+к / р)

Заключение

В работе впервые с помощью логистических моделей с учетом предыстории описываются внутриэтнические эволюционные процессы на стадии подъема. Приведенные модели могут представлять несомненный интерес хотя бы для ориентировочной оценки численности этноса в периодах подъема.

Список литературы References

1. Вернадский В.И. 1988. Философские мысли натуралиста. М., Наука.

Vernadsky V.I. 1988. Philosophskie mysli naturalista [Philosophical thoughts naturalist]. Moscow, Nauka.

2. Коробицын В.В. 2001. Моделирование этнических полей. Омск, Вестник Омского государственного университета.

Korobitsin V.V. 2001. Modelirovanie etnicheskikh poley [Modelling of ethnic fields]. Omsk. Bulletin of the Omsk State University.

3. Вольтерра В. 1976. Математическая теория борьбы за существование. М., Наука. Matematicheskaya teoriya borby za sushestvovanie [The mathematical theory of the struggle for existence].

Moscow, Nauka.

4. Хакен Г. 1980. Синергетика. М., Мир.

Haken G. 1980. Sinergetika[Synergetics]. Moscow, Mir.

5. Нахушев А.М. 1995. Уравнения математической биологии. М., Высшая школа.

Nakhushev A.M. 1995. Uravneniya matematicheskoy biologii [The equations of mathematical biology]. Moscow, Graduate School.

6. Нахушев А.М. 2003. Дробное исчисление и его применение. М., ФИЗМАТЛИТ.

Nahushev A.M. 2003. Drobnoe ischislenie i ego primenenie [The fractional calculus and its application]. Moscow, FIZMATLIT.