Journal influence
Bookmark
Next issue
An algorithm and software implementation of test object model synthesis based on the solution of the nonparametric identification equation
Abstract:The paper considers the development of the theory of testing in general and the experimental-theoretical method in particular. In the aspect of this issue, the authors have developed an algorithm for synthesizing a model of a test object based on solving the equation of nonparametric identification of a dynamic system using hyperdelta approximation and the Laplace transform. Unlike the existing ones, the algorithm is applicable to input and output signals of arbitrary shape and physical quantities. In addition, it does not require large computing resources. Taking into account these features, the algorithm enables formalizing a multidimensional relationship between factors and performance characteristics of the test object through repeated use for different input and output signals. The authors have implemented a mathematical library for identifying a test object model and an application with a graphical user interface for automating calculations using the C++ and Python programming languages. The presented software solution is made similar to classical machine learning models. To substantiate the possibility of using the developed algorithm, the authors carried out a computational experiment that involved various types of input and output signals (periodic, non-periodic and random) with different hyperdelta approximation accuracy. Based on the results of the computational experiment, the authors have made recommendations on using the algorithm. In particular, they recommended to increase the number of initial moments of the hyperdelta approximation at high amplitudes of the output signal.
Аннотация:Настоящая работа посвящена развитию теории испытаний в целом и опытно-теоретического метода в частности. Авторами разработан алгоритм синтеза модели объекта испытаний, основанный на решении уравнения непараметрической идентификации динамической системы с использованием гипердельтной аппроксимации и преобразования Лапласа. В отличие от существующих данный алгоритм применим для входных и выходных сигналов произвольной формы и физических величин. Кроме того, он не требует больших вычислительных ресурсов. Алгоритм позволяет формализовать многомерную зависимость между факторами и тактико-техническими характеристиками объекта испытаний. С помощью языков программирования C++ и Python реализованы математическая библиотека идентификации модели объекта испытаний и приложение с графическим пользовательским интерфейсом для автоматизации расчетов. Представленное программное решение выполнено по аналогии с классическими моделями машинного обучения. Для обоснования возможности применения разработанного алгоритма проведен вычислительный эксперимент на различных типах входных и выходных сигналов (периодических, непериодических и случайных) с разной точностью гипердельтной аппроксимации. По результатам вычислительного эксперимента получены рекомендации по использованию алгоритма, в частности, при высоких амплитудах выходного сигнала следует увеличить количество начальных моментов гипердельтной апроксимации.
Authors: Gusenitsa, Ya.N. (yaromir226@gmail.com ) - Military Innovative Technopolis «ERA» (Head of Testing Laboratory), Anapa, Russia, Ph.D, Mingachev, E.R. (era_otd1@mil.ru) - Military Innovative Technopolis "ERA" (Senior Operator), Anapa, Russia, Iskhakov, N.U. (era_otd1@mil.ru) - Military Innovative Technopolis "ERA" (Senior Operator), Anapa, Russia, Kolokolov, M.I. (era_otd1@mil.ru) - Military Innovative Technopolis "ERA" (Operator), Anapa, Russia | |
Keywords: stochastic processes, dynamic system modeling, non-parametric identification, mathematical model, test object |
|
Page views: 1792 |
PDF version article |
Для современных испытаний характерны значительная продолжительность, сложность, трудоемкость, высокая стоимость, наличие ограничений по воспроизведению всего комплекса условий функционирования испытываемого объекта [1]. Например, стоимость проведения испытаний ракетно-космической техники составляет от 50 до 80 % общих затрат на ее разработку, ракетно-артиллерийского вооружения – от 45 до 60 % [2]. Поэтому в настоящее время наблюдается тенденция к сокращению объема испытаний. Одним из подходов, хорошо зарекомендовавших себя на практике при проведении ограниченного объема испытаний, является опытно-теоретический метод. Он предполагает проведение расчетов на математических моделях объекта испытаний в сочетании с натурными экспериментами [3, 4]. При этом результаты натурных экспериментов используются в качестве исходных данных для моделирования, а также для проверки правильности объекта испытаний [5, 6]. Основным достоинством метода является то, что он позволяет обеспечить высокую полноту результатов испытаний [7]. Однако, поскольку достоверность результатов испытаний зависит от адекватности используемых моделей, для его применения необходима идентификация параметров этих моделей. Целью настоящей работы являются разработка и программная реализация алгоритма синтеза модели объекта испытаний на основе решения уравнения непараметрической идентификации. Структура модели объектаиспытаний При реализации опытно-теоретического метода важное место отводится разработке математической модели объекта испытаний, которая в общем виде представляет собой вектор тактико-технических характеристик: , где yi(t) – i-я тактико-техническая характеристика; K – общее количество тактико-технических характеристик. Каждая i-я тактико-техническая характеристика yi(t) зависит от различных параметров: yi(t) = f(xj, t), где xj(t) – j-й параметр, определяющий значение i-й тактико-технической характеристики объекта испытаний. Данная зависимость может быть определена на основе интегрального уравнения Дюамеля: , (1) где h(t) – неизвестная весовая функция, позволяющая определять y(t) при произвольном значении x(t). Учитывая, что значение i-й тактико-технической характеристики зависит от нескольких параметров, выражение (1) можно представить следующим образом: , (2) где L – общее количество параметров. Модель объекта испытаний представляет собой матрицу H[KL](t), элементами которой являются весовые функции hij(t) (табл. 1). Данная матрица является исчерпывающим описанием объекта испытаний (рис. 1) и позволяет при произвольном значении j-го параметра xj(t) определять i-ю тактико-техническую характеристику yi(t). Важное место в синтезе модели объекта испытаний отводится решению уравнения непараметрической идентификации, которое при наличии по одному входному и выходному случайным сигналам в общем виде будет , (3) где Kyx(t) – взаимно-корреляционная функция входного x(t) и выходного y(t) сигналов; Kxx(t) – автокорреляционная функция входного x(t) сигнала. Для решения уравнения (3) можно использовать различные временные и частотные (спектральные) методы и алгоритмы, как, например, в работах [8–10]. Вместе с тем существующие методы непараметрической идентификации имеют ограниченные условия применения. Так, например, одни методы оказываются непригодными для идентификации динамических систем, входные и выходные сигналы которых имеют вероятностный характер [11, 12], другие требуют исчерпывающей информации о входных и выходных сигналах и существенных затрат вычислительных ресурсов [13, 14]. Поэтому предлагается использовать алгоритм синтеза модели объекта испытаний, основанный на методе решения уравнения непараметрической идентификации динамической системы, который лишен указанных недостатков [15]. Содержание алгоритма Для нахождения элементов матрицы H[KL](t) (табл. 1) необходимо решить уравнение непараметрической идентификации относительно соответствующих входного и выходного сигналов. Решение проходит в несколько этапов. 1. Вычисление автокорреляционной функции Kxx(t) входного x(t) сигнала и взаимно-кор-реляционной функции Kyx(t) входного x(t) и выходного y(t) сигналов. 2. Нормализация функций Kxx(t) и Kyx(t), то есть их перенос в положительную ось и приведение области значений от 0 до 1. В результате получаются нормализованные функции Fxx(t) и Fyx(t). 3. Вычисление гипердельтной аппроксимации функций Fxx(t) и Fyx(t) по формуле , где Сi – вероятности и Ti – параметры, удовлетворяющие приведенной ниже системе уравнений; δ(∙) – дельта-функция Дирака. Для получения вероятностей Сi и параметров Ti на основе нормализованных функций Fxx(t) и Fyx(t) рассчитываются начальные моменты случайной величины νi. Далее решается система уравнений относительно переменных Сi и Ti: (4) где νi – i-й начальный момент случайной величины, распределенной с плотностью вероятности f(t); N – степень гипердельтной аппроксимации (точность приближения). Для степеней гипердельтной аппроксимации N > 2 система решается с использованием численных методов. В данной работе применяется метод Ньютона–Рафсона. В результате получаются аппроксимированные плотности распределения вероятностей fΔxx(t) и fΔyx(t). 4. Нахождение функций распределения на основе вычисленных плотностей распределения вероятностей: , . 5. Нахождение изображений по Лапласу полученных функций: , и весовой функции: . 6. Нахождение весовой функции с использованием приближенного способа обращения преобразования Лапласа: h(t) = sh*(s) при s = 1/t. В результате получается весовая функция h(t) для входного x(t) и выходного y(t) сигналов (тактико-технических характеристик). Данная процедура повторяется для всех пар входных и выходных сигналов – элементов матрицы весовых функций H[KL](t). Таким образом, происходит идентификация модели объекта испытаний. Полученная матрица (модель) позволяет при произвольном значении l-го параметра xl(t) предсказывать k-ю тактико-техническую характеристику yk(t). Для этого применяются формулы интегрального уравнения Дюамеля (3) и его модификации для нескольких входных параметров (4). Программная реализация алгоритма На основе предложенного алгоритма с использованием языков программирования C++ (стандарта C++17) и Python (версии 3.11) разработана математическая библиотека идентификации модели объекта испытаний. Как и классические модели машинного обучения, модель объекта испытаний является C++/Python объектом с методами fit для обучения модели на выборке и predict для вычисления предсказания значений параметров [16]. Для демонстрации работы модели и автоматизации процесса идентификации модели объекта испытаний разработано приложение с графическим пользовательским интерфейсом. Приложение написано на языке Python (версии 3.11) с использованием библиотеки PyQt5 (5.15.7) для реализации графического интерфейса. Приложение состоит из двух частей – пользовательского интерфейса (frontend), реализующего взаимодействие с пользователем, и расчетного блока (backend), обеспечивающего автоматизацию идентификации модели объекта испытаний. Расчетный блок включает модули Hyperdelta для решения уравнения непараметрической идентификации по алгоритму и NewtonRaphson для численного решения системы уравнений (4) по методу Ньютона–Рафсона. В пользовательский интерфейс входят модули, реализующие различные графические элементы и их логику взаимодействия с пользователем, а также обработку входных данных. Пользовательский интерфейс представлен на рисунке 2. Приложение позволяет выполнять следующие функции. 1. Задавать наборы входных xi(t) (параметров) и выходных yj(t) сигналов (тактико-техни-ческих характеристик) и их взаимосвязи (рис. 2a–в). 2. Идентифицировать модель объекта испытаний в виде матрицы H[KL](t) весовых функций hij(t) с необходимой точностью аппроксимации (рис. 2в). 3. Вычислять предсказания относительно характеристик объекта испытаний – значения выходных y'j(t) сигналов для произвольных входных x'i(t) (рис. 2г). Результаты Для обоснования возможности применения разработанного алгоритма проведен вычислительный эксперимент на различных типах входных и выходных сигналов (периодических, непериодических и случайных) с разной точностью гипердельтной аппроксимации. Примеры сигналов и результатов идентифика-ции представлены на рисунке 3, где синим цветом выделен выходной сигнал, желтым – идентифицированный выходной сигнал на основе использования гипердельтной аппроксимации с двумя начальными моментами, красным – идентифицированный выходной сигнал с восьмью начальными моментами. Данный пример иллюстрирует следующие тенденции. Достаточно точно аппроксимируются простые непериодические (a) и периодические (б) сигналы. Случайные сигналы без четкого распределения (в) упрощаются до средних значений. При этом с ростом количества начальных моментов гипердельтной аппроксимации увеличивается амплитуда результирующих значений. Для получения оценки быстродействия разработанного алгоритма нормировалось общее время расчетов и количество итераций. Расчеты производились по 10 раз для различного количества начальных моментов гипердельтной аппроксимации. Результаты (средние, минимальные и максимальные значения измеряемых показателей) представлены в таблице 2. На основе данных таблицы можно сделать следующий вывод: среднее время идентификации растет экспоненциально с увеличением количества начальных моментов гипердельтной аппроксимации, как растет и разброс минимального и максимального времени расчета. Аналогичную тенденцию можно заметить и в отношении количества итераций метода Ньютона–Рафсона при решении системы уравнений (4). Как показывает эксперимент, основную часть времени расчетов занимает итеративное решение системы уравнений (4). Время расчета остальных элементов уравнения непараметрической идентификации в целом не изменяется и по сути зависит только от объема входных данных. Заключение В настоящей работе представлен алгоритм синтеза модели объекта испытаний. Алгоритм основан на методе решения уравнения непараметрической идентификации динамической системы, в отличие от существующих не требующем исчерпывающей информации об объекте испытаний и существенных затрат вычислительных ресурсов. Математическая библиотека идентификации модели объекта испытаний и приложение с графическим пользовательским интерфейсом для автоматизации расчетов реализованы с помощью языков C++ и Python. Полученные результаты позволяют обеспечить многофакторную оценку динамических параметров объектов испытаний и могут быть применены при проверке соответствия тактико-технических характеристик предъявляемым требованиям. Список литературы
Reference List
|
Permanent link: http://swsys.ru/index.php?page=article&id=5007&lang=en |
Print version |
The article was published in issue no. № 2, 2023 [ pp. 320-326 ] |
Perhaps, you might be interested in the following articles of similar topics:
- Использование алгоритмов оптимизации с самообучением для управления динамически изменяющимися системами
- Использование программного обеспечения для определения и прогнозирования показателей качества экструдированной продукции
- Модели планирования производства изделий, основанных на нанотехнологиях
- Математические модели реограмм состояния в программах Table Curve 2d/3d как основа интеллектуальной системы управления процессами структурирования многокомпонентных эластомерных композитов
- Информационно-расчетные системы для определения лазерно-локационных характеристик объектов
Back to the list of articles