Научная статья на тему 'ПАРАБОЛИЧЕСКАЯ МОДЕЛЬ В ФОРМЕ ПРОСТРАНСТВА СОСТОЯНИЙ ДИНАМИКИ НАКОПЛЕНИЙ'

ПАРАБОЛИЧЕСКАЯ МОДЕЛЬ В ФОРМЕ ПРОСТРАНСТВА СОСТОЯНИЙ ДИНАМИКИ НАКОПЛЕНИЙ Текст научной статьи по специальности «Математика»

CC BY
60
5
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ / СТОХАСТИЧЕСКАЯ ДИНАМИЧЕСКАЯ СИСТЕМА / ВХОДНЫЕ И ВЫХОДНЫЕ ДАННЫЕ НАБЛЮДЕНИЙ / ПАССИВНАЯ ИДЕНТИФИКАЦИЯ / АДДИТИВНЫЙ БЕЛЫЙ ШУМ / ЛИНЕЙНОЕ КОНЕЧНО-РАЗНОСТНОЕ УРАВНЕНИЕ / УРАВНЕНИЕ ПАРАБОЛИЧЕСКОГО ТИПА / МОДЕЛЬ В ПРОСТРАНСТВЕ СОСТОЯНИЙ / ФИЛЬТР КАЛМАНА

Аннотация научной статьи по математике, автор научной работы — Абденова Гаухар Амирзаевна, Базикова Карлыгаш Манаповна, Кенжегалым Жангул Медеткызы

Важное место в теории уравнений с частными производными и ее приложениях занимает уравнение теплопроводности, представитель класса так называемых параболических уравнений. Известно, что для проверки корректности математической модели, основанной на параболическом уравнении, очень важно существование ее решения, так как математическая модель не всегда адекватна конкретному явлению, и из существования решения реальной прикладной задачи не следует существование решения соответствующей математической задачи. Поэтому методы решения дифференциальных уравнений в частных производных - как аналитических, так и численных - всегда актуальны. В настоящее время вычислительный эксперимент стал мощным средством теоретических исследований. Он проводится над математической моделью исследуемого объекта, но при этом по одним параметрам модели вычисляются другие параметры и делаются выводы о свойствах изучаемого объекта или явления. В работе рассматривается задача пассивной параметрической идентификации систем с распределенными параметрами для динамики аккумуляции ресурсов множества домохозяйств, а также преобразование ее в модель в форме пространства состояний с учетом белых шумов модели динамики исследуемого объекта и измерительной системы линейного типа. Используя метод конечных разностей, решение уравнений с частными производными параболического типа свели к решению системы линейных конечно-разностных и алгебраических уравнений. Для более точного оценивания поведения объекта было предложено использование алгоритма фильтрации по схеме Калмана. Осуществлены расчеты с помощью математической системы Matlab на основе данных наблюдений за пять лет, взятых с сайта «Бюро национальной статистики Агентства по стратегическому планированию и реформам Республики Казахстан». Оценивание коэффициентов уравнений аккумуляции ресурсов домохозяйств в форме пространства состояний с использованием данной методики в достаточной степени универсально и может быть применено и в других областях науки и техники.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Абденова Гаухар Амирзаевна, Базикова Карлыгаш Манаповна, Кенжегалым Жангул Медеткызы

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

A PARABOLIC MODEL IN THE FORM OF SPACE STATES OF THE DYNAMICS OF SAVINGS

An important place in the theory of partial differential equations and its applications is occupied by the heat equation, a representative of the class of the so-called parabolic equations. It is known that to check the correctness of a mathematical model based on a parabolic equation, the existence of its solution is very important since a mathematical model is not always adequate to a specific phenomenon and the existence of a solution to a corresponding mathematical problem does not follow from the existence of a solution to a real applied problem. Therefore, methods for solving partial differential equations, both analytical and numerical, are always relevant. Nowadays, a computational experiment has become a powerful tool for theoretical research. It is carried out over a mathematical model of the object under study, but at the same time, other parameters are calculated using one of the parameters of the model and conclusions are drawn about the properties of the object or phenomenon under study. The problem of passive parametric identification of systems with distributed parameters for resource accumulation dynamics of many households using a stochastic distributed model in the form of a state space with regard to the white noise of the dynamics model of the object under study and the white noise of the model of a linear-type measuring system is considered in the paper. The use of the finite difference method allowed us to reduce the solution of partial differential equations of a parabolic type to the solution of a system of linear finite difference and algebraic equations represented by models in the form of a state space. It was also proposed to use a filtering algorithm based on the Kalman scheme for reliable estimation of the object behavior. Calculations were carried out using the Matlab mathematical system based on statistical data for five years, taken from the site “Agency for Strategic Planning and Reforms of the Republic of Kazakhstan Bureau of National Statistics”. Estimation of the coefficients of the equations for the household resource accumulation in the form of a state space using this technique is sufficiently universal and can be applied in other fields of science and technology.

Текст научной работы на тему «ПАРАБОЛИЧЕСКАЯ МОДЕЛЬ В ФОРМЕ ПРОСТРАНСТВА СОСТОЯНИЙ ДИНАМИКИ НАКОПЛЕНИЙ»

ISSN 2782-2001 Системы анализа и обработки данных том 82, № 2, 2021, с. 7-18

http://journals.nstu.ru/vestnik Analysis and data processing systems Vol. 82, No. 2, 2021, pp. 7-18

ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И УПРАВЛЕНИЕ

INFORMATICS, COMPPUTER ENGINEERING AND CONTROL

УДК 517.95: 519.24 Б01: 10.17212/2782-2001-2021-2-7-18

Параболическая модель в форме пространства

*

состояний динамики накоплений

ГА. АБДЕНОВА", КМ. БАЗИКОВА4, Ж.М. КЕНЖЕГАЛЫМС

010008, Республика Казахстан, г. Нур-Султан, ул. Сатпаева, 2, Евразийский национальный университет имени Л.Н. Гумилева

а gauhar.phd@gmail.com ь kmbazikova@mail.ru с zhangul_flower@mail.ru

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

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

Осуществлены расчеты с помощью математической системы МаНаЪ на основе данных наблюдений за пять лет, взятых с сайта «Бюро национальной статистики Агентства по стратегическому планированию и реформам Республики Казахстан». Оценивание коэффициентов уравнений аккумуляции ресурсов домохозяйств в форме пространства состояний с использованием данной методики в достаточной степени универсально и может быть применено и в других областях науки и техники.

Статья получена 18 января 2021 г. Исследование выполнено в Евразийском национальном университете имени Л.Н. Гумилева Республики Казахстан.

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

ВВЕДЕНИЕ

В настоящее время построение динамических моделей для трудно формализуемых областей, использующих дифференциальные уравнения с частными производными [1] и учитывающих случайные факторы и аддитивные шумы, вызывает особый интерес. Так, для моделирования аккумулирования ресурсов исследуемого объекта будем использовать стохастическую параболическую распределенную модель, описывающую динамику накоплений до-мохозяйств, в виде уравнений с частными производными, а модель измерительной системы - в виде линейных стохастических алгебраических моделей с аддитивными белыми шумами.

Анализ современных научных исследований показывает эффективность решения задач идентификации с использованием моделей пространства состояний (ПС) [2-4]:

й (г +1) = Фй (г) + Бх(г), й (1) = й1, (1)

У (г +1) = Нй (г + 1) + т(г + 1), г - 1, 2,..., (2)

где г - время; й(г) - «-вектор состояния; Ф - переходная матрица состояния, размерность п х п; х(г) - г-вектор управления; Б - переходная матрица управления, размерность п х г; й (1) - «-вектор начального состояния; У (г) - га-вектор выходных данных; Н - матрица наблюдения, размерность т х «; м>(г) - «-вектор белого шума измерительной системы, имеющий нулевое математическое ожидание и неизвестную ковариационную матрицу Я.

1. ПОСТАНОВКА ЗАДАЧИ

