Научная статья на тему 'ПОСТРОЕНИЕ НА ОСНОВЕ ИНТЕРПОЛЯЦИИ МОДЕЛЕЙ РАЗЛИЧНЫХ ФИЗИЧЕСКИХ И ХИМИЧЕСКИХ СИСТЕМ МЕТОДОМ МАТЕМАТИЧЕСКОГО ПРОТОТИПИРОВАНИЯ ЭНЕРГЕТИЧЕСКИХ ПРОЦЕССОВ'

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

CC BY
4
2
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
метод математического прототипирования энергетических процессов / цифровые двойники / интерполяция / линейная идентификация / машинное обучение / method of mathematical prototyping of energy processes / digital twins / interpolation / linear identification / machine learning

Аннотация научной статьи по математике, автор научной работы — Игорь Евгеньевич Старостин

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

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

Похожие темы научных работ по математике , автор научной работы — Игорь Евгеньевич Старостин

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

BUILDING, BASED ON INTERPOLATION, MODELS OF VARIOUS PHYSICAL AND CHEMICAL SYSTEMS BY METHOD OF MATHEMATICAL PROTOTYPING OF ENERGY PROCESSES

Background. Solving practical problems (design and operation of systems) implies the construction of mathematical models of systems. For mathematical modeling of systems of various physical and chemical natures, the authors proposed a method of mathematical prototyping of energy processes, based on modern nonequilibrium thermodynamics, mechanics and electrodynamics. The mentioned method provides adequate models of systems, i.e. that do not contradict general physical laws, as well as the peculiarities of processes in a particular system. However, to determine the controlled parameters of a system from its measured parameters from a system of differential equations obtained by the method of mathematical prototyping, it is necessary to solve the very labor-intensive task of identifying a large number of parameters of these equations. One of the ways to combat the complexity of the mentioned identification problems is the use of interpolation methods, which determines the relevance of the task of developing a methodology for constructing system models by mathematical prototyping of energy processes using interpolation methods. Matherials and methods. To synthesize equations for the dynamics of physical and chemical processes, the method of mathematical prototyping of energy processes is used. In order to simplify the identification of the parameters of equations obtained by the method of mathematical prototyping, special methods for integrating systems of differential equations are used, reducing the integration of differential equations to solving algebraic equations. In order to simplify the solution of the resulting algebraic equations, interpolation methods are used. To determine the constant coefficients of the resulting model from experimental data, methods of identification theory are used. Results. The system model constructed by the methods proposed in this article is correct, i.e. does not contradict general physical laws and incorporates the peculiarities of the processes in a particular system. Also, model training is block-by-block and can be reduced to linear identification, which makes it possible to process experimental data as it arrives. Conclusions. The proposed architecture of transformed system models makes it possible to use the resulting models as part of the mathematical core of digital twins. The architectural features of the proposed models make it possible to passively identify them.

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

ПРОЕКТИРОВАНИЕ И ТЕХНОЛОГИИ ПРИБОРОСТРОЕНИЯ И РАДИОЭЛЕКТРОННОЙ АППАРАТУРЫ

DESIGN AND TECHNOLOGY OF INSTRUMENTATION AND ELECTRONIC EQUIPMENT

УДК 519.622.2

doi: 10.21685/2307-4205-2024-1-6

ПОСТРОЕНИЕ НА ОСНОВЕ ИНТЕРПОЛЯЦИИ МОДЕЛЕЙ РАЗЛИЧНЫХ ФИЗИЧЕСКИХ И ХИМИЧЕСКИХ СИСТЕМ МЕТОДОМ МАТЕМАТИЧЕСКОГО ПРОТОТИПИРОВАНИЯ ЭНЕРГЕТИЧЕСКИХ ПРОЦЕССОВ

И. Е. Старостин

Московский государственный технический университет гражданской авиации, Москва, Россия

starostinigo@yandex.ru

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

© Старостин И. Е., 2024. Контент доступен по лицензии Creative Commons Attribution 4.0 License / This work is licensed under a Creative Commons Attribution 4.0 License.

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

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

Для цитирования: Старостин И. Е. Построение на основе интерполяции моделей различных физических и химических систем методом математического прототипирования энергетических процессов // Надежность и качество сложных систем. 2024. № 1. С. 49-58. doi: 10.21685/2307-4205-2024-1-6

