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

Экономические и гуманитарные науки

УДК 517.938

ИССЛЕДОВАНИЕ УПРАВЛЯЕМОЙ МАТЕМАТИЧЕСКОЙ МОДЕЛИ МНОГОСЕКТОРНОЙ ЭКОНОМИКИ С ПОСТОЯННЫМ УРОВНЕМ ПОТРЕБЛЕНИЯ И.С. Потапова

Кафедра математики и методики преподавания математических дисциплин,

ГОУ ВПО «Рязанский государственный университет им. С.А. Есенина»;

Irina00000@yandex.ru

Представлена членом редколлегии профессором В.И. Коноваловым

Ключевые слова и фразы: вектор-функция; внешние воздействия; допустимое управление; кредитные вложения; множество сочетаний; объем основных фондов; система дифференциальных уравнений; фундаментальная матрица.

Аннотация: Построена математическая модель изменения объемов основных производственных фондов. Найдены необходимые и достаточные условия существования управления, при котором модель имеет решение, определяющее развитие отраслей экономики с заранее заданным уровнем потребления.

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

Цель работы - исследовать управляемую математическую модель и найти условия существования управления, при котором модель имеет решение, определяющее процесс развития отраслей экономики с заранее заданным уровнем потребления.

Пусть Х1 (/), Х2 (/), ..., хп (/) - объемы основных производственных фондов (станки, предприятия, здания, денежные средства и др.) отраслей экономики в момент времени /. Будем предполагать, что темп изменения объема фондов Х1 (/) ^й отрасли при / е [0, Т], Т > 0 - некоторое число, пропорционален объему основных производственных фондов i-й отрасли xi (/) и зависит от объемов производственных фондов остальных отраслей. На темп изменения фондов i-й отрасли

могут оказать влияние вложения собственных средств отрасли, кредитные вложения, определяемые вектором и (Г) = (^(Г), U2(t), ..., um (Г)), и внешние воздействия 'цi (Г). Математическая модель изменения объемов фондов многоотраслевой экономики может быть записана в виде системы дифференциальных уравнений

х = А(Г)х + В(Г )и + у(Г), (1)

в которой А(Г) = (Яу (Г))п, а^ (Г) - доля фонда ]-й отрасли, оказывающая влияние

на темп изменения фонда /-й отрасли; В(Г) = (Ь^ (г))”!”, (г) - часть кредитного

вложения и}- в развитие фонда /-й отрасли; у (г) = (^(Г), у2(Г), к, (Г)), (г) -

характеризует внешние воздействия на темп изменения фондов /-й отрасли.

Выпуск продукции z(Г), идущей на потребление в момент времени г е[0, Т], определится равенством z(Г) = С(Г)х(Г), где О(Г) - к х п-матрица, непрерывная на сегменте Г е [0, Т]. Уровень потребления за период времени [0, Т] задается в виде [1]

Т Т

| z(t)dt =| G(t)x(t)dt = б, (2)

о о

в котором б - к-мерный постоянный вектор.

Предположим, что на сегменте [0, Т ] матрицы А(Г), В(Г), 0(Г) определены

и непрерывны, и = {и(Г)} - множество допустимых управлений, и(Г)- т-мерная

вектор-функция, определенная на сегменте [0, Т].

Определение. Система (1) называется управляемой на сегменте [0, Т] во множестве допустимых управлений и, если для любых векторов а е Еп и б е Ек существует управление и(Г) е и, при котором математическая модель (1) имеет решение х(Г), х(0) = а, определенное на сегменте [0, Т], удовлетворяющее равенству (2).

Пусть б - произвольный, но фиксированный к-мерный вектор.

Ставится задача: найти условия управляемости модели (1) на сегменте [0, Т]. Пусть Х(Г) - фундаментальная матрица системы х = А(Г)х, Х(0) = Е, и(Г) е и. Тогда решение х(Г) (х(0) = а) математической модели (1) определится равенством

Г

х(Г) = Х(Г)а + Х(Г) | X-1 ©(В©и© + у©) ^.

0

Следовательно, для того чтобы система (1) была управляемой во множестве допустимых управлений и, необходимо и достаточно, чтобы для любых векторов а е Еп и б е Ек существовало управление и(Г) е и, удовлетворяющее равенству

Т г

| С(Г )Х(Г) | X-1 фВфиф^Г = у, (3)

