На правах рекламы:
ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Авторитетность издания

ВАК - К1
RSCI, ядро РИНЦ

Добавить в закладки

Следующий номер на сайте

2
Ожидается:
16 Июня 2024

Алгоритм и программная реализация синтеза модели объекта испытаний на основе решения уравнения непараметрической идентификации

An algorithm and software implementation of test object model synthesis based on the solution of the nonparametric identification equation
Дата подачи статьи: 20.12.2022
Дата после доработки: 17.02.2023
Дата принятия к публикации: 21.02.2023
УДК: 004.42
Статья опубликована в выпуске журнала № 2 за 2023 год. [ на стр. 320-326 ]
Аннотация:Настоящая работа посвящена развитию теории испытаний в целом и опытно-теоретического метода в частности. Авторами разработан алгоритм синтеза модели объекта испытаний, основанный на решении уравнения непараметрической идентификации динамической системы с использованием гипердельтной аппроксимации и преобразования Лапласа. В отличие от существующих данный алгоритм применим для входных и выходных сигналов произвольной формы и физических величин. Кроме того, он не требует больших вычислительных ресурсов. Алгоритм позволяет формализовать многомерную зависимость между факторами и тактико-техническими характеристиками объекта испытаний. С помощью языков программирования C++ и Python реализованы математическая библиотека идентификации модели объекта испытаний и приложение с графическим пользовательским интерфейсом для автоматизации расчетов. Представленное программное решение выполнено по аналогии с классическими моделями машинного обучения. Для обоснования возможности применения разработанного алгоритма проведен вычислительный эксперимент на различных типах входных и выходных сигналов (периодических, непериодических и случайных) с разной точностью гипердельтной аппроксимации. По результатам вычислительного эксперимента получены рекомендации по использованию алгоритма, в частности, при высоких амплитудах выходного сигнала следует увеличить количество начальных моментов гипердельтной апроксимации.
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.
Авторы: Гусеница Я.Н. (yaromir226@gmail.com ) - Военный инновационный технополис «ЭРА» (начальник испытательной лаборатории информатики и вычислительной техники), Анапа, Россия, кандидат технических наук, Мингачев Э.Р. (era_otd1@mil.ru) - Военный инновационный технополис «ЭРА» (старший оператор научной роты), Анапа, Россия, Исхаков Н.У. (era_otd1@mil.ru) - Военный инновационный технополис «ЭРА» (старший оператор научной роты), Анапа, Россия, Колоколов М.И. (era_otd1@mil.ru) - Военный инновационный технополис «ЭРА» (оператор научной роты), Анапа, Россия
Ключевые слова: случайные процессы, динамическая система, непараметрическая идентификация, математическая модель, объект испытаний
Keywords: stochastic processes, dynamic system modeling, non-parametric identification, mathematical model, test object
Количество просмотров: 1002
Статья в формате PDF

Размер шрифта:       Шрифт:

Для современных испытаний характерны значительная продолжительность, сложность, трудоемкость, высокая стоимость, наличие ограничений по воспроизведению всего комплекса условий функционирования испытываемого объекта [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) сигнала.

Таблица 1 
Модель объекта испытаний
Table 1 
The test object model

Тактико-техническая 
характеристика 
объекта испытаний	Параметр
	x1(t)	x2(t)	…	xj(t)	…	xL(t)
y1(t)	h11(t)	h12(t)	…	h1l(t)	…	h1L(t)
y2(t)	h21(t)	h22(t)	…	h2l(t)	…	h2L(t)
…	…	…	…	…	…	…
yi(t)	hk1(t)	hk2(t)	…	hkl(t)	…	hkL(t)
…	…	…	…	…	…	…
yK(t)	hK1(t)	hK2(t)	…	hKl(t)	…	hKL(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), обеспечивающего автоматизацию идентификации модели объекта испытаний.

