CC BY
34
Вестник СамГУ — Естественнонаучная серия. 2013. № 9/2(110)
УДК 517.9
215
ТРАЕКТОРИИ-УТКИ В ДИНАМИЧЕСКОЙ МОДЕЛИ ЧИСЛЕННОСТИ ПОПУЛЯЦИИ ВРЕДИТЕЛЕЙ ЛЕСА
© 2013 Е.В. Трошкина1
Работа посвящена исследованию модели динамики популяции гусеницы листовертки-почкоеда елового (Choristoneura fumiferana). Изучено поведение решений системы в зависимости от значений параметров, определены условия возникновения уточных решений.
Ключевые слова: популяционная модель, медленно-быстрая система, траектории-утки.
1. Математическая модель
Динамика популяции гусеницы листовертки-почкоеда елового (СЬоп81опеига ГитИегапа) и периоды наибольшей активности моделируются системой из трех обыкновенных дифференциальных уравнений, выведенных Д. Людвигом и др. в 1978 году [1]. Позже эта модель была упрощена Р.М. Мэйем [2]. Эта сокращенная модель имеет следующий вид:
йВ = гВ(1 - КБ) - в
Ш = - ^) - £В-
йТ—'^У1- К'Б> Г (а'Б)2 + В2- (1 1)
йБ _ „а/л Б \ ^п V1'1/
Здесь:
В, Б — плотность популяции вредителя и средняя площадь листвы, приходящаяся на одно дерево соответственно;
г, р — удельная скорость роста популяции и листвы соответственно; К', Бтах — пропускная способность численности популяции и средней площади листвы соответственно (пропускная способность означает емкость среды или максимально возможное значение переменной среды);
а'Б — продуктивность атак хищников, то есть значение, при котором начинается "хищничество" (чем меньше это значение, тем быстрее возрастает эффект "хищничества" при росте плотности популяции почкоеда); в — уровень насыщенности для хищника; £ — коэффициент пропорциональности.
Изменение площади поверхности листвы описывается вторым уравнением системы (1.1), где первое слагаемое отражает натуральный прирост листвы, описанный логистическим соотношением, а второе — потери, которые линейно пропорциональны плотности популяции вредителя.
хТрошкина Елена Вячеславовна ([email protected]), кафедра технической кибернетики Самарского государственного аэрокосмического университета (Национального исследовательского университета), 443001, Российская Федерация, г. Самара, ул. Молодогвардейская, 151.
2. Уточные решения
2.1. Предварительные данные
Пока в системе (1.1) выполняется соотношение для ее параметров р ^ r, данная система будет являться медленно-быстрой или сингулярно возмущенной системой.
Систему (1.1) можно записать в виде
eb = f (b,s), ( )
s = g(b,s), (2.1)
где £ = р — малый положительный параметр, s — быстрая переменная, b — медленная переменная системы.
При £ ^ 0 система (2.1) принимает вид
0 = f (b,s), (2.2)
S = g(b,s). (2.3)
На фазовом портрете (b, s) график (2.2) изображается кривой, которая называется медленной кривой Г. Части кривой, для которых выполняется неравенство
f (b,s) < 0, (2.4)
называются устойчивыми участками Г, если же выполняется противоположное неравенство
f (b,s) > 0, (2.5)
то такие участки называются неустойчивыми. Если положительный параметр £ мал, тогда траектория системы проходит вдоль устойчивого участка кривой, и когда она достигает £-окрестности точки (b, s), где
f (b, s) = 0, (2.6)
то совершает скачок. Как следствие, формируется предельный цикл классических релаксационных колебаний (рис. 1, а).
Если же после достижения точки P траектория системы еще некоторое расстояние проходит вдоль неустойчивой медленной кривой и только потом срывается на устойчивую часть Г, то такое решение называют уточным (рис. 1, b) [3].
Предположим, что обе правые части в (2.1) зависят от дополнительного параметра р. Тогда на параметрическом портрете (£,p) можно отметить область существования уточного решения. Ширина этой области имеет порядок O(e-1/E)
2.2. Результаты численного исследования
Введем безразмерные переменные
B s
b = ~j7i , s = —,t = рТ (2.7)
а smax smax
и преобразуем систему (1.1) к виду
£f = b а - ¿) - Y= f(b,s,a,1),
ds = s(1 - s) - \b = g(b,s,5). '
(0 р
\ Л.
V .
, \ *
0> г ■ р Л
Л | 1 Д,
у V -
V. \ -
Рис. 1. Траектории системы: а — график релаксационного колебания; Ь — график уточного решения
Система (2.8) содержит четыре положительных безразмерных параметра:
К' га!Бт„Т „ Р Р ,
а = = -¿Т^'6 = = (2.9)
а' р а г
Для качественного исследования получившейся системы уравнений (2.8) определим особые точки. Наиболее интересным оказывается случай, когда особая точка является точкой экстремума медленной кривой Г. Такая особая точка описывается набором из трех уравнений:
/ = О,
д(Ь,в,5) = 0, (2.10)
Э! = 0
эъ = 0
Заметим, что при изменении значения параметра а меняется положение особой точки на медленной кривой, а также при прохождении через бифуркационное значение меняется характер ее устойчивости и происходит бифуркация Андро-нова-Хопфа.
Эволюция решения в зависимости от значения параметра а представлена на рис. 2-4. Первый график в каждом отдельном случае представляет собой развертку решения по временной координате, второй же — медленная кривая и фазовый портрет системы, построенные при фиксированных значениях параметров.
При достижении параметром бифуркационного значения происходит бифуркация Андронова-Хопфа, то есть особая точка становится неустойчивым фокусом и от нее отделяется замкнутая траектория.
Медленная кривая Г состоит из трех частей. Нижняя и верхняя дуги являются устойчивыми ((2.4) выполняется), а средняя дуга — неустойчивая ((2.5) выполняется). Таким образом, замкнутая траектория, изображенная на рис. 3, имеет одновременно устойчивые и неустойчивые участки медленного движения, то есть является траекторией-уткой.
При дальнейшем увеличении параметра формируется релаксационный цикл (рис. 4).
Рис. 2. График решения системы (2.8). Особая точка — устойчивый фокус: а = = 8,1546; 6 =10; е = 0,1; 7 = 0, 66
а 10 20 19 40 50 1
Рис. 3. График уточного решения системы (2.8). Особая точка — неустойчивый фокус: а = 8,1565; 6 = 10; е = 0,1; 7 = 0, 66
О 10 20 }0 -И Я
Рис. 4. График решения системы (2.8). Предельный цикл: а = 9,1; 6 = 10; е =
= 0, 1; 7 = 0, 66
Заключение
Исследование показало, что модель динамики популяции гусеницы листовертки-почкоеда елового (Choristoneura fumiferana) обладает интересными харак-
теристиками. Главным образом, было исследовано существование уточного решения системы. Это явление объясняется тем, что особая точка, для которой происходит бифуркация Андронова-Хопфа, лежит на медленной кривой очень близко к точке срыва, разделяющей устойчивую и неустойчивую части кривой. Рождение замкнутой траектории делает возможным соединение устойчивого и неустойчивого интегральных многообразий в окрестности медленной кривой. Через такого рода точку и проходит уточное решение.
Литература
[1] Ludwig D., Jones D.D., Holling C.S. Qualitative analysis of insect outbreak systems: the spruce budworm and forest // Journal of Animal Ecology. 1978. № 47. P. 315-332.
[2] May R.M. Thresholds and breakpoints in ecosystems with a multiplicity of stable states // Nature. 1977. № 269. P. 471-477.
[3] Соболев В.А., Щепакина Е.А. Редукция моделей и критические явления в макрокинетике. M.: Физматлит, 2010. 319 c.
[4] Методы анализа нелинейных динамических моделей / М. Холодниок [и др.]. М.: Мир, 1991. 363 c.
Поступила в редакцию 18/XI/2013; в окончательном варианте — 19/XII/2013.
220
Е.В. TpowKuua
CANARD CURVES IN THE DYNAMIC MODEL OF FOREST INSECT PESTS
© 2013 E.V. Troshkina2
The paper is concerned with the investigation of the model of dynamics of population of caterpillar budworm (Choristoneura fumiferana). The system
solutions behavior depending on the parameters values has been investigated, conditions of origination of canard curves have been found.
Key words: population model, slow-fast system, canard curves.
Paper received 18/XI/2013. Paper accepted 19/XH/2013.
2Troshkina Elena Vyacheslavovna ([email protected]), the Dept. of Technical Cybernetics, Samara State Aerospace University, Samara, 443001, Russian Federation.