00 Т Т г

в котором у = б -10(Г)Х(Г^Г -1С(Г)Х(Г)|Х~l(£>)у(£>)d£>dt.

0 0 0

Определим множество допустимых управлений равенством

U = { u(t) = S(t )x + R(t )v}, (4)

где S(t) - т х «-матрица, непрерывная на сегменте [0, Т]; Щі) = г ^ )т,

к ___________________________________________________________

Гу (і) = ^ уФА, (t), г(^ - действительные числа, при любом Хєі, к (t) -

Х=1

известные функции, определенные и кусочно-непрерывные на сегменте [0, Т], V - «-мерный постоянный вектор.

Тогда система (1) преобразуется в систему

X = Н(1)х + B(t )Щ0 v + у(1), (5)

в которой

Н(0 = A(t) + ВО1 )S(t). (6)

Для простоты записи для фундаментальной матрицы системы X = Н(1)х сохраним прежние обозначения. Тогда равенство (3) примет вид

Т і

| 0(1)Х(1) | X-1 ©В©Я©^ Л V = у. (7)

0 0

Из равенства (7) вытекает следующая теорема.

Теорема 1. Система (1) тогда и только тогда управляема во множестве допустимых управление и, определенном равенством (4), когда существует матрица Я(|),

Т і

удовлетворяющая неравенству det Р Ф 0, где Р = | 0(1)Х(1)| X" -1©в©я©^.

00

Пусть г - тпк-мерный вектор, определенный равенством г = (гі, г2, ...,

гт, ^, гтп, ^, гтпкX в котором Г1 = Г111, г2 = г211, •••, гт = ^тЪ •••, гтп = ГПт, •••,

гктп = гпт.

Непосредственным вычислением устанавливаем, что

Л* Р = X а(1, І2, к, іп )г1 г2 к \ ,

(іЬ i2,к, іп)Є-0

где В - множество сочетаний из натуральных чисел 1, кшп по п.

Теорема 2. Для того чтобы система (1) была управляемой на сегменте [0, Т] во множестве допустимых управлений, определенных равенством (4), необходимо и достаточно, чтобы существовало хотя бы одно сочетание (-1, -2, ..., 1п) е В, при

К°ТОр°м Я(/1, /2, к, /п ) ^ ^

Доказательство. Необходимость. Так как система (1) управляема во множестве допустимых управлений, определенных равенством (4), то согласно теореме 1 существует матрица И(§) такая, что выполненяется неравенство detР Ф 0. Это значит, что существует сочетание (/1, /'2, ..., /п) е В, удовлетворяющее неравенству г к, /п) Ф 0.

Достаточность. Пусть сочетание (-1, -2,к,гп) е В таково, что выполняется неравенство а^,-2,.,-п) Ф0. Тогда, положив гг1 = г/2 =... = г-п = 1, г. = 0 при

у е {-1, -2,., гп }, получим матрицу Я* (/), удовлетворяющую неравенству

Т I

det J G(t)X(t)J X-1©B©R* ©d^dt

ф 0.

Следовательно, на основании теоремы 1 система (1) управляема на сегменте [0, Т] во множестве управлений, определенном равенством (4).

Теорема доказана.

Согласно теореме 2, для решения проблемы управляемости конкретных линейных систем необходимо иметь явное представление фундаментальной матрицы системы х = Н(/ )х.

Рассмотрим систему уравнений

Н(0 = Л(0 + Б8(/), (8)

в котором Л(/), Н(/) - известные п х п-матрицы; Б - известная постоянная ненулевая п х п-матрица; 8(/) - искомая матрица.

Будем предполагать, что матрицы Л(/), Н(/) представимы равенствами

I I

Л(/)=X Л1(1 ^ Н(/)= X Н Л (/), в которых при любых це1,1 Лц, Нц -

Ц=1 Ц=1

постоянные матрицы; (/) - некоторая известная функция, непрерывная на сег-

I

менте [0, Т ]. Тогда матрицу 8(/) будем искать в виде 8(/) = X 8Л (/).

Ц=1

Для того чтобы выполнялось равенство (8) достаточно, чтобы при любом це1,1 существовала матрица 8 ц, для которой выполняется равенство

Нц = Лц + Б8ц. (9)

Цустъ Бц =(Б,НЦ- Лц).

Теорема 3. Если при любом ц е 1,1 матрица Н ц такова, что

га^(Б) = га^(Бц), то существует матрица 8(/), удовлетворяющая равенству (8).

Доказательство теоремы следует из теоремы Кронекера-Капелли. Теорема 4. Пусть ш = п -1. Если любой минор (п - 1)-го порядка матрицы Б не равен нулю, то существуют матрица 8ц и диагональная матрица Нц, удовлетворяющие равенству (9).

Доказательство. Пусть ц - произвольное фиксированное число,

це1,1. Рассмотрим систему уравнений: (со1оп(&21,...,Ьп{), со1оп(^2,Ьп2),..,

со1оп(Й2ш, ьпш)) со1оп(.(ц), з^, ..., Ш) = со1оп(з - , 0 - а^,..., 0 - ).

Так как всякий минор (п - 1)-го порядка матрицы Б не равен нулю, то можно найти вектор со1оп(.[1),я^,...,.Щ?), удовлетворяющий системе (9). Аналогично

можно отыскать остальные столбцы матрицы 8 ц. Из равенства (9) определим диагональные элементы матрицы Нц. Из произвольности ц следует справедливость теоремы.

Теорема доказана.

Теорема 5. Если при любом t е [0, Т] det Б(/) Ф 0, то для любых матриц Л(/), Н(0 существует матрица 8(0, удовлетворяющая равенству (6). Доказательство очевидно.

Следовательно, теоремы 3-5 определяют условия, при которых матрица Н^) может быть выбрана диагональной, а фундаментальная матрица системы

х = Н(0х записана в явном виде.

Заключение. Таким образом, получены необходимые и достаточные условия управляемости модели (1) на сегменте [0, Т] .

Установлено посредством теорем 1-5, что кредитные вложения в развитие фондов можно выбрать так, что за промежуток времени [0, Т] будет достигнут

заранее заданный уровень потребления.

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

х = Л(0х + Би, (10)

~ w ч I , I 2 „ 1 1 1 1 ) , (

в которой A(t) =1 colonl--------, - 3----------------, 1 +------1--I, colon

1 4 t(t+2) t + 2 t(t+2) t + 2 t(t+2)' ^

-2 +

ч

3 2 ,3 5.5 1 0.5 \ л L 4 1

- 6 +----------------------,----------------------I, colonl 1 --

t + 2 t (t + 2) t + 2 t (t + 2) t + 2 t(t + 2) J Ч t + 2 t (t + 2)

10 - +------------5—,3 +—5-9— ||; B = (colon(1, 2,0), colon(0,1, -1)); u(t)= (u1(t),

t + 2 t (t + 2) t + 2 t ^ + 2)^

42^)), где Щ^) - кредитные вложения собственных средств, ^(0 - вложения внешних заимствований.

Выпуск продукции z(t), идущей на потребление в момент времени t е [0,1],

определим равенством z(t) = О(0х^), в котором ) = (со1оп(1, 0, 0),

со1оп(0, t, 0), со1оп(0,0, t(t + 2))) - матрица, определяющая часть объема фонда х(0, идущего на потребление продукции. Уровень потребления за период времени [0,1] зададим равенством

11

| z(t)dt =| )x(t)dt = б, (11)

0 0

где б - 3-мерный постоянный вектор.

Пусть б = со1оп(1,1,1), а = со1оп(0,1,0.5). Найдем условия существования управления и^), при котором математическая модель (10) имеет решение, удовлетворяющее равенству (11).

Для решения поставленной задачи достаточно найти условия управляемости системы (10) на сегменте [0,1].

+

Определим множество допустимых управлений равенством (4), в котором 8(0 - 2х3-матрица, непрерывная на сегменте [0,1]; Я(0 = ^ + Я2t + Rзt2,

Я2, Я3 - постоянные 2х3-матрицы, Я ^ =Г? > Г , V - 3-мерный постоянный вектор.

Тогда система (10) преобразуется в систему

х = H(t )х + БR(t) V. (12)

Поскольку для матрицы А^) выполняется равенство А^) = А1 + А2 —+

+ А3 —1—, то матрицу Н^) возьмем в виде Н^) = Н + Н 2 —1—+ Н3 —1—, 3 ^ + 2) 1 2 t + 2 3 t(t + 2)

где Н1 = diag(1, - 2,1), Н2 = diag(0, - 2, 3), Н1 = diag(2, - 2, - 2). Пусть Бц = (Б,Нц - Ац), це1, 3.

Можно убедиться, что выполняются условия теоремы 3. Это значит, что существует матрица 8(0, удовлетворяющая равенству (6), и управление можно искать в виде и(0 = 8(0х + )у.

Фундаментальная матрица Х(/) системы х = H(t)x примет вид

= diag(exp(t - 1п^ +1) + 1п(0), ехр(- 2t - 1n(t2 + 2t)), ехр(1 + 1п^ + 2) - 1п^))).

Из определения матрицы ) следует, что г = (Г1, Г2, к, Гб, к, Г(8), в кото-

(1) (1) (1) (2) (2) (3)

ром Г1 = ГlУ, Г2 = Г\—, • ••, Г6 = Г2з , Г7 = Г11 , Г8 = Г21 % •••, Г18 = Г23

Учитывая, что в рассматриваемом случае

det Р = Е °(П, <2, <3)Г1 Г 2 Г<3>

(^1, <2, <3)еБ

где Б - множество сочетаний из натуральных чисел 1,18 по 3, непосредственным вычислением устанавливаем, что detР = -1.332^7^7 +.... Поскольку а(4 7 17) = -1.332, то, положив Г4 = Г7 = Г17 = 1, остальные Г< = 0 при < й {4, 7,17}, получим, что detР = -1.332 при Я! = (со1оп(0,0), со1оп(0,1), со1оп(0, 0)), Я2 = (со1оп(1, 0), со1оп(0,0), со1оп(0, 0)), Яз = (со1оп(0, 0), со1оп(0,0), со1оп(1, 0)), при этом матрица Я*^) определяется равенством Я*^) = (co1on(t, 0),

со1оп(0,1), co1on(t2, 0)).

