МУЗЫКА (алгоритм)

редактировать

Радиопеленгования с помощью алгоритма MUSIC

MUSIC ( MUltiple SIgnal Classification) - это алгоритм, используемый для оценки частоты и радиопеленгации.

СОДЕРЖАНИЕ

1 История
2 Теория
3 Сравнение с другими методами
4 Другие приложения
5 См. Также
6 Ссылки
7 Дальнейшее чтение

История

Во многих практических задачах обработки сигналов цель состоит в том, чтобы оценить на основе измерений набор постоянных параметров, от которых зависят принимаемые сигналы. Было несколько подходов к таким проблемам, включая так называемый метод максимального правдоподобия (ML) Капона (1969) и метод максимальной энтропии Бурга (ME). Хотя эти методы часто успешны и широко используются, они имеют определенные фундаментальные ограничения (особенно смещение и чувствительность в оценках параметров), в основном потому, что они используют неверную модель (например, AR, а не специальную ARMA ) измерений.

Писаренко (1973) был одним из первых, кто использовал структуру модели данных в контексте оценки параметров сложных синусоид в аддитивном шуме с использованием ковариационного подхода. Шмидт (1977), работая в Northrop Grumman и независимо Бьенвену и Копп (1979), был первым, кто правильно использовал модель измерения в случае массивов датчиков произвольной формы. Шмидт, в частности, добился этого, сначала получив полное геометрическое решение в отсутствие шума, а затем ловко расширив геометрические концепции, чтобы получить разумное приближенное решение в присутствии шума. Полученный алгоритм получил название MUSIC (MUltiple SIgnal Classification) и широко изучался.

В результате детальной оценки, основанной на тысячах симуляций, лаборатория Линкольна Массачусетского технологического института в 1998 году пришла к выводу, что среди принятых в настоящее время алгоритмов высокого разрешения MUSIC является наиболее многообещающим и ведущим кандидатом для дальнейшего изучения и фактической реализации оборудования. Однако, хотя преимущества в производительности MUSIC существенны, они достигаются за счет затрат на вычисления (поиск в пространстве параметров) и хранение (данных калибровки массива).

Теория

Метод MUSIC предполагает, что вектор сигнала, состоит из комплексных экспонент, частоты которых неизвестны, в присутствии гауссовского белого шума, как задано линейной моделью ${\ displaystyle \ mathbf {x}}$ $\ mathbf {x}$ ${\ displaystyle p}$ $п$ ${\ displaystyle \ omega}$ $\омега$ ${\ Displaystyle \ mathbf {п}}$ $\ mathbf {n}$

{\ displaystyle \ mathbf {x} = \ mathbf {A} \ mathbf {s} + \ mathbf {n}.}

{\ displaystyle \ mathbf {x} = \ mathbf {A} \ mathbf {s} + \ mathbf {n}.}

Вот это матрица Вандермонда из управляющих векторов и амплитуда вектора. Ключевым предположением является то, что количество источников, меньше, чем количество элементов в векторе измерения, т. Е. ${\ Displaystyle \ mathbf {A} = [\ mathbf {a} (\ omega _ {1}), \ cdots, \ mathbf {a} (\ omega _ {p})]}$ ${\ Displaystyle \ mathbf {A} = [\ mathbf {a} (\ omega _ {1}), \ cdots, \ mathbf {a} (\ omega _ {p})]}$ ${\ displaystyle M \ times p}$ ${\ displaystyle M \ times p}$ ${\ displaystyle \ mathbf {a} (\ omega) = [1, e ^ {j \ omega}, e ^ {j2 \ omega}, \ ldots, e ^ {j (M-1) \ omega}] ^ { T}}$ ${\ displaystyle \ mathbf {a} (\ omega) = [1, e ^ {j \ omega}, e ^ {j2 \ omega}, \ ldots, e ^ {j (M-1) \ omega}] ^ { T}}$ ${\ Displaystyle \ mathbf {s} = [s_ {1}, \ ldots, s_ {p}] ^ {T}}$ ${\ Displaystyle \ mathbf {s} = [s_ {1}, \ ldots, s_ {p}] ^ {T}}$ ${\ displaystyle p}$ $п$ ${\ displaystyle M}$ $M$ ${\ displaystyle p lt;M}$ ${\ displaystyle p lt;M}$

