专利名称：实波束形成器的设计方法
技术领域：
本发明属信号处理、涉及阵列信号波束形成技术，尤其涉及高分辨率 实波束形成器的设计方法。
技术背景在阵列信号波束形成技术中，其主动/被动测向系统为了利用孔径有限 的阵列接收的信号进行多目标测向，需要使用高分辨率技术将信号在方向上分开。但是，经典的空间傅里叶(Fourier)变换方法的方向分辨率受瑞 利限(Rayleigh Limit)的限制，而且瑞利限仅由阵列孔径的尺寸决定， 无法通过增加阵列接收信号的观测时间或提高信噪比的方法突破瑞利限。 在现有阵列测向技术中，通常采用卡彭(Capon)的波束形成方法和典型的 子空间方法，如多信号分类方法(MUSIC, Multiple Signals Classification)实现高分辨率测向。这类方法虽然能够突石皮瑞利限对方 向分辨率的限制，但是在数据量较少、信噪比较低的情况下，仍然无法分 开方向上靠得很近的信号，并使测向精度产生较大程度的恶化。在阵列信号处理技术广泛应用的今天，高分辨率的波束形成器通常是 采用线性约束下的最小方差来实现。在卡彭的波束形成方法中，采用线性 约束保证一个方向的信号无失真通过，采用最小方差优化技术最大程度地 抑制其它方向的干拔zf言号。采用这种设计方法，虽然波束形成器的分辨率有提高。但是，由于理 想的波束形成器与自相关矩阵的真实值有关，当阵列快摄数较小、信噪比 较低时，样本自相关矩阵将不可避免地受噪声的影响较大，导致波束形成 器的方向分辨率恶化。因此，在阵列孔径较小的前提下，尽可能降低对阵列快摄数和信噪比的要求，是有效提高波束形成器性能的关键。 发明内容本发明的目的是针对现有技术的不足，提供一种适用于阵列孔径和阵 列快摄数较小、信噪比较低情况的高分辨率波束形成器的设计方法。本发明的目的是这样达到的设计出实波束形成器的实数权向量，利 用两种不同的正交关系，使实波束形成器的输出在对应的信号方向上出现 明显的局部最大值，从而完成旁瓣水平较低、对弱信号和相干信号具有较 高方向分辨率的实波束形成。其步骤是首先，由阵列接收数据确定样本自相关矩阵及其逆矩阵；然后进行第一次正交确定对应搜索方向的酉变换矩阵，对样本自相 关矩阵的逆矩阵进行与搜索方向有关的酉变换，并只取酉变换结果的实数 部分；根据酉变换结果的实数部分的矩阵的列向量之和，该向量在对应的 信号方向上出现较小的和值，得到实波束形成器的权向量。再进行第二种正交确定样本自相关矩阵的实数部分； 实波束形成器的权向量作用于样本自相关矩阵的实数部分，确定实波束形 成器的空间语。最后，根据空间语估计，通过搜索空间语的峰值位置得到对应信号的 方向估计，做出波束形成。所述由阵列接收数据确定的样本自相关矩阵为其中，X(z)为阵列在,时刻的快摄向量，[f表示向量共轭转置，AT为 快摄次数。所述通过搜索空间语的峰值位置得到对应信号的方向估计，是指实波 束形成器的实数权向量对样本自相关矩阵的实数部分做波束形成的输出结果乘以任何非负数时的情况。当其输出结果乘以任何非正值时，则是搜索 空间谱估计的谷值位置得到对应信号的方向。
对样本自相关矩阵进行与搜索方向有关的酉变换矩阵为 0(6>) = diag(fl(^))其中，d(0为阵列的方向向量；
2(0 = 0'(^ -^p的其中，()'表示共轭变换，(广表示矩阵的逆； 只取酉变换结果的实数部分，得到的实数矩阵是A(e)-Re(e^)) 其中，Re()表示实数部分。
所述根据酉变换结果的实数部分的矩阵的列向量之和确定待设计实波束 形成器的实数权向量为
其中，厶为元素都等于1的M维列向量。实数权向量对样本自相关矩 阵的实数部分做波束形成，其实波束形成器的输出结果为
gRBF(^>= r,。、D L、，。、 方程式的右边可以乘以任何非零值。
所述由阵列接收数据确定样本自相关矩阵的形态是均匀线阵列或非均匀线 阵列或圆阵列或随机分布的阵列。
所述对样本自相关矩阵进行与搜索方向有关的酉变换矩阵的形态是均匀线 阵列或非均匀线阵列或圓阵列或随机分布的阵列。
本发明的突出优点是在阵列孔径较小的前提下，能尽量降低对阵列 快摄数和信噪比的要求，设计出高分辨率的实波束形成器。本实波束形成 器的旁瓣水平较低，在低信噪比、多信号、弱信号情况下分辨率高，对相 干信号的分辨率高，同时，实数计算量大幅减少，对阵列信号处理具有重 要意义。将本实波束形成器设计方法应用于阵元数为16、相邻阵元间隔为 半波长的均勻线阵，在快摄数等于64的条件下，即可以清晰的分辨信噪比分别为8. OdB和6. OdB、方向间隔等于2度的两个信号。而在相同条件下 卡彭的波束形成方法和基于正交子空间分析的MUSIC方法已经无法分辨这 两个信号。
将本发明应用于工业化阵列系统的方向高分辨率技术的设计中，以获 得较窄的主瓣和较低的旁瓣水平，从而满足智能天线阵列、声纳阵列、无 线电成像阵列等阵列信号处理系统对波束形成器的性能要求。
图1示出本发明实波束形成器的设计方法流程图； 图2示出本发明的参考阵列结构图3示出1个信号情况下已有的空间谱估计方法和采用本发明方法设 计的实波束形成器之间归一化空间镨估计比较的图；其中"FT"表示空间 傅立叶变换方法，"MUSIC"表示多信号分类方法，"CBF"表示卡彭的波 束方法，"RBF"表示本发明的设计方法，下同。
图4示出2个信号情况下已有的空间谱估计方法和采用本发明的设计 方法设计的实波束形成器之间归 一化空间语估计比较的图。
图5示出4个等功率信号情况下已有的空间谱估计方法和采用本发明 的设计方法设计的实波束形成器之间归一化空间i普估计比较的图。
图6示出4个不等功率信号情况下已有的空间镨估计方法和采用本发 明的设计方法设计的实波束形成器之间归一化空间谱估计比较的图。
图7示出2个相干信号情况下已有的空间i普估计方法和采用本发明的 设计方法设计的实波束形成器之间归 一化空间谱估计比较的图。
具体实施例方式
以下结合附图对本发明的具体实施方式
进行说明。 由阵列接收数据确定样本自相关矩阵，及-丄堂AT(,)AT"(,) 该阵列不受形状的限制，包括均匀线阵列、非均匀线阵列、圆阵列、随^/L分布的阵列 等。以如图2所示的均匀线阵为例子，阵列的阵元数为M,相临阵元之间 的间距为"假设空间有D个非相干的窄带信号A(,), iH7，A…，A到达此 线阵，与线阵的夹角为《。以阵列的第一阵元为参考点，f时刻阵列的输出
信号向量为 AT")-
柳 x2(,)
=J>(《K(0 + v(O = , + v(f)
"1
其中，
*(t) = [4)
1
a)
(2)
(3)
々",
2;r
2"柳
y,(A/-l)rfsin(6>)
(4);
信噪比(dB ): [8. 0]
用最高旁瓣水平(犯)表示各种方法的比较结果是
FT; -12.3; MUSIC: -26.4; CBF: -19.2; RBF : —32.8 。
图4示出了本实波束形成器在较低信噪比情况下的分辨率4^高。
其比较条件是
阵元数16;信号数2; 快摄数64;波达方向(度)[3,5]; 信噪比UB): [8.0,6.0]。
图5示出了实波束形成器在多信号情况下的分辨率较高。 其比较条件是
阵元数32;信号数4; 快摄数64;波达方向(度)[2,4,6,8]; 信噪比(dB): [8.0， 8.0， 8.0, 8.0]。
图6示出了实波束形成器对弱信号的分辨率较高。 其比较条件是
阵元数32;信号数4; 快摄数64;波达方向(度)[2,5， 7， 9];信噪比(dB ): [16.5， 16.5, 6.0， 16.5]。
图7示出了实波束形成器对相干信号的分辨率较高。 其比较条件是
阵元数16;信号数2; 快摄数64;波达方向(度)[2， 6]; 信噪比(dB): [14.0, 14.0]
根据本发明的设计方法进行的实际设计表明，本发明设计的实波束形成 器具有比空间傅里叶(Fourier)变换方法、多信号分类方法(MUSIC)和 卡彭(Capon )的波束形成方法更高的方向分辨率。可以—险证，卡彭(Capon) 的波束形成方法只是利用了第一种正交关系。如图3、图4、图5、图6和 图7的结果所示，与已有的空间傅里叶变换方法(FT)、多信号分类方法 (MUSIC)和卡彭(Capon)的波束形成方法(CBF )相比，利用本发明的设 计方法设计出的实波束形成器(RBF)的旁瓣水平较低、在较低信噪比情况 下的分辨率较高、在多信号情况下的分辨率较高、对弱信号的分辨率较高、 对相干信号的分辨率也较高。在图4、图5、图6和图7所示的情况下，已 有的空间傅里叶变换方法(FT)、多信号分类方法(MUSIC )和卡彭(Capon) 的波束形成方法(CBF)都存在无法分辨的信号，而利用本发明的设计方法 设计出的实波束形成器(RBF)能够清晰的分辨所有信号。
由于都需要计算及-1 ，所以只比较CBF计算(13)式和RBF计算(7) 、 (8)、 (IO)式所需要的计算量，结果如表一所示。
表一:
实k计算量CBFRBF
乘法次数4M2+4M2M2+2M
加法此凄史2M2-22. 5M2-0. 5M-1
可见，本发明设计的实波束形成器(RBF)所需的乘法次凄t比卡彭的波束形成方法少一半。由于本发明设计的实波束形成器不需要先估计信号个数，
明显比多信号分类方法(MUSIC )的计算量少。由于空间傅里叶变换方法(FT) 不适用于高分辨率应用，所以不与该法比较计算量。
虽然已经参考附图对本发明的方法以举例方式进行了描述，但是本发 明不限于上述这些细节，应该包含由本方法衍生的各种变型或改变。
权利要求
1、一种实波束形成器的设计方法，其特征在于设计出实波束形成器的实数权向量，利用两种不同的正交关系，使实波束形成器的输出在对应的信号方向上出现明显的局部最大值，从而完成旁瓣水平较低、对弱信号和相干信号具有较高方向分辨率的实波束形成；其步骤是首先，由阵列接收数据确定样本自相关矩阵及其逆矩阵；然后，进行第一次正交确定对应搜索方向的酉变换矩阵，对样本自相关矩阵的逆矩阵进行与搜索方向有关的酉变换，并只取酉变换结果的实数部分；根据酉变换结果的实数部分的矩阵的列向量之和，该向量在对应的信号方向上出现较小的和值，得到实波束形成器的权向量；再进行第二种正交确定样本自相关矩阵的实数部分；实波束形成器的权向量作用于样本自相关矩阵的实数部分，确定实波束形成器的空间谱；最后，根据空间谱估计，通过搜索空间谱的峰值位置得到对应信号的方向估计，做出波束形成。
2、 根据权利要求1所述的实波束形成器的设计方法，其特征在于 所述由阵列接收数据确定的样本自相关矩阵为其中，AT(0为阵列在^时刻的快摄向量，[f表示向量共轭转置，JV为 快摄次数。
3、 根据权利要求1或2所述的实波束形成器的设计方法，其特征在于 所述通过搜索空间谱的峰值位置得到对应信号的方向估计，是指实波束形成器的实数权向量对样本自相关矩阵的实数部分做波束形成的输出结果乘 以任何非负数时的情况，当其输出结果乘以任何非正值时，则是搜索空间 i瞽估计的谷值位置得到对应信号的方向。
4、 根据权利要求1或2所述的实波束形成器的设计方法，其特征在于 对样本自相关矩阵进行与搜索方向有关的酉变换矩阵为承(0 = diag("的) 其中，fl(0)为阵列的方向向量。对样本自相关矩阵的逆矩阵进行与搜索方向^有关的酉变换结果为^^-^(^r-^W) 其中，()'表示共轭变换，(广表示矩阵的逆； 只取酉变换结果的实数部分，得到的实数矩阵是《(e) = Re^W》 其中，Re()表示实数部分。
5、 根据权利要求1或2所述的实波束形成器的设计方法，其特征在于 所述根据酉变换结果的实数部分的矩阵的列向量之和确定待设计实波束形 成器的实数权向量为其中，/M为元素都等于1的M维列向量，实数权向量对样本自相关矩 阵的实数部分做波束形成，实波束形成器的输出结果为g(0= r//n ) 、 ，n、 方程式的右边可以乘以任何非零值。 < (P)Re(及
6、 根据权利要求2所述的实波束形成器的设计方法，其特征在于所 述由阵列接收数据确定样本自相关矩阵的形态是均匀线阵列或非均匀线阵 列或圆阵列或随^L分布的阵列。
7、根据权利要求4所述的实波束形成器的设计方法，其特征在于所述 对样本自相关矩阵进行与搜索方向有关的酉变换矩阵的形态是均匀线阵列 或非均匀线阵列或圆阵列或随^L分布的阵列。
全文摘要
一种实波束形成器的设计方法。利用两种不同的正交关系设计出实波束形成器的实数权向量，使实波束形成器的输出在对应的信号方向上出现明显的局部最大值，从而完成旁瓣水平较低、对弱信号和相干信号具有较高方向分辨率的实波束形成。步骤是由阵列接收数据确定样本自相关矩阵及其逆矩阵；确定对应搜索方向的酉变换矩阵并进行酉变换，只取酉变换结果的实数部分；根据酉变换结果的实数部分的矩阵的列向量之和得到实波束形成器的权向量。实波束形成器的权向量作用于样本自相关矩阵的实数部分，确定实波束形成器的空间谱通过搜索空间谱的峰值位置得到对应信号的方向估计，完成波束形成。本发明为阵列信号处理提供了全新方法，效果显著，具有重要意义。
文档编号G01S7/02GK101295018SQ20071004893
公开日2008年10月29日 申请日期2007年4月23日 优先权日2007年4月23日
发明者群 万, 杨万麟, 杨瑞明, 沈晓峰, 衡 窦 申请人:电子科技大学