
Математические модели в экономике.-5
.pdf
C K a 2bd3 .
Рассмотрим случай когда одна из фирм зная, что конкурент выбрал стратегию Курно, действует не по стратегии Курно, а используюя эту ин-
формацию выбирает себе другую стратегию. Такая стратегия называется стратегией Стакельберга. Для простоты будем считать, что d1 d2 d .
Пусть первая фирма дает возможность узнать свой ход x1, тогда вторая фирма выбирает стратегию
x |
(d x1 ) |
. |
(5.9) |
|
|||
2 |
2 |
|
|
|
|
|
Но первая фирма будет действовать тайно по другому с учетом знания стратегии второй фирмы. Подставим в формулу для прибыли первой фир-
мы значение x2 , определяемое формулой (5.9), в результате получим:
W1 bx1(d (x1 x2 )) bx1(d (x1 d x1 )) bx1(d x1 ) . 2
Тогда решив уравнение:
dW1 0 , dx1
Найдем оптимальное значение x1 , соответствующее стратегии Стакельбер-
га, в результате получим:
x1S d2 .
Подставив это значение формулу (8) получим:
x2S d4 .
Общий выпуск товаров по стратегии Стакельберга равен:
xS x1S x2S 34d .
Вычислим прибыли фирм, в результате получим:
W S bxS (d (xS xS )) |
bd 2 |
W K , |
|||
|
|||||
1 |
1 |
1 |
2 |
8 |
1 |
|
|
|
|
|
41

W S bxS (d (xS xS )) |
bd 2 |
W K . |
|||
|
|||||
2 |
2 |
1 |
2 |
16 |
2 |
|
|
|
|
|
Суммарная прибыль по Стакельбергу равна:
W S 316bd .
Установившаяся цена на товар равна
CS a b(x1S x2S ) a 3bd4 .
Как видно, цена на товар в случае Стратегии Стакельберга меньше,
чем в стратегии Курно, т.е. стратегия Стакельберга для потребителя вы-
годней, чем стратегия Курно. Итак, мы видим, что прибыль первой фирмы возрастает. Но здесь обе фирмы могут попасть в ловушку, если обе начнут действовать по стратегии Стакельберга. Тогда их прибыль будет умень-
шаться. Эта стратегия опасна для фирм.
Фирмы могут объединиться в монополию или образовать картель (это тайный сговор двух фирм с целью поддержания заданной цены). Рассмот-
рим его подробней в предположении, что d1 d2 d . В этом случае сум-
марная прибыль двух фирм равна:
W W1 W2 bx1(d x) bx2 (d x) bx(d x) .
Максимум прибыли достигается при выпуске:
xM d2 .
Равновесная цена будет следующей:
CM a b(x1M x2M ) a bd2 .
Для объединения фирм в форме картеля имеем:
x1C x2C d4 , W1C W2C bd8 2 .
42
Итак, из расчетов видно, для потребителя наиболее выгодна страте-
гия Стакельберга и наименее выгодной является стратегия монополии или картеля.
ЗАДАНИЕ
1. Для заданных значений параметров a , b , C01 , C02 , x1(0) , x2 (0) полу-
чить графики динамики изменения объемов выпуска фирм, динамики из-
менения прибылей фирм и динамики изменения цены для стратегии Кур-
но. Построить фазовый портрет, найти точку Курно, установившиеся зна-
чения прибылей фирм, объемов выпуска и установившуюся цену. Исход-
ные данные приведены в таблице 5.1.
Таблица 5.1
n/n |
a |
b |
C01 |
C02 |
x1(0) |
x1(0) |
1 |
15 |
1,1 |
1,1 |
1,4 |
7,0 |
4,5 |
|
|
|
|
|
|
|
2 |
17 |
1,2 |
1,05 |
1,5 |
5,0 |
4,7 |
|
|
|
|
|
|
|
3 |
20 |
1,1 |
1,2 |
1,3 |
6,5 |
4,2 |
|
|
|
|
|
|
|
4 |
19 |
1,3 |
1,2 |
1,1 |
5,6 |
5,5 |
|
|
|
|
|
|
|
5 |
21 |
1,3 |
1,2 |
1,3 |
6,5 |
4,7 |
|
|
|
|
|
|
|
6 |
13 |
1,4 |
1,1 |
1,7 |
7,6 |
4,5 |
|
|
|
|
|
|
|
7 |
24 |
1,0 |
1,2 |
1,5 |
4,5 |
6,2 |
|
|
|
|
|
|
|
8 |
20 |
1,3 |
1,08 |
1,2 |
6,2 |
6,2 |
|
|
|
|
|
|
|
9 |
15 |
1,0 |
1,4 |
1,2 |
6,0 |
4,5 |
|
|
|
|
|
|
|
10 |
20 |
1,1 |
1,2 |
1,3 |
6,2 |
4,8 |
|
|
|
|
|
|
|
11 |
15 |
1,6 |
1,3 |
1,4 |
5,0 |
6,5 |
|
|
|
|
|
|
|
12 |
25 |
1,35 |
1,4 |
1,3 |
6,5 |
5,2 |
|
|
|
|
|
|
|
43

