Имитационное статистическое моделирование рынка одного товара с оптимальной детерминированной стратегией поставки товара в условиях стохастичности спроса Текст научной статьи по специальности «Экономика и бизнес»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2012 Управление, вычислительная техника и информатика № 1(18)

УДК 519.865

В.В. Поддубный, О.В. Романович

ИМИТАЦИОННОЕ СТАТИСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РЫНКА ОДНОГО ТОВАРА С ОПТИМАЛЬНОЙ ДЕТЕРМИНИРОВАННОЙ СТРАТЕГИЕЙ ПОСТАВКИ ТОВАРА В УСЛОВИЯХ СТОХАСТИЧНОСТИ СПРОСА

Рассматривается имитационная модель инерционного рынка одного товара в условиях запаздывания поставок и флуктуаций покупательского спроса при оптимальной детерминированной стратегии поставки товара на рынок. Математическое описание рынка представляется рестриктивной (подчиняющейся ограничениям типа неравенств) динамической моделью с запаздывающим управлением, учитывающей случайные флуктуации покупательского спроса. Показано, что оптимальная в смысле максимума прибыли продавца стратегия поставки товара на рынок определяется условиями состояния рынка (товарного дефицита, затоваривания рынка, динамического равновесия рынка). Построен алгоритм нахождения оптимальной детерминированной стратегии поставки товара на рынок на основе упреждающей детерминированной модели рынка. Проведено имитационное моделирование и статистический анализ функционирования рынка с одновременным использованием стохастической и упреждающей детерминированной моделей.

Ключевые слова: динамическая система, ограничения типа неравенств, запаздывающее управление, оптимизация, рынок одного товара, флуктуации спроса, статистическое имитационное моделирование.

1. Постановка задачи

В работе [1] рынок одного товара, функционирующий в условиях запаздывания поставок товара на рынок, рассматривался как детерминированная инерционная нелинейная управляемая динамическая система, подчиняющаяся ограничениям типа неравенств. Получено точное решение задачи оптимального управления поставкой товара на рынок при детерминированном спросе на товар в указанных условиях. Было показано, что оптимальная стратегия поставки товара на рынок может быть реализована с помощью детерминированной математической модели, позволяющей точно прогнозировать состояние рынка вперёд на время запаздывания. Проведено имитационное моделирование инерционного рынка одного товара при оптимальном управлении поставкой товара на рынок в условиях запаздывания поставок. В работе [2] описывается математическая модель рынка одного товара, в которую включено влияние флуктуаций покупательского спроса на переменные, описывающие функционирование рынка в условиях запаздывания поставки товара на рынок при той же детерминированной (но уже не оптимальной) стратегии поставки товара. В настоящей работе проводится статистическое имитационное моделирование рынка с флуктуирующим покупательским спросом и исследуются его статистические характеристики.

Рассмотрим модель рынка одного товара, функционирующую в дискретном времени ґ = 0,1,2,... . Пусть Р(ґ) - цена товара в момент времени ґ, QO(t) - остаток товара на рынке. Стратегия заказа товара определяется переменной QZ(t - т) -

объём товара, заказываемого в момент времени t - т для поставки на рынок в момент t (предполагается запаздывание т в поставках товара). Пусть в момент t дискретного времени спрос QD(t) на товар при цене P(t) имеет вид простейшей линейной зависимости с аддитивным членом ^(t), характеризующим случайные флуктуации спроса:

Qd (t) = Qm -aP(t) + £(t), (1)

где Qm > О и a > О - заданные константы.

Будем считать, что ^(t) - стационарный некоррелированный случайный процесс с нулевым средним значением и произвольным стандартным отклонением с. Неотрицательность спроса обеспечим нелинейной рестриктивной зависимостью спроса от цены при любом значении с путём усечения процесса £(t):

QD (t ) = {Qm - aP (t ) + W ^ если Qm - aP (t ) + W )> 0, (2)

Q ltj \ О, если Qm - aP (t) + 5(t)< О. (2)

Объём товара Q(t), предлагаемого к продаже в момент времени t, складывается из остатка товара в объёме Q0(t) и товара в объёме QZ(t - т), заказанного продавцом в момент времени t - т (с учётом запаздывания поставки) для поставки его на рынок к моменту времени t:

Q(t) = Q0 (t) + Qz (t-т).

Тогда объём продаж QS(t) на интервале t дискретного времени равен

QS (t) = min(Qd (t), Q° (t) + Qz (t -т)), (3)

а остатки товара удовлетворяют рекуррентному соотношению

Q0 (t+1) = Q0 (t)+Qz (t-т)-QS (t). (4)

Будем считать целевой функцией рыночного механизма ценообразования прибыль продавца J(t) в момент времени t. В рассматриваемых условиях прибыль продавца равна разности между выручкой от продажи товара и затратами на его приобретение и хранение. Если P1 - цена закупки товара (на оптовом рынке или у производителя), P2 - цена хранения единицы товара, не проданного на предыдущем интервале дискретного времени, то прибыль продавца в момент времени t составит величину

J(t) = QS (t)P(t)-Qz (t)P -Q° (t)P2 -R2(P(t)-P(t-1))2. (5)

Последнее слагаемое (с коэффициентом R > О) «штрафует» продавца за изменение цены товара и определяет инерционность рынка - за резкое повышение цены могут последовать санкции законодательного характера, за резкое снижение цены - «санкции» конкурентов, выражающиеся в нанесении ущерба продавцу в размере, эквивалентном этой штрафной функции.

Решим задачу нахождения такой цены товара P(t) (устанавливаемой рынком) и такой дополнительной поставки QZ(t - т) товара на рынок (устанавливаемой продавцом), чтобы прибыль продавца при заданной линии спроса на t-ом интервале дискретного времени была максимальной:

J (t) ^ sup . (6)

P(t ),QZ (t-т)

При решении этой задачи должны выполняться ограничения на величину цены товара P(t):

P1 < Pmin < P(t) < Pmax(t) = (Qm + Ф)) / a и на величину дополнительного заказа товара Qz(t - т) > О.

2. Стохастическая и детерминированная модели.

Оптимизация цены и поставки товара

Принцип решения поставленной задачи (6) без учета возможных флуктуаций спроса описан в работе [1]. При наличии флуктуаций спроса ^(t) схема решения остаётся аналогичной и впервые представлена в работе [2]:

1) Находится оптимальная (обеспечивающая максимум прибыли продавца (5)) цена P(t) товара при фиксированных значениях предыдущей цены P(t - 1) товара и объёма предложения Q(t) (следовательно, при фиксированном объёме заказа QZ(t - т)):

J (t) ^ max .

P(t )IQ(t)

2) Находится оптимальное значение объёма Q(t) товара, предлагаемого на рынке, и заказа Qz(t - т) товара в зависимости от его остатка Q0(t) на рынке.

Однако имеется важное различие между стохастической и детерминированной задачами оптимизации. В детерминированной постановке задачи в каждый момент дискретного времени t на рынок поступает оптимальный объём Qz(t - т) товара, заказанного за т шагов до момента t в условиях полной предсказуемости поведения рынка в течение этих шагов. А в стохастической постановке рынок идёт по стохастической траектории, так что предсказать точно будущий спрос на товар и найти точное значение оптимального объёма Qz(t - т) поставки товара на рынок в момент t - т невозможно (будущие значения флуктуаций спроса £(t) ещё не известны).

Поэтому решать стохастическую задачу приходится с использованием двух моделей. Одна модель (стохастическая) является имитационной. Она моделирует стохастическую динамику рынка при известном на каждом шаге t объёме поставки Qz(t - т). Другая модель (детерминированная) - скользящая упреждающая. Она запускается на каждом текущем шаге t на т шагов вперёд с начальными условиями, выражающими состояние стохастической модели в этот момент времени, для предсказания будущего состояния рынка на момент t + т. По предсказанному состоянию производится расчёт требуемой к этому моменту времени оптимальной детерминированной поставки товара Qz(t + т), которая и принимается к реализации в имитационной модели (или на настоящем рынке).

На рис. 1 приведена схема взаимодействия стохастической и детерминированной моделей рынка.