Известно, что исследуемые объекты в каждый момент времени распределены неравномерно. В работах [5, 6] на основе применения принципа сплошных сред и введения функции, которая опишет распределение домохо-зяйств по аккумуляции ресурсов, получено параболическое уравнение с частными производными вида (3), которому удовлетворяет плотность исследуемых объектов на пространстве и:

ди д ,, „ ч 1 д2 . ... — = -—((с + Р )и) + -—(Ьи) + #, (3)

дг дх 2 дх2

с начальными условиями

х(г)|г-0 - х0 ,

(4)

где х(:)| : _0 - гауссовская величина с математическим ожиданием X} и дисперсией р.

В уравнении (3) u (х, :) - искомые значения состояний эволюции аккумуляции ресурсов домохозяйств распределенного типа; ^ (х, :) - заданная функция распределенного типа; с, Ь, d - постоянные коэффициенты; / (х, :) - число домохозяйств, которые попадут из других множеств домохо-зяйств на отрезок единичной длины пространства накоплений за единичный интервал времени в окрестности х и :. Уравнение (3) называется дифференциальным уравнением параболического типа, описывающим накопления множества домохозяйств. С учетом аддитивных шумов процесса аккумуляции ресурсов уравнение (3) можно записать в виде

((с + + 1-^т(bu) + df + р^, (5)

от ох 2 7х

где Р1 = Q(Xk, ) - шумовой коэффициент в уравнении состояния; w - аддитивные белые шумы измерителей.

Предположим, что за исследуемыми объектами можно вести наблюдение, поэтому уравнение можно записать в виде

у^, :3) = , :3) + Р2в^, :3),

_ _ (6) k = 1, п, 5 = 1, m,

где u (х, :) - пространственно-временная функция состояния, измеряемая в дискретных точках Xk и : ^ы(х, ?) « u(Xk, ) = Uk 5, k = 1, п, 5 = 1, mj; у (Xk, ) - выходные наблюдения; h - весовой коэффициент; |в(Xk, ) = вks, k = 1, п, 5 = 1, mj - белый шум системы с неизвестной ко-

вариацией р2 = Л(Xk, ).

Учитывая указанные условия на основе дискретного выхода наблюдаемой системы у (Xk, ) распределенного типа и начальное условие (4), можно сформулировать задачу оценивания параметров с, Ь, d, Р1, Р2 .

2. АЛГОРИТМ ОЦЕНИВАНИЯ ПАРАМЕТРОВ

В ПАРАБОЛИЧЕСКОМ УРАВНЕНИИ АККУМУЛЯЦИИ РЕСУРСОВ

Для реализации решения задачи пассивной идентификации используем дискретную форму выражения (5) в частных производных [7, 8]:

8u 8 ,82 ...

