阵列信号基础：天线阵列模型

天线阵列模型

0. 前言

阵列信号处理基础本质上属于参数估计问题，和信道估计知识基本上别无二致。末学在这里整理了阵列信号处理的基础知识。

3. 天线阵列模型¹

天线阵元之间的摆放位置影响着阵列接收信号的数学模型，不同的天线阵列模型有着不一样的应用场景，并影响着阵列信号处理的方法。

3.1 均匀线阵

假设 $M$ 个阵元等距离排列成一条直线，阵元间距为 $d$。假设选取最左边的天线阵元作为参考点，则第 $m$ 个天线阵元相对于参考点的时间延迟可以表示为
$${\tau }_m({\theta }_i)=\frac{(m-1)d\sin {\theta }_i}{v}$$

• 各天线阵元之间的互耦效应忽略不计，且阵元间距为最高频率源信号的半波长，即 $d=\frac{\lambda }{2}$。

• 传播距离远大于阵列大小，即信号在介质中以平面波的形式到达阵列。

• 接收机各个通道拥有完全相同的特性。

根据以上假设，方向向量可以写为
$$\mathbf{a}({\theta }_i)={\left[1,e^{j\pi \sin {\theta }_i},e^{j2\pi \sin {\theta }_i},\cdots ,e^{j(M-1)\pi \sin {\theta }_i}\right]}^T$$

在这里很尴尬地说，我用了好久的模型都是 $e^{-j(k-1)\pi {\beta }_i}$，今天才发现如果参考点为第一个阵元在最左边的话，指数应该是没有符号的，也就是上面的式子是正确的！
%%
当然，如果是如下图这种，参考阵元第一个在最右边，从右往左建模，则需要加符号！
%%
updated by 2019/7/19

3.2 L 型阵

如图为 L 型阵列模型，在 $x$ 轴和 $y$ 轴均 $M$ 个阵元，阵元间距为 $d$。 $K$ 个信源，二维入射角为 $({\theta }_i,{\phi }_i)$。$x$ 轴上的接收信号模型可以视为线性阵列，$y$ 轴同理：
$$\begin{array}{rl}X& =A_{x}S+N_x\\ Y& =A_{y}S+N_y\end{array}$$

其中，方向矩阵为
$$\begin{array}{rl}{A}_x& =\left[a_x({\theta }_1,{\phi }_1),\cdots ,a_x({\theta }_K,{\phi }_K)\right]\\ A_y& =\left[a_y({\theta }_1,{\phi }_1),\cdots ,a_y({\theta }_K,{\phi }_K)\right]\\ a_x({\theta }_i,{\phi }_i)& ={\left[1,e^{-j\frac{2\pi }{\lambda }d\cos {\theta }_i\sin {\phi }_i},\cdots ,e^{-j\frac{2\pi }{\lambda }(M-1)d\cos {\theta }_i\sin {\phi }_i}\right]}^T\\ a_y({\theta }_i,{\phi }_i)& ={\left[1,e^{-j\frac{2\pi }{\lambda }d\sin {\theta }_i\sin {\phi }_i},\cdots ,e^{-j\frac{2\pi }{\lambda }(M-1)d\sin {\theta }_i\sin {\phi }_i}\right]}^T\end{array}$$

3.3 均匀平面阵

假设有 $K$ 个信号源，均匀平面阵列指的是 $M$ 个天线阵元等间距的排列成一个正方形或者矩形，如上图所示由 $M\times N$ 个阵元组成。$x$ 轴方向有 $N$ 个间距为 $d$ 的均匀线阵，$y$ 轴方向有 $M$ 个间距为 $d$ 的均匀线阵。如果选取原点为参考点，另外某个阵元的坐标可以写为 $(x_n,y_m)$，于是时延差可以写为

$$\begin{array}{rl}{\tau }_{n,m}({\phi }_k,{\theta }_k)& =\frac{x_n\sin {\phi }_k\cos {\theta }_k+y_m\sin {\phi }_k\sin {\theta }_k}{v}\\ a_{n,m}({\phi }_k,{\theta }_k)& =\exp \left(-j{w}_0{\tau }_{n,m}({\phi }_k,{\theta }_k)\right)\\ X& =AS+N\\ S& ={\left[s_1(t),\cdots ,s_K(t)\right]}^T\in {\mathbb{C}}^{K\times 1}\\ A& =\left[a({\theta }_1,{\phi }_1),\cdots ,a({\theta }_K,{\phi }_K)\right]\in {\mathbb{C}}^{NM\times K}\end{array}$$
简单解释下这里的时延差公式怎么来的：
由之前的 ULA 基础可知，沿着入射波方向的长度差等于阵元 $(x_n,y_m)$ 到原点直连距离 $r$ 乘以入射波与平面阵法线的夹角角度的正弦 $\sin {\phi }_k$。
即这段沿着入射波方向的长度差为 $r\cdot \sin {\phi }_k$。而在 $xoy$ 平面上，
$$\begin{array}{rl}{x}_n^2+y_m^2& =r^2\\ x_n& =r\cos {\theta }_k\\ y_m& =r\sin {\theta }_k\end{array}$$
于是就有
$$\begin{array}{rl}r& =r({\cos }^2{\theta }_k+{\sin }^2{\theta }_k)\\ & =r\cos {\theta }_k\cdot \cos {\theta }_k+r\sin {\theta }_k\cdot \sin {\theta }_k\\ & =x_n\cos {\theta }_k+y_m\sin {\theta }_k\end{array}$$

从 $x$ 轴对应的 $N$ 个阵元的方向矩阵和从 $y$ 轴对应的 $M$ 个阵元的方向矩阵为
$$\begin{array}{rl}{A}_x& =\left[a_x({\theta }_1,{\phi }_1),\cdots ,a_x({\theta }_K,{\phi }_K)\right]\\ A_y& =\left[a_y({\theta }_1,{\phi }_1),\cdots ,a_y({\theta }_K,{\phi }_K)\right]\\ a_x({\theta }_i,{\phi }_i)& ={\left[1,e^{-j\frac{2\pi }{\lambda }d\cos {\theta }_i\sin {\phi }_i},\cdots ,e^{-j\frac{2\pi }{\lambda }(N-1)d\cos {\theta }_i\sin {\phi }_i}\right]}^T\\ a_y({\theta }_i,{\phi }_i)& ={\left[1,e^{-j\frac{2\pi }{\lambda }d\sin {\theta }_i\sin {\phi }_i},\cdots ,e^{-j\frac{2\pi }{\lambda }(M-1)d\sin {\theta }_i\sin {\phi }_i}\right]}^T\end{array}$$

则可以得到方向矩阵
$$\begin{array}{rl}X& =AS+N\\ S& ={\left[s_1(t),\cdots ,s_K(t)\right]}^T\in {\mathbb{C}}^{K\times 1}\\ A& =\left[a({\theta }_1,{\phi }_1),\cdots ,a({\theta }_K,{\phi }_K)\right]\\ & =\left[\begin{array}{c}{A}_x{D}_1(A_y)\\ A_x{D}_2(A_y)\\ ⋮\\ A_x{D}_M(A_y)\end{array}\right]\in {\mathbb{C}}^{NM\times K}\end{array}$$
其中 $D_i(\cdot )$ 表示取矩阵的第 $i$ 行作为构成对角矩阵的对角元素。

3.4 均匀圆阵

如图所示为均匀圆阵模型，$M$ 个阵元均匀分布在圆周上，假设 $K$ 个信源，二维入射角为 $({\theta }_i,{\phi }_i)$，一般取圆周上两个阵元的间距为 $\lambda /2$，对应的圆半径取为 $R=\frac{\lambda }{4}/\sin (\frac{\pi }{M})$，阵列的第 $m$ 个阵元与 $x$ 轴的角度用 $\frac{2\pi }{M}\cdot (m-1)$ 表示。

以原点为参考点，则位于 $x$ 轴正方向的阵元视为沿着半径方向为轴的参考系。则有时延差：
$${\tau }_{m=1,i}=\frac{R\sin {\phi }_i\cos {\theta }_i}{v}$$
其中 ${\theta }_i$ 是入射投影与半径方向的轴之间的夹角。同理第 $m$ 个阵元即把阵元与原点之间的半径方向作为参考轴，此时入射投影与半径方向的轴之间的夹角为 ${\theta }_i-\frac{2\pi (m-1)}{M}$。则有时延差：
$${\tau }_{m,i}=\frac{R\sin {\phi }_i\cos \left({\theta }_i-\frac{2\pi (m-1)}{M}\right)}{v}$$
于是阵列方向向量为
$$\begin{array}{rl}A& =\left[a({\theta }_1,{\phi }_1),\cdots ,a({\theta }_K,{\phi }_K)\right]\\ a({\theta }_i,{\phi }_i)& =\left[\begin{array}{c}\exp \left(-j2\pi R\sin {\phi }_i\cos \left({\theta }_i\right)\right)\\ \exp \left(-j2\pi R\sin {\phi }_i\cos \left({\theta }_i-\frac{2\pi }{M}\right)\right)\\ ⋮\\ \exp \left(-j2\pi R\sin {\phi }_i\cos \left({\theta }_i-\frac{2\pi (M-1)}{M}\right)\right)\end{array}\right]\end{array}$$

3.5 任意三维阵列天线模型

如图所示为任意阵列模型。假设 $M$ 个阵元任意分布在空间中，二维入射角为 $({\theta }_i,{\phi }_i)$，假设有 $K$ 个信号源。则入射的方向矢量为：
$$\mathbf{V}={\left[\begin{array}{ccc}\sin {\phi }_i\cos {\theta }_i& \sin {\phi }_i\sin {\theta }_i& \cos {\phi }_i\end{array}\right]}^T$$
如果第 $m$ 个阵元的坐标位置为 ${\mathbf{r}}_m=(x_m,y_m,z_m)$ ，波速为 $v$。那么

$${\tau }_m({\theta }_i,{\phi }_i)=\frac{1}{v}\left(x_m\sin {\phi }_i\cos {\theta }_i+y_m\sin {\phi }_i\sin {\theta }_i+z_m\cos {\phi }_i\right)$$

因此可得导向矢量：
$$\mathbf{a}({\theta }_i,{\phi }_i)=\left[\begin{array}{c}{p}_1({\theta }_i,{\phi }_i)\exp (j\frac{2\pi }{\lambda }{\mathbf{r}}_1\cdot \mathbf{V})\\ p_2({\theta }_i,{\phi }_i)\exp (j\frac{2\pi }{\lambda }{\mathbf{r}}_2\cdot \mathbf{V})\\ ⋮\\ p_M({\theta }_i,{\phi }_i)\exp (j\frac{2\pi }{\lambda }{\mathbf{r}}_M\cdot \mathbf{V})\end{array}\right]$$
对于传统阵列，通常极化矢量 $p_k(\theta ,\phi )=1$ 一般省略。但是对于极化敏感阵列，不同载体的影响产生屏蔽效应，所以 $p_k(\theta ,\phi )$ 不能省略。

定义传统阵列方向矩阵为：$A=\left[a({\theta }_1,{\phi }_1),\cdots ,a({\theta }_K,{\phi }_K)\right]$。

参考文献

---

1. 毫米波低复杂度 DOA 估计与波束成形技术的研究 ↩︎