Научная статья на тему 'Эколого-экономическая система эксплуатации биоресурсов с векторными платежами'

Эколого-экономическая система эксплуатации биоресурсов с векторными платежами Текст научной статьи по специальности «Математика»

CC BY
95
12
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЗАДАЧА УПРАВЛЕНИЯ БИОРЕСУРСАМИ / BIORESOURCE MANAGEMENT PROBLEM / МНОГОКРИТЕРИАЛЬНЫЕ ИГРЫ / MULTICRITERIA GAMES / ДИНАМИЧЕСКИЕ ИГРЫ / DYNAMIC GAMES / АРБИТРАЖНОЕ РЕШЕНИЕ НЭША / NASH BARGAINING SOLUTION

Аннотация научной статьи по математике, автор научной работы — Реттиева Анна Николаевна

Исследована теоретико-игровая модель эколого-экономической системы в дискретном времени. В игре участвуют агенты (фирмы или рыболовецкие артели), производящие вылов биоресурсов на конечном промежутке времени. Агенты эколого-экономической системы хотят добиться нескольких целей, поэтому используются векторные функции выигрышей. Предложены концепции определения оптимального некооперативного и кооперативного поведения в многокритериальной динамической игре. Для построения некооперативного равновесия использована конструкция арбитражной схемы Нэша (произведения Нэша), а для определения кооперативного арбитражная схема Нэша для всего периода продолжения игры. Целью работы является определение кооперативного поведения участников и применение разработанной схемы для рационального использования биоресурсов водоемов Республики Карелия.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ECOLOGY-ECONOMIC SYSTEM OF BIORESOURCE EXP-LOITATION WITH VECTOR PAYOFFS

A discrete time game-theoretic model of an ecology-economic system is considered. Agents (firms or fishermen’s artels) that exploit the fish stock on a finite planning horizon are the participants of the game. Agents wish to achieve several goals, hence the vector functions for players’ payoffs are considered. The approaches to determining the optimal noncooperative and cooperative behavior in dynamic multicriteria games are constructed. The multicriteria Nash equilibrium is obtained via the Nash bargaining design (Nash products), and the cooperative equilibrium is determined by the Nash bargaining procedure for the entire planning horizon. The main goal of the paper is to define the cooperative behavior of the participants and to apply the procedure for sane exploitation of Karelian bioresources.

Текст научной работы на тему «Эколого-экономическая система эксплуатации биоресурсов с векторными платежами»

Труды Карельского научного центра РАН № 7. 2018. С. 83-91 DOI: 10.17076/mat847

УДК 519.83

ЭКОЛОГО-ЭКОНОМИЧЕСКАЯ СИСТЕМА ЭКСПЛУАТАЦИИ БИОРЕСУРСОВ С ВЕКТОРНЫМИ ПЛАТЕЖАМИ

А. Н. Реттиева

Институт прикладных математических исследований КарНЦ РАН, ФИЦ «Карельский научный центр РАН», Петрозаводск, Россия

Исследована теоретико-игровая модель эколого-экономической системы в дискретном времени. В игре участвуют агенты (фирмы или рыболовецкие артели), производящие вылов биоресурсов на конечном промежутке времени. Агенты эколого-экономической системы хотят добиться нескольких целей, поэтому используются векторные функции выигрышей. Предложены концепции определения оптимального некооперативного и кооперативного поведения в многокритериальной динамической игре. Для построения некооперативного равновесия использована конструкция арбитражной схемы Нэша (произведения Нэша), а для определения кооперативного — арбитражная схема Нэша для всего периода продолжения игры. Целью работы является определение кооперативного поведения участников и применение разработанной схемы для рационального использования биоресурсов водоемов Республики Карелия.

Ключевые слова: задача управления биоресурсами; многокритериальные игры; динамические игры; арбитражное решение Нэша.

A. N. Rettieva. ECOLOGY-ECONOMIC SYSTEM OF BIORESOURCE EXPLOITATION WITH VECTOR PAYOFFS