6. Динамические модели фирмы
Рассмотрим модель производства n видов товаров в условиях рынка.
Вектор состояния x(k) представлен компонентами:
|
z1 (k ) |
|
|
|
||||
|
v (k ) |
|
|
|
||||
|
|
1 |
|
|
|
|
|
|
|
z2 |
(k ) |
|
|
||||
|
|
|
|
|
|
|
|
|
x(k ) |
v2 |
(k ) |
|
|
||||
|
|
|
|
|
, |
|
(6.1) |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
zn |
(k ) |
|
|
||||
|
v |
n |
(k ) |
|
|
|||
|
|
|
|
|
|
|
||
|
w(k ) |
|
|
|
||||
где zi (k) количество товаров i -го вида на рынке; vi (k) |
количество то- |
|||||||
|
|
|
прибыль. |
|
||||
варов i -го вида у потребителя, i 1, n ; w(k) |
|
Математическая модель динамики изменения количества товаров у потребителей и на рынке, а также прибыли может быть записана в следу-
ющем виде:
z1(k 1) (1 k11)z1(k) s1(k) u1(k) , v1(k 1) (1 k21)v1(k) s1(k) ,
z2 (k 1) (1 k12 )z2 (k) s2 (k) u2 (k) , v2 (k 1) (1 k22 )v2 (k) s2 (k) ,
………………….
………………….
zn (k 1) (1 k1n )zn (k) sn (k) un (k) ,
vn (k 1) (1 k2n )vn (k) sn (k) |
|
, |
|
w(k 1) w(k) c1s1(k) c2s2 (k) ... cnsn (k) k31z1(k) k32z2 (k) ... |
|
... k3n zn (k) c01u1(k) c02u2 (k) ... c0nun (k) , |
(6.2) |
44

