Везде
Экономические и эволюционные аспекты оптимального рыбного промысла
Оглавление
Аннотация Оценить Содержание публикации
Библиография Комментарии
Поделиться
QR
Метрика
Экономические и эволюционные аспекты оптимального рыбного промысла
Аннотация
Код статьи
S042473880005778-4-1
DOI
10.31857/S042473880005778-4
Тип публикации
Статья
Статус публикации
Опубликовано
Авторы
Ильичев Виталий Григорьевич 
Должность: главный научный сотрудник
Аффилиация: Южный Научный Центр РАН
Адрес: Ростов-на-Дону, РФ
Выпуск
Том 55 Номер 3
Страницы
86-99
Аннотация

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

Ключевые слова
многолетний промысел, оптимизация, пространственная адаптация, внутренние цены.
Источник финансирования
Статья подготовлена в рамках реализации ГЗ ЮНЦ РАН (проект АААА-А18-118122790121-5), а также гранта Российского фонда фундаментальных исследований (проект 18-01-00453).
Классификатор
Получено
13.08.2019
Дата публикации
22.08.2019
Всего подписок
92
Всего просмотров
1813
Оценка читателей
0.0 (0 голосов)
Цитировать   Скачать pdf
1

ВВЕДЕНИЕ

Проблема разработки стратегии промысла является популярной в теоретических и прикладных исследованиях (Скалецкая и др., 1979; Абакумов, Израильский, 2017; Clark, 2010; и т.д.). В основном здесь обсуждаются вопросы максимизации многолетнего вылова или устойчивости равновесий в рамках жестких математических моделей. В работе (Ильичев и др., 2000) авторы обнаружили эффект стабилизации хаотической динамики (подробнее см. в (Якобсон, 1976)) в модели Риккера (Ricker, 1954) при оптимальном промысле. Иногда при обсуждении проблем промысла используется и теоретико-игровой подход (Mazalov, Rettieva, 2004).
2 В традиционных моделях оптимального вылова, как правило, не учитывается эволюционный ответ рыбной популяции на промысел. Так, пространственно-временной характер вылова может вызвать изменение маршрута миграции рыбы. Требование учета механизма адаптации вносит новизну в привычную задачу оптимизации и вызывает новые содержательные постановки задач. Основная трудность здесь заключается в эффективной реализации эволюционного процесса.
3 Напомним, что в известных работах (например, (Свирежев, Пасеков, 1982)) предлагается прямолинейное классическое копирование биологических факторов. А именно, сначала задается достаточно большой набор близких субпопуляций (мутантов) с несколько различными параметрами, а затем в результате конкуренции выявляется победитель — субпопуляция с наилучшим параметром. Такой пассивный подход требует задания тысячи и более переменных, и а компьютерная реализация — значительных затрат машинного времени. В работе (Ильичев, Ильичева, 2014) был предложен быстрый метод адаптации, когда задается малое число мутантов, а параметры исходной популяции активно изменяются. Вероятно, эти модельные механизмы можно использовать и в моделях эволюционной экономики (Экономическая трансформация…, 2004).
4 Далее будет описана мягкая эколого-эволюционная модель, в которой маршрут миграции популяции рыб может целесообразно изменяться в зависимости от распределения кормовых ресурсов и вылова. Оказывается, во многих вопросах пространственной адаптации ключевую роль играет не финальная матрица миграции М, а ее положительный собственный вектор . Этот (перроновский) вектор удовлетворяет соотношению , и при положительной М оказывается единственным (Робертс, 1986). Удобно считать, что сумма всех компонент равна 1, тогда является относительным временем пребывания популяции в районе i. Установлено (Ильичев, Ильичева, 2018), что всякое адаптивное поведение сводится, по сути, к перестройке времен пребывания в некотором районе. Таким образом, матрица является лишь одной из многих форм конкретной миграции, а истинным содержанием миграции будет перроновский вектор.
5 Оказывается, задачу оптимального промысла можно дополнить экономической начинкой, основанной на идее внутренних цен рыбной популяции для различных районов. Это позволяет задать эффективные экономические механизмы: налогообложение, торговые обмены и др. В данной статье приведена иллюстрация указанных подходов.
6

  • ВОДОЕМ С ЗАПОВЕДНИКОМ. МИГРАЦИЯ И АДАПТАЦИЯ РЫБ

    		•
    7 Пусть водоем разбит на две зоны (района) U и Z. В зоне U допускается вылов, а Z расположен заповедник (неприкосновенный запас). Предполагается, что рыбная популяция мигрирует между этими районами. Формально, пусть xt и yt — текущие численности данной популяции в зонах U и Z соответственно. Тогда ее перемещение за единицу времени задается линейной моделью:
    8 (1)
    9 где . Поскольку общая численность популяции при перемещениях не меняется, сумма элементов по каждому столбцу равна 1. Такие матрицы называют марковскими. Удобно представлять двумерные марковские матрицы в виде:
    10 (2)
    11 где . Внедиагональные элементы и характеризуют интенсивность движения популяции. Так, — это доля популяции, которая может перейти из района Z в U. Аналогично, определяет долю популяции, переходящую из U в Z. Воспользуемся представлением матриц вида (2) как некоторых точек плоскости:
    12 (3)
    13 В совокупности они заполняют весь квадрат с вершинами . Каждая точка является выпуклой комбинацией этих вершин. По сути, они составляют своеобразный базис. Приведем базисные матрицы:
    14 . Считаем, что исходная популяция с матрицей миграции может порождать мутантов с маршрутами , близкими к . А именно:
    15 (4)
    16 где и мало, — одна из базисных матриц.
    17 В процессе конкуренции исходной популяции с семейством таких мутантов может выявиться сильный мутант, численность которого превосходит численность исходной популяции в 10 и более раз. Если сильных мутантов несколько, выбираем самого многочисленного. А для всякого более слабого мутанта порождаем новый маршрут (4) с другими и .
    18 Далее, считаем, что сильный мутант задает направление естественного отбора. Формально пусть — матрица миграции сильного мутанта. Тогда деформируем матрицу исходной популяции следующим образом:
    19 (5)
    20 где и мало, параметр задает скорость адаптации.
    21 Пусть и — две марковские матрицы. Тогда к можно сколь угодно близко подойти с помощью композиции формул вида (4) и (5) при подходящем выборе базисных матриц (Ильичев, 2012). Содержательно, это означает практическую достижимость любой марковской матрицы с помощью композиции формул вида (4) и (5).
    22

  • ВОДОЕМ С ЗАПОВЕДНИКОМ. ОПТИМИЗАЦИЯ И ВНУТРЕННИЕ ЦЕНЫ

    • 23 Введем параметр , который характеризует долю зоны U в водоеме: при район промысла мал, а при велик. Обозначим через K общую величину корма в водоеме. В этом случае и — это величина корма в зоне U и в зоне Z соответственно.
    24 Рассмотрим простейший случай, когда запасы корма фиксированы и не изменяются в каждом районе.
    25 В оптимизационном блоке модели рассматривается промысел в зоне U на многолетний период. Обозначим через — оптимальный доход за лет, x и y — начальные запасы рыбы равны в зоне U и в зоне Z соответственно. Полагаем на всем допустимым
    26 (6)
    27 Здесь — коэффициент дисконтирования (полагаем ); — величина вылова в момент t; ; p — вогнутая и монотонно возрастающая, гладкая функция полезности с , например, . Для определенности пусть
    28 Динамику численности неподвижной популяции формально представим в виде . Используем простейший вариант нелинейного роста: пусть f — гладкая, строго возрастающая и вогнутая функция от x; и . Последнее неравенство гарантирует положительную численность популяции для всех t. Введем требование на ограничение роста численности популяции, например, пусть популяция растет в медленном темпе, т.е. при некотором имеет место для каждого K и для всех достаточно больших z. В частности, ниже будет использована зависимость Контуа (Contois, 1959): где d = 0,25 и r = 3 — скорости смертности и роста.
    29 Обозначим и . При заданной матрице миграции и последовательности ежегодных выловов динамика фазовых переменных в водоеме с заповеднике задается системой:
    30 (7)
    31 Компьютерные расчеты показали, что траектория (7) при оптимальной стратегии вылова стремится к некоторому равновесию. В облегченной версии вылова, когда при фиксированном , это явление можно строго доказать.
    32 Утверждение 1. Пусть существует положительное равновесие в модели (7). Тогда оно оказывается глобально устойчивым.
    33 Обоснование данного утверждения приведено в Приложении, п. 1.
    34 Маршрут исходной популяция является эволюционно устойчивой матрицей (ЭУ-матрица), если данная популяция не вытесняется своими мутантами с близкими к маршрутами.
    35 В природе могут реализоваться только ЭУ-матрицы миграции. В наиболее сильном варианте этого понятия число мутантов может быть любым. В простейшем варианте можно ограничиться и одним мутантом. В приведенных ниже компьютерных расчетах использовались два мутанта, в которых перебор базисных матриц происходит в противоположных направлениях.
    36 Актуальна задача, насколько зависят получаемые результаты от числа генерируемых мутантов? Предварительные расчеты показали, что найденные здесь закономерности слабо обусловлены выбранным семейством мутантов. Строго говоря, это пока открытая проблема.
    37 Для краткости, запишем модель (7) в форме:
    38 (8)
    39 Очевидно, что все функции в (8) являются строго возрастающими и вогнутыми по каждой переменной. Теперь согласно известной схеме динамического программирования при построении целесообразно воспользоваться рекурсией Беллмана (Беллман, 1960):
    40 (9)
    41 где . Каждая является непрерывной функцией. Более того, рекурсия (9) сходится к некоторой предельной непрерывной функции (обоснование см. в Приложении, п. 2).
    42 Отметим, что при бесконечном горизонте планирования промысла оптимальный вылов зависит только от значений фазовых переменных. Напротив, при конечном горизонте планирования оптимальный вылов существенно определяется моментом времени t. Поэтому исследование бесконечномерной версии данной задачи удобнее и проще.
    43 Далее, наследует хорошие свойства функции , поэтому функция Беллмана строго возрастает и вогнута по обеим переменным. Согласно известной теореме Лебега (Колмогоров, Фомин, 1972), всякая монотонная функция почти всюду дифференцируема. Известно, что возможные точки недифференцируемости вогнутой функции устроены просто. Так, всегда существуют и при этом выполняется неравенство Имеются веские основания полагать, что здесь добавка является почти лишней (Рохлин, 2000). Ниже для краткости также будем опускать эту добавку. Будем говорить: из вогнутости по переменной x следует, что убывает по x. Отметим, что из строгой вогнутости и непрерывности функции Беллмана вытекает непрерывность управления (Ланкастер, 1972).
    44 Математическую функцию можно трактовать как большую экономическую функцию полезности, когда начальные запасы рыбной популяции по районам равны x и y. Поэтому частные производные и можно трактовать как (внутренние) цены рыбы в первом и во втором районах, соответственно. Любопытно было бы строго выяснить, какая из этих цен выше? Расчеты показали: промысловая рыба существенно дороже заповедной.
    45 Важное значение имеет следующее утверждение.
    46 Утверждение 2. Пусть в точке оптимальный вылов , положителен, тогда существует и реализуется соотношение
    47 (10)
    48 Обоснование данного утверждения о связи большой и малой функций полезности приведено в Приложении, п. 3.
    49 В частности, при оптимальном вылове имеет место Разумеется, верно и обратное: если в точке выполняется неравенство то в ней оптимальный вылов положителен. Значит, при оптимальный вылов равен нулю.
    50 Величину можно использовать в качестве единицы налога на выловленную рыбу. Так, рассмотрим решение простой задачи оптимизации
    51 по всем (11)
    52 и попутно выясним, насколько сильно отличаются управления в (11) и (6). Для задачи (11) существует два случая.
    53
    1. Для целевая гладкая функция H:
    54 (12)
    55 является строго убывающей, поскольку ее производная отрицательна . Поэтому здесь решение оказывается оптимальным.
    56
    1. Для в силу строгой вогнутости H и искомое управление должно удовлетворять уравнению , т.е. . Очевидно, что данное соотношение совпадает с выражением (10). Таким образом, решение проблемы многошаговой оптимизации (6) свелось к решению задачи одношаговой оптимизации (11).
    57 Приведем ряд численных расчетов оптимального вылова, представляющих определенный теоретический интерес. Так, при одновременно малых вылов будет равен нулю. А если значение переменной y фиксировано, то функция возрастает по x.
    58 Пусть ежегодный вылов выбран оптимальным в соответствие с решением задачи (9). Оказывается, что и в этом случае асимптотические значения фазовых переменных не зависят от задания начального условия и совпадают с равновесием модели (7).
    59 Любопытно было бы ввести в модель некоторые глобальные экономические характеристики водоема с заповедником. Например, вычислить средние значения внутренних цен или по всем . Однако оказалось, что такие значения не очень чувствительны к изменению параметров водоема. Более эффективными показали себя характеристики и
    60

  • КОАДАПТАЦИЯ ПРОМЫСЛА И МИГРАЦИИ

    		•
    61 Поиск стратегии оптимального вылова и построение ЭУ-матриц миграции включает два последовательных этапа (блока).
    62 1. Блок «Оптимизация». При заданной матрице миграции M на основе метода динамического программирования производится поиск оптимального вылова u. Пусть x и y — соответствующие численности популяции в двух районах. Тогда можно построить отображение
    63 (13)
    64 2. Блок «Эволюция». При фиксированной функции вылова реализуется конкурентное взаимодействие исходной популяции (с матрицей миграции M) с несколькими мутантами (с матрицами, близкими к M). Когда в результате отбора выявляется сильный мутант, тогда согласно (5) производим некоторую деформацию M до некоторой матрицы . Поэтому реализуется отображение
    65 (14)
    66 Если в результате многократного повторения этапов 1 и 2 этот процесс сходится, то определяется эволюционно-устойчивая матрица M* и вылов u* при заданном значении корма в районах U и Z.
    67 Методология расчетов. Рассматривалось несколько типов заповедника Z: от самого малого (10% водоема) до большого (90%). В каждом случае произведены соответствующие расчеты при вариации начальных значений на некоторой сетке. В качестве итоговых значений определялись средние значения и .
    68 Вычисления показали, что финальные значения матрицы миграции могут зависеть от выбора начального маршрута рыбной популяции. Это явление множественности ЭУ-состояний характерно для многих проблем эволюционной экологии (Ильичев, Ильичева, 2018).
    69 Установление эволюционной устойчивости иногда противоречит наивным представлениям об адаптации маршрутов популяции. Так, например, рассмотрим ситуацию, когда и , а в качестве промысла используется хищная функция . Казалось бы, полный вылов в зоне U должен заставить популяцию полностью покинуть этот район. Однако расчеты показали, что — финальная матрица, и .
    70 Напомним, что — это вероятность перехода рыбных особей из района j в i. В данной матрице имеет место неравенство 0,26 < 0,47. Значит, район U менее предпочтителен, чем Z, но все же не настолько, чтобы его окончательно оставить. В целом, ЭУ-матрицы имеют скорее оборонное значение и обеспечивают исходной популяции устойчивое сосуществование со своими мутантами. Здесь во многих случаях достигается и максимум численности популяции.
    71

  • ВОДОЕМ С ЗАПОВЕДНИКОМ. ЭВОЛЮЦИЯ И ЭКОНОМИКА ПРОМЫСЛА

    		•
    72 Пусть кормовой ресурс промысловой зоны характеризуется параметром , который принимает дискретные значения от 0,1 до 0,9. Этот коэффициент является долей промысловой зоны в общей акватории водоема. Представляют интерес анализ поведения экономических и экологических элементов при вариации . Основные результаты компьютерных расчетов представлены в табл. 1.
    73 Таблица 1. Динамика эколого-экономической системы с заповедником
    74
    Пара-метры 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
    0,6 0,6 0,6 0,6 0,6 1,9 3,0 4,8 5,2
    10,1 8,9 7,8 6,8 5,4 4,6 3,0 2,1 0,9
    0,5 0,5 0,5 0,5 0,5 1,2 1,9 2,5 2,7
    5 5 5 5 5 2,1 1,3 0,9 0,7
    0,7 0,7 0,7 0 0,2 0,9 0,4 0,4 0,3
    75 Примечание. В таблице выделено критическое значение , для которого слева стоят низкопродуктивные зоны вылова, а справа — высокопродуктивные.
    76 Обсудим поведение данных в табл. 1 в свете роста параметра [[[image104]]][[[image104]]]. Для него имеют место следующие ясные закономерности.
    77
    1. Численность рыбной популяции x в промысловой зоне мала и слабо изменяется при . Поэтому величина вылова почти совпадает с численностью рыбы в данном районе. Здесь промысловая зона существует только за счет дотаций заповедника. Однако при численность рыбы и ее вылов в промысловой зоне значительно возрастает, и этот район становится самодостаточным.
    2. В заповеднике численность рыбной популяции x равномерно уменьшается без особых скачков.
    3. Напомним, π-перроновский вектор финальной матрицы миграции. С ростом увеличивается кормовой ресурс промысловой зоны, и, несмотря на вылов, растет и время пребывания рыбной популяции в ней.
    4. При запасы рыбы в промысловой зоне незначительны, поэтому цена на дефицитную рыбу будет довольно высока. Но после критического значения λ численность рыбы здесь возрастает, и происходит снижение цены .
    78 Для цены рыбы в заповеднике не обнаружено явных закономерностей. Вероятно, здесь требуются расчеты на более мелкой сетке. Однако можно утверждать: в рамках данной модели рыбной популяции цена существенно меньше . Представляет интерес построение аналога формулы (10), связывающего оптимальный вылов с внутренней ценой заповедной рыбы. Интуитивно ясно, что заповедник — это резерв будущей промысловой рыбы.
    79 Приведем один из таких вариантов. При условии, что и цена является линейной комбинацией величин и с коэффициентами, зависящими от фазовых переменных. Это выражение весьма сложно, но при один из коэффициентов обращается в нуль и данная формула приобретает простой вид:
    80

    81 Обоснование этой формулы громоздко и опускается.
    82

  • ВОДОЕМ БЕЗ ЗАПОВЕДНИКА.

    • ЭКОНОМИКА ПРОМЫСЛОВОЙ КОНКУРЕНЦИИ
    83 Рассмотрим задачу, когда акватория поделена на две зоны и в каждой из них ведет промысел свой собственник. Такими участниками конкурентного промыслового процесса могут являться, в частности, два соседних государства.
    84 Перейдем к формализации данной задачи. Пусть — численности рыбной популяции в первой и второй зонах водоемах в момент времени t; u и v — соответствующие этому состоянию выловы:
    85 (15)
    86 p и q — функции полезности для первой и второй зон соответственно. Будем считать, что конкуренция не злонамерена, поэтому возможна общая целевая функция, например, максимизация суммарной полезности. Тогда рекурсия Беллмана приобретает вид
    87 (16)
    88 с условием (15). Разумеется, (обычно полагают ); с учетом вылова реализуются соотношения
    89 (17)
    90 Итерационный процесс (16) сходится к некоторой непрерывной функции B(x, y), которая, по сути, является большой функцией полезности для всей популяции. Она наследует основные свойства малых функций полезности p и q: монотонность, вогнутость и др.
    91 Пусть, например, значение переменной фиксировано, тогда функция является монотонно возрастающей, и следовательно, по теореме Лебега оказывается дифференцируемой почти для всех x. Здесь имеются веские основания полагать, что B(x, y) является дифференцируемой всюду (Ильичев и др., 2000). Поэтому, как и ранее, в дальнейшем слово почти будем опускать. Далее, в силу вогнутости, производная одноместной функции по x убывает. Разумеется, аналогичные результаты справедливы и для относительно переменной y.
    92 Зададим внутренние цены единицы рыбной популяции в районах:
    93 (18)
    94 Такие цены зависят от значений фазовых переменных. Здесь имеют место результаты исследования связи оптимальных выловов с функцией Беллмана.
    95 Утверждение 3. Пусть тогда выполняется равенство
    96 (19)
    97 поэтому имеет место .
    98 И обратно, пусть выполняется неравенство , тогда обязательно должно быть . Иначе, согласно (19), получаем .
    99 Утверждение 4. Пусть тогда выполняется равенство
    100 . (20)
    101 И неравенство является критерием условия .
    102 На данных соотношениях основаны простые экономические механизмы построения оптимальных выловов. Так, пусть в состоянии задана цена рыбы , которую будем использовать в качестве налога на единицу выловленной рыбы в первом районе. В этом случае максимизация ежегодной прибыли для собственника данного района сводится к решению задачи
    103 по всем (21)
    104 Обсудим решение (21). Для этого рассмотрим вспомогательную (целевую) функцию . Она является строго вогнутой и .
    105 При имеет место для всех . Значит, аргумент доставляет максимум убывающей функции H. При имеем . Несколько сложнее оценить знак . Здесь согласно формуле (19) существует из промежутка , при котором , поэтому, в силу убывания получаем . Значит, в точке максимума функции H на J выполняется условие , т.е. Следовательно, решение одношаговой задачи (21) совпадает с решением многошаговой задачи оптимизации (16) при .
    106 Аналогичная конструкция налогообложения применима и к поиску оптимального вылова во втором районе. Здесь при следует рассмотреть одношаговую задачу по всем .
    107 Разумеется, и — локальные цены, зависящие от конкретного состояния . Представляет интерес построение глобальных экономических показателей. Так, обозначим через T — достаточно большой период оптимального промысла (16). Зададим некоторое множество возможных начальных значений . Например, такой естественной областью является дискретная решетка в достаточно большом прямоугольнике вида . Под действием динамической системы (17) образуется соответствующее множество конечных значений . Пусть — осреднение семейства финальных векторов. Отметим, что переменные в системе (17) при оптимальных выловах стремятся к некоторому равновесию, которое не зависит от выбора начальных условий. Поэтому точка является таким равновесным состоянием.
    108 Теперь положим и . Это глобальные ценовые характеристики отдельных районов водоема в зависимости, например, от величины кормовых ресурсов (рис. 1). Графики локальных цен представляются вполне естественными. В частности, при первый район располагает малым кормом, поэтому здесь численность рыбы мала. И значит, ее цена будет велика. Отметим, что, например, при возможна выгодная продажа дешевой рыбы (собственника второго района) собственнику первого района.
    109

    Рис. 1. Динамика средних внутренних цен (сплошная линяя) и зависимости их от параметра λ (пунктирная)
    110

  • ОПТИМАЛЬНЫЙ ПРОМЫСЕЛ КАК ИГРА В ПОДДАВКИ

    		•
    111 В рамках предыдущей задачи о конкурентном промысле двух государств обсудим возможность увеличения дохода, например, страны 1 путем временного добровольного уменьшения ее величины вылова. Будем действовать в духе иерархического управления с уточненными прежними соглашениями о стратегиях вылова. Считаем, что центр на основе решения задачи (16) выдает рекомендации этим государствам о квоте вылова, которая не должна превышать величин и , зависящих от состояния рыбных запасов . Разумеется, можно в какие-то годы вылавливать и меньше (т.е. держать свою акваторию “под паром”).
    112 Рассмотрим следующую симметричную ситуацию. Пусть в начальный момент времени матрица миграции имеет вид , и кормовые ресурсы акваторий данных стран одинаковы, т.е. .
    113 Если страны действуют согласно рекомендациям центра (с учетом вида и решения задачи (16)), матрица миграции сохраняется, и устанавливаются значения:
    114 (22)
    115 Пусть теперь страна 1 не вылавливает рыбу (), а страна 2 ведет промысел в прежнем объеме . Таким образом, происходит своеобразное заманивание рыбной популяции в акваторию 1. В итоге матрица миграции преобразуется к виду Иными словами, акватория 1 стала более предпочтительной, чем акватория 2 (в соотношении 3 : 2). С учетом вида матрицы центр согласно решению задачи (16) выдает новые рекомендации оптимального вылова:
    116 (23)
    117 Теперь ситуация в акватории 1 улучшилась. Более того, если страны станут придерживаются вылова (23), то матрица не портится. Следовательно, является эволюционно-устойчивой при управлениях (23).
    118 После процедуры заманивания вылов в стране 1 всегда будет больше, чем в стране 2. Эффективность этой процедуры обнаружена для других начальных матриц миграции:
    119

    120 А поскольку , у коммерческих структур страны 1 возникает возможность дополнительной торговой выгоды: продать стране 2 свою дешевую рыбу (в смысле ее цена стала ниже, чем в другой стране). Вероятно, чтобы восстановить прежний статус-кво страна 2 должна, в свою очередь, применить процедуру заманивания.
    121 Отметим, что процедура заманивания является эффективной и при незначительном уменьшении вылова (75% от рекомендуемого).
    122

    ЗАКЛЮЧЕНИЕ

    123 Перечислим основные результаты, полученные в данной статье.
    124
    1. Сравнительные расчеты показали, что с учетом поведенческой адаптации рыб их вылов снижается. По сути, данный процесс адаптации — это поиск компромисса между наличием корма и интенсивности вылова.
    2. В случае когда водоем разделен на две зоны (промысловую и заповедник) и промысловая зона невелика, она становится дотационной от рыбных запасов заповедника, и тогда оптимальный вылов заключается в практически полном изъятии рыбы. А если данная зона велика, она становится самодостаточной, и здесь вылов не превышает половины рыбных запасов.
    3. В случае двух промысловых зон, каждая из которых облавливается своим собственником, возможна парадоксальная стратегия долгосрочного вылова. А именно, первый собственник может на некоторое время сильно сократить свой вылов, что приведет к эволюционной деформации маршрута миграции рыб с предпочтением первой зоны. После этого он может перейти к стратегии оптимального вылова, который будет превышать оптимальный вылов второго собственника.
    4. В статье введено естественное понятие внутренних цен для рыбных запасов в том или ином районе. Внутренние цены могут быть использованы в качестве налога на единицу выловленной рыбы. В этом случае проблема многолетней оптимизации сводится к решению задачи максимизации одногодичного вылова. Пространственная неоднородность внутренних цен позволяет конструировать разнообразные торговые механизмы обмена выловленной рыбой.
    125 Приложение
    126 1. Доказательство утверждения 1.
    127 На точках фазового пространства (7) определим отношение частичного порядка. А именно, для точек и положим , если одновременно выполняются неравенства и .
    128 В функциональном анализе принято считать (Опойцев, 1986), что данное отношение порождено конусом (положительный квадрант).
    129 Обозначим через — равновесие (7). Пусть выбран некоторый параметр . Зададим точки (слабая точка) и (сильная точка). Очевидно, что . Теперь построим λ-параметрическое семейство вложенных прямоугольников окружающих равновесие : Внутренность можно трактовать, как промежуточное множество (конусный отрезок) между слабой точкой и сильной точкой .
    130 Разумеется, совпадает с равновесием , а при допустимой вариации эти прямоугольники полностью заполняют положительный ортант.
    131 Обозначим через — отображение (за единицу времени), задаваемое моделью (7). В силу строгого роста правых частей (7), выполняется свойство . И более того, из вытекает .
    132

    Рис. 2. Образ (выделен фоном) куба [[[image209]]][[[image209]]] под действием отображения [[[image1]]][[[image1]]]
    133 Покажем, что образ строго вложен в (рис. 2). Здесь достаточно установить два строгих предпочтения:
    134 — слабая точка идет вперед;
    135 — сильная точка идет назад.
    136 Ниже при всех i решающее значение имеет следующая лемма.
    137 Лемма 1. Для всех и справедливы неравенства:
    138 (П1)
    139 (П2)
    140 Доказательство. Исходя из строгой вогнутости при получаем . Поскольку , то последнее соотношение равносильно (П1). Положим , тогда из неравенства (П1) вытекает (П2). ▄
    141 Теперь подставим слабую точку в уравнения расширенной модели. С учетом леммы 1 получаем оценки для новой точки :
    142

    143 Аналогично устанавливается соотношение Следовательно, слабая точка идет вперед. Также показываем, что сильная точка идет назад. Отсюда сразу вытекает глобальная в устойчивость и единственность равновесия .
    144 2. Доказательство сходимости рекурсии (9) к непрерывной функции
    145 Лемма 2. При достаточно большом S на границе прямоугольника векторное поле (7) не направлено наружу.
    146 Доказательство. Обозначим через и две точки из П. Разумеется, П можно трактовать как промежуточное семейство точек для этих вершин. Поэтому достаточно установить, что точки и не покидают П. Действительно, сначала имеем .
    147 Далее, для -компоненты точки с учетом допущения о медленном росте получаем оценку А далее, для y-компоненты точки аналогично получаем ▄
    148 Обозначим через семейство всех непрерывных функций от двух переменных с чебышовой метрикой по всем из П. Известно, что это пространство является полным (Колмогоров, Фомин, 1972).
    149 C учетом обозначений (8) на определим оператор : по всем . Согласно (Ланкастер, 1972) и лемме 2 функция непрерывна, и поэтому отображает в себя. Более того, покажем, что — сжатие.
    150 Так, пусть и — две функции из данного пространства. Формально имеем при своем оптимальном .
    151 В рамках предыдущих обозначений справедлива оценка . После вычитания последних соотношений находим . Подобным образом, используя управление , устанавливаем . Окончательно, находим .
    152 Поскольку , оператор является сжатием. И следовательно, для него существует единственная неподвижная точка (непрерывная функция двух переменных). Она является решением уравнения Беллмана с бесконечным горизонтом управления.
    153
    1. Доказательство утверждения 2
    154 В данной точке согласно обозначениям (8) и (9) имеет место соотношение
    155 . (П3)
    156 Пусть достаточно малое число, тогда в точке вылов является допустимым независимо от знака . Поскольку он не может быть лучше оптимального, то справедливо неравенство
    157 (П4)
    158 Вычитая из (П4) равенство (П3), находим . Отсюда при и имеет место
    159 (П5)
    160 а при и —
    161 .(П6)
    162 Поскольку для вогнутой по переменной x функции всегда имеет место неравенство , существует , а c учетом (П5) и (П6) выполняется требуемое равенство .

    Библиография

    1. Абакумов А.И., Израильский Ю.Г. (2017). Стабилизирующая роль структуры рыбной популяции в условиях промысла при случайных воздействиях среды обитания // Компьютерные исследования и моделирование. Т. 9. № 4. C. 609–620.

    2. Беллман Р. (1960). Динамическое программирование. М.: Издательство иностранной литературы [Пер. Bellman R. (1957). Dynamic Programming. New Jersey: Princeton University Press.]

    3. Ильичев В.Г., Рохлин Д.Б., Угольницкий Г.А. (2000). Об экономических механизмах управления биоресурсами // Известия РАН. Теория и системы управления. № 4. С. 104—110.

    4. Ильичев В.Г., Ильичева В.В. (2014). Пространственная адаптация и оптимальный промысел рыбных популяций // Экономика и математические методы. Т. 50. № 3. С. 119–129.

    5. Ильичев В.Г., Ильичева О.А. (2018). Пространственная адаптация популяций в моделях экологии // Биофизика. Т. 63. Вып. 2. С. 373–381.

    6. Колмогоров А.Н., Фомин С.В. (1972). Элементы теории функций и функционального анализа. М.: Наука.

    7. Ланкастер К. (1972). Математическая экономика. М.: Советское радио.

    8. Опойцев В.И. (1986). Нелинейная системостатика. М.: Наука.

    9. Робертс Ф.С. (1986). Дискретные математические модели с приложением к социальным, биологическим и экологическим задачам. М.: Наука.

    10. Рохлин Д.Б. (2000). Производная решения функционального уравнения Беллмана и цена биоресурсов // Сибирский журнал индустриальной математики. Т. 3. № 1. С. 169–181.

    11. Свирежев Ю.М., Пасеков В.П. (1982). Основы математической генетики. М.: Наука.

    12. Скалецкая Е.И., Фрисман Е.Я., Шапиро А.П. (1979). Дискретные модели динамики численности популяций и оптимизация промысла. М.: Наука.

    13. Экономическая трансформация и эволюционная теория Й. Шумперта (2004). В сб. статей V международного симпозиума по эволюционной экономике. Россия, г. Пущино, 25–27 сентября 2003 г. М.: Институт экономики РАН.

    14. Якобсон М.В. (1976). О свойствах однопараметрического семейства динамических систем // Успехи математических наук. T. 31. № 2 (188). С. 239–240.

    15. Clark C.W. (2010). Mathematical Bioeconomics. The Mathematics of Conservation. New Jersey: J. Wiley and Sons Publ.

    16. Contois D.E. (1959). Kinetics of Bacterial Growth Relationship between Population Density and Specific Growth Rate of Continuous Culture // Journal of General Microbiology. No. 1–2. P. 40–50.

    17. Mazalov V.V., Rettieva A.N. (2004). A Fishery Game Model with Migration Reserved Territory Approach //Game Theory and Applications. Vol. 10. P. 97–108.

    18. Ricker W.E. (1954). Stock and Recruitment // Journal Fisheries Research Board of Canada. Vol. 11. No. 5. P. 569–623.

    Комментарии

    Сообщения не найдены

    Написать отзыв
    Перевести