CC BY
7
Вестник КРАУНЦ. Физ.-мат. науки. 2018. № 3(23). C. 124-130. ISSN 2079-6641 DOI: 10.18454/2079-6641-2018-23-3-124-130
УДК 517.927
ОБ ОДНОЙ МОДЕЛИ ДИНАМИКИ ЧИСЛЕННОСТИ НАСЕЛЕНИЯ С УЧЕТОМ ПОЛОВОЙ СТРУКТУРЫ
Ф. М. Лосанова, Р. О. Кенетова
Институт прикладной математики и автоматизации, 360000,
г. Нальчик, ул. Шортанова, 89А
E-mail: losanovaf@gmail.com, raisa.kenetova@mail.ru
В данной работе проведен анализ математической модели, которая описывает динамику популяции с учетом половой структуры. Для определения плотности семейных пар найдено нелинейное интегральное уравнение в свертках, обсуждены подходы к решению.
Ключевые слова: модель с половой структурой, динамика численности населения, плотность семейных пар, нелинейное интегральное уравнение, корень свертки.
@ Лосанова Ф.М., Кенетова Р. О., 2018
MSC 34L99
ON ONE MODEL OF POPULATION DYNAMICS WITH REGARD TO SEX STRUCTURE
F. M. Losanova, R. O. Kenetova
Institute of Applied Mathematics and Automation, 360000, Nalchik, Shortanova st., 89A, Russia
E-mail: losanovaf@gmail.com, raisa.kenetova@mail.ru
In this paper, we consider a mathematical model describing population dynamics in view of the sexual structure. We obtain a nonlinear convolution equation for determining the density of family pairs, and discuss approaches to its solution.
Key words: sexual structure model, population dynamics, density of family pairs, nonlinear integral equation, convolution root.
© Losanova F.M., Kenetova R.O., 2018
Введение
В настоящее время теория структурированных биологических популяций является достаточно активно развивающейся областью науки. Она находит важные приложения в биофизике и медицине. Значительный прикладной интерес представляют также те задачи теории структурированных биологических популяций, которые касаются связи внутрипопуляционной дифференциации, обусловленной внутренней структурой их элементов, с пространственно-временными аспектами их эволюции.
Рождаемость и смертность в определенных пределах могут быть описаны линейными процессами, а образование семейных пар является существенно нелинейным феноменом. Задача моделирования процесса формирования пары была названа "двухполой задачей"[1]. Построением и исследованием систем обыкновенных дифференциальных уравнений, а также интегральных уравнений, описывающих возрастную структуру занимались многие авторы, такие как D.G. Kendall [2], N. Keyfitz [3], D.D McFarland [4], B. Parlett [5], J.H. Pollard [1], A.G. Fredrickson [6].
0.В. Староверов [7] предложил общую модель, учитывающую возрастную структуру и предпочтительный возраст вступления в брак.
Недостатком моделей, которые описываются обыкновенными дифференциальными уравнениями без возрастной структуры является то, что они подразумевают, что даже новорожденные могут незамедлительно образовывать пары. Период созревания можно моделировать только путем введения в дальнейшем полу-дискретных моделей интегральных уравнений в бесконечных пространствах [8].
В работе [8] представлена система трех связанных обыкновенных дифференциальных уравнений, которая включает в себя многие модели из более ранних, в виде особых случаев. Исследуемая в [8] модель нелинейна, стационарные решения представлены в экспоненциальном виде, которые описывают популяции с постоянной социальной структурой и соотношением полов, которые растут (или распадаются) экспоненциально по времени. "Рождение"следует интерпретировать как возраст особи достигшего совершеннолетнего возраста.
1. Описание модели
В данной работе проведен анализ математической модели, в основе которой лежит модель, предложенная в работе [8]
x = az — ßx — F (x, y), (1)
y = bz — vy — F (x, y), (2)
Z = F (x, y) — cz, (3)
x(0) = xo, y(0) = yo, z(0) = zo. (4)
Здесь x = x(t), y = y(t) - численность неженатых мужчин и незамужних женщин; z = z(t) - число семейных пар; 0 < t < T; ß, v — коэффициенты интенсивности
смертности мужчин и женщин (состоящих и не состоящих в браке), a, b, c — параметры модели, F(x,y) — интенсивность (скорость) образования пар в результате контактов между неженатыми мужчинами и незамужними женщинами. Причем функция F(x,y) должна удовлетворять следующим свойствам:
1) F(0, y)= F(x, 0)= 0, для всех x > 0, y > 0;
2) если и > 0, & > 0, то F(x + и,y + &) > F(x,y), и, & = const.
В нашей модели предлагается учитывать интенсивность образования семейных пар путем введения нелинейного оператора
F (x, y) = x * y * a,
где
(x * y)(t) = / x(t — s)y(s)ds J 0
— свертка Лапласа функций x(t) и y(t), a(t) — заданная неотрицательная функция. Очевидно, что оператор F(x,y) обладает свойствами 1) и 2). Кроме того, в предложенном виде модель более проста для анализа, а выбор плотности a(t) предоставляет более широкие возможности при описании конкретных данных.
2. Анализ модели
С помощью уравнения (3), (1) и (2) можно переписать в виде
x + jix = (a — c)z — z, (5)
y + vy = (b — c)z — Z. (6)
Предполагая правые части известными, уравнения (5) и (6) перепишем в виде
г t
x(t) = x0e—'ш + / ((a — c)z — Z)e—'u(t—n, (7)
0
y(t) = y0e—vt + f ((b — c)z — Z)e—v(t—n . (8)
0
Интегрируя по частям уравнения (7) и (8), получим
x(t) = (x0 + z0)e—M 1 + (a — c + д )z * e—д — z, (9)
y(t) = (y0 + z0)e—vt + ((b — c + v )z * e—vt — z. (10)
Уравнение (3) также можно переписать в виде
z(t)= z0e—ct + F(x, y) * e—ct. (11)
2.1. Случай д = V
Будем считать, что д = V. Вычислим ^(х,у)
^(х,у) = [(хо + го)е-дг + (а - с + д)г * е-дг - г]* *[(уо + zо)e-vг + ((Ь - с + V)г * е-^ - г] * а =
е-Дг ¿-дг
= (хо + го)(уо + го)-* а + [(хо + го)(Ь - с + V)--(12)
д - V 1 д - V
¿-V; - е-дг
- (хо + го)е-дг + (а - с + д)(уо + го)----(уо + го)е-* г * а+
д ^
+ [(а - с + д)(Ь - с + V)—д—V--(а - с + д)е-дг - (Ь - с + V+ 1] * г * г * а.
2.2. Случай д = V
Учитывая, что е-дг * е-дг = ге-дг, формула (12) запишется в виде
^(х,у) = [г * г * (а^ге-дг - (а + Ь^е-дг + 1)+
+г* (((уо + го)а + (хо + го)^)^-дг - (хо + уо + 2го)е-дг+ + (хо + го) (уо + го)ге-дг] * а.
3. Интегральное уравнение для плотности семейных пар
Подставляя (12) в (11) получим
г(г) = (г * г * В)(г) + (г * С) (г) + ), (13)
где
) = [Bie-(v+c)t - В2е-(д+c)t + e-ct] * а, (а - c + v )(b - c + v) (а - c + д)(Ь - c + д)
B1 =-, B2 = -
д - v д - v
C(t) = [Cie-(v+c)t - С2е-(д+c)t] * а,
(а - c + v )(yo + zo) + (b - c + v )(xo + zo) д - v ,
(b - c + д)xo - (а - c + д )yo + zo (b - а)
Ci = C2 =
D(t)= Di[e-(v+c)t - е-(д+c)t] * а + zoe-ct, Di =
д - v
(xo + zo)(yo + zo)
д - V
Перенеся слагаемое (г * С)(г) уравнения (13) в левую часть и считая оставшуюся правую часть известной, приходим к уравнению
г(г )= ^ (г ) + (я * ^ )(г), (14)
где ^(*) = (г * г * )+ ), ) = £ С*и(?) - резольвента ядра С(*). Здесь С*и(?)
п=0
означает п-ую степень свертки: С*0 = С(?), С*и(?) = (С*Си-1)(?). Подставив (14) в (13), приходим к интегральному уравнению
) = (г * г * р)(,) + Q(í), (15)
где Р(*) = ) + (Л * 5) (*), Q(í) = ) + (Л * Л) (*).
4. О методах решения уравнения (15)
Таким образом, для определения плотности семейных пар г(0 мы пришли к нелинейному интегральному уравнению (15). Далее обсудим некоторые подходы к решению этого уравнения.
4.1. Сведение к вопросу о нахождении корня свертки
Обозначим через = (г*Р)(?). Далее, подействуем на обе части уравнения (15) оператором Р*, то есть возьмем свертку с функцией Р(?), получим
) = (м * ) + ^ * Р)(*). (16)
Или, после простых преобразований
МО -1) *(М0 -1) = 1 - Ш * Р)(0. (17)
Таким образом, вопрос нахождения решения уравнения (15) сведен к нахождению квадратного корня свертки от правой части (17). То есть, необходимо найти функцию ) такую, что
(Л*Л)(г) = 1 - (Q*Р)(г).
Определив Л(?) получим, что
М(г ) = 1 ± Л(г).
Далее, ) ищется как решение уравнения (г * Р)(?) = ), которое является интегральным уравнением первого рода типа свертки, и может быть обращено, например, если Р(0) = 0 и Р'(?) е С[0, Т].
4.2. Линеаризация
Для поиска приближенного решения (15) можно применить метод линеаризации, предложенный в работе [9, стр. 68]. Для этого вместо исходного уравнения рассматривается линеаризованное нагруженное уравнение
) = (г * г * р)(,) + Q(í), (18)
1 Т
где г = 1 /0 г^О^у — усреднение искомого решения. Полученное уравнение является линейным и его решение можно найти (например методом последовательных приближений), считая г известным. Далее, усредняя найденное решение, получаем функциональное уравнение для нахождения г. Данный метод позволяет находить приближенное решение, которое в определенных случаях оказывается достаточным для проведения качественного анализа модели.
Заключение
В данной работе исследована математическая модель, описывающая динамику популяции с учетом половой структуры, которая представляет собой систему трех связанных обыкновенных дифференциальных уравнений. Решая их мы приходим к нелинейному интегральному уравнению (15) в свертках для определения плотности семейных пар z(t). В работе обсуждены методы решения данного уравнения. Один из них сводит вопрос о решении нелинейного уравнения к вопросу о нахождении квадратного корня свертки. Второй метод позволяет находить приближенные решения для проведения качественного анализа модели.
Список литературы
[1] Pollard J.H., Stochastic processes and population growth, Cambridge University Press, 1973.
[2] Kendall D.G., "Models for pair formation in bisexual populations", Roy. Statist. Soc., Ser B 2 (1949), 230-264.
[3] Keyfitz N., "The mathematics of sex and marriage", Proceedings of the Sixth Berkeley Symposion on Mathematical Statistics and Probability. V. IV, Biology and health, 1972, 89-108.
[4] McFarland D. D., "Comparison of alternative marriage models", Population dynamics, Academic Press, New York-London, 1972, 89-106.
[5] Parlett B., "Can there be a marriage function?", Population Dynamics, Academic Press, New York-London, 1972, 107-135.
[6] Fredrickson A. G., "A mathematical theory of age structure in sexual populations: Random mating and monogamous marriage models", Math. Biosci., 10 (1971), 117-143.
[7] Staroverov, O. V., "Reproduction of the structure of the population and marriages", Ekonomika i matemati6eskije metody, 13 (1977), 72-82.
[8] Hadeler K. P., Waldstatter R., Worz-Busekros A., "Models for pair formation in bisexual populations", J. Math. Biol., 26:6 (1988), 635-639.
[9] Нахушев А.М., Нагруженные уравнения и их применение, Наука, М., 2012, 232 с. [Nahushev A.M., Nagruzhennye uravneniya i ih primenenie, Nauka, M., 2012, 232 pp.]
Список литературы (ГОСТ)
[1] Pollard J. H. Stochastic processes and population growth. Cambridge University Press, 1973.
[2] Kendall D. G. Models for pair formation in bisexual populations // Roy. Statist. Soc. 1949. Ser B 2. pp. 230-264.
[3] Keyfitz N. The mathematics of sex and marriage // Proceedings of the Sixth Berkeley Symposion on Mathematical Statistics and Probability. 1972. IV: Biology and health. pp. 89-108.
[4] McFarland D. D. Comparison of alternative marriage models. Population dynamics. New York-London: Academic Press, 1972. C. 89-106.
[5] Parlett B. Can there be a marriage function?. Population Dynamics. New York-London: Academic Press, 1972. pp. 107-135.
[6] Fredrickson A. G. A mathematical theory of age structure in sexual populations: Random mating and monogamous marriage models // Math. Biosci. 1971. vol. 10. pp. 117-143.
[7] Staroverov O. V. Reproduction of the structure of the population and marriages // Ekonomika i matematicheskije metody. 1977. vol. 13. pp. 72-82.
[8] Hadeler K. P., Waldstatter R., Worz-Busekros A. Models for pair formation in bisexual populations // J. Math. Biol. 1988. vol. 26. no. 6. pp. 635-639.
[9] Нахушев А. М. Нагруженные уравнения и их применение. М.: Наука, 2012. 232 с.
Для цитирования: Лосанова Ф.М., Кенетова Р. О. Об одной модели динамики численности населения с учетом половой структуры // Вестник КРАУНЦ. Физ.-мат. науки. 2018. № 3(23). C. 124-130. DOI: 10.18454/2079-6641-2018-23-3-124-130
For citation: Losanova F. M., Kenetova R. O. On one model of population dinamics with regard to sex structure, Vestnik KRAUNC. Fiz.-mat. nauki. 2018, 23: 3, 124-130. DOI: 10.18454/2079-6641-2018-23-3-124-130
Поступила в редакцию / Original article submitted: 08.06.2018
