《阵列信号处理及MATLAB实现》信源和噪声模型、阵列天线统计模型2.5 信源和噪声模型2.6 阵列天线的统计模型

68 阅读 0 评论 45 点赞

我是靠谱客的博主不安酸奶，最近开发中收集的这篇文章主要介绍《阵列信号处理及MATLAB实现》信源和噪声模型、阵列天线统计模型2.5 信源和噪声模型2.6 阵列天线的统计模型，觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

PS：文章内容为本人读书笔记，如想阅读更详细内容请购买正版书籍

2.5 信源和噪声模型

2.5.1 窄带信号

如果信号带宽远小于其中心频率，则该信号为窄带信号，即：

$W_{B}/f_0< 1/10$

其中， $W_B$ 为信号带宽， $f_0$ 为中心频率。通常将正弦信号和余弦信号统称为正弦型信号，正弦型信号是典型的窄带信号。若无特殊说明，本书中所提及的窄带信号表示为：

$s(t)=a(t)e^{j[w_0t+theta (t)]}$

式中， $a(t)$ 为慢变幅度调制函数（或称实包络）， $theta(t)$ 为慢变相位调制函数， $w_0=2pi f_0$ 为载频。一般情况下， $a(t)$ 和 $theta (t)$ 包含了全部有用信息。

2.5.2 相关系数

对于接受到的多个信号，一般可利用相关系数（或称互相关系数）来衡量信号之间的关联程度。对于两个平稳信号 $s_i(t)$ 和 $s_j(t)$ ，其相关系数定义为：

$rho _{ij}=frac{E[(s_i(t)-E[s_i(t)])(s_j(t)-E[s_j(t)])]}{sqrt{E[s_i(t)-E[s_i(t)]]^2E[s_j(t)-E[s_j(t)]]^2}}$

当 $rho _{ij}=0$ 时，称 $s_i(t)$ 和 $s_j(t)$ 不相关（或称不相干）；当 $0<left | rho _{ij} right |<1$ 时，称 $s_i(t)$ 和 $s_j(t)$ 部分相关；当 $rho _{ij}=1$ 时，称 $s_i(t)$ 和 $s_j(t)$ 完全相关（或称相干）

2.5.3 噪声模型

本书中，无特殊说明，阵元接收到的噪声均假设为平稳零均值高斯白噪声，方差为 $sigma^2$ 。各阵元间噪声互不相关，且与目标源不相关。这样，噪声向量 $vec{n}(t)$ 的二阶矩满足：

$Ebegin{Bmatrix} vec n(t_1),vec n^H(t_2) end{Bmatrix}=sigma ^2I delta_{t_1,t_2}$

$Ebegin{Bmatrix} vec n(t_1),vec n^T(t_2) end{Bmatrix}=0$

2.6 阵列天线的统计模型

2.6.1 前提及假设

关于接收天线阵的假设：接收阵列由位于空间已知坐标处的无源阵元按一定的形式排列而成。假设阵元的接收特性仅与位置有关而与尺寸无关（认为是一个点），并且阵元都是全向阵元，增益均相等，相互之间互耦忽略不计。阵元接收信号时将产生噪声，假设其为加性高斯白噪声，各阵元上的噪声相互统计独立，且噪声与信号是统计独立的。

关于空间源信号的假设：假设空间信号的传播介质是均匀且各向同性的，这是空间信号在介质中按直线传播；同时假设阵列处于空间信号辐射的远场中，所以空间源信号到达阵列时可以看做平行平面波，则空间源信号到达阵列各阵元在时间上的不同时延，可由阵列的几何结构和空间波来向决定。空间波的来向在三维空间中常用仰角和方位角来表征。

2.6.2 阵列的基本概念

令信号载波为 $e^{jwt}$ ，并以平面波形式在空间沿波数向量 $vec k$ 的方向传播，设基准点处的信号为 $s(t)e^{jwt}$ ，则距离基准点r处的阵元接收信号为：

$s_r(t)=s(t-frac{1}{c}vec r ^T vec alpha)exp[j(wt-vec r^Tvec k)]$

式中， $vec k$ 为波数向量； $vec{alpha}=vec{k}/ left | vec kright |$ 为电磁波传播方向单位向量； $left | vec k right |=w/c=2pi /lambda$ 为波数，其中c为光速， $lambda$ 为电磁波的波长； $frac{1}{c}vec r ^T vec alpha$ 为信号相对于基准点的延迟时间； $vec r^Tvec k$ 为电磁波传播到离基准点r处的阵元相对于电磁波传播到基准点的滞后相位。 $theta$ 为电磁波传播方向角，他是相对于x轴的逆时针旋转方向定义的，显然，波数向量可表示为：

$vec k =k[costheta,sintheta]^T$

电磁波从点辐射源以球面波向外传播，只要距离足够远，在接收的局部区域，球面波可以近似为平面波。

设在空间有M个阵元组成阵列，阵元从1到M编号，并以阵元1（也可以选择其他阵元）作为基准或参考点。设各阵元无方向性（全向），相对于基准点的位置向量分别为 $r_i(i=1,...,M;r_1=0)$ 。若基准点接收信号为 $s(t)e^{jwt}$ ，则各阵元上的接收信号分别为：

$s_r(t)=s(t-frac{1}{c}vec r_i ^T vec alpha)exp[j(wt-vec r_i^Tvec k)]$

在通信中，信号的频带B比载波值w小得多，所以s(t)的变化相对缓慢，延时 $frac{1}{c}vec r_i ^T vec alpha<<(1/B)$ ，故有： $s(t-frac{1}{c}vec r_i ^T vec alpha)approx s(t)$ ，即信号包络在各阵元上的差异可忽略，称为窄带信号。

此外，阵列信号总是变换到基带再进行处理，因而可将阵列信号用向量形式表示为：

$s(t)=[s_1(t),s_2(t),...,s_M(t)]^T=s(t)[e^{-jvec r_1^T},e^{-jvec r_2^T},...,e^{-jvec r_M^T}]^T$

上式中的向量部分称为方向向量，因为当波长和阵列的几何结构确定时，该向量只与到达波的空间角向量 $theta$ 有关，方向向量记作 $vec a(theta)$ 。例如，若选择第一个阵元为基准点，则方向向量为：

$vec a(theta)=[1,e^{-jbar r_2vec k},...,e^{-jbar r_Mvec k}]$

公式中： $bar{r_i}=vec r_i-vec r_1(i=2,...,M)$

实际使用的阵列结构要求方向向量 $vec a(theta)$ 必须与空间角向量 $theta$ 一一对应，不能出现模糊现象。当有多个（例如K个）信源时，到达波的方向向量可分别用 $vec a(theta)$ 表示。这K个方向向量组成的矩阵 $A=[vec a(theta_1),vec a (theta_2),...,vec a(theta_K)]$ 称为阵列的方向矩阵或响应矩阵，它表示所有信号源的方向。

改变空间角 $theta$ ，使方向向量 $vec a (theta)$ 在M维空间内扫描，所形成的曲面称为阵列流行。

阵列流行常用符号A表示，即有：

$A=begin{Bmatrix} vec a(theta|theta in Theta ) end{Bmatrix}$

其中， $Theta=[0,2pi)$ 是波达方向 $theta$ 所有可能取值的集合。因此阵列流行A是阵列方向向量的集合。阵列流行A包含了阵列几何结构、阵元模式、阵元间的耦合、频率等影响。

2.6.3 天线阵模型

假设有一个天线阵列，由M具有任意方向性的阵元按任意排列构成。同时设有K个具有相同中心频率 $w_0$ 、波长 $lambda$ 的空间窄带平面波（M>K）分别以来向角 $Theta_1,Theta_2,...,Theta_K$ 入射到该阵列。如下图所示。

其中 $Theta_i=(theta_i,phi_i),i=1,2,...K$ 。 theta_i 和 $phi_i$ 分别是第i个入射信号的仰角和方位角， $0leq theta_i <90^{circ},0leq phi_i<360^{circ}$

则阵列第m个阵元的输出可以表示为：

$x_m(t)=sum_{i=1}^{K}s_i(t)e^{jw_0tau_m(Theta_i) }+n_m(t)$

其中， $s_i(t)$ 为入射到阵列的第i个源信号， $n_m(t)$ 为第m个阵元的加性噪声， $tau_m(Theta_i)$ 为来自 $Theta_i$ 方向的源信号投射到第m个阵元时，相对于选定参考点的时延。并记：

$X(t)=[x_1(t),x_2(t),...,x_M(t)]^T$

$N(t)=[n_1(t),n_2(t),...,n_M(t)]^T$

另外，S(t)为Kx1维列向量

$S(t)=[s_1(t),s_2(t),...,s_K(t)]$

$A(Theta)$ 为MxK矩阵的方向矩阵： $A(Theta)=[vec a(Theta_1),vec a(Theta_2),...,vec a(Theta_K)]$

矩阵 $A(Theta)$ 中任一列向量 $vec a(Theta_i)$ 是阵列在空间源信号中一个来向为 $Theta_i$ 的方向向量，且是Mx1维列向量： $vec a(Theta)=[e^{jw_0tau_1(Theta_1)},e^{jw_0tau_2(Theta_2)},...,e^{jw_0tau_M(Theta_M)}]^T$

因此，如果用矩阵描述，即使在最一般化的情况下，阵列信号模型可简单的表示为：

$X(t)=A(Theta)S(t)+N(t)$

很显然，方向矩阵 $A(Theta)$ 与矩阵的形状、信号源的来向有关，而一般在实际应用中，天线阵的形状一旦固定就不会再改变。所以，矩阵 $A(Theta)$ 中任意一列总是和某个空间源信号的来向是紧密联系着的。

2.6.4 阵列的方向图

阵列输出的绝对值与来波方向之间的关系称为天线的方向图。方向图一般有两类，一类是阵列输出的直接相加（不考虑信号及其来向），即静态方向图；另一类是带指向的方向图（考虑信号指向），当然信号的指向是通过控制加权的相位来实现的。从前面的信号模型可知，对于某一确定的 $M$ 元空间阵列，在忽略噪声的条件下，第 $l$ 个阵元的复振幅为：

$x_l=g_0e^{-jwtau_l},l=1,2,...,M$

式中， $g_0$ 为来波的复振幅， $tau_l$ 为第 $l$ 个阵元与参考点之间的延迟。设第 $l$ 个阵元的权值为 $w_l$ ，那么所有阵元加权的输出结果为：

$Y_0=sum_{l=1}^{M}w_lg_0e^{-jwtau_l},l=1,2,...,M$

对上式取绝对值并归一化之后可得到空间阵列的方向图 $G(theta)$ 为：

$G(theta)=frac{begin{vmatrix} Y_0 end{vmatrix}}{ max(begin{vmatrix} Y_0 end{vmatrix}) }$

如果式中 $w_l=1,l=1,2,...,M$ 上式则成为静态方向图 $G(theta)$

下面考虑均匀线列阵的方向图。假设均匀线列阵的间距为 $d$ ，且以最左边的阵元为参考点；另假设信号入射方位角为 $theta$ ，其中方位角表示信号入射方向与线阵法线方向的夹角。则第 $l$ 个阵元与参考点的波程差(时延)为 $tau_l=(x_ksintheta)/c=(l-1)dsintheta/c$ ，则阵列的输出为：

$Y_0=sum_{l=1}^{M}w_lg_0e^{-jwtau_l}=sum_{l=1}^Mw_lg_0e^{-jw(frac{(l-1)dsintheta}{c})}=sum_{l=1}^Mw_lg_0e^{-jfrac{w}{c}(l-1)dsintheta}=sum_{l=1}^Mw_lg_0e^{-jfrac{2pi}{lambda}(l-1)dsintheta}=sum_{l=1}^Mw_lg_0e^{-j(l-1)beta}$

上式中， $beta=2pi d sintheta/lambda$ , $w=2pi f$ , $lambda$ 为入射信号的波长。

当上式中 $w_l=1,l=1,2,...,M$ 时，公式可以进一步化简为：

$Y_0=Mg_0e^{j(M-1)beta /2}frac{sin(Mbeta/2)}{Msin(beta/2)}$

可得均匀线阵的静态方向图为（公式2.6.17）：

$G_0(theta)= |frac{sin(Mbeta/2)}{Msin(beta/2)}|$

当上式中 $w_l=e^{j(l-1)beta_d},beta_d = 2pi d sintheta_d/lambda,l=1,2,...,M$ 时，公式可以简化为：

$Y_0=Mg_0e^{j(M-1)(beta-beta_d)/2}frac{sin[M(beta-beta_d)/2]}{Msin[(beta-beta_d)/2]}$

于是可得指向为 theta_d 的阵列方向图为（公式2.6.19）：

$G(theta)=| frac{sin[M(beta-beta_d)/2]}{Msin[(beta-beta_d)/2]}|$

2.6.5 波束宽度

线阵的测向范围为[-90°，90°]，而一般的面阵（如圆阵）的测向范围为[-180°，180°]。为了说明波束宽度，下面只考虑线阵。

由公式2.6.17可知，M个阵元的均匀线阵的静态方向图为：

$G_0(theta)= |frac{sin(Mbeta/2)}{Msin(beta/2)}|$

式中，空间频率 $beta = (2pi d sintheta)/lambda$

这对于天线静态方向图主瓣的零点，由 $|G_0(theta)|^2=0$ 可得零点波束宽度 $BW_0$ 为：

$BW_0 = 2$arcsin$(lambda/Md)$

而由 $|G_0(theta)|^2=1/2$ 可得到半功率点波束宽度 $BW_{0.5}$ ，在 Mdgg lambda 的条件下有：

$BW_{0.5}approx 0.886lambda/Md$

本书中一般考虑的是静态方向图的半功率点波束宽度，即对于均匀线阵而言，其波束宽度为：

$BW_{0.5}approx frac{51^{circ}}{D/lambda}=frac{0.89}{Dlambda}rad$

式中，D为天线的有效孔径， $lambda$ 为信号的波长，rad表示弧度单位。对于M阵元的等距均匀线阵，阵元间距为 $lambda/2$ ，则天线的有效孔径为 $D=(M-1)lambda/2$ ，所以对于ULA阵，阵列波束宽度的近似计算公式为：

$BWapprox 102^{circ}/M$

关于波束宽度，以下几点需要注意。

（1）波束宽度与天线孔径成反比，一般情况下天线的半功率点波束宽度与天线孔径之间关系为：

$BW_{0.5}=(40sim 60)frac{lambda}{D}$

（2）对于某些阵列（如线阵），天线的波束宽度与波束指向有关系，如波束指向为 theta_d 时，均匀线阵的波束宽度为：

$BW_0=2$arcsin$(frac{lambda}{Md}+sintheta_d)$

$BW_{0.5}approx 0.886frac{lambda}{Md}frac{1}{costheta_d}$

（3）波束宽度越窄，阵列的指向性越好，也就说明阵列分辨空间信号的能力越强。

2.6.6 分辨率

阵列测向中，在某方向上对信源的分辨率与在该方向附近阵列方向矢量的变化率直接相关。在方向矢量变化较快的方向附近，随信源角度变化阵列快拍数据变化也越大，相应的分辨率也高。定义一个表征分辨率 $D(theta)= ||frac{dvec a(theta)}{dtheta}||propto ||frac{dtau}{dtheta} ||$

$D(theta)$ 越大表示在该方向的分辨率越高。

对于均匀线阵，有 $D(theta)propto costheta$ ，说明信号在0°方向分辨率最高，而在60°方向分辨率已经降了一半，所以一般线阵的测向范围为-60°到60°。