FDA-MIMO雷达参数估计

引言

FDA-MIMO雷达通过在发射端阵元间引入一个远小于信号载频的频偏，使得发射导向矢量与距离和角度相关。利用距离维的额外自由度，FDA-MIMO可以无模糊地估计出目标距离和角度参数。

FDA-MIMO的信号模型

第m个阵元的发射信号：
$$s_m\left(t\right)={\varphi }_m\left(t\right)e^{j2\pi f_{m}t},0\le t\le T_p,m=1,2,...,M,$$
其中，$T_p$为脉冲持续时间， ${\varphi }_m\left(t\right)$为第m个天线的基带调制信号，M个天线发射的信号正交。

在接收端经过混频、匹配滤波等操作后，第n个天线单元输出的m路信号：
$$y_n={\left[y_{1,n},y_{2,n},\dots ,y_{M,n}\right]}^{\text{T}}={\xi }_k\left[\begin{array}{c}1\\ e^{-j4\pi \frac{\Delta f}{c}{r}_k+j2\pi \frac{d_{\text{T}}}{{\lambda }_0}\sin {\theta }_k}\\ ⋮\\ e^{-j4\pi \frac{\Delta f}{c}(M-1)r_k+j2\pi \frac{d_{\text{T}}}{{\lambda }_0}(M-1)\sin {\theta }_k}\end{array}\right]e^{j2\pi \frac{d_{\text{R}}}{{\lambda }_0}(n-1)\sin {\theta }_k}.$$
将N个接收天线输出的MN路信号排成MN×1的列向量：
$$\begin{array}{rl}& x(l)={\left[y_1^{\text{T}},y_2^{\text{T}},\dots ,y_N^{\text{T}}\right]}^{\text{T}}\\ & =As(l)+n(l)\\ & =\left(A_r\circ A_t\right)s(l)+n(l),l=1,2,\dots L,\end{array}$$
其中，$A$为联合收发域导向矢量，$A_r,A_t$分别为发射域和接收域导向矢量。

MUSIC算法

首先计算采样协方差矩阵：
$${\hat{R}}_x=\frac{1}{L}\sum _{l=1}^{L}x(l)x^{\text{H}}(l)=\sum _{i=1}^K{\lambda }_i{e}_i{e}_i^{\text{H}}+\sum _{i=K+1}^{MN}{\lambda }_i{e}_i{e}_i^{\text{H}}={\hat{E}}_s{\Lambda }_s{\hat{E}}_s^{\text{H}}+{\sigma }_n^2{\hat{E}}_n{\hat{E}}_n^{\text{H}},$$
根据噪声子空间计算二维空间谱函数：
$$P(\theta ,r)=\frac{1}{a^{\text{H}}(\theta ,r){\hat{E}}_n{\hat{E}}_n^{\text{H}}a(\theta ,r)}=\frac{1}{{\Vert {\hat{E}}_n^{\text{H}}a(\theta ,r)\Vert }^2}.$$
遍历距离和角度的搜索区间，通过寻找谱函数的峰值来估计目标距离和角度:
$$\left(\hat{\theta },\hat{r}\right)=\underset{\theta ,r}{\arg }\text{ max}P(\theta ,r)$$

仿真实验

仿真条件：载频$f_0=10\mathrm{GHz}$，频率增量$\Delta f=30\mathrm{KHz}$，发射天线数和接收天线数都设置为8，四个远场不相关目标的距离分别为[8 8 10 10]*1e3m，角度为[-20 20 -20 20]${}^{\circ }$，以下分别是MIMO雷达和FDA-MIMO雷达的距-角度二维空间谱：


可以看出，MIMO雷达无法识别角度相同但分布在不同距离处的目标，而FDA-MIMO雷达由于导向矢量中包含目标的真实距离和角度，可以实现距离和角度的联合估计。

// 2D-MUSIC算法代码
//MIMO空间谱
for i=1:length(range) //range,angle分别为距离角度搜索区间
    for j=1:length(angle)
        at = exp(1i*2*pi*f0/c*dr*sind(angle(j))*(0:N-1)');  
        ar = exp(1i*2*pi*f0/c*dr*sind(angle(j))*(0:M-1)');
        atr = kron(ar,at);  //克罗内克积
        P(i,j) = 1/abs(atr'*En*En'*atr); //En为噪声子空间
    end
end
P = P/max(max(P));
figure(); 
imagesc(range,angle,10*log10(P.'));
xlabel('range'); ylabel('angle'); 
axis xy; axis tight; colorbar;colormap(jet)

//FDA-MIMO空间谱
for i=1:length(range)
    for j=1:length(angle)
        at_d = exp(-1i*2*pi*2*delta_f/c*range(i)*(0:N-1)' );
        at_r = exp(1i*2*pi*f0/c*dt*(0:N-1)'*sind(angle(j)));
        at = at_d.*at_r;
        ar = exp(1i*2*pi*f0/c*dr*sind(angle(j))*(0:M-1)');
        atr = kron(ar,at);
        P(i,j) = 1/abs(atr'*En*En'*atr);
    end
end
P = P/max(max(P));
figure(); 
imagesc(range,angle,10*log10(P.'));
xlabel('range'); ylabel('angle'); 
axis xy; axis tight; colorbar; colormap(jet)