Автокорреляционная матрица затем дается ${\ displaystyle M \ times M}$ $M \ раз M$ ${\ displaystyle \ mathbf {x}}$ $\ mathbf {x}$

{\ displaystyle \ mathbf {R} _ {x} = \ mathbf {A} \ mathbf {R} _ {s} \ mathbf {A} ^ {H} + \ sigma ^ {2} \ mathbf {I},}

{\ displaystyle \ mathbf {R} _ {x} = \ mathbf {A} \ mathbf {R} _ {s} \ mathbf {A} ^ {H} + \ sigma ^ {2} \ mathbf {I},}

где - дисперсия шума, - единичная матрица и - матрица автокорреляции. ${\ displaystyle \ sigma ^ {2}}$ $\ sigma ^ {2}$ ${\ displaystyle \ mathbf {I}}$ $\ mathbf {I}$ ${\ displaystyle M \ times M}$ $M \ раз M$ ${\ displaystyle \ mathbf {R} _ {s}}$ ${\ displaystyle \ mathbf {R} _ {s}}$ ${\ displaystyle p \ times p}$ $р \ раз р$ ${\ displaystyle \ mathbf {s}}$ $\ mathbf {s}$

Матрица автокорреляции традиционно оценивается с использованием выборочной корреляционной матрицы. ${\ displaystyle \ mathbf {R} _ {x}}$ $\ mathbf {R} _ {x}$

{\ displaystyle {\ widehat {\ mathbf {R}}} _ {x} = {\ frac {1} {N}} \ mathbf {X} \ mathbf {X} ^ {H}}

{\ displaystyle {\ widehat {\ mathbf {R}}} _ {x} = {\ frac {1} {N}} \ mathbf {X} \ mathbf {X} ^ {H}}

где - количество векторных наблюдений и. Учитывая оценку, MUSIC оценивает частотный состав сигнала или автокорреляционной матрицы, используя метод собственного подпространства. ${\ displaystyle Ngt; M}$ ${\ displaystyle Ngt; M}$ ${\ Displaystyle \ mathbf {X} = [\ mathbf {x} _ {1}, \ mathbf {x} _ {2}, \ ldots, \ mathbf {x} _ {N}]}$ ${\ Displaystyle \ mathbf {X} = [\ mathbf {x} _ {1}, \ mathbf {x} _ {2}, \ ldots, \ mathbf {x} _ {N}]}$ ${\ displaystyle \ mathbf {R} _ {x}}$ $\ mathbf {R} _ {x}$

Поскольку матрица является эрмитовой, все ее собственные векторы ортогональны друг другу. Если собственные значения отсортированы в порядке убывания, собственные векторы, соответствующие наибольшим собственным значениям (то есть направлениям наибольшей изменчивости), охватывают подпространство сигнала. Остальные собственные векторы соответствуют собственному значению, равному шумовому подпространству, которое ортогонально сигнальному подпространству, и охватывают его. ${\ displaystyle \ mathbf {R} _ {x}}$ $\ mathbf {R} _ {x}$ ${\ displaystyle M}$ $M$ ${\ Displaystyle \ {\ mathbf {v} _ {1}, \ mathbf {v} _ {2}, \ ldots, \ mathbf {v} _ {M} \}}$ ${\ Displaystyle \ {\ mathbf {v} _ {1}, \ mathbf {v} _ {2}, \ ldots, \ mathbf {v} _ {M} \}}$ ${\ displaystyle \ mathbf {R} _ {x}}$ $\ mathbf {R} _ {x}$ ${\ Displaystyle \ {\ mathbf {v} _ {1}, \ ldots, \ mathbf {v} _ {p} \}}$ ${\ Displaystyle \ {\ mathbf {v} _ {1}, \ ldots, \ mathbf {v} _ {p} \}}$ ${\ displaystyle p}$ $п$ ${\ displaystyle {\ mathcal {U}} _ {S}}$ ${\ displaystyle {\ mathcal {U}} _ {S}}$ ${\ displaystyle Mp}$ $Mp$ ${\ displaystyle \ sigma ^ {2}}$ $\ sigma ^ {2}$ ${\ displaystyle {\ mathcal {U}} _ {N}}$ ${\ displaystyle {\ mathcal {U}} _ {N}}$ ${\ Displaystyle {\ mathcal {U}} _ {S} \ perp {\ mathcal {U}} _ {N}}$ ${\ Displaystyle {\ mathcal {U}} _ {S} \ perp {\ mathcal {U}} _ {N}}$

Обратите внимание, что для, MUSIC идентичен гармоническому разложению Писаренко. Общая идея метода MUSIC состоит в том, чтобы использовать все собственные векторы, охватывающие подпространство шума, для повышения производительности оценщика Писаренко. ${\ displaystyle M = p + 1}$ $М = р + 1$

Поскольку любой вектор сигнала, который находится в подпространстве сигнала, должен быть ортогонален подпространству шума, он должен быть таким для всех собственных векторов, которые охватывают подпространство шума. Чтобы измерить степень ортогональности по отношению ко всем, алгоритм MUSIC определяет квадрат нормы ${\ displaystyle \ mathbf {e}}$ $\ mathbf {e}$ ${\ displaystyle \ mathbf {e} \ in {\ mathcal {U}} _ {S}}$ ${\ displaystyle \ mathbf {e} \ in {\ mathcal {U}} _ {S}}$ ${\ Displaystyle \ mathbf {е} \ perp {\ mathcal {U}} _ {N}}$ ${\ Displaystyle \ mathbf {е} \ perp {\ mathcal {U}} _ {N}}$ ${\ Displaystyle \ mathbf {е} \ perp \ mathbf {v} _ {я}}$ ${\ Displaystyle \ mathbf {е} \ perp \ mathbf {v} _ {я}}$ ${\ Displaystyle \ {\ mathbf {v} _ {я} \} _ {я = p + 1} ^ {M}}$ ${\ Displaystyle \ {\ mathbf {v} _ {я} \} _ {я = p + 1} ^ {M}}$ ${\ displaystyle \ mathbf {e}}$ $\ mathbf {e}$ ${\ displaystyle \ mathbf {v} _ {i} \ in {\ mathcal {U}} _ {N}}$ ${\ displaystyle \ mathbf {v} _ {i} \ in {\ mathcal {U}} _ {N}}$

{\ displaystyle d ^ {2} = \ | \ mathbf {U} _ {N} ^ {H} \ mathbf {e} \ | ^ {2} = \ mathbf {e} ^ {H} \ mathbf {U} _ {N} \ mathbf {U} _ {N} ^ {H} \ mathbf {e} = \ sum _ {i = p + 1} ^ {M} | \ mathbf {e} ^ {H} \ mathbf { v} _ {i} | ^ {2}}

{\ displaystyle d ^ {2} = \ | \ mathbf {U} _ {N} ^ {H} \ mathbf {e} \ | ^ {2} = \ mathbf {e} ^ {H} \ mathbf {U} _ {N} \ mathbf {U} _ {N} ^ {H} \ mathbf {e} = \ sum _ {i = p + 1} ^ {M} | \ mathbf {e} ^ {H} \ mathbf { v} _ {i} | ^ {2}}

где матрица - это матрица собственных векторов, охватывающих шумовое подпространство. Если, то как следует из условия ортогональности. Использование выражения, обратного квадрату нормы, создает резкие пики на частотах сигнала. Функция оценки частоты для МУЗЫКИ (или псевдоспектра): ${\ Displaystyle \ mathbf {U} _ {N} = [\ mathbf {v} _ {p + 1}, \ ldots, \ mathbf {v} _ {M}]}$ ${\ Displaystyle \ mathbf {U} _ {N} = [\ mathbf {v} _ {p + 1}, \ ldots, \ mathbf {v} _ {M}]}$ ${\ displaystyle {\ mathcal {U}} _ {N}}$ ${\ displaystyle {\ mathcal {U}} _ {N}}$ ${\ displaystyle \ mathbf {e} \ in {\ mathcal {U}} _ {S}}$ ${\ displaystyle \ mathbf {e} \ in {\ mathcal {U}} _ {S}}$ ${\ displaystyle d ^ {2} = 0}$ $d ^ 2 = 0$

{\ displaystyle {\ hat {P}} _ {MU} (e ^ {j \ omega}) = {\ frac {1} {\ mathbf {e} ^ {H} \ mathbf {U} _ {N} \ mathbf {U} _ {N} ^ {H} \ mathbf {e}}} = {\ frac {1} {\ sum _ {i = p + 1} ^ {M} | \ mathbf {e} ^ {H } \ mathbf {v} _ {i} | ^ {2}}},}

{\ displaystyle {\ hat {P}} _ {MU} (e ^ {j \ omega}) = {\ frac {1} {\ mathbf {e} ^ {H} \ mathbf {U} _ {N} \ mathbf {U} _ {N} ^ {H} \ mathbf {e}}} = {\ frac {1} {\ sum _ {i = p + 1} ^ {M} | \ mathbf {e} ^ {H } \ mathbf {v} _ {i} | ^ {2}}},}

где - собственные векторы шума и ${\ Displaystyle \ mathbf {v} _ {я}}$ $\ mathbf {v} _ {i}$

{\ displaystyle \ mathbf {e} = {\ begin {bmatrix} 1 amp; e ^ {j \ omega} amp; e ^ {j2 \ omega} amp; \ cdots amp; e ^ {j (M-1) \ omega} \ end {bmatrix}} ^ {T}}

{\ displaystyle \ mathbf {e} = {\ begin {bmatrix} 1 amp; e ^ {j \ omega} amp; e ^ {j2 \ omega} amp; \ cdots amp; e ^ {j (M-1) \ omega} \ end {bmatrix}} ^ {T}}

- потенциальный вектор управления. Расположение наибольших пиков функции оценки дает оценки частоты для компонентов сигнала. ${\ displaystyle p}$ $п$ ${\ displaystyle p}$ $п$

{\ displaystyle {\ hat {\ omega}} = \ arg \ max _ {\ omega} \; {\ hat {P}} _ {MU} (e ^ {j \ omega}).}

{\ displaystyle {\ hat {\ omega}} = \ arg \ max _ {\ omega} \; {\ hat {P}} _ {MU} (e ^ {j \ omega}).}

МУЗЫКА - это обобщение метода Писаренко, которое сводится к методу Писаренко, когда. В методе Писаренко только один собственный вектор используется для формирования знаменателя функции оценки частоты; а собственный вектор интерпретируется как набор коэффициентов авторегрессии, нули которых могут быть найдены аналитически или с помощью алгоритмов поиска полиномиального корня. Напротив, MUSIC предполагает, что несколько таких функций были добавлены вместе, поэтому нули могут отсутствовать. Вместо этого есть локальные минимумы, которые могут быть обнаружены путем вычислительного поиска пиков в функции оценки. ${\ displaystyle M = p + 1}$ ${\ displaystyle M = p + 1}$

Сравнение с другими методами

MUSIC превосходит простые методы, такие как выбор пиков спектров DFT в присутствии шума, когда количество компонентов известно заранее, потому что он использует знание этого числа, чтобы игнорировать шум в своем окончательном отчете.

В отличие от DFT, он может оценивать частоты с точностью выше одной выборки, потому что его функция оценки может быть оценена для любой частоты, а не только для интервалов DFT. Это форма сверхразрешения.

Его главный недостаток состоит в том, что он требует, чтобы количество компонентов было известно заранее, поэтому исходный метод не может быть использован в более общих случаях. Существуют методы оценки количества исходных компонентов исключительно на основе статистических свойств автокорреляционной матрицы. См., Например, Кроме того, MUSIC предполагает, что сосуществующие источники некоррелированы, что ограничивает ее практическое применение.

Недавние итерационные полупараметрические методы предлагают надежное сверхразрешение, несмотря на сильно коррелированные источники, например SAMV.

Другие приложения

Модифицированная версия MUSIC, обозначенная как Time-Reversal MUSIC (TR-MUSIC), недавно была применена для вычислительной визуализации с обращением времени. Алгоритм MUSIC также был реализован для быстрого обнаружения частот DTMF ( двухтональная многочастотная передача сигналов ) в виде библиотеки C - libmusic.

Смотрите также

использованная литература

^ Хейс, Монсон Х., Статистическая обработка и моделирование цифровых сигналов, John Wiley amp; Sons, Inc., 1996. ISBN 0-471-59431-8.
^ Шмидт, Р. О., " Определение местоположения нескольких излучателей и параметров сигнала ", IEEE Trans. Антенны распространения, Vol. AP-34 (март 1986 г.), стр. 276-280.
^ Barabell, AJ (1998). «Сравнение производительности алгоритмов обработки массивов сверхвысокого разрешения. Пересмотрено» (PDF). Массачусетский технологический институт Lexington Lincoln Lab.
^ Р. Рой и Т. Кайлат, « ESPRIT-оценка параметров сигнала с помощью методов инвариантности вращения », в IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, нет. 7, pp. 984-995, июль 1989 г.
^ Фишлер, Эран и Х. Винсент Бедный. « Оценка количества источников в несбалансированных массивах с помощью критериев теории информации ». IEEE Transactions по обработке сигналов 53.9 (2005): 3543-3553.
^ Abeida, Habti; Чжан, Цилинь; Ли, Цзянь; Мерабтин, Наджим (2013). «Итерационные подходы на основе разреженной асимптотики с минимальной дисперсией для обработки массивов». Транзакции IEEE по обработке сигналов. Институт инженеров по электротехнике и радиоэлектронике (IEEE). 61 (4): 933–944. arXiv : 1802.03070. Bibcode : 2013ITSP... 61..933A. DOI : 10.1109 / tsp.2012.2231676. ISSN 1053-587X.
^ Чжан, Цилинь; Абейда, Хабти; Сюэ, Мин; Роу, Уильям; Ли, Цзянь (2012). «Быстрая реализация разреженной итерационной ковариационной оценки для локализации источника». Журнал акустического общества Америки. 131 (2): 1249–1259. Bibcode : 2012ASAJ..131.1249Z. DOI : 10.1121 / 1.3672656. PMID 22352499.
^ Devaney, AJ (2005-05-01). «Визуализация с обращением времени скрытых целей из мультистатических данных». Транзакции IEEE по антеннам и распространению. 53 (5): 1600–1610. Bibcode : 2005ITAP... 53.1600D. DOI : 10.1109 / TAP.2005.846723. ISSN 0018-926X.
^ Ciuonzo, D.; Romano, G.; Солимен, Р. (2015-05-01). "Анализ производительности МУЗЫКИ с обращением времени". Транзакции IEEE по обработке сигналов. 63 (10): 2650–2662. Bibcode : 2015ITSP... 63.2650C. DOI : 10.1109 / TSP.2015.2417507. ISSN 1053-587X.
^ «Данные и сигнал - ИТ-решения, быстрое обнаружение частоты сверхвысокого разрешения с использованием алгоритма MUSIC». Архивировано из оригинала на 2019-06-26. Проверено 14 июля 2018. Цитировать журнал требует |journal=( помощь )

дальнейшее чтение

Оценка и отслеживание частоты, Куинн и Ханнан, Cambridge University Press 2001.