CN113422629A - Covariance matrix reconstruction self-adaptive beam forming method and system - Google Patents
Covariance matrix reconstruction self-adaptive beam forming method and system Download PDFInfo
- Publication number
- CN113422629A CN113422629A CN202110673467.2A CN202110673467A CN113422629A CN 113422629 A CN113422629 A CN 113422629A CN 202110673467 A CN202110673467 A CN 202110673467A CN 113422629 A CN113422629 A CN 113422629A
- Authority
- CN
- China
- Prior art keywords
- covariance matrix
- data
- array
- virtual
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 138
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000013598 vector Substances 0.000 claims abstract description 68
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 17
- 238000000605 extraction Methods 0.000 claims abstract description 10
- 230000003044 adaptive effect Effects 0.000 claims description 24
- 238000001228 spectrum Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000011068 loading method Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 4
- 230000001629 suppression Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B15/00—Suppression or limitation of noise or interference
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a covariance matrix reconstruction self-adaptive beam forming method and a covariance matrix reconstruction self-adaptive beam forming system based on a virtual linear array.A uniform linear array is used for receiving time domain data, forward and backward linear prediction coefficients are respectively calculated, and the received data at the position of a virtual array element are fitted in sequence according to a recursion mode, so that snapshot data of each virtual array element are predicted; then, calculating an interference and noise covariance matrix with higher dimensionality by using data containing actual array elements and virtual array elements, removing reconstruction errors by adopting a mode of characteristic decomposition and large characteristic value extraction, and finally adding an ideal noise covariance matrix to obtain an improved interference and noise covariance matrix corresponding to the virtual array; and finally, reconstructing an expected signal covariance matrix corresponding to the virtual array in the same way, and extracting the eigenvector corresponding to the maximum eigenvalue as a guide vector, thereby calculating to obtain a beam former weighting vector and beam output corresponding to the virtual array.
Description
Technical Field
The invention relates to the field of array signal processing, in particular to a covariance matrix reconstruction self-adaptive beam forming method and a covariance matrix reconstruction self-adaptive beam forming system.
Background
Adaptive beamforming techniques for sensor arrays have evolved rapidly over the past several decades, and related efforts have been widely applied in various fields. As the best known adaptive beamformer, Capon beamformer can achieve the best of its theoretical performance when the training data does not contain the desired signal. However, in practical application environments, a desired signal often exists in received data of an array, and particularly when the signal-to-noise ratio is high, the classical adaptive beamforming method is more sensitive to mismatching of a steering vector of the desired signal.
In order to improve the robustness of the adaptive beam former, in recent years, scholars at home and abroad propose a plurality of robust beam formers, which are summarized into the following two types according to the used covariance matrix:
1) the diagonal loading method and the indeterminate set constraint method disclosed in the document 1 "Robust Adaptive Beamforming," IEEE Transactions on optics Speech and Signal Processing, vol.35, No.10, pp.1365-1376, Oct,1987 "and the document 2" Robust Caption Beamforming, "IEEE Signal Processing Letters, vol.10, No.6, pp.172-175, Jun, 2003" can improve the robustness well, but the diagonal loading amount and the indeterminate set size are difficult to determine according to actual conditions and difficult to use practically.
2) The reconstructed Interference-plus-noise Covariance Matrix (INCM) algorithm disclosed in document 3 "Robust Adaptive Beamforming Based on Interference Covariance Matrix Reconstruction and engineering Vector Estimation, IEEE trans signal process, vol.60, No.7, pp.3881-3885, jul.2012" can effectively remove the desired signal component in the Covariance Matrix, thereby reducing the sensitivity of the beamformer to the desired signal Steering Vector error, but at the same time, there may be interfering Steering Vector errors in the reconstructed INCM, thereby reducing the Interference suppression capability.
To sum up, the method of the type 1) directly utilizes a received data covariance matrix (SCM), and greatly sacrifices the interference suppression capability in order to ensure the robustness; the method of the type 2) can effectively ensure the robustness, but when the model mismatch exists, the interference suppression capability is reduced. Both algorithms are based on array element received data, and the difference value of the output signal-to-interference-and-noise ratio and the input signal-to-noise ratio does not exceed the array gain.
Disclosure of Invention
In order to avoid the defects of the existing adaptive beam forming, the invention provides a covariance matrix reconstruction adaptive beam forming method based on a virtual linear array.
In order to achieve the purpose, the invention adopts the technical scheme that: a covariance matrix reconstruction self-adaptive beam forming method comprises the following specific processes:
acquiring time domain data received by a uniform linear array, respectively calculating a forward fitting coefficient and a backward fitting coefficient, and sequentially fitting the data received by the virtual array elements in a recursion mode, thereby predicting snapshot data of each virtual array element;
calculating an interference-plus-noise covariance matrix with higher dimensionality by utilizing snapshot data containing actual array elements and virtual array elements, eliminating reconstruction errors by adopting a mode of characteristic decomposition and extraction of larger characteristic values, and finally introducing an ideal noise covariance matrix to obtain an improved interference-plus-noise covariance matrix corresponding to the virtual array;
and reconstructing an expected signal covariance matrix corresponding to the virtual array in the same way of calculating and improving the interference and noise covariance matrix, extracting a characteristic vector corresponding to the maximum characteristic value of the expected signal covariance matrix to be used as a guide vector of an expected signal, calculating a beam former weighting vector corresponding to the virtual array according to the guide vector, and further obtaining beam output according to the weighting vector and a receiving data matrix of the virtual linear array.
Acquiring time domain data received by a uniform linear array, and respectively calculating forward and backward linear prediction coefficients, wherein the specific process comprises the following steps:
the uniform linear array receives time domain snapshot data of signals and noise from space;
calculating forward sample covariance matrixes corresponding to other array element data except the 1 st vibration element in the uniform linear array, and calculating a fitting vector of forward fitting;
calculating a forward fitting coefficient according to the fitting vector of the forward fitting;
calculating forward sample covariance matrixes corresponding to the data of the residual vibration elements except the last vibration element in the uniform linear array, and calculating a fitting vector of forward fitting;
and calculating to obtain a backward fitting coefficient according to the fitting vector.
Sequentially fitting the data received at the positions of the virtual array elements in a recursion mode, so that the process of predicting the snapshot data of each virtual array element is as follows:
taking the data of the vibration elements except the last vibration element in the virtual vibration array, calculating the forward fitting data of the vibration element No. 0, and sequentially calculating to the nth forward fitting data;
and taking the data of the other vibration elements except the first vibration element in the virtual vibration array, calculating the data of the maximum vibration element with the backward number, and sequentially and respectively calculating the nth backward fitting data.
Calculating an interference and noise covariance matrix with higher dimensionality by utilizing snapshot data containing actual array elements and virtual array elements, eliminating reconstruction errors by adopting a mode of characteristic decomposition and extraction of larger characteristic values, and finally adding an ideal noise covariance matrix to obtain an improved interference and noise covariance matrix corresponding to a virtual array and an expected signal covariance matrix corresponding to a reconstructed virtual array, wherein the specific process is as follows:
combining the forward fitting data, the original array element receiving data and the backward fitting data together to obtain a receiving data matrix of the virtual linear array;
calculating a sample covariance matrix of the virtual linear array according to the received data matrix;
according to the array element number of the virtual linear array, the guide vector of the corresponding angle is obtained,
calculating diagonal loading Capon space spectrums corresponding to all angles based on the guide vectors;
dividing an expected signal sector according to the Capon space spectrum, and then dividing interference and noise sectors except the expected sector;
preliminarily calculating an interference-plus-noise covariance matrix in the interference-plus-noise sector;
performing characteristic decomposition on the interference and noise covariance matrix;
obtaining a corrected interference and noise covariance matrix according to the characteristic decomposition result;
the desired signal covariance matrix is reconstructed within the desired signal sector.
The method comprises the following steps of obtaining uniform linear array received time domain data, calculating forward and backward linear prediction coefficients respectively, fitting the data received by virtual array elements in sequence according to a recursion mode, and predicting snapshot data of each virtual array element, wherein the method specifically comprises the following steps:
a uniform linear array of M transducers receives time domain snapshot data x (k) x ═ x of signals and noise from space1(k),x2(k),…,xM(k)]TWhere the signal is uncorrelated with noise, xm(k) Representing the time domain data received by the mth array element, and calculating the data x of the 2 nd to M array elementsF(k)=[x2(k),…,xM(k)]TCorresponding forward sample covariance matrix
Calculating a fitting vector d of the forward fitF:
Finally obtaining a forward fitting coefficient alphaF:
Calculating the 1 st to M-1 st array metadata XB(k)=[x1(k),…,xM-1(k)]TThe corresponding forward sample covariance matrix,
recalculating the fitting vector d of the forward fitB:
Finally obtaining a backward fitting coefficient beta:
taking the array element data x from 1 st to M-1 sttemp(k)=[x1(k),…,xM-1(k)]TCalculating the forward fitting data of the array element No. 0,
according to the rule, the nth forward fitting data are respectively calculated in sequence:
taking the 2 nd to M array element data xtemp(k)=[x2(k),…,xM(k)]TCalculating the data of the backward M +1 number array element,
according to the rule, the nth backward fitting data are sequentially and respectively calculated:
combining the forward fitting data, the original array element receiving data and the backward fitting data together to obtain a receiving data matrix of the virtual linear array,
xLP(k)=[x1-N(k),…,x1(k),…,xM(k),…,xM+N(k)]T,
then, a sampling covariance matrix is calculated,
calculating an interference-plus-noise covariance matrix with higher dimensionality by utilizing snapshot data containing actual array elements and virtual array elements, eliminating reconstruction errors by adopting a mode of characteristic decomposition and large characteristic value extraction, and finally combining an ideal noise covariance matrix to obtain an improved interference-plus-noise covariance matrix corresponding to the virtual array; the method comprises the following specific steps:
according to the array element number of the virtual linear array, the corresponding guide vector a of the theta angle is obtainedLP(theta), calculating diagonal loading Capon space spectrums p (theta) corresponding to all the angles,
where I is the unit diagonal matrix,
according to Capon space spectrum, dividing expected signal sectorWhere B is the set region width, and then demarcating interference plus noise sectors other than the desired sectorFor the complement of Θ, the interference-plus-noise co-ration is initially calculated in the interference-plus-noise sectorVariance matrix
Wherein: lambda [ alpha ]m(M is 1, …, M) is a feature value arranged in descending order, v ismIs the corresponding feature vector; dintRepresenting a diagonal matrix containing all the larger eigenvalues, VintIncluding eigenvectors corresponding to larger eigenvalues, ReRepresenting an error covariance matrix consisting of smaller eigenvalues and eigenvectors, the larger eigenvalues being based onEstimating the number of interference signals as Q, adding 1 signal number as Q +1, wherein Q +1 is the minimum positive integer for ensuring the above formula to be satisfied; the number of Q interference signals arranged from large to small plus 1 signal number Q +1 forms a larger characteristic value; at this time, the interference plus noise covariance matrix is modified as:
Wherein theta islIs the l-th discrete angle within Θ, estimates the steering vector of the desired signal as,
computing beamformer weighting vectors w corresponding to virtual arraysLP:
Calculating beam output y (k) from the weight vector and the received data matrix of the virtual linear array:
a covariance matrix reconstruction adaptive beam forming system comprises a snapshot data acquisition module, a covariance matrix calculation module,
The snapshot data acquisition module is used for acquiring time domain data received by the uniform linear array, respectively calculating forward and backward linear prediction coefficients, and sequentially fitting the data received by the virtual array elements in a recursion mode, so as to predict snapshot data of each virtual array element;
the covariance matrix calculation module is used for calculating an interference-plus-noise covariance matrix with higher dimensionality according to snapshot data containing actual array elements and virtual array elements, eliminating reconstruction errors in a mode of characteristic decomposition and large characteristic value extraction, introducing an ideal noise covariance matrix, and obtaining an improved interference-plus-noise covariance matrix corresponding to the virtual array;
the output beam calculation module is used for reconstructing an expected signal covariance matrix corresponding to the virtual array in the same way as the interference and noise covariance matrix is calculated and improved, extracting a characteristic vector corresponding to the maximum characteristic value as a guide vector of an expected signal, calculating a beam former weighting vector corresponding to the virtual array according to the guide vector, and obtaining beam output according to the weighting vector and a receiving data matrix of the virtual linear array.
A computer device comprises a processor and a memory, wherein the memory is used for storing a computer executable program, the processor reads part or all of the computer executable program from the memory and executes the computer executable program, and the processor can realize the covariance matrix reconstruction adaptive beamforming method when executing part or all of the computer executable program.
A computer readable storage medium having stored thereon a computer program which, when executed by a processor, is capable of implementing the covariance matrix reconstruction adaptive beamforming method of the present invention.
Compared with the prior art, the invention has at least the following beneficial effects: according to the invention, linear prediction and covariance matrix reconstruction technology are combined, and received data of forward and backward virtual array elements are fitted by linear prediction, so that the number of channels of the array is increased. After the linear prediction and the covariance matrix reconstruction technology are combined, not only can the array gain of the adaptive beam former be improved, but also the complexity of the adaptive beam former can be improved, so that the beam former can form deeper depression in the interference direction when the array error is small, and can suppress interference through lower sidelobes when the array error is large, and finally the interference and noise suppression capability of the adaptive beam former is improved.
Drawings
Fig. 1 shows the output SINR of the proposed virtual array beamforming method as a function of the input SNR.
Fig. 2 shows the variation of the output SINR of the proposed virtual array beamforming method with the input fast beat number.
Fig. 3 shows the variation of the output SINR of the proposed virtual array beamforming method with the forward-backward fitting number N.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The invention utilizes the existing uniform linear array to receive data, and fits the virtual array elements in a linear prediction mode, thereby forming the virtual array, improving the number of channels of the array and improving the performance of self-adaptive beam forming.
Examples
Simulation 1: the underwater sound velocity is 1500m/s, a 10-wire array is considered, the assumed array element spacing is 0.5m, the expected signal is assumed to be located at 15 degrees, the signal-to-noise ratio is changed from-10 dB to 20dB, the interference is assumed to be located at-25 degrees and 35 degrees, the dry-to-noise ratio is 20dB, and the noise is Gaussian uniform white noise. The desired signal and the interfering signal are both narrow-band incoherent signals with a center frequency of 1500 Hz. N is 5 virtual array elements are fitted in the front and back directions, and the fast beat number is K is 30.
In order to illustrate the statistical performance of the method, the simulation is performed for 200 times, and finally the statistical result of the output new interference-to-noise ratio changing along with the input signal-to-noise ratio is obtained. For each simulation, the following operation steps of the method are repeated to obtain a weight vector:
step 1: a uniform linear array of 10 hydrophones receives time domain snapshot data of signal and noise from space x (k) ═ x1(k),x2(k),…,x10(k)]TWhere the signal is uncorrelated with noise, xm(k) Representing the time domain data received by the mth array element. Virtual vibration element and actual vibration elementForming a virtual linear array together, and calculating the 2 nd to 10 th array element data xF(k)=[x2(k),…,x10(k)]TThe corresponding forward sample covariance matrix,
then calculating the fitting vector of the forward fitting,
Next, the 1 st to 9 th array metadata X are calculatedB(k)=[x1(k),…,x9(k)]TThe corresponding forward sample covariance matrix,
then calculating the fitting vector of the forward fitting,
Step 2: taking the 1 st to 9 th array element data xtemp(k)=[x1(k),…,x9(k)]TCalculating the forward fitting data of the array element No. 0,
according to the rule, 5 th forward fitting data are sequentially and respectively calculated:
taking 2 nd to 10 th array element data xtemp(k)=[x2(k),…,x10(k)]TAnd calculates the data of the backward 11 number array element,
according to the rule, 5 th backward fitting data are sequentially and respectively calculated:
and step 3: 5 data are fitted in the forward direction and the backward direction, the forward fitting data, the original array element receiving data and the backward fitting data are combined together to obtain a receiving data matrix of the virtual linear array,
xLP(k)=[x-4(k),…,x1(k),…,x10(k),…,x15(k)]T,
then, a sampling covariance matrix is calculated,
setting the discrete angle interval to 1 deg., then-90 deg., 90 deg]A total of 181 discrete angles. According to the array element number of the virtual linear array, the corresponding guide vector a of the theta angle is obtainedLP(theta), which contains 20 elements, calculates the diagonal loaded Capon space spectrum corresponding to each angle,
where I is the unity diagonal matrix.
And 4, step 4: estimating the angle of the expected signal according to the Capon space spectrumIf the set region width B is 10, the desired signal sector Θ is divided into [10 °,20 °]Then define interference plus noise sectors other than the desired sectorThere are a total of 170 discrete angles, 181-11. Preliminary computation of interference-plus-noise covariance matrix in interference-plus-noise sectors
Wherein theta islIs thatThe ith discrete angle. Then, the characteristic decomposition is carried out,
wherein: lambda [ alpha ]m(m-1, …,10) is a feature value arranged in descending order, v ismIs the corresponding feature vector; dintRepresenting a diagonal matrix containing all the larger eigenvalues, VintIncluding eigenvectors corresponding to larger eigenvalues, ReRepresenting an error covariance matrix consisting of smaller eigenvalues and eigenvectors, the larger eigenvalues being based onEstimating the number of interference signals as Q, adding 1 signal number as Q +1, wherein Q +1 is the minimum positive integer for ensuring the above formula to be satisfied; the greater characteristic is formed by adding 1 signal number Q +1 to the number of Q interference signals arranged from large to smallA value; at this point, the interference plus noise covariance matrix may be modified to,
and 5: there are 11 discrete angles in the desired signal sector, and the desired signal covariance matrix is reconstructed as described in step 4, resulting in,
wherein theta islIs the l-th discrete angle within Θ, estimates the steering vector of the desired signal as,
step 6: the corresponding beamformer weight vector for the virtual array is calculated,
calculating a beam output based on the weight vector and the received data matrix of the virtual linear array,
fig. 1 shows a comparison of the proposed LP-REB with REB and RCB beamforming methods when the fast beat number is 30, it can be seen that RCB performs poorly at high signal-to-noise ratio, whereas REB suffers from below the optimal upper limit (optimal) due to the presence of reconstruction errors. And due to the use of the virtual array, the LP-REB increases the total array element number (10 real array elements and 10 virtual array elements), can effectively improve the output signal-to-interference-and-noise ratio, and is even better than the optimal upper limit of a 10-element array.
Simulation 2: the input signal-to-noise ratio is kept at 10dB, the snapshot number K is changed from 10 to 100 at an interval of 10. the rest of simulation conditions are the same as that of simulation 1, 200 Monte Carlo experiments are carried out, and a graph of the change of the average value of the output signal-to-interference-and-noise ratio along with the snapshot number is obtained. The steps are the same as the corresponding steps of simulation 1.
Fig. 2 shows the variation of the output snr of each beamformer with the snapshot number, and it can be seen that the RCB algorithm has gradually enhanced performance due to the estimation error of the sample covariance matrix at low sample number, but still much lower than 20dB, and REB and LP-REB are insensitive to the snapshot number due to the reconstructed covariance matrix. The output signal-to-noise ratio of the LP-REB is far higher than that of the REB due to the adoption of more virtual array metadata.
Simulation 3: the input signal-to-noise ratio is fixed to be 10dB, the snapshot number is fixed to be 30, the forward and backward fitting number N is changed from 0 to 9 at intervals of 1, the rest simulation conditions are the same as those of simulation 1, 200 Monte Carlo experiments are carried out, and a change graph of the average value of the output signal-to-interference-and-noise ratio along with the forward and backward fitting number N is obtained.
Fig. 3 shows the change of the SINR of the LP-REB with the number of forward and backward fits, and it can be seen that as the number of forward and backward fits increases, the number of virtual array elements also increases, thereby effectively improving the gain of the virtual array and finally improving the SINR.
The invention provides a covariance matrix reconstruction self-adaptive beam forming system, which comprises a snapshot data acquisition module, a covariance matrix calculation module,
The snapshot data acquisition module is used for acquiring time domain data received by the uniform linear array, respectively calculating forward and backward linear prediction coefficients, and sequentially fitting the data received by the virtual array elements in a recursion mode, so as to predict snapshot data of each virtual array element;
the covariance matrix calculation module is used for calculating an interference-plus-noise covariance matrix with higher dimensionality according to snapshot data containing actual array elements and virtual array elements, eliminating reconstruction errors in a mode of characteristic decomposition and large characteristic value extraction, introducing an ideal noise covariance matrix, and obtaining an improved interference-plus-noise covariance matrix corresponding to the virtual array;
the output beam calculation module is used for reconstructing an expected signal covariance matrix corresponding to the virtual array in the same way as the interference and noise covariance matrix is calculated and improved, extracting a characteristic vector corresponding to the maximum characteristic value as a guide vector of an expected signal, calculating a beam former weighting vector corresponding to the virtual array according to the guide vector, and obtaining beam output according to the weighting vector and a receiving data matrix of the virtual linear array.
The invention can also provide a computer device, which comprises a processor and a memory, wherein the memory is used for storing a computer executable program, the processor reads part or all of the computer executable program from the memory and executes the computer executable program, and when the processor executes part or all of the computer executable program, the covariance matrix reconstruction adaptive beamforming method can be realized.
In another aspect, the present invention provides a computer readable storage medium, in which a computer program is stored, and the computer program, when being executed by a processor, can implement the covariance matrix reconstruction adaptive beamforming method according to the present invention.
The computer equipment can be an onboard computer, a notebook computer, a tablet computer, a desktop computer, a mobile phone or a workstation.
The processor may be a Central Processing Unit (CPU), a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), or an off-the-shelf programmable gate array (FPGA).
The memory of the invention can be an internal storage unit of a vehicle-mounted computer, a notebook computer, a tablet computer, a desktop computer, a mobile phone or a workstation, such as a memory and a hard disk; external memory units such as removable hard disks, flash memory cards may also be used.
Computer-readable storage media may include computer storage media and communication media. Computer storage media includes volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data. The computer-readable storage medium may include: a Read Only Memory (ROM), a Random Access Memory (RAM), a Solid State Drive (SSD), or an optical disc. The Random Access Memory may include a resistive Random Access Memory (ReRAM) and a Dynamic Random Access Memory (DRAM).
Claims (9)
1. A covariance matrix reconstruction adaptive beam forming method is characterized by comprising the following specific processes:
acquiring time domain data received by a uniform linear array, respectively calculating a forward fitting coefficient and a backward fitting coefficient, and sequentially fitting the data received by the virtual array elements in a recursion mode, thereby predicting snapshot data of each virtual array element;
calculating an interference-plus-noise covariance matrix with higher dimensionality by utilizing snapshot data containing actual array elements and virtual array elements, eliminating reconstruction errors by adopting a mode of characteristic decomposition and extraction of larger characteristic values, and finally introducing an ideal noise covariance matrix to obtain an improved interference-plus-noise covariance matrix corresponding to the virtual array;
and reconstructing an expected signal covariance matrix corresponding to the virtual array in the same way of calculating and improving the interference and noise covariance matrix, extracting a characteristic vector corresponding to the maximum characteristic value of the expected signal covariance matrix to be used as a guide vector of an expected signal, calculating a beam former weighting vector corresponding to the virtual array according to the guide vector, and further obtaining beam output according to the weighting vector and a receiving data matrix of the virtual linear array.
2. The covariance matrix reconstruction adaptive beamforming method of claim 1, wherein the uniform linear array received time domain data is obtained, and forward and backward linear prediction coefficients are calculated respectively, the specific process is as follows:
the uniform linear array receives time domain snapshot data of signals and noise from space;
calculating forward sample covariance matrixes corresponding to other array element data except the 1 st vibration element in the uniform linear array, and calculating a fitting vector of forward fitting;
calculating a forward fitting coefficient according to the fitting vector of the forward fitting;
calculating forward sample covariance matrixes corresponding to the data of the residual vibration elements except the last vibration element in the uniform linear array, and calculating a fitting vector of forward fitting;
and calculating to obtain a backward fitting coefficient according to the fitting vector.
3. The covariance matrix reconstruction adaptive beamforming method of claim 1, wherein the data received at the position of the virtual array element is sequentially fitted in a recursive manner, so as to predict the snapshot data of each virtual array element as follows:
taking the data of the vibration elements except the last vibration element in the virtual vibration array, calculating the forward fitting data of the vibration element No. 0, and sequentially calculating to the nth forward fitting data;
and taking the data of the other vibration elements except the first vibration element in the virtual vibration array, calculating the data of the maximum vibration element with the backward number, and sequentially and respectively calculating the nth backward fitting data.
4. The covariance matrix reconstruction adaptive beam forming method according to claim 1, wherein a higher-dimensional interference-plus-noise covariance matrix is calculated by using snapshot data including actual array elements and virtual array elements, reconstruction errors are eliminated by means of feature decomposition and extraction of larger eigenvalues, and an ideal noise covariance matrix is finally added to obtain an improved interference-plus-noise covariance matrix corresponding to a virtual array and an expected signal covariance matrix corresponding to a reconstructed virtual array, and the specific process is as follows:
combining the forward fitting data, the original array element receiving data and the backward fitting data together to obtain a receiving data matrix of the virtual linear array;
calculating a sample covariance matrix of the virtual linear array according to the received data matrix;
according to the array element number of the virtual linear array, the guide vector of the corresponding angle is obtained,
calculating diagonal loading Capon space spectrums corresponding to all angles based on the guide vectors;
dividing an expected signal sector according to the Capon space spectrum, and then dividing interference and noise sectors except the expected sector;
preliminarily calculating an interference-plus-noise covariance matrix in the interference-plus-noise sector;
performing characteristic decomposition on the interference and noise covariance matrix;
obtaining a corrected interference and noise covariance matrix according to the characteristic decomposition result;
the desired signal covariance matrix is reconstructed within the desired signal sector.
5. The covariance matrix reconstruction adaptive beamforming method of claim 1, wherein the method comprises the steps of obtaining time domain data received by a uniform linear array, calculating forward and backward linear prediction coefficients, fitting the data received by the virtual array elements in turn according to a recursive method, and predicting snapshot data of each virtual array element, and specifically comprises the steps of:
a uniform linear array of M transducers receives time domain snapshot data x (k) x ═ x of signals and noise from space1(k),x2(k),…,xM(k)]TWherein the signal and the noiseAcoustic uncorrelated, xm(k) Representing the time domain data received by the mth array element, and calculating the data x of the 2 nd to M array elementsF(k)=[x2(k),…,xM(k)]TCorresponding forward sample covariance matrix
Calculating a fitting vector d of the forward fitF:
Finally obtaining a forward fitting coefficient alphaF:
Calculating the 1 st to M-1 st array metadata XB(k)=[x1(k),…,xM-1(k)]TThe corresponding forward sample covariance matrix,
recalculating the fitting vector d of the forward fitB:
Finally obtaining a backward fitting coefficient beta:
taking the array element data x from 1 st to M-1 sttemp(k)=[x1(k),…,xM-1(k)]TCalculating the forward fitting data of the array element No. 0,
according to the rule, the nth forward fitting data are respectively calculated in sequence:
taking the 2 nd to M array element data xtemp(k)=[x2(k),…,xM(k)]TCalculating the data of the backward M +1 number array element,
according to the rule, the nth backward fitting data are sequentially and respectively calculated:
combining the forward fitting data, the original array element receiving data and the backward fitting data together to obtain a receiving data matrix of the virtual linear array,
xLP(k)=[x1-N(k),…,x1(k),…,xM(k),…,xM+N(k)]T,
then, a sampling covariance matrix is calculated,
6. the covariance matrix reconstruction adaptive beam forming method according to claim 1, wherein a higher-dimensional interference-plus-noise covariance matrix is calculated by using snapshot data including actual array elements and virtual array elements, reconstruction errors are eliminated by adopting a mode of characteristic decomposition and large eigenvalue extraction, and finally an improved interference-plus-noise covariance matrix corresponding to the virtual array is obtained by combining an ideal noise covariance matrix; the method comprises the following specific steps:
according to the array element number of the virtual linear array, the corresponding guide vector a of the theta angle is obtainedLP(theta), calculating diagonal loading Capon space spectrums p (theta) corresponding to all the angles,
where I is the unit diagonal matrix,
according to Capon space spectrum, dividing expected signal sectorWhere B is the set region width, and then demarcating interference plus noise sectors other than the desired sectorFor the complement of Θ, the interference-plus-noise covariance matrix is initially computed in the interference-plus-noise sector
Wherein: lambda [ alpha ]m(M is 1, …, M) is a feature value arranged in descending order, v ismIs the corresponding feature vector; dintRepresenting a diagonal matrix containing all the larger eigenvalues, VintIncluding eigenvectors corresponding to larger eigenvalues, ReRepresenting an error covariance matrix consisting of smaller eigenvalues and eigenvectors, the larger eigenvalues being based onEstimating the number of interference signals as Q, adding 1 signal number as Q +1, wherein Q +1 is the minimum positive integer for ensuring the above formula to be satisfied; the number of Q interference signals arranged from large to small plus 1 signal number Q +1 forms a larger characteristic value; at this time, the interference plus noise covariance matrix is modified as:
Wherein theta islIs the l-th discrete angle within Θ, estimates the steering vector of the desired signal as,
computing beamformer weighting vectors w corresponding to virtual arraysLP:
Calculating beam output y (k) from the weight vector and the received data matrix of the virtual linear array:
7. the covariance matrix reconstruction self-adaptive beam forming system is characterized by comprising a snapshot data acquisition module, a covariance matrix calculation module,
The snapshot data acquisition module is used for acquiring time domain data received by the uniform linear array, respectively calculating forward and backward linear prediction coefficients, and sequentially fitting the data received by the virtual array elements in a recursion mode, so as to predict snapshot data of each virtual array element;
the covariance matrix calculation module is used for calculating an interference-plus-noise covariance matrix with higher dimensionality according to snapshot data containing actual array elements and virtual array elements, eliminating reconstruction errors in a mode of characteristic decomposition and large characteristic value extraction, introducing an ideal noise covariance matrix, and obtaining an improved interference-plus-noise covariance matrix corresponding to the virtual array;
the output beam calculation module is used for reconstructing an expected signal covariance matrix corresponding to the virtual array in the same way as the interference and noise covariance matrix is calculated and improved, extracting a characteristic vector corresponding to the maximum characteristic value as a guide vector of an expected signal, calculating a beam former weighting vector corresponding to the virtual array according to the guide vector, and obtaining beam output according to the weighting vector and a receiving data matrix of the virtual linear array.
8. A computer device, comprising a processor and a memory, wherein the memory is used for storing a computer executable program, the processor reads part or all of the computer executable program from the memory and executes the computer executable program, and the processor can realize the covariance matrix reconstruction adaptive beamforming method according to any one of claims 1 to 6 when executing part or all of the computer executable program.
9. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, is adapted to carry out the method of covariance matrix reconstruction adaptive beamforming according to one of claims 1-6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110673467.2A CN113422629B (en) | 2021-06-17 | 2021-06-17 | Covariance matrix reconstruction self-adaptive beam forming method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110673467.2A CN113422629B (en) | 2021-06-17 | 2021-06-17 | Covariance matrix reconstruction self-adaptive beam forming method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113422629A true CN113422629A (en) | 2021-09-21 |
CN113422629B CN113422629B (en) | 2023-02-21 |
Family
ID=77788895
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110673467.2A Active CN113422629B (en) | 2021-06-17 | 2021-06-17 | Covariance matrix reconstruction self-adaptive beam forming method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113422629B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114779176A (en) * | 2022-04-19 | 2022-07-22 | 四川大学 | Low-complexity robust adaptive beam forming method and device |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050259006A1 (en) * | 2004-05-24 | 2005-11-24 | Samsung Electronics Co., Ltd. | Beam forming apparatus and method using interference power estimation in an array antenna system |
US20050271016A1 (en) * | 2004-05-07 | 2005-12-08 | Byoung-Yun Kim | Beam forming apparatus and method for an array antenna system |
US20160014393A1 (en) * | 2014-07-09 | 2016-01-14 | Massachusetts Institute Of Technology | Methods and Apparatus for Virtual Sensor Array |
CN105306117A (en) * | 2015-10-15 | 2016-02-03 | 哈尔滨工程大学 | Para-virtual antenna array beamforming method based on covariance matrix extending |
US20160157828A1 (en) * | 2014-06-05 | 2016-06-09 | Chikayoshi Sumi | Beamforming method, measurement and imaging instruments, and communication instruments |
CN106405487A (en) * | 2016-10-17 | 2017-02-15 | 重庆邮电大学 | General spatial spectrum estimation method based on extended ESPRIT |
CN106569181A (en) * | 2016-11-15 | 2017-04-19 | 大连大学 | Algorithm for reconstructing robust Capon beamforming based on covariance matrix |
CN109600152A (en) * | 2018-12-17 | 2019-04-09 | 西北工业大学 | A kind of Adaptive beamformer method based on the transformation of subspace base |
US10574320B1 (en) * | 2019-07-15 | 2020-02-25 | Southwest Research Institute | Matrix methods to speed processing for MVDR beamforming |
CN111049556A (en) * | 2019-12-27 | 2020-04-21 | 中国科学技术大学 | Mutual prime matrix robust self-adaptive beam forming method based on interference covariance matrix reconstruction |
CN111123192A (en) * | 2019-11-29 | 2020-05-08 | 湖北工业大学 | Two-dimensional DOA positioning method based on circular array and virtual extension |
CN111651718A (en) * | 2020-05-25 | 2020-09-11 | 中国科学技术大学 | Robust adaptive beam forming method for non-Gaussian signals in Gaussian noise |
WO2020212569A1 (en) * | 2019-04-17 | 2020-10-22 | Université Du Luxembourg | Method and device for beamforming in a mimo radar system |
-
2021
- 2021-06-17 CN CN202110673467.2A patent/CN113422629B/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050271016A1 (en) * | 2004-05-07 | 2005-12-08 | Byoung-Yun Kim | Beam forming apparatus and method for an array antenna system |
US20050259006A1 (en) * | 2004-05-24 | 2005-11-24 | Samsung Electronics Co., Ltd. | Beam forming apparatus and method using interference power estimation in an array antenna system |
US20160157828A1 (en) * | 2014-06-05 | 2016-06-09 | Chikayoshi Sumi | Beamforming method, measurement and imaging instruments, and communication instruments |
US20160014393A1 (en) * | 2014-07-09 | 2016-01-14 | Massachusetts Institute Of Technology | Methods and Apparatus for Virtual Sensor Array |
CN105306117A (en) * | 2015-10-15 | 2016-02-03 | 哈尔滨工程大学 | Para-virtual antenna array beamforming method based on covariance matrix extending |
CN106405487A (en) * | 2016-10-17 | 2017-02-15 | 重庆邮电大学 | General spatial spectrum estimation method based on extended ESPRIT |
CN106569181A (en) * | 2016-11-15 | 2017-04-19 | 大连大学 | Algorithm for reconstructing robust Capon beamforming based on covariance matrix |
CN109600152A (en) * | 2018-12-17 | 2019-04-09 | 西北工业大学 | A kind of Adaptive beamformer method based on the transformation of subspace base |
WO2020212569A1 (en) * | 2019-04-17 | 2020-10-22 | Université Du Luxembourg | Method and device for beamforming in a mimo radar system |
US10574320B1 (en) * | 2019-07-15 | 2020-02-25 | Southwest Research Institute | Matrix methods to speed processing for MVDR beamforming |
CN111123192A (en) * | 2019-11-29 | 2020-05-08 | 湖北工业大学 | Two-dimensional DOA positioning method based on circular array and virtual extension |
CN111049556A (en) * | 2019-12-27 | 2020-04-21 | 中国科学技术大学 | Mutual prime matrix robust self-adaptive beam forming method based on interference covariance matrix reconstruction |
CN111651718A (en) * | 2020-05-25 | 2020-09-11 | 中国科学技术大学 | Robust adaptive beam forming method for non-Gaussian signals in Gaussian noise |
Non-Patent Citations (2)
Title |
---|
J. LI等: "Versatile robust Capon beamforming: theory and applications", 《PROCESSING WORKSHOP PROCEEDINGS, 2004 SENSOR ARRAY AND MULTICHANNEL SIGNAL》 * |
鲁欢等: "基于IAA的协方差矩阵重构稳健波束形成方法", 《太赫兹科学与电子信息学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114779176A (en) * | 2022-04-19 | 2022-07-22 | 四川大学 | Low-complexity robust adaptive beam forming method and device |
CN114779176B (en) * | 2022-04-19 | 2023-05-05 | 四川大学 | Robust self-adaptive beam forming method and device with low complexity |
Also Published As
Publication number | Publication date |
---|---|
CN113422629B (en) | 2023-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108809398B (en) | Robust adaptive beam forming method based on information source number constraint | |
CN108181507B (en) | Robust adaptive beam forming method | |
CN111049556B (en) | Mutual prime matrix robust self-adaptive beam forming method based on interference covariance matrix reconstruction | |
CN105302936B (en) | The Adaptive beamformer method reconstructed based on correlation computations and covariance matrix | |
CN110045321B (en) | Robust DOA estimation method based on sparse and low-rank recovery | |
CN108663668B (en) | IAA-based interference plus noise covariance matrix reconstruction robust beam forming method | |
CN109597046B (en) | Metric wave radar DOA estimation method based on one-dimensional convolutional neural network | |
CN109639333B (en) | Beam forming method based on effective reconstruction covariance matrix | |
CN109245814B (en) | Self-adaptive beam forming method based on maximum likelihood resampling | |
CN113422630B (en) | Self-adaptive focusing broadband beam forming method and system | |
CN110988854A (en) | Robust self-adaptive beam forming algorithm based on alternative direction multiplier method | |
CN111651718B (en) | Robust adaptive beam forming method for non-Gaussian signals in Gaussian noise | |
CN113422629B (en) | Covariance matrix reconstruction self-adaptive beam forming method and system | |
CN113673419B (en) | Beam domain rapid sparse Bayesian azimuth estimation method suitable for strong interference environment | |
CN111698011B (en) | Steady self-adaptive beam forming method for asymmetric signals in symmetric noise | |
CN111431574B (en) | Multiple steady self-adaptive beam forming method | |
CN115808659A (en) | Robust beam forming method and system based on low-complexity uncertain set integration | |
CN114047481A (en) | Robust adaptive beam forming method based on subspace orthogonality | |
CN114779176B (en) | Robust self-adaptive beam forming method and device with low complexity | |
CN107135026B (en) | Robust beam forming method based on matrix reconstruction in presence of unknown mutual coupling | |
CN114996653A (en) | Two-dimensional robust self-adaptive beam forming method based on atomic norm minimization | |
CN114563764A (en) | Self-adaptive beam generation method and system | |
CN117478186B (en) | Beam forming method based on covariance matrix estimation and subspace construction | |
KR101206297B1 (en) | Signal-blocking-based robust adaptive beamforming apparatus by interference null space projection | |
CN115494468B (en) | Clutter suppression method based on small sample sub-aperture smooth sparse recovery |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |