plotSpectrum

Системный объект: поэтапный. MUSICEstimator2D
Пакет: поэтапный

Постройте 2D спектр MUSIC

Синтаксис

plotSpectrum(estimator)
output_args = method(estimator,Name,Value)
lh = plotSpectrum(___)

Описание

plotSpectrum(estimator) строит 2D MUSIC пространственный спектр, вычисленный новым осуществлением метода step для phased.MUSICEstimator2D, estimator.

output_args = method(estimator,Name,Value) строит 2D MUSIC пространственный спектр с дополнительными опциями, заданными одним или несколькими аргументами пары Name,Value.

lh = plotSpectrum(___) возвращает указатель строки на фигуру.

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

развернуть все

2D средство оценки MUSIC, заданное как Системный объект phased.MUSICEstimator2D.

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

Укажите необязательные аргументы в виде пар ""имя, значение"", разделенных запятыми. Имя (Name) — это имя аргумента, а значение (Value) — соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Модули используются для графического вывода, заданные как пара, разделенная запятой, состоящая из 'Unit' и 'db', 'mag' или 'pow'.

Пример:

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

Постройте нормированный спектр, заданный как пара, разделенная запятой, состоящая из 'NormalizedResponse' и false или true. Нормализация устанавливает значение самого большого значения спектра одному.

Пример: true

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

Заголовок графика, заданного как пара, разделенная запятой, состоящая из 'Title' и вектора символов.

Пример: true

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

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

развернуть все

Указатель строки графика.

Примеры

развернуть все

Примите, что две синусоидальных волны частот 450 Гц и 600 Гц ударяют URA от двух различных направлений. Сигналы прибывают от азимута на-37 °, повышения на 0 ° и азимута на 17 °, повышения на 20 °. Используйте 2D MUSIC, чтобы оценить направления прибытия двух сигналов. Массив рабочая частота составляет 150 МГц и частота дискретизации сигнала, составляет 8 кГц.

f1 = 450.0;
f2 = 600.0;
doa1 = [-37;0];
doa2 = [17;20];
fc = 150e6;
c = physconst('LightSpeed');
lam = c/fc;
fs = 8000;

Создайте URA с изотропными элементами по умолчанию. Установите область значений частотной характеристики элементов.

array = phased.URA('Size',[11 11],'ElementSpacing',[lam/2 lam/2]);
array.Element.FrequencyRange = [50.0e6 500.0e6];

Создайте два сигнала и добавьте случайный шум.

t = (0:1/fs:1).';
x1 = cos(2*pi*t*f1);
x2 = cos(2*pi*t*f2);
x = collectPlaneWave(array,[x1 x2],[doa1,doa2],fc);
noise = 0.1*(randn(size(x))+1i*randn(size(x)));

Создайте и выполните 2D средство оценки MUSIC, чтобы найти направления прибытия.

estimator = phased.MUSICEstimator2D('SensorArray',array,...
    'OperatingFrequency',fc,...
    'NumSignalsSource','Property',...
    'DOAOutputPort',true,'NumSignals',2,...
    'AzimuthScanAngles',-50:.5:50,...
    'ElevationScanAngles',-30:.5:30);
[~,doas] = estimator(x + noise)
doas = 2×2

   -37    17
     0    20

Предполагаемые DOAs точно совпадают с истинным DOAs.

Постройте 2D пространственный спектр

plotSpectrum(estimator);

Введенный в R2017b