При решении как стохастической, так и детерминированной задачи (с учётом её рестриктивности в силу соотношения (3)) следует выделить области, соответствующие дефициту товара на рынке (область 1, в которой Q(t) < QD(t)), затовариванию рынка (область 2, в которой Q(t) > QD(t)) и балансу спроса и предложения (область 3, область динамического равновесия, в которой Q(t) = QD(t)). Соотношения, получаемые для стохастической модели, являются более общими. В детерминированной модели во всех формулах следует положить £(t) = О.

Детерминированная модель

г 1й&т

Стохастическая модель функционирования рынка Рис. 1. Взаимодействие детерминированной и стохастической моделей рынка

1) В области товарного дефицита Q(t) < QD(t), и в соответствии с соотноше-

нием (3) имеем Q (t) = Q(t), так что

-\0

R

J(t) = Q(t)P(t)-Q(t)P + Qw (t)(P -P2)--(P(t)-P(t-1))2 ^ sup . (7)

2 P(t)IQ(t)

Это вогнутая квадратичная функция переменной P(t), достигающая условного максимума при

P (t) = P (t -1) + = P(1)(t). (8)

R

Как видим, P(1)(t) растет с ростом Q(t) по линейному закону. Выражение (8) справедливо не при любом Q(t), а лишь при Q(t), удовлетворяющем условию Q(t) < Q°(t) принадлежности к области 1. Это условие с учетом (1) и (8) имеет вид

Q (t )< Qm +P(t)-»P»l(t )= Qm +P(t)-aP (t -1)-Rq (t),

Q (t )< R (Q„ + Kt )-aP (--1)) = Q,,)

a + R

Поскольку в силу (2) QD (t) = Qm + £ (t) - aP(1) (t) > 0, возможная (допустимая) флуктуация ^(t), не приводящая к отрицательности спроса, подчиняется неравенству

£(t )>-(Qm - aP (t -1)), (10)

что обеспечивает выполнение неравенства Q(1) (t) > 0. Таким образом, в области 1 P(t) линейно растёт с ростом Q(t) от значения

P(1)(t )Q(t )=о = P (t -1)

до значения

откуда

(9)

э0)(-) | ... = Qm +^(t)+RP(t-1) = Pm1 (t)

P(1)(t )|

q(-)=Q(1)(t)

a + R

причем условие принадлежности Q(t) к области 1 выражается неравенством (9) с

учётом (10), то есть 0 < Q(t) < б(1)(/). После подстановки Р(0 = Р(1)(0 в выражение (7) для ^) имеем

3 (t )=^2т+(р ( -1}-р )Q ^)+Q0 (t )(р - р2 )=3 (1)^) •

Дальнейшая оптимизация J(1)(t) по Q(t) производится только в детерминированной модели, не содержащей в течение т предыдущих шагов флуктуации спроса 4(0, и нужна для определения оптимального объёма поставки QZ(t - т) товара продавцом.

Видно, что ./^(О монотонно растет с ростом Q(t) по линейно-квадратичному закону, достигая максимального значения на границе области при Q(t) = Q(1)(t):

Jmax (t )=j (і)(-)q

0

Q(t)=Q(1)(t) ■

Если остаток товара QO(t) от продаж предыдущего интервала дискретного времени не превышает величину Q(1)(t), то дополнительный заказ и поставка товара на рынок в объеме QZ(t - т) = Q(1)(t) - QO(t) (в частности, QZ(t - т) = 0 при

QO(t) = Q(1)(t)) обеспечивает максимум прибыли. Иначе, при QO(t) > Q(1)(t), следу-

ет искать решение задачи оптимизации прибыли в областях 2 или 3.

2) В области затоваривания рынка Q(t) > QD(t), и в соответствии с выражением (3) имеем QS(t) = QD(t), так что с учетом (1)

J(t) = (Qm + 4(t)-aP(t))P(t)-Q(t)P + QO (t)(P -P2)-

-R(P(t)-P(t-!))2 ^ sup . a1

2 P(t )IQ(t)

Это вогнутая квадратичная функция переменной P(t) с точкой максимума

P(t) = Qm + £(t) + RP(t-1) = P(2)(t). (12)

V 7 2a + R W

При выполнении условия (10) эта величина неотрицательна и не меньше P(1)(t). Как видим, P(2)(t) не зависит от Q(t) (остаётся постоянной при любом Q(t) в этой области). Выражение (12) справедливо лишь при условии, что Q(t) > QD(t), то есть при условии

,(2)^)= (a + R XQm + 4(t)-aRP (t -1)

Q (t )> Qm + %(t)- aP (t ) =

2a + R

/ч р(Qm + 4(ґ)-аР(ґ-1)) + а(Qm + 4(ґ)) (2), ч

откуда Q(ґ)> У^т )------------V----Ц----^ >’ =Q(2)(t). (13)

2а + Я

Последнее неравенство определяет условие принадлежности Q(t) к области 2. Нетрудно показать, что Q(2)(t) > Q(1)(t). В области 2 Q(t) > Q(2)(t). После подста-

новки в З(ґ) Р(ґ) = Р(2)(0, не зависящего от Q(t), имеем

З (ґ) = ( + 4 (ґ) - аР(2) (ґ))Р(2) (ґ) - Q (ґ)Р + Qo (ґ) (Р - Р2) -

- Я2 (Р(2) (ґ)-Р ( -1))2 = ) (2)(ґ).

Дальнейшая оптимизация ./2)(0 по Q(t) производится только в детерминированной модели, не содержащей в течение т предыдущих шагов флуктуации спроса 4(0, и нужна для определения оптимального объёма поставки QZ(t - т) продавцом.

Как видим, J(2)(t) монотонно убывает с ростом Q(t) по линейному закону, так что достигает в этой области наибольшего значения при Q(t) = Q(2)(t):

jmax (t )=J(2)(t )Iq(, ,=q,)(, ).

Если QO(t) < Q(2)(t), то дополнительный заказ и поставка товара на рынок в объеме QZ(t - т) = Q(2)(t) - QO(t) обеспечивает получение этого максимума прибыли. Если же QO(t) > Q(2)(t), то QZ(t - т) = 0 и достигается лишь значение прибыли

J(4(t v нл, )< jmx (').

3) В области баланса спроса и предложения (то есть в области динамического равновесия рынка) Q(t) = QD(t) и в соответствии с выражением (3) имеем, как и в области 2, объём продаж, равный спросу, то есть QS(t) = QD(t), и прибыль J(t) в виде (10). Но Q(t) = Qm + 4(t) - aP(t), откуда

P (, ) = Qm +£(t )-Q (t) = p(3) (,). (14)

a

Границами области 3 по Q(t) являются точки Q(1)(t) и Q(2)(t):

Q(1)(t )< Q (t )< Q(2) (t).

Как видим из (13), в этой области P(3)(t) линейно убывает с ростом Q(t) от

P(3)(t) | = Qm + RP (t - 1) = P(1) (t)

P (t)IQ(t)=Q(1)(t)= a + R ~ Pmax (t)

до р(3)(, ) | = Qm +£(t )+RP (t-1)= p(2) (,)

до P (t Q(t)=Q(2)(,)= 2a + R ~ P (t).

После подстановки P(t) = P(3)(t) в J(t) для этой области имеем

J (t) = Q (t) Qm +£(t )-Q (t) - Q (t) р + QO (t) (P - P2 ) -

a

R ^ Qm +£(t )-Q (t)

Дальнейшая оптимизация J'3)(t) по Q(t) производится только в детерминированной модели, не содержащей в течение т предыдущих шагов флуктуации спроса

4(t), и нужна для определения оптимального объёма поставки QZ(t - т). J(3)(t) - вогнутая линейно-квадратичная функция переменной Q(t). В точке

Q(t) = R(Qm +£(t)-ap(t-l)) + a(Qm +4(t)-ap) = Q(3) (,)

^ ’ 2a + R У ’

прибыль J(3)(t) достигает максимального значения. Нетрудно показать, что Q(1)(t) < Q(3)(t) < Q(2)(t), то есть точка максимума Q(3)(t) лежит в области 3. Максимальное значение прибыли в области 3 (при QO(t) < Q(3)(t))

jmax (, )=jmc )q„ )=q(3)(, ).

Это значение является и глобально максимальным, потенциально возможным при QO(t) < Q(3)(t). Если же Q(3)(t) < QO(t) < Q(2)(t), то глобально максимальное значение прибыли не может быть достигнуто, и достигается только условномаксимальное значение (при фиксированном QO(t)) внутри области 3, лежащее

между jmax (,) и j(2)(t).

Qm+^(?) Q(t) -p(t-1)^j = j(3)(t).

3. Статистический анализ стохастической динамики рынка

Статистический анализ функционирования имитационной модели рынка в условиях стохастических флуктуаций спроса проводился в различных ситуациях, в том числе в ситуации, когда сначала, до некоторого момента времени t = 0 включительно, рынок находился в состоянии равновесия и флуктуации спроса отсутствовали. Равновесная цена товара P* = ^т + aP1) / (2а), спрос 0°* = Qm - aP' равнялся предложению, запасы товара на складе отсутствовали (<2° = 0). Таким образом, в момент времени t = 0 цена товара P0 оставалась равной P , остаток товара Q0 = 0. Начиная со следующего момента, спрос начинал флуктуировать, вследствие чего рынок выводился из состояния равновесия. Флуктуации спроса моделировались последовательностью независимых нормально распределённых случайных величин с нулевым средним и одинаковой дисперсией с2. При этом для обеспечения неотрицательности спроса производилось его усечение в соответствии с рестриктивным преобразованием (2). Расчёты проводились при Я = 50, Qm = 4, a = 0,4, Т = 2000, Po = P, Pl = 3, P2 = 0,1, Pmln = Pl + P2, Pmax = Qm / а, P = 6,5, Q0 = 0, т = 7, с = 0,1. Эксперимент повторялся N = 2000 раз. При каждой п-й реализации имитационного моделирования рынка (прогонке модели) рекуррентно уточнялись средние значения М((), дисперсии D(t) = с2(0 и функции автокорреляции К(^ t+т) = с(0с^+т)г(^ t+т) всех переменных модели по схеме:

Мп+1 (t)=~П~Л Мп ^)+~Ц-хп+1, М1 = Х^ п = 1 N - 1,