BUILDING, BASED ON INTERPOLATION, MODELS OF VARIOUS PHYSICAL AND CHEMICAL SYSTEMS BY METHOD OF MATHEMATICAL PROTOTYPING OF ENERGY PROCESSES

I.E. Starostin

Moscow State Technical University of Civil Aviation, Moscow, Russia starostinigo@yandex.ru

Abstract. Background. Solving practical problems (design and operation of systems) implies the construction of mathematical models of systems. For mathematical modeling of systems of various physical and chemical natures, the authors proposed a method of mathematical prototyping of energy processes, based on modern nonequilibrium thermodynamics, mechanics and electrodynamics. The mentioned method provides adequate models of systems, i.e. that do not contradict general physical laws, as well as the peculiarities of processes in a particular system. However, to determine the controlled parameters of a system from its measured parameters from a system of differential equations obtained by the method of mathematical prototyping, it is necessary to solve the very labor-intensive task of identifying a large number of parameters of these equations. One of the ways to combat the complexity of the mentioned identification problems is the use of interpolation methods, which determines the relevance of the task of developing a methodology for constructing system models by mathematical prototyping of energy processes using interpolation methods. Matheri-als and methods. To synthesize equations for the dynamics of physical and chemical processes, the method of mathematical prototyping of energy processes is used. In order to simplify the identification of the parameters of equations obtained by the method of mathematical prototyping, special methods for integrating systems of differential equations are used, reducing the integration of differential equations to solving algebraic equations. In order to simplify the solution of the resulting algebraic equations, interpolation methods are used. To determine the constant coefficients of the resulting model from experimental data, methods of identification theory are used. Results. The system model constructed by the methods proposed in this article is correct, i.e. does not contradict general physical laws and incorporates the peculiarities of the processes in a particular system. Also, model training is block-by-block and can be reduced to linear identification, which makes it possible to process experimental data as it arrives. Conclusions. The proposed architecture of transformed system models makes it possible to use the resulting models as part of the mathematical core of digital twins. The architectural features of the proposed models make it possible to passively identify them.

Keywords: method of mathematical prototyping of energy processes, digital twins, interpolation, linear identification, machine learning

For citation: Starostin I.E. Building, based on interpolation, models of various physical and chemical systems by method of mathematical prototyping of energy processes. Nadezhnost' i kachestvo slozhnykh sistem = Reliability and quality of complex systems. 2024;(1):49-58. (In Russ.). doi: 10.21685/2307-4205-2024-1-6

Введение

Для построения математических моделей (ММ), предназначенных для решения практических задач (ПЗ) проектирования и эксплуатации различных систем [1-4], авторами был предложен в рамках современной неравновесной термодинамики, механики и электродинамики метод математического прототипирования энергетических процессов (ММПЭП) [5]. Построенные ММПЭП модели динамики процессов различной физической и химической природы не противоречат общим законам физики (законам сохранения, началам термодинамики, и т.д.), а также физическим особенностям конкретной системы [5, 6].

Для определения в соответствии с дифференциальными уравнениями (ДУ) динамики системы, полученных ММПЭП, контролируемых параметров (КП) системы по ее измеряемым параметрам

(ИП) необходимо идентифицировать большое число постоянных параметров, входящих в эти ДУ [5, 6], что представляет собой трудоемкую задачу [5-7]. Одним из путей борьбы с такой трудоемкостью является использование специальных методов интегрирования систем ДУ [8], позволяющих свести задачу к решению алгебраических уравнений [8] наряду с локальным упрощением ДУ [6, 8, 9]. Для снижения трудоемкости (в общем случае большой [8]) решения алгебраических уравнений могут быть применены методы интерполяции [8].

Настоящая работа посвящена построению ММ систем, применимых для решения упомянутых ПЗ [1-4], на базе методов ММПЭП [5, 6, 10] и методов интерполяции [8] и идентификации [7].

Материалы и методы

В соответствие с ММПЭП состояние системы характеризуется независимо от ее предыстории параметрами состояния (ПС), изменяющимися в результате протекания физических и химических процессов в ней [5, 6]. Причиной протекания процессов в системе являются динамические силы, определяемые через частные производные свободной энергии по ПС, расходуемой на протекание процессов в системе [5, 6]. Динамические силы и кинетические свойства однозначно определяют протекание процессов в системе [5, 6]. Динамика ПС в системе в свою очередь определяет динамику ИП и КП (рис. 1) [5, 6]. Приведенные на рис. 1 факторы определяют общую структуру ДУ динамики процессов, а также динамики ИП и КП, получаемых ММПЭП [5, 6].

