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