Сопоставление некоторых способов моделирования при описании нелинейной динамики рядами Вольтерра

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

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

В качестве вариантов развития экономики предполагается использовать заложенные сценарные условия, применяемые Министерством экономического развития России. Помимо этого модель позволяет оценить влияние на реальный сектор проводимой Центральным банком России денежно-кредитной политики и изменений в налоговой политике государства.

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

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

БЛАГОДАРНОСТИ: Работа поддержана ЗАО «Прогноз».

Simonov P.M., Shultz M.N. APPLICATION GENERAL EQUILIBRIUM THEORY FOR MODELING RUSSIAN ECONOMY

The possibility of use of the general equilibrium approach for modelling of economy of Russia at its modern stage of development, the basic stages of construction of model, possibility which open at application of the given method, is considered.

Key words: general equilibrium theory; computable general equilibrium model.

УДК 517.968

СОПОСТАВЛЕНИЕ НЕКОТОРЫХ СПОСОБОВ МОДЕЛИРОВАНИЯ ПРИ ОПИСАНИИ НЕЛИНЕЙНОЙ ДИНАМИКИ РЯДАМИ ВОЛЬТЕРРА

© С.В. Солодуша

Ключевые слова: уравнения Вольтерра I рода; система автоматического регулирования. Работа посвящена сопоставлению алгоритмов моделирования нелинейной динамики с помощью квадратичных полиномов Вольтерра. Приведены результаты вычислительных экспериментов для эталонной модели теплообмена, связанных с задачей автоматического регулирования.

2678

В теории математического моделирования нелинейных динамических систем хорошо известен подход, дающий представление отклика у(Ь) системы типа «вход-выход» на входной сигнал х(Ь) = (х1(Ь), х2(Ь), в виде полинома Вольтерра N -ой степени:

N г г

у

^=11^1^... ^0 0

г £ [о, т ], причем у(о) = о, у (г) £ СрТ]. Построить математическую модель в виде (1) — значит решить задачу идентификации ядер Вольтерра Кг1,...,ги, симметричных по тем переменным, которые соответствуют совпадающим индексам *1,..., ги. В монографии [1] разработана методика построения PN {х(г)), основанная на задании многопараметрических семейств кусочно-постоянных тестовых входных сигналов. Эта методика была реализована в программном обеспечении [2], использующем эталонную модель переходного процесса в элементе теплообменного аппарата с независимым подводом тепла:

г * *

ЛлЛ2 I ( Оп \ ( -^1 /'Р(3)й3 -^2 /'Л(з)4з\

АШ = Л—Л;] [А2

0

В (2) Ь — время; Л1 и Л2 — некоторые константы; индексами «0» обозначены параметры начального стационарного режима; А — приращение, например '0(Ь)= Рп + &Ю(Ь).

В работе сравниваются некоторые способы моделирования нелинейной динамики в наиболее важном для приложений случае, когда N = 2 в (1). Исследуется стационарная динамическая система с входом х(Ь) = (А'О(Ь), АQ(t)) и выходом у(Ь) = Агтоа(Ь). Предположим, что задача построения Р2(х(Ь)) решена. На примере эталонной модели (2) рассмотрим задачу автоматического регулирования, связанную с поиском управляющего воздействия х(Ь), поддерживающего сигнал у(Ь) на заданном уровне у* = 0. Выберем в качестве управляющего воздействия сигнал х1(Ь)=АВ(Ь), воздействие х2(Ь) = АQ(t) считаем заданным. В этом случае (1) является полиномиальным интегральным уравнением Вольтерра I рода, непрерывное решение которого, вообще говоря, носит локальный характер. Пусть АQ(t) = = 7QпI(Ь), 7 ^ 75 %.

Один из исследуемых в работе способов связан с использованием уравнения

г г г

У К 1

г г

_)П(Ь — в1 — К)П(Ь — в2 — п п п

г г г

+ ! J К 12(в1,в2)п(Ь — в1 — h)АQ(t — в2)йв1 йв2 = /(Ь) — J К2(в)АQ(t — в)йв,

п п п

возникающего в задаче автоматического регулирования замкнутой нелинейной динамической системы. В (3) К 1(0) = 0, х1 (Ь) = п (Ь — К), п (£) =0, £ £ [—h, 0] , h — известное запаздывание, /(Ь)= е(Ь) — е(Ь — К), сигнал е(Ь) = у* — Аг1тоЛ(Ь), е(£)=0, £ £ [—К, 0], считаем ошибкой управления.

В другом подходе используется полиномиальное уравнение, возникающее в задаче автоматического регулирования динамическим объектом при отсутствии обратной связи: г г г

/т,в)Щв)Ь + ! J К 11

п п п

2679

Ь Ь Ь

+ У У к 12(81, в2)АВ(Ь — в1)АЯ(г — 82^8^82 = А%2тоЛ(Ь) — ^ в)АЯ(в)йв, Ь € [0,Т],

0 0 о

где К\(0, 0) = 0. В основе (4) — комбинация линейной нестационарной и билинейной стационарной составляющих. Расчетные схемы решения (3), (4), приведены в [3], [4] соответственно. Пусть Ьг = гН, г = 1,п, пН = Т. Сеточные решения уравнений (3), (4) обозначим через А'Р^(Ьг), А'и^(Ьг). Ядра Вольтерра в (3), (4) были настроены на сигналы с амплитудами |а| ^ 60%. Вычисления проводились с помощью программного обеспечения [2]. Было получено, что специфика (3), (4) проявляется при входных воздействиях АQ(t) с 7 ^ 71тах и 7 ^ 72тах соответственно, что означает возможную потерю управляемости. Значения 71тах и 72тах (в %), при которых существуют АТ'^^г) и АВ^^) для 0 ^ и ^ Т — Н, приведены в табл. 1. Из таблицы видно, что существуют такие Т, |а|, при которых 71тах и 72тах совпадают с точностью 8 = 10_1. Численный эксперимент показал также важность выбора тестовых сигналов, используемых для идентификации ядер Вольтерра. Анализ исследуемых подходов позволил выделить области предпочтительности того или иного алгоритма для эталонной модели (2).

Таблица 1: Результаты вычислительного эксперимента.

Н (в %) 40 45 50 55 60

T (с) 7lmax 72max 71max 72max 71max 72max 71max 72max 71max 72max

30 55,62 55,59 56,17 56,12 57,23 57,17 58,88 58,83 61,30 61,28

35 54,27 54,25 54,94 54,98 56,16 56,12 58,03 58,00 60,78 60,76

40 53,40 53,38 54,15 54,13 55,49 55,47 57,52 57,51 60,51 60,50

45 52,80 52,78 53,61 53,60 55,04 55,03 57,21 57,20 60,37 60,36

50 52,36 52,35 53,24 53,23 54,74 54,73 57,00 56,99 60,28 60,28

ЛИТЕРАТУРА

1. Апарцин А.С. Неклассические уравнения Вольтерра I рода: теория и численные методы. Новосибирск: Наука, 1999.

2. Солодуша С.В. Программное обеспечение и алгоритмы для моделирования нелинейной динамики полиномами Вольтерра // Программные продукты и системы. 2012. № 4. С. 137-141.

3. Солодуша С.В. Численное моделирование нелинейных динамических систем с векторным входом квадратичными полиномами Вольтерра // Вестн. Южно-Урал. гос. ун-та. Серия: Математика. Механика. Физика. 2012. № 34. С. 53-59.

4. Солодуша С.В. Численное моделирование динамики теплообмена модифицированным квадратичным полиномом Вольтерра // Вычислительные технологии. 2013. Т.18. № 2. С. 83-93.

БЛАГОДАРНОСТИ: Работа поддержана грантом РФФИ, проект № 12-01-00722-а.

Solodusha S.V. COMPARISON OF SOME MODELING APPROACHES FOR DESCRIPTION OF NONLINEAR DYNAMICS BY VOLTERRA SERIES

A comparison of methods of modeling nonlinear dynamics using quadratic Volterra polynomials is considered. The results of computational experiments for reference models of heat exchange associated with the of automatic control problem are given.

Key words: Volterra equations of the first kind; system of automatic control.

2680