где ui (k) количество товаров, выпускаемых за один такт, i 1, n , k1i ко-
эффициенты потерь; k2i |
коэффициенты потребления; k3i стоимость |
||
хранения единицы товаров; |
с0i себестоимости; si (k) количество про- |
||
|
|
|
|
данных товаров i -го вида в один такт, i 1, n . Функции продаж определя- |
ются по формулам:
s (k) n exp( c )(1 v (k)Y |
1)z (k) , |
(6.3) |
||||||
i |
i |
i |
i |
i |
i |
|
||
ni коэффициенты продаж; |
|
цены на товары, i |
|
; Yi |
потенциаль- |
|||
сi |
1, n |
ный спрос для i -го вида товара (объем рынка для i -го вида товара).
Модель (6.2), (6.3) может быть представлена в следующем векторно-
матричном виде:
x(k 1) A (x(k)) Bu(k) , |
(6.4) |
||||
где вектор (x(k)) следующий: |
|
|
|
|
|
|
|
z1 (k ) |
|
||
|
|
v1 (k ) |
|
||
|
|
|
|||
|
|
|
|
||
|
|
|
|
|
|
|
|
zn (k ) |
|
||
|
|
|
|||
|
|
vn (k ) |
|
||
( x(k )) |
|
w(k ) |
. |
||
|
|
|
|||
v (k )z (k ) |
|
||||
|
1 |
1 |
|
|
|
v2 (k )z2 (k ) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
v |
n |
(k )z |
n |
(k ) |
|
|
|
|
|
Матрица динамики A для данного объекта имеет вид:
|
a11 |
0 |
0 |
|
0 |
|
0 |
0 |
a1,2n 2 |
0 |
|
a21 |
a22 |
0 |
|
0 |
|
0 |
0 |
a2,2n 2 |
0 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
A |
0 |
0 |
0 |
a2n 1,2n 1 |
|
0 |
0 |
0 |
0 |
|
|
|
|||||||||
|
0 |
0 |
0 |
a |
2n,2n 1 |
a |
2n,2n |
0 |
0 |
0 |
|
|
|
|
|
|
|
|
|
||
|
0 |
a2n 1,3 |
a2n 1,2n 1 |
|
0 |
1 |
a2n 1,2n 2 |
a2n 1,2n 3 |
||
|
|
|
||||||||
a2n 1,1 |
|
где элементы матрицы определяются по формулам:
0 |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
a2n 1,3n 1 |
|
, |
|
|
|
a2n,3n 1 |
|
|
|
|
|
a2n 1,3n 1 |
|
|
|
|
45

