空间平滑doa估计

 
 
clear all;close all;clc;
lamda = 1
d=0.5;
m=10;
trial = 100
for t = 1:1:trial
theta=[-25,30];
snum = length(theta)
f=400;
snapshot=2000;
S1=2.72*exp((j*2*pi*f*[1:snapshot]/snapshot)+pi/10); 
S2=4.48*exp((j*2*pi*f*[1:snapshot]/snapshot)+pi/2); 
S=[S1;S2];
 
A1=exp(-j*2*pi*d*[0:m-1]'*sin(theta(1)*pi/180)/lamda);
A2=exp(-j*2*pi*d*[0:m-1]'*sin(theta(2)*pi/180)/lamda);
A=[A1,A2];
N= wgn(10,2000,0);
s_power1=10*log(2.72^2/2);
s_power2=10*log(4.48^2/2);
snr1=s_power1;
snr2=s_power2;
X=A*S+N;
X1 = X(1:m-1,:)
X2 = X(2:m,:)
Rx1x1 = X1*X1'/snapshot;
Rx2x2 = X2*X2'/snapshot
Rxx = 1/2*(Rx1x1+Rx2x2 )