一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法

摘要

发明提供的是一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法。(1)光纤组合阵的三条子阵分别接收入射信号，所述光纤组合阵包含三条子阵1、2、3依次排列，各子阵为均匀等间距直线阵，各子阵内部阵元间距彼此互质，三条子阵彼此间距为接收信号半波长；(2)各子阵所接收的入射信号，利用各互质子阵进行波束形成，获取各扫描方位的波束数据输出；(3)利用各子阵间阵元间距互质的特性，对三个子阵波束域信息进行综合处理，得到抑制栅瓣后的空间谱和目标方位估计结果。本发明能够提高对空间目标方位估计的可靠性和准确度，获得更高的分辨力；相比传统均匀面阵，能够扩大等效孔径，减少阵元数量，具有较强的工程实用价值。

权利要求

1.一种基于光纤组合阵的栅瓣抑制方法，包括如下步骤：(1)光纤组合阵的三条子阵分别接收入射信号，所述光纤组合阵包含三条子阵1、2、3依次排列，各子阵为均匀等间距直线阵，各子阵内部阵元间距彼此互质，三条子阵彼此间距为接收信号半波长；

(2)各子阵所接收的入射信号，利用各互质子阵进行波束形成，获取各扫描方位的波束数据输出；

(3)利用各子阵间阵元间距互质的特性，对三个子阵波束域信息进行综合处理，得到抑制栅瓣后的空间谱和目标方位估计结果；

其特征是步骤(2)具体包括：

所述波束形成的扫描范围为：方位角：θ∈(0°,180°)；俯仰角： 目标方位为 对于每个子阵i＝1,2,3，子阵i在时刻t接收的信号 表示成一个列矢量：式中，N为入射信号数量，ω0为接收信号频率，si(1～N)(t)表示子阵i在t时刻的接收信号复包络， 表示子阵i各阵元在t时刻的噪声， 表示第N个信号到达子阵i各阵元相对于参考阵元的时延；

将上述接收信号列矢量表示为矢量形式如下：

Xi(t)＝AiSi(t)+Ni(t),i＝1,2,3式中，

Si(t)＝[si1(t),si2(t),...,siN(t)]T，Ai为子阵i的Mi×N维流型矩阵，表示为：Ai＝[ai1(ω0) ai2(ω0) … aiN(ω0)],i＝1,2,3其中，导向矢量ai(ω0)为：

将子阵2的第一个阵元定义为参考阵元，位于坐标原点(0,0,0)，空间任意一个阵元的位置为(x,y,z)，由几何关系推导出空间两阵元间的时延差为：将上述公式以 表示，设子阵i的权值为wi，则子阵i的波束输出为：常规波束形成中取 在扫描范围内进行搜索，进而得到空间功率谱，即：其中Rxi为数据协方差矩阵，采用有限快拍数的阵列接收信号对数据协方差矩阵Rxi进行估计，即：其中p为快拍数，根据子阵i的空间谱结果 以及 获得关于目标的方位信息。

2.根据权利要求1所述的基于光纤组合阵的栅瓣抑制方法，其特征是步骤(3)具体包括：首先对各子阵i进行波束输出做乘积处理，在每个扫描方位对不同子阵的波束输出进行互相关运算，并将其作为该方位的空间谱输出，即：其中， 为子阵i,j的波束输出， 为乘积处理空间谱，进行乘积处理后，对各子阵间的互相关输出进行最小处理，即选取各方位上的最小值作为该方位的空间谱结果：之后，按照步骤(2)二对各子阵端阵元组成的半波长均匀线阵进行波束形成，得到其空间功率谱，表示为Pula；对Pmin1,2,3和Pula再进行一步最小处理，得到的综合处理结果表示为PCom，即：PCom＝min(Pmin1,2,3,Pula)。

3.根据权利要求1或2所述的基于光纤组合阵的栅瓣抑制方法，其特征是步骤(1)中：三个子阵的互质因子分别为C1、C2、C3为互为质数的正整数，且C1,C2,C3≥2，子阵间距为λ/2，其中λ为接收信号波长，子阵1阵元数量为M1＝C2*C3，间距为C1*λ/2；子阵2阵元数量为M2＝C1*C3，间距为C2*λ/2；子阵3阵元数量为M3＝C1*C2，间距为C3*λ/2。

说明书

一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法

技术领域

[0001] 本发明涉及的是一种光纤组合阵。本发明申请也涉及一种基于光纤组合阵的栅瓣抑制方法。

背景技术

