数字多波束相控阵基础知识

均匀线阵

对于一个均匀线阵，假设每个阵元都是理想全向点源，布局如上图所示， 第 1 个阵元位于坐标原点，如果阵元间距为d ，空间窄带远场平面波到达原点的信号可以表示为
$$s_1(t)=s(t)=a(t)coswt-b(t)sinwt$$
第$l$个阵元接收信号表示为
$$s_l(t)=s(t-{\tau }_l)$$
由于信号是窄带信号，相对于载波相位变化，基带信号相位变化可以忽略不计，则
$$s_l(t)=a(t)cosw(t-{\tau }_l)-b(t)sinw(t-{\tau }_l)$$
第$l$个阵元接收信号相对于第1个阵元接收信号的延时
$${\tau }_l=-(l-1)\frac{dsin\theta }{c}$$
注意${\tau }_l$符号为负，因为第$l$个阵元接收信号的等相位面落后于第1个阵元接收信号，故需要加上正的延时。
其中$\theta$为信号入射角度，范围为$[-\frac{\pi }{2},\frac{\pi }{2}]$。
不同的阵元接收信号，延时不同，信号间存在由延时引起的相位差，该相位差的大小等于
$$w{\tau }_l=w(l-1)\frac{dsin\theta }{c}=\frac{2\pi (l-1)dsin\theta }{\lambda }$$
常取阵元间距$d=\frac{\lambda }{2}$，则$w{\tau }_l=\pi (l-1)sin\theta$

相控阵基本结构图

相控阵天线的基本原理，是通过对每个阵元接收信号进行加权处理，校正信号由于波程差所导致的相位差，该过程如上图所示 。
如果相控阵采用模拟方式实现，则加权移相是一个移相器。当相控阵采用数字化方式实现时，移相功能则通过一个复数乘法器实现。由于一个实数乘法器无法完成移相功能，所以，在数字相控阵当中，需要首先将阵元接收的实信号，
通过一定的正交变换，转换为复信号，也就是该实信号的解析信号。完成正交化后的复信号，通过复数加权，就可以实现移相功能。
常用的相控阵信号表示方式是采用复数矢量的方式，此时的阵列接收信号为一个列矢量，每个元素对应一个阵元接收信号，并且信号也是解析表示。此时第l 个阵元接收解析信号可以表示为
$$x_l(t)=s(t)e^{-jw{\tau }_l}$$
则，阵列接收信号矢量表示为
$$\mathbf{x}(t)=\left[\begin{array}{c}{x}_1(t)\\ ...\\ x_L(t)\end{array}\right]=s(t)\left[\begin{array}{c}{e}^{-jw{\tau }_1}\\ ...\\ e^{-jw{\tau }_L}\end{array}\right]+\mathbf{n}(t)=s(t)\mathbf{v}(\theta )+\mathbf{n}(t)$$
其中$s(t)=[a(t)+jb(t)]e^{jwt}=m(t)e^{jwt}$，$m(t)$称为信号复包络或复基带，信号为窄带信号时，$m(t-\tau )\approx m(t)$。
噪声列矢量为$\mathbf{n}(t)$为
$$\mathbf{n}(t)=\left[\begin{array}{c}{n}_1(t)\\ ...\\ n_L(t)\end{array}\right]$$
而矢量$\mathbf{v}(\theta )$称为信号的方向矢量，或导引矢量，它是由信号到达方向决定的列矢量，即信号入射角度$\theta$。
$$\mathbf{v}(\theta )=\left[\begin{array}{c}{e}^{-jw{\tau }_1}\\ ...\\ e^{-jw{\tau }_L}\end{array}\right]$$
定义相控阵的权矢量为
$$\mathbf{w}=\mathbf{v}({\theta }_0)$$
其中${\theta }_0$为期望信号入射角度。
所以，窄带相控阵的权矢量就是信号的方向矢量。相控阵天线加权后的输出信号为
$$y(t)={\mathbf{w}}^H\mathbf{x}(t)=\mathbf{v}({\theta }_0{)}^{H}s(t)\mathbf{v}(\theta )+\mathbf{v}({\theta }_0{)}^H\mathbf{n}(t)=Ls(t)+{\mathbf{w}}^H\mathbf{n}(t)$$
上述加权求和的操作，在相控阵中通常称为波束合成。符号“$H$”为共轭转置操作。故$${\mathbf{w}}^H\mathbf{v}(\theta )={\left[\begin{array}{c}{e}^{-jw{\tau }_1}\\ ...\\ e^{-jw{\tau }_L}\end{array}\right]}^H\left[\begin{array}{c}{e}^{-jw{\tau }_1}\\ ...\\ e^{-jw{\tau }_L}\end{array}\right]=1+...+1=L$$

