✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
⛄ 内容介绍
与相控阵发射波束只有角度依赖性不同,频控阵(FDA)通过在阵元上附加一个频偏增量,使其具有距离-角度二维依赖特性,不仅解决了波束发射时的角度指向问题,还可以解决波束发射时的距离指向问题.通过在FDA中引入MIMO,可以将发射信息整合到接收端,从而实现更加灵活的距离维自由度利用,提升距离-角度二维联合信号处理的能力,这使得FDA-MIMO在抗距离维主瓣欺骗干扰方面具有独特的应用优势.
⛄ 部分代码
function [DOA_esprit,r_esprit]=esprit(X,K,M,N,delta_f,c,J)
R_x=X*X'/J; %the covariance matrix
[V,D]=eig(R_x); %eigenvalue decomposition
[D,index]=sort(diag(D));
en=V(:,index(1:M*N-K));es=V(:,index(M*N-K+1:M*N)) ; % noise and signal subspace
Es3=[];Es4=[];
for hh=1:N
Es3=[Es3;es((1:M-1)+M*(hh-1),:)];
Es4=[Es4;es((2:M)+M*(hh-1),:)];
end
RzzEs3Es4 = [Es3'; Es4'] * [Es3 Es4]; %%% TLS-ESPRIT
[EVec1, EVal1] = eig(RzzEs3Es4);
[EVal1, Indx1] = sort(diag(EVal1).');
for Col = 1 : 2*K
E(:,Col) = EVec1(:,Indx1(2*K+1-Col));
end;
E12 = E(1:K,K+1:2*K);
E22 = E(K+1:2*K,K+1:2*K);
[Vt,THETA_T] = eig(-1*E12/E22);
THETA_T=diag(THETA_T);
Es1=es((1:M*N-M),:);%%% TLS-ESPRIT
Es2= es(M+1:M*N,:);
RzzEs1Es2 = [Es1'; Es2'] * [Es1 Es2];
[EVec, EVal] = eig(RzzEs1Es2);
[EVal, Indx] = sort(diag(EVal).');
for Col = 1 : 2*K
E(:,Col) = EVec(:,Indx(2*K+1-Col));
end;
E12 = E(1:K,K+1:2*K);
E22 = E(K+1:2*K,K+1:2*K);
[Vr,THETA_R] = eig(-1*E12/E22);
THETA_R=diag(THETA_R);
R_pair=inv(Vt)*Vr; %%% a new simple pairing method
for ppaa=1:K
iirr=find(R_pair(ppaa,:)==max(R_pair(ppaa,:)));
THETA_R1(ppaa)= THETA_R(iirr)*180/pi;
end
DOA_esprit=asin(angle(THETA_R)/pi)*180/pi; %%% estimated DOA
% DOD_esprit=asin(angle(THETA_T)/pi)*180/pi;
%r_esprit=(angle(THETA_R(2,:))-angle(THETA_T(2,:)))*c/(4*delta_f);
r_esprit=(angle(THETA_R1).'-angle(THETA_T))*c/(4*pi*delta_f);
⛄ 运行结果
⛄ 参考文献
[1] 刘润东, 薛峰涛, 杨赟秀, et al. FDA-MIMO雷达实值降维求根MUSIC参数估计算法[J]. 现代雷达, 2022, 44(4):7.
[2] 朱圣棋, 王可, 王宇卓,等. 机载FDA-MIMO双基地雷达距离模糊杂波抑制及降维搜索方法:, CN113253222A[P]. 2021.
[3] 陈浩, 马建朝, 吕明久,等. 基于压缩感知的FDA-MIMO雷达波束形成算法[J]. 空军预警学院学报, 2020(004):034.
[4] 李星星, 王党卫, 范甘霖,等. 一种改进的FDA-MIMO雷达波束形成算法[C]// 中国电子学会. 中国电子学会, 2017.