step
Системный объект: фазированный. MUSICEstimator
Пакет: поэтапный
Оценить направление прибытия используя MUSIC
Синтаксис
spectrum = step(estimator,X)
[spectrum,doa]
= step(estimator,X)
Описание
Примечание
Вместо использования step метод для выполнения операции, заданной Системной object™, можно вызвать объект с аргументами, как если бы это была функция. Для примера, y = step(obj,x) и y = obj(x) выполнять эквивалентные операции.
spectrum = step(estimator,X)X.
[ также возвращает широкие направления прихода сигнала, spectrum,doa]
= step(estimator,X)doa. Чтобы использовать этот синтаксис, установите DOAOutputPort свойство к true.
Примечание
Объект выполняет инициализацию при первом выполнении объекта. Эта инициализация блокирует нетронутые свойства и входные спецификации, такие как размерности, сложность и тип данных входных данных. Если вы изменяете свойство nontunable или спецификацию входа, системный объект выдает ошибку. Чтобы изменить нетронутые свойства или входы, необходимо сначала вызвать release метод для разблокировки объекта.
Входные параметры
Выходные аргументы
Примеры
Постройте спектр MUSIC двух сигналов, поступающих в ULA
Оцените DOA двух сигналов, принятых стандартным ULA с 10 элементами, имеющим интервал между элементами 1 метр. Затем постройте график спектра MUSIC.
Примечание: Можно заменить каждый вызов функции эквивалентным step синтаксис. Для примера замените myObject(x) с step(myObject,x).
Создайте массив ULA. Рабочая частота антенны составляет 150 МГц.
fc = 150.0e6; array = phased.ULA('NumElements',10,'ElementSpacing',1.0);
Создайте поступающие сигналы в ULA. Истинное направление прихода первого сигнала составляет 10 ° по азимуту и 20 ° по повышению. Направление второго сигнала составляет 60 ° по азимуту и -5 ° по повышению.
fs = 8000.0; t = (0:1/fs:1).'; sig1 = cos(2*pi*t*300.0); sig2 = cos(2*pi*t*400.0); sig = collectPlaneWave(array,[sig1 sig2],[10 20; 60 -5]',fc); noise = 0.1*(randn(size(sig)) + 1i*randn(size(sig)));
Оцените DOAs.
estimator = phased.MUSICEstimator('SensorArray',array,... 'OperatingFrequency',fc,... 'DOAOutputPort',true,'NumSignalsSource','Property',... 'NumSignals',2); [y,doas] = estimator(sig + noise); doas = broadside2az(sort(doas),[20 -5])
doas = 1×2
9.5829 60.3813
Постройте график спектра MUSIC.
plotSpectrum(estimator,'NormalizeResponse',true)
Вычисление DOA двух ближайших сигналов с использованием MUSIC
Сначала оцените DOA двух сигналов, принятых стандартным ULA с 10 элементами, имеющим интервал между элементами половину длины волны. Затем постройте график пространственного спектра.
Примечание: Можно заменить каждый вызов функции эквивалентным step синтаксис. Для примера замените myObject(x) с step(myObject,x).
Рабочая частота антенны составляет 150 МГц. Направления прихода двух сигналов разделены на 2 °. Направление первого сигнала составляет 30 ° азимут и 0 ° повышения. Направление второго сигнала составляет 32 ° азимута и 0 ° повышения. Оцените количество сигналов с помощью критерия минимальной длины описания (MDL).
Создайте сигналы, поступающие в ULA.
fs = 8000; t = (0:1/fs:1).'; f1 = 300.0; f2 = 600.0; sig1 = cos(2*pi*t*f1); sig2 = cos(2*pi*t*f2); fc = 150.0e6; c = physconst('LightSpeed'); lam = c/fc; array = phased.ULA('NumElements',10,'ElementSpacing',0.5*lam); sig = collectPlaneWave(array,[sig1 sig2],[30 0; 32 0]',fc); noise = 0.1*(randn(size(sig)) + 1i*randn(size(sig)));
Оцените DOAs.
estimator = phased.MUSICEstimator('SensorArray',array,... 'OperatingFrequency',fc,'DOAOutputPort',true,... 'NumSignalsSource','Auto','NumSignalsMethod','MDL'); [y,doas] = estimator(sig + noise); doas = broadside2az(sort(doas),[0 0])
doas = 1×2
30.0000 32.0000
Постройте график спектра MUSIC.
plotSpectrum(estimator,'NormalizeResponse',true)