Вычислите направления прибытия 3 некоррелированых сигналов, прибывающих в ULA с 11 элементами с интервалом полудлины волны. Примите, что сигналы прибывают из поперечных углов 0 °,-12 ° и 85 °. Шум в каждом элементе является Гауссовым белым шумом и является некоррелированым между элементами. ОСШ составляет 5 дБ.

Задайте количество элементов ULA и интервала элемента (в длинах волн).

nelem = 11;
d = 0.5;
snr = 5.0;
elementPos = (0:nelem-1)*d;

Задайте количество сигналов и их поперечных углов падения.

nsig = 3;
angles = [0.0 -12.0 85.0];

Создайте ковариационную матрицу датчика.

covmat = sensorcov(elementPos,angles,db2pow(-snr));

Оцените поперечные углы падения.

doas = musicdoa(covmat,nsig)

doas = 1×3

    85     0   -12

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

Вычислите направления прибытия 3 некоррелированых сигналов, прибывающих в ULA с 11 элементами с интервалом полудлины волны. Примите, что сигналы прибывают из поперечных углов 0 °,-12 ° и 85 °. Шум в каждом элементе является Гауссовым белым шумом и является некоррелированым между элементами. ОСШ составляет 2 дБ.

Задайте количество элементов ULA и интервала элемента (в длинах волн).

nelem = 11;
d = 0.5;
snr = 2.0;
elementPos = (0:nelem-1)*d;

Задайте количество сигналов и их поперечных углов падения.

nsig = 3;
angles = [0.0 -12.0 85.0];

Создайте ковариационную матрицу датчика.

covmat = sensorcov(elementPos,angles,db2pow(-snr));

Вычислите спектр MUSIC и оцените поперечные углы падения.

[doas,spec,specang] = musicdoa(covmat,nsig);

Постройте спектр MUSIC.

plot(specang,10*log10(spec))
xlabel('Arrival Angle (deg)')
ylabel('Magnitude (dB)')
title('MUSIC Spectrum')
grid

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

Вычислите направления прибытия 4 некоррелированых сигналов, прибывающих в ULA с 11 элементами. Интервал элемента является 0,5 длинами волн. Примите, что сигналы прибывают из поперечных углов-60.2 °,-20.7 °, 0,5 ° и 84,8 °. Шум в каждом элементе является Гауссовым белым шумом и является некоррелированым между элементами. ОСШ составляет 0 дБ.

Задайте количество элементов ULA и интервала элемента (в длинах волн).

nelem = 11;
d = 0.5;
snr = 5.0;
elementPos = (0:nelem-1)*d;

Задайте количество сигналов и их поперечных углов падения.

nsig = 4;
angles = [-60.2 -20.7 0.5 84.8];

Создайте ковариационную матрицу датчика.

covmat = sensorcov(elementPos,angles,db2pow(-snr));

Вычислите спектр MUSIC и оцените поперечные углы падения в диапазоне от-70 ° до 90 ° с шагом на 0,1 °.

[doas,spec,specang] = musicdoa(covmat,nsig,'ScanAngles',[-70:.1:90]);

Постройте спектр MUSIC.

plot(specang,10*log10(spec))
xlabel('Arrival Angle (deg)')
ylabel('Magnitude (dB)')
title('MUSIC Spectrum')
grid

disp(doas)

   84.8000    0.5000  -60.2000  -20.7000

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

Вычислите направления прибытия 4 некоррелированых сигналов, прибывающих в ULA с 11 элементами. Интервал элемента является 0,4 интервалами длин волн. Примите, что сигналы прибывают из поперечных углов-60 °,-20 °, 0 ° и 85 °. Шум в каждом элементе является Гауссовым белым шумом и является некоррелированым между элементами. ОСШ составляет 0 дБ.

Задайте количество элементов ULA и интервала элемента (в длинах волн).

nelem = 11;
d = 0.4;
snr = 0.0;
elementPos = (0:nelem-1)*d;

Задайте количество сигналов и их поперечных углов падения.

nsig = 4;
angles = [-60.0 -20.0 0.0 85.0];

Создайте ковариационную матрицу датчика.

covmat = sensorcov(elementPos,angles,db2pow(-snr));

Вычислите спектр MUSIC и оцените поперечные углы падения.

[doas,spec,specang] = musicdoa(covmat,nsig,'ElementSpacing',d);

Постройте спектр MUSIC.

plot(specang,10*log10(spec))
xlabel('Arrival Angle (deg)')
ylabel('Magnitude (dB)')
title('MUSIC Spectrum')
grid

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