CC BY
14УДК 534.6:534.843.12
идентификация импульсных шумов
в СЛУЖЕБНОМ МОДУЛЕ РОССиЙСКОГО СЕГМЕНТА МКС © 2017 г. Зайцев К.и., Половнев А.Л.
Ракетно-космическая корпорация «Энергия» имени С.П.Королёва (РКК «Энергия») Ул. Ленина, 4А, г. Королёв, Московская обл., Российская Федерация, 141070, e-mail: post@rsce.ru
Представлены результаты обработки бортовых импульсных шумов, зарегистрированных научной аппаратурой «Средства оперативного определения координат точки пробоя» в космическом эксперименте «Пробойд с помощью малогабаритных микрофонов, установленных внутри служебного модуля Российского сегмента Международной космической станции (РС МКС) во время работы экипажей МКС-41/42, МКС-43/44 и МКС-45/46.
Разработан и реализован в программном коде алгоритм распознавания источников импульсных акустических шумов, действующих на борту служебного модуля РС МКС.
Выполнены параметрические исследования разрешающей способности коэффициентов линейных корреляций для различных представлений спектров. Выполнен обоснованный выбор параметров спектрального и корреляционного анализов.
Сформирована корреляционная матрица 1 /3-октавных спектров, которая содержит четко выделенные группы импульсных шумов, включая группу, соответствующую спектрам сигналов имитатора пробоя, который размещался более чем в 100 точках объема служебного модуля, в т. ч. и в запанельном пространстве.
Разработанный алгоритм планируется использовать при доработке программного обеспечения научной аппаратуры по выделению в реальном масштабе времени импульсных сигналов пробоя с передачей телеметрической информации на Землю в экспедиции МКС-53/54.
Результаты исследований показывают принципиальную возможность идентификации импульсных источников шума на борту РС МКС.
Ключевые слова: идентификация импульсных шумов, космический эксперимент «Пробой», алгоритм распознавания источников импульсных акустических шумов, коэффициент линейной корреляции, 1 /3-октавный спектр, корреляционная матрица, импульсные сигналы пробоя.
IDENTIFYING PULSE NOISES IN THE SERVICE MODULE OF THE ISS RUSSIAN SEGMENT Zaytsev K.I., Polovnev A.L.
S.P. Korolev Rocket and Space Public Corporation Energia (RSC Energia) 4A Lenin str., Korolev, Moscow region, 141070, Russian Federation, e-mail:post@rsce.ru
The paper presents the results of processing the onboard pulsed noises recorded by the scientific equipment «Means of Real-time Determination of the Puncture-point Location» in the space experiment Proboi («Puncture») using small microphones installed inside the Service Module of the Russian Segment of the International Space Station (ISS RS) during the work of the crews of ISS-41/42, ISS-43/44 and ISS-45/46.
Developed and implemented in the form of a software code was an algorithm for recognizing sources of pulsed acoustic noises acting onboard the Service Module of the ISS RS.
Parametric studies were done on the resolution of coefficients of linear correlations for various spectral representations. An informed choice of parameters for spectral and correlation analyses was performed.
A correlation matrix of 1/3-octave spectra was generated, which contains clearly defined groups of pulsed noises, including a group corresponding to signal spectra of the puncture simulator, which was placed at more than 100 points within the volume of the Service Module, including spaces behind the panels.
There are plans to use the algorithm in updating the software for the scientific equipment for identifying in real-time pulsed signals of a puncture with transmission of telemetry data to Earth in Expedition ISS-53/54.
The results of the studies demonstrate that in principle it is feasible to identify sources of pulsed noise onboard the ISS RS.
Key words: identification of pulsed noises, space experiment «Proboi», algorithm for recognizing sources of pulsed acoustic noises, linear correlation coefficient, 1/3-octave spectrum, correlation matrix, pulsed signals of puncture.
Ж
ПОЛОВНЕВ А.Л.
ЗАЙЦЕВ Константин Иванович — аспирант, инженер 2 категории РКК «Энергия», e-mail: post@rsce.ru
ZAYTSEV Konstantin Ivanovich — Post-graduate, Engineer 2 category at RSC Energia, e-mail: post@rsce.ru
ПОЛОВНЕВ Антон Леонидович — кандидат физико-математических наук, старший научный сотрудник РКК «Энергия», e-mail: post@rsce.ru
POLOVNEV Anton Leonidovich — Candidate of Science (Physics and Mathematics), Senior research scientist at RSC Energia, e-mail: post@rsce.ru
ЗАЙЦЕВ К.И.
Введение
Космический эксперимент «Пробой» проводится с целью верификации метода оперативного определения координат точки пробоя гермооболочки пилотируемого модуля высокоскоростной частицей. Определение координат происходит путем анализа времени прихода акустического импульса (бегущей звуковой волны, имеющей характер резкого кратковременного изменения давления) на датчики (малогабаритные микрофоны), расположенные по определенной схеме [1, 2].
Импульсные сигналы, наблюдающиеся в воздушной среде пилотируемого модуля, могут быть вызваны работой различных бортовых систем и аппаратуры, жизнедеятельностью экипажа станции и, возможно, реальными ударами микрочастиц.
Целью данной работы является идентификация источников импульсных сигналов, т. е. выделение среди всего множества сигналов групп, соответствующих одному источнику. Для этого необходимо решить задачу нахождения схожих импульсов.
Методика расчета
Для выявления схожести импульсных сигналов используется взаимный спектрально-корреляционный анализ.
В качестве критерия схожести импульсов принята величина коэффициента корреляции Пирсона, большая некоторого порогового значения г , как это применяется при построении корреляционных плеяд Терен-тьева и вроцлавской таксономии [3]. Основным критерием выполнения задачи является выделение в отдельную группу импульсов от срабатывания имитатора пробоя, входящего в состав научной аппаратуры «Средства оперативного определения координат точки пробоя» и предназначенного для создания импульсной акустической волны, близкой по параметрам к волне, возникающей в воздушной среде пилотируемого модуля при пробое гермооболочки высокоскоростной частицей [4].
Пусть имеются K спектров импульсов {а} от некоторого количества источников, где г = 1, 2, ..., К. Каждый спектр представляет
собой массив уровней звукового давления в Т спектральных полосах ||ай||, где Ь = 1, 2, ..., Т.
Спектры сигналов, относящихся к одному и тому же источнику, обладают попарно высоким линейным коэффициентом корреляции
-а- а)(а*
■ а)
т.. =
ч
\1ц= (аи - а-) 2 1 ( а
а)2
где I, . = 1, 2, ..., К. Верхняя черта означает усреднение по уровням спектра.
Корреляционная матрица ||т..|| является симметричной, а все ее диагональные элементы равны 1.
При рассмотрении коэффициентов корреляции спектров сигналов важны только близкие к 1 значения. Близкие к -1 значения линейных коэффициентов корреляции ничего не говорят о взаимосвязи спектров сигналов, в отличие от близких к -1 линейных коэффициентов корреляции самих сигналов.
Выбрав пороговый уровень -1 < гп < 1, мы можем одновременной перенумерацией строк и столбцов привести матрицу ||т..|| к виду
1
т.,
г1
т..
У
К
т
'к
1
где в незаштрихованных областях каждый коэффициент т.. < тп, а вдоль главной диагонали располагаются одна или несколько квадратных симметричных подматриц, таких, что каждая либо состоит из единственного элемента, либо для любых двух столбцов . и & (^ &) существует такая строка г, что т.. > т , т., > т .
г. п ш п
Существование такого преобразования очевидно в предельных случаях: при т = -1 в общем случае вдоль диагонали будет располагаться одна подматрица, совпадающая с матрицей ||г.||, а при тп = 1 вдоль диагонали будут располагаться К матриц ||1||. Очевидно, что с ростом тп количество подматриц, располагаемых вдоль диагонали, не может уменьшаться. Соответственно, при -1 < тп < 1 количество таких подматриц будет лежать на отрезке [1; К].
В ходе эксперимента на борту служебного модуля (СМ) РС МКС с декабря 2014 г. по ноябрь 2015 г. (экспедиции МКС-41/42, МКС-43/44 и МКС-45/46) блоком преобразования акустических сигналов (БПАС) было зафиксировано 1 112 700 акустических импульсов. Каждый импульс представляет собой запись сигналов с частотой опроса ^ = 50 кГц длительностью Ь = 80,192 мс
^ опроса ^ ^ порции 7
по каждому датчику (4 096 опросов).
Для корреляционного анализа были проведены исследования параметров представления спектров сигналов:
• октавные, 1/3-октавные или 1/12-ок-тавные спектры;
• представление в физических единицах или в логарифмических (дБ);
• диапазон анализируемых частот.
Было установлено, что наилучшие результаты дает представление 1/3-октавных спектров в децибелах.
На рис. 1 и 2 представлены узкополосный (в полосах 1 Гц) и 1/3-октавный спектры хлопков имитатора пробоя на фоне прочих импульсов, а на рис. 3 — развертка этого импульса по времени.
Частота,
Рис. 1. Узкополосный спектр акустического импульса от имитатора пробоя (красный) на фоне спектров прочих импульсов
Рис. 2. 1/3-октавный спектр акустического импульса от имитатора пробоя (красный) на фоне спектров прочих импульсов
т
1
1
1
1
1
т
367 300
200
100
б -Н)0
и
е -200
г
£-300
-408
[ 1
л Л' и', АДА ■Л' \Ь \ ЛлЛ А -Л л.
г 1 1 /Ц V у/ V V V V V > *\)
1 I ( I У 1
II
18,4330 18,4340 18,4350 18,4360 18,4370 18,4380 18,4390 18,4400 18,4410 18,4420
I, с
Рис. 3. Импульс от имитатора пробоя
Из рис. 2 видно, что в 1/3-октавных полосах с центральными частотами 315 Гц и ниже уровни звукового давления имитатора пробоя сравнимы с фоновыми. Поэтому для корреляционного анализа были выбраны 1/3-октавные спектры с центральными полосами частот от 400 Гц до 16 кГц. Наибольшее превышение спектров импульсов имитатора над спектрами фонового шума наблюдается на частоте 8 кГц. Далее рассмотрим те сигналы, для которых уровни звукового давления в 1/3-октавной полосе с центральной частотой 8 кГц превышали 67,5 дБ по всем каналам. Это гарантировало, что все сигналы имитатора, которых по протоколам проведения космического эксперимента «Пробой» было зарегистрировано 144, в рабочем отсеке СМ попадут в выборку. Всего в выборку попала 2 941 запись. Оставшиеся сигналы также обрабатывались, но ввиду больших объемов в данной статье их анализ не проводился.
Из рис. 3 видно, что длительность импульса от имитатора пробоя составляет не более t =10 мс. Для учета в анализе
импульса ^ ^
только спектра самого импульса, а не фоновых шумов или отраженных волн, возьмем от каждого сигнала порцию N = 512 опро-
^ 1 ^ порции 1
сов. Длительность порции определяется по формуле:
T =(1// Ш
порции опроса порции
и составляет 10,24 мс. Результаты
Построение корреляционной матрицы происходило итерационно, с постепенно растущим пороговым значением гп от -1 до 1 с шагом 0,001.
На рис. 4 изображена преобразованная корреляционная матрица, в которой точками
обозначены значения г.. > 0,970: цветными — запротоколированные сигналы, черными — все прочие сигналы.
Горизонтальными серыми линиями в верхней правой части матрицы показаны сноски на номера сигналов, приведенных в нескольких колонках справа от матрицы. Номера сигналам присвоены в порядке убывания их уровней. Так как корреляционная матрица симметрична, сноски на рис. 4 не мешают разглядеть структуру преобразованной матрицы.
Из рис. 4 видно, что большинство цветных сигналов сгруппировалось в одном квадрате (увеличенное изображение этого квадрата показано в левом нижнем углу рисунка). В него попали 75 выделенных цветом импульсов, соответствующих хлопкам имитатора пробоя, зарегистрированным в протоколах, а также 21 «черный» импульс, которые соответствуют хлопкам имитатора, не вошедшим в протоколы испытаний. Последнее обусловлено особенностями процедуры верификации системы с использованием имитатора пробоя и чувствительностью системы.
Чувствительность системы связана с установленным пороговым уровнем АЬ срабатывания системы по превышению уровня сигнала над уровнем фонового шума (по умолчанию, АЬ = 20 дБ). При установке слишком низкого порогового уровня чувствительность системы повышалась, и после срабатывания системы на звук имитатора пробоя оператор, работающий за лэптопом, не всегда успевал занести данные о координатах места пробоя и времени события, выданные системой в протокол эксперимента до следующего срабатывания системы. Наоборот, при установке слишком высокого порогового уровня система могла не среагировать на сигнал имитатора пробоя при его размещении в заглушенном запанельном пространстве.
Рис. 4. Преобразованная корреляционная матрица
Примечание. Точками на рисунке обозначены значения коэффициентов взаимной корреляции Пирсона ту, такие, что ту > 0,970, при этом различным цветам соответствуют запротоколированные сигналы имитатора пробоя в различных зонах СМ: В — в запанельном пространстве рабочего отсека; В — на открытых участках рабочего отсека; I — в переходном отсеке; В — в запанельном пространстве переходной камеры; — на открытых участках переходной камеры; В — все прочие сигналы. Горизонтальными серыми линиями показаны сноски на номера сигналов (правая колонка) в порядке убывания их уровней.
В то же время пороговый уровень АЬК для регистрации импульсных шумов на внутренний носитель системы был значительно ниже порогового уровня срабатывания системы (АЬК = 6 дБ), поэтому все пропущенные оператором сигналы были зарегистрированы и вошли в наш квадрат «без протокола».
При этом 21 запротоколированный сигнал имитатора не попал в вышеупомянутый квадрат ввиду того, что экипаж МКС-41/42 зачастую использовал нештатные средства
генерации акустических импульсов, спектральные характеристики которых существенно отличались от характеристик имитатора пробоя.
Непопадание оставшихся 27 из 144 запротоколированных событий в корреляционную матрицу вызвано срабатыванием системы БПАС, размещенной в рабочем отсеке СМ, не на все сигналы от имитатора пробоя при его установке в смежных отсеках СМ: переходной камере или переходном отсеке. Для регистрации и протоколирования
сигналов в этих отсеках использовался переносной автономный регистратор, обработка данных с которого в данной статье не рассматривается.
Заключение
Разработан и реализован в программном коде алгоритм распознавания источников импульсных акустических шумов, действующих на борту служебного модуля РС МКС.
Выполнены параметрические исследования разрешающей способности коэффициентов линейных корреляций для различных представлений спектров. Выполнен обоснованный выбор параметров спектрального и корреляционного анализов.
Сформирована корреляционная матрица 1/3-октавных спектров, которая содержит четко выделенные группы импульсных шумов, включая группу, соответствующую спектрам сигналов имитатора пробоя, который размещался более чем в 100 точках объема служебного модуля, в т. ч. и в запа-нельном пространстве.
Разработанный алгоритм планируется использовать при доработке программного обеспечения научной аппаратуры по выделению в реальном масштабе времени импульсных сигналов пробоя с передачей телеметрической информации на Землю в экспедиции МКС-53/54.
Результаты исследований показывают принципиальную возможность идентификации импульсных источников шума на борту РС МКС.
Список литературы
1. Половнев А.Л. Определение координат точки пробоя высокоскоростной частицей на борту служебного модуля МКС // Научный вестник МГТУ ГА. 2009. № 56. С. 198-203.
2. Патент 2387966 С1. Российская Федерация. Устройство определения координат места пробоя гермооболочки пилотируемого космического объекта и способ определения координат места пробоя. Авершьев С.П., Болотин В.А., Дементьев В.К., Дядькин А.А., Макаревич Г.А., Михайлов А.В., Пелипен-ко Л.Ф., Половнев А.Л., Рыбак С.П., Тихомиров Н.А.; заявители и патентообладатели — ОАО РКК «Энергия» и ФГУП ЦНИИмаш; заявка 2008140955/28; приоритет от 16.10.2008 г., опубликовано 27.04.2010 г.
3. Шитиков В.К., Розенберг Г.С., Зин-ченко Т.Д. Количественная гидроэкология: методы системной идентификации. Тольятти: ИЭВБ РАН, 2003. 463 с.
4. Авершьев С.П., Будаев В.С., Макаревич Г.А., Михайлов А.В., Пелипенко Л.Ф., Тихомиров Н.А., Дементьев В.К., Рыбак С.П., Половнев А.Л., Сычев А.В. Акустические волны в гермоотсеке космического аппарата при его пробое высокоскоростной частицей // Космонавтика и ракетостроение. 2011. № 1(62). С. 12-17.
Статья поступила в редакцию 18.08.2017 г.
Reference
1. Polovnev A.L. Opredelenie koordinat tochki proboya vysokoskorostnoi chastitsei na bortu sluzhebnogo modulya MKS [Determining coordinates of a leakage point punctured by a high-velocity particle onboard the ISS Service Module]. Nauchnyi vestnik MGTU GA, 2009, no. 56, pp. 198 -203.
2. Patent 2387966 S1. Rossiiskaya Federatsiya. Ustroistvo opredeleniya koordinat mesta proboya germoobolochki pilotiruemogo kosmicheskogo ob»ekta i sposob opredeleniya koordinat mesta proboya [A device for determining coordinates of a puncture in the pressure shell of a manned space vehicle and a method for determining coordinates of the puncture]. Aversh'ev S.P., Bolotin V.A., Dement'ev V.K., Dyad'kin A.A., Makarevich G.A., Mikhailov A.V., Pelipenko L.F., Polovnev A.L., Rybak S.P., Tikhomirov N.A.; the applicants and the patent owners — OAO RKK «Energiya» and FGUP TsNIImash; application 2008140955/28; priority of 16.10.2008, published 27.04.2010.
3. Shitikov V.K., Rozenberg G.S., Zinchenko T.D. Kolichestvennaya gidroekologiya: metody sistemnoi identifikatsii [Quantitative hydroecology: system identification methods]. Tol'yatti: Institute of Ecology of Volga basin RAS publ., 2003. 463 p.
4. Aversh'ev S.P., Budaev V.S., Makarevich G.A., Mikhailov A.V., Pelipenko L.F., Tikhomirov N.A., Dement'ev V.K., Rybak S.P., Polovnev A.L., Sychev A.V. Akusticheskie volny v germootseke kosmicheskogo apparata pri ego proboe vysokoskorostnoi chastitsei [Acoustic waves in the pressure shell of a spacecraft when punctured by a high-velocity particle]. Kosmonavtika i raketostroenie, 2011, no. 1(62), pp. 12-17.