A discrete time game-theoretic model of an ecology-economic system is considered. Agents (firms or fishermen's artels) that exploit the fish stock on a finite planning horizon are the participants of the game. Agents wish to achieve several goals, hence the vector functions for players' payoffs are considered. The approaches to determining the optimal noncooperative and cooperative behavior in dynamic multicriteria games are constructed. The multicriteria Nash equilibrium is obtained via the Nash bargaining design (Nash products), and the cooperative equilibrium is determined by the Nash bargaining procedure for the entire planning horizon. The main goal of the paper is to define the cooperative behavior of the participants and to apply the procedure for sane exploitation of Karelian bioresources.

Keywords: bioresource management problem; multicriteria games; dynamic games; Nash bargaining solution.

Введение

Задачи рационального природопользования являются актуальными для Республики Карелия, обладающей большими запасами возобновляемых ресурсов. Одной из таких за-

дач является определение оптимального поведения в случае несимметричности агентов эколого-экономической системы, т. к. возобновляемые ресурсы Республики Карелия подвергаются совместной эксплуатации различными экономическими субъектами.

®

Математические модели, учитывающие наличие нескольких целевых функций у участников конфликтно-управляемых процессов, более приближены к реальности. Зачастую игроки хотят достичь одновременно нескольких целей, которые могут быть несравнимы. Такая ситуация типична для водоемов Республики Карелия, где хозяйствующие субъекты несут различные затраты на эксплуатацию, что влечет за собой различие в ценах для конечного потребителя. При этом участники хотят одновременно и увеличить свою прибыль от продажи эксплуатируемого ресурса, и уменьшить издержки. Такая постановка задачи влечет введение вектор-функций выигрышей участников процесса эксплуатации ресурса и исследование многокритериальных игр.

Ллойд Шепли [7] в 1959 г. ввел понятие многокритериальной игры, т. е. игры с векторными функциями выигрышей участников, и расширил для таких игр концепцию равновесия по Нэшу, получив, таким образом, оптимальность по Парето (сильную и слабую). В последние годы много работ посвящено играм с векторными выигрышами и концепциям их решения. Например, в [8] предложено понятие идеального равновесия по Нэшу, а в [5] предложена концепция Е-равновесия. Классическим способом решения векторных игр является их скаляризация путем оптимизации взвешенной суммы критериев, который неприменим для несравнимых критериев. Тем не менее понятие равновесия по Парето является самой распространенной концепцией решения таких игр. В кооперативных многокритериальных играх для распределения общего кооперативного выигрыша используется естественное расширение вектора Шепли.

Однако все предложенные концепции решений используются только в статических многокритериальных играх. Малоисследованной проблемой является построение равновесий в динамических многокритериальных играх. В работе [6] было формализовано понятие многокритериального равновесия по Нэшу с использованием конструкции арбитражной схемы Нэша и предложены три варианта построения гарантированных выигрышей.

Главной целью представленной работы является исследование кооперативной динамической многокритериальной игры и формализация оптимального решения. Необходимость исследования кооперации связана с тем, что кооперативное поведение благоприятно влияет на эксплуатируемый ресурс. Для определения кооперативного поведения используется разработанный ранее подход построения ко-

оперативного равновесия в теоретико-игровых моделях с несимметричными участниками [2]. А именно - кооперативные стратегии и выигрыши участников определяются из решения арбитражной схемы Нэша для всего периода продолжения игры. При этом точками статус-кво выступают некооперативные выигрыши, полученные при использовании игроками многокритериальных равновесных по Нэ-шу стратегий. Предложенная схема позволяет спрогнозировать поведение агентов эколого-экономической системы при заключении между ними кооперативного договора. При этом многокритериальный подход позволяет учесть специфику хозяйствующих субъектов и их влияние на эксплуатируемый возобновляемый ресурс.

Предложенные концепции решения апробированы для популяции сига озера Сямозе-ро. По имеющимся данным матричным методом оценки запаса была восстановлена численность популяции и параметры функции развития [3]. Построены некооперативное и кооперативное равновесия. Проведено сравнение стратегий игроков и размера эксплуатируемого ресурса для различных вариантов построения гарантированных выигрышей. Показано, что вариант определения гарантированных выигрышей с помощью равновесных по Нэшу решений является наилучшим для экологической ситуации.

Проведено сравнение состояния экологической системы и прибыли агентов при кооперативном и эгоистическом поведении. Численное моделирование показало, что использование арбитражной схемы для определения кооперативного поведения не только выгодно агентам эколого-экономической системы, что позволит мотивировать хозяйствующих субъектов заключать кооперативные договоры, но и, что самое главное, благотворно влияет на экологическую систему.

ДИНАМИЧЕСКИЕ МНОГОКРИТЕРИАЛЬНЫЕ ИГРЫ

Формализуем понятие динамической многокритериальной игры. Рассмотрим бикрите-риальную динамическую игру с двумя игроками в дискретном времени. Игроки эксплуатируют некоторый общий возобновляемый ресурс и стремятся достигнуть двух различных целей. Динамика развития возобновляемого ресурса имеет вид

Хг+1 = /(хг,Пц,П2г) , Хо = х, (1)

где х^ ^ 0 — размер ресурса в момент времени £ ^ 0, /{х1,иц,П2г) — функция развития

©

возобновляемого ресурса, пц € и — стратегия (интенсивность эксплуатации) игрока г в момент времени Ь ^ 0, г = 1, 2.

Вектор-функции выигрышей игроков на бесконечном промежутке планирования имеют вид

Т =

Т =

/

/те \

Л1 = Е 5*5 !(п 1 г ,П24)

4=0 те

Л2 = Е ¿г£2(и 1 г,П2г)

\ 4=0

/те \

Т = Е (п 1 4, «24) 4=0 те

Л2 =Е 5^ («14, «24) \ 4=0

(2)

/

где 5(п1г,п2г) ^ 0 — функции «мгновенного» выигрыша, г, ^ = 1, 2, 5 € (0,1) — общий коэффициент дисконтирования.

статус-кво:

Я1(щг,П2г) = = (Л1(п1г,п2г) - С1)(Л2(п1г,П2г) - С?),

Н2(п1г,п2г) = = (^(«14, П2г) - С2)(Т|(«14, «24) - О2).

Следующее определение дает предложенную концепцию некооперативного решения динамической многокритериальной игры [6].

Определение 1. Профиль стратегий К ,«") является многокритериальным равновесием по Нэшу в игре (1), (2), если

Н1(пи) ^ Я1(п14,п^) € и

Многокритериальное равновесие по Нэшу

Для построения некооперативного равновесия в многокритериальной динамической игре используется подход, предложенный в [6]. А именно — применяется конструкция арбитражной схемы Нэша (произведения Нэша). Поэтому сначала определяются гарантированные выигрыши, которые играют роль точек статус-кво. Предложены три варианта построения гарантированных выигрышей.

В первом из них все гарантированные выигрыши определяются из решений антагонистических игр. В частности, первый гарантированный выигрыш С1 — это решение антагонистической игры, в которой первый игрок стремится максимизировать свой первый критерий, а второй — минимизировать его. Остальные гарантированные точки О1, С1, О^ строятся аналогично.

Во втором способе гарантированные точки определяются из решений антагонистических игр с суммами критериев. То есть оба гарантированных выигрыша первого игрока С11 и С21 — это решение антагонистической игры, в которой первый игрок стремится максимизировать сумму своих критериев, а второй — минимизировать его. Аналогично и для второго игрока.

В третьем варианте гарантированные выигрыши строятся как равновесия по Нэшу в играх с первыми и вторыми критериями игроков соответственно.

Для построения выигрышей игроков в динамической многокритериальной игре используем произведения Нэша, при этом гарантированные выигрыши играют роль точек

#2«,«£) ^ #2«,П2г) УП24 € и .

Следовательно, как и в классическом равновесии по Нэшу, игрокам невыгодно отклоняться от оптимальной стратегии. При этом они максимизируют произведение отклонений оптимальных выигрышей от гарантированных.

кооперативное

Многокритериальное равновесие

Теперь предположим, что игроки хотят действовать кооперативно. Для определения кооперативного поведения используется арбитражная схема Нэша для всего периода продолжения игры. При этом в качестве точки статус-кво выступают некооперативные выигрыши, полученные при использовании игроками многокритериальных равновесных по Нэ-шу стратегий. Первая точка — это сумма некооперативных выигрышей по первому критерию, вторая — сумма по второму критерию обоих участников соответственно.

Выигрыши игроков в многокритериальном равновесии по Нэшу (см. определение 1) примут вид

Т ^

Т2

/ = Е 5г£1С< ,п£

4=0 те

тГ = Е 5г£2« ,«N4)

\ г=0 /

Т = Е 5г51(п^г ,«N4) х

4=0 те

ТГ = Е ^(«й ,«£)

\ 4=0

(3)

/

Для определения кооперативного поведения используется арбитражная схема Нэша, следовательно, необходимо решить следую-

щую задачу:

(V/ c + c - J1 N - n ) ■

(y2c + V22c - J2N - J2N) =

E^E gl (ui t,«2t) - J1 N - J21 N t=0 i=1

те 2 :t

Гарантированные выигрыши и многокритериальное равновесие по Нэшу

Первый вариант

Используя результаты [6], заключаем, что гарантированные выигрыши для первых критериев обоих игроков равны нулю: С} = С2 = 0, а для вторых критериев принимают вид

те 2

■ (Е * Е gi2(u1t, u2t) - J2N - J|N) ^ mmax , (4

V n . ' U1t,U2t

t=0 i=1

TjN

xt+1 = ext - U1t - U2t, xo = x,

(5)

J =

J2 =

/те \

J1 = E 5*Р1«14 t=0

те

J2 = - E РоЩ^ \ t=0

/

J1 = E

\

t=0 те

J| = - ^ 5*С«14U2t \ t=0

(6)

/

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

G2 — G-i — G — — c—--— X0 21 45 0

где — некооперативные выигрыши, определенные в (3), г,^ = 1,2, г = j.

Определение 2. Профиль стратегий (и}4,и24) является многокритериальным кооперативным равновесием в игре (1), (2), если является решением задачи (4).

