本发明专利技术涉及声纳信号处理领域，主要是一种Ｕ型阵波束形成加权方法，根据“Ｕ”型阵阵元在与波束预成方向垂直的波阵面方向的投影坐标位置关系确定插值点，通过对窗函数在插值点的插值，计算出与“Ｕ”型阵阵元位置相匹配的修正窗函数权值。本发明专利技术提供了一种新型的，计算方法简单，计算量小，易于实现的“Ｕ”型阵波束形成加权方法，把基于均匀阵的窗函数加权推广到非均匀阵加权，有效地抑制了旁瓣，提高了“Ｕ”型阵波束形成的性能。
【技术实现步骤摘要】
本专利技术涉及声纳信号处理领域，具体涉及阵列信号处理领域，主要是一种U型阵波束形成加权方法。
技术介绍
波束形成是声纳信号处理的一个热点研究问题。无论是被动声纳还是主动声纳都要有波束形成系统。波束形成系统是现代声纳系统的核心部件，是声纳具有良好的战术、技术性能的基础。常规波束形成方法(Coventional beamforming，简称CBF)是阵列信号处理中的经典方法，它的优点是结构简单，易于实现，而且对各种恶劣的环境有较好的适应能力。虽然现在各种高分辨算法不断涌现，但是常规波束形成方法(CBF)仍然是目前工程上运用最广泛的方法。 对各个阵元的输出信号在幅度上乘以一个实数，用这种方法来改善基阵的指向性或抑制某种特定的干扰，这就是基阵的加权。工程中运用最普遍的窗函数加权方法，是基于均匀阵提出的。对于像“U”型阵这样的非均匀阵，窗函数的权值与阵元位置不匹配，因而权值不是最佳的。
技术实现思路
本专利技术所要解决的技术问题是提供一种适合“U”型水听器基阵波束形成的新型加权方法，使窗函数的权值与阵元位置匹配，使“U”型阵获得比常规窗函数加权更低的旁瓣结构，提高波束形成的性能。 为解决以上技术问题，本专利技术是提出以下技术方案实现的 其包括以下步骤 步骤一采用M元“U”型阵接收信号，通过前置预处理机对基阵接收信号进行时间滤波、A/D变换、放大等前置预处理。令阵元编号分别为1，2，…，M，则阵元输出信号为 x(n)＝T (1) 步骤二信号处理机对时间信号进行快速傅立叶变换(FFT)，将信号从时域变换到频域，频域输出为 X(f)＝T (2) 其中xm(f)为第m个阵元在频率f的值，m＝1，2，…，M。设接收信号为带限信号，则fl≤f≤fh，fl，fh分别为信号频率的下限和上限； 步骤三考虑到“U”型阵波束形成时存在的遮挡问题，如果要形成θ方位的波束，需要通过计算斜率为智tgθ的直线ll与“U”型阵的切点位置P，如图2所示。根据P点位置，确定第1，2，…，N号阵元参加波束形成，由(2)式选出这N个阵元的频域输出为 则波束形成的驾驶向量为 W(f，θ)＝(4) 其中wn(f，θ)＝exp(j2πfτn(θ))，τn(θ)为第n号阵元到参考点(一般设为坐标原点0点)的时延补偿，n＝1，2，…，N； 步骤四把“U”型阵中参加波束形成的N个阵元投影到与θ垂直的波阵面(直线qq)方向，得到投影点A1，A2，…，AN。先假定在投影点左右两端最边缘的两点(设为点A1和点AN)之间形成了一个均匀阵，对该均匀阵进行常规窗函数加权，再对形成的“窗口”进行插值。插值点由各投影点在投影连线(线A1AN)的位置确定，投影点An(其中，n＝1，2，...，N)对应的插值点计算公式如下 其中 表示线段A1An的长度， 表示线段A1AN的长度。根据插值点对常规窗函数在范围内进行插值，插值后的值即为窗函数的修正权值。设所得到的修正窗函数权值为 B＝ (6) 其中bn为对应第n号阵元的修正窗函数权值，n＝1，2，…，N； 步骤五结合步骤三和步骤四的结果，进行频域波束形成得到θ方位的频域波束输出为 本专利技术能带来以下有益效果 本专利技术根据“U”型阵阵元在与波束预成方向垂直的波阵面方向的投影坐标位置关系确定插值点，通过对窗函数在插值点的插值，计算出与“U”型阵阵元位置相匹配的修正窗函数权值。本专利技术提供了一种新型的，计算方法简单，计算量小，易于实现的“U”型阵波束形成加权方法，把基于均匀阵的窗函数加权推广到非均匀阵加权，有效地抑制了旁瓣，提高了“U”型阵波束形成的性能。 附图说明 图1本专利技术的波束形成流程框图； 图2“U”型水听器基阵θ方位加权波束形成示意图； 图3等弧长“U”型水听器基阵结构图； 图4阵元随机布放“U”型水听器基阵结构图； 图5等弧长U阵修正的海明加权窗函数与海明窗比较图； 图6阵元随机布放U阵修正的海明加权窗函数与海明窗比较图； 图7等弧长“U”型阵加权波束图； 图8阵元随机布放“U”型阵加权波束图。 具体实施例方式 下面结合具体实施例和附图对本专利技术做进一步的描述 本专利技术的波束形成流程框图如图1所示，这种U型阵波束形成加权方法，具体包括以下步骤 步骤一采用M元“U”型阵接收信号，通过前置预处理机对基阵接收信号进行时间滤波、A/D变换、放大等前置预处理。令阵元编号分别为1，2，…，M，则阵元输出信号为 x(n)＝T(1) 步骤二信号处理机对时间信号进行快速傅立叶变换(FFT)，将信号从时域变换到频域，频域输出为 X(f)＝T(2) 其中xm(f)为第m个阵元在频率f的值，m＝1，2，…，M。设接收信号为带限信号，则fl≤f≤fh，fl，fh分别为信号频率的下限和上限； 步骤三考虑到“U”型阵波束形成时存在的遮挡问题，如果要形成θ方位的波束，需要通过计算斜率为tgθ的直线ll与“U”型阵的切点位置P，本专利技术形成θ方位加权波束示意图如图2所示。根据P点位置，确定第1，2，…，N号阵元参加波束形成，由(2)式选出这N个阵元的频域输出为 则波束形成的驾驶向量为 W(f，θ)＝(4) 其中wn(f，θ)＝exp(j2πfτn(θ))，τn(θ)为第n号阵元到参考点(坐标原点0点)的时延补偿，n＝1，2，…，N； 步骤四把“U”型阵中参加波束形成的N个阵元投影到与θ垂直的波阵面(直线qq)方向，得到投影点A1，A2，…，AN。先假定在投影点左右两端最边缘的两点(设为点A1和点AN)之间形成了一个均匀阵，对该均匀阵进行常规窗函数加权，再对形成的“窗口”进行插值。插值点由各投影点在投影连线(线A1AN)的位置确定，投影点An(其中，n＝1，2，...，N)对应的插值点计算公式如下 其中 表示线段A1An的长度， 表示线段A1AN的长度。根据插值点对常规窗函数在范围内进行插值，插值后的值即为窗函数的修正权值。设所得到的修正窗函数权值为 B＝(6) 其中，bn为对应第n号阵元的修正窗函数权值，n＝1，2，…，N； 步骤五结合步骤三和步骤四的结果，进行频域波束形成得到θ方位的频域波束输出为 具体实例分别采用192元等弧长间距和阵元随机布放的“U”型水听器基阵，结构分别如图3和图4所示。信号采用200kHz单频信号，波束形成的幅度加权方式以窗函数中的海明窗(hamming)加权为例。 假定要形成90度方向的波束，在这个方向上直线ll与“U”型阵的切点在第1号和第192号阵元，因此192个阵元都要参与波束形成。把“U”型阵上各个阵元投影到90度垂直方向(x轴上)，得到点A1，A2，…，A192。假设点A1和点A192为左右两段最边缘的两个点，根据点A1，A2，…，A192在线段A1A192的位置关系利用(4)式生成插值点，再对192元的海明窗进行插值，即得到海明窗的修正权值。 等弧长间距和阵元随机布放的修正海明加权窗函数与海明窗比较分别如图5和图6所示。由各个修正的海本文档来自技高网...
【技术保护点】
一种Ｕ型阵波束形成加权方法，其特征在于：该方法包含以下步骤：步骤一：采用Ｍ元“Ｕ”型阵接收信号，对基阵接收信号进行前置预处理，令阵元编号分别为１，２，…，Ｍ，则阵元输出信号为：ｘ（ｎ）＝［ｘ↓［１］（ｎ），ｘ↓［２］（ｎ），…，ｘ↓［Ｍ］（ｎ）］↑［Ｔ］（１）；步骤二：信号处理机对阵元输出信号进行快速傅立叶变换（ＦＦＴ），将信号从时域变换到频域，频域输出为：Ｘ（ｆ）＝［ｘ↓［１］（ｆ），ｘ↓［２］（ｆ），…，ｘ↓［ｍ］（ｆ）］↑［Ｔ］（２）其中ｘ↓［ｍ］（ｆ）为第ｍ个阵元在频率ｆ的值，ｍ＝１，２，…，Ｍ，设接收信号为带限信号，则ｆ↓［ｌ］≤ｆ≤ｆ↓［ｈ］，ｆ↓［ｌ］，ｆ↓［ｈ］分别为信号频率的下限和上限；步骤三：通过计算斜率为ｔｇθ的直线ｌｌ与“Ｕ”型阵的切点位置Ｐ，形成θ方位的波束，根据Ｐ点位置，确定第１，２，…，Ｎ号阵元参加波束形成，由（２）式选出这Ｎ个阵元的频域输出为：＊（ｆ）＝［ｘ↓［１］（ｆ），ｘ↓［２］（ｆ），…，ｘ↓［Ｎ］（ｆ）］↑［Ｔ］（３）则波束形成的驾驶向量为：Ｗ（ｆ，θ）＝［ｗ↓［１］（ｆ，θ），ｗ↓［２］（ｆ，θ），…，ｗ↓［Ｎ］（ｆ，θ）］（４）其中ｗ↓［ｎ］（ｆ，θ）＝ｅｘｐ（ｊ２πｆτ↓［ｎ］（θ）），τ↓［ｎ］（θ）为第ｎ号阵元到参考点（坐标原点０点）的时延补偿，ｎ＝１，２，…，Ｎ；步骤四：把“Ｕ”型阵中参加波束形成的Ｎ个阵元投影到与θ垂直的波阵面，即直线ｑｑ方向，得到投影点Ａ↓［１］，Ａ↓［２］，…，Ａ↓［Ｎ］；假定在投影点左右两端最边缘的两点，即点Ａ↓［１］和点Ａ↓［Ｎ］之间形成了一个均匀阵，对该均匀阵进行常规窗函数加权，再对形成的“窗口”进行插值；插值点由各投影点在投影连线，即线Ａ↓［１］Ａ↓［Ｎ］的位置确定，投影点Ａ↓［ｎ］（其中，ｎ＝１，２，…，Ｎ）对应的插值点计算公式如下：＊＊＊（５）其中＊表示线段Ａ↓［１］Ａ↓［ｎ］的长度，＊表示线段Ａ↓［１］Ａ↓［Ｎ］的长度；根据插值点对常规窗函数在［０，Ｎ－１］范围内进行插值，插值后的值即为窗函数的修正权值；设所得到的修正窗函数权值为：Ｂ＝［ｂ↓［１］，ｂ↓［２］，…，ｂ↓［Ｎ］］（６）其中，ｂ↓［ｎ］为对应第ｎ号阵元的修正窗函数权值，ｎ＝１，２，…，Ｎ；步骤五：结合步骤三和步骤四的结果，进行频域波束形成得到θ方位的频域波束输出为：Ｙ（ｆ，θ）＝（Ｂ．＊Ｗ（ｆ，θ））′＊＊（ｆ）（７）。...
【技术特征摘要】
【专利技术属性】
技术研发人员：胡鹏，丁烽，李然威，
申请(专利权)人：中国船舶重工集团公司第七一五研究所，
类型：发明
国别省市：86[中国|杭州]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。