Файл: Приц, А. К. Принцип стационарных состояний открытых систем и динамика популяций.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 30.10.2024

Просмотров: 46

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
< 0 только при нестационарных решениях (1—8—3),т. е. при

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)

эквивалентны уравнениям

(19—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