Рис. 1. Факторы, определяющие динамику физических и химических процессов

В общем случае модель объекта (рис. 2) представляет собой алгоритм определения из ИП измеряемых характеристик объекта начального состояния и индивидуальных параметров системы с последующим определением КП системы (с учетом внешних условий протекания процессов) (видно из рис. 1) [5, 7]. Следует отметить, что в КП системы могут войти и ее ИП в последующие моменты времени, а также ее ИП, измерение которых в процессе эксплуатации осуществить невозможно. Для получения модели системы, пригодной для решения ПЗ (рис. 2), необходимо получить из уравнений ММПЭП аналитическое выражение динамики состояния системы, а также функции для определения ее начального состояния и ее индивидуальных параметров [10].

Рис. 2. Математическая модель системы для решения практических задач. В параллелограммах приведены обучаемые параметры модели

Аналитическое выражение общего решения системы уравнений ММПЭП задается в виде [6]

Ах() = Ах* (е(е0, Г), йи), я^, I()), и = й[и()], (1)

(2) (3)

где [6, 9]:

х() = х (Ах(), и( и), qTon, s (е(е0, i), ii( и), s(), qy

i(t) = s(е(ео,t), и(t, и), qупр), ео = x*-1 (xo,u(0, Iй), qTOn, qaHCC, q

е(ео,0) = ео, е(ео,t + т) = е(е(ео,т),t), dim(е) = dim(х),

Ve0, и, qдИСС, S 3 lim (Ах* ((ео, t), и, qДИСС, s)) = Ах** ( U, qДИСС, s), ( йх (t)

из

S 0 сдадует s (u,s,q^) s (u,s,q^);

V Zext

(4)

(5)

(6) (7)

из U (t )= U = const следует и (,и )= и (U),

где х(t) - динамика ПС системы; Ах(t) - динамика величин, характеризующих независимые ПС (не связанные законами сохранения с параметрами баланса [6]); U(t) - динамика характеристик системы,

не изменяющихся в результате протекания процессов внутри нее, а изменяющихся только в результате внешних воздействий на систему; (dx(t)/dt) t - внешние потоки. Как видно из (1)-(7),

функции X (Ах,и, qтоп,8 ) и Ах (и, qдисс,8) характеризуют топологическую и диссипативную составляющие системы соответственно [6], а функции § (е,и, 8, qупр) и £ (е,и, qупр) - ее управляющую

составляющую [6]. Для любого задания динамики ПС в виде (1)-(7) найдется система уравнений ММПЭП, общим решением которого будет эта динамика [6]. Отсюда (1)-(7) является способом задания динамики системы, альтернативным уравнениям ММПЭП. Квадратные скобки в уравнении (1) означают взятие функционала от и(*).

Измеряемые у () и контролируемые г (*) параметры системы [6]:

у(*) = У[*,х(*),и(*)], г(*) = Z[*,х(*),и(*)] . (8)

Для преобразования системы уравнений ММПЭП к виду, пригодному для решения ПЗ, используется локальное упрощение этих ДУ с последующим заданием соответствующих упрощенных аналитических выражений динамик х() в виде (1)-(7), из которых формируется полное решение [5, 6,

9, 11, 12].

Локальные упрощения системы уравнений ММПЭП осуществляются путем локальных упрощений функций состояния (ФС) для свойств веществ и процессов (СВП) с сохранением соответствующих ограничений на упрощенные ФС [5, 6]. Благодаря этому упрощенные уравнения ММПЭП также не противоречат законам физики [5, 6]. Запишем уравнения ММПЭП в виде [11]

^ = {(х (* ),и(' ),Р,а^ (х (*),и(* ),р,а)), (9)

где р - индивидуальные параметры системы; а - параметры ФС (задаваемых в виде функциональных разложений) для СВП [6, 10, 11]. При этом зависимость правой части (9) к w меньше, чем к х, и, р , а [9, 11]. Это дает возможность представить (9) в итерационном виде [11]:

dxk+1,r

— = f(xk+i,r(t),U(t),p,a,wKr(t)), wKr(t) = w(xk,r(t),U(t),p,a),

&

хКг е Хг (10)

или в эквивалентном виде [11]:

хк+1,г ( )= х" (х0,г, *, и г , q, ^^ к ,г ), ^^ к ,г = ^^ г [w к ,г ( )] , г = й г [и ( )] ,

w к г (* ) = w (,г (* ),и(* ),q), * е*, ^ ], х0,г = х (*г), (11)

где Xг - области локального упрощения (9). Правая часть (11) получается из (10) [11, 12]. Упрощения (10) могут быть такими, что система (9) при фиксированных wkr может распасться на независимые друг от друга подсистемы уравнений [11, 12]. При выборе достаточно малых областей Xг итерации (10) равномерно приближаются к искомой динамике х (*) [11] (один из специальных методов интегрирования систем ДУ [8]). Критерием остановки итераций (11) является заданная близость хкг () к х (*) из класса динамик (1)-(7), включающих в себя всю физику рассматриваемой системы [6, 10]. Это обеспечивает корректность модели системы, показанной на рис. 2 [10].

Синтез итераций построения преобразованной модели системы

Итерации (11) могут быть получены из уравнения (10) либо аналитически (в случае достаточно простой системы (10) [8]), либо с использованием специальных методов интегрирования систем ДУ [8].

В соответствии со специальными методами интегрирования систем ДУ динамика ПС х () рассматриваемой системы задается через опорные динамики хкг1 (е(е0,1),йгйг1.)), ■ = 1,поп ПС х(), полученные из уравнения (10) для соответствующих значений {(йг1, р1, а1, кг1)} . " в виде (1)-(7), в соответствии с работой [8]:

хк,г (х0,г ,*) = Гх,к,г ^{(хк,г,.■ (.г(*), Тг,.■ (*)), ск)} . 1 ,Iйг,Р,™к,Ьк^ , (12)

Тк,г. (*)= (е(е0, *), ск,г■,иг,Р,^к,г,Ьк,г ) '= 1Пп , (13)

и г ,.■ (* ) = и г ,.■ (1а г (, и г), Ск , II г , Р, ^ к ,г , Ь к , ) , . = П п , (14)

е0 = Гх,к ,г ( Х0,{Ск ,г} й г , Р, ^ к ,г , Ь к, ) , (15)

где постоянные коэффициенты Ск,г= (йг,р.,а.,йкг), . = 1,поП, Ьк,г = Ьк,г (а); функция гх,к,г ({(,ск,г,г )} ,иг,Р'^^к,г,Ьк,г) - непрерывна и ограничена; функции е*^ (е,ск,г, иг,Р,^к г,Ьк г),

. = 1, поп таковы, что функции

{хк ,г ,,■ ((,г ,,■ (е (е0, *) С к ,г,, )и г ■■ (и г (, и г ,,■),с к ,г,, ))}

Ск,г=(Ск,г,<,иг,Р,^к,г,Ьк,г ) ■ = 1поп

с учетом (13) представимы в виде (1)-(6), а функции (1 г (*,иг) удовлетворяют выражению (7). В таком случае нетрудно видеть, что итерации (11), полученные в силу (12)-(15), также удовлетворяют (1)-(7). Опорными динамиками хкг1 (е (е0, *), и г (*, и г1)), ■ = 1, поп могут быть общие аналитические решения

еще более локально упрощенной системы ДУ (10), качественное поведение которых такое же, как и качественное поведение общего решения системы ДУ (10). Отсюда упомянутые итерации (11)—(15) равномерно сходятся к общему решению некоторой системы ДУ, полученной ММПЭП, при любой непрерывной ограниченной функции Ех,к,г ({(хк,г■, Ск,гл)} 1, йг, Р, к г, Ьк,г), а значит, являются корректными [6].

Это дает возможность выбора структуры упомянутой интерполяционной функции (с точностью до исследуемых экспериментально входящих в уравнения (12)-(15) коэффициентов Ьк г = Ьк г (а)) [8].

Правые части (12)-(15) получаются из моделирования динамик на основе уравнений ММПЭП при различных значениях параметров иг, Р, ^кг [10, 13]. Правые части (13) и (14) при этом определяются таким образом, чтобы особые точки опорных динамик совпадали с особыми точками динамики хк г (х0 г, *), соответствующей иг, Р , IVк г (так как особые точки упомянутых динамик обусловливаются диссипативной и управляющей составляющими (видно из работы [6])). Это дает возможность задавать правые части (13) и (14) также в виде интерполяционных функций [8]. Что в свою очередь дает возможность одновременного синтеза правых частей (12)-(14) на динамиках, полученных ММПЭП, причем определяемые параметры входят в уравнения (12)-(14) линейно [8], что существенно упрощает синтез правых частей уравнений (12)-(14) [7].

Структура и опорные динамики правых частей уравнений (12)-(14) выбираются на основе качественного анализа динамик ПС х () системы, получаемых из уравнений ММПЭП. В качестве интерполяционных функций, через которые задаются правые части уравнений (12)—(14), может быть взята интерполяция сплайнами [8], интерполяционный полином Лагранжа [8], нейронные сети как универсальный интерполятор [14, 15], а также используя вейвлеты [16].

Идентификация преобразованной модели

В общем случае динамики ИП у () (8) являются динамиками медленно меняющихся составляющих динамик у () некоторых характеристик рассматриваемой системы [17]:

у (* ) = У (х (*), и (*)), у () = У [*, у ()] , (16)

причем динамика у () полностью восстанавливается по динамике у () [17]:

у () = У[*,у()] . (17)

Приведенный подход обусловливается ограниченной пропускной способностью измерительных датчиков [17]. Из уравнения (16) видно, что, выбрав у в качестве части параметров состояния х, мы, доведя итерации (11)—(15) до полной модели

х ()=Гх ({(х, (((), й (*)), с, )}\Ь ^, £0=г-1 ({с К!;, Ь), (18)

^(*) = к ((е0,*),с,.,Ь), й-() = й, (й(и,),с,,Ь), г = 1,поп, (19)

Ь = Ь(Ь,й,р), Ь = Ь(а), (20)

представимой в виде (1)-(7), можем вышеописанным путем обучение модели (18)-(20) свести к обучению (по данным у ()) правых частей (18), (19).

Далее на основе полученных коэффициентов Ь мы определяем, используя (20), значения параметров Ь, й, р , на основе которых в силу (1)-(8) определяются КП г(*). Пусть имеется для некоторой совокупности экземпляров рассматриваемой системы некоторого бренда множество данных ИП

Г ^ 1 пэкз _

{{(й,,],у,,] (*))} } , для каждой совокупности у(), ] = 1,, г = 1, пэкз из них описанным выше

I- ]=^ ,=1

путем определим соответствующие 1Ь {] , е0 {], ] = 1, Ni, г = 1, пэкз. Здесь пэкз - число экземпляров системы рассматриваемого бренда, а Ni, г = 1, пэкз - число испытаний для каждого экземпляра. Затем значения параметров Ь, р,, , = 1, пэкз с использованием (20) определим из работы [10]:

Ь,]. = Ь(ь,й,,],р,.), ]=, г = 1,Пэкз, (21)

получив при этом из уравнения (21) функционалы для Ь, р,, г = 1, пэз [10]: Ь = Ь

Û j 'b

, Ь,,j = b [y,, (t)] , j =1,N , i = 1, nx3, (22)

p, = p (b,{ûi.j, bi.j }j= ), b,,j = F[ y i, j (t )] , j = Щ, i = 1, nx3, (23)

где функционал Ь представляет собой вышеописанный алгоритм определения параметров Ь, входящих в уравнения (18) и (19), из динамики у(). Функции (22) и (23) могут быть получены путем аналитического либо численно-аналитического (с использованием методов [8, 14-16, 18]) преобразования (21) в уравнения (22) и (23) [10]. Соответственно, правые части уравнений (12)-(14) желательно так выбирать, чтобы уравнения (22) и (23) были получены аналитически из формулы (21).

На основе выражения (23) получим функцию для индивидуальных параметров системы, определяемых через контрольные динамики ИП системы у(кошр) (), г = 1, N при соответствующих кон-

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

трольных воздействиях й(контр), , = 1, N на нее, входящую в архитектуру модели (см. рис. 2) [10]:

Р = Р(ь, {й(контр), Ь(контр)}^), Ь(контр) = Ь[у(контр) (*)], ■ = 1, ^конр . (24)

Аналогично (24) вводится функционал Ье, для определения е0 через текущую динамику ИП у ():

е0 = К [у (*)] , (25)

представляющий собой функцию для начального состояния системы, входящую в архитектуру модели (см. рис. 2) [10].

В частном случае в силу системы (18)-(20) динамика ПС х() может быть выражена через контрольные динамики с точностью до параметров Ь . Это облегчает определение параметров Ь из статистики ИП, что, например, и было сделано в работе [19].

Результаты

Выражения (17)-(20), (24), (25) представляют собой ММ системы (см. рис. 2) [10], пригодную для решения ПЗ [1-4, 10]. Выражения (17) и (22) представляют собой алгоритм обучения упомянутой ММ системы [10].

ММ системы (17)-(20), (24), (25), структура которой показана на рис. 2, строится в соответствии со следующим алгоритмом:

1) синтезируем систему ДУ динамики физических и химических процессов в рассматриваемой системе на базе ММПЭП;

2) на основе особенностей датчиков ИП рассматриваемой системы формируем соотношение (17);

3) выполняем локальные упрощения ФС для СВП рассматриваемой системы, получив тем самым упрощенные ДУ динамики системы;

4) для локально упрощенной системы ДУ динамики системы рассчитываем различные динамики ее ПС, ИП и КП при различных случайно сгенерированных значениях экспериментально исследуемых параметров ДУ;

5) на основе полученных динамик ПС, ИП и КП рассматриваемой системы формируем итерации (11)—(15) с использованием методов интерполяции;

6) на основе сформированных итераций (11)—(15) получаем выражения (18)—(20);

7) методами теории идентификации, исходя из преобразованной ММ (18)—(20) и из особенностей испытаний рассматриваемой системы, формируем выражения (22), (24), (25);

8) на основе (22) обучаем полученную ММ (17)—(20), (24), (25);

9) тестируем полученную ММ (17)—(20), (24), (25) и при необходимости возвращаемся к предыдущим пунктам предложенного алгоритма с целью корректировки синтезированной ММ (17)—(20), (24), (25).

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

Обсуждение

Предложенная в настоящей работе методика синтеза ММ (17)—(20), (24), (25) является развитием методов, изложенных в работах [7—10, 13, 14, 16, 20]. Основным отличием предложенного метода от методов, изложенных в работах [14, 16], является использование базисов, вбирающих в себя физику конкретной рассматриваемой системы, что гарантирует корректность построенной модели [5, 10, 20]. Использование локального упрощения уравнений ММПЭП с последующим синтезом сложных динамик ПС из простых динамик ПС, полученных из упрощенных уравнений ММПЭП, в предложенном методе позволяет существенно сократить по сравнению с методами [13] объем вычислений [20].

