2D MVDR (Каплун) пространственное средство оценки спектра
MVDREstimator2D объект вычисляет 2D минимальное отклонение ответ без искажений (MVDR) пространственная оценка спектра. Это средство оценки DOA также упоминается как средство оценки Каплуна.
Оценить пространственный спектр:
Задайте и настройте свое 2D пространственное средство оценки спектра MVDR. Смотрите Конструкцию.
Вызовите step оценить пространственный спектр согласно свойствам phased.MVDREstimator2D. Поведение step характерно для каждого объекта в тулбоксе.
Примечание
Запуск в R2016b, вместо того, чтобы использовать step метод, чтобы выполнить операцию, заданную Системой object™, можно вызвать объект с аргументами, как будто это была функция. Например, y = step(obj,x) и y = obj(x) выполните эквивалентные операции.
H = phased.MVDREstimator2D создает 2D пространственный Системный объект средства оценки спектра MVDR, H. Объект оценивает пространственный спектр сигнала с помощью узкополосного формирователя луча MVDR.
H = phased.MVDREstimator2D( создает объект, Name,Value)H, с каждым заданным набором имени свойства к заданному значению. Можно задать дополнительные аргументы пары "имя-значение" в любом порядке как (Name1, Value1..., NameN, ValueN).
|
Обработайте к сенсорной матрице Задайте сенсорную матрицу как указатель. Сенсорная матрица должна быть объектом массивов в Значение по умолчанию: |
|
Скорость распространения сигнала Задайте скорость распространения сигнала, в метрах в секунду, как положительная скалярная величина. Можно задать это свойство как одинарную или двойную точность. Значение по умолчанию: Скорость света |
|
Система рабочая частота Задайте рабочую частоту системы в герц как положительная скалярная величина. Значение по умолчанию соответствует 300 МГц. Можно задать это свойство как одинарную или двойную точность. Значение по умолчанию: |
|
Количество битов квантования фазовращателя Количество битов раньше квантовало компонент сдвига фазы формирователя луча или регулирующий векторные веса. Задайте количество битов как неотрицательное целое число. Значение нуля указывает, что никакое квантование не выполняется. Можно задать это свойство как одинарную или двойную точность. Значение по умолчанию: |
|
Выполните прямое обратное усреднение Установите это свойство на По умолчанию: false |
|
Углы сканирования азимута (степени) Задайте углы сканирования азимута (в градусах) как вектор действительных чисел. Углы должны быть между –180 и 180, включительно. Необходимо задать углы в порядке возрастания. Можно задать это свойство как одинарную или двойную точность. Значение по умолчанию: |
|
Углы сканирования вертикального изменения Задайте углы сканирования вертикального изменения (в градусах) как вектор действительных чисел или скаляр. Углы должны быть между –90 и 90, включительно. Необходимо задать углы в порядке возрастания. Можно задать это свойство как одинарную или двойную точность. Значение по умолчанию: |
|
Enable DOA выход Чтобы получить направление прибытия (DOA) сигнала, установите это свойство на По умолчанию: false |
|
Количество сигналов Задайте количество сигналов для оценки DOA как положительное скалярное целое число. Это свойство применяется, когда вы устанавливаете Значение по умолчанию: |
| plotSpectrum | Постройте пространственный спектр |
| сброс | Сбросьте состояния 2D пространственного объекта средства оценки спектра MVDR |
| шаг | Выполните пространственную оценку спектра |
| Характерный для всех системных объектов | |
|---|---|
release | Позвольте изменения значения свойства Системного объекта |
Оцените DOAs двух сигналов, полученных URA с 50 элементами с прямоугольной решеткой. Антенна рабочая частота составляет 150 МГц. Фактическое направление первого сигнала является −37 ° в азимуте и 0 ° в вертикальном изменении. Направление второго сигнала составляет 17 ° в азимуте и степени на 20 ° в области вертикального изменения. Затем постройте пространственный спектр.
Создайте прибывающие сигналы.
fs = 8000; t = (0:1/fs:1).'; x1 = cos(2*pi*t*300); x2 = cos(2*pi*t*400); array = phased.URA('Size',[5 10],'ElementSpacing',[1 0.6]); array.Element.FrequencyRange = [100e6 300e6]; fc = 150e6; x = collectPlaneWave(array,[x1 x2],[-37 0;17 20]',fc);
Добавьте шум.
noise = 0.1*(randn(size(x))+1i*randn(size(x)));
Создайте средство оценки DOA MVDR и оцените DOAs.
estimator = phased.MVDREstimator2D('SensorArray',array,... 'OperatingFrequency',fc,... 'DOAOutputPort',true,'NumSignals',2,... 'AzimuthScanAngles',-50:50,... 'ElevationScanAngles',-30:30); [~,doas] = estimator(x + noise);
Постройте спектр.
plotSpectrum(estimator)
[1] Деревья фургона, H. Оптимальная обработка матриц. Нью-Йорк: Wiley-межнаука, 2002.
