О применении Z-преобразования к анализу многокритериальных линейных моделей экономической динамики
Автор: Медведев Алексей Викторович
Журнал: Сибирский аэрокосмический журнал @vestnik-sibsau
Рубрика: Математика, механика, информатика
Статья в выпуске: 3 (16), 2007 года.
Бесплатный доступ
Предлагается обобщение подхода к анализу линейных многошаговых задач оптимального управления, связанного с использованием z-преобразования, для доказательства существования их решения и получения достаточных условий неэффективности инвестиционных проектов на примере двухкритериалъной модели развития региональной экономики.
Короткий адрес: https://sciup.org/148175554
IDR: 148175554
Текст научной статьи О применении Z-преобразования к анализу многокритериальных линейных моделей экономической динамики
Процесс оптимизации реальных инвестиций производственных отраслей в регионе можно представить в виде следующей многокритериальной многошаговой задачи линейного программирования (ММЗЛП) [1]:
X k ( t + 1) = X k ( t ) + U k ( t )
( k = 1,..., n ; t = 0,..., T - 1 ) , x n + 1( t + 1) = x n + 1( t ) +
+ L Uk(t) (t = 0,..., T2 -1), k=1
X n + 1( t + 1) =
= -:Lxk(t)/Tk + Xn+1(t) + Luk(t) x k=1 k=1
x(t = T2,..., T-1), xn+2 (t+ 1) = -a2 xn+1(t) + xn+2 (t) -
Un + k ( t ) < ^ kXk ( t ) x x( k = 1,..., n ; t = T 2,..., T - 1 ) ; Xn + 3 ( T 1 ) < I 0 , U 2 n + 2(0) < K 0 ,
U k ( t ) > 0 ( k = 1,..., n ; t = 0,..., T - 1 ) , U 2 n + 1( t ) > 0 ( t = 0,..., T 1 - 1 ) , u 2 n + 2 (0) > 0;
J = { J 1 , J 2 } ^ max,
r T 1 - 1 U 2 n + 1( t )
гд е J 1 = - L----- r t = 0 (1 + r ) t
T - 1 + l l t = T 2
V Xk(t) “3 L Л k=1 Tk
-
U 2 n + 2(0) +
n
- 6 Xn + 1 ( t ) + г L Un + k ( t )
-------------k = 1-------J +
(1 + r ) t
Д x n + 1 ( T ) о T T - 1 ,
n
- L u k ( t ) + U 2 n + 1( t ) + U 2 n + 2( t ) ( t = 0 ) , k = 1
T - 1
J2 = L2
t = T 2
- “ 3 L n X k^ t ) + 6 X n + 1 ( t ) + c L U n + k ( t ) k = 1 T k k = 1
( 1 + r ) t
соот-
x n + 2 ( t + 1) = -a 2 x n + 1( t ) + x n + 2 ( t ) - n
- L u k ( t ) + U 2 n + 1( t ) ( t = 1,..., T 2 - 1 )
k = 1
x n + 2 ( t + 1) = “ 3 L Xk ^) - 6 x n + 1 ( t ) + x n + 2 ( t ) k = 1 T k
ветственно дисконтированная сумма собственных средств
^^^^^B
-^L u k ( t ) + k = 1
+ г L U n + k ( t ) + u 2 n + 1 ( t ) ( t = T 2,..., T 1 - 1 )
k = 1
xn+2 (t+ 1) = “3 L ^k^ - 6xn+1(t) + xn+2 (t) - k=1 Tk
- L U k ( t ) + г L U n + k ( t ) ( t = T 1 ,..., T - 1 ) ; k = 1 k = 1
xn+3(t+ 1) = xn+3(t) + U2n+1(t) (t = 0,..., T1 -1) , xn+3(t+1) = xn+3(t) (t = t 1,...,t-1);
X k (0) = 0 ( k = 1,..., n + 3);
X n + 2 ( t ) > 0 ( t = 1,..., T ),
X k ( t ) ..
-L^--“2 Xn+1( t) + k=1 Tk
n
+ ( i -B ) L Un + k ( t ) ^ 0 ( t = t 2,..., t - 1 ) ;
k = 1
0 < Un + k ( t ) ^ q k ( t + 1),
производственного сектора и регионального центра (налоговых поступлений в регион); U k ( t ) ( t = 0,..., T - 1 ) , Un + k ( t ) ( k = 1,..., n ; t = T 2,..., T - 1 ) , U 2 n + 1( t ) ( t = 0,..., T ^ -1 ) и u 2 n + 2 (0) - стоимость приобретаемых основных про
изводственных фондов (ОПФ), выручка от реализации продукции k-го типа, внешние и внутренние инвестиции соот
ветственно; X k ( t ) , Xn + 1( t ) , Xn + 2 ( t ) , Xn + 3( t )( k = 1,..., n ; t = 0,..., T ) - соответственно накопленная стоимость всех ОПФ к - го
типа, остаточная стоимость всех ОПФ, текущие денежные средства предприятия и накопленные суммы внешних инвестиций в момент t; q k ( t + 1) ( t = T 2,..., T - 1 ) ,, V k , T k , c k и P k - соответственно прогнозный спрос в стоимостном выражении для момента t + 1, производительность, срок службы, стоимость единицы ОПФ и стоимость единицы продукции k-го типа; 1 0 , K 0 - суммы внешних и внутренних инвестиций, выделяемых на весь срок действия инвестиционных проектов (ИП); “ 1 , “ 2 , “ 3 , “ 4 - ставки налогов на добавленную стоимость (НДС), на имущество (НИ), на прибыль (НП) и единый социальный налог (ЕСН) соответственно (НДС включается в цену продукции, поэтому можно считать, что “ 1 = 0 ); в - доля выручки от реализации, выделяемая на фонд оплаты труда (ФОТ); T 1 , T 2, T (1 < T 2 < T 1 < T )-соответ-
ственно моменты завершения внешнего инвестирования, начала производства и срок действия ИП; 6 = (1 - “ 3 )“ 2 , д k = PkVk / ck ( k = 1,-, n ), г = (1 - “3 )(1 - b), c = (1 - b)“ 3 + “ 4 B, r - ставка доходности ИП; д (0 < 8 < 1) - доля остаточной стоимости всех ОПФ на момент t = Т от ее балансовой стоимости, определяемая в общем случае экспертно.
В работе [2] рассмотрен пример применения z-преобразования для анализа однокритериальной МЗЛП. Обобщим указанный прием на случай ММЗЛП. Согласно источникам [3; 4], ММЗЛП (1), (2) равносильна однокритериальной МЗЛП, в которой векторный целевой критерий (2) заменен на скалярный: J(м) = м^J। + м2J2 , где me M = {(м1;м2)е E2 | мi > 0;м1 + М2 = 1} -вектор параметров, E2 - двумерное евклидово пространство. Таким об разом, учитывая, что м2 = 1 - М1 и полагая def м = М1,
мож
но свести исследование ММЗЛП (1), (2) к исследованию однопараметрической МЗЛП (1) при условии
J (м) = m J 1 + (1 - м) J 2 ^ max (3) где 0<ц<1.
Рассмотрим задачу оценки стоимости ИП, описанного моделью (1), (2), когда 8=0, т. е. продажа ОПФ регионом не предполагается. Для упрощения дальнейшего анализа доопределим управления u n + k ( t ) ( k = 1,..., n ; t = 0,..., T 2 - 1 ) , u 2 n + 1( t ) ( t = T 1,..., T - 1 ) , u 2 n + 2 ( t ) ( t = 1,-, T — 1 ) , полагаяих равными нулю и сопоставив им в исходных данных нулевые элементы. Дополнив задачи (1), (3) неравенствами
X n + 3 ( t ) ^ 1 0 ( t = 1,..., T );
un+k (t) - 0 , un+k(t) > 0 (k = 1, ..., n; t = 0, ..., T2 - 1) , u 2 n+1( t) - 0 , u 2 n+1( t) > 0( t = T V", T - 1) ;
u 2 n + 2 ( t ) - 0 , u 2 n + 2 ( t ) > 0( t = 1,..., T - 1) , (4)
перейдем к однокритериальной задаче с условием
-(9 + z -1) ^ELUk (z) + k=1
+r(z -1) EE Un+k (z) + k=1
+ ( z - 1)( U 2 n + 1( z ) + U 2 n + 2 ( z )) > 0 ,
« 2 ^Lut ( z ) ,
--k =L---+ (1 - в) E U n + k ( z ) > 0 , (6) z - 1 k = 1
U n + k ( z ) ^ Q k ( z ) ,
U n + k ( z ) ^ ^ k Uk^T ( k = 1,..., n ), z - 1
U 2 n + 1( z ) - I 0 , U 2 n + 2 ( z ) - K 0 ,
Uk ( z ) > 0( k = 1,..., 2 n + 2),
[1 - 2 m]9 E^ U k ( z )
j (m, z ) =--------k^1------+ z - 1
+[мг + (1 - м)с] EE Un+k(z) - k=1
- H U 2 n + 1 ( z ) - U 2 n + 2 ( z )] ^ max, (7)
^
где U j ( z ) = E u j ( t ) z - t ; (j = 1,..., 2n + 2) - z-изображения t = 0
управляющих переменных задач (1), (4), (5), J ( m, z ) = lim J (m) ,
T -^^
где j 1 =
T - 1 Е t = 0
T - 1
j 2 = E t = 0
9,
( t ) = •
J ( м) = м J 1 + (1 - м) J 2 ^ max , (5)
nn
- « 3 Е T k ( t ) xk ( t ) + 0 0 ( t ) xn + 1 ( t ) + г Е un + k ( t ) - u 2 n + 1 ( t ) - u 2 n + 2 ( t )
. k = 1 k = 1 :
(1 + r ) t
nn
« 3 Е T k ( t ) x k ( t ) - 6 0 ( t ) xn + 1 ( t ) + с Е un + k ( t )
. k = 1 k = 1:
(1 + r ) t
def ^
Q k ( z ) = E q k ( t + 1) z - t ( k = 1,..., n ). (8)
t = T 2
Докажем, что в задачах (1), (2) существует решение, используя для этого, в отличие от работы [1], z-преобразование. Покажем сначала, что в ЗЛП (6), (7) допустимое множество непусто. Так как функциональный ряд (8) при постоянных и одинаковых q k ( t + 1) есть геометрическая прогрессия, то, очевидно, справедливы неравенства вида
0, t = 0,..., T 2 - 1
-9, t = T 2..... T -1’
T k ( t ) = {
0, t = 0,..., T 2 - 1
- 1/ T k , t = T 2,...., T - 1,
а управляющий вектор u(t) имеет постоянную на каждом шаге размерность 2и + 2. Здесь и далее значком (*) обозначены оптимальные значения переменных и целевых функций, чертой сверху - соответствующая целевая функция и параметры задачи с управлениями постоянной размерности.
Применим для исследования задач (1), (4), (5) z-преобразование. С этой целью устремим Т к бесконечности и положим z = 1 + г > 1, T 2 = 1 . Учитывая соотношение Z ( x ( t + 1)) = z [ X ( z ) - x (0) ] и импликацию T ^+тс^ T k ^+~ ( k = 1,..., n ), а также исключая
^
xk ( z ) = Е X k ( t ) z - t ( k = 1,..., n + 3) - z-изображения фазо- t = 0
вых переменных указанной задачи, получим двухпараметрическую (по параметрам м и z) статическую задачу линейного программирования (ЗЛП), аналогичную приведенной в [2]:
Qk ( z ) ----- kprx ( k = 1,..., n ), r>0, откуда, в частности, r (1 + r ) T - 1
при T 2 = 1 имеем
Q k ( z ) - q k ( k = 1,..., n ), (9)
r
где q k =
max q k ( t + 1) ( k = 1, t = T 2,...
n ) - наибольший прогноз
ный спрос на продукцию k-го типа за весь период произ
водства.
Будем предполагать, что прогнозный спрос по всем видам продукции конечен в течение всего периода производства, описываемого моделью (1), (2), т. е. справедливы условия q k < +~ ( k = 1,..., n ), из которых, с учетом условия (9), следует, что
Q k ( z ) <+“ ( k = 1,..., n ; z > 1). (10)
Найдем условия, при которых существует решение задач (1), (2) на бесконечном временном интервале. Ограниченность переменных U 2 n + 1 ( z ) и U 2 n + 2 ( z ) следует непосредственно из пятого, шестого и седьмого условий в выражении (6). Из условия (10), третьего и седьмого ограничений в выражении (6) следует ограниченность
переменных U n + k ( z ) ( k = 1,..., n ) приz > 1. Из второго неравенства в (6) следует, что n , х ( z - 1)(1 - в) n , х
У Uk(z) <---------У Un+k(z). Учитывая седьмое ус- k=1 “2 k=1
ловие в выражении (6), дляк = 1, ..., и получимследую-
. . (z -1)(1 - в) П _ щую систему оценок: 0 < Uk (z) < --------- У Qk (z), из “2 k=1
которой следует ограниченность переменных U k ( z ) при “ 2 > 0 . Заметим, что если “ 2 = 0 , то 6 = 0, и первое неравенство в (6) преобразуется к виду y^ Uk ( z ) <-г У U n + k ( z ) + U 2 n + 1 ( z ) + U 2 n + 2 ( z ) . Следова- k = 1 k = 1
тельно, и в этом случае переменные U k ( z )( k = 1,..., n ) являются ограниченными. Поскольку изображения вида U k ( z ) = 0 ( k = 1,..., 2 n + 2; z > 1) удовлетворяют ограничениям ЗЛП (6), (7), то допустимое множество D ( z ) ( z > 1) указанной задачи непусто. В силу нестрогости ограничений (6) указанное множество замкнуто, а значит, компактно. Легко проверить, что управление вида u k ( t ) = 0 ( k = 1,..., n ; t = 0,..., T - 1) ; U n + k ( t ) = 0( k = 1,..., n ; t = T 2,..., T - 1) ; u 2 n + 1 ( t ) = 0( t = 0,..., T 1 - 1), u 2 n + 2 ( t ) = 0( t = 0) является допустимым в ММЗЛП (1), (2), поэтому допустимое множество D(z) этой задачи при T ^ +~ является непустым. Кроме того, очевидно, что при переходе в задаче (1) к пределу при T ^ +~ имеет место соотношение
D ( z ) с Пр D ( z )( z > 1) ,
где Пр D ( z ) - проекция множества d ( z ) на оси, соответствующие определенным в модели (1), (2) управлениям u j ( t ) ( j = 1 —, 2 n + 2) при T ^ +~ (4). Непустота множества d ( z ) следует из (11) и непустотыD(z). Множество D ( z ) ограничено, поэтому и D(z) также является ограниченным в силу формулы (11). По теореме Вейерштрасса, в силу непрерывности J ( м, z ) и J ( м) по своим аргументам (при фиксированных параметрах z > 1 и 0< ц <1) в задачах (1), (2) и (6), (7) существует решение, т. е. теорема доказана.
Теорема 1. Если выполняются условия qk < +“ (k = 1,..., n);T ^ +~;r > 0;T2 = 1,(12)
то задачи (1), (2) и (6), (7) имеют решение, причем справедливы следующие оценки:
0 < Uk (z) <(z -1)(1 - в) У Qk (z), “2
0 < U n + k ( z ) < Q k ( z ) ( k = 1 .... n ) ,
0 < U2n+1(z) < 10 , 0 < U2n+2(z) < K0 .(13)
Заметим, что, зная диапазоны выражения (13), изменения z-изображений Uj ( z ) ( j = 1, -, 2 n + 2) , можно получить интервал изменения соответствующих управлений задачи (1), (2). Покажем это, например, для
U k ( t )( k = 1,..., n ; t = 0,..., T - 1). Пусть П(г) - множество изображений, удовлетворяющих ограничениям (13). Тогда по теореме 1 D ( z ) с П( z ) ( z > 1) , откуда, в силу (11), имеем
D ( z ) с П( z )( z > 1) , (14)
где ПрП(г) - соответствующая проекция множества П(/).
Пусть uk(t)(k = 1,...,n; t = 0,...,T-1) - допустимое управление в задаче (1), (2). Тогда с учетом формулы (14) и того, что все указанные выше управляющие перемен ные определены и в задаче (6), (7), для uk (t) выполняется первое условие в (13):
^
0 <У ukz-t < k=0
< ( z - 1 )( 1 - в) j^ Q k ( z )( k = 1,..., n ). (15)
“ 2 k = 1
С другой стороны, def
D t с lim Пр D t 0 = Пр D ( z ) ( z > 1) ,(16)
T ° ^~ u ( t ) = 0( t = T ,...) u ( t ) = 0( t = T ,...)
где ПрDT0 - проекция множества управлений DT", u (t) = 0( t = T,...)
определяемого ограничениями выражение (1) для t 0
шагов, на оси uk ( t )( k = 1,..., n ; t = 0,..., T - 1), u 2 n + 1( t ) = 0( t = 0,..., T 1 - 1), u 2 n + 2( t ) = 0( t = 0) при
T0 ^ +^ , Dt - множество управлений, заданное условиями (1) для Тшагов. Следовательно, в силу (16), дополнительно к ограничениям (1) можем полагать, что u0(t) = 0(k = 1,..., n; t = T,...). (17)
k
Тогда из выражения (15) получим v 0 —t / (z -1)(1 - в) n
У ukz = У ukz < -------- У Qk(z) , откуда при k=0 k=0 “2
0 -1 ^( z -1)(1 - в), k^ 1, „., и имеем ukz <----------У Qk(z) и, следова-
“2
0 (z -1)(1 - в)zt v тельно, uk <----------У Qk(z). Поскольку
“2
u 0 ( t ) ( k = 1,..., n ; t = 0,..., T - 1) произвольно, то последнее ограничение имеет место для любого допустимого в задачах (1), (2) управления u k ( t ) ( t = 0,..., T - 1) .Для остальных управлений оценки сверху получаются аналогично. Таким образом, справедливо следствие теоремы 1.
Следствие 1. Для управляющих переменных задачи (1), (2) имеют место следующие неравенства:
0 < u k ( t ) < ( z 1)(1 в) z t У Q k ( z ) “ 2 k = 1
( k = 1,..., n ; t = 0,..., T - 1),
0 < u 2 n + 1( t ) < 1 0 z t ( t = 0,..., T - 1),
0 < u 2 n + 2(0) < K 0 . (18)
Несмотря на то, что некоторые оценки на управления в МЗЛП можно получить и из уравнений движения и ограничений задачи, следует отметить, что предлагаемый подход к определению диапазонов управляющих переменных, основанный на идеях операционного исчисления, является универсальным и может применяться при увеличении количества переменных и неравенств модели и использоваться для быстрой оценки поискового пространства переменных многокритериальных МЗЛП и численной реализации решения указанного класса задач на основе соответствующих алгоритмов. Рассмотрим важные следствия и приложения используемого подхода. Для перенесения получаемых с помощью z-преобразования результатов на конечный временной интервал сформулируем теорему 2 и следствие из нее, справедливое приц^О, ц= 1.
Теорема 2. Функция j T ( м) является неубывающей функцией переменной Т при фиксированных m e [ 0;1 ] .
Следствие. J * (i = 1, 2) есть неубывающие функции от параметра Т.
По построению для целевых критериев J ( м) и j ( м) соответственно задач (1), (3) и (1), (4), (5) имеет место неравенство: J ( m) < J (м) .Вчастности,
J ( m) < J (м), 0 < м < 1 . (19)
При T ^ +~ из теорем 1,2 получим: j *(ц) < J *( р , z ). В свою очередь, из (19) и последнего неравенства следует выражение вида
J ( m) < J (м z ) (0 < м < 1) . (20)
Переходя от задачи (6), (7) к эквивалентной МЗЛП (зависящей от двух параметров миг), аналогично тому, как это сделано в источнике [2], найдем ее решение в следующем виде:
Uk ( z ) = rrQk ( z )/9; U * + k ( z ) = Qk ( z ), m e (0;1/2);
Uk ( z ) = r [ P kQk ( z )/Д k + [1 - P k +1Kr Qk ( z ) - ° k № U»+k ( z ) = Qk ;p k + 1 e [0;1],m = 1/2;
Uk i z ) = 0; U* + k ( z ) = 0, 8 k < Щ И 2 /(1 - =),m e (1/2;1); (21)
uk ( z ) = c krQk ( z )/д k ; U*n+k ( z ) = c kQk ( z );c k e [0;1], д k = Ш И 2/(1 - b), Pe (1/2;1);
Uk ( z ) = rQk ( z )/ д k ; иП + k ( z ) = Qk ( z ),д k > Щ « 2 /(1 - =),m e (1/2;1)( k = 1,..., n ),
U * n+1( z) = 0; U 2 n +2( z) = 0, где щ = [2м—1]r- . Подставляя (21), например для мг + (1 - м)с me (0;1/2), в выражение J(м,z), с учетом условия (7), получим
*
J ( м, z ) = ({[1 -2м] r + м}г + n def
+ (1 - m)c) £ Q k ( z ) = r(M, z )(m e (0;1/2>> , (22)
k = 1
и учитывая (3), (20):
J ( m) = m J i +
+ (1 - m) J 2 <Г(м z ) (m e (0;1/2)) . (23)
Анализ двухкритериальной модели (1), (2) показал, что значения критериев J* ( м) ( г = 1,2) , соответствующие оптимальному значению функции j * ( м) в МЗЛП (1),(3) при фиксированном m e (0;1) , неотрицательны:
J* ( м) > 0( г = 1,2;m e (0;1)) , (24)
где J * ( m) ( I = 1,2) -вкладкритерия Ji (м) ( г = 1,2) вопти-мальное значение целевой функции j * ( м) из условия (3) при фиксированном m e (0;1) . Из выражений (23) и (24) следует, что:
J * ( м) <Г * (м, z )( i = 1,2;m e (0;1/2)) , (25)
def def где Г1(р,z) = Г(р,z)/р;Г1(р,z) = F(p,z)/(1 -р)(ре (0;1/2)).
Отметим, что теорема 2 позволяет перенести условие (23), а значит и (25) на конечный интервал времени. Условие (25) означает, что множество Парето в пространстве критериев ММЗЛП (1), (2) мажорируется линией, координаты ( Г 1 (м, z ); Г 2 (м, z )) точек которой зависят от параметров м e (0;1/2) иг> 1. Таким образом, указанная линия может использоваться для оценки этого множества в критериальном пространстве.
Приведем ниже еще одно приложение z-преобразования к оценке эффективности ИП, описываемого моделью (1), (2). Пусть i-е лицо, принимающее решение (ЛПР)
рассчитывает получить в результате реализации некоторого ИПзначение своего целевого критерия не ниже величины J * ( * = 1,2), т. е.
J * ( m) > J * ( г = 1,2) . (26)
С учетом (23) и (26) имеем mJ1 + (1 - m) J2 < mJ1 + (1 - m) J2 < Г(м, z), откуда, учитывая (22) и группируя слагаемые при р, по- дП
J 1 - J 2 + ( [1 - 2 r ]г + c ) ^ Q k ( z ) ^ < ( r r + c ) ^ Q k ( z ) - J 2 . k = 1 J k = 1
Обозначая выражения при р и правую часть последнего неравенства соответственной и В, получим следу ющие оценки на параметр р:
м < B / A , A > 0
м>- B / A , A < 0
. Если,
кроме того,
0 < B / A < 1/2
, то получим сужение ин-
0 <- B / A < 1/2
тервалов изменения параметра р вида
"m e ( 0; B / A ) с ( 0;1/2 ) , B/A > 1/2;
m e ( - B / A ;1/2 ) c ( 0;1/2 ) , - B/A > 0.
Таким образом, доказана теорема.
Теорема 3. Если имеют место условия
B / A < 1/2 9
- B / A > 0 9
A > 0
A < 0
, то в ММЗЛП (1), (3) параметр р изменяется соответственно в диапазонах, определяемых условиями (27).
Теорема 3 позволяет исключить из Парето-множества ММЗЛП (1), (2) те точки, для которых условия (26) не выполняются, и сузить его до множества точек в (1), (2), потенциально приемлемых для ЛПР.
Таким образом, основанный на z-преобразовании подход позволяет обосновывать ограниченность множества переменных, а значит доказывать существование решения в многокритериальных линейных моделях экономической динамики; проводить быструю предварительную оценку поискового пространства переменных многокритериальных МЗЛП для дальнейшей численной реализации решения задач указанного класса на основе соответствующих алгоритмов; строить явно задаваемые линии (поверхности) в критериальном пространстве, ограничивающие Парето-множество исходной ММЗЛП путем перехода к эквивалентной однокритериальной МЗЛП со сверткой критериев; исключать из Парето-множества ММЗЛП те точки, для которых значения целевых критериев ЛПР ниже требуемых им значений, сужая его до множества точек, потенциально приемлемых для ЛПР.
За счет простоты, быстроты и комплексности применения предложенного подхода для получения вышеописанного спектра результатов даже не математиками, он может быть рекомендован инвестиционным аналитикам, управляющим органам предприятий, регионов. при многоплановом экономико-математическом анализе линейных моделей экономической динамики. В настоящее время разработано математическое и программное обеспечение для реализации описанного подхода в виде пакета инвестиционного и финансового анализа [5].