Заключение

Методика построения ММ (17)—(20), (24), (25) произвольной системы дает ММ системы, которые могут обучаться в том числе на лету (в процессе функционирования объекта). Благодаря этому ММ (17)—(20), (24), (25) целесообразно использовать в составе цифровых двойников.

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

1. Юревич Е. И. Основы проектирования техники. СПб. : Санкт-Петербургский государственный политехнический университет, 2012. 135 с.

2. Барзилович Е. Ю. Модели технического обслуживания сложных систем. М. : Высш. шк., 1982. 231 с.

3. Колодежный Л. П., Чернодаров А. В. Надежность и техническая диагностика. М. : Изд-во ВВА им. проф. Н. Е. Жуковского и Ю. А. Гагарина, 2010. 452 с.

4. Бессекерский В. А., Попов Е. П. Теория систем автоматического управления. СПб. : Профессия, 2003. 768 с.

5. Khalyutin S. P., Starostin I. E., Agafonkina I. V. Generalized Method of Mathematical Prototyping of Energy Processes for Digital Twins Development // Energies. 2023. Vol. 16. P. 1933-1958. doi: 10.3390/en16041933

6. Старостин И. Е., Дружинин А. А. Аналитическое приближение решений уравнений метода математического прототипирования энергетических процессов путем качественного анализа этих уравнений // Надежность и качество сложных систем. 2023. № 2 (42). С. 22-31. doi: 10.21685/2307-4205-2023-2-3

