Добавил:
Студент Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ААВасин Математические_модели_рынков_и_аукционов

.pdf
Скачиваний:
66
Добавлен:
29.05.2024
Размер:
3.15 Mб
Скачать
Q Qij , i, j L

130

Рынки с сетевой структурой

рынка. Таким образом, для решения задачи (10.2) могут быть использованы известные численные методы (см. Гасников А.В., 2015).

Далее рассмотрим задачу поиска оптимального множест-

ва L*. В ряде случаев для ее решения удается использовать свойства субмодулярности или супермодулярности функции общественного благосостояния (см. Хачатуров В.Р., 1989). Функция

W R , определенная для каждого подмножества R конечного

множества L, является субмодулярной (соответственно супермодулярной) на L , если

L1 , L2 L W L1 W L2 W L1 L2 W L1 L2 .

Для любого вектора пропускных способностей рассмотрим исходный рынок без затрат на строительство

и при ограничениях

 

qij

 

Qij ,

i, j L . Обозначим

через

 

 

pi Q , i N , равновесные цены, соответствующие

такому

рынку.

В работе Vasin A., Dolmatova M. (2016) получен следующий результат относительно супермодулярности функции благосос-

тояния для рынка типа цепи, где L i,i 1 ,i 1, ,n 1 .

Теорема 11.2. Пусть начальные цены pi Q0 ,i 1,.., n , мо-

нотонно убывают по i . Тогда функция W R является супер-

модулярной. Сложность поиска оптимального множества L* при Q0 0 не превышает n 1 n / 2 .

Данная оценка сложности связана с тем, что при Q0 0 не нужно перебирать при добавлении несвязанные подмножества L. В общем случае (при Q0 0 ) задача остается NP-трудной.

Часть III

131

Рассмотрим алгоритм, который позволяет в среднем эффективно находить оптимальный набор L* для случайно выбранной супермодулярной функции W R . Алгоритм работает с нижней

оценкой Lmin и верхней оценкой Lmax множества L* . На каждом шаге производится попытка либо расширить текущее множество Lmin , добавляя очередное подмножество S Lmax \ Lmin , либо сузить текущее множество Lmax , исключая очередное подмножество S Lmax \ Lmin . Если на текущем шаге количество исключае-

мых элементов k2 больше, чем количество k1 добавляемых, то выполняется операция исключения. В противном случае выпол-

няется

добавление.

Если

k1 k2

 

L

\ L

 

,

то

 

 

 

 

 

 

 

max

min

 

 

 

L* Arg max W Lmin ,W Lmax , решение найдено. Алгоритм основан на следующем утверждении.

Теорема 11.3. Пусть для S Lmax

\ Lmin выполняются ус-

ловия: 1)W Lmin S W Lmin

и 2)

W Lmin R W Lmin

R S. Тогда

Lmin S – уточненная нижняя оценка. Если же

выполняются

условия:

3)W Lmax \ S W Lmax и

4)W Lmax \ R W Lmax R S , то Lmax \ S – уточненная верхняя оценка.

Наихудший случай – когда L* L / 2 и R : R L / 2, R L*,

выполнено W R W и W L \ R W L . В этих условиях

алгоритм может осуществить полный перебор всех подмножеств R L. В то же время для широкого круга задач со случайными

параметрами средняя эффективность данного алгоритма оказывается достаточно высокой. Приведем результаты вычислительного эксперимента для рынка типа цепи с m ребрами, функциями

спроса вида Di max 0,dif pci / 2 и предложения

132

Рынки с сетевой структурой

Si p pci / 2, если p 2dif / ci , Si p dif pci , .

если p 2dif / ci ,i 1,.., m 1

Их разность – линейная функция Si p dif pci . Коэффици-

енты

c

,

d f

выбираются случайно

при условиях

 

i

 

 

i

 

 

pi 1 Q

0 pi Q

0 ,i 1,.., m . Для каждого ребра

