【通信】基于增广矩阵束的L型阵列的二维DOA估计附matlab完整代码

1 内容介绍

基于增广矩阵束的L型阵列的二维DOA估计附matlab完整代码

2 完整代码

%

clear all
close all
derad = pi/180;
radeg = 180/pi;
twpi = 2*pi;
kelmx = 8;   
kelmy = 10;             % 
dd = 0.5;               %
dx=0:dd:(kelmx-1)*dd;   %
dy=0:dd:(kelmy-1)*dd;
iwave = 3;              % number of DOA
L=iwave;
theta1 = [10 30 50];  
theta2 = [15 35 55];
snr = 20;% input SNR (dB)           
n = 200;                % 
Ax=exp(-i*twpi*dx.'*(sin(theta1*derad).*cos(theta2*derad)));
Ay=exp(-i*twpi*dy.'*(sin(theta1*derad).*sin(theta2*derad)));%%%% direction matrix
S=randn(iwave,n);
X0=Ax*S;
X=awgn(X0,snr,'measured');
Y0=Ay*S;
Y=awgn(Y0,snr,'measured');
Rxy=X*Y';
P=5;
Q=6;
Re=[];
for kk=1:kelm