Разработка алгоритма идентификации - Программирование, компьютеры и кибернетика дипломная работа

Состав и принцип работы аппаратуры. Выбор параметров корреляционного анализа и Фурье-анализа. Разработка и применение алгоритма корреляционного анализа. Реализация алгоритма Фурье-анализа на языке С++ и алгоритма корреляционного анализа на языке С#.


посмотреть текст работы


скачать работу можно здесь


полная информация о работе


весь список подобных работ


Нужна помощь с учёбой? Наши эксперты готовы помочь!
Нажимая на кнопку, вы соглашаетесь с
политикой обработки персональных данных

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Определения, обозначения и сокращения
1. Научная аппаратура СОКП и исходные данные для анализа
1.1 Состав и принцип работы аппаратуры
1.2 Исходные научные данные для анализа
2.1 Выбор параметров корреляционного анализа
3. Алгоритм корреляционного анализа
3.1 Разработка алгоритма корреляционного анализа
3.2 Применение алгоритма корреляционного анализа
Приложение 1. Реализация алгоритма Фурье-анализа на языке С++
Приложение 2. Реализация алгоритма корреляционного анализа на языке С#
Определения, обозначения и сокращения
АСУ - ассенизационное санитарное устройство
БПАС - блок преобразования акустических сигналов
ВГК - внешний гидравлический контур
ИУС - информационно-управляющая система
МКС - Международная Космическая Станция
НПО ИТ - Научно-производственное объединение измерительной техники
РКК - Ракетно-космическая корпорация
СКВ - система кондиционирования воздуха
СОКП - средства оперативного определения координат точки пробоя
СОТР - система обеспечения теплового режима
СУБА - система управления бортовой аппаратурой
ЦПК - Центр подготовки космонавтов им. Ю.А. Гагарина
ФГУП ЦНИИмаш - Федеральное государственное унитарное предприятие Центральный научно-исследовательский институт машиностроения
В настоящее время всё более актуальной становится проблема возрастающей засоренности околоземного космического пространства космическим мусором. В связи с чем, вероятность столкновения космической станции с различными техногенными частицами также возрастает. Кроме этого, всегда существует вероятность столкновения с природными небесными телами - метеороидами, движущимися в межпланетном пространстве, и по размеру занимающими промежуточное положение между межпланетной пылью и астероидом.
Контроль за космическим мусором вокруг МКС идет постоянно, и при необходимости специалисты ЦУП «заставляют» МКС делать маневр «Уклонение от осколков».
В результате столкновения космического объекта с корпусом станции может произойти разгерметизация модулей станции, которая представляет собой реальную опасность, в результате которой необходимо мгновенное принятие решения, чтобы устранить течь.
Печальный опыт разгерметизации модулей станции в истории пилотируемой космонавтики был, когда произошла разгерметизация модуля «Спектр» космической станции «Мир». Тогда найти место разгерметизации модуля и восстановить его работу не удалось.
В связи с этим при проектировании новых пилотируемых космических модулей и эксплуатации уже существующих должна учитываться возможность пробоя гермооболочки в результате столкновения с техногенной частицей или микрометеороидом. Вероятность непробоя гермооболочки МКС даже при усиленной экранной защите модулей оценивается величиной на уровне 0,85 в течение 15 лет эксплуатации.
Н а Российском сегменте МКС в настоящее время имеются течеискатели, которые позволяют лишь при сканировании всей поверхности пилотируемого модуля определять точку пробоя. П олученный на станции «Мир» при разгерметизации модуля «Спектр» опыт , показал необходимость повышения эффективности способов обеспечения живучести станции при её разгерметизации. Особенно это касается оперативности определения места пробоя.
При образовании небольших отверстий (площадью не более нескольких единиц и десятков мм 2 ) существует техническая перспектива заделки отверстия силами экипажа. Основным условием при этом является оперативное определение координат места пробоя (удара), позволяющее свести к минимуму время для принятия решения о необходимых действиях экипажа.
В настоящее время на РС МКС проводятся космические эксперименты «Отклик» и «Пробой». В КЭ «Отклик» исследуется возможность определения места соударения гермооболочки космического модуля путём анализа волн вибраций, распространяющихся по самому гермокорпусу модуля, а в КЭ «Пробой» - посредством анализа акустических волн, распространяющихся в воздушной среде модуля.
Космический эксперимент «Пробой» проводится с целью верификации метода оперативного определения координат точки пробоя. Подзадачами КЭ «Пробой» являются изучения собственной акустической обстановки и идентификация постоянно или периодически действующих на борту РС МКС источников импульсного акустического шума. Данная идентификация необходима в связи с потребностью отсеивать известные сигналы работы бортовой аппаратуры от потенциально возможных сигналов реального пробоя. Изучение акустической среды в РС МКС необходимо для учёта всех параметров при проектировании штатной системы, максимизации точности и устойчивости решений.
Шумящее оборудование на борту РС МКС включает в себя насосы СПН и ГЖТ ВГК СОТР, клапаны БВК, элементы системы АСУ, циркуляционные и комфортные вентиляторы, компрессор и вентилятор СКВ , различная научная аппаратура и многое другое. Также дополнительными источниками акустических импульсов являются действия самих космонавтов: занятия на беговой дорожке, приготовление пищи, работа с приборами, открытие и интерьерных панелей и дверей кают, разговоры с Землёй и между собой. Спектры звуков защёлок интерьерных панелей были отдельно исследованы на УТМ в ЦПК им. А.Ю. Гагарина.
Распознавание сигналов широко используется в различных сферах человеческой деятельности: в сейсмологии [3, 4], океанологии [5], кардиологии, сфере IT [6, 7, 8] и многих других.
При цифровом распознавании сигналов используются различные методы: сингулярные разложения [9], метод главных компонент [10], спектральный анализ и анализ фурье-методами [11, 12], корреляционные методы [12], вейвлет-преобразования [13, 14].
Целью данной работы является идентификация источников импульсного акустического шума, действующих на борту РС МКС, для дальнейшей доработки научной аппаратуры «СОКП» и создания штатной системы.
Задачей работы является разработка алгоритма, позволяющего эффективно различать импульсные шумы от различной бортовой аппаратуры, его реализация в программном коде и применение для анализа научной информации, полученной в ходе КЭ «Пробой» .
Объектом данного исследования являются действующие на борту источ ники импульсных шумов, а предметом исследования - алгоритм их идентификации.
В первом разделе работы описан научная аппаратура «СОКП», её устройство и принцип работы, а также формируемые ею в ходе КЭ «Пробой» данные для последующего анализа.
Во втором разделе проанализированы критерии, по которым следует идентифицировать источники импульсного акустического шум а, произведён выбор параметров анализа.
В третьем разделе описан алгоритм, который позволяет разделять импульсные сигналы на группы, различного происхождения, и приведены результаты его применения к научной информации, полученной в ходе КЭ «Пробой» в период с декабря 2014 года по ноябрь 2015 года.
В приложениях приведены листинги реализаций двух составных частей алгоритма в программном коде.
1 . Научная аппаратура СОКП и исходные данные для анализа
1.1 Состав и принцип работы аппаратуры
Подзадачами космического эксперимент а «Пробой» являются получение характеристик постоянно действующих помех в реальных условиях полета для настройки пороговых уровней срабатывания системы оперативного определения координат точки пробоя, идентификация источников импульсных акустических сигналов , действующих на борту СМ, определение характеристик воздушной среды в СМ, влияющих на скорость распространения звуковых волн и верификация метода оперативного определения координат точки пробоя с помощью имитатора звуковой волны пробоя в натурных условиях полета.
Результаты космического эксперимента планируется использовать при создании штатной системы оперативного определения координат места пробоя.
Принцип действия системы оперативного определения координат пробоя основан на измерении разновременности прихода акустической волны от точки удара к акустическим датчикам (микрофонам), расположенным в различных точках гермоотсека по определенной схеме (на способ и схему размещения датчиков имеются патенты [1, 2]), и расчете координат источника акустической волны на основании этих данных.
Возможность регистрации акустических волн пробоя и определения координат пробоя подтверждена в наземных испытаниях, в которых исследовались акустические волны при натурном пробое элементов конструкции КА для скоростей соударения до 7 км/с.
Важнейшим элементом аппаратуры КЭ является имитатор акустической волны пробоя. Характеристики акустической волны, получаемой в имитаторе, близки к характеристикам акустической волны пробоя, что подтверждено в наземных экспериментах по натурному пробою элементов конструкции КА. В качестве источника импульсного звука в имитаторе использован проц есс разрыва латексной мембраны.
Значительное место в наземной отработке заняло также решение проблемы электромагнитной совместимости научной аппаратуры космического эксперимента с бортовыми техническими средствами. Проблема была успешно решена, причем на основе элементов отечественной электронной компонентной базы.
Метод определения координат пробоя был успешно подтвержден в экспериментах, проведенных на КС СМ в РКК «Энергия» и на УТМ СМ в ЦПК им . Ю.А. Гагарина.
Постановщиком эксперимента «Пробой» и разработчиком научной аппаратуры является ФГУП ЦНИИмаш. В разработке и изготовлении имитатора пробоя принимало участие ОАО «Р КК «Энергия», а в изготовлении вторичных источников питания - ОАО «НПО ИТ».
Научная аппаратура КЭ «Пробой» состоит из двух независимых регистрирующих блоков: БПАС ( рисунок 1 ), АР ( рисунок 2 ) с вторичными источниками питания, комплекта микрофонов в адаптерах, имитатора пробоя (рисунок 3) .
БПАС размещается в РО СМ, АР используется как переносной регистратор и размещается поочередно в ПхО СМ или ПрК СМ.
К каждому из регистрирующих блоков подключено по 6 датчиков . В качестве чувствительных элементов используются фазированные ICP -микрофоны массива типа 4958 фирмы “ Bruel & Kjaer ”. Микрофоны закреплены в специальных адаптерах, которые защищают их от случайных повреждений. Адаптеры крепятся на панелях интерьера или вблизи гермооболочки с помощью велкро (рисунок 4) .
Рис. 4. Микрофон в сборе закреплен на панели интерьера
Кабели от микрофонов в РО прокладываются на стыках панелей интерьера и закрепляются с помощью липкой серой ленты, при этом предусматривается возможность свободного открытия панелей. Кабели в ПхО и ПрК прокладываются вдоль гермооболочки.
На рисунке 5 представлена структурная схема НА «СОКП» взаимодействия с системами СМ.
Размещено на http://www.allbest.ru/
Рис. 5. Структурная схема НА «СОКП»
Электрическая схема подключения НА «СОКП» приведена на рисунке 6.
Рис. 6. Схема подключения НА «СОКП»
Основные технические характеристики НА «СОКП» приведены в таблице 1.
Технические характеристики НА «СОКП»
Полоса частот регистрируемых акустических сигналов, Гц
Динамический диапазон измерений, дБ
Уровень импульсного звука на расстоянии 3 м от имитатора пробоя, дБ (при длительности импульса не более 70 мс)
Пробой гермооболочки пилотируемого космического модуля высокоскоростной частицей сопровождается импульсным энерговыделением, приводящим к образованию в атмосфере гермоотсека ударной волны, которая по мере распространения внутри модуля вырождается в звуковую.
Принцип работы системы оперативного определения координат точки пробоя корпуса космического аппарата высокоскоростной частицей основан на измерении сдвигов по времени моментов прихода акустической волны, возникающей в точке пробоя, к измерительным датчикам и последующем определении координат точки пересечения шаровых фронтов обращенных виртуальных воздушных волн, генерируемых в точках расположения датчиков, с учётом измеренных временных сдвигов.
Исходными данными для решения обратной задачи нахождения координат источника звука являются координаты микрофонов, моменты времени прихода волны к микрофонам и скорость звука.
Основные затруднения при определении момента прихода переднего фронта волны от гермооболочки связаны с прохождением звука в рабочее пространство модуля через приборный отсек, отделенный от рабочей зоны обитания космонавтов панелями интерьера толщиной от 8 до 20 мм, изготовленными из вспененного полистирола, обшитого тонким слоем алюминиевой фольги толщиной ~ 0.3 мм и оклеенными ворсовой тканью типа «велкро». Наземные эксперименты показали, что при прохождении акустической волны через панель интерьера ее амплитуда снижается на 30-40 дБ.
Учитывая, что уровень звука, генерируемого при пробое, составляет не менее 140 дБ на расстоянии 1 м от точки пробоя, а уровень фоновых шумов в рабочем отсеке составляет 60-70 дБА, можно достаточно надежно фиксировать акустическую волну пробоя после ее прохо ждения сквозь панели интерьера.
Предполагая, что звук от источника распространяется к микрофонам прямолинейно с постоянной скоростью, можно записать систему уравнений, решение которой позволяет определить координаты источника акустической волны.
Реальная погрешность определения момента прихода волны к датчику увеличивается за счет спектральных характеристик импульсных шумов, дифракции волны на элементах оборудования и т.д.
В космическом эксперименте для выявления переднего фронта импульсной звуковой волны на фоне шума оборудования реализованы различные алгоритмы обработки сигналов.
При срабатывании системы сигнал записывается регистратором в файл, который оперативно передается на бортовой лэптоп по беспроводному интерфейсу WiFi. При этом на бортовом лэптопе запускается вычислительная часть ПО, которая обрабатывает полученный файл, находит моменты времени прихода переднего фронта волны к разным датчикам, вычисляет координаты точки пробоя и выводит на экран математические координаты и проекцию этой точки на гермооболочку модуля. При определении координат точки пробоя используется проецирование расчетной точки пробоя на гермооболочку.
В сеансах верификации, периодически проводимых с участием экипажа, в заранее выбранную точку на гермооболочке устанавливается излучатель имитатора пробоя. Скорость звука в воздушной среде уточняется по измерениям те мпературы внутренней атмосферы.
При срабатывании системы на экран бортового лэптопа выдается изображение проекции модуля, на которой мигающей точкой п оказано расчетное место пробоя.
Излучатель имитатора пробоя достаточно компактен, что позволяет размещать его вблизи гермооболочки. При этом имеется возможность достаточно плотно закрывать панель интерьера, что позволяет наиболее полно моделировать условия распространения звуковой волны при реальном пробое.
Кроме коротких сеансов верификации, в которых определяются координаты источника акустической волны, имитирующего пробой, в космическом эксперименте также проводится непрерывный сбор сигналов фоновых шумов в звуковом диапазоне при превышении амплитуды сигнала фонового уровня шумов на определенную величину.
В длительном сеансе измерений регистрирующие блоки включаются для непрерывного сбора информации об импульсных и постоянных шумах в течение месяца.
Все результаты коротких и длительных сеансов верификации записы ваются в файлы на флэш-накопители, которые возвращаются на Землю с экипажем.
1.2 Исходные научные данные для анализа
Акустический импульс (от лат. impulsus - удар, толчок) - это бегущая звуковая волна, имеющая характер резкого кратковременного изменения давления. Такие звуковые волны создаются, например, взрывом, искровым разрядом, соударением тел. Каждый такой импульс содержит как область повышенного, так и область пониженного давления. Спектр такого акустического импульса сплошной, с максим умом в области частот, период которых близок к длительности импульса .
В 2008 году проводил и сь исследования пробоя фрагмента гермооболочки модуля на установке ЦНИИмаш с помощью газокумулятивной пушки, в стволе которой снаряд из стали или алюминия массой до 0,2 г разгонялся до 5 км/c. На рисунках 7 и 8 показан импульс в разных масштабах времени, а на рисунке 9 - спектр данного импульса.
Рис. 7. Импульс пробоя гермооболочки высокоскоростным снарядом
Рис. 8. Развёртка импульса пробоя гермооболочки
Рис. 9. 1/3-октавный спектр импульс пробоя гермооболочки
Из рисунка 8 видно, что длительность основного импульса составляет порядка t импульса = 0,01 с.
Метод определения координат точки пробоя (равно как и любого другого акустического импульса) основан на конечности скорости звука и определении разницы во времени прихода фронта акустической волны на разные датчики.
Первый этап КЭ «Пробой» проводился на борту СМ РС МКС в трёх экипажами трёх миссий:
- миссия 41 -4 2 с декабря 2014 по февраль 2015;
- миссия 43 - 44 с мая 2015 по июнь 2015;
- миссия 45-46 с октября 2015 по ноябрь 2015.
В каждой миссии проводилось несколько коротких сеансов верификации и длительные сеансы непрерывного включения аппаратуры СОКП. Всего суммарная длительность длительных сеансов составила 91,66 суток.
В ходе первого этапа КЭ «Пробой» в РО СМ РС МКС непрерывно находилось 6 микрофонов, подключенных к БПАС . Ещё 6 микрофонов, подключенных к АР попеременно располагались в ПхО и ПрК. Микрофоны располагались согласно запатентованной схеме [1, 2], а также указаниям специалистов по бортовой аппаратуре для обеспечения комфортного взаимодействия с приборами и системами служебного модуля. В миссии 45-46 эксперимент проводился с проводной синхронизацией, для чего на 7 и 8 каналы БПАС заводились данные с АР. Координаты микрофонов приведены в таблице 2 и на рисунк ах 10 и 11 . На рисунках 10 и 11 также представлен внешний вид программы [3].
Рис. 10. Координаты датчиков (вид СМ сверху) и внешний вид программы.
Рис. 11. Координаты микрофонов (вид на правый борт СМ сверху).
Данные с микрофонов снимались с частотой опроса f опроса = 50 кГц синхронно по всем каналам. Из этих данных формировались буферы по N буфера = 1024 показания по каждому каналу. Значит, каждый буфер соответствует сигналу длительностью t буфера = 0,02048 с. Как видно, t буфера примерно вдвое превышает длительность t импульса .
В каждом буфере вычислялся суммарный уровень звукового давления. Эти суммарные уровни записывались в отдельные файлы. Циклограмма типовых рабочих суток экипажа МКС (12 февраля 2015 года) по датчику М4 представлена на рисунке 12.
Рис. 12. Зависимость уровня звука (в дБА) от времени суток (в ч)
Как видно из рисунка, максимального уровня в 110 дБА шум достигает в часы работы экипажа. В ночное же время над уровнем фона превышают в основном только импульсы от срабатывания клапанов БВК каждые 15 минут.
Если суммарный уровень превышал уровень фоновых шумов минимум на 6 дБ хотя бы по двум датчикам, этот буфер, а также один предыдущий и два последующих записывались в файл импульса. Кроме самого сигнала из 4Ч N буфера = 4096 точек в файлы записывались дата время регистрации этого импульса с точностью до 1 мс. На рисунке 13 приведен вид программы LoadSignal [12] с загруженной записью импульса.
Рис. 13. Порция сигнала в программе LoadSignal
С 15:13 13 ноября 2015 года по 03:56 14 ноября 2015 года производилась регистрация импульсов по превышению на 6 дБ не по двум, а по одному каналу. Это привело к аномальному количеству записей импульсов: за t 2 = 91,66 суток длительных сеансов было зарегистрировано m 2 = 387 000 импульсов, в то время как за t 1 = 0,53 суток 13-14 ноября 2015 зарегистрировано m 1 = 725 700 импульсов. Индексы в данном случае обозначают количество датчиков, по которым происходила регистрация импульса.
Аномальность числа импульсов m 1 выражается в том, что за время t 1 суммарное время всех записанных импульсы составляет
Подставив значения, получим T 1 = 0,69 суток , что превышает t 1 . Это можно объяснить тем, что фактом, что дискретность регистрации импульсов равна длительности одного t буфера , а размер записи превосходит в 4 раза это время. Таким образом, может происходить перекрытие соседних записей.
Для углубления знаний об акустической обстановке на борту РС МКС эти избыточные данные после удаления зон перекрытий были соединены в большие wav -файлы. Таким образом, удалось собственным слухом почувствовать все события, произошедшие на борту за тот период. В частности, ночью удалось зафиксировать периодические изменения в акустической обстановке, вызываемые попеременной работой одного или двух роторов, обеспечивающих циркуляцию воздуха. Циклограммы ночных записей в 1/3-октавных спектральных полосах представлены на рисунке 14.
Рис. 14. Циклограммы ночных часов в 1/3-октавных полосах частот
Из рисунка видно, что на некоторых частотах между П-образными подъёмами, соответствующими срабатываниям клапанов БВК, уровни звукового давления в некоторых 1/3-октавных полосах частот заметно меняются. Особенно это касается 1/3-октавной полосы частот с центральной частотой 3150 Гц, в которой изменение уровня звукового давления достигает 7 дБ.
Для более подробного анализа были построены узкополосные спектры с разрешением 1 Гц двух последовательных периодов работы роторов, и произведено вычитание одного из другого. Таким образом, удалось исключить общие фоновые шумы из спектров каждого интервала и более наглядно продемонстрировать, что же именно изменилось. Результат приведён на рисунке 15.
Рис. 15. Разность спектров последовательных интервалов работы роторов БВК
Верхняя (выше 0) часть графика соответствует работе одиночного ротора, а нижняя (ниже 0) - работе пары роторов. Легко заметить, что узкие (тональные) пики, присутствующие в обоих случаях, меняют своё положение и амплитуду, в зависимости от того, сколько роторов работает одновременно. Частоты этих пиков являются кратными к собственной частоте вращения ротора. По графику можно определить, что в случае работы одного ротора частота вращения составляет 110,0 Гц, а при работе двух - 107,7 Гц. Эту разницу можно объяснить тем, что двум роторам, работающим совместно, не нужно так же много усилий, для обеспечения необходимых расходо-напорных характеристик, как одному.
Также из графика виден высокий резонанс на частоте 3300 Гц, который объясняется собственной частотой воздуховода. Не смотря на то, что суммарный уровень при работе двух роторов несколько выше, чем при работе одного, наличие узкого тонального пика делает звук одного ротора более раздражающим для человеческого уха, в чём в очередной раз удалось убедиться благодаря результатам эксперимента.
Таким образом мы показали, что НА СОКП помимо своей основной задачи оперативного определения координат точки пробоя может выполнять функции мониторинга, показывая изменение в режимах работы шумящей аппаратуры, связанные, возможно, в том числе и с её неисправностью, а также поиске резонансов, соответствующих собственным частотам акусти ческих объёмов, чтоб использовать эти данные при проектировании шумящей аппаратуры и элементов звукопоглащения.
Скорость звука в течение эксперимента могла меняться в зависимости от параметров окружающей среды (температуры и влажности, которые определялись согласно показаниям термогидрометра «ИВА-6» непосредственно экипажем станции и заносились в электронные протоколы эксперимента) в пределах c = 340±15 м/с. Учитывая максимальное расстояние между микрофонами БПАС М1 и М5 ? L = 5,258 м, максимальная разница во времени прихода фронта акустического импульса определяется по формуле
Таким образом, независимо от места расположения источника акустического импульса разница во времени прихода его волнового фронта на все датчики составляет не более длительности одного буфера t буфера . Это означает, что мы сможем зафиксировать момент прихода импульса по всем датчиками в рамках одной четырёхбуферной порции.
За время проведения длительных сеансов средняя плотность регистрации импульсов с
и составляет около 4 2 2 0 импульсов в сутки или 0,05 импульсов / с. Положим, что источники каждого импульса независимы (хотя это не совсем так, некоторые приборы имеют характерную циклограмму). В таком случае вероятность наложения двух импульсов определяется по формуле:
и равна 2 , 5 •10 - 7 . Учитывая, что всего было зарегистрировано m 2 = 387 000 импульсов, из них двойными были в среднем
что составляет 0, 1 . Таким образом, приходим к заключению, что количество наложенных импульсов практически равно нулю.
Тем самым мы приходим к следующим выводам:
1) Все импульсы, вызвавшие превышения по двум датчикам - одиночные, а значит, нет необходимости искать методы их разделения.
2) Записей сигналов импульсов в 4096 точек более чем достаточно для того, чтоб зафиксировать импульс по всем датчикам.
2.1 Выбор параметров корреляционного анализа
Идентифицировать источник импульсного сигнала можно несколькими способами: по суммарному уровню импульса, по спектру импульса и по координатам источника импульса. Для лучшего результата будем использовать комбинацию всех этих методов.
Для сопоставления двух и более сигналов по спектру и уровню проще и нагляднее сравнивать не сами сигналы, а их спектры. Выделим в каждом сигнале спектральные полосы частот, например 1/3 октавные. Пусть есть спектры сигналов и , где i = 1, 2, …, I ( I - количество спектральных полос частот) . Классически коэффициент линейной корреляции определяется по формуле:
Такой коэффициент лежит в диапазоне [-1; 1]. Случай, когда с = 1, соответствует абсолютно коррелирующим си гналам. П ри с = -1 сигналы анти коррелирующие. Эта формула совпадает с выражением для косинуса угла между двумя векторами и в многомерном ортонормированном пространстве.
Именно такая форма коэффициента линейной корреляции использовалась в работе [15]. Однако она имеет существенны е недостат ки. Во-первых, она учитывает только разницу в спектрах, а не в суммарных уровнях, так как на значение косинуса угла между векторами не оказывает влияние их длина. То есть нельзя различить источники с похожим спектром и разными суммарными уровнями сигнала. Во вторых, вычитание математических ожиданий и может сделать спектры похожими, даже если они таковыми не являлись.
Рассмотрим в таком пространстве радиус-вектора , где k = 1, 2, …, K ( K - количество векторов, k - порядковый номер вектора). Концы векторов образуют некоторое множество (облако) точек в этом пространстве. Причём вектора, соответствующие спектрам сигналов от одного источника, будут указывать на близкорасположенные точки. Действительно, если источники отличаются только по уровню сигнала, но схожи по спектру, тогда радиус-вектора их импульсов будут направленны в одну сторону, но различаться по длине, и чем больше будет различие по уровню, тем дальше будут разнесены концы радиус-векторов. Если же источники близки по суммарному уровню, но различаются по спектру, тогда радиус-векторы их будут направленны в разные стороны, и чем больше будет различие по спектрам, тем дальше будут разнесены концы радиус-векторов.
Теперь передвинем начало координат в центр облака с координатами
Далее будем работать с полученным множеством радиус-векторов . Такое преобразование даст б у льшую разрешающую способность при определении косинуса угла, так как на всё том же множестве (облаке) вершин радиус-векторов будут заведомо присутствовать те, которые отличаются на углы, близкие к р. То есть, если вы стоите в центре рощи, то в какую сторону бы вы не смотрели, за вашей спиной всегда будут деревья.
Запишем выражение для коэффициентов линейной корреляции - косинусов угла между этими новыми векторами:
Значения уровней звукового давления в полосах частот может быть представлены в линейных ( физических ) величинах (паскалях) или в логарифмических (децибелах). Сравним корреляционные матрицы, составленные из коэффициентов линейной корреляции, для множеств векторов , , в линейных и логарифмических величинах для каждого случая. Итого получим шесть корреляционных матриц. Ввиду всех изложенных выше выкладок можно предположить, что наилучший результат даст корреляция на множестве , а наихудший - на множестве . Кроме того, корреляция величин в логарифмическом масштабе должна дать лучший результат в сравнении с корреляцией в линейном масштабе, так как и датчики, и источники звука имеют мультипликативную, а не аддитивную природу.
Теперь выберем множество, на котором будем сравнивать корреляционные матрицы. Возьмём в качестве первого и второго сигналов спектр ы от имитатора пробоя снятые в реальных условиях космического эксперимента по двум разным каналам. В качестве векторов возьмём пять спектров различных сигналов в 1/3-октавных полосах частот с центральными частотами от 400 Гц до 16 кГц.
В качестве третьего сигнала возьмём второй , преобразованный так, чтоб сохранились его форма спектра и суммарный уровень, но значительно уменьшилось среднеквадратичное отклонение от среднего уровня.
В качестве четвёртого спектра возьмём так называемый «розовый шум» , спектральная плотность которого обратно пропорциональна частоте. Этот выбор обусловлен тем, что в шуме бортовой аппаратуры СМ РС МКС присутствуют шумы, близкие к «розовому» . На самом деле, чистого «розового шума» не наблюдается, т ак как «розовый шум» - это лишь математическая модель, к тому же обраща ющая знаменатель коэффициента линейной корреляции для . Чтоб избежать этого, мы введём в него малое возмущение - добавим 0,01 дБ в 1/3-октавной полосе частот с центральной частотой 400 Гц.
Ну и, наконец, в качестве пятого сигнала возьмём копию первого , сдвинутую на 20 дБ вниз .
Спектры всех пяти сигналов представлены на рисунке 1 6 .
Рис. 1 6. Спектры тестовых сигналов
В таблице 3 приведены спектры сигналов в числовом виде, а также корреляционные матрицы для каждого из 6 случаев. Корреляционные матрицы для удобства раскрашены в цветовые градации: чем больше коэффициент линейной корреляции (ближе к 1) , тем фон зеленее, чем он меньше (ближе к -1), том фон краснее.
Из таблицы 3 видно, что методы вычисления коэффициента линейной корреляции дают высокие значения (близкие к 1) для корреляции сигналов 1 и 2, что и следовало ожидать.
Корреляции для логарифмического масштаба на множестве не удалось разрешить ни одну группу сигналов, что, очевидно, не удовлетворяет нас. Корреляция для логарифмического масштаба на множестве причислила сигнал 3 к группе с сигналами 1, 2, а не к группе с сигналом 4, что не удовлетворяет условиям различия источников п
Разработка алгоритма идентификации дипломная работа. Программирование, компьютеры и кибернетика.
Реферат: Целеполагание в системе мотивации. Скачать бесплатно и без регистрации
Дезидерий Эразм Роттердамский
Курсовая работа по теме Розробка основних техніко-економічних показників виробничо-господарської діяльності підприємства
Реферат: Кабеля для компьютерной сети
Реферат: Высокочастотная связь по ВЛЭП. Скачать бесплатно и без регистрации
Реферат по теме Способы борьбы с наркотической зависимостью
Контрольные Работы Рейнбоу 5 Класс
Реферат: Adobe Photoshop. Монтаж. Скачать бесплатно и без регистрации
Курсовая работа: Инновационный менеджмент
Дипломная работа по теме Преодоление феномена 'выученной беспомощности' в процессе обучения младших подростков иностранному языку
Дипломная Работа На Тему Проектирование Системы Организации Труда Персонала
Реферат На Тему Научно-Технический Прогресс
Информационное Обеспечение Оперативно Розыскной Деятельности Контрольная Работа
Дипломная работа по теме Разработка сайта колледжа
Курсовая работа по теме Статистические ряды распределения в изучении структуры рынка
Доклад: Творчество Николая Клюева. Скачать бесплатно и без регистрации
Многие Истины Были Сначала Кощунством Эссе
Курсовые Работы По Фармакологии Готовые
Решебник По Русской Литературе 5 Класс Сочинения
Отчет по практике по теме Организация деятельности дошкольного образовательного учреждения
Анализ использования трудовых ресурсов на предприятии - Бухгалтерский учет и аудит курсовая работа
Уреаплазмоз: симптомы и последствия - Медицина презентация
The place of “Macbeth” among Shakespeare’s tragedies - Иностранные языки и языкознание курсовая работа