Многокритериальная задача управления возобновляемыми ресурсами

Исследуется динамическая многокритериальная модель управления биоресурсами. Два игрока (фирмы или рыболовецкие артели) эксплуатируют ресурс на протяжении бесконечного промежутка времени. Динамика развития популяции имеет вид

Вследствие того, что первые гарантированные выигрыши равны нулю, для определения многокритериального равновесия по Нэшу в задаче (5), (6) необходимо решить следующую задачу:

J]L(u1t,U2t)(J12(u1t,U2t) - G) ^ max,

«it

J21(u1t,u2i)(j|(u1i,u2i) - G) ^ max .

Используя принцип Беллмана и строя оптимальные стратегии в линейном виде, получим многокритериальное равновесие по Нэшу

uit = uN

(fe2 - i)(g - 1)

г52 - 1 + 5г(г - 1)

xt,

(7)

где Х£ ^ 0 — размер популяции в момент времени £ ^ 0, е ^ 1 — коэффициент естественного роста, иц ^ 0 — стратегия (вылов) игрока г в момент времени £ ^ 0, г = 1, 2.

Игроки стремятся достичь двух целей — максимизировать доход от продажи ресурса и минимизировать затраты на эксплуатацию. Предполагается, что цена на рынке для участников различна, а затраты равны и зависят от интенсивности эксплуатации обоих игроков. Тогда вектор-функции выигрышей на бесконечном промежутке планирования примут вид

и размер ресурса при некооперативном поведении примет вид

xt

г2 5 + г - 2

U52 - 1 + 5г(г - 1)J

t

x0

(8

Второй вариант

Используя результаты [6], заключаем, что гарантированные выигрыши принимают вид

р,(5е2 - 1) х р2(5е2 - 1) С, — "ггг;-тг- Хо

25(г - 1) 2с5(г - 1)2 с(5г2 - 1) 2 , р2(5г2 - 1) ,

G2 —--

45

x0 +

4с5(г - 1)2

i = 1,2 . (9)

Для определения многокритериального равновесия по Нэшу в задаче (5), (6) необходимо решить задачу

где р, ^ 0 — рыночная цена за единицу ресурса для игрока г, г = 1,2 с ^ 0 — затраты на эксплуатацию и 5 € (0,1) — общий коэффициент дисконтирования. Предполагается, что параметры таковы, что 5е ^ 1.

Е - G1 j Е 5*си14и24 - G2j ^ max,

t=0 t=0 Ult те те

Е 5*Р2и24 - G2) (- Е 5*си14и24 - G^ ^ max

t=0

t=0

Аналогично, используя принцип Беллмана и строя оптимальные стратегии в линейном

виде, получим многокритериальное равновесие по Нэшу

uNt =

(¿e2 - 1)(е -1}

xt +

+

¿e2 - 1 + ¿e(e - 1) ¿£P2G}(£ - 1) - pi(¿e2 - 1) 2pip2(¿e2 - 1 + ¿e(e - 1))

(¿e2 - 1)(e - 1)

uNt

Xt +

+

¿e2 - 1 + ¿e(e - 1) ¿epiGl(e - 1) - Glp2(¿e2 - 1) 2pip2(¿e2 - 1 + ¿e(e - 1))

(10)

а размер ресурса при некооперативном поведении принимает вид

Xt =

¿e2 + e - 2

L¿e2 - 1 + ¿e(e - 1)J

í P2GI + piG2 xo + i 2

Gi _ Pi x Gi _ P2 x

Gi — — Xo , G2 — — Xo . ¿¿

P2

(12)

G2 — G2 — c¿e2 1 x2 Gi — g2 — -c 4¿ x0

(13)

Утверждение 1. Если x0 < — {Pi,P2}, mo

c(e - 1)

xvav2 - xvavi - xvav3

xxx,

t

uVtar2 - uvtavi - uvtav3 , i — 1, 2.

t

va-^ít

—ax(pi,p2>

Если x0 > ----—, то

c(e - 1)

^vavi ^ rpvav2 ^ /y)vav3

uvavi - uvav2 - uv

t

vav3 ít

, i — 1, 2.

2р1Р2(е - 1) ' (11)

Третий вариант

В третьем варианте построения гарантированных выигрышей С11 и С12 являются равновесием по Нэшу в игре (1,11, и, и2, Т|, Т1) и принимают вид

Аналогично, определяя равновесие по Нэ-шу в игре со вторыми критериями ^(«^,«24) и Т|(«14, «24), получены вторые гарантированные выигрыши

Многокритериальное равновесие по Нэшу в данном случае имеет такой же вид, как и во втором варианте (10),(11), но с соответствующими гарантированными выигрышами (12),(13).

Сравнение некооперативных решений

Сравним стратегии игроков и размер популяции для различных вариантов построения гарантированных выигрышей. Оказывается, что существуют два случая в зависимости от параметров модели: в первом их них наихудшим для экологии вариантом является второй (когда гарантированные выигрыши определяются из решения антагонистической игры с суммами критериев), во втором — первый (когда гарантированные выигрыши определяются из решения антагонистических игр). При этом в обоих случаях наилучшим для состояния эксплуатируемой системы и выгодным для участников является третий вариант построения гарантированных выигрышей (равновесные по Нэшу решения).

Доказательство. Сравним размер популяции для всех вариантов построения гарантированных выигрышей. Обозначим х^1 =

е25 + е - 2 и*

-:—;-т х0 размер популяции при

е52 - 1 + 5е(е - 1) ] использовании первого варианта построения

гарантированных выигрышей (8). Тогда во втором и третьем вариантах (11) размер популяции примет вид

xt

xvavi +

P2^i + G2 Pi 2pip2(e - 1)

с соответствующими гарантированными выигрышами, определенными в (9) и (12).

Для третьего варианта рассмотрим числитель дроби 1), получим

Р201 + О1Р1 = —-— Х0 > 0, 5

следовательно, xvav3

> xvav i

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

, г = 1,2.

Для второго варианта также рассмотрим числитель, получим

П1 . П1 Р1Р2(5е2 - 1) / Р1 + Р2 \

Р2О1 + О1Р1 = 2*(£-1) 12Х0 - ¿ГТ)),

следовательно, при x0 ^

Pi + P2 2m(e - 1)

x

vav2

xvavi, а при x0 - Pi +_P2 xvav2 - xrvi.

2C(e 1)

Сравним стратегии игроков для всех вариантов построения гарантированных выигрышей. Обозначим uvavi — uvavi —

(¿e2 - 1)(e - 1) vavi

—;-—--xv i некооперативные

e¿2 - 1 + ¿e(e - 1) 1

стратегии при использовании первого варианта построения гарантированных выигрышей (7). Тогда во втором и третьем вариантах (10) стратегии примут вид

uit

U2t

_vavi + ui t +

(¿e - 1)Gi

