Построить 10-элементную полуволновую однородную линейную матрицу, принимающую две плоские волны, поступающие от азимута 0 ° и -25 °. Оба угла возвышения равны 0 °. Предположим, что два сигнала частично коррелированы. SNR для каждого сигнала составляет 5 дБ. Шум является пространственным и временным гауссовым белым шумом. Сначала создайте пространственную ковариационную матрицу из сигнала и шума. Затем решите для количества сигналов, используяrootmusicdoa. Затем выполните пространственное сглаживание ковариационной матрицы, используя spsmoothи снова решить для углов прихода сигнала с помощью rootmusicdoa.

Настройте массив и сигналы. Затем генерируют пространственную ковариационную матрицу для матрицы из сигналов и шума.

N = 10;
d = 0.5;
elementPos = (0:N-1)*d;
angles = [0 -25];
ac = [1 1/5];
scov = ac'*ac;
R = sensorcov(elementPos,angles,db2pow(-5),scov);

Решите углы прихода с помощью исходной ковариационной матрицы.

Nsig = 2;
doa =  rootmusicdoa(R,Nsig)

doa = 1×2

    0.2099   41.3166

Решаемые углы прихода неверны - они не согласуются с известными углами прихода, используемыми для создания ковариационной матрицы.

Затем решите углы прихода с помощью сглаженной ковариационной матрицы. Выполнить пространственное сглаживание для обнаружения L-1 когерентных сигналов. Выберите L = 3.

Nsig = 2;
L = 2;
RSM = spsmooth(R,L);
doasm = rootmusicdoa(RSM,Nsig)

doasm = 1×2

  -25.0000   -0.0000

В этом случае вычисленные углы согласуются с известными углами прихода.