Matlab实现基于ESPRIT算法的单基地MIMO雷达测向算法

%LCESPRITMONO% clear all close all clc %basic parameter% pt=5;%tagart number mt=8;%transmitter number mr=8;%receiver number fc=10000;%carrier frequency for n=1:pt ft(n)=300+n*100; end lamda=1; dt=lamda/2; dr=lamda/2; SNR=10; ns=1000;%snapshot number t=0:1/fc:(ns-1)/fc; %theta=-50+rand(1,pt)*100;%transmit angle theta=[-40 0 5 -25 25]; %phi=-50+rand(1,pt)*100;%reveive angle phi=theta; dtt=0:dt:(mt-1)*dt;dtt=dtt'; drr=0:dr:(mr-1)*dr;drr=drr'; A1=zeros(mr,pt); A2=zeros(mt,pt); A=zeros(mr*mt,pt); s=zeros(pt,ns); for n=1:pt for m=1:mr A1(m,n)=exp(-j*2*pi*drr(m)/lamda*sin(phi(n)/180*pi)); end end for n=1:pt for m=1:mt A2(m,n)=exp(-j*2*pi*dtt(m)/lamda*sin(theta(n)/180*pi)); end end for n=1:pt A(:,n)=kron(A1(:,n),A2(:,n)); end noise=(randn(mt*mr,ns)+j*randn(mt*mr,ns))/sqrt(2); for n=1:pt amp(n)=sqrt(cov(noise(1,:))*10^(SNR/10)); end for n=1:pt s(n,:)=amp(n)*sin(2*pi*ft(n)*t); end X=A*s+noise; G=zeros(mt*mr,mt+mr-1); for n=1:mr G((1+(n-1)*mt):n*mt,n:(mt+n-1))=eye(mt); end W=G'*G; for n=1:(mt+mr-1) W(n,n)=W(n,n)^(-1/2); end Y=W*G'*X; ry=(Y*Y')/ns; [V,D]=eigs(ry,pt,'LM'); E1=V(1:(mt+mr-2),:); E2=V(2:(mt+mr-1),:); E=inv(E2'*E2)*E2'*E1; [VEC,VALUE]=eigs(E); monophi=zeros(pt,1); for n=1:pt monophi(n)=asin(angle(VALUE(n,n))/pi)*180/pi; %?????? end plot(monophi,5,'*r');hold on xlabel('DOAs(degree)'); grid on; title('LS-ESPRIT Algorithm DOA Estimation For MONOSTATIC MIMO Radar');