2pi(¿e2 - 1 + ¿e(e - 1)):

u2avi +

(¿e - 1)G2

2p2(¿e2 - 1 + ¿e(e - 1))'

с соответствующими гарантированными выигрышами, определенными в (9) и (12).

Для третьего варианта из (12) заметим, что

(¿£-1)01

числители выражений 2р(¿£2-1+ге(£-1)), г 1, 2, положительны, следовательно, и

или

-иагЗ

>

и

,4

, г = 1, 2.

Для второго варианта построения гарантированных выигрышей аналогичные числите-

Р1

ли зависят от знака выражений хо

и х0--—Р2—-, следовательно, и^4аг2 ^ и™г1

с(е — 1)

с(е — 1)

а

г = 1, 2, при х0 ^ -^—и и|4аг2 ^ и™г1

при х0 ^

шт{р1,р2}

с(е — 1)

и

заметим, что при этом выпол-

с(е — 1)

Теперь рассмотрим случай х0 ^

шах{р1,р2}

с(е — 1) ,

^ Р1 + Р2 п

нено х0 ^ —---. Сравним стратегии и раз-

2с(е — 1)

мер популяции во втором и третьем вариантах построения гарантированных выигрышей. Заметим, что при выполнении данных условий гарантированные выигрыши в третьем варианте (12) больше, чем во втором (9). Следовательно, <4аг3 ^ <4аг2, г = 1, 2, и хГг3 ^ х£аг2

откуда и получаем утверждение.

Многокритериальное кооперативное равновесие

Используя первый вариант построения гарантированных выигрышей, определим выигрыши в многокритериальном равновесии по Нэшу

3

12

4=0

оо

5е2 — 1 5е — 1

х,

тШ & N 5е — 1

32 = 2_^ 5р2ии = р2——г х,

т22 31

4=0

т22 32

5е — 1

-сУ^54

4=0

"и14 и24

= _ с (5е2 — 1)2(е — 1) х2

С (5е — 1)2(3 — е) х '

Для построения многокритериального кооперативного равновесия необходимо решить следующую задачу:

+ ^21с — 312 — 3212)

(V2- + VIе — 322 — 3222).

шах

+ +

(14)

4=0

