阅读说明：本技术 一种基于分块矩阵的快速自适应波束形成算法 (Rapid self-adaptive beam forming algorithm based on block matrix ) 是由 谷新禹 孙微 于 2020-12-14 设计创作，主要内容包括：本发明公开了一种基于分块矩阵的快速自适应波束形成算法,主要包括矩阵分块处理,将高维矩阵分拆为低维矩阵,利用自适应信号处理的特殊性以及递归处理方法,对低维矩阵运算过程中的重叠部分进行整理合并,降低计算复杂度。本发明不需要对数据协方差矩阵进行求逆,同时矩阵维度只有原矩阵的一半,降低了工程实现难度,并且降维运算节省了大量运算资源,对于大规模基阵信号处理效果显著,经验证本发明较常规自适应算法计算效率提高20％以上。(The invention discloses a fast self-adaptive beam forming algorithm based on a block matrix, which mainly comprises matrix block processing, a high-dimensional matrix is divided into low-dimensional matrices, and overlapping parts in the operation process of the low-dimensional matrices are sorted and combined by utilizing the specificity of self-adaptive signal processing and a recursion processing method, so that the calculation complexity is reduced. According to the invention, the data covariance matrix does not need to be inverted, meanwhile, the matrix dimension is only half of that of the original matrix, the engineering realization difficulty is reduced, the dimension reduction operation saves a large amount of operation resources, the large-scale matrix signal processing effect is obvious, and the calculation efficiency is improved by more than 20% compared with that of the conventional adaptive algorithm through verification.)

1. A fast adaptive beamforming algorithm based on a blocking matrix is characterized in that: according to the method, a high-dimensional matrix is divided into low-dimensional matrices through matrix blocking processing, and overlapping parts in the operation process of the low-dimensional matrices are sorted and combined by using the specificity of adaptive signal processing and a recursive processing method; the method comprises the following specific steps:

1) based on the idea of a divide-and-conquer method, partitioning high-dimensional array data: assuming that the number of array elements is N, and the input data per snapshot is X ═ X₁ x₂ … x_N]^H，x_iFor the sampled data of the ith array element, let the array be uniform, with the array element spacing d, and the steering vector w ═ exp (-j ω τ), where ω is the angular frequency, τ is_i＝d_isin (theta)/c is the delay of the ith array element relative to the first array element in the theta direction, and the spatial energy spectrum formed by the MVDR adaptive beam isAnd R is a data covariance matrix, partitioning R, equally dividing into four parts, and considering the covariance matrix as a Hermite matrix to obtain:

wherein, A, B and C are N/2 dimensional square matrixes, and R is a Hermite matrix, so that the transposes of the upper right corner matrix and the lower left corner matrix are the same, namely, the A, B and C can represent R, and the w is made to be [ w ] by combining the block matrix inversion principle and the amplification matrix inversion formula₁,w₂]And then:

wherein S is B-C^HA^-1C，

2) Using the cholesky decomposition, a matrix of variables with fixed patterns is decomposed into the product of two matrices, and a simple expression for 1/P is obtained by the joint cancellation between the variables: let w₂＝αw₁Where α ═ exp (-j ω τ)_M/2) Is w₁And w₂1/P is simplified to:

wherein S^-1＝s*s^HAnd s is obtained by matrix decomposition, then:

whereinTo matrix I-YY^HPerforming cholesky decomposition to obtain:

I-YY^H＝q*q^H

wherein

Therefore, it is not only easy to uses_n+1＝s_nq, same as A^-1Decomposition is also carried out:

whereinThenAnd calculating to obtain the spatial energy spectrum of the MVDR algorithm.

Technical Field

The invention relates to the field of sonar array signal processing, in particular to a fast self-adaptive beam forming algorithm based on a block matrix.

Background

In the sonar detection process, a beam forming technology is needed to be used for processing array receiving signals, enhancing target signals and suppressing interference and noise signals. With the development of science and technology and the increasing complexity of underwater environment, the requirements on the beam forming quality are higher and higher. Among various beamforming algorithms, the adaptive algorithm is drawing attention due to its high resolution and high output signal-to-noise ratio, and is currently widely applied to various types of ocean exploration equipment. However, the core of the adaptive algorithm is to perform update iteration on the steering vector by using data feedback, and an inversion operation needs to be performed on a data covariance matrix. Because matrix multiplication and inversion operation in engineering application need to consume a large amount of resources, the calculation complexity is exponentially increased along with the increase of the matrix scale, the array scale of the underwater acoustic detection equipment is larger and larger at present, for example, underwater imaging sonar equipment can reach the scale of hundreds of array elements, and the complicated calculation of a high-dimensional matrix seriously restricts the performance of the underwater acoustic detection equipment.

