Асимптотический анализ двухкритериальной модели реальных инвестиций на основе Z-преобразования

Автор: Победаш Павел Николаевич

Журнал: Сибирский аэрокосмический журнал @vestnik-sibsau

Рубрика: Математика, механика, информатика

Статья в выпуске: 2 (23), 2009 года.

Бесплатный доступ

Предлагается приложение z-преобразования к асимптотическому анализу линейных многошаговых задач оптимального управления на примере двухкритериальной модели оптимизации реальных инвестиций предприятия с неограниченным количеством видов производимой продукции и/или спросом на нее на бесконечном горизонте планирования.

Z-преобразование, двухкритериальная модель оптимального управления, планирование инвестиций

Короткий адрес: https://sciup.org/148175894

IDR: 148175894

Текст научной статьи Асимптотический анализ двухкритериальной модели реальных инвестиций на основе Z-преобразования

В данной работе рассматривается задача, обобщаю- тие имеет капитал K0 . При этом государственный орган щая задачу из работы [1] на случай двух критериев, ко- (ГО) для реализации инвестиционного проекта (ИП) торую сформулируем следующим образом. Предприя- выделяет инвестиции не более величины I0 на приоб- ретение активных основных производственных фондов (ОПФ) n видов. Необходимо определить стоимость (количество) всех приобретаемых в моменты t = 1, _., T ОПФ каждого вида, при которых дисконтированные суммы собственных средств предприятия и его налоговых поступлений в ГО максимальны за все время действия ИП. При этом выполнены следующие основные предпосылки:

  • 1)    учитываются налоги, составляющие большую часть затрат предприятия: налог на добавленную стоимость (НДС), налог на прибыль (НП), налог на имущество (НИ), единый социальный налог (ЕСН) и отчисления в фонд оплаты труда (ФОТ);

  • 2)    предприятие имеет достаточные запасы сырья;

  • 3)    срок T действия ИП меньше сроков Tk службы единицы ОПФ каждого типа: T T k ( к = 1,..., n );

  • 4)    на ОПФ каждого вида производится лишь один тип продукции. С учетом приведенных предпосылок сформулированная выше задача имеет вид двухкритериальной многошаговой задачи линейного программирования (МЗЛП), которую, согласно работе [2], назовем моделью B 1:

x k ( t + 1) = x k ( t ) + u k ( t )

( к = 1,..., n ; t = 0,..., T - 1),

n

X n +1 ( t + 1) = - S X k ( t ) / T k + X n +1 ( t ) + к =1 n

+ S U k ( t ) ( t = 0,..., T - 1)

k =1

xn+2 (t + D = -a2 xn+1 (t) + xn+2 (t) - n

  • - S u k ( t ) + u 2 n +1( t ) + u 2 n +2( t ) ( t = 0) ,            (1)

k =1

nL X. ( t )

xn+2 (t + D = a3 S --6xn+1 (t) + xn+2 (t ) - k=1 Tk nn

  • - S uk ( t ) + y S u n + k ( t ) ( t = U-, T - 1);

k =1                k =1

X k (0) = 0( k = 1,..., n + 2);

X n +2 ( t ) 0 ( t = 1,..., T ), Xk ( t )             , X

S v    a 2 Xn+1 (t)+ k=1   Tk

+ (1 - р£> + k ( t ) > 0 ( t = 1,..., t - 1);

k =1

U n + k ( t ) s k X k ( t ) ( k = 1,..., n ; t = 1,..., t - 1),

U 2 n +1 (0) I 0 , U 2 n +2 (0) K 0 , uk ( t ) > 0 ( k = 1,..., n ; t = 0,..., T - 1), U n + k ( t ) 0 ( k = 1,..., n ; t = 1,..., T - 1),

U 2 n +1 (0) 0, u 2 n +2 (0) 0,

J = { J 1 , J 2} ^ max , где J 1 =- u 2 n +1(0) - u 2 n +2(0) +

T -1 S t =1

a 3 S ^^ - 6 Xn +1 ( t ) + y S u n + k ( t )

k =1 T k _____________________ k =1 __________ J +

(1 + r ) t

5 Xn +1 ( T )

(1 + r ) T -1

T - 1

J2 =S t = T2

- Из] n: x Д( t ) + 6 X i( t ) + PlLU n + k ( t ) k = 1 T k           n + 1            k = 1

(1 + Г ) t соот-

ветственно дисконтированная сумма собственных средств предприятия и налоговых поступлений в ГО; Uk(t)(t = 0,..., T -1), Un+k(t) (k = 1,..., n; t = 1,..., T -1), u2n+1(0) и u2n+2(0) - стоимость приобретаемых ОПФ, вы- ручка от реализации продукции k -го типа, внешние и внутренние инвестиции соответственно; Xk(t) (k = 1,..., n), Xn+1 (t), Xn+2( t) (t = 0,..., T) - соответственно накопленная стоимость всех ОПФ k-го типа, остаточная стоимость всех ОПФ, текущие денежные средства предприятия в момент t; Vk, Tk, ck и Pk - соответственно производительность, срок службы, стоимость единицы ОПФ и стоимость единицы продукции k-го типа; 10, K0 - суммы внешних и внутренних инвестиций, выделяемых на весь срок действия ИП; а1, а2, а3, а4 - ставки НДС, НИ, НП и ЕСН соответственно (НДС включается в цену продукции, поэтому можно считать, что а, = 0); Р - доля выручки от реализации, выделяемая на ФОТ; 6 = (1 - а3)а2, 5 k = PkVk / ck (k = 1,..., n), Y = (1 - а3)(1 - в), p = (1 - в)а3 + а4в, r - ставка доходности ИП; 5 (0 < 5 < 1) - доля остаточной стоимости всех ОПФ на момент t = T от ее балансовой стоимости, определяемая в общем случае экспертно.

Согласно исследованию [3], многокритериальная МЗЛП (ММЗЛП) (1), (2) равносильна однокритериальной задаче с условиями (1) и максимизацией свертки критериев J (р) = р 1 J 1 + p2 J 2, где p е M = {(р 1 ; р2) е E 21 р i 0 ( i = 1, 2); р 1 + р2 = 1} - вектор параметров, E 2 - двумерное евклидово пространство. Поскольку р2 = 1 - р 1 , то, обозначая р = р 1 , перейдем от ММЗЛП (1), (2) к эквивалентной ей однокритериальной задаче (1) при следующем условии:

J (р) = p J 1 + (1 - р) J 2 ^ max (р е (0;1)). (3)

Для задачи (1), (3) имеет место лемма [2], аналогичная лемме из публикации [1].

Лемма. Для оптимальных значений переменных U n + k ( t ) и X k ( t )( k = 1,..., n ; t = 1,..., T - 1) модели (1), (3) справедливо равенство

U n + k ( t ) = 5 k X k ( t )( k = 1,..., n ; t = 1,..., T - 1). (4)

Здесь и далее * обозначены оптимальные значения переменных и критериев. Условие (4) означает следующее: в оптимуме выручка от реализации в модели (1), (3) (а значит, и в задаче (1), (2)) равна стоимости произведенной по каждому виду продукции.

Рассмотрим задачу двухкритериальной оценки ИП, описанного моделью (1), (2), когда продажа ОПФ предприятием не предполагается, т. е. 5 = 0. Для применения z -преобразования к анализу модели (1), (3), доопределим управления, Uk ( t )( k = 2 n + 1,2 n + 2; t = 1,..., T - 1) до вектора постоянной размерности n + 2 (поскольку переменные Un + k ( t ) ( k = 1,..., n ; t = 1,..., T - 1) можно исключить согласно (4)), полагая отсутствующие переменные равными нулю, т. е. дополняя указанную задачу условиями U 2 n +1 ( t ) < 0 , U 2 n +1 ( t ) >  0 ; U 2 n +2 ( t ) <  0 , u 2 n + 2 ( t ) 0 ( t = 1,..., T - 1). Учитывая (4) и начальные условия, запишем ее ограничения и критерии единообразно, перейдя к МЗЛП:

xk ( t + 1) = x k ( t ) + u k ( t )

(k = 1,..., n; t = 0,..., T-1), xn+1 (t + D = -S xk (t) / Tk + xn+1 (t) + k=1

+ £X( t ) ( t = 0,..., T - 1)

k =1

x n +2 ( t + 1) = S Y k x k ( t ) - 9 x n +1 ( t ) + x n +2 ( t ) - k =1

- S U ( t ) + u 2 n +1 ( t ) + u 2 n +2 ( t ) ( t = 0 T - 1) ;

X k (0) = 0 ( k = 1,..., n + 2);               (5)

X n +2 ( t + 1) 0 ( t = 0,..., T - 1),

S о k X k ( t ) - a 2 X n +1 ( t ) > 0 ( t = 0,..., T - 1), k =1

U 2 n +1 ( t ) 1 0 , U 2 n +2 ( t ) K a ( t = 0);

U 2 n+1( t) < 0, U 2 n+2 ( t) < 0( t = 1,..., T-1), uk(t) > 0 (k = 1,..., n; 2n +1, 2n + 2; t = 0,..., T -1);

J(ц) = ц J1 + (1 - ц) J2 > max (ц e (0; 1)), где

_    T -1

J1 =S t=0

n

S Y k x k ( t ) - 9 x n +1 ( t ) - U 2 n +1 ( t ) - U 2 n +2 ( t )

_ k =1

(1 + r ) t

_    T -1

J 2 = S t=0

n

S to k x k ( t ) + 9 x n +1( t )

. k =1

Y k = аз / Tk + Y 5 k ;

( 1 + r ) t

to k = - а з / T k + p5 k ; о k = (1 - P)5 k - 1/ T k ( k = 1^..., n ). По построению для целевых критериев J (ц) и J (ц) соответственно задач (1), (3) и (5) имеет место неравенство: J (ц) J (ц). В частности,

*

J *(ц) J (ц)(ц e (0;1)). (6)

Полагая z = 1 + r > 1, перейдем в задаче (5) к пределу при T ^да . В силу предпосылки (3), T ^ +да fi Tk ^ +да , откуда следует, что y k ^ Y^ k , to k ^ p6 k , о k ^ (1 - в)3 k ( k = 1,..., n ). Применяя к последней МЗЛП z -преобразование и учитывая формулу Z ( x ( t + 1)) = z [ X ( z ) - x (0) ] , получим двухпараметрическую (по параметрам ц и z ) статическую задачу линейного программирования (ЗЛП):

zXk ( z ) = X k ( z ) + Uk ( z )( k = 1,..., n ),

n zXn+1( z) = Xn+1( z) + S Uk (z), k=1

n zXn+2( z) = yS 5 kXk (z) - 9Xn+1( z) + k=1

n

+ X n +2 ( z ) - S U k ( z ) + U 2 n +1 ( z ) + U 2 n +2 ( z ) , k =1

zX n +2 ( z ) 0,

(1 - e) S 5 k X k ( z ) - a 2 X n +1 ( z ) 0,         (7)

k =1

U 2 n +1 ( z ) 1 0 , U 2 n +2 ( z ) K 0 ;

Uk ( z ) 0 ( k = 1,..., n ;2 n + 1,2 n + 2);

J (ц, z ) = ц J 1 ( z ) + (1 - ц) J 2 ( z ) ^ max (ц e (0;1); z 1),

n где J 1(z) = yS 5kXk (z) - 9Xn+1 (z) - U2n+1 (z) - U2n+2 (z) , k=1

n

J 2( z) = pS 5 kXk (z) + 9Xn+1( z), k=1

Uj(z) = j^Uj(t)z-t (j = 1,..., n;2n +1,2n + 2)             и t = 0 да

Xk ( z ) = S x k ( t ) z - t ( k = 1,..., n + 2) - z -изображения уп- t =0

равляющих_и фазовых переменных задачи (5), J (ц, z ) = T im J (ц), причем при T ^ +да по построению J * (ц) J * (ц, z ) (ц e (0; 1)), откуда в силу выражения (6) получим

*

J (ц) J (ц, z )(ц e (0;1); z 1).          (8)

Выражая из уравнений в ЗЛП (7) переменные Xk ( z ) ( k = 1,..., n + 2), имеем

X k ( z ) = Uk^z ) ( k = 1,..., n ), z - 1

n

S U k ( z )

Xn+1( z) =            - z -1

X n +2 ( z ) =

n

S (y5 k - 9 - z + 1) U k ( z ) + ( z - 1)( U 2 n +1 ( z ) + U 2 n +2 ( z ))

_ k =1 ________________________________________________________________ ( z - 1)2                                .

Подставляя последние формулы в оставшиеся выражения указанной задачи и учитывая, что z = 1 + r >  1, запишем параметрическую статическую задачу, совпадающую с приведенной в работе [2] и эквивалентную модели ZB 1:

n

S (9 + r - y5 k )Uk ( z ) - r(U 2 n +1 ( z ) + U 2 n +2 ( z )) 0, k =1

n

S (9 - y5 k )Uk ( z ) 0, k =1

U 2 n +1 ( z ) 1 0 , U 2 n +2 ( z ) K 0 ,

Uk ( z ) 0( k = 1,..., n ;2 n + 1,2 n + 2),         (9)

J ( ц, z ) = 1 S Ф kUk ( z ) - Ц [ U 2 n +1( z ) + r k

+ U 2 n +2( z )] ^ тах( ц e (0; 1 ))

где ф k = [1 - 2ц]9 + [цY + (1 - ц)р]5 k ( k = 1,..., n ; ц e (0; 1)).

Заметим, что в указанной работе ММЗЛП (1), (2) получена формально из более общей модели A , отличающейся от нее дополнительным ограничением

Un+k (t) < qk (t +1) (k = 1,..., n; t = T 2,._ T-1), (10) при условиях qk(t+1) >+ /

( k = 1,..., n ; t = T 2,._ T - 1);                (11)

T 1 = T 2 = 1.

Здесь qk ( t + 1) ( k = 1,..., n ; t = T2, ..., T - 1) - прогнозный спрос на продукцию k -го типа в момент t + 1, T 1 и T 2 - моменты окончания инвестирования и начала производства. При этом, согласно работе [2], в модели ZA , полученной после применения z -оператора к задаче А , появляются агрегированные ограничения

U n + к ( z ) < Q k ( z ), Un + к ( z ) <a U 1 z 1 ( к = 1,..., n ), (12) z - 1

def ”'

где Qk ( z ) = ^ qk ( t + 1) z t ( к = 1,..., n ). Первое из них в t = T1

силу условия (11) примет вид Un + к ( z ) < +да ( к = 1,..., n ), и исключается как избыточное, а второе, являющееся единственным ограничением на переменные Un + к ( z ) ( к = 1,..., n ; z 1) сверху, в оптимуме выполняется как равенство, поскольку целевая функция в ZA максимизируется, а ее коэффициенты при указанных переменных положительны:

С к ( z ) = "^ 'z ' ( к = 1,..., n ). z -1

Указанное выше совпадение ЗЛП (9) с задачей из монографии [2] подтверждает правильность выводов, полученных в отмеченной работе формальным многократным предельным переходом согласно 1-го из условий (11). При этом рассмотренный выше подход является более строгим (поскольку не использует указанного предельного перехода, а значит, не требует обоснования его правомерности). Отметим, что, если не учитывать «динамические» соотношения (4) (достаточно рутинное доказательство которых приведено в том же источнике), то получим ЗЛП, в которой, в отличие от задачи ZA , отсутствует 1-е из неравенств (12). Учитывая «статические» равенства (13), вновь получим ЗЛП (9). Авторская методика из указанной монографии [2] может рассматриваться в качестве правдоподобных и простых наводящих рассуждений (требующих последующего обоснования), позволяющих трактовать задачи (1), (2) и (9) как частные предельные случаи при условиях (11) моделей A и ZA соответственно. Решение z -задачи (9), найденное там же, определяется следующей теоремой.

Теорема 1. Если найдется такой номер к 0 е {1,..., n }, что б к = max б к , причем — б. 6 + r , то в ЗЛП (9)

0 к=1,..., n                        Y 0 Y существует нетривиальное решение

U k о ( z ) = r ( I о + K о )/(6 + r - Y6 к о );

U n + к о ( z ) = б к о ( I о + K о )/(6 + r - y8 к о ),

U k ( z ) = о; U n + к ( z ) = о( к = 1,..., n ; к # к о );

U 2 n + 1 ( z ) = I о , U 2 n + 1 ( z ) = K о , которому в пространстве критериев соответствует единств енная ненулевая Парето-точка с координатами =*      ( 2уб - 26 - r )          =•      ( 6+рб^ )

J ( z ) = I --------l [ Io + K ]; J 2 ( z ) = I---- k-2— l [ Io + K ]

IV/                              со о 2                                Н~о о

I 6 + r - y6l l                      I 6 + r - y6l l

V             kо J                         V            k о J для всех значений ц, задаваемых соотношениями "б к о > (26 + r )/(2y), ц е (о;1)

б^ < (26 + r) / (2y), ц е (о; (6 + рб^ )/[36 + r + (р - 2y)6ц ]). о                                                                        о

При этом оптимальное значение свертки J (ц, z )

в

указанной

ЗЛП

равно

=*

J ( ц, z ) =

U21 , (z) = Uk(z); U t (z) = Un+к(z) (к = 1,..., n); 1 k I1

V1(z) = U2n+1(z); Vг(z) = U2n+1(z)(14)

и переходя формально к пределу при n ^ да , переформулируем указанную теорему в следующем виде.

Теорема 2. Если б, = max бк таково, что к =1,...

6 . 6 + r

- < б1 <----- ,                   (15)

γγ то в ЗЛП (9) даже при неограниченном количестве видов продукции существует нетривиальное решение

U ( z ) = r ( I о + K о )/(6 + r - y6 1 );

Ц ; ( z ) = б 1 ( I о + K о )/(6 + r - y6 1 );

U k ( z ) = о( к = 3,...);

V( z) = Iо, Vi( z) = Kо, которому в пространстве критериев соответствует единственная ненулевая Парето-точка с координатами

_ *

J 1 ( z ) =

2y6 1 - 26 - r )

—----J— l[ I о + K о] ;

6 + r - y61 J

_ *

J 2 ( z ) =

6+рб 1

6 + r - y61

для всех значений ц , задаваемых соотношениями "8 1 > (26 + r )/(2у), ц е (о;1)

8 1 < (26 + r ) / (2у), ц е (о; (6 + рб 1 )/[36 + r + (р - 2у)б 1 ]).

=*

При этом предел оптимального значения свертки J (ц, z ) в         указанной         задаче         равен

= *.    . (6+рб, + [<2у-р)б, -36 -r]ц 1, lim J (ц, z) = I —t-1---1---------H l[Iо + Kо]. n^”           V         6 + r - уб1          J

Справедливы следующие теоремы, доказательство которых также приведено в монографии [2].

Теорема 3. В ЗЛП вида f (x) = (c, x) ^ max;

n gj(x) = Zajixi + aj < о(j =1’-,k)       (17)

i =1

решение существует тогда и только тогда, когда множество X , задаваемое ее ограничениями, не пусто и функция f ( x ) ограничена сверху на этом множестве.

Теорема 4. Оптимальное значение свертки J * (ц) в проекте, описываемом моделью A , есть неубывающая функция от параметров T , n , T 1 , Y, б; б к , q k ( t ) ( к е {1,..., n }; t е { T 2 + 1,..., T }); I о , K о и невозрастающая от T 2 ,6 и r (T , n , T 1 , T 2 е {1, 2,...}) при неизменных значениях остальных параметров и ц е [ о;1 ] .

Отметим, что доказательство теоремы 3, предложенное в работе [2], не использует конечность числа переменных и ограничений задачи (17), задающих множество X . Поэтому данная теорема выполняется и в том случае, когда число переменных и/или неравенств бесконечно, т. е. полагаем в частности, что в указанной задаче n ^ да и/или к ^ да . _, _*

Поскольку J (ц, z ) = J (ц, z ) (ц е (о;1); z 1), то при n ^ да из теоремы 2 и неравенства (8) получим условие, аналогичное следствию 3.5.1 из работы [2] для конечного n :

lim lim J * (ц) <

T ^m n

6 + р5 , + [(2y - р)5 .

6 + r - y^x

36 r I [ 1 0 + K 0 ],(18)

из которого следует, что при условиях (15) и неограниченных n , T целевая функция J (ц) (ц е (0; 1)) МЗЛП (1), (3) ограничена сверху. Поскольку управление (полученное перенумерацией исходных пер еме нных указанной задачи, со-

Ji | A , Ji ( z )| ZA – значения i -го критерия в Парето-точке задач A и ZA соответственно, из которого в силу (16) следует, что даже при условиях теоремы 5 (следствия) координаты Парето-точек критериального пространства задачи A ограничены сверху:

J ,4 <

о тве тствующей заменам (14)) u

2 к - 1

u M ,( t ) = 0( к = 1,...; t = 1,...);

2 кI v2 n+2( t) = 0( t = 0)     является

( t ) = 0( к = 1,...; t = 0,...);

v 2 n +i ( t ) = 0( t = 0), допустимым   при

J * L <

2y61-26 - r 1

6 + r - y6 , J1 1 0 + K 0] ;

.( p I [ + ]

( 6 + r - y61 )

T ^ m; n ^ m , а задачи (1), (2) и (1), (3) эквивалентны, то по теореме 3 доказали следующую теорему.

Теорема 5. Если 61 = max 6 к удовлетворяет условиям к = 1, ...

(15), то в ММЗЛП (1), (2) даже при неограниченном количестве видов продукции на бесконечном временном интервале существует нетривиальное решение, причем имеет место неравенство (18).

Для содержательной корректности предельного перехода при n ^ м в теореме 2 номер к 0 е {1,..., n } должен перенумеровываться так, чтобы он не зависел от n .

При условиях (11) из теоремы 4 и теоремы Вейершт-расса о существовании предела ограниченной монотонной последовательности следует справедливость теоремы 5 и для модели A с конечным спросом и на конечном интервале времени. Более того, принимая во внимание, что предложенный в этой работе подход строго обосновывает правомерность предельного перехода согласно 1-го из условий (11) и получаемых с его помощью результатов, а значит, позволяет с формальной и экономической точки зрения интрепретировать модели B 1 и ZB 1 как частные асимптотические варианты моделей A и ZA соответственно при указанных условиях, из теоремы 5 получим следствие, являющееся по сути ее переформулировкой.

Следствие. Если 61 = max 6 к удовлетворяет условиям (15), то в задаче А даже пр инеограниченном количестве видов продукции и неограниченном спросе в каждый момент времени в период производства и по каждому виду на бесконечном временном интервале существует нетривиальное решение, причем имеет место неравенство (18).

Заметим, что теорема 5 (или следствие) носит содержательно парадоксальный характер: хотя основные экономические характеристики ИП, описываемого моделью B 1 (или A ) – число различных видов ОПФ, горизонт планирования (и спрос на продукцию) бесконечны, значение свертки ее критериев, соответствующее Парето-точке критериального пространства, является конечным при условиях (15) и любых значениях остальных параметров проекта. Согласно численных экспериментов, представленных в главе 2 работы [2], справедливо эмпирическое соотношение J * | . J * (z ) | ZA ( i = 1,2; ц е (0;1)), где

Этот же факт можно обосновать строго от противного, учитывая, что ц е (0; 1), вид свертки в (3) и условие (18). Содержательно указанный результат объясняется тем, что если максимальная фондоотдача ОПФ каждого вида (характеризующая их экономическую эффективность) заключена в диапазоне (15), то при условиях теоремы 5 (следствия) в проекте оптимизации реальных инвестиций предприятия, описываемом моделью B 1 ( A ), даже в оптимуме целевые критерии ограничены сверху.

Отметим, что если не делать замен (14), то теорему 5 сложно интерпретировать полностью содержательно. Аналогично можно показать, что если в модели A отсутствует ограничение (1 0) для некоторого

к е {1,..., n }; t е { T 2 ,

..

., T - 1}, то в соответствующей z -за-

даче отсутствует 1-е из ограничений (12) для указанного номера k , что формально можно получить из упомянутого ограничения при q k ( t ) ^ +м , т. е. при Q k ( z ) ^ +м .

На основе z-преобразования можно дать более простой метод доказательства леммы. Допустим, что существуют к е {1,..., n }; t е {1,..., T - 1}, при которых равенство (4) в оптимуме не выполняется, т е. u n+ к ( t ) 6 kx*k ( t ). Поскольку в ЗЛП оптимум достигается на границе допустимого множества, то последнее неравенство можно исключить, поскольку оно является строгим. Тогда, аналогично предыдущему замечанию, 2-е из условий (12) отсутствует, что противоречит (13), а значит, указанное строгое неравенство невозможно, откуда следует (4).

Статья научная