Е 54(р1и14 + Р2иу — ■

4=0 те

■ (—2с V 54и14и24 — ^х2) ^ шах , (15)

(5е2 — 1)(р1 + р2)

где С2 = —2с

С1 =

(5е2 - 1)2(е - 1)

5е 1

(5е — 1)2(3 — е) ' Используя принцип Беллмана и строя стратегии как и ранее в линейном виде, получим следующее утверждение.

Утверждение 2. Кооперативные стратегии в многокритериальной игре (55), (6) имеют вид

14

24

2сС1 + 125еА + К 6с(р1 + р2) + 245А

х4 ,

(16)

где

К = (2сС1(сС1 + 125еА) —

— 12(сС2р2 + 5Ар(3се2 + 2^)))2 ,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

а параметр А является решением следующего уравнения:

К3

27А(с(р1 + р2) + 45А)2 = — +

+ 45АС1С2(1085А + 45с(р1 + Р2)) +

+ 18с2С1(4е5АС1 + С2(р1 + Р2)2) —

— 27е5А(с(р1 + Р2) — 45А) ■

■ ((се2 + 3)(р1 + р2) + 3сеС1).

Динамика развития ресурса при кооперативном поведении принимает вид

"3се(р1 + Р2) — 2сС1 — К14

х4 =

3с(р1 + р2) + 125А

х0.

(17)

Моделирование популяции сига озера Сямозеро

Для моделирования использовались фактические данные о многотычинковом сиге озера Сямозеро [4]. Для восстановления численности сига был использован матричный метод оценки запаса [1].

Полученные оценки размера популяции за длительный период позволили оценить параметры функции ее развития. Наиболее адекватной ситуации оказалась степенная функция развития [3]. В данной работе будем использовать линейную аппроксимацию функции развития популяции с параметром е = 1, 3.

Численное моделирование было проведено со следующими параметрами:

Хо = 100000 , S = 0,8 , c = 50 , pi = 100 , p2 = 150 .

На рисунке 3 представлена динамика развития популяции, а на рисунке 4 — вылов первого игрока при кооперативном и некооперативном поведении. Заметим, что при кооперации экологическая обстановка лучше, и при этом интенсивность эксплуатации игрока выше.

На рисунке 1 показана динамика развития популяции, на рисунке 2 — вылов первого игрока для различных вариантов построения гарантированных выигрышей. Заметим, что, как это и было получено аналитически, наихудшим для экологической обстановки вариантом является первый, он ведет к перелову популяции. А самым лучшим вариантом построения гарантированных выигрышей является третий, когда гарантированные выигрыши — это равновесие по Нэшу, здесь и размер популяции, и интенсивности эксплуатации участников больше.

x(t)

Рис. 1. Размер популяции: темная линия -1 вариант, пунктир - 2 вариант, светлая линия -3 вариант

Fig. 1. Size of population: dark line shows 1st variant, dotted line - 2nd variant, light line - 3rd variant

u1(t)

40000 35000 30000 25000 20000 15000 10000

20 30 _Time t_

Рис. 2. Вылов первого игрока: темная линия -1 вариант, пунктир - 2 вариант, светлая линия -3 вариант

Fig. 2. Catch of player 1: dark line shows 1st variant, dotted line - 2nd variant, light line - 3rd variant

x(t)

Рис. 3. Размер популяции: темная линия - кооперация, светлая - равновесие по Нэшу Fig. 3. Size of population: dark line shows cooperation, light line - Nash equilibrium

u1(t)

Рис. 4. Вылов первого игрока: темная линия - кооперация, светлая - равновесие по Нэшу Fig. 4. Catch of player 1: dark line shows cooperation, light line - Nash equilibrium

Заключение

В работе исследованы многокритериальные динамические игры в дискретном времени и предложены подходы к определению равновесных стратегий. Для определения многокритериального некооперативного равновесия используются произведения Нэша, а для построения точек статус-кво предложены три варианта. Для иллюстрации предложенной концепции решения исследована динамическая многокритериальная модель управления эколого-экономической системой эксплуатации возобновляемого ресурса. Проведено

сравнение стратегий агентов и размера эксплуатируемого ресурса для различных вариантов построения гарантированных выигрышей. Представлены два случая в зависимости от параметров модели: в первом случае наихудшим для экологической ситуации вариантом является второй (когда гарантированные выигрыши определяются из решения антагонистической игры с суммами критериев), во втором случае наихудшим для экологии вариантом является первый (когда гарантированные выигрыши определяются из решения антагонистических игр).

Исследовано кооперативное поведение в многокритериальных динамических играх. Для определения кооперативных стратегий и выигрышей игроков была использована арбитражная схема Нэша для всего периода продолжения игры. При этом в качестве точки статус-кво выступает многокритериальное равновесие по Нэшу. Концепция решения применена для исследования кооперативного поведения в динамической многокритериальной модели управления эколого-экономической системой эксплуатации возобновляемого ресурса. Предложенная схема позволяет спрогнозировать поведение агентов при заключении между ними кооперативного договора. При этом многокритериальный подход позволяет учесть специфику хозяйствующих субъектов и их влияние на эксплуатируемый возобновляемый ресурс.

Для численного моделирования были использованы данные о популяции сига озера Сямозеро. Построены некооперативное и кооперативное равновесия. Показано, что применение арбитражной схемы для определения кооперативного поведения выгодно обоим агентам и при этом благоприятно влияет на состояние экологической системы. Исследованный многокритериальный подход позволяет спрогнозировать состояние эксплуатируемой популяции на длительный промежуток времени, что позволит управляющим структурам контролировать состояние экологических систем Республики Карелия и принимать ме-

References

1. Abakumov A. I., Kol'ev N. V., Maksimenko V. P., Gorr S. V. Matrichnii metod otsenki zapasa i prognozirovanie vylova populatsii morskikh organizmov [A matrix method for marine populations stock assessment and catching prediction]. Voprosy ikhtiologii [Journal of Ichthyology]. 1994. Vol. 34, no. 3. P. 400-407.

ры во избежание перелова и гибели популяции.

Финансовое обеспечение исследований осуществлялось из средств федерального бюджета на выполнение государственного задания КарНЦ РАН (Институт прикладных математических исследований КарНЦ РАН) и при финансовой поддержке РФФИ (проекты 16-01-00183_а, 16-41-100062 р_а).

Литература

1. Абакумов А. И., Кольев Н. В., Максимен-ко В. П., Горр С. В. Матричный метод оценки запаса и прогнозирования вылова популяций морских организмов // Вопросы ихтиологии. 1994. Т. 34, № 3. С. 400-407.

2. Реттиева А. Н. Задача управления биоресурсами с асимметричными игроками // Математическая теория игр и ее приложения. 2013. Т. 5, вып. 3. С. 72-87.

3. Реттиева А. Н. Эколого-экономическая система эксплуатации биоресурсов с асимметричными участниками // Труды КарНЦ РАН. Сер. Математическое моделирование и информационные технологии. 2016. № 8. С. 91-97. doi: 10.17076/mat441

4. Стерлигова О. П., Павлов В. Н., Иль-маст Н. В., Павловский С. А., Комулай-нен С. Ф., Кучко Я. А. Экосистема Сямозера (биологический режим, использование). Петрозаводск: КарНЦ РАН, 2002. 119 с.

5. Pusillo L., Tijs S. E-equilibria for multicriteria games // Annals of ISDG. 2013. Vol. 12. P. 217228.

6. Rettieva A. N. Multicriteria dynamic games // International Game Theory Review. 2017. Vol. 1(19). P. 1750002.

7. Shapley L. S. Equilibrium points in games with vector payoffs // Naval Research Logistic Quarterly. 1959. Vol. 6. P. 57-61.

8. Voorneveld M., Grahn S., Dufwenberg M. Ideal equilibria in noncooperative multicriteria games // Mathematical Methods of Operations Research. 2000. Vol. 52. P. 65-77.

Поступила в редакцию 23.04.2018

2. Rettieva A. N. Zadacha upravleniya bioresursami s asimmetrichnymi igrokami [Discrete-time bioresource management problem with asymmetric players]. Matematicheskaya teoriya igr i ee prilozheniya [Mathematical Game Theory and its Applications]. 2013. Vol. 5, iss. 3. P. 72-87.

3. Rettieva A. N. Ecologo-economicheskaya systema ekspluatatsii bioresursov c asymmetrichnymi uchastnikami [Environmental-economic system of

bioresource use with assymetric agents]. Trudy KarNTs RAN [Trans. KarRC RAS]. 2016. No. 8. P. 91-97. doi: 10.17076/mat441

4. Sterligova O. P., Pavlov V. N., Il'mast N. V., Pavlovskii S. A., Komulainen S. F., Kuchko Ya. A. Ecosistema Syamozera (biologicheskii rezhim, ispol'zovanie) [Ecosystem of Lake Syamozero: biological regime and use]. Petrozavodsk: KarRC RAS, 2002. 119 p.