п +1 п + 1

^+1 (t) = Dn ()+ “ (Хп+1 (t)-Mn+1 (t))2, D1 (t) = (x1 (t)-М1 (t))2 = ^ п = 1 N - 1,

п п

Кп+1 (t,t + Т) = —Кп (^t + Т) +1 (п+1 (t)-Мп+1 (t))Хп+1 (t + Т)-Мп+ (t + Т)) ,

п п

К1 (, t + т) = 0, п = 1, N -1.

На рис. 2 - 7 изображена реализация стохастической динамики основных переменных состояния рынка вместе с их математическими ожиданиями М и границами полос разброса значений (шириной 2с(0). Видно, что средние значения всех параметров рынка в условиях флуктуаций спроса постепенно отклоняются от начальных равновесных значений параметров детерминированного (невозмущённого) рынка и асимптотически приближаются к значениям, характеризующим статистическое равновесие. Это смещение объясняется нелинейным (именно, рестриктивным, связанным с ограничениями вида неравенств (2), (3)) характером рассматриваемой модели рынка.

Как видно на рис. 2 и 3, хаотические флуктуации спроса приводят к сравнительно медленным флуктуациям оптимальной цены и предложения товара, определяемого оптимальной детерминированной стратегией поставки товара на рынок. Остальные переменные флуктуируют достаточно быстро.

На рис. 8 представлено поведение нормированной автокорреляционной функции цены + т), полученное по N = 2000 прогонкам модели рынка в одних и тех же условиях. Верхняя и нижняя горизонтальные плоскости, показанные на рисунке, указывают границы 95 %-го интервала статистической незначимости функции автокорреляции, то есть статистической значимости нулевой гипотезы

6.8

6.6

6.4

6.2

5.8 5.6 5.4

О 500 1000 1500 ^

Рис. 2. Динамика цены товара

1.75

1.7 1.65

1.6 1.55

1.5 1.45

14 1 35

' 0 500 1000 1500 I

