довой комбинаторной оптимизации // Экономика и ма-тем. методы. 2000. Т. 36, №2. С. 141-144. 11.Емец О.А., Недобачий С.И., Колечкина Л.Н. Неприводимая система ограничений комбинаторного многогранника в дробнолинейной задаче оптимизации на перестановках // Дискретная математика. 2001. Т. 13. Вып. 1. С. 110 - 118. 12.Ємець О О., Роскладка О.В., Недобачій С.І. Незвідна система обмежень для загального многогранника розміщень // Укр. мат. журн. 2003. Т. 55, №1. С. 3-11. 13.Черненко О.О. Дослідження множини допустимих розв’язків задачі оп-тимізації дробово-лінійної функції на евклідовій комбінаторній множині розміщень // В кн. X міжн. наук. конф. ім. ак. М.Кравчука (13-15 травня 2004 р., Київ): Матеріали конференції. К., 2004. С. 548. 14.Ємець О.О., Барболіна Т.М. Розв’язування задач нелінійної умовної оптимізації на розміщеннях методом відсікання // Укр. мат. журн. 2003. Т. 55, №5. С. 604-612. 15.Емец О.А., Барболина Т.Н. Решение линейных задач оптимизации на размещениях методом отсечения // Кибернетика и системный анализ. 2003. №6. С. 131-141. 16. Емец О.А., Барболина Т.Н. Решение задач евклидовой комбинаторной оптимизации методом построения лексикографической эквивалентности // Кибернетика и системный анализ. 2004. №5. С. 115-
УДК517.9 "
УСТОЙЧИВОСТЬ СТАЦИОНАРНЫХ СОСТОЯНИЙ СМЕШАННОЙ КУЛЬТУРЫ С ВИДОМ ВЗАИМОДЕЙСТВИЯ ТИПА КОММЕНСАЛИЗМ
ЯЛОВЕГА И.Г.______________________________
Проводится анализ устойчивости стационарных состояний смешанной культуры с видом взаимодействия типа комменсализм с помощью метода Ляпунова. Показывается динамика системы во времени. Строятся графики областей устойчивости при изменении скорости протока и концентрации лимитирующего субстрата.
1. Введение
Поведение смешанных культур, смесей организмов различных типов, имеет большое значение для экологии микроорганизмов в почве, воде, при изучении болезней и порчи продуктов [2,5]. Смешанные культуры имеют также большое значение в приготовлении пищевых продуктов брожения и изготовления микробных продуктов. В качестве примеров промышленных процессов с участием смешанных культур микроорганизмов можно указать на сыроварение и биологическую очистку сточных вод. Более того, в естественных системах смешанные популяции микроорганизмов являются скорее правилом, чем исключением. Из-за сложности поведения смешанных культур использование математических моделей различных систем для описания и предсказания поведения культуры прио бретает осо бенное значение.
Одной из областей применения дифференциальных уравнений в моделировании микробиологических сообществ является анализ возможных стационарных состояний смешанной культуры [6,7].
BE, 2005, 1 4
125. 17. Емец О.А., Колечкина Л.Н. Решение задач оптимизации с дробно-линейными целевыми функциями и дополнительными линейными ограничениями на перестановках // Кибернетика и системный анализ. 2004. №3. С. 30-43.
Надійшла до редколегії 20.12.2005
Рецензент: д-р фіз.-мат. наук Лагно В.І.
Ємець Олег Олексійович, д-р фіз.-мат. наук, професор, завідувач кафедри математичного моделювання та соціальної інформатики Полтавського університету споживчої кооперації України. Наукові інтереси: комбінаторна оптимізація. Захоплення: колекціонування марок. Адреса: Україна, 36003, Полтава, а/с 1671, тел. 8-050-30-47-156, роб. (8-0532)509-204, e-mail: slemets@e-mail.pl.ua
Черненко Оксана Олексіївна, лаборант кафедри економічної кібернетики Полтавського університету споживчої кооперації України. Наукові інтереси: дослідження задачі з дробово-лінійною функцією цілі на розміщеннях та її розв’язування. Захоплення: математика, психологія. Адреса: Україна, Полтава, вул. Г.Сталінграда, 9, кв.66, тел. 8-097-452-62-59, роб. (8-0532)509-205, дом. 66-81-54, e-mail: chemenko7@ukr.net
Имеющийся в литературе фактический материал отражает огромное разнообразие конкретных типов взаимодействий (конкуренция, амменсализм, комменсализм, симбиоз и т. д.), осуществляемых через выделение или потребление различных веществ (метаболитов) в среду и из среды [1,2,5]. Положительное или отрицательное влияние этих метаболитов, а также лимитирующих рост субстратов приводит к различным исходам взаимодействий: вытеснению, сосуществованию или доминированию. Эти исходы дают представление об обширной гамме возможных вариантов в динамике популяций одного трофического уровня.
При комменсализме в случае смешанных культур, состоящих из двух штаммов микроорганизмов, определенные преимущества получает второй вид. Широко распространен вариант комменсализма, когда один вид продуцирует соединения, ускоряющие рост организмов другого вида. Конечные продукты жизнедеятельности одного вида, обеспечивающие эффект комменсализма, довольно многочисленны. В зависимости от конкретной системы продуцируемое соединение может служить источником энергии или углерода для другого вида [2,5].
Проточное, или непрерывное, культивирование широко применяется в промышленном производстве, так как обеспечивает непрерывный вывод из установки однородной массы клеток. Кроме того, существует возможность в некоторых пределах управлять параметрами системы для получения оптимальных условий выращивания. Проточный культиватор, построенный на принципе хемостата, обладает устойчивым состоянием равновесия, когда масса вымываемых клеток компенсируется нарождающимися. Изменение скорости протока или концентрации питательного субстрата приводит к новому состоянию равновесия.
149
2. Постановка задачи
Рассмотрим математическую модель, описывающую производство смешанной культуры, которая состоит из двух видов популяций, в проточной системе, а именно в хемостате. Предполагается, что виды лимитируются одним субстратом, а также то, что в присутствии вида 1 вид 2 растет быстрее, т. е. наблюдается взаимодействие типа комменсализм. Исходя из указанных взаимоотношений между видами, можно составить систему уравнений, описывающую динамику компонентов в проточной системе:
( .. с Л
2. Х2 = 0, S2 =
3. уз = 0, S3 =
DKy у 2 Sq - S
P у _ D a у
DKx x3 So _ S:
P x _ D a x
(2)
4. S4
о DKx
S4 = Х
(
Х4 =-Р x - D’ у
D -
Р yS4
V
S + K
у;
Х =
^ Xs
s+Kx
Ц yS
- D
У 4 =
Х,
d(Sq -S4 ~аxx4) a y(D-yx4)
S + K
- D
■yx
у
(1)
S = d(Sq -s)_а?
Ц Xs
s + Kx
x-а
у
Ц yS
S + K
■у.
у
Здесь S - концентрация субстрата; х, у - концентрации биомассы 1 и 2 в культиваторе;
D - скорость протока, или ско- pxS
рость разбавления; S0 - концен- s + Kx трация субстрата, поступившего
Для исследования устойчивости особых точек рассматривают линеаризованную систему дифференциальных уравнений, которая описывает движение вблизи положения равновесия. Характер устойчивости особых точек установим, пользуясь методом Ляпунова [4,6,7]. Характеристический определитель системы (1) имеет вид:
в культиватор; аЛ а,,1- эко-
- D-X
ТУ
а Х Р XS
Р yS
S + K
Х ’ "У
номические коэффициенты биомассы 1 и 2, показывающие, какая часть поглощенного субстра- S + Kx
та идет на приращение биомассы; р.Х, Ру - максимальные удельные скорости роста для двух видов; Kx, K, - константы Михаэлиса по субстрату S; у - коэффициент воздействия биомассы 1 на биомассу 2, в нашем случае комменсализма.
В задачах прикладного характера важно найти условия для совместного сосуществования в субстрате двух культур. Исследование автономной системы на устойчивость при изменении скорости протока или концентрации питательного субстрата дает возможность для оптимального управления процессом производства смешанных культур. Для исследования поведения модели во времени проведем анализ возможных стационарных состояний.
3. Анализ устойчивости стационарных состояний
Система (1) имеет четыре особые точки:
1. Х1 = 0, у1 = 0, S1 = Sq.
PxDK у
- D + yx - X
PxKxx
(S + Kx)
P yKyy
2
у
2
^ у P уС
S + K
- D -
у
(S + K,)
axPxKxx ^“уРуКуу
(S + Kx)2 (S + Ky)2
-X
= 0 . (3)
Для четырех особых точек ( 2 ) соответствующие им характеристические определители имеют вид:
1.
P xS0
S0 + Kx
- D-X
0
axPxS0 S0 + Kx
P yS0 S0 + K,
■-D-X
0у ^ у ^ YS0 S0 + K
- D-X
у
P xS0
V S0 + Kx
- D-X
P yS0
VS0 + Ky
У
- D-X
(-D-X).
(4)
DKy +p ,Kx-DKx У(PyS0 _ S0D _ DKy) ^ у(Ру _ D) a x P xDKy DKy +p ,Kx —DKx
--D-X Q
-X
(pySQ - SqD - DKy )(цу - D)
-ayD -D-
^ yP yKy
(pySQ - SqD - DKy )(цу - D)
P yK
-X
у
2
3
1
у
0
0
0
0
2
Q
1 50
BE, 2005, 1 4
Ц xDK
Л
y
DKy +p, yKx-DKx
■-D-X
-X
- D -
y • y x x у
(цySo - SqD - DKy)(py - D)
a yD
Ц yKy
(ЦySo - SoD - DKy )(py - D)
a
y Ц yKy
-X
\
(5)
-і/
-X
0
-a xD
p yDKx
DKx + pxKy-DKy
^ y P yDKx
-D + y| So--DK^-\-X
Ц x D
w
0
DKx + p xKy-DKy
Ц yDKx
(
- — D +y
So -
DKx
V
H- x D J
\ \
-X
DKX + pxKy-DKy
x(-X[-D-ax(pxSo -SoD-DKxX^x -D)-\l+ (6)
+ axD(pxSo - SoD -dKxXh-x - D)),
здесь w = (p xSo - SoD - DKx)(px - d).
4.
-X
У y
-a xD
o
-X
y Ц yDKx
DKx + pxKy-DK
u
v
-D-a xu-a yV-X
где u =-
D(H-x -D2(DKx +ЦxKy _DKy -pyKJ
y(dkx + P xKy - dkJ :
V =
H"yKy (цx D) y
(DKx + p xKy - DKy)'
Е) ^1,2 = ІЮ , 8 < o , ^1,2 є C, X з > o , X з є R -
Ж) ^1,2 = 8 +iro , S > o, ^1,2 є C, Xз < o , Xз є R -
З) X12 =8 +k> , 8> o , ^1,2 є C, X 3 > o , X3 є R - (8)
Из всех восьми случаев только два (А и Д) определяют устойчивые состояния. В случае А особая точка является устойчивым узлом, в случае Д - устойчивым фокусом.
4. Динамика развития системы и области устойчивости численно заданной модели при изменении параметров
Прикладной интерес в изучении смешанной культуры представляют устойчивые состояния нетривиальных особых точек- Для анализа системы (1) разработано программное обеспечение, позволяющее показать как динамику развития системы, так и области устойчивости численно заданной модели при изменении параметров. Возможность сосуществования и управления равновесным составом популяции непосредственно следует из анализа положительных значений стационарных компонентов системы.
Рассмотрим поведение системы (1) при следующих параметрах:
D = o,5, у = o,1, So = o,5, рx = 1,2, Цy = °,8,
Kx = o,oo4, Ky = o,o2 , ax = o,o5 , ay = o,3 . Входные данные:
-D-axw-X
x = 2, y = 1, S = 3.
(9)
|2 , y = y 4 . (7)
Значения X і , полученные из решений уравнений, определяют характер движения вблизи особых точек исходной нелинейной системы [4,6,7]. Возможны следующие варианты:
A) X1 <o, X2 <o, Xз <o, Xі єR.
Б) Х1 > o , X 2 < o, X з < o , X і є R.
B) Х1 >o, X2 >o, Xз <o, Xі єR.
Г) ^1 >o, X2 >o, Xз >o, Xі єR.
Д) ^1,2 = i® , 8 < o , ^1,2 є C, Xз < o , Xз є R .
Система (1) с данными параметрами имеет четыре особые точки (2):
1. x1 = o , y1 = o , S1 = o,5 .
2. x2 = o, y 2 = 1,556, S2 = o^.
3. x3 = 9,94з, уз = o, S3 = o^. (Ю)
4. x4 = 4, У4 = 4,952, S4 = o^.
Для каждой особой точки найдены Xі из (з) - (7):
1. Х1 = o,69 , X 2 =-o,5, X з = o,269 .
2. Х1 = o,571, X2 = -o,5 , Xз = -2,625 .
3. Х1 = o,594 , X2 =-o,5, Xз =-5o,75 .
4. Х1 =-o,111 + o,282i, X2 =-o,111-o,282i,
Xз = -66,195 .
Из полученных значений видно, что только одна особая точка является устойчивой. Четвертая особая точ-
3
o
y
BE, 2oo5, 1 4
1 5 1
ка, которая является нетривиальной, отвечает случаю Д из (8) и является устойчивым фокусом.
Полученные результаты были найдены с помощью теории устойчивости. При помощи численных методов найдено решение системы нелинейных уравнений (1), а именно, динамика состояний. На рис.1 представлено окно, в котором при заданных параметрах (таблица System parameters) и начальных условиях (Initial conditions) выведены решения системы (1) в таблице под Current state, полученные с помощью численных методов. Фазовый портрет системы в проекции на плоскость (x,y) представлен на рис.2.
5MR: D5P Control/Debug window
Sterns control I Stationery slate diagram | Output properties ]
*]
System selection: ] Unknown dynamic system (ODE system, 34h order)
pSystem parameters: 1
D 0.500 ■
Gamma 0.100
SO 0.500
Д1рЬз_к 0.050
АІрЬз_у 0 300
M_ к 1.200
К_и 0004
М_у 0.000
К_у 0.020
X Cared I
V AppV
Initial conditions і Current state:
к 2.000 14.000
У 1.000 4.962
$ Зі 000 0.003
Рис. 1. Решение системы (1) при численно заданных параметрах (9)
Рис. 2. Фазовый портрет системы (1) в плоскости (x,y) при численно заданных параметрах (9)
На фазовом портрете системы видно, что в точке (4;4,952) она входит в стационарное состояние и
данная точка является устойчивым фокусом. Из рис. 1 и 2 видно, что результаты, полученные численно, а именно стабилизация системы и значения устойчивого состояния, совпадают с результатами и значениями, полученными с помощью теории устойчивости.
При проточном культивировании в хемостате легко регулируемыми величинами являются скорость протока D и концентрация субстрата S0 .Варьируя значения этих параметров, возможно управлять процессом выращивания смешанной культуры в производственных целях [2]. При изменении параметров возможно изменение характера устойчивости, а также переход системы в качественно новое состояние, например, из устойчивого в неустойчивое.
С помощью программного обеспечения построены диаграммы, показывающие изменение поведения модели при вариации входных параметров. Так, на рис.3 представлена диаграмма в плоскости параметров D и S0 , показывающая области устойчивости и неустойчивости нетривиальной особой точки x4 = 4,
View Options
4 I %L 0 D
Stationary states |x, (D) |y, (SO) ^
0.3 0.5 0.7
Г- - г-
О О
LЛ Ю
О п 0.
ГЭ
О
0.3 0.5 0.7
y 4 = 4,952, S4 = 0,03. D и S0 изменяются в промежутке от 0,25 до 0,75.
Рис. 3. Диаграмма областей устойчивости и неустойчивости в плоскости параметров D и S0 для системы (1) при численно заданных параметрах (9)
На рис.3 видно, что при такой вариации параметров существуют три области: верхняя область (С—), где 7-1,2 =8+ iro , 5<0, 7-1,2 є C, X3 <0, X3 є R, что удовлетворяет случаю Д, т. е. в этой области нетривиальная особая точка является устойчивым фокусом; средняя область (R—), где 7-1 < 0 , X 2 < 0, X 3 < 0 , X і є R, что удовлетворяет случаю А, т. е. в этой области нетривиальная особая точка является устойчивым узлом; и нижняя область (R—+), где 7-1 > 0 ,
1 52
BE, 2005, 1 4
X2 <0, Xз <0, Xі єR, что удовлетворяет случаю Б, т. е. в этой области нетривиальная особая точка является неустойчивой.
Таким образом, на одном численном примере показано, насколько различным может быть поведение системы, описывающей производство смешанной культуры, которая состоит из двух видов популяций, в проточной системе при взаимодействии типа комменсализм при изменении некоторых характеристик.
5. Выводы
Научная новизна. Впервые для общей системы нелинейных дифференциальных уравнений, описывающей производство смешанной культуры, которая состоит из двух видов популяций, с видом взаимодействия типа комменсализм, в проточной системе, а именно в хемостате, выписаны в общем виде характеристические определители для найденных четырех особых точек. На численном примере показана возможность различных исходов производства описанной смешанной культуры при вариации скорости протока D и концентрации субстрата S0.
Практическая значимость. Полученные результаты дают возможность управлять процессом выращивания смешанной культуры, создавать условия для совместного сосуществования в субстрате двух культур в производственных целях с помощью вариации некоторых параметров. Исследование автономной системы на устойчивость при изменении скорости протока или концентрации питательного субстрата позволяет оптимально управлять процессом производства смешанных культур.
Сравнение с аналогами. Проблема сосуществования двух конкурирующих популяций была впервые поднята Г аузе. В соответствии с его принципом конкурентного исключения, в каждой экологической нише выживает только один вид. Однако в дальнейшем
понятие экологической ниши исследователи стали толковать по-разному, а по поводу самого принципа нередко высказывали крайне противоречивые суждения. Не так давно к проблеме поддержания разнообразия бактерий, конкурирующих за один и тот же субстрат, с совершенно другой стороны подошел Т. Чарань из исследовательской группы по теоретической биологии при Венгерской академии наук и его коллеги из Нидерландов - Р. Хукстра и Л. Паги. Предложенная ими математическая модель основывается на хорошо известной способности некоторых штаммов бактерий вырабатывать токсические вещества, угнетающие рост других штаммов. Такая модель показывает, что сосуществование нескольких видов, конкурирующих за один и тот же ресурс, теоретически возможно благодаря выработке конкурентами средств нападения, средств защиты. Предложенная в работе модель также показывает возможность сосуществования двух видов, конкурирующих за один субстрат.
Литература: 1. Gause G. F. The struggle for existence. Baltimore, 1934. 2. Бейли Дж., Оллис Д. Основы биохимической инженерии. Ч. 2. М.: Мир, 1989. 590 с. З.ГаузеГ. Г. Борьба за существование. 1999. //www.ggause.com/. 4. Йосс Ж., Джозеф Д. Элементарная теория устойчивости и бифуркаций. М.: Мир, 1983. 304 с. 5.Перт С. Дж. Основы культивирования микроорганизмов и клеток. М.: Мир, 1978. 334 с. 6. Романовский Ю. М., СтепановаН. В., Чернавский Д. С. Математическое моделирование в биофизике. М.: Наука, 1975. 344 с. 7.Романовский Ю. М., Степанова Н. В., Чернавский Д. С. Математическая биофизика. М.: Наука, 1984.
Поступила в редколлегию 27.11.2005
Рецензент: д-р тех. наук, проф. Удовенко С. Г.
Яловега Ирина Георгиевна, аспирантка кафедры ПМ ХНУРЭ. Научные интересы: математическое моделирование в биофизике и экологии. Адрес: Украина, 61000, Харьков, ул. Новгородская, 20, кв. 14, тел. 720-44-41.
BE, 2005, 1 4
1 5 3