— = —[—u + bl—2u + ^ + Рlw, (7)

8: 8х 8х2

где «1 = с + ^, ¿1 = Ь /2, ^ = й, Р1 - коэффициенты; ? - время, ? > 0; х - координата, 0 < х < Ь; w(х, ?) - белый шум с неизвестной ковариацией Q и нулевым математическим ожиданием; и (х, ?) - состояние эволюции аккумуляции ресурсов домохозяйств распределенного типа с граничным и начальным условиями:

ди (х, t)

дх

или

и (х t )1 х=0 =

(8а)

х=0

и (х, t)1х=Ь = иЬ , и(х, t)|?=0 = и(х, 0) = и0(х) .

(8б) (9)

Для ? и х интервалы квантования приравняем: Д? = 1 и Дх = 1, получим следующее выражение:

1 ¿1

(ик,5 +1 " ик5 ) =-Т(ик,5 +1 " 2ик,5 + ик5-1) "

Дх2

Д?

- ~Г(ик,5 +1 - ик,5 ) + й1/к,5 + Р^к,5 ,

Дх

где к = 1, п, 5 = 1, т.

Значения функций и (х, ?) объединим по координатам к, 5 узлов сетки при интервалах Д? = 1 и Дх = 1.

Получим систему уравнений при к = 1, 2,..., п :

и1,5 +1 = с1и0,5 - с2и1,5 + с3и2,5 + и1,5 + ¿1/1,5 + Р1Щ,5 , и2,5+1 = с1и1,5 - с2и2,5 + с3и3,5 + и2,5 + й1/2,5 + Рlw2,

ип,5+1 = с1ип,5 - с2ип,5 + с3ип +1,5 + и1,5 + й1/П,5 +

(10)

Сделаем следующие обозначения:

( и1,5+1 ^ и2,5+1

V ип,5+1 Упх1

( и1,5 Л

и2,5

и

V^ Упх1

/1

/2.

V/n,s Упх1

ж =

( П 5 ^ w2,s

V у

( С1 >

Е2 =

пх1

V 0 у

( 0 >

Во =

пх1

с3,

пх1

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

05 = =0; л* = ч

5 _ П+1,5 ■

(11)

(12)

Уравнение (7) в векторно-матричной форме примет вид (-с2 с3 0 0 ••• 0 0 0 ^ с1 - с2 с3 0 ••• 0 0 0

А

м=

0 0 0 0 0 0 0 0

С1 - с2 с3

0 С1 - С2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(13)

^пхп

ОА = d1E; ПА = Р1Е,

(14)

где Е - единичная матрица размерностью п х п. Используя (11)—(14), выражение (10) представим в векторно-матричной форме:

^+1 = ми8 + аи8 + НЕ + ^ + (Е205 + В2Ц3 ).

(15)

В [9, 10] приведены обозначения, которые позволяют представить соотношение (15) в виде модели в ПС.

В итоге уравнение (10) запишется в виде

и5+1 = ми5 + ОН5 + ИЕ + .

(16)

Соотношение (16) по виду соответствует выражению (1), где время может иметь значения: 5 = {1, 2,..., m}, а начальные условия можно записать следующим образом:

и0 =

(щл ^

^,1

(^,0 ^

^,0

(17)

V пЛ У V п,° Упх1

А

Далее, при условии h = const, введем обозначения для уравнения (6):

Ys =

'Л, ^

У2,.

Vy",s Jnx1

Us =

(u1,s ^

u2,

u

Vйn,s Уих1

ь1, s S2,s

V n,s Jnxl

Y0 = У0 / h •

В итоге уравнение (6) запишем в векторном виде

Ys = hUs + p2Ss , s = 0, 1, 2,..., m .

(18)

(19)

Полученные соотношения (16), (17), (19) - модель стохастической системы с распределенными параметрами, представленные как модель ПС (1)-(2).

При значениях п и 5 сгруппируем в (16) дискретные данные наблюдений ик j за четыре года и выведем четыре уравнения с неизвестными

сЪ -c2, c3, ¿1:

un-4,s +1 = c1un-5,s - c2un-4,s + c3un-3,s +

+ un-4,s + d1fn-4,s + p1wn-4,s;

un-3,s+1 = c1un-4,s - c2un-3,s + c3un-2,s +

+ un-3,s + d1fn-3,s + p1wn-3,s;

un-2,s+1 = c1un-3,s - c2un-2,s + c3un-1,s +

+ un-2,s + d1fn-2,s + p1 wn-2,s;

un-1,s+1 = c1un-2,s - c2un-1,s + c3un,s +

+ un-1,s + d1 fn-U + P1 wn-1„

(20)

В уравнениях (16) и (19) Р1, Р2 есть коэффициенты при шумовых неконтролируемых переменных Ж и в, которые дают количественную характеристику дисперсии аддитивных шумов динамики и измерительной системы, являются белыми гауссовскими последовательностями с нулевым средним и дисперсиями, равными единице. В работах [11-13] предлагается использование алгоритма расчета всех дисперсионных характеристик модели динамики

£ s =

и измерительной системы с учетом совокупности выходных данных {Ук,5, к = 1,«, 5 = 1,от).

Можно найти средние оценки дисперсий для Р1 « Q = Q и Р2 ® Щ = Щ. Данный алгоритм можно использовать для каждой строки к = 1, п . Оценки дисперсий используются при параметрической идентификации состояния объекта с применением калмановской фильтрации [2, 11, 13, 14].

В дальнейшем для расчетов оценок коэффициентов С1, - С2, С3, ^ можно использовать уравнения фильтра Калмана. Если рассматривать систему (21) как систему четырех уравнений с четырьмя неизвестными, то можно рассчитать четверку оценок С1, - С2, С3, ^ [15] на основе следующих формул:

а1 =

Г и«-1,5+1 - и«-1,5 и«-2,5 +1 - и«-2,5 1 /Г и«-2,5 - 2ип-1,5 + ип

1«-1,5

ип-3,5 - 2ип-2,5 + ип-1,5

П—1,5

(21)

ь =

*п-2,5+1 ип-2,5

Лп-1,5

«1

ип-3,5 - 2ип-2,5 + ип-1,5

Лп-1,5

(22)

3. ЧИСЛЕННЫЕ РЕЗУЛЬТАТЫ АЛГОРИТМА

ИДЕНТИФИКАЦИИ КОЭФФИЦИЕНТОВ УРАВНЕНИЯ

Для примера возьмем данные аккумуляции ресурсов домохозяйств за четыре года (2014-2017 гг.).

Пример.

а) Чтобы осуществить расчеты для проведения исследований на основе вычисления коэффициентов уравнений накоплений домохозяйств и в последующем для получения наиболее достоверных оценок поведения исследуемого объекта, по схеме фильтра Калмана промоделируем значения при следующих значениях индексов узлов сетки модели: п = {1, 2, 3, 4, 5, 6}, 5 = {1, 2, 3, 4, 5, 6, 7, 8}. Данный пример рассчитывался в математической системе МайаЬ [16].

б) Задача численного оценивания коэффициентов распределенного уравнения параболического типа аккумуляции ресурсов множества до-мохозяйств при заданных данных наблюдений и краевых, начальных условиях. Получили значения коэффициентов уравнений: с = -0.0006477; ?2 = 1.000606 ; с3 = -0.018023 ; « = 74.6055, использовав выражения (21), (22). Смоделированы значения наблюдений расходов Щ и Щ исследуемого объекта за 2014-2017 гг.

в) Более точно была решена задача оценивания состояния динамики объекта с применением калмановской фильтрации [11, 13]. Для расчета оце-

нок предсказаний и оценок фильтрации использованы данные наблюдений: п = 24, у = [76200 76097 82343 82285 82339 79187 79111 79013 82485 82393 82307 83346 83286 83161 90188 90049 89992 86508 86425 86299 91191 91025 90960 94894]; дисперсия динамики объекта Q = 1282; дисперсия измерительной системы Я = 33177; расчетные коэффициенты дифференциального уравнения: а = 0.7012; Ь = 21997.

Результаты расчета и графики поведений наблюдений, оценок предсказания и фильтрации представлены на рисунке.

хр = [90132 85198 79357 76812 75670 76806 77289]; (Пр) - % оценки предсказания.

х/ = [90132 81803 78173 76544 78164 78853 79176]; (Ф) - % оценки фильтрации.

у^ = [90132 76291 76200 76097 82343 82285 82339]; (Н) - % наблюдения.

Графики наблюдений, оценок предсказания и фильтрации Plots of Observation and Prediction and Filtering Estimates

г) Обратная задача оценивания коэффициентов уравнения аккумуляции ресурсов домохозяйств решена с применением фильтрационных оценок поведения объекта по схеме Калмана. Расчетные оценки С, С2, £3, di были использованы для наиболее точного моделирования выходных данных системы.

ЗАКЛЮЧЕНИЕ

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

СПИСОК ЛИТЕРАТУРЫ

1. Самарский А.А., Гулин А.В. Численные методы математической физики. - М.: Научный мир, 2003. - 316 с.

2. Горский В.Г., Адлер Ю.П., Талалай А.М. Планирование промышленных экспериментов (модели динамики). - М.: Металлургия, 1978. - 112 с.

3. АнисимовА.С. Идентификация объектов управления: учебное пособие. - Новосибирск: НЭТИ, 1985. - 80 с.

4. Mehra R.R. Optimal input signals for parameter estimation in dynamic systems - Servey and new results // IEEE Transactions on Automatic Control. - 1974. - Vol. 19, iss. 6. - P. 753-768. -DOI: 10.1109/TAC.1974.1100701.

5. Ерофеенко В.Т., Козловская И.С. Уравнения с частными производными и математические модели в экономике: курс лекций. - Минск: Изд-во БГУ, 2004. - 245 с.

6. ЧернавскийД.С., ПопковЮ.С., РахимовА.Х. Математические модели типологии семейных накоплений // Экономика и математические методы. - 1994. - Т. 30, вып. 2. -С. 98-106.

7. Араманович И.Г., Левин В.И. Уравнения математической физики. - М.: Наука, 1969. -

288 с.

8. Абденова Г.А. Структурно-параметрическая идентификация систем с распределенными параметрами с использованием модели типа «вход-состояние-выход» // Научный вестник НГТУ. - 2006. - № 1 (38). - С. 9-16.

9. Абденов А.Ж., Абденова Г.А. Методика пассивной идентификации коэффициентов уравнения теплопроводности с учетом ошибок оценок состояния объекта и измерительной системы // Автометрия. - 2016. - Т. 52, № 2. - С. 43-51. - DOI: 10.15372/AUT20160205.

10. Оценивание коэффициентов уравнения теплопроводности с учетом шумов измерительной системы / А.Ж. Абденов, Ж.К. Нурбекова, В.Б. Уткин, Г.А. Абденова // Вестник ЕНУ им. Л.Н. Гумилева. - 2015. - № 4 (107). - С. 5-13.

11. Mehra R.R. Identification and adaptive Kalman filtering // Mechanics. - 1971. - N 3. -P. 34-52.

12. Абденова Г.А., Воевода А.А. Оценивание параметров и характеристик шумов нестационарных процессов в стохастических системах, описываемых в пространстве состояний // Сборник научных трудов НГТУ. - 2010. - № 3 (61). - C. 11-18.

13. Абденова Г.А. Прогнозирование значений уровня временного ряда на основе уравнений фильтра Калмана // Ползуновский вестник. - 2010. - № 2. - С. 4-6.

14. Эйкхофф П. Основы идентификации систем управления. - М.: Мир, 1975. - 680 с.

15. Abdenov A., Abdenova G., KulbaevD. Estimation of equation coefficients of free and forced string vibrations in continuous medium with consideration of dynamic noise and measurement // Materials Today: Proceedings. - 2019. - Vol. 16. - P. 336-342.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

16. ДьяконовВ.П. MATLAB. Полный самоучитель. - М.: ДМК Пресс, 2012. - 768 с.

Абденова Гаухар Амирзаевна, доцент кафедры «Математическое и компьютерное моделирование» механико-математического факультета Евразийского национального университета им. Л.Н. Гумилева, кандидат технических наук, PhD. Область научных интересов: математическое моделирование в экономике, параметрическая идентификация систем. Автор более 40 научных работ. E-mail: gauhar.phd@gmail.com

Базикова Карлыгаш Манаповна, докторант специальности «Математическое и компьютерное моделирование» Евразийского национального университета им. Л.Н. Гумилева. Направление научных интересов - математическое и компьютерное моделирование. Имеет 2 печатные работы. E-mail: kmbazikova@mail.ru

Кенжегалым ЖангYл Медетцызы, магистрант специальности «Математическое и компьютерное моделирование» Евразийского национального университета им. Л.Н. Гу-

милева. Направление научных интересов - математическое и компьютерное моделирование. Научных работ и учебных пособий пока не имеет. E-mail: zhangul_flower@mail.ru

Abdenova Gaukhar A, associate professor at the Department of Mathematical and Computer Modeling, Faculty of Mechanics and Mathematics, L.N. Gumilyov Eurasian National University, PhD (Eng.). Her research interests are focused on mathematical modeling in economics and parametric identification of systems. She has more than 40 scientific publications. E-mail: gauhar.phd@gmail.com

Bazikova Karlygash M., a postdoctoral student specializing in "Mathematical and Computer Modeling", L.N. Gumilyov Eurasian National University. Her research interests are focused on mathematical and computer modeling. She has 2 scientific publications. E-mail: kmbazikova@mail.ru

Kenzhegalym Zhangul M., a master degree student specializing in "Mathematical and Computer Modeling", L.N. Gumilyov Eurasian National University. His research interests include mathematical and computer modeling. He has not published any papers yet. E-mail: zhangul_flower@mail.ru

DOI: 10.17212/2782-2001-2021-2-7-18 A parabolic model in the form of space states of the dynamics of savings

G.A. ABDENOVAa, K.M. BAZIKOVAb, Zh.M. KENZHEGALYM0

L.N. Gumilyov Eurasian National University, 2 Satpayev Street, Nur-Sultan, 010008, Republic of Kazakhstan

a gauhar.phd@gmail.com b kmbazikova@mail.ru c zhangul_flower@mail.ru

Abstract

An important place in the theory of partial differential equations and its applications is occupied by the heat equation, a representative of the class of the so-called parabolic equations. It is known that to check the correctness of a mathematical model based on a parabolic equation, the existence of its solution is very important since a mathematical model is not always adequate to a specific phenomenon and the existence of a solution to a corresponding mathematical problem does not follow from the existence of a solution to a real applied problem. Therefore, methods for solving partial differential equations, both analytical and numerical, are always relevant.

Nowadays, a computational experiment has become a powerful tool for theoretical research. It is carried out over a mathematical model of the object under study, but at the same time, other parameters are calculated using one of the parameters of the model and conclusions are drawn about the properties of the object or phenomenon under study. The problem of passive parametric identification of systems with distributed parameters for resource accumulation dynamics of many households using a stochastic distributed model in the form of a state space with regard to the white noise of the dynamics model of the object under study and the white noise of the model of a linear-type measuring system is considered in the paper. The use of the finite difference method allowed us to reduce the solution of partial differential equations of a parabolic type to the solution of a system of linear finite difference and algebraic equations represented by models in the form of a state space. It was also proposed to use a filtering algorithm based on the Kalman scheme for reliable estimation of the object behavior.

Calculations were carried out using the Matlab mathematical system based on statistical data for five years, taken from the site "Agency for Strategic Planning and Reforms of the Republic of Kazakhstan Bureau of National Statistics". Estimation of the coefficients of the equa-

*

Received 18 January 2021.

tions for the household resource accumulation in the form of a state space using this technique is sufficiently universal and can be applied in other fields of science and technology.

Keywords: partial differential equation, stochastic dynamical system, input and output observational data, passive identification, additive white noise, linear finite-difference equation, parabolic equation, model in the form of a state space, Kalman filter

REFERENCES

1. Samarskii A.A., Gulin A.V. Chislennye metody matematicheskoi fiziki [Numerical methods of mathematical physics]. Moscow, Nauchnyi mir Publ., 2003. 316 p.

2. Gorskii V.G., Adler Yu.P., Talalai A.M. Planirovanie promyshlennykh eksperimentov (modeli dinamiki) [Planning of industrial experiments. Models of dynamics]. Moscow, Metallurgiya Publ., 1978. 112 p.

3. Anisimov A.S. Identifikatsiya ob"ektov upravleniya [Identifying control objects]. Novosibirsk, NETI Publ., 1985. 80 p.

4. Mehra R.R. Optimal input signals for parameter estimation in dynamic systems - Servey and new results. IEEE Transactions on Automatic Control, 1974, vol. 19, iss. 6, pp. 753-768. DOI: 10.1109/TAC.1974.1100701.

5. Erofeenko V.T., Kozlovskaya I.S. Uravneniya s chastnymi proizvodnymi i matematicheskie modeli v ekonomike [Equations with private derivatives and mathematical models in economics]. Minsk, BSU Publ., 2004. 245 p.

6. Chernavskii D.S., Popkov Yu.S., Rakhimov A.Kh. Matematicheskie modeli tipologii se-meinykh nakoplenii [Mathematical models of family savings typology]. Ekonomika i matematicheskie metody = Economics and Mathematical Methods , 1994, vol. 30, iss. 2, pp. 98-106.

7. Aramanovich I.G., Levin V.I. Uravneniya matematicheskoi fiziki [Equations of mathematical physics]. Moscow, Nauka Publ., 1969. 288 p.

8. Abdenova G.A. Strukturno-parametricheskaya identifikatsiya sistem s raspredelennymi s ispol'zovaniem modeli tipa "vkhod-sostoyanie-vykhod" [Structural and parametric identification of the systems with distributed parameters of using the model of the type "entry-state-exit"]. Nauchnyi vestnik Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta = Science bulletin of the Novosibirsk state technical university, 2006, no. 1 (38), pp. 9-160.

9. Abdenov A.Zh., Abdenova G.A. Metodika passivnoi identifikatsii koeffitsientov uravneniya teploprovodnosti s uchetom oshibok otsenok sostoyaniya ob"ekta i izmeritel'noi sistemy [Passive identification of heat equation coefficients with account for errors in estimating the state of the object and measuring system]. Avtometriya = Optoelectronics, Instrumentation and Data Processing, 2016, vol. 52, no. 2, pp. 43-51. DOI: 10.15372/AUT20160205. (In Russian).

10. Abdenov A.Zh., Nurbekova Zh.K., Utkin V.B., Abdenova G.A. Otsenivanie koeffitsientov uravneniya teploprovodnosti s uchetom shumov izmeritel'noi sistemy [Coefficients evaluation of heat equation with account of noise measuring system]. Vestnik ENU im. L.N. Gumileva = Herald of the L.N. GumilyovEurasian National University, 2015, no. 4 (107), pt. 1, pp. 5-13.

11. Mehra R.R. Identification and adaptive Kalman filtering. Mechanics, 1971, no. 3, pp. 34-52.

12. Abdenova G.A., Voevoda A.A. Otsenivanie parametrov i kharakteristik shumov nes-tatsionarnykh protsessov v stokhasticheskikh sistemakh, opisyvaemykh v prostranstve sostoyanii [The parameters estimation and noise characteristics of time-varying processes in the stochastic systems in the state space]. Sbornik nauchnykh trudov Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta = Transaction of scientific papers of the Novosibirsk state technical university, 2010, no. 3 (61), pp. 11-18.

13. Abdenova G.A. Prognozirovanie znachenii urovnya vremennogo ryada na osnove uravnenii fil'tra Kalmana [Prediction of time series level values based on Kalman filter equations]. Polzunovskii vestnik = Polzunov Bulletin, 2010, no. 2, pp. 4-6.

14. Eykhoff P. System identification. London, New York, Wiley-Interscience, 1974 (Russ. ed.: Eikkhoff P. Osnovy identifikatsii sistem upravleniya. Moscow, Mir Publ., 1975. 680 p.).

15. Abdenov A., Abdenova G., Kulbaev D. Estimation of equation coefficients of free and forced string vibrations in continuous medium with consideration of dynamic noise and measurement. Materials Today: Proceedings, 2019, vol. 16, pp. 336-342.

16. D'yakonov V.P. MATLAB. Polnyi samouchitel' [MATLAB. Complete tutorial]. Moscow, DMK Press Publ., 2012. 768 p.

Для цитирования:

Абденова Г.А., Базикова К.М., Кенжегалым Ж.М. Параболическая модель в форме пространства состояний динамики накоплений // Системы анализа и обработки данных. - 2021. -№ 2 (82). - С. 7-18. - DOI: 10.17212/2782-2001-2021-2-7-18.

For citation:

Abdenova G.A., Bazikova K.M., Kenzhegalym Zh.M. Parabolicheskaya model' v forme pros-transtva sostoyanii dinamiki nakoplenii [A parabolic model in the form of space states of the dynamics of savings]. Sistemy analiza i obrabotki dannykh = Analysis and Data Processing Systems, 2021, no. 2 (82), pp. 7-18. DOI: 10.17212/2782-2001-2021-2-7-18.

ISSN2782-2001, http://journals.nstu.ru/vestnik Analysis and data processing systems Vol. 82, No 2, 2021, pp. 7-18

i Надоели баннеры? Вы всегда можете отключить рекламу.