7. Eykhoff P. Systems identification: parametrs and state estimation. Eindhoven, Netherlands : University of technology, 1975. 680 p.

8. Калиткин Н. Н. Численные методы. СПб. : БХВ-Петербург, 2011. 592 с.

9. Ланцов В. Н. Методы понижения порядка моделей сложных систем. Владимир : Изд-во ВлГУ, 2017. 84 с.

10. Старостин И. Е., Дружинин А. А., Гавриленков С. И. Использование машинного обучения с учителем для построения математических моделей систем методом математического прототипирования энергетических процессов // Труды Международного симпозиума Надежность и качество. 2023. Т. 1. С. 66-72.

11. Старостин И. Е. Программная реализация решения потенциально-потоковым методом задач построения моделей систем из результатов испытаний этих систем // Надежность и качество сложных систем. 2020. № 3 (33). С. 128-136. doi: 10.21685/2307-4205-2020-3-15

12. Starostin I. E., Khalyutin S. P., Bykov V. I. Simplification of potential-flow equations of physical and chemical processes in dynamics for obtaining a mathematical model of a system // The Complex Systems. 2019. № 1. P. 74-87.

13. Starostin I. E., Khalyutin S. P. Identification of system models from potential-stream equations on the basis of deep learning on experimental data // Civil Aviation High Technologies. 2020. Vol. 23, № 2. P. 47-58.

