基于Matlab实现FDA-MIMO子空间算法

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，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.

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料