k 1, 2,.., m пере-

менные затраты на увеличение пропускной способности составляют ek Q2 . Коэффициенты ek , а также затраты на передачу ekt ,

фиксированные издержки ekf и исходные пропускные способности Qk0 , k 1,.., m , являются случайными.

Рисунок 11.1

На Рис. 11.1 указаны среднее и максимальное число подмножеств R , для которых решалась вспомогательная задача (11.2)

вычисления V R для определения L*. Средняя сложность решения в диапазоне 10 25 линий не превышает 8m . В более широ-

Часть III

133

ком диапазоне 10 50 линий средняя сложность хорошо приближается полиномом f3 x 0.014x3 0.16x2 2.75x.

В общем случае, однако, даже рынки типа цепи не обладают ни одним из указанных свойств (см. Vasin A., Dolmatova M., 2016). Тем не менее, в довольно общих условиях существует возможность разработки эффективных алгоритмов, подобных описанному. Рассмотрим соответствующие понятия для рынка с транспортной структурой общего вида.

Определение 11.1. Ребро q называется дополнительным (соответственно конкурентным) к ребру l , если M L \ l, q

выполнено W M q,l W M q W M l W M .

Теорема 11.4. Пусть M1 l и M2 l - множества дополнительных и конкурентных ребер к ребру l. Тогда

W L1 L2 l W L1 L2 L1 M1 l и L2 M2 l .

Алгоритмы, подобные предложенным для суб- и супермодулярных функций, применимы, если для любого l L можно оп-

ределить M1 l , M2 l

и M1 l M 2 l L \ l . Возможность ус-

тановить полные отношения дополнительности и конкурентности тесно связана со следующим условием.

Определение 11.2. Рассматриваемая модель удовлетворяет

условию инвариантности структуры потока (УИСП), если

Q Q0 , i, j L sign pi Q pj Q sign pi Q0 pj Q0 .

Рассмотрим наиболее общий случай, когда удается определить указанные отношения. Пусть паре N, L соответствует граф

типа дерева. Тогда для любых ребер l1,l2 существует единственный путь L l1,l2 без самопересечений, начинающийся c ребра l1 и заканчивающийся ребром l2. Назовем l2 исходно дополнительным (соответственно конкурентным) к l1, если при начальных пропускных способностях Q0 потоки на этих ребрах имеют

134

Рынки с сетевой структурой

одинаковые (соответственно противоположные) направления относительно пути L l1,l2 . Обозначим L10 l , L02 l множества ис-

ходно дополнительных и конкурентных ребер к l .

Теорема 11.5. Рынок типа дерева удовлетворяет УИСП то-

гда и только тогда, когда

 

 

 

 

 

 

 

 

l i, j

 

sign pi Q0 pj Q0

 

 

sign p

 

Q0

 

,Q

p

j

 

Q0

,Q

 

.

i

0

l

0

 

0

0

 

 

 

L1

L2 l

 

 

L1 l

L2

l

При этом условии l

M1 l

L10 l , M2 l

 

L02 l при лю-

бом увеличении пропускных способностей.

Для графов, включающих цикл, ввести полные отношения конкурентности и дополнительности в общем случае не удается. Так, для графа на Рис. 11.2 в зависимости от конкретных значений параметров модели ребро (3,4) может быть как конкурентным, так и дополнительным для ребра (1,2).

Рисунок 11.2. Граф без полных отношений дополнительности

иконкурентности

12.Двухузловой рынок

вмодели олигополии Курно

Рассмотрим два локальных рынка, соединенных линией передачи. Каждый локальный рынок l 1, 2 характеризуется конеч-

ным множеством Al производителей, | Al | nl , функциями затрат

Часть III

135

Ca (v), a Al , и функцией спроса

Dl ( p) . Обозначим через

k (0, L) коэффициент потерь, который показывает долю утра-