14. Круглов В. В., Дли М. И., Голунов Р. Ю. Нечеткая логика и искусственные нейронные сети. М. : Физмат-лит, 2001. 221 с.

15. Галкин В. А., Гавриленко Т. В., Смородинов А. Д. Некоторые аспекты аппроксимации и интерполяции функций искусственными нейронными сетями // Вестник КРАУНЦ. Физ.-мат. науки. 2022. Т. 38, № 1. С. 54-73. doi: 10.26117/2079-6641-2022-38-1-54-73

16. Дьяконов В. П. Вейвлеты. От теории к практике. М. : СОЛОН-Пресс, 2010. 440 с.

17. Шишкин И. Ф. Теоретическая метрология. Ч. 2. Обеспечение единства измерений. СПб. : Питер, 2012. 240 с.

18. Flach P. Machine Learning. The Art and Science of Algorithms that Make Sense of Data. Cambridge : Cambridge University Press, 2015. 400 p.

19. Starostin I. E., Khalyutin S. P., Davidov A. O. [et al.]. Obtaining a model for the voltage and temperature of the US18650VTC6 Series lithium-ion battery in constant current discharge mode from the analysis of physical and chemical processes in the accumulator // 18th Technical Scientific Conference on aviation dedicated to the memory of N.E. Zhukovsky, TSCZH 2021. 2021. P. 109-117.

20. Малинецкий Г. Г. Математические основы синергетики: Хаос, структуры, вычислительный эксперимент. М. : Ленанд, 2017. 312 с.

References

1. Yurevich E.I. Osnovy proektirovaniya tekhniki = Fundamentals of engineering design. Saint Petersburg: Sankt-Peterburgskiy gosudarstvennyy politekhnicheskiy universitet, 2012:135. (In Russ.)

2. Barzilovich E.Yu. Modeli tekhnicheskogo obsluzhivaniya slozhnykh system = Models of maintenance of complex systems. Moscow: Vyssh. shk., 1982:231. (In Russ.)

3. Kolodezhnyy L.P., Chernodarov A.V. Nadezhnost' i tekhnicheskaya diagnostika = Reliability and technical diagnostics. Moscow: Izd-vo VVA im. prof. N.E. Zhukovskogo i Yu.A. Gagarina, 2010:452. (In Russ.)

4. Bessekerskiy V.A., Popov E.P. Teoriya sistem avtomaticheskogo upravleniya = Theory of automatic control systems. Saint Petersburg: Professiya, 2003:768. (In Russ.)

5. Khalyutin S.P., Starostin I.E., Agafonkina I.V. Generalized Method of Mathematical Prototyping of Energy Processes for Digital Twins Development. Energies. 2023;16:1933-1958. doi: 10.3390/en16041933

6. Starostin I.E., Druzhinin A.A. Analytical approximation of solutions of equations of the method of mathematical prototyping of energy processes by qualitative analysis of these equations. Nadezhnost' i kachestvo slozhnykh system = Reliability and quality of complex systems. 2023;(2):22-31. (In Russ.). doi: 10.21685/2307-4205-2023-2-3

