CC BY
12
Математические структуры и моделирование 2019. №2(50). С. 79-84
УДК 004.9 : 631.4 + 519.6 DOI: 10.25513/2222-8772.2019.2.79-84
РАВНОВЕСНАЯ ДИНАМИКА НЕРАЗЛОЖИВШЕГОСЯ ОПАДА В ЛЕСНЫХ ЭКОСИСТЕМАХ
Л.А. Володченкова
к.б.н., доцент, e-mail: volodchenkova2007@yandex.ru А.К. Гуц
д.ф.-м.н., профессор, e-mail: guts@omsu.ru Омский государственный университет им. Ф.М. Достоевского, Омск, Россия
Аннотация. В статье исследуются равновесные состояния Нэша для нераз-ложившегося опада в лесных экосистемах в рамках теории дифференциальных игр. Показано, что в равновесии разложение опада идёт с меньшей скоростью.
Ключевые слова: Равновесие Нэша, лесная экосистема, опад, растительность, дифференциальные игры.
Введение
Динамику различных компонент х (биомасса, опад, грибы, гумус и пр.) лесных экосистем часто описывают системой дифференциальных уравнений
^ = $(г,г,иг,...,им), г = (х,у) € К2. аъ
с управляющими внешними факторами иг,...,и^,.
В качестве внешних управляющих факторов могут браться самые различные характеристики лесных компонент г. Динамика показателя г зависит от того, в каких границах изменяются управляющие факторы, которыми распоряжаются не только люди, но и то, что можно назвать природой.
Наши знания ограничены, поэтому любые математические модели являются упрощённым взглядом на изучаемые зависимости. Чаще всего принимается постоянным, неизменным во времени то, что на языке математиков именуется коэффициентами, входящими в выводимые дифференциальные уравнения.
Как выявить, насколько неточна изучаемая модель? Один из возможных способов предлагается в данной статье. Суть состоит в том, что коэффициенты рассматриваются как управляющие факторы, а управляют ими неизвестные игроки, которые являются агентами природы. В таком случае мы ищем такое оптимальное управление, которое трудно отнести к людям в силу его идеалистичности, неприемлемости для людей. Таким оптимальным управлением являются равновесия Нэша — в них каждый игрок бережно относится к интересам другого.
Затем сравним решение, предлагаемое системой с постоянными коэффициентами, и решение, получаемое в случае оптимального управления.
Для реализации нашего замысла мы найдём равновесие Нэша для системы уравнений, описывающих опад1, используя теорию дифференциальных игр [2].
1. Уравнения, описывающие опад
Пусть Ь — это свежий неразложившийся опад, Р — частично разложившийся опад. Тогда их динамика описывается системой дифференциальных уравнений [1, с. 86]
^ = -(Ь + к3)Ь,
М (1)
— = к3Ь - к2Р. ■ №
где коэффициенты к\,к2,к3 описывают функционирование комплекса грибов и микрофауны. Они оценены на основе известных экспериментов литературных данных по лабораторным экспериментам по определению скоростей разложения органических остатков, а также на основе анализа многочисленных экспериментов по разложению опада в контролируемых условиях в лаборатории биохимии почв Биологического института СПбГУ [1, с. 83].
2. Алгоритм нахождения равновесий Нэша
Рассматриваем уравнение (1) как дифференциальную игру с тремя игроками
и = к\ + к3, и2 = к3, и3 = к2:
¿Ь
-ж = -и1Ь,
¿Р
— = и2Ь - изР, №
(2)
Будем рассматривать игру с ненулевой суммой, поскольку по нашему замыслу «выигрыши» наших игроков слабо связаны.
Если игрок формирует «своё» управляющее воздействие в виде только функции времени и(Ь) на всю продолжительность игры, то и(Ь) — это программное управление игрока. Ранее мы называли его, используя термин «управление». Однако игрок может выбирать своё управление в зависимости от того, в каком положении (Ь,Р) в момент времени Ь находится система. В таком случае игрок конструирует управляющее воздействие в виде функции и( , Ь, Р), зависящей уже от позиции {Ъ, (Ь,Р)}, и для и(Ь,Ь,Р) используется термин позиционное управление игрока [3]. Часто пишут просто и(Ь,Р).
Мы будем искать позиционное управление, позиционное равновесие Нэша.
1Опад — отмершие части растений (ветки, листья и др.), опавшие на поверхность почвы или дно водоёма.
Математические структуры и моделирование. 2019. №2(50)
81
Для дифференциальной игры Ж-игроков
^ = !(*) + ¿Л- , т = 0,
3 = 1
г=(х, у) € К2, щ € К,
N
^,щ,...,им) = / + (щ)2](И, (г=1,...,Ж),
о *=г
где числа
< > 0, Кгг > 0, Пг1 ^ 0, существование равновесий Нэша
Зг(и\,и*2,и*^, ...,и%) ^ Зг(и**,и*2, ...,и*-г,иг,и*+1, ...,П*Ы), Ущ, г=1,...,Ы (3)
сводится к крайне сложной задаче отыскания решения Уг(г) нелинейного уравнения Гамильтона-Якоби
1 М
(ЪУг)Т!(г) + <&) - (?№з)-1 (Зз(*))ТФЪ)+
3 = 1
1 М
+ 1 Е^ Ъ )Т9э [( Язз )-1]2(9з (?))Т (V Ъ) = 0, (4)
3 = 1
v V = ) , (V V?' = ((V)x, (Ъ)у),
по которому строится равновесие Нэша [2, Theorem 10.4-2, утверждение b.]:
u*(z) = Ui(Vi(z)) = -Ir-l(g*(z))T(vV), i = l,..., N. (5)
3. Нэшевское равновесие системы для опада
В нашем случае N = 3, и рассматриваем
Rll = R22 = R33 = 1, Rij =0 (i = j), z=(L,F).
Полагая, что
Vi(z) = V2(z) = V3(z) = 1(L2 + F2) > 0
и подставляя эти функции в уравнения Гамильтона-Якоби, получаем
< = \ьА + \ь2Р2 + \р4 > 0,
< = -ь4 + -ь2р2 + -р4 > о, (6)
/2 = -Ь4 + -Ь2Р2 + -/З 2 4 2
<3 = 2ь4 + \ь2р2 + 4^4 > 0.
Следовательно, если < выбрать именно такими, то уравнения Гамильтона-Якоби выполняются.
Поэтому по теореме 10.4-2 из [2] имеем равновесие Нэша
«I = 2ь2, и* = -2ьр, из = 2р2, (7)
найденное по формулам (5).
Выигрышные / проигрышные функции поэтому имеют вид:
^(Ь, Р, ь + кз, кз, к2) = I [Я1(х) + (к1 + кз)2](И,
о
32(1, Р, к1 + кз, кз, к2) = I Шх) + кЗ^сН, (8)
о
3з(Ь, Р1 + кз, кз, кз) = J [<<з(х) + кз]сИ.
о
Эволюция опада Ь в случае равновесия Нэша (7) находится посредством подстановки (7) в первое уравнение системы (2) с последующим интегрированием:
= -«1Ь = --ьз
(И 1 2
или
Ь(1) = 1 =. (9)
\/2(Ь + с о п в Ь)
Если непосредственно проинтегрировать первое уравнение системы (1), то получим, что
Ь(Ь) = со'пМ^ е-(к1+кз)\ (10)
4. Заключение
Из уравнений (9), (10) видим, что в случае равновесного управления Нэша разложение опада идёт по степенному закону, а в случае постоянных коэффициентов по экспоненте, т. е. гораздо быстрее. Другими словами, учёт «интересов» всех игроков:
Математические структуры и моделирование. 2019. №2(50)
83
- игрока к\ — скорость минерализации свежего органического материала (разложение в органическом горизонте);
- игрока к2 — скорость минерализации КГВ (комплекс гумусовых веществ);
- игрока к3 — скорость образования КГВ,
ведёт к замедлению скорости разложения опада. Природа более «милостлива» к опаду — она в каждый момент времени выбирает коэффициенты2 к\, к2, к3 так, что они коррелируют свои значения с текущим состоянием неразложившейся и разложившейся частей опада.
Литература
1. Моделирование динамики органического вещества в лесных экосистемах / отв. ред.
B.Н. Кудеяров: Ин-т физ.-хим. и биолог. проблем почвоведения РАН. М. : Наука, 2007. 380 с.
2. Lewis F.L., Vrabie D.L., Syrmos V.L. Optimal control. New Jersey : John Wiley & Sons, Inc., 2012. 540 p.
3. Тынянский Н.Т., Жуковский В.И. Дифференциальные игры с ненулевой суммой (кооперативный вариант) // Итоги науки и техн. Сер. Мат. анал. 1979. Т. 17.
C. 3-112.
EQUILIBRIUM DYNAMICS OF UNDECOMPOSED LITTER IN FOREST
ECOSYSTEMS
L.A. Volodchenkova
Ph.D. (Biology), Associate Professor, e-mail: volodchenkova2007@yandex.ru
A.K. Guts
Dr.Sc. (Phys.-Math.), Professor, e-mail: guts@omsu.ru Dostoevsky Omsk State University
Abstract. In the article the Nash equilibrium states for undecomposed litter in forest ecosystems within the framework of the theory differential games are investigated. It is shown that in equilibrium the decay of litter goes at a lower speed.
Keywords: The Nash equilibrium, forest ecosystem, litter, vegetation, differential games.
References
1. Modelirovanie dinamiki organicheskogo veshchestva v lesnykh ekosistemakh. Otv. red. V.N. Kudeyarov, In-t fiz.-khim. i biolog. problem pochvovedeniya RAN, Moscow, Nauka Publ., 2007., 380 p. (in Russian)
2В [1, S 2.3.4] косвенным образом это подтверждается в ходе наблюдения зависимости скоростей минерализации и гумификации от содержания азота и зольности в опаде.
