Файл: Приц, А. К. Принцип стационарных состояний открытых систем и динамика популяций.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 30.10.2024
Просмотров: 46
Скачиваний: 0
rjl |
' ^ |
‘ dr,* |
|
a=i |
|
-О » J J I (oi-O i)Siil {/i7/—/i/j}duidQ= 0 ( 1 - 8 - 3 )
имеет место:
dH
;0, если ( 1 - 8 - 4 )
dt
dH
=0, если ( 1 - 8 - 5 )
dt
dH dt
функциях распределения f явно зависящих от времени. Стацио- dН
нарные решения (1—8—3) возможны при |
0- |
Легко показать, что единственным стационарным решением является распределение Максвелла, ибо определяющие аргу
менты фуНКЦИЙ |
Vi', t), |
V/ t), fi{rU Vi, t), |
fj{rh |
Vj, |
t) |
|
связаны между собой уравнениями: |
|
|
|
|||
|
|
al/+ 'o / = U i+ vjt |
( 1- |
8- |
6) |
|
|
|
|
|
|||
|
(o/)2+ ,(o/)2= (»i)2+ ( » i) 2- |
|
|
|
||
Распределение Максвелла является частным случаем распре |
||||||
деления Гиббса. Из сравнения (1—8— 1) и (1—3—9) |
видно, что |
|||||
|
|
Sm= - H |
( 1 - 8 - 7 ) |
|||
при соответствующем подборе значения х. Поэтому |
(1—8—4) |
и |
||||
(1—в—5) можно записать в виде: |
|
|
|
|||
|
d5 |
> 0 , |
при |
(1—8—8) |
||
|
d^ |
|
|
|
|
|
|
dS |
= 0, |
при f / f / ^ f i f j . |
( 1 - 8 - 9 ) |
||
|
dt |
|
|
|
|
|
§9. Стационарные состояния, модель Вольтерра
иколебательный характер численности популяций
Уравнение Вольтерра является одним из наиболее употреби тельных и наиболее изученных уравнений при описании динами ки популяций, например, [25, 40, 41, 42].
25
Рассмотрим |
систему, состоящую из двух компонент. ;В каче |
||
стве |
функции, |
характеризующей первую компоненту, возьмем |
|
численность |
рыбной популяции определенного вида, обитаю |
||
щей в определенном ареале. |
В качестве функции, характеризу |
||
ющей |
вторую |
компоненту, |
численность N2 рыбодобывающего |
флота, базирующегося на данный ареал.
Механизм взаимного влияния этих двух компонент в ареале обитания компоненты Ny может быть описан с помощью уравне
ний Вольтерра |
|
|
|
dNj_ |
aiNi— |
( 1 - 9 - 1 ) |
|
dt |
|||
|
|
||
dN2 |
0-2^2“Ь^2^1Л^21 |
(1—9—2) |
|
dt |
|||
|
|
где KiNiNz — есть скорость убывания численности первой компо ненты при ее взаимодействии со второй, т. е. убывание числен ности рыбной популяции при встрече рыбных стай с рыболовны ми судами, а сцМ — скорость прироста численности Ni; %2NiN2— дает скорость увеличения численности рыбодобывающего флота, д2Л^2— скорость убывания этой численности.
Выражение (1—9—2) показывает, что численность флота увеличивается при успешном лове, т. е. при достаточном рыбном запасе, который и обуславливает частоту встреч косяков рыб с рыболовными судами.
Уравнения (1—9— 1) и (1—9—2), отражающие связь запаса
рыбы (первая компонента) |
и численность |
рыбодобывающего |
флота (вторая компонента), |
могут быть записаны в виде: |
|
dNi |
V i |
(1—9—3) |
— ±- =/г. ^ + рг 1 2 , aijNiNj . |
j=i
Первый член суммы представляет поведение i-ой компоненты в отсутствии второй. При k C > 0, Nj= ,0 происходит экспоненци альный рост численности /-ой компоненты.
При йг<0 наблюдается экспоненциальное уменьшение чис ленности.
Постоянная ац может быть положительной, отрицательной или равной нулю. Если a ij> 0 , то это означает насколько часто та встреч судна с косяками рыб ведет к увеличению числа судов, отрицательное же значение aij отвечает на вопрос о влиянии ча стоты встречи рыбы и рыболовного судна на уменьшение чис ленности рыб. Если же ciij равна нулю, то это соответствует от сутствию взаимосвязи между числом судов и численностью ры бы (суда не ловят рыбу).
26
Таким образом, в результате взаимодействия t-ой и /-ой ком понент происходит увеличение численности t-ой компоненты и. уменьшение /-ой. Поэтому ciij и ац имеют противоположные
знаки.
В результате «встреч» t-ой и /-ой компонент произведение р£—ipj дает уменьшение (или увеличение) численности t-ой ком поненты в единицу времени.
Исключая особые случаи, будем иметь:
|
(1—9—4) |
ац — 0. |
(1—9—5) |
■Стационарное состояние системы определяется совокупно
стью значений ЛД для которых |
|
~~77~ — 0 |
( 1 - 9 - 6 ) |
сп |
|
при всех значениях /.
Те же значения ЛД которыеудовлетворяют уравнению
(1—9—6), обозначим через {gy}. Отсюда |
|
|
gi 1&Р*+ 2 |
= 0 - |
(1 —9—7) |
3 |
|
|
Если gi=£0, то получаем |
|
|
k $ i + 2 |
аб‘Я з=°- |
(1—9—8) |
3 |
|
|
Уравнение i(l—9—8) было впервые исследовано Д. Анкона [40] при введении интегралов движения в динамику развития попу ляций.
Введем обозначения
о ,= 1 п - ^ - . N j = g j e x p (vj) • |
( 1 - 9 - 9 ) |
gj |
|
Совершенно очевидно, что Vj стремится к нулю при стремлении Nj к gj и может служить характеристикой отклонения значений Nj от своего стационарного значения gj.
Подставляя в (1—9—3) вместо Nj его значение из (1—9—9), получим:
Р ; ^ = = 2 ^ М ^ - 1 ) . |
(1 - 9 - Ю ) |
i |
|
Умножив обе части уравнения '(1—9— 10) на gy[e®i— 1] и про суммировав по всем /, получим
27
- ^ |
2 |
Pi&te,,,— °*1=! |
|
|
i |
|
|
= 2 ^ |
^ |
. - D (ev3—1) = 0 , |
(1—9— 11) |
ji |
|
|
|
ибо двойная сумма антисимметрична относительно i и /. 'Находим
|
|
G= 2 |
6 } — У ,! Piffle”;—Щ-] = const ■ |
(1—9— 12) |
|
|
3 |
31 |
|
где G — тот |
интеграл движения, который надо было найти. |
|||
Значения Gj больше нуля, т. к. |
|
|||
а) |
V j>0 |
и evi>-Vj, |
|
|
б) |
vjCO |
и —ол-> 0 , экспонента же evj будет положительной, |
||
независимо |
от знака vj. |
|
||
Поэтому |
G > 0 . |
Ввиду того, что ' 2 Big,- есть |
величина по |
|
стоянная, можно также найти, что |
|
iY
Ga= — “ Р ^ (1 —V— evs) — const. |
(1—9— 13) |
Воспользуемся константой движения G для выяснения перио дичности состояний системы, состоящей из двух взаимодейству
ющих компонент в соответствии |
с (1—9— 1) и |
(1—9—2). |
||||||
|
Уравнения |
(1—9— 1) |
и |
(1—9—2) |
эквивалентны уравнениям |
|||
(1—9—3) если |
|
|
|
|
|
|||
|
|
ai=& i, |
Ol= |
—ki, |
(1—9— 14) |
|||
или |
A.lPi= —#12, |
taP2= 021 |
(1—9— 15) |
|||||
(3i= Х2Р2 =—— |
$21- |
(1—9— 16) |
||||||
г> |
И" |
|||||||
|
|
|
|
|
|
|||
Отсюда, используя определение gi, находим: |
|
|||||||
|
|
„ |
“2 |
|
ai |
(1—9— 17) |
||
|
|
^ |
*2 |
’ |
|
Я, ' |
||
|
|
|
|
|||||
Вводя функции |
|
|
|
|
|
|||
|
|
f |
а д |
’ |
f |
N2(t) |
(1—9— 18) |
|
|
|
/1= |
гг |
/2 = |
S2 |
|||
|
|
|
Si |
|
|
|
||
найдем,что |
|
|
|
|
|
( 1 - 9 - 1 9 ) |
||
|
|
|
|
|
|
|
||
Подставляя |
значения |
(1—9— 19), |
(1—9— 17), |
(1—9— 16) и |
||||
(1—9— 15) в (1—9— 12) получим: |
|
|
28
1 |
i |
|
{/че_А} “i {fie-ty аг = const. |
(1—9—20) |
|
Введя обозначения |
* |
|
i |
|
|
Ui= { f le~h}'^, |
U2= {f2e - h } ^ 2 , |
( 1 - 9 —21) |
получим гиперболу |
|
(1—9—22) |
UiUz= const, . |
определяющую набор стационарных состояний в пространстве переменных Ux и Uz (рис. 2а).
Рис. 2. График периодичности решения уравнений (1—9—3)
Рисунки 2b и 2с описывают зависимость Ut и Uz от ft и f2. Осо бенностью фигур 2b и 2с является то, что они имеют максимумы в точках М2 и Mi. Поэтому гипербола имеет границы в точках Л и В. Одному значению точки 0, лежащей между Л и В, .соот ветствуют два значения f\{al) и fi{b ]) и два значения f2(c) и fz(b). На рисунке 2d устанавливается связь между функциями fi, f2При перемещении точки 0 от Л к В траектория этой точки представится в виде окружности «а рис. 2d. Значения Л и В соответствуют двум экстремальным значениям функции f i и f%со ответственно.
29