Рис. 3. Стратегия заказа товара

2

1.9

1.8 1.7

1.6 1.5 1.4 1.3

0 500 1000 1500 ь

Рис. 4. Покупательский спрос

О 500 1000 1500

Рис. 5. Остаток непроданного товара

500 1000 1500

Рис. 6. Динамика продаж

500 1000 1500

Рис. 7. Динамика прибыли

Н0: г(/, /+т) = 0 о некоррелированности флуктуаций. Они определяются /-критерием Стьюдента и равны ± \/ф + (Ы _ 2))^_2 0 975 = ± 0,0438, где (ы_2 0 975 -

квантиль уровня 0,975 распределения Стьюдента с N - 2 степенями свободы. Как видно на рис. 8, при критическом уровне значимости 5 %-й корреляции значений цены товара статистически значимо отличаются от нуля для всех значений параметра т сдвига в интервале от т = 1 до т = 200, не менее. Аналогично ведёт себя и функция автокорреляции ге2(/,/ + т) предложения (заказа) товара. Все остальные переменные имеют слабо коррелированные во времени значения. Наименее авто-коррелированы остатки 0° непроданного товара, функция автокорреляции гд°(/,/ + т) которых приведена на рис. 9. Поверхности функций автокорреляции объёма О продаж, спроса 0° и прибыли 3 имеют вид, промежуточный между поверхностями, изображёнными на рис. 8 и 9, но ближе к рис. 9. Почти всюду их автокорреляции практически статистически незначимы.

12-.

1

0 8 гР 0 6'

04020

-0 2 2000

1000

Рис. 8. Нормированная автокорреляционная функция цены и границы 95 %-го интервала её незначимости

1.2-г" 1 -■■■■■■'

0.8---""

1000 200 / 0 300 Т

Рис. 9. Нормированная автокорреляционная функция остатков и границы 95 %-го интервала её незначимости

Заключение

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

ЛИТЕРАТУРА

1. Поддубный В.В., Романович О.В. Рестриктивная динамическая модель инерционного рынка одного товара с оптимальной поставкой товара на рынок в условиях запаздывания // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. 4(17). С. 16-24.

2. Поддубный В.В., Романович О.В. Имитационная модель рынка одного товара с оптимальной детерминированной стратегией поставки товара в условиях стохастичности спроса // Информационные технологии и математическое моделирование (ИТММ-2011): Материалы Х Всероссийской научно-практической конференции с международным участием (25-26 ноября 2011 г.). Томск: Изд-во Том. ун-та, 2011. Ч. 2. - С. 47-53.

Поддубный Василий Васильевич Романович Ольга Владимировна Томский государственный университет

E-mail: vvpoddubny@gmail.com; njkm@ngs.ru Поступила в редакцию 1 декабря 2011 г.

Poddubny Vasily V., Romanovich Olga V. (Tomsk State University). Simulation statistical modeling of the market of one goods with the optimal deterministic strategy of the delivery of goods under condition of stochastic demand.

Keywords: dynamic system, restrictions of an inequality type, time-lag control, optimization, market of one goods, casual fluctuations of demand, statistical simulation modeling.

The simulation model of the inertial market of one goods is considered under conditions of delivery time-lag and casual fluctuations of the purchase requirement under the optimal deterministic strategy of the supply of goods to the market. The mathematical description of the market is given by the restrictive (because of restrictions of inequalities type) dynamic model with the time-lag of control taking into account casual fluctuations of the purchase requirement. It was shown that the optimal (in the sense of the maximum profit of the seller) strategy of the goods delivery to the market is determining by the market conditions (commodity deficit, overstocking of the market, or the dynamic balance state of the market). The algorithm of the finding of optimum deterministic strategy of the delivery of goods to the market is constructed with using of the forecasting deterministic model of the market. The simulation modeling and statistical analysis of the market dynamics are organized with simultaneous use of the stochastic model and of the forecasting deterministic model. The statistical simulation of the model behavior is performed; average values, variances, and auto-correlation functions for all variables of the model are investigated. It is shown that in view of non-linearity of the model the casual fluctuations of demand shift the statistical equilibrium of the market from its deterministic equilibrium.