plotSpectrum
Системный объект: поэтапный. MUSICEstimator2D
Пакет: поэтапный
Постройте 2D спектр MUSIC
Синтаксис
plotSpectrum(estimator)
output_args = method(estimator,Name,Value)
lh = plotSpectrum(___)
Описание
plotSpectrum( строит 2D MUSIC пространственный спектр, вычисленный новым estimator)step осуществление метода для phased.MUSICEstimator2D, estimator.
output_args = method(estimator,Name,Value)Name,Value парные аргументы.
lh = plotSpectrum(___) возвращает указатель линии на фигуру.
Входные параметры
Выходные аргументы
Примеры
Оцените DOAs двух сигналов
Примите, что две синусоидальных волны частот 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);