[0002] 阵列信号处理利用阵列对信号进行空间和时间采样，采用适当方法进行运算，加工和处理，并提取有用信息。阵列信号处理有两大基本任务：一是波束形成，二是空间谱估计。阵列信号处理理论发展至今，出现了如CBF、MVDR、MUSIC、ESPRIT、最大似然法、稀疏重构等算法。上述算法大多基于均匀阵列而言，均匀阵列的内部阵元间距都必须小于或等于入射信号的半波长，导致阵列孔径受到阵元数的制约。对于均匀阵列而言，要获得高分辨率，必须要增加阵元数，由此导致了硬件成本高、阵列设计困难等问题。

[0003] 光纤水听器具有灵敏度高，传输损耗小，信号串扰小等优点，大大降低了阵列的体积和重量，使阵列能够搭载于多类型平台，便于建立大范围、全方位、立体化的水声传感网络，提高对水下空间的实时监视能力。然而，由于光纤水听器自身的独特结构，在进行水声高频信号的方位估计时，无法满足阵列半波长间距布阵的要求，易受到方位估计模糊问题的影响。

发明内容

[0004] 本发明的目的在于提供一种稀疏且可以抑制栅瓣干扰的光纤组合阵。本发明的目的还在于提供一种能够解决利用光纤水听器阵列进行方位估计时由于栅瓣造成的方位模糊问题的基于光纤组合阵的栅瓣抑制方法。

[0005] 本发明的光纤组合阵包含三条子阵1、2、3依次排列，各子阵为均匀等间距直线阵，各子阵内部阵元间距彼此互质，三条子阵彼此间距为接收信号半波长。

[0006] 本发明的光纤组合阵中三个子阵的互质因子分别为C1、C2、C3为互为质数的正整数，且C1,C2,C3≥2，子阵间距为λ/2，其中λ为接收信号波长，子阵1阵元数量为M1＝C2*C3，间距为C1*λ/2；子阵2阵元数量为M2＝C1*C3，间距为C2*λ/2；子阵3阵元数量为M3＝C1*C2，间距为C3*λ/2。

[0007] 本发明的基于光纤组合阵的栅瓣抑制方法包括如下步骤：

[0008] (1)光纤组合阵的三条子阵分别接收入射信号，所述光纤组合阵包含三条子阵1、2、3依次排列，各子阵为均匀等间距直线阵，各子阵内部阵元间距彼此互质，三条子阵彼此间距为接收信号半波长；

[0009] (2)各子阵所接收的入射信号，利用各互质子阵进行波束形成，获取各扫描方位的波束数据输出；

[0010] (3)利用各子阵间阵元间距互质的特性，对三个子阵波束域信息进行综合处理，得到抑制栅瓣后的空间谱和目标方位估计结果。

[0011] 本发明的基于光纤组合阵的栅瓣抑制方法还可以包括：

[0012] 1.步骤(2)具体包括所述波束形成的扫描范围为：方位角：θ∈(0°,180°)；俯仰角：目标方位为 对于每个子阵i＝1,2,3，子阵i在时刻t接收的信号
表示成一个列矢量：

[0013]

[0014] 式中，N为入射信号数量，ω0为接收信号频率，si(1～N)(t)表示子阵i在t时刻的接收信号复包络， 表示子阵i各阵元在t时刻的噪声， 表示第N个信号到达子阵i各阵元相对于参考阵元的时延；

[0015] 将上述接收信号列矢量表示为矢量形式如下：

[0016] Xi(t)＝AiSi(t)+Ni(t),i＝1,2,3

[0017] 式中，

[0018] Si(t)＝[si1(t),si2(t),...,siN(t)]T，Ai为子阵i的Mi×N维流型矩阵，表示为：

[0019] Ai＝[ai1(ω0) ai2(ω0) … aiN(ω0)],i＝1,2,3

[0020] 其中，导向矢量ai(ω0)为：

[0021]

[0022] 将子阵2的第一个阵元定义为参考阵元，位于坐标原点(0,0,0)，空间任意一个阵元的位置为(x,y,z)，由几何关系推导出空间两阵元间的时延差为：

[0023]

[0024] 将上述公式以 表示，设子阵i的权值为wi，则子阵i的波束输出为：

[0025]

[0026] 常规波束形成中取 在扫描范围内进行搜索，进而得到空间功率谱，即：

[0027]

[0028] 其中Rxi为数据协方差矩阵，采用有限快拍数的阵列接收信号对数据协方差矩阵Rxi进行估计，即：

[0029]

[0030] 其中p为快拍数，根据子阵i的空间谱结果 以及 获得关于目标的方位信息。