Учитывая, что б = со1оп(1,1,1), а = со1оп(0,1, 0.5), из равенства Ру = у будем иметь V = со1оп(7.439, - 0.92, - 5.111). Следовательно, управление, при котором система (10) имеет решение, удовлетворяющее равенству (11), определяется соотношением и^) = 8(0х(0 + Я* (0 у, где х(0 (х(0) = а) - решение системы (12), а значит и системы (10) при условии, что ) = Я* ^).

Список литературы

1. Максимов, В.П. О некоторых обобщениях некоторых дифференциальных уравнений, краевых задач и их приложения к задачам экономической динамики / В.П. Максимов // Вестн. Перм. техн. гос. ун-та. - 1997. - № 4. - С. 103-120.

Research into the Controlled Mathematical Model of the Multi-Sector Economy with the Constant Level of Consumption

I.S. Potapova

Department «Mathematics and Teaching Methods of Mathematical Disciplines», Ryazan State University Named after S.A. Esenin;

Irina00000@yandex.ru

Key words and phrases: acceptable control; external influence; fundamental matrix; provision of crediting; set of combinations; size of main assets; system of the differential equations; vector function.

Abstract: The mathematical model of change of the capital stock size is constructed. The necessary and sufficient conditions for control when the model has the solution, defining development of the branches of economy with the given level of consumption are specified.

Untersuchung des steuernden matematischen Modells der mehrsektoralen Wirtschaft mit dem standigen Verbrauchsniveau

Zusammenfassung: Es ist das matematische Modell der Veranderung der Volumen der Hauptbetriebsfonds gebaut. Es sind die notigen und genugende Bedingungen des Existierens der Steuerung, bei dem das Modell den Beschluft hat, der die Entwicklung der Wirschaftszweige mit dem fruher vorgegebenen Verbrauchsniveau bestimmt, gefunden.

Etude du modele mathematique controle de l’economie multisectorielle avec un niveau constant de la consommation

Resume: Est construit le modele mathematique des changements des volumes des fonds principaux de production. Sont trouvees les conditions necessaires et admissibles de la gestion lors de laquelle le modele a une solution determinant le developpement des branches de l’economie avec un niveau de la consommation calcule d’avance.

Автор: Потапова Ирина Сергеевна - аспирант кафедры математики и методики преподавания математических дисциплин, ГОУ ВПО «РГУ им. С.А. Есенина».

Рецензент: Терехин Михаил Тихонович - доктор физико-математических наук, профессор кафедры математики и методики преподавания математических дисциплин, ГОУ ВПО «РГУ им. С. А. Есенина».