УДК 681.11.054
DOI: 10.12737/article_5a02f9fe97a5a7.40674404
С.И. Досько, В.В. Киренков, А.А. Молчанов
МЕТОД ФОРМИРОВАНИЯ ОПТИМАЛЬНОГО ТЕСТИРУЮЩЕГО ВОЗДЕЙСТВИЯ ПРИ ИСПЫТАНИЯХ ДИНАМИЧЕСКИХ СИСТЕМ
Описан метод формирования оптимальных тестирующих воздействий при проведении проверок и экспериментов с динамическими системами при их доводке, диагностике, приемосдаточных испытаниях. В качестве критерия оптимальности используется известный принцип максимума следа информационной матрицы Фишера, который
представляет собой осредненный по оцениваемым параметрам критерий достоверности (точности) проведенных оценок.
Ключевые слова: идентификация, входное воздействие, принцип двойственности, теория оптимального управления, след информационной матрицы, динамические системы.
S.I. Dosko, V.V. Kirenkov, A.A. Molchanov
FORMATION METHOD OF OPTIMUM TESTING IMPACT AT DYNAMIC SYSTEMS TESTS
The paper reports the description of the formation method of optimum testing influences in the course of carrying out tests and experiments with dynamic systems, diagnostics, acceptance tests. A problem is posed in the following way: there is assigned a time interval intended for system tests before a mode of its normal operation, a simulator of the system and input effects maximum allowable according to a module. In what way in this time interval, often limited, one should excite a system in order to ensure, in a certain sense, optimum conditions for system identification?
One of possible methods, for the substantiated solution such a problem, based, from the mathematical point of view, on the duality principle of problems of experimental data processing (identification) and problems of optimum management is considered. As a criterion of optimality there is used a well-known principle of
Введение
Задачам обработки и интерпретации опытных данных, полученных при испытаниях динамических систем, всегда уделялось большое внимание как в теоретическом, так и в практическом плане. Среди этих задач все большее место занимают так называемые обратные задачи [1; 2], одним из возможных путей совершенствования методологии решения которых является привлечение принципа двойственности - методов теории оптимального управления. Основанием для такого подхода является то, что обе задачи (идентификации и управления) преследуют
maximum of Fisher information matrix trace, which is a criterion of assessments carried out averaged on parameters assessed of reliability.
Particular attention is paid to the identity of methods for the solution of some inverse problems arising at the assessment of test results of complex dynamic systems with classical problems of optimum management theory. In particular, it is shown, in what way a problem of the optimum testing effect formation in the course of active experiments can be solved with the use of the methodology for the solution of optimum management standard problems based on the well-known Pontryagin maximum.
Key words: identification, input action, duality principle, optimum management theory, information matrix trace, dynamic systems.
с математической точки зрения идентичную конечную цель - нахождение экстремума квадратичного функционала от каких-то характеристик (обычно невязок) решаемой задачи.
В статье [10] показано, как использование принципа двойственности может способствовать решению одной из типовых обратных задач, возникающих при оценке результатов штатной эксплуатации динамической системы (в режиме пассивного эксперимента). В настоящей статье с этой точки зрения рассматривается один принципиальный вопрос,
возникающим и при проведении активного эксперимента. Речь идет о формировании
оптимального тестирующего воздействия.
Математические модели объекта контроля
Традиционная форма математических моделей объектов или процессов, подлежащих испытаниям и последующей идентификации, обычно представляется в следующем виде (форма Коши):
Х=ДХ,и^), (1)
где X(t) - я-мерный вектор состояния; U(t) -т-мерный вектор управления; F(t) - р-мерный вектор возмущений.
Линеаризация этой модели посредством ее представления в отклонениях от опорного программного движения может быть записана в следующем стандартном, в частности для теории оптимальных процессов, виде:
dX
— = A(t) X+B(t)U+D(t )F, dt
(2)
где A(n x n) определяющая
квадратная матрица, изменение состояния
Выбор тестирующего воздействия и решение
Основной задачей, которая решается в процессе активного эксперимента, является определение фактических значений и отклонений от расчетных значений какого-то комплекса наиболее важных параметров регулятора системы а/, в частности параметров закона управления.
В такой постановке, помимо решения задачи собственно идентификации, всегда встает и другая упомянутая выше задача: как выбрать тестирующие воздействия, которые с учетом естественных ограничений могли бы обеспечить наилучшие (в каком-то смысле) условия идентифицируемости системы? Эта задача может решаться с привлечением методов теории оптимальных процессов и теории чувствительности [2; 4; 5]. Пусть имеется динамическая система, уравнение которой в матрично-векторной форме по аналогии с (1) может быть представлено в следующем виде:
Х=ДХ,а,и^), (4)
и измерении
нестационарной линейной системы; В(п x m), D (п xp) - матрицы управляющих U и возмущающих F воздействий
соответственно.
Целесообразность перехода от формы (1) к форме (2) обусловлена тем, что для линейных систем разработаны методы решения задач оптимального управления [24], охватывающие различные, в том числе и квадратичные, критерии качества.
Уравнение измерений также представляется в традиционной для теории оптимального управления форме:
Y(t) = C(t) X + L(t)U + V(t), (3)
где Y(l x 1), С (l x n), L(p x r) - вектор и матрица измерений соответственно; V (l x 1) - случайный процесс, характеризующий ошибки измерений.
задачи идентификации
где а[т х 1] - упомянутый выше вектор оцениваемых параметров.
В заданном интервале времени (0; Т) система подвергается испытаниям, в процессе которых на ее вход подается внешнее тестирующее воздействие F(t). В ходе испытаний производится регистрация части фазовых координат в соответствии с уравнением измерений (3):
У (?) = С(1) X + У(г). (5)
Решение задачи идентификации подразумевает в общем случае оценку параметров, входящих в правую часть соотношения (4), в частности неизвестных элементов вектора а и входных воздействий F. Условия проведения активного эксперимента методически облегчают эту задачу, ибо предполагают задание в качестве F(t) известных, заранее выбранных тестирующих воздействий. С учетом этого задача идентификации элементов вектора а может решаться следующим образом.
Первый этап - ее линеаризация место при относительно опорного решения, имеющего
X[a0 +Aa, F(t)} = X[a0 F(t)] + W(t)Aa, где Aa[m x 1] - вектор отклонений оцениваемых параметров; W(t)(n x m) - так называемая матрица коэффициентов
номинальных
значениях
элементов а, а именно: 0 (6) чувствительности, определяемая по правилам дифференцирования матриц и векторов следующим образом:
W (t) =
дх
да
дхг дхг дхг
дах да2 да m
дх2 дх2 дх2
дах да2 да m
дх n дх n дх n
дах да да m
f dx
(7)
Элементы матрицы чувствительности W(t) могут быть получены после этого [5] путем дифференцирования по а соотношения (4):
щол = л. (8)
дх да да дх да
Таким образом, соотношения (4) при а = ао и (8) при W(0)=0 определяют значения функций чувствительности при внешнем воздействии В частности, при условии
= 0 они могут быть сведены к типовому
операторному уравнению,
рассматриваемому при классических методах решения обратных задач [1].
Полученные соотношения позволяют свести исходную задачу в разряд линейных задач, для которых разработаны и исследованы методы идентификации, дающие несмещенные и эффективные оценки. В предположении гауссовского характера распределения ошибок измерений V(t) решение для вектора Аа ищется исходя из условия min квадратичного функционала, имеющего вид
I = | [(7 (г) - С (г) X (/, а0 0 - С (г )Ш (г )Аа)Т (У (г) - С (г) X (/, а0) - С (г )Ш (г)Аа)]Уг
0
гд е Т (в подынтегральном выражении и далее) - символ транспонирования. Минимизация функционала (9) по Аа приводит к следующему матрично-
Ш(г)Т С (г)ТС(г)Ж(г)]Уг Аа Т(г)Т С (г)АУ у = 0:
(9)
векторному уравнению для определения составляющих этого вектора:
(10)
где АУ = У (г) - С (г) X (/, аи) - вектор отклонений фактических результатов измерений от модельных значений, соответствующих Аа =0.
Уравнение (10) представляет собой алгебраическое уравнение т-го порядка относительно составляющих вектора Аа (Аа1, Аа2 ... Аам), решение которого в общем виде записывается следующим образом:
Аа = Б'АУ, (11)
где B(m х m) = T [W(t)T C(t)TC(t)W(t)]dt - так
0
называемая матрица идентификации (иногда ее называют информационной матрицей Фишера).
Элементы этой матрицы определяются характером изменения по времени и величиной функций чувствительности W(t), определяемых входными воздействиями F(t). В условиях активного эксперимента имеется возможность так формировать эти воздействия, чтобы обеспечить в каком-то
смысле оптимальные условия для определения элементов вектора Аа. В качестве этих оптимальных условий целесообразно выбрать такие, которые обеспечивают минимально возможные ошибки определения элементов этого вектора.
Как следует из соотношения (11), эти погрешности могут характеризоваться диагональными элементами матрицы В'1. Критерий, к которому нужно стремиться, состоит в том, чтобы сумма всех диагональных членов этой матрицы была минимальной в рамках дозволенных пределов изменения входных воздействий
Известно, что сумма диагональных элементов квадратной матрицы (или ее след равна сумме собственных чисел этой матрицы. Если обозначить собственные
числа матрицы В через Х1, Х2 ... то след матрицы В'1 и условия оптимальности запишутся следующим образом:
т |
8рБ4 = ^ — = шт.
1 К
(12)
При решении аналогичных задач условие минимизации следа обратной матрицы В'1 заменяют более простым условием - максимизации следа матрицы В, т.е:
БрБ = ^К = шах.
( 13)
Условие (13) при решении задач идентификации трактуется как условие максимума следа информационной матрицы Фишера [2].
Меандровый сигнал и принцип максимума Л.С.
Рассмотренный подход позволяет свести исходную задачу идентификации к задаче выбора оптимального управления процессами, использовав классический принцип максимума Л.С. Понтрягина [4]. Из этого принципа следует, что оптимальное входное воздействие для
рассматриваемой задачи имеет кусочно-постоянный (меандровый) вид, который может быть представлен следующим образом:
F (0 = ( 14)
где С - модуль предельно допустимого входного воздействия; signA(щ) - знак функции вспомогательных переменных ¥ метода Л.С. Понтрягина А(^).
Общий вид оптимального
тестирующего воздействия (в соответствии с
Модельный пример
Ниже приведен модельный пример формирования оптимального тестирующего воздействия для некой механической системы. Пусть аналоговая передаточная функция между входным воздействием (сила, момент) X и выходом -
Понтрягина
(14)) при продолжительности испытаний Т приведен на рисунке.
Аналитических методов определения времени изменения знака вспомогательных переменных ¥ не существует. Поэтому даже для сравнительно простых систем это возможно только применением
приближенных численных методов.
Рис. Оптимальное (по Л.С. Понтрягину) тестирующее воздействие
перемещением (угловым или линейным) (Ж=У/Х) принята в форме
к
Ж(Р) =Т 2 2+ ОТ* + 1 . (15) 1 р +21 фр +1
Расчетные значения параметров к, Т и реализуемые этим набором, составляли 1,0; 0,051с и 0,256. Используя методологию
2-преобразования, перейдем от аналоговой
У}] = 2хх[1 - 1] + 22х[ со значениями коэффициентов zl=0,379, Z2=0,319, Z5=0,908, Z4=-0,605. При моделировании амплитуда тестирующего воздействия вида (1 4) принималась равной С = ±1, шаг по времени АТ= 0,05с, интервал тестирования с учетом ограничения п=30
(1,5 с).
Эффективность того или иного воздействия оценивалась, как указано во введении, посредством сравнения коэффициентов
Данные таблицы свидетельствуют о том, что для выбранного примера - системы с передаточной функцией (15) и интервалом тестирования 1,5 с - оптимальным входным сигналом является меандровый сигнал с числом переключений N=3. Достигаемые
Заключение
В статье обращается внимание на идентичность методов решения некоторых обратных задач, возникающих при оценке результатов испытаний сложных динамических систем, с классическими задачами теории оптимального управления. В частности, показано, как задача формирования оптимальных тестирующих воздействий в процессе активных экспериментов может решаться с использованием методологии решения типовых задач оптимального управления на основе известного принципа максимума
формы (15) к системе разностных уравнений
+ У[ -1] + ¿4 УУ - 2] (16)
Аг.
влияния ошибок измерений D (К = ) на
точности оценок. Результаты расчетов для случая классического гармонического воздействия с собственной частотой & = 3 Гц, меандровых воздействий различного числа переключений N и ступенчатого N=0 приведены в таблице. Во всех случаях амплитуда входного воздействия ограничивалась значением С = ±1.
Таблица
при этом точности оценок параметров закона управления (коэффициентов zi ) существенно выше традиционных тестирующих сигналов - как в виде гармонической функции, так и в виде единичного ступенчатого воздействия.
Л.С. Понтрягина. Не исключено, что найдутся и другие обратные задачи, где окажется возможным эффективное использование методов оптимального управления. Рассмотренная в статье методология позволяет сделать объективное суждение о близости теоретически оптимальной и находящейся в эксплуатации реальной системы контроля, включая её тестирующие воздействия, и
целесообразности в этом отношении её доработки.
Коэффициенты влияния ошибок измерений K
Коэффициенты разностных уравнений
Гармоническая функция Q = 3 Гц Меандровое воздействие
N=0 N=1 N=2 N=3 N=4
Zi 0,938 1,001 0,454 0,346 0,269 0,316
Z2 1,312 1,538 0,693 0,525 0,399 0,443
Z3 0,602 1,732 0,734 0,537 0,537 0,586
Z4 0,709 1,086 0,491 0,363 0,369 0,372
СПИСОК ЛИТЕРАТУРЫ
1. Тихонов, А.Н. Методы решения некорректных задач / А.Н. Тихонов, В.Я. Арсенин. - М.: Наука, 1974.
2. Брайсон, А. Прикладная теория оптимальных процессов / А. Брайсон, Хо Ю-ши. - М.: Мир, 1972.
3. Ройтенберг, Я.И. Автоматическое управление / Я.И. Ройтенберг. - М.: Наука, 1978.
4. Понтрягин, Л.С. Математическая теория оптимальных процессов / Л.С. Понтрягин [и др.]. -М.: Физматгиз, 1961.
5. Розенвассер, Е.Н. Чувствительность систем управления / Е.Н. Розенвассер, Р.М. Юсупов. - М.: Наука, 1981.
6. Киренков, В.В. Типовые обратные задачи и методы их решения при оценке результатов испытаний изделий ракетно-космической техники / В.В. Киренков, С.И. Досько // Ракетно-космическая техника: сб. науч. тр. РКК «Энергия»
1. Tikhonov, A.N. Methods for Incorrect Problem Solution / A.N. Tikhonov, V.Ya. Arsenin. - M.: Science, 1974.
2. Braison, A. Applied Theory of Optimum Processes / A. Braison, Ho Yu-shi. - M.: Mir, 1972.
3. Roitenberg, Ya.I. Automated Management / Ya.I. Roitenberg. - M.: Science, 1978.
4. Pontryagin, L.S. Mathematical Theory of Optimum Processes / L.S. Pontryagin [et al.]. - M.: Physmath-giz, 1961.
5. Rosenwasser, E.N. Management System Sensitivity / E.N. Rosenwasser, R.M. Yusupov. - M.: Science, 1981.
6. Kirenkov, V.V. Standard inverse problems and methods for their solution at test results assessment of missile-space engineering / V.V. Kirenkov, S.I. Dosko // Missile-Space Engineering; Proceedings of Korolev
Сведения об авторах:
Досько Сергей Иванович, к.т.н., ст. науч. сотрудник Института конструкторско-технологической
информатики РАН, тел. 8-977-494-71-15, е-mail: [email protected].
Киренков Вениамин Васильевич, к.т.н, вед. науч. сотрудник РКК «Энергия» им. С.П. Королева, тел.: 8916-352-45-03.
Dosko Sergey Ivanovich, Can. Eng., Senior researcher, Institute of Design-Technological Informatics of RAS, е-mail: [email protected]
Kirenkov Veniamin Vasilievich, Can. Eng., Leading researcher of Korolev MSC "Energy", phone: 8-916-35245-03.
им. С.П. Королева. - 2014. - Сер. XII. - Вып. 3. - С. 99.
7. Калман, Р. Очерки по математической теории систем / Р. Калман [и др.]. - М.: Мир, 1971.
8. Тихонов, А.Н. Численные методы решения некорректных задач / А.Н. Тихонов, А.В. Гончарский [и др.]. - М.: Наука, 1990.
9. Киренков, В.В. Решение обратных задач оценки испытаний изделий ракетно-космической техники с использованием методов теории оптимального управления / В.В. Киренков, В.Г. Микитенко // Космическая техника и технологии. - 2015. - № 4 (11). - С. 74-81.
10. Досько, С.И. Использование принципа двойственности при решении обратных задач оценки результатов испытаний / С.И. Досько, В.В. Киренков, В.Г. Микитенко // Вестник Брянского государственного технического университета. - 2017. - № 2 (55).
MSC "Energy". - 2014. - Series XII. - Edition #. - pp. 99.
7. Kalman, R. Sketches on Mathematical Theory of Systems / R. Kalman [et al.]. - M.: Mir, 1971.
8. Tikhonov, A.N. Numetical Methods of Incorrect Problems Solution / A.N. Tikhonov, A.V. Gon-charsky [et al.]. - M.: Science, 1990.
9. Kirenkov, V.V. Inverse problems solution of missile-space product tests using methods of optimum management theory / V.V. Kirenkov, V.G. Mikitenko // Space Engineering and Technologies. - 2015. - No.4 (11). - pp. 74-81.
10. Dosko, S.I. Duality principle use at solution of test result inverse problems / S.I. Dosko, V.V. Kirenkov, V.G. Mikitenko // Bulletin of Bryansk State Technical University. 2017. - No.2 (55).
Статья поступила в редколлегию 5.10.17.
Рецензент: д.т.н., профессор ИКТИРАН
Шептунов С.А.
Молчанов Александр Александрович, аспирант Московского государственного технического университета им. Н.Э. Баумана, НТИМИ, научный сотрудник, тел.: 8-925-053-63-52, е-шаП: А1ехап[email protected].
Molchanov Alexander Alexandrovich, Post graduate student of Bauman State Technical University of Mos-cow,e-mail: [email protected].