ченного товара (в частности, электроэнергии) при передаче с одного рынка на другой, Q – пропускную способность линии. Рас-

смотрим конкуренцию по Курно в этой модели. Тогда каждый производитель выбирает объем va [0 V a ]

Условия первого порядка равновесия Курно для исхода пер-

вого типа с ценами p*,

p* , такими что

 

1

2

 

 

 

1 p*

p* ,

(12.1)

 

2

1

 

в целом аналогичны условиям (3.1), (3.2) для локального рынка: va * ( pi* Ca ' (va *)) | Di '( pi* ) |

для любого a Ai

такого, что Ca '(0) p* ,

(12.2)

 

i

 

va * 0

при Ca '(0) p* ,

(12.3)

 

i

 

где Ca '(v) [Ca '(v),Ca '(v)] в точках скачка функции предельных

затрат.

Кроме того,

 

va * D( pi* ),

i 1, 2 .

(12.4)

 

Ai

 

 

 

 

 

Получим условия первого порядка для равновесия Курно вто-

рого типа с

 

 

 

 

 

 

 

q (0,Q), p* p* .

 

(12.5)

 

 

 

1

2

 

 

Отметим, что при любом малом изменении стратегии va *

производитель

a A1 остается

на

рынке

с функцией спроса

D1( p1(v )) (D2 ( p1(v )) va ) ,

где цена

 

p1(v )

удовлетворяет

 

A2

 

 

 

 

 

условию va

D1( p1) (D2 ( p1) va )) .

Таким образом, для

A1

 

 

A2

 

 

 

a A1

136Рынки с сетевой структурой

va * ( p1* Ca '(va *)) | D1 '( p1*) 2 D2 '( p1*) | ,

если Ca '(0) p*

,

(12.6)

 

1

 

 

va * 0

при Ca '(0) p* .

(12.7)

 

 

1

 

Аналогично, для производителей на рынке 2 спрос равен

D2 ( p1) 1 (D1( p1) va ) и

 

 

A1

 

 

 

va * ( p* Ca '(va *)) | D2 '( p* ) D1 '( p* ) / 2

|

 

1

1

1

 

 

при Ca '(0) p2* ,

 

(12.8)

 

va * 0 при

Ca '(0) p2* .

(12.9)

Наконец, если ограничение на пропускную способность ак-

тивно, то

 

 

 

 

q Q,

p1* p2* , va* D( p1* ) Q , va* D( p2* ) Q , (12.10)

 

A1

A2

 

 

условия первого порядка –

 

 

 

 

va* ( p* Ca ' (va* )) | Di '( p* ) | при Ca ' (0) p* ,

 

i

i

i

 

 

va* 0 при Ca' (0) p* , a Ai ,

i 1, 2.

(12.11)

 

 

i

 

 

Утверждение 12.1. Каждое равновесие по Нэшу в модели конкуренции Курно для двухузлового рынка принадлежит к одному из описанных трех типов, то есть, либо удовлетворяет условиям (12.1)–(12.4), либо условиям (12.5)–(12.9), либо условиям

(12.10)–(12.11).

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

Часть III

137

Допустим, что ограничение пропускной способности не является существенным (формально Q ), и рассмотрим локальное

равновесие типа 1 (с нулевым перетоком между рынками). При достаточно большом росте объема производства игрока a цена

на рынке 1 уменьшается до уровня p1 p2* . Дальнейшее увели-

чение объема позволяет игроку продавать продукцию также на рынке 2. Рисунок 12.1 показывает функцию спроса для этого производителя.

D

D1( p1) vb* (D2 ( p1) vb* )

A1\a

A2

D1( p1) vb*

A1\a

p*

 

p*

p

2

 

1

1

 

 

Рисунок 12.1

 

Отметим, что

va*

оптимальная

стратегия агента a для

функции спроса D1( p1) vb* . Его оптимальная стратегия при

 

A1\a

 

 

 

спросе D1( p1) vb* (D2

( p1) vb* )

соответствует цене p1

A1\a

 

A2

 

 

такой, что

 

 

 

 

D1( p1) vb* (D2 ( p1) vb* )

(12.12)

 

A1\a

 

A2

v a ( p

Ca (v a )) | D1 ( p ) 2 D2 ( p ) |

 

1

 

1

1

 

138

Рынки с сетевой структурой

Так как левая часть удовлетворяет условиям Теоремы 3.1, оптимальные значения v a , p1 определяются единственным образом из (12.12).

Ниже мы предполагаем, что функции спроса D1( p) , D2 ( p) удовлетворяют условиям Теоремы 3.1:

либо a) D( p) 0 и e( p) p при p ( p, M ), D( p) 0 при p M ,

либо b) D( p) 0

и e( p) p при

p p, lim e( p) L 1/ m ,

 

 

p

где m

общее число производителей на рынке. Эластичность

функции

D1( p) D2 ( p)

возрастает

при

D1( p) 0

и D2 ( p) 0.

Утверждение 12.2. Точка локального равновесия типа a) (которая удовлетворяет условиям (12.1)–(12.4)) не является равновесием Нэша тогда и только тогда, когда для некоторого производителя a A1 выигрыш при спросе D1( p1) vb* (D2 ( p1) vb* ) по цене p1 , заданной согласно

 

A1\a

A2

 

 

 

(12.12),

превышает f a (v * ) ,

или когда для некоторого произво-

дителя

a A2

выигрыш

при

спросе

D2 ( p2 ) vb* (D1( p2 ) vb*)

по оптимальной цене превы-

 

A2 \a

A1

 

 

 

шает f a (v * ) .

 

 

 

 

Теперь рассмотрим локальное равновесие типа b). Чтобы из-

бежать путаницы, обозначим через

p1 , v

цену локального равно-

весия и

производственный

объем

для

этого случая,

а через

pi* ,i 1, 2 – равновесные по Нэшу цены изолированных рынков.

В этом случае функция спроса для производителя a A1 показана на Рисунке 12.2 и аналогична показанной на Рисунке 12.1.

Отличие заключается в том, что цена p2 определяется урав-

нением D2 ( p2 ) vb и p1 p2 / . Уменьшая свое предложе-

A2

Часть III

139

ние, производитель a A1 может увеличить рыночную цену до p2 / . При дальнейшем уменьшении рынки распадаются. Опти-

мальная цена p1

для функции спроса D1( p1) vb

удовлетво-

 

 

 

 

 

 

 

 

 

 

A1\a

 

 

 

 

 

 

 

 

ряет отношению

D1( p

)

 

vb

 

a

 

Ca

 

a

 

D1

 

 

 

 

.

 

 

 

 

 

 

 

v

p

v

p

 

1

 

1

 

 

 

 

1

 

 

 

D

 

 

A1\a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D1( p1) vb (D2 ( p1) vb )

 

 

 

 

 

 

 

 

 

 

A1\a

 

 

 

 

A2

 

 

 

 

 

 

 

 

D1( p1) vb

A1\a

p1

 

p2

p1

 

Утверждение 12.3. Точка локального равновесия типа b) (которая удовлетворяет условиям (12.5)–(12.9)) не является равновесием Нэша тогда и только тогда, когда для некоторого произ-

водителя a A1 выигрыш по оптимальной цене p1 при спросе

D1( p1) vb превышает f a (v) .

A1\a

Отметим, что агенту a A2 невыгодно отклоняться от локального равновесия типа b) с положительным перетоком из рынка 1 на рынок 2. При фиксированных стратегиях других игроков функция спроса выглядит как на Рисунке 12.3.

Путем достаточно большого сокращения объема производства агент a может уменьшить цену и разделить рынки в некоторых случаях. Однако его выигрыш в этом случае понизится, поскольку функция спроса для p p лежит ниже кривой спроса

объединенного рынка.