7. Eykhoff P. Systems identification: parametrs and state estimation. Eindhoven, Netherlands: University of technology, 1975:680.

8. Kalitkin N.N. Chislennye metody = Numerical methods. Saint Petersburg: BKhV-Peterburg, 2011:592. (In Russ.)

9. Lantsov V.N. Metody ponizheniya poryadka modeley slozhnykh system = Methods of lowering the order of models of complex systems. Vladimir: Izd-vo VlGU, 2017:84. (In Russ.)

10. Starostin I.E., Druzhinin A.A., Gavrilenkov S.I. Using machine learning with a teacher to build mathematical models of systems by mathematical prototyping of energy processes. Trudy Mezhdunarodnogo simpoziuma Nadezhnost' i kachestvo = Proceedings of the International Symposium Reliability and Quality. 2023;1:66-72. (In Russ.)

11. Starostin I.E. Software implementation of the solution of problems of building models of systems from the test results of these systems by a potentially streaming method. Nadezhnost' i kachestvo slozhnykh system = Reliability and quality of complex systems. 2020;(3):128-136. (In Russ.). doi: 10.21685/2307-4205-2020-3-15

12. Starostin I.E., Khalyutin S.P., Bykov V.I. Simplification of potential-flow equations of physical and chemical processes in dynamics for obtaining a mathematical model of a system. The Complex Systems. 2019;(1):74-87.

13. Starostin I.E., Khalyutin S.P. Identification of system models from potential-stream equations on the basis of deep learning on experimental data. Civil Aviation High Technologies. 2020;23(2):47-58.

14. Kruglov V.V., Dli M.I., Golunov R.Yu. Nechetkaya logika i iskusstvennye neyronnye seti = Fuzzy logic and artificial neural networks. Moscow: Fizmatlit, 2001:221. (In Russ.)

15. Galkin V.A., Gavrilenko T.V., Smorodinov A.D. Some aspects of approximation and interpolation of functions by artificial neural networks. Vestnik KRAUNTs. Fiz.-mat. nauki = Vestnik KRAUNTS. Phys.-mat. nauki. 2022; 38(1):54-73. (In Russ.). doi: 10.26117/2079-6641-2022-38-1-54-73

16. D'yakonov V.P. Veyvlety. Ot teorii k praktike = Wavelets. From theory to practice. Moscow: SOLON-Press, 2010:440. (In Russ.)

17. Shishkin I.F. Teoreticheskaya metrologiya. Ch. 2. Obespechenie edinstva izmereniy = Theoretical metrology. Part. 2. Ensuring the uniformity of measurements. Saint Petersburg: Piter, 2012:240. (In Russ.)

18. Flach P. Machine Learning. The Art and Science of Algorithms that Make Sense of Data. Cambridge: Cambridge University Press, 2015:400.

19. Starostin I.E., Khalyutin S.P., Davidov A.O. et al. Obtaining a model for the voltage and temperature of the US18650VTC6 Series lithium-ion battery in constant current discharge mode from the analysis of physical and chemical processes in the accumulator. 18th Technical Scientific Conference on aviation dedicated to the memory of N.E. Zhukovsky, TSCZH 2021. 2021:109-117.

20. Malinetskiy G.G. Matematicheskie osnovy sinergetiki: Khaos, struktury, vychislitel'nyy eksperiment = Mathematical foundations of synergetics: Chaos, structures, computational experiment. Moscow: Lenand, 2017:312. (In Russ.)

Информация об авторах / Information about the authors

Игорь Евгеньевич Старостин

кандидат технических наук,

доцент кафедры электротехники

и авиационного электрооборудования,

Московский государственный технический

университет гражданской авиации

(Россия, г. Москва, Кронштадтский бульвар, 20)

E-mail: starostinigo@yandex.ru

Igor E. Starostin

Candidate of technical sciences, associate professor of the sub-department of electrical engineering and aviation electrical equipment, Moscow State Technical University of Civil Aviation (20 Kronshtadtskiy boulevard, Moscow, Russia)

Автор заявляет об отсутствии конфликта интересов / The author declares no conflicts of interests.

Поступила в редакцию/Received 12.02.2024 Поступила после рецензирования/Revised 22.02.2024 Принята к публикации/Accepted 29.02.2024

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