Расчетный блок включает модули Hyper­delta для решения уравнения непараметрической идентификации по алгоритму и Newton­Raphson для численного решения системы уравнений (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.

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

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

  1. Гусеница Я.Н., Хайбуллов О.Д., Загрутдинов Ю.А., Скорик Е.А. Опытно-теоретический метод исследования степени защиты специальных сооружений // Изв. ТулГУ. Технич. науки. 2022. № 12. С. 264–267.
  2. Любарчук Ф.Н. Методическое обеспечение научно-методических положений сокращенных испытаний образцов ракетно-артиллерийского вооружения // Изв. РАРАН. 2016. № 1. С. 80–86.
  3. Арсеньев В.Н., Дубинин Д.П. Обоснование метода оценивания характеристик сложных систем при ограниченном числе натурных испытаний // ОКНТПР. 2020. № 2. С. 30–36.
  4. Козлов Н.Н., Красный В.П., Решетников А.В. Особенности современной методологии испытаний систем вооружения воздушно-космической обороны // Военная мысль. 2015. № 6. С. 42–50.
  5. Александровская Л.Н., Кириллин А.В., Кербер О.Б., Иосифов П.А. Один подход к сокращению объема натурных испытаний при отработке сложных технических систем // Тр. ФГУП «НПЦАП». Системы и приборы управления. 2018. № 3. С. 59–63.
  6. Кофман В.М. Математическая модель обратной термогазодинамической задачи для экспериментально-расчетной оценки показателей эффективности камеры сгорания и турбины опытного газогенератора ГТД // Вестн. ПНИПУ. Аэрокосмическая техника. 2018. № 53. С. 97–108. doi: 10.15593/2224-9982/2018.53.09.
  7. Найденов В.Г., Крупский К.А., Бочкарев А.В. Методический подход к оценке потребного количества натурных экспериментов при проведении испытаний сложных образцов вооружения, военной и специальной техники // Вооружение и экономика. 2015. № 1. С. 4–11.
  8. Антонова Т.В. Методы идентификации параметра в ядре уравнения первого рода типа свертки на классе функций с разрывами // СибЖВМ. 2015. Т. 18. № 2. C. 107–120. doi: 10.15372/SJNM20150201.
  9. Воскобойников Ю.Е., Боева В.А. Алгоритмы непараметрической идентификации сложных технических систем // Науч. вестн. НГТУ. 2020. № 4. С. 47–64. doi: 10.17212/1814-1196-2020-4-47-64.
  10. Данилов А.М., Гарькина И.А. Ретроспективная идентификация динамической системы // Приднепровский науч. вестн. 2018. Т. 6. № 4. С. 62–65.
  11. Шишкина А.В. Об анализе непараметрических алгоритмов идентификации // Решетневские чтения: сб. Междунар. конф. 2018. Т. 2. № 22. С. 180–181.
  12. Ярещенко Д.И. О непараметрической идентификации частично-параметризованного дискретно-непрерывного процесса // СибЖНТ. 2020. Т. 21. № 1. С. 47–53. doi: 10.31772/2587-6066-2020-21-1-47-53.
  13. Корнеева А.А., Чернова С.С., Шишкина А.В. Непараметрические алгоритмы восстановления взаимно неоднозначных функций по наблюдениям // СибЖНТ. 2017. Т. 18. № 3. С. 510–519.
  14. Вожов С.С., Чимитова Е.В. Параметрическая и непараметрическая идентификации закона распределения по интервальным данным // Метрология. 2018. № 1. С. 3–16.
  15. Гусеница Я.Н. Решение уравнения непараметрической идентификации динамической системы на основе гипердельтной аппроксимации // Проблемы машиностроения и надежности машин. 2022. Т. 21. № 1. С. 96–102. doi: 10.31857/s0235711921060092.
  16. Wang H. Machine learning basics. In: Goodfellow I., Bengio Y., Courville A. (eds.) Deep Learning, 2016, pp. 98–164.

 

Reference List

  1. Gusenitsa, Ya.N., Khaybullov, O.D., Zagrutdinov, Yu.A., Skorik, E.A. (2022) ‘Experimental-theoretical method for studying the degree of protection of special structures’, Proc. of the TSU. Tech. Sci., (12), pp. 264–267 (in Russ.).
  2. Lubarchuk, F.N. (2016) ‘Methodical ensuring scientific and methodical provisions the reduced tests of samples of rocket and artillery arms’, Izv. of the RARAS, (1), pp. 80–86 (in Russ.).
  3. Arsenev, V.N., Dubinin, D.P. (2020) ‘Justification of the method of estimation of characteristics of complex systems with a limited number field tests’, S&T Center of the Defense Complex "Kompas", (2), pp. 30–36 (in Russ.).
  4. Kozlov, N.N., Krasny, V.P., Reshetnikov, A.V. (2015) ‘The features of modern methodology for testing aerospace defense weapon systems’, Military Thought, (6), pp. 42–50 (in Russ.).
  5. Aleksandrovskaya, L.N., Kirillin, A.V., Kerber, O.B., Iosifov, P.A. (2018) ‘An approach to reducing the number of field tests for developing complex technical systems’, Proc. FSUE NPTsAP. Systems and Control Devices, (3), pp. 59–63 (in Russ.).
  6. Kofman, V.M. (2018) ‘Mathematical model of the inverse thermogasodynamic task for experimental and calculating estimation of the indicators of the efficiency of the combustion chamber and the turbine of the experimental GTE gas generator’, PNRPU Aerospace Engineering Bull., (53), pp. 97–108 (in Russ.).
  7. Naydeonov, V.G., Krupsky, K.A., Bochkarev, A.V. (2015) ‘Methodical approach to the estimation of the required amount of natural experiments testing complex samples weapons, military and special equipment’, Armament and Economics, (1), pp. 4–11 (in Russ.).
  8. Antonova, T.V. (2015) ‘Methods of identifying a parameter in the kernel of the first kind equation of the convolution type on the class of functions with discontinuities’, Numerical Analysis and Applications, 18(2), pp. 107–120. doi: 10.15372/SJNM20150201 (in Russ.).
  9. Voskoboynikov, Yu.E., Boeva, V.A. (2020) ‘Non-parametric identification algorithms for complex engineering systems’, Sci. Bull. of the NSTU, 80(4), pp. 47–64. doi: 10.17212/1814-1196-2020-4-47-64 (in Russ.).
  10. Danilov, A.M., Garkina, I.A. (2018) ‘Retrospective identification of a dynamic system’, Pridneprovsky Sci. Bull., 6(4), pp. 62–65 (in Russ.).
  11. Shishkina, A.V. (2018) ‘On the analysis of nonparametric identification algorithms’, Proc. Reshetnev Readings, 2(22), pp. 180–181 (in Russ.).
  12. Yareshchenko, D.I. (2020) ‘About non-parametric identification of partial-parametred discretecontinuous process’, Siberian J. of Sci. and Tech., 21(1), pp. 47–53. doi: 10.31772/2587-6066-2020-21-1-47-53 (in Russ.).
  13. Korneeva, A.A., Chernova, S.S., Shishkina, A.V. (2017) ‘Non-parametric algorithms of reconstruction of mutually ambiguous functions from observations’, Siberian J. of Sci. and Tech., 18(3), pp. 510–519 (in Russ.).
  14. Vozhov, S.S., Chimitova, E.V. (2018) ‘Parametric and non-parametric identification of the distribution law by interval data’, Metrologiya, (1), pp. 3–16 (in Russ.).
  15. Gusenitsa, Ya.N. (2022) ‘Solution of the nonparametric identification equation in a dynamic system based on the hyperdelta approximation’, Problems of Mechanical Engineering and Machine Reliability, 51(1), pp. 80–85 (in Russ.).
  16. Wang, H. (2016) ‘Machine learning basics’, in Goodfellow, I., Bengio, Y. and Courville A. (eds.) Deep Learning. pp. 98–164.

Постоянный адрес статьи:
http://swsys.ru/index.php?page=article&id=5007
Версия для печати
Статья опубликована в выпуске журнала № 2 за 2023 год. [ на стр. 320-326 ]

Возможно, Вас заинтересуют следующие статьи схожих тематик: