rootmusicdoa

Направление прибытия с помощью Root MUSIC

Описание

пример

ang = rootmusicdoa(R,nsig) оценивает направления прибытия, ang, из набора плоских волн, полученных на равномерной линии массива (ULA). Оценка использует алгоритм root MUSIC. Входными параметрами являются оцененная пространственная ковариационная матрица между элементами датчика, Rи количество поступающих сигналов, nsig. В этом синтаксисе элементы датчика разнесены на половину длины волны.

пример

ang = rootmusicdoa(___,'Name','Value') позволяет вам задать дополнительные входные параметры в виде пар "имя-значение". Этот синтаксис может использовать любой из входных параметров в предыдущем синтаксисе.

Примеры

свернуть все

Предположим, что равномерная матрица с интервалом линии половину длины волны с 10 элементами. Три плоские волны прибывают из направления 0 °, -25 ° и 30 ° азимута. Углы возвышения 0 °. Шум является пространственно и временно белым Гауссовым шумом.

Установите ОСШ для каждого сигнала равным 5 дБ. Найдите углы прибытия.

N = 10;
d = 0.5;
elementPos = (0:N-1)*d;
angles = [0 -25 30];
Nsig = 3;
R = sensorcov(elementPos,angles,db2pow(-5));
doa = rootmusicdoa(R,Nsig)
doa = 1×3

   -0.0000  -25.0000   30.0000

Эти углы согласуются с известными входными углами.

Предположим, что равномерный массив линии 10 элементами, как в предыдущем примере. Но теперь интервал между элементами меньше половины длины волны. Три плоские волны прибывают из направления 0 °, -25 ° и 30 ° азимута. Углы возвышения 0 °. Шум является пространственно и временно белым Гауссовым шумом. ОСШ для каждого сигнала составляет 5 дБ.

Установите интервал между элементами равным 0,4 длины волны, используя ElementSpacing Пара "имя-значение". Затем найдите углы прибытия.

N = 10;
d = 0.4;
elementPos = (0:N-1)*d;
angles = [0 -25 30];
Nsig = 3;
R = sensorcov(elementPos,angles,db2pow(-5));
doa = rootmusicdoa(R,Nsig,'ElementSpacing',d)
doa = 1×3

   -0.0000  -25.0000   30.0000

Решение согласуется с известными углами.

Входные параметры

свернуть все

Пространственная ковариационная матрица, заданная как комплексная, положительно-определенная, N -by - N матрица. В этой матрице N представляет количество элементов в массиве ULA. Если R не является Эрмитовым, Эрмитова матрица формируется путем усреднения матрицы и ее сопряженного транспонирования, (R+R')/2.

Пример: [4.3162, -0.2777 -0.2337i; -0.2777 + 0.2337i, 4.3162]

Типы данных: double
Поддержка комплексного числа: Да

Количество поступающих сигналов, заданное в виде положительного целого числа. Количество сигналов должно быть меньше, чем количество элементов в массиве ULA.

Пример: 2

Типы данных: double

Аргументы в виде пар имя-значение

Задайте необязательные разделенные разделенными запятой парами Name,Value аргументы. Name - имя аргумента и Value - соответствующее значение. Name должны находиться внутри кавычек. Можно задать несколько аргументов в виде пар имен и значений в любом порядке Name1,Value1,...,NameN,ValueN.

Пример: «ElementSpacing», 0,4

Интервал между элементами ULA, заданный как действительный, положительный скаляр. Модули измерения положения измеряются в терминах длины волны сигнала.

Пример: 0.4

Типы данных: double

Выходные аргументы

свернуть все

Направления угла прихода, возвращенные как действительный, 1-бай- M вектор. Размерность M является количеством поступающих сигналов, заданным в аргументе nsig. Модулями измерения угла являются степени, а значениями угла являются значения от -90 ° до 90 °.

Ссылки

[1] Van Trees, H.L. Optimum Array Processing. Нью-Йорк: Wiley-Interscience, 2002.

Расширенные возможности

.
Введенный в R2013a