a11 1 k11 n1 exp( c1 ),
a |
n1 exp( c1 ) , |
||
1,2n 2 |
|
Y1 |
|
|
|
a21 n1 exp( c1 ),
a22 1 k21,
a2,2n 2 n1 exp( c1 ) ,
Y1
a2n 1,2n 1 1 k1n nn exp( cn ),
a |
2n |
1,3n |
1 |
nn exp( cn ) |
, |
||
|
|
|
Yn |
|
|||
|
|
|
|
|
|
|
|
a2n,2n 1 nn exp( cn ) , |
|||||||
|
|
a2n,2n 1 k2n , |
|
||||
a |
|
|
|
nn exp( cn ) |
, |
||
2n,3n 1 |
|
||||||
|
|
|
|
Yn |
|
||
|
|
|
|
|
|
|
a2n 1,1 k31 c1n1 exp( c1 ),
a2n 1,3 k32 c2n2 exp( c2 ),
a2n 1,2n 1 k3n cnnn exp( cn ), ,
|
|
|
cini |
exp( ci ) |
|
|
|
|
|
|
|
|||
a |
2n 1, j |
, , |
j 2n 2,3n 1 , |
i 1,n . |
||||||||||
|
Yi |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Матрица B имеет вид: |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
1 |
0 |
0 |
|
0 |
0 |
|
|||||
|
|
|
0 |
0 |
0 |
|
0 |
0 |
|
|||||
|
|
|
|
|
||||||||||
|
|
|
0 |
1 |
0 |
|
0 |
0 |
|
|||||
|
|
|
0 |
0 |
0 |
|
0 |
0 |
|
|||||
|
|
|
|
|
||||||||||
|
B |
0 |
0 |
1 |
|
0 |
0 |
. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
0 |
0 |
0 |
|
0 |
1 |
|
|||||
|
|
|
0 |
0 |
0 |
|
0 |
0 |
|
|||||
|
|
|
|
|
||||||||||
|
|
|
c |
c |
c |
|
c |
c |
|
|||||
|
|
|
01 |
02 |
03 |
04 |
|
0n |
46
Рассмотрим модель производства двух видов товаров в условиях рын-
ка. Вектор состояния x(k) состоит из пяти компонент:
z1 (k ) |
x1 (k ) |
|
|
|
||||
v (k ) |
|
x |
|
(k ) |
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
x(k ) z2 (k ) |
x3 (k ) |
|
, |
|
(6.5) |
|||
|
|
|
|
|
|
|
|
|
v2 (k ) |
x4 (k ) |
|
|
|
|
|||
w(k ) |
|
x (k ) |
|
|
|
|
||
|
|
|
5 |
|
|
|
|
|
где z1(k) , z2 (k) количество товаров 1-го и 2-го вида на рынке; v1(k) , |
v2 (k) |
|||||||
количество товаров 1-го и 2-го вида у потребителя, w(k) |
прибыль. |
|
Математическая модель динамики изменения количества товаров у потребителей и на рынке, а также прибыли может быть записана в следу-
ющем виде:
z1(k 1) (1 k11)z1(k) s1(k) u1(k) , |
|
v1(k 1) (1 k21)v1(k) s1(k) , |
|
z2 (k 1) (1 k12 )z2 (k) s2 (k) u2 (k) , |
|
v2 (k 1) (1 k22 )v2 (k) s2 (k) , |
|
w(k 1) w(k) c1s1(k) c2s2 (k) k31z1(k) k32 z2 (k) |
|
c01u1(k) c02u2 (k) , |
(6.6) |
где u1(k) , u2 (k) количество товаров, выпускаемых за один такт; k11 , k12
коэффициенты потерь; k21 , k22 коэффициенты потребления; k31 , k32 сто-
имость хранения единицы товаров; с01 , с02 себестоимости; s1 (k) , s2 (k)
количество проданных товаров 1-го и 2-го вида в один такт (функции про-
даж). Формулы для s1 (k) , s2 (k) имеют вид:
s (k) n |
exp( c )(1 v (k)Y 1)z (k) , |
(6.7) |
||||
1 |
1 |
1 |
1 |
1 |
1 |
|
s (k) n |
exp( c |
)(1 v |
(k)Y |
1)z (k) , |
(6.8) |
|
2 |
2 |
2 |
2 |
2 |
2 |
|
47
n1 , n2 коэффициенты продаж; с1 , |
с2 цены на товары; Y1 , |
Y2 потенци- |
||
альный спрос на товар 1-го вида и 2-го вида. |
|
|||
В векторно-матричном виде модель следующая: |
|
|||
x(k 1) A (x(k)) Bu(k) , x(0) x0 , |
(6.9) |
|||
В (6.9) вектор (x(k)) представляется в виде: |
|
|||
|
x1(k ) |
|
|
|
|
x2 (k ) |
|
|
|
|
|
|
||
|
|
|
||
|
x3 (k ) |
|
|
|
|
|
|
|
|
( x(k )) |
x4 (k ) |
. |
(6.10) |
|
|
x5 (k ) |
|
|
|
|
|
|
|
|
x1 |
(k )x2 (k ) |
|
||
x |
(k )x |
(k ) |
|
|
3 |
4 |
|
|
|
Матрица динамики А для данного объекта имеет вид:
|
|
k11 |
n1 exp( c1) |
|||
|
1 |
|||||
|
|
|
|
|
|
|
|
|
|
|
n1 exp( c1 ) |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
A |
|
|
|
|
0 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
31 |
c n exp( c ) |
|||
|
|
|
|
1 1 |
1 |
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
0 |
0 |
|
|
n1 exp( c1 ) |
|
||
|
|
|
|
|
|
|
Y1 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|||
1 k |
|
|
|
0 |
|
0 |
0 |
|
n1 exp( c1 ) |
|
|||
|
21 |
|
|
|
|
|
|
|
|
|
Y1 |
||
|
|
|
|
|
|
|
|
|
|
|
|||
0 |
1 k12 n2 exp( c2 ) |
0 |
0 |
0 |
|
|
|||||||
0 |
|
|
n2 exp( c2 ) |
1 k22 |
0 |
0 |
|
|
|||||
0 |
k |
|
c n |
exp( c ) |
0 |
1 |
|
c1n1 exp( c1 ) |
|||||
32 |
|
||||||||||||
|
|
2 |
2 |
2 |
|
|
|
|
|
Y1 |
|||
|
|
|
|
|
|
|
|
|
|
|
Матрица В и вектор управления следующие:
|
1 |
0 |
|
|
|
|
0 |
0 |
|
|
|
|
|
u1 |
(k ) |
||
B 0 |
|
, |
|||
1 |
u(k) |
. |
|||
|
0 |
0 |
|
u2 |
(k ) |
|
|
|
|
||
|
c |
c |
|
|
|
|
01 |
02 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n2 exp( c2 ) |
|
|
|
|
|
|
. |
|||
|
|
Y2 |
|
|
|
|
|
|
|
||
|
n2 exp( c2 ) |
|
|
||
Y2 |
|
|
|||
|
|
|
|||
|
|
|
|
||
c2n2 exp( c2 ) |
|||||
|
|
Y2 |
|
|
|
|
|
|
|
48
ЗАДАНИЕ
1. Для модели фирмы, производящей два вида товаров (6.6) (6.8) вы-
полнить моделирование для следующих значений параметров:
u1 60, u2 65 количество товаров, выпускаемых фирмой за один такт; n1 1,95, n2 1,8 коэффициенты продаж; c1 2,5 у.е., c2 1,5 у.е. цены на товары; c0,1 1,0 у.е., c0,2 0,9 у.е. себестоимости; Y1 Y2 1000 потенци-
альный спрос (объем рынка); k1,1 0,15, k1,2 0,13 коэффициенты потерь; k2,1 0,1, k2,2 0,055 коэффициенты потребления; k3,1 0,002 у.е., k3,2 0,001
у.е. стоимости хранения единицы товара за один день.
Моделирование выполнить на интервале времени от 0 до 140 (один такт соответствует 1 дню) для следующих начальных условий:
z1(0) 150, z2 (0) 300, v1(0) 250, v2 (0) 170, w(0) w0 у.е.
Построить графики процессов (величина w0 приведена в таблице 6.1).
2. Исследовать влияние различных стратегий управления фирмой на полученную прибыль.
Стратегия 1. Увеличить цену 2-го вида товара c2 до величины 2,3 у.е.
Привести в отчет графики изменения прибыли. Определить прибыль в по-
следний день исследуемого периода ( w140 ). Оценить возможность реальной реализации этой стратегии. Сделать выводы.
Стратегия 2. Увеличить коэффициент продаж n2 до величины 3,2 (увеличение этого коэффициента можно осуществить, реализовав реклам-
ную компанию). В модели учесть затраты на рекламу в 2у.е. в течении первых 10 дней. Затем этот коэффициент должен уменьшаться по линей-
ному закону в течение 60 дней до первоначальной величины n2 1,8 . Затем опять провести рекламную компанию в течение 10 дней.
49
Построить графики изменения прибыли. Определить прибыль в по-
следний день исследуемого периода ( w140 ). Сделать выводы.
Стратегия 3. Увеличить потенциальный спрос (объем рынков для 1-го и 2-го вида товаров). В модели учесть затраты на расширение рынка в 8у.е.
в течении первых 60 дней. По окончании этого периода значения Y1 и Y2
принять равными 2000 (увеличение этих параметров осуществляется по-
средством расширения рынка в течении первых 60 дней, например, создав новые торговые точки в новом регионе).
Построить графики изменения прибыли. Сделать выводы.
3. Применить метод покоординатного спуска для максимизации кри-
терия J (u1, u2 ) w140 (прибыли фирмы в последний день), применив метод деления шага пополам. Начальное значение шага принять равным 10. оп-
тимизацию осуществить сначала по переменной u2 , затем по переменной u1 .
Промежуточные результаты оформить в виде таблицы. Привести в отчете оптимальные значения объемов производства и оптимальное значе-
ние прибыли.
4. Найти оптимальные значения объемов производства и прибыли с учетом ограничений (величина uмах приведена в таблице 6.2):
u1 u2 uмах .
|
|
|
|
|
|
Таблица 6.1 |
|
|
|
|
|
|
|
|
|
n/n |
1 |
2 |
3 |
4 |
5 |
|
6 |
|
|
|
|
|
|
|
|
w0 |
10 |
25 |
30 |
0 |
45 |
|
50 |
|
|
|
|
|
|
|
|
50