Вычислите направления прихода 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

   -12     0    85

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

Вычислите направления прихода 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)

    0.5000   84.8000  -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

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