确定信号功率和随机信号功率

对于确知信号$s(t)$
$$P=\underset{T\to \infty }{\lim }\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}{s}^2(t)dt$$
若$s(t)$为复信号，则
$$P=\underset{T\to \infty }{\lim }\frac{1}{T}{\int }_{-\frac{T}{2}}^{\frac{T}{2}}s(t)s^{\ast }(t)dt$$

对于一个随机过程$\xi (t)$，其均值或数学期望定义
$$E\left[\xi (t)\right]={\int }_{-\infty }^{\infty }x{f}_1(x,t)dx$$
常记为$a(t)$，其中$f_1(x,t)$为$\xi (t)$的一维概率密度函数。
随机过程的方差定义为
$$D\left[\xi (t)\right]=E\left\{{\left[\xi (t)-a(t)\right]}^2\right\}=E\left[{\xi }^2(t)\right]-a^2(t)$$
其常记为${\sigma }^2(t)$
$$D\left[\xi (t)\right]={\int }_{-\infty }^{\infty }{x}^2(t)f_1(x,t)dx-[a(t){]}^2$$
故方差等于均方值与均值平方之差。前者为交流功率，后者为直流功率。通常$x(t)$为交流信号，故其信号功率就是方差。
另外有
$$P={\int }_{-\infty }^{\infty }P(f)df，P(f)为功率谱密度$$

数字波束合成 （Digital Beam Forming，DBF）

相控阵方向图

根据阵列天线方向图乘积定理，方向图等于阵元因子与阵因子的乘积。
$$P(\theta ,\varphi )=P_E(\theta ,\varphi )P_A(\theta ,\varphi )$$
通常假设阵元都是相同的全向天线，故$P_E(\theta ,\varphi )=1$
由输出信号与接收信号的关系为$y(t)={\mathbf{w}}^H\mathbf{x}(t)$
$$y(t)y^{\ast }(t)dt={\mathbf{w}}^{H}s(t)\mathbf{v}(\theta )s^{\ast }(t)[{\mathbf{w}}^H\mathbf{v}(\theta ){]}^{\ast }=|s(t){|}^2\cdot |{\mathbf{w}}^{H}s(t)\mathbf{v}(\theta ){|}^2$$
其中$\ast$对于矩阵或矢量而言即为$H$，反之$H$对于复数而言即是$\ast$。
若$P_i=1$，则输出信号功率表示为
$$P(\theta )={\Vert \begin{array}{c}{\mathbf{w}}^H\mathbf{v}(\theta )\end{array}\Vert }^2$$
则场强方向图表示为
$$B(\theta )=\Vert \begin{array}{c}{\mathbf{w}}^H\mathbf{v}(\theta )\end{array}\Vert$$
式中，符号$\Vert \begin{array}{}\end{array}\Vert$为2-范数。

相控阵的信噪比

输出信号$y(t)=Ls(t)$的功率
$$P_{so}=\frac{1}{T}\int y(t)y^{\ast }(t)dt=\frac{1}{T}\int Ls(t)L{s}^{\ast }(t)dt=L^2\frac{1}{T}\int s(t)s^{\ast }(t)dt=L^2{P}_{si}$$
噪声功率
由$y(t)=)=Ls(t)+{\mathbf{w}}^H\mathbf{n}(t)$，移相不改变噪声功率。
$$P_{no}=L{P}_{ni}$$
其可用等效噪声温度描述。
故
$$SN{R}_o=\frac{P_{so}}{P_{no}}=L\frac{P_{si}}{P_{ni}}$$

相控阵的优点

抗干扰、多波束、宽带。