CC BY
77
УДК 519.2
с1о1:10.15217/1ззп1684-8853.2015.5.55
метод обработки
V V
неоднородной статистическом информации
о характеристиках точности системы управления
A. А. Ардашова, канд. техн. наук, старший научный сотрудник
B. Н. Арсеньева, доктор техн. наук, профессор
C. Б. Силантьева, канд. техн. наук, профессор
аВоенно-космическая академия им. А. Ф. Можайского, Санкт-Петербург, РФ
Введение: в процессе опытной отработки систем управления перспективных объектов невозможно обеспечить полную идентичность условий испытаний отдельных образцов из-за проведения доработок, изменения граничных условий и т. д. Одним из путей устранения неоднородности информации, получаемой в процессе испытаний, является приведение результатов испытаний отдельных образцов к некоторым заранее заданным условиям. Качество оценок характеристик системы управления, получаемых по объединенной таким образом выборке, существенно зависит от точности используемых операторов приведения. Цель: повышение точности определения операторов приведения по сравнению с известными методами решения этой задачи. Результаты: предложен новый подход к объединению неоднородных опытных данных, позволивший повысить точность приведения результатов испытаний к единым условиям и качество оценок характеристик системы. В основу определения оператора приведения положены условия полного совпадения математических ожиданий и максимальной близости ковариационных матриц, характеризующих точность системы в различных условиях. Таким образом обеспечен наиболее полный учет ограниченных опытных данных, полученных в процессе опытной отработки системы управления. Приведен пример оценивания характеристик точности системы управления предложенным методом. Практическая значимость: применение полученных результатов позволяет повысить точность оценок характеристик системы управления.
Ключевые слова — система управления, характеристики точности, опытные образцы, условия испытаний, неоднородная статистическая информация.
Введение
Заключительным этапом процесса создания системы управления (СУ) любого объекта являются натурные испытания опытных образцов (ОО). В ходе их проведения ведутся доработки СУ, корректируются алгоритмы управления, могут изменяться граничные условия, полезные нагрузки и т. д. [1-8]. Из-за отличий условий испытаний отдельных опытных образцов получаемые данные о характеристиках СУ являются неоднородными [9-13]. Одним из путей устранения неоднородности является приведение результатов испытаний отдельных образцов к некоторым заранее заданным условиям испытаний [2, 8-10]. Качество оценок характеристик СУ, получаемых по объединенной таким образом выборке, существенно зависит от точности используемых операторов приведения. Целью работы является повышение точности определения операторов приведения по сравнению с известными методами решения этой задачи [8-10].
Постановка задачи
Имеется модель СУ в виде системы нелинейных дифференциальных уравнений
л
л X л л
— = Р(Х,и, х, г), Х(^) = Хо, (1)
где значок «л» используется для отличия случайной величины от детерминированной;
л л
Х(г) = Хн (г) + Д Х(г) е К" — вектор фазовых координат СУ в момент времени г; Хн(г) — его нол
минальное значение; Д Х(г) — вектор случайных вариаций фазовых координат СУ относительно номинального значения Хн(г); Р(-) — в общем случае нелинейная вектор-функция; и е ИД —
л
вектор-функция программ управления; X е Кт —
вектор случайных параметров системы; г е (г0, гк); t0 и £к — моменты времени начала и окончания функционирования СУ.
Полагается [4], что распределение вектора вал
риаций фазовых координат СУ ДХ(гк) в конечной
точке траектории движения является многомерным нормальным ЩМдХ, КдХ) с математическим ожиданием МдХ и ковариационной матрицей КдХ. Величины МдХ и КдХ характеризуют область рассеивания фазовых координат СУ в момент времени £к и называются точностными характеристиками (характеристиками точности) СУ. Они зависят от условий испытаний ОО, которые могут отличаться граничными значениями, полезными нагрузками, значениями параметров СУ и т. д.
Пусть в некоторых первых условиях испытано Ь1 опытных образцов, а во вторых, отличных от первых, 12. Результаты испытаний представлены
соответственно множествами АХ1у-, - = и ДХ2у, ]' = 1,¿2 значений векторов вариаций фазовых координат СУ в конечных точках траекторий движения.
Необходимо получить оценки M12 и К12 математического ожидания МдХ и ковариационной матрицы КдХ, характеризующие точность СУ в первых условиях, по результатам испытаний ОО в первых и вторых условиях испытаний.
Определение оператора приведения
Полагается, что в первых и вторых условиях
Л Л
испытаний векторы А Х1 и А Х2 вариаций фазовых координат СУ в конечных точках траекторий движения связаны линейной зависимостью:
А Х1 = Р12Д XX 2, (2)
где P12 — матрица размерности пхп, которая в дальнейшем называется оператором приведения результатов испытаний ОО во вторых условиях к первым условиям.
В силу (2) матрица P12 должна удовлетворять двум уравнениям:
М л = Р12М Л ;
А Х1 А Х2
■»Т
К
А Х1
Р^К Л Р1Т2.
А Х2
(3)
(4)
Оценки параметров М Л , М Л , К Л и К Л
АХ1 АХ2 АХ1 АХ2
могут быть получены путем многократных статистических испытаний модели (1). Поскольку число испытаний может быть сколь угодно большим, то оценки обозначаются так же, как сами параметры.
Задача определения оператора приведения P12 сводится к совместному решению линейного векторного уравнения (3) и нелинейного матричного уравнения (4). В качестве показателей точности полученного решения используются величины
8М =
М Л А Х1 " Р12М Л А Х2
М Л АХ1
100% и
=
К Л " Р12К РТ Л Р12
А Х1 А Х2
К Л
АХ1
100%.
Можно показать [4], что система уравнений (3), (4) приводится к виду
м = ^2; (5)
QQT = I,
(6)
где
= б:1^2, м
А Х;
Q = Б"12
8тР1282Б
12. 2 ;
8;Б;8; = К Л ; S¿ и D¿ — ортогональная и диа-
А Х;
гональная матрицы, состоящие из собственных векторов и собственных значений матрицы К Л соответственно; I = 1, 2. А Х;
В общем случае точного решения системы уравнений (3), (4) [и соответствующей системы уравнений (5), (6)] может не существовать. Это может иметь место, когда реальная зависимость
ЛЛ
между векторами АХ1 и АХ2 отличается от линейной (2). В связи с этим предлагается следующий метод приближенного решения системы уравнений (3), (4).
Из (5) видно, что существует бесчисленное множество матриц Q, удовлетворяющих этому уравнению.
Пусть Q0 — любая матрица той же размерности, что и Q. Тогда общее решение уравнения (5) можно найти, минимизируя функционал J = - Q) х - Q)T] при условии ц1 = Qц2, где 1г[] — функция определения следа матрицы.
Решение этой задача, полученное методом неопределенных множителей Лагранжа, имеет вид
(
Q = Qc
л
Т
МзМ2
Т
Т
МзМ2
Т
М^2 .
(7)
Легко проверить, что при любой матрице Q0 выполняется уравнение (5). В частном случае можно взять в качестве Q0 единичную матрицу I.
Матрица (7) подставляется в уравнение (6). Очевидно, что в общем случае она не обращает это уравнение в тождество. Однако путем поочередного изменения элементов матрицы Q0 с небольшим шагом (0,001-0,01) в сторону уменьшения величины 5к можно получить достаточно точное с точки зрения показателей 5М и 5К решение поставленной задачи.
При этом любая полученная таким образом матрица Q обеспечивает строгое выполнение условия (5), а соответствующая ей матрица
Р12 = 81Б1^
12
8Т
(8)
является решением уравнения (3). Оценивание характеристик точности СУ
Результаты испытаний ОО во вторых условиях, приведенные к первым условиям испытаний, Р12АХ2У, ]' = 1,¿2, и результаты испытаний ОО в первых условиях АХ^-, - = 1^, рассматриваются как выборка из одной генеральной совокупности с законом распределения ЩМ^, К АХ 1). Тогда оценки М12 и К12 математического ожидания М^ 1 и ковариационной матрицы К АХ 1, ха-
рактеризующие точность СУ в первых условиях испытаний и учитывающие результаты испытаний ОО в первых и вторых условиях, определяются по формулам
(
М12 =-
I! +12 1
1 ДХц + Р12 1 ДХ2У 1=1 1=1
\
11 +12 К12 =
¿1Й1 л + ¿2Р12М
Д Х1
Д Х2
(9)
1 [ДХц- - М
1=1'
Д Х1
¿1 + ¿2 - 2
ДХ1у - М
Д Х1
+Р
12 •
1 |ДХ21 - М
1=1
Д Х2
ДХ2у - М
Д Х2
Рт
Р12
(10)
где мдХ! =11ДХи; дХ2 =11ДХ21 •
¿1 |=1 ¿2 I=1
Можно показать, что дисперсии элементов оценок М12 и К12 меньше дисперсий соответствующих элементов оценок, полученных только по результатам испытаний ОО в первых условиях.
Пример. В качестве исследуемого объекта рассматривается ракета-носитель (РН), предназначенная для выведения космических аппаратов (КА) на две различные орбиты. Условия пусков РН отличаются координатами точки выведения.
л
Трехмерный вектор ДХ состоит из отклонений фазовых координат РН от расчетных значений в конечных точках траекторий выведения. Отклонения измеряются в метрах.
Априорные оценки характеристик точности СУ РН и другие исходные данные заимствованы из работы [10]:
МДХ1 = [25 6 37Г; мдХ2 = [24 8 123
КДХ1
КДХ 2
814081 561 1403298 561 1497 1031 1403298 1031 2428612
776354 548556 547578" 548556 388292 387339 547578 387339 387377
Соответствующий этим данным оператор приведения, полученный по формуле (8), имеет вид
Р1
12
1,22 0,23 -0,51 0,38 -0,83 0,30 0,62 -0,65 2,27
Величины, характеризующие точность данного решения: 5М = 0, 5К = 0,04 %. Для сравнения оператор приведения, найденный путем решения уравнения (4) (в соответствии с [8]), дает 5М = 43 % и 5К = 0, а в соответствии с [9] — 5М = 18 % и 5К =0,06 %.
Поскольку реальные пуски РН не проводились, то необходимые статистические данные были получены путем испытаний модели (1). В первых условиях испытаний проведено 6 модельных пусков, результаты которых ДХ11 = = [-1276 -56 -2290]т; ДХ12 = [26 -24 51]т; ДХ13 = [1427 36 2456]т; ДХ14 = [-729 11 -1332]т; ДХ15 = [684 17 1184]т; ДХ16 = [54 39 189]т, а во вторых — 8 пусков с результатами ДХ21 = = [-1257 -933 -910]т; ДХ22 = [32 10 29]т; ДХ23 = [1388 999 974]т; ДХ24 = [-730 -521 -540]т; ДХ23 = [667 482 471]т; ДХ26 = [70 8 8 81]т; ДХ25 = [-446 -308 -316]т; ДХ28 = [429 258 289]т.
По этим данным рассчитаны оценки
ММДХ1 = [31 4 43]т; 1МДХ2 = [19 9 9^ математических ожиданий МДХ и МдХ2 соответственно.
Тогда оценки М12 и К12 характеристик точности выведения РН в первых условиях пусков, учитывающие результаты пусков РН в первых и вторых условиях, в соответствии с формулами (9), (10) примут вид
М12 =[26 3 34^;
К12 =
810956 6662 1419402 6662 820 12178 1419402 12178 2486901
Очевидно, что дисперсии элементов оценок М12 и К12, учитывающих результаты 14 пусков РН в первых и вторых условиях, будут меньше дисперсий соответствующих элементов оценок МдХх и КдХполученных только по результатам шест и пусков РН в первых условиях.
Заключение
Приведение разнородных опытных данных к единым условиям испытаний позволяет повысить качество оценок характеристик точности СУ. Результат решения этой задачи существенно зависит от ошибок определения оператора приведения. Предложенный в данной статье подход, как показали многочисленные исследования, во всех случаях обеспечивает нулевую ошибку 5М и, как правило, малое значение ошибки 5К. Он может быть использован при решении ряда смежных задач. Так, результаты испытаний ОО в первых условиях могут быть пересчитаны на вторые. Или же, в более общем случае, все
х
испытания в первых и вторых условиях могут быть использованы для повышения качества оценивания характеристик точности СУ в некоторых, отличных от первых двух, третьих условиях испытаний. Возможны и другие, достаточно очевидные постановки задач, связанных с обработкой неоднородной статистической информации.
Литература
1. Поляков А. П. и др. Справочник по эксплуатации космических средств / под ред. А. П. Полякова. 2-е изд., перераб. и доп. — СПб.: ВКА им. А. Ф. Можайского, 2006. — 758 с.
2. Миронов В. И. Эффективность, надежность и испытания систем управления: учеб. пособие / МО СССР. — М., 1981. — 200 с.
3. Элементы теории испытаний и контроля технических систем / под ред. Р. М. Юсупова. — Л.: Энергия, 1978. — 192 с.
4. Эльясберг П. Е. Определение движения по результатам измерений. — М.: Наука, 1976. — 415 с.
5. Буренок В. М., Найденов В. М. Испытательная база: выход из кризиса // Воздушно-космическая оборона. 2009. № 1(44). С. 18-25.
6. Шаракшанэ А. С., Халецкий А. К., Морозов И. А.
Оценка характеристик сложных автоматизированных систем. — М.: Машиностроение, 1993. — 272 с.
7. Кринецкий Е. И., Александровская Л. Н., Шаронов А. В., Голубков А. С. Летные испытания ракет и космических аппаратов. — М.: Машиностроение, 1979. — 464 с.
8. Миронов В. И. Задача приведения вариаций фазовых координат нелинейных динамических систем к заданным условиям испытаний // Изв. АН СССР. Техн. кибернетика. 1970. № 3. С. 31-38.
9. Арсеньев В. Н. Определение характеристик точности системы по результатам ее испытаний в различных условиях испытаний // Изв. АН РФ. Техн. кибернетика. 1992. № 2. С. 118-121.
10. Арсеньев В. Н. Оценивание характеристик точности системы управления ракеты-носителя по результатам пусков в различных условиях // Изв. вузов. Приборостроение. 2015. № 1. С. 27-32.
11. Половко А. М., Гуров С. В. Основы теории надежности. — СПб.: БХВ-Петербург, 2006. — 704 с.
12. Вентцель Е. С. Теория вероятностей. — М.: Академия, 2003. — 576 с.
13. Пугачев В. С. Теория вероятностей и математическая статистика. — М.: Физматлит, 2002. — 496 с.
UDC 519.2
doi:10.15217/issn1684-8853.2015.5.55
Processing Miscellaneous Statistical Information about Control System Accuracy Characteristics
A. A. Ardashova, PhD, Tech., Senior Research Associate, avgust.ar.@yandex.ru V. N. Arsenieva, Dr. Sc., Tech., Professor, vladar56@mail.ru S. B. Silantyeva, PhD, Tech., Professor, silantev2008@yandex.ru
aA. F. Mozhaiskii Military Space Academy, 13, Zhdanovskaia St., 197198, Saint-Petersburg, Russian Federation
Introduction: When testing control systems, it is impossible to provide that the test conditions for various samples are fully identical, because of modifications, changes in the boundary conditions, etc. A way to remove the information heterogeneity is reducing the test results for different samples to certain conditions determined in advance. The estimates of the control system characteristics obtained from the selection united in this way significantly depends on the accuracy of the reduction operators. Purpose: Increasing the definition accuracy of the reduction operators as compared to the known methods of solving this task. Results: A new approach to amalgamate miscellaneous test data has been proposed which helps to increase the accuracy of reducing the test results to unified conditions and improve the quality of estimating the system characteristics. The reduction operator is determined on the base of full coincidence of population means and the maximum proximity of the correlation matrixes characterizing the accuracy of the system under various conditions. This provides the best possible consideration of limited test data obtained in the course of test functioning of the control system. An example of estimating the accuracy characteristics of the control system is given. Practical relevance: Applying the results obtained in this paper can help to increase the accuracy of estimating the characteristics of the control system.
Keywords — Control System, Accuracy Characteristics, Prototypes, Test Condition, Miscellaneous Statistical Information.
References
1. Polyakov A. P., et al. Spravochnik po ekspluatatsii kosmicheskikh sredstv [The Reference Book on Operation of Space Means]. Saint-Petersburg, VKA im. A. F. Mo-zhaiskogo Publ., 2006. 758 p. (In Russian).
2. Mironov V. I. Effektivnost', nadezhnost' i ispytaniia sistem upravleniia [Efficiency, Reliability and Tests of Control Systems]. MO SSSR Publ., 1981. 200 p. (In Russian).
3. Elementy teorii ispytanii i kontrolia tekhnicheskikh sistem [Elements of the Theory of Tests and Control of Technical Systems]. Ed. R. M. Yusupov. Saint-Petersburg, Energiia Publ., 1978. 192 p. (In Russian).
4. Eliasberg P. E. Opredelenie dvizheniia po rezul'tatam izme-renii [Definition of the Movement on Results of Measurements]. Moscow, Nauka Publ., 1976. 415 p. (In Russian).
5. Burenok V. M., Naidyonov V. M. Test Base: Recovery from the Crisis. Vozdushno-kosmicheskaia oborona, 2009, no. 1(44), pp. 18-25 (In Russian).
6. Sharakshane A. S., Haletsky A. K., Morozov I. A. Otsenka kharakteristik slozhnykh avtomatizirovannykh sistem [Estimation of Characteristics of the Difficult Automated Systems]. Moscow, Mashinostroenie Publ., 1993. 272 p. (In Russian).
7. Krinetsky E. I., Aleksandrovskaia L. N., Sharonov A. V., Golubkov A. S. Letnye ispytaniia raket i kosmicheskikh ap-paratov [Flight Tests of Rockets and Spacecraft]. Moscow, Mashinostroenie Publ., 1979. 464 p. (In Russian).
8. Mironov V. I. The Task of Reduction of Nonlinear Dynamic Systems Phase Coordinates Variations to the Set Conditions of Test. Izvestiia AN SSSR. Tekhn. kibernetika, 1970, no. 3, pp. 31-38 (In Russian).
9. Arseniev V. N. Opredeleniye harakteristik tochnosti systemi po rezulytatam ee ispitaniy v razlichnih usloviyah [Definition of System Accuracy Characteristics by Results of Its Tests In Various Conditions]. Izvestiia AN RF. Tekhn. kibernetika, 1992, no. 2, pp. 118-121 (In Russian).
10. Arseniev V. N. Estimation of Carrier Rocket Control System Accuracy Characteristics by Results of Start-up in Various Conditions. Izvestiia vuzov. Priborostroenie, 2015, no. 1, pp. 27-32 (In Russian).
11. Polovko A. M., Gurov S. V. Osnovy teorii nadezhnosti [Bases of The Theory of Reliability]. Saint-Petersburg, BKhV-Pe-terburg Publ., 2006. 704 p. (In Russian).
12. Ventcel E. S. Teoriia veroiatnostei [Probability Theory]. Moscow, Akademiia Publ., 2003. 576 p. (In Russian).
13. Pugachev B. C. Teoriia veroiatnostei i matematicheskaia statistika [Probability Theory and Mathematical Statistics]. Moscow, Fizmatlit Publ., 2002. 496 p. (In Russian).
Уважаемые авторы!
При подготовке рукописей статей необходимо руководствоваться следующими рекомендациями.
Статьи должны содержать изложение новых научных результатов. Название статьи должно быть кратким, но информативным. В названии недопустимо использование сокращений, кроме самых общепринятых (РАН, РФ, САПР и т. п.).
Объем статьи (текст, таблицы, иллюстрации и библиография) не должен превышать эквивалента в 20 страниц, напечатанных на бумаге формата A4 на одной стороне через 1,5 интервала Word шрифтом Times New Roman размером 13, поля не менее двух сантиметров.
Обязательными элементами оформления статьи являются: индекс УДК, заглавие, инициалы и фамилия автора (авторов), ученая степень, звание (при отсутствии — должность), полное название организации, аннотация и ключевые слова на русском и английском языках, электронные адреса авторов, которые по требованию ВАК должны быть опубликованы на страницах журнала. При написании аннотации не используйте аббревиатур и не делайте ссылок на источники в списке литературы.
Статьи авторов, не имеющих ученой степени, рекомендуется публиковать в соавторстве с научным руководителем, наличие подписи научного руководителя на рукописи обязательно; в случае самостоятельной публикации обязательно предоставляйте заверенную по месту работы рекомендацию научного руководителя с указанием его фамилии, имени, отчества, места работы, должности, ученого звания, ученой степени — эта информация будет опубликована в ссылке на первой странице.
Формулы набирайте в Word, не используя формульный редактор (Mathtype или Equation), при необходимости можно использовать формульный редактор; для набора одной формулы не используйте два редактора; при наборе формул в формульном редакторе знаки препинания, ограничивающие формулу, набирайте вместе с формулой; для установки размера шрифта никогда не пользуйтесь вкладкой Other..., используйте заводские установки редактора, не подгоняйте размер символов в формулах под размер шрифта в тексте статьи, не растягивайте и не сжимайте мышью формулы, вставленные в текст; в формулах не отделяйте пробелами знаки: + = -.
Для набора формул в Word никогда не используйте Конструктор (на верхней панели: «Работа с формулами» — «Конструктор»), так как этот ресурс предназначен только для внутреннего использования в Word и не поддерживается программами, предназначенными для изготовления оригинал-макета журнала.
При наборе символов в тексте помните, что символы, обозначаемые латинскими буквами, набираются светлым курсивом, русскими и греческими — светлым прямым, векторы и матрицы — прямым полужирным шрифтом.
Иллюстрации в текст не заверстываются и предоставляются отдельными исходными файлами, поддающимися редактированию:
— рисунки, графики, диаграммы, блок-схемы предоставляйте в виде отдельных исходных файлов, поддающихся редактированию, используя векторные программы: Visio 4, 5, 2002-2003 (*.vsd); Coreldraw (*.cdr); Excel (*.xls); Word (*.doc); Adobelllustrator (*.ai); AutoCad (*.dxf); Matlab (*.ps, *.pdf или экспорт в формат *.ai);
— если редактор, в котором Вы изготавливаете рисунок, не позволяет сохранить в векторном формате, используйте функцию экспорта (только по отношению к исходному рисунку), например, в формат *.ai, *.esp, *.wmf, *.emf, *.svg;
— фото и растровые — в формате *.tif, *.png с максимальным разрешением (не менее 300 pixels/inch).
Наличие подрисуночных подписей обязательно (желательно не повторяющих дословно комментарии к рисункам в тексте статьи).
В редакцию предоставляются:
— сведения об авторе (фамилия, имя, отчество, место работы, должность, ученое звание, учебное заведение и год его окончания, ученая степень и год защиты диссертации, область научных интересов, количество научных публикаций, домашний и служебный адреса и телефоны, e-mail), фото авторов: анфас, в темной одежде на белом фоне, должны быть видны плечи и грудь, высокая степень четкости изображения без теней и отблесков на лице, фото можно представить в электронном виде в формате *.tif, *.png с максимальным разрешением — не менее 300 pixels/inch при минимальном размере фото 40x55 мм;
— экспертное заключение.
Список литературы составляется по порядку ссылок в тексте и оформляется следующим образом:
— для книг и сборников — фамилия и инициалы авторов, полное название книги (сборника), город, издательство, год, общее количество страниц;
— для журнальных статей — фамилия и инициалы авторов, полное название статьи, название журнала, год издания, номер журнала, номера страниц;
— ссылки на иностранную литературу следует давать на языке оригинала без сокращений;
— при использовании web-материалов указывайте адрес сайта и дату обращения.
Список литературы оформляйте двумя отдельными блоками по образцам lit.dot на сайте журнала (http://i-us.ru/paperrules) по разным стандартам: Литература — СИБИД РФ, References — один из мировых стандартов.
Более подробно правила подготовки текста с образцами изложены на нашем сайте в разделе «Оформление статей».
Контакты
Куда: 190000, Санкт-Петербург, Б. Морская ул., д. 67, ГУАП, РИЦ Кому: Редакция журнала «Информационно-управляющие системы» Тел.: (812) 494-70-02 Эл. почта: ius.spb@gmail.com Сайт: www.i-us.ru
