BMP
基本步骤:首先估计空间中主瓣干扰的入射角度(需要精确估计入射角,但在工程中很难实现),再根据估计出的角度来构造阻塞矩阵B,接着对雷达接收到的数据进行变换处理以消除主瓣干扰分量,最后再做一般自适应波束形成抑制旁瓣干扰。(BMP 方法将会损失一个有效的阵列单元自由度)
详细步骤:
- 实际数据的协方差矩阵: R x = E [ X X H ] = A R s A H + σ n 2 I R_{x}=E[XX^{H}]=AR_{s}A^{H}+\sigma _{n}^{2}I Rx=E[XXH]=ARsAH+σn2I一般无法得到,所以我们采用样本协方差矩阵 R ^ x = 1 N X X H \hat{R}_{x}=\frac{1}{N}XX^{H} R^x=N1XXH来作为实际协方差矩阵的估计。阵列接收信号中干扰信号的强度比噪声信号的强度大很多,并且假设在主瓣内只有一个干扰信号,所以我们可以选取比较简单的Capon 法来估计,由于仅仅需要估计主瓣干扰的方向,仅需要在方向图的主瓣内进行角度搜索即可,采用Capon估计主瓣干扰的入射角度 α 1 \alpha_{1} α1: α 1 = max α , a l p h a ∈ m a i n l o b e 1 a H ( α ) R ^ x − 1 a ( α ) \alpha_{1}=\max \limits_{\alpha,alpha\in mainlobe}\frac{1}{a^{H}(\alpha)\hat{R}{_{x}}^{-1}a(\alpha)} α1=α,alpha∈mainlobemaxaH(α)R^x−1a(α)1。
- 用估计得到的
α
1
\alpha_{1}
α1构造阻塞矩阵:
B = [ 1 − e − j u 1 0 ⋯ 0 0 0 1 − e − j u 1 ⋯ 0 0 ⋮ ⋮ ⋮ ⋮ 0 0 ⋯ 1 − e − j u 1 0 0 ⋯ 0 1 − e − j u 1 ] B=\begin{bmatrix} 1& -e^{-ju_{1}}& 0& \cdots & 0& 0\\ 0& 1& -e^{-ju_{1}}& \cdots & 0& 0\\ \vdots & \vdots & & & \vdots & \vdots \\ 0& 0& \cdots & 1& -e^{-ju_{1}}& \\ 0& 0& \cdots & 0& 1& -e^{-ju_{1}} \end{bmatrix} B=⎣⎢⎢⎢⎢⎢⎡10⋮00−e−ju11⋮000−e−ju1⋯⋯⋯⋯1000⋮−e−ju1100⋮−e−ju1⎦⎥⎥⎥⎥⎥⎤
其中 u 1 = 2 π ( d / λ ) s i n α 1 u_{1}=2\pi (d/\lambda) sin\alpha _{1} u1=2π(d/λ)sinα1, α 1 \alpha _{1} α1为主瓣干扰的入射角,这里的阻塞矩阵是 ( N − 1 ) × N (N-1)\times N (N−1)×N的。 N N N为阵元数,通过这个方式,可以在不损失期望信号的前提下由邻近的传感器消除主瓣干扰。由于引入阻塞矩阵会导致一个自由度的损失,但是预处理数据 X B l o c k X_{Block} XBlock中不再含有主瓣干扰。 - 对数据进行预处理:
X
B
l
o
c
k
=
B
X
X_{Block}=BX
XBlock=BX
可以得到预处理之后的数据的协方差矩阵的估计: R ^ y = B R ^ x B H \hat{R}_{y}=B\hat{R}_{x}B^{H} R^y=BR^xBH,这里的 R ^ y \hat{R}_{y} R^y为 ( N − 1 ) × ( N − 1 ) (N-1)\times (N-1) (N−1)×(N−1)维。 - 进行常规的自适应波束形成:构造权矢量
w
a
d
p
t
=
R
^
x
−
1
a
ˉ
(
α
0
)
a
ˉ
H
(
α
0
)
R
^
x
−
1
a
ˉ
(
α
0
)
w_{adpt}=\frac{\hat{R}{_{x}}^{-1}\bar{a}(\alpha_{0})}{\bar{a}^{H}(\alpha_{0})\hat{R}{_{x}}^{-1}\bar{a}(\alpha_{0})}
wadpt=aˉH(α0)R^x−1aˉ(α0)R^x−1aˉ(α0)
其中, a ˉ ( α 0 ) \bar{a}(\alpha_{0}) aˉ(α0)为 ( N − 1 ) × 1 (N-1)\times 1 (N−1)×1维。