[0031] 2.步骤(3)具体包括：

[0032] 首先对各子阵i进行波束输出做乘积处理，在每个扫描方位对不同子阵的波束输出进行互相关运算，并将其作为该方位的空间谱输出，即：

[0033]

[0034] 其中， 为子阵i,j的波束输出， 为乘积处理空间谱，

[0035] 进行乘积处理后，对各子阵间的互相关输出进行最小处理，即选取各方位上的最小值作为该方位的空间谱结果：

[0036]

[0037] 之后，按照步骤(2)二对各子阵端阵元组成的半波长均匀线阵进行波束形成，得到其空间功率谱，表示为Pula；对Pmin1,2,3和Pula再进行一步最小处理，得到的综合处理结果表示为PCom，即：

[0038] PCom＝min(Pmin1,2,3,Pula)。

[0039] 针对已有技术中存在的问题，本发明了一种新的稀疏阵列结构及可以抑制栅瓣干扰的信号处理方法。不同于均匀阵列，其阵元间距可以突破半波长的限制。在阵元孔径相等的情况下相比于均匀阵列，稀疏阵列能够降低阵元数量，减小硬件成本。所提出的综合处理方法，不仅能够准确地估计入射信号方位，还能够抑制栅瓣干扰，且谱峰更加尖锐。因此，在利用光纤水听器进行方位估计时，本发明具有较高的现实意义，能够有效的推动相关阵列信号处理算法在现实中的应用，节约成本，易于实现。

[0040] 相比于半波长间距布阵的均匀阵列，本发明利用更少的阵元数，获得了更尖锐的谱峰，并能够消除栅瓣干扰对目标方位估计造成的方位模糊。在采用光纤阵列进行方位估计或在高频背景下，基于光纤组合阵列的栅瓣抑制方法具有非常高的现实意义。

附图说明

[0041] 图1为光纤组合阵阵列模型示意图；

[0042] 图2为空间任意两阵元的几何关系图；

[0043] 图3为乘积处理器结构图；

[0044] 图4为最小处理器结构图；

[0045] 图5为具有栅瓣干扰抑制能力的综合处理器结构图；

[0046] 图6a-图6b为半波长间距均匀平面阵常规波束形成方位谱；

[0047] 图7a-图7b为整阵常规波束形成方位谱；

[0048] 图8a-图8b为基于光纤组合阵的栅瓣抑制综合处理方位谱。

具体实施方式

[0049] 下面举例对本发明做更详细的描述。

[0050] 本发明设计光纤传感器阵列布阵结构，包含三条子阵1、2、3依次排列，子阵为均匀等间距直线阵，且各子阵内部阵元间距彼此互质。三条子阵彼此间距为接收信号半波长，分别接收入射信号。根据各子阵所接收的入射信号，利用各互质子阵进行波束形成(包括CBF、MVDR算法等)，本发明以常规波束形成算法为例，获取各扫描方位的波束数据输出。依据所设计的阵列结构，利用各子阵间阵元间距互质的特性，对三个子阵波束域信息进行综合处理，得到抑制栅瓣后的空间谱和目标方位估计结果。

[0051] (1)结合图1所示，三个子阵的互质因子分别为C1、C2、C3为互为质数的正整数，且C1,C2,C3≥2，子阵间距为λ/2，其中λ为接收信号波长。子阵1阵元数量为M1＝C2*C3，间距为C1*λ/2；子阵2阵元数量为M2＝C1*C3，间距为C2*λ/2；子阵3阵元数量为M3＝C1*C2，间距为C3*λ/2。

[0052] (2)结合图1、图2所示，所述常规波束形成扫描范围为：方位角：θ∈(0°,180°)；俯仰角： 目标方位为 对于每个子阵i＝1,2,3，子阵i在时刻t接收的信号 可以表示成一个列矢量：

[0053]

[0054] 式中，N为入射信号数量，ω0为接收信号频率，si(1～N)(t)表示子阵i在t时刻的接收信号复包络， 表示子阵i各阵元在t时刻的噪声， 表示第N个信号到达子阵i各阵元相对于参考阵元的时延。

[0055] 将上述接收信号列矢量表示为矢量形式如下：

[0056] Xi(t)＝AiSi(t)+Ni(t),i＝1,2,3

[0057] 式中，

[0058] Si(t)＝[si1(t),si2(t),...,siN(t)]T，Ai为子阵i的Mi×N维流型矩阵，表示为：

[0059] Ai＝[ai1(ω0) ai2(ω0) … aiN(ω0)],i＝1,2,3

[0060] 其中，导向矢量ai(ω0)为：

[0061]

[0062] 将子阵2的第一个阵元定义为参考阵元，位于坐标原点(0,0,0)，空间任意一个阵元的位置为(x,y,z)，由几何关系推导出空间两阵元间的时延差为：

[0063]

[0064] 将上述公式以 表示，设子阵i的权值为wi，则子阵i的波束输出为：

[0065]

[0066] 常规波束形成(CBF)中取 在扫描范围内进行搜索，进而得到空间功率谱，即：

[0067]

[0068] 其中Rxi为数据协方差矩阵。实际应用中，一般采用有限快拍数的阵列接收信号对数据协方差矩阵Rxi进行估计，即：

[0069]

[0070] 其中p为快拍数。可以根据子阵i的空间谱结果 以及 获得关于目标的方位信息。

[0071] (3)结合图1、图2、图3、图4和图5所示，由于空间欠采样，步骤(2)中的每个子阵利用常规波束形成扫描得到的空间谱均受到栅瓣干扰的影响，严重影响信号方位估计效果。步骤(3)利用各子阵间阵元间距互质的特性，对三个子阵获取的目标方位信息进行综合处理。利用各子阵间阵元间距互质的特性，对三个子阵获取的目标方位信息进行综合处理。

[0072] 首先对各子阵i进行波束输出做乘积处理，在每个扫描方位对不同子阵的波束输出进行互相关运算，并将其作为该方位的空间谱输出，即：

[0073]

[0074] 其中， 为子阵i,j的波束输出， 为乘积处理空间谱。通过乘积处理在栅瓣位置处形成一定抑制能力，但与半波长间距布阵的均匀阵相比，所得空间谱在非目标方位仍具有较大起伏。

[0075] 进行乘积处理后，对各子阵间的互相关输出进行最小处理，即选取各方位上的最小值作为该方位的空间谱结果：

[0076]

[0077] 进行最小处理后，对栅瓣位置处的干扰抑制能力进一步降低。之后，按照步骤二对各子阵端阵元组成的半波长均匀线阵进行常规波束形成，得到其空间功率谱，表示为Pula。对Pmin1,2,3和Pula再进行一步最小处理，得到的综合处理结果表示为PCom，即：

[0078] PCom＝min(Pmin1,2,3,Pula)

[0079] 所提出的综合处理兼具乘积处理和最小处理的优点，成功地消除了栅瓣对目标方位估计造成的方位模糊，能准确估计到目标的方位。

[0080] 本发明的基于光纤组合阵的栅瓣抑制方法主要实现步骤包括：

[0081] (1)设计光纤传感器阵列布阵结构，包含三条子阵1、2、3依次排列，子阵为均匀等间距直线阵，且各子阵内部阵元间距彼此互质。三条子阵彼此间距为接收信号半波长，分别接收入射信号。

[0082] (2)根据各子阵所接收的入射信号，利用各互质子阵进行波束形成(包括CBF、MVDR算法等)，本发明以常规波束形成算法为例，获取各扫描方位的波束数据输出。

[0083] (3)依据所设计的阵列结构，利用各子阵间阵元间距互质的特性，对三个子阵波束域信息进行综合处理，得到抑制栅瓣后的空间谱和目标方位估计结果。

[0084] 下面通过仿真来验证本发明的效果。仿真条件与结果如下：

[0085] 结合图1所示，所构成光纤组合阵列的3个均匀稀疏互质线阵的阵元数分别为12、15、20，目标信号中心频率设为1500Hz，声速设为1500m/s，对应的阵元间距分别为2.5m(5倍半波长)、2.0m(4倍半波长)、1.5m(3倍半波长)，信噪比SNR＝0dB。采用常规波束形成方法，扫描角度范围为θ∈(0°,180°)， 目标方位角设为

[0086] 结合图6a、图6b、图7a、图7b、图8a、图8b所示，分别为与光纤组合阵相同孔径的均匀半波长间距面阵DOA估计得到的二维方位谱、光纤组合阵整体DOA估计得到的二维方位谱以及利用基于光纤组合阵的栅瓣抑制方法得到的二维方位谱。与相同孔径均匀半波长间距面阵相比，本发明方法所估计的目标方位谱峰更窄，能够达到较好的目标方位估计效果；与光纤组合阵整体DOA估计相比，本发明方法对栅瓣干扰的抑制作用显著，成功地消除了栅瓣对目标方位估计造成的方位模糊，能准确估计到目标的真实方位。