5. Pusillo L., Tijs S. E-equilibria for multicriteria games. Annals of ISDG. 2013. Vol. 12. P. 217-228.

СВЕДЕНИЯ ОБ АВТОРЕ:

Реттиева Анна Николаевна

зам. директора по научной работе, д. ф.-м. н. Институт прикладных математических исследований КарНЦ РАН, Федеральный исследовательский центр «Карельский научный центр РАН» ул. Пушкинская, 11, Петрозаводск, Республика Карелия, Россия, 185910 эл. почта: annaret@krc.karelia.ru тел.: (8142) 766312

6. Rettieva A.N. Multicriteria dynamic games. International Game Theory Review. 2017. Vol. 1(19). P. 1750002.

7. Shapley L.S. Equilibrium points in games with vector Payoffs. Naval Research Logistic Quarterly. 1959. Vol. 6. P. 57-61.

8. Voorneveld M., Grahn S., Dufwenberg M. Ideal equilibria in noncooperative multicriteria games. Mathematical Methods of Operations Research. 2000. Vol. 52. P. 65-77.

Received April 23, 2018

CONTRIBUTOR:

Rettieva, Anna

Institute of Applied Mathematical Research, Karelian Research Centre, Russian Academy of Sciences 11 Pushkinskaya St., Petrozavodsk, Karelia, 185910, Russia

e-mail: annaret@krc.karelia.ru tel.: (8142) 766312

i Надоели баннеры? Вы всегда можете отключить рекламу.