Disclosure of Invention

The invention aims to overcome the defects in the prior art and provide a fast self-adaptive beam forming algorithm based on a block matrix.

The purpose of the invention is achieved by the following technical scheme: the fast self-adaptive beam forming algorithm based on the block matrix divides a high-dimensional matrix into low-dimensional matrices through matrix block processing, and arranges and combines overlapped parts in the operation process of the low-dimensional matrices by utilizing the particularity of self-adaptive signal processing and a recursion processing method; the method mainly comprises the following steps:

1) based on the idea of a divide-and-conquer method, high-dimensional array data is partitioned, and matrix inversion operation is avoided by using a recursion mode: assuming that the number of array elements is N, and the input data per snapshot is X ═ X₁ x₂ … x_N]^H，x_iFor the sampled data of the ith array element, let the array be uniform, with the array element spacing d, and the steering vector w ═ exp (-j ω τ), where ω is the angular frequency, τ is_i＝d_isin (theta)/c is the delay of the ith array element relative to the first array element in the theta direction, and the spatial energy spectrum formed by the MVDR adaptive beam isR is a data covariance matrix, when the magnitude of N reaches hundreds or even thousands, no matter the inversion or matrix multiplication of R, a large amount of resources are consumed, the practicability is reduced, in order to solve the problem, the matrix can be partitioned from the aspect of a divide and treat method, the data covariance matrix R is partitioned, the data covariance matrix R is divided into four parts, and the covariance matrix is considered to be Hermite momentAnd (3) obtaining:

wherein A, B and C are N/2 dimensional square matrix, and w is [ w ] by combining block matrix inversion principle and amplification matrix inversion formula₁,w₂]And then:

wherein S is B-C^HA^-1C，

2) By using cholesky decomposition, a variable matrix with a fixed mode is decomposed into a product of two matrixes, and a simple expression about 1/P is obtained through the combination and cancellation among the variables, so that the calculation efficiency is improved: let w₂＝αw₁Where α ═ exp (-j ω τ)_M/2) Is w₁And w₂1/P is simplified to:

wherein S^-1＝s*s^HS is obtained by matrix decomposition, then

WhereinFor the matrix I-YY^HPerforming cholesky decomposition to obtain:

I-YY^H＝q*q^H

whereinTherefore, it is not only easy to uses_n+1＝s_nq is calculated. Same principle A^-1Decomposition is also carried out:

whereinThen 1/P ═ w₁a)(w₁a)^H+(w₁(αI-A^-1C)s)(w₁(αI-A^-1C)s)^H＝P₁P₁ ^H+P₂P₂ ^HAnd calculating to obtain the spatial energy spectrum of the MVDR algorithm. .

The invention has the beneficial effects that: the invention provides a divide and conquer method thought, a high-dimensional matrix is divided into low-dimensional matrices through matrix block processing, and overlapped parts in the operation process of the low-dimensional matrices are sorted and combined by utilizing the particularity of adaptive signal processing and a recursion processing method, so that the calculation complexity is reduced; the data covariance matrix does not need to be inverted, meanwhile, the matrix dimension is only half of that of the original matrix, the engineering realization difficulty is reduced, a large amount of operation resources are saved through dimension reduction operation, the large-scale matrix signal processing effect is obvious, and the calculation efficiency is improved by more than 20% through the verification compared with that of a conventional adaptive algorithm.

Drawings

FIG. 1 is a schematic structural diagram of the present invention.

FIG. 2 is a time-consuming comparison between the present invention and a conventional method.

FIG. 3 is a schematic diagram of the improvement of the present invention over the conventional method.

Detailed Description

The invention will be described in detail below with reference to the following drawings:

as shown in fig. 1-3, the fast adaptive beamforming algorithm based on the block matrix mainly includes the following steps:

1) let the array element number be N, and input data per snapshot be X ═ X₁ x₂ … x_N]^H,x_iThe sampled data of the ith array element is the data covariance matrix R

Wherein a, b and s are weighting factors representing the proportion of new and old components in the covariance matrix, and are all set to 1 for convenient operation. Assuming the array is uniform and the array elements are spaced apart by d, the steering vector is w ═ exp (-j ω τ), where ω is the angular frequency and τ is_i＝d_isin (theta)/c is the delay of the ith array element relative to the first array element in the theta direction, and taking the MVDR adaptive beam forming example, the spatial energy spectrum is:

where R is the data covariance matrix. The scale of the sonar basic array influences the calculation efficiency, and when the number of array elements reaches hundreds or even thousands, no matter the inversion or the matrix multiplication of R, a large amount of resources are consumed, and the practicability is reduced. In order to solve the problem, from the aspect of divide and conquer method, the matrix can be processed in a blocking way, so that the dimension of the matrix can be reduced, and the calculation efficiency is improved.

2) Partitioning the R into four parts, and considering the covariance matrix as a Hermite matrix to obtain the final product

Wherein A, B and C are N/2 dimensional square matrixes, and can be obtained according to the inversion principle of the block matrix

Wherein S is reversible and S ═ B-C^HA^-1C, as can be seen from the formula (4), since A, B and C are obviously known, only two N/2-dimensional matrixes A need to be calculated^-1,S^-1To obtain R^-1。

For matrix A, there are

A_n+1＝A_n+X₁X₁ ^H (5)

Wherein X₁＝[x₁ x₂ … x_N/2]^H

Recombination of augmented matrix inversion formula

(M+NPQ)^-1＝M^-1-M^-1N(QM^-1N+P^-1)^-1QM^-1 (6)

Let M be A_n,N＝X₁,When P is 1, the inverse matrix of a can be obtained by recursion of the following equation

WhereinFor S, the same principle

WhereinObviously E is a one-dimensional column vector, so we will bring (8) into equation (6):

wherein

3) An inverse matrix of a and S is obtained, and then the spatial energy spectrum P can be calculated by combining equation (2) with w ═ w₁,w₂]，w₁,w₂Respectively, the guide vectors of the block matrixes, then:

1/P＝[w₁,w₂]R^-1[w₁,w₂]^H

can be brought into formula (4)

In equation (10), a large number of matrix multiplication operations exist, and simplification of processing is required. In the case of a uniform array far field, the steering vector is linear, w₂＝αw₁Where α ═ exp (-j ω τ)_M/2) Is w₁,w₂Is observed, and the presence of A in (10)^{- 1}C,S^-1The factor term, then the formula (10) can be simplified to

Wherein S^-1＝s*s^HS can be obtained by matrix decomposition, but is computationally complex and requires the computation of S^-1The calculation amount is large, and a recursion form is considered here. According to equation (9), assume S^-1Can be written as s^HIn the form of

WhereinFor the matrix I-YY^HCholesky decomposition proceeds due to I-YY^HIs given by the particularity of

I-YY^H＝q*q^H (13)

Wherein

Therefore, it is not only easy to uses_n+1＝s_nq。

Same principle A^-1Decomposition can also be performed, avoiding complex calculations of equation (7):

wherein

Then, the formula (11) can be simplified to

Wherein P is₁＝w₁a, can be obtained by recursive calculation using the formulae (15), (16) and (17), P₂＝w₁(αI-A^-1C) s, which can be obtained by recursion using the equations (12), (13) and (14), the spatial energy spectrum of the MVDR algorithm can be obtained by calculation.

Fig. 1 is a schematic block diagram of the present invention, which takes MVDR beam forming method as an example. By using the augmented matrix inversion formula and the block matrix inversion principle, the inverse matrix A of the covariance matrix subarray A can be obtained^-1Recursion decomposition of (A) and inverse matrix S of Schur' S complement S of B^-1By using factorization, R can be simplified^-1Each element of (1) to R^-1Can be written as P₁P₁ ^H+P₂P₂ ^HThe characteristic of linear relation of the uniform array guide vectors is utilized, and a dimension reduction expression of the self-adaptive beam forming space spectrum can be obtained.

The MVDR spatial energy spectrum is calculated for arrays of different sizes, and the average value is calculated 200 times, and the result is shown in figures 2 and 3. FIG. 2 shows the time consumption of the method compared with the conventional method, and FIG. 3 shows the improvement of the method compared with the conventional method. It can be seen from the figure that the calculation efficiency of the method of the invention is obviously better than that of the traditional method along with the increase of the array element scale, and when the array element scale is larger than 256 elements, the improvement degree is larger than 20%. It is seen from the figure that, when the array element size is smaller than 256 elements, the calculation efficiency of the method of the present invention is inferior to that of the conventional method, because the debugging machine performs simulation operation by using Matlab, the matrix inversion function carried by Matlab is an optimized result, the calculation efficiency is extremely excellent in the case of low dimension, and in practical engineering application, software and hardware means such as C language and FPGA are mostly used for calculation without the support of Matlab function library, so the method of the present invention has superior effect in such a case.

It should be understood that equivalent substitutions and changes to the technical solution and the inventive concept of the present invention should be made by those skilled in the art to the protection scope of the appended claims.