CN113189570B - Array signal processing method and system based on complex domain compressed sensing - Google Patents
Array signal processing method and system based on complex domain compressed sensing Download PDFInfo
- Publication number
- CN113189570B CN113189570B CN202110442884.6A CN202110442884A CN113189570B CN 113189570 B CN113189570 B CN 113189570B CN 202110442884 A CN202110442884 A CN 202110442884A CN 113189570 B CN113189570 B CN 113189570B
- Authority
- CN
- China
- Prior art keywords
- complex
- array
- domain
- processing
- array element
- 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.)
- Active
Links
- 238000003672 processing method Methods 0.000 title abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 113
- 238000012545 processing Methods 0.000 claims abstract description 62
- 239000011159 matrix material Substances 0.000 claims abstract description 48
- 238000009825 accumulation Methods 0.000 claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 27
- 238000001514 detection method Methods 0.000 claims abstract description 18
- 230000008447 perception Effects 0.000 claims abstract description 17
- 238000001228 spectrum Methods 0.000 claims description 42
- 238000004458 analytical method Methods 0.000 claims description 19
- 238000010276 construction Methods 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 6
- 230000017105 transposition Effects 0.000 claims description 5
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000003786 synthesis reaction Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000015556 catabolic process Effects 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000004304 visual acuity Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002620 method output Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention relates to the field of sonar signal processing, in particular to an array signal processing method and system based on complex domain compressed sensing. The method comprises the following steps: carrying out complex analytic transformation on array signals received by the sonar array to obtain complex analytic data of each array element, carrying out time delay compensation, correlation and accumulation processing on the complex analytic data of each array element in a complex domain according to the estimated orientation, constructing an observation sequence and a complex domain perception matrix, and realizing array signal processing by adopting a complex domain compressed sensing method. Compared with a frequency domain compressed sensing method, the method provided by the invention has the advantages that the minimum requirement on the input signal-to-noise ratio is reduced by nearly 10lgMdB (M is the number of channels) under the same detection probability, and the detection capability of a weak target is improved.
Description
Technical Field
The invention relates to the field of sonar signal processing, in particular to an array signal processing method and system based on complex domain compressed sensing.
Background
In array signal processing, in order to realize an efficient and high-performance spatial spectrum synthesis technology, researchers have conducted intensive research on space from the aspects of beam forming, subspace decomposition, compressed sensing and the like, and have obtained certain research results. However, as a matter of concern for researchers: background level, spatial resolution, requirements on input signal-to-noise ratio and the like in the spatial spectrum synthesis technology are not well solved. Although the requirement on the input signal-to-noise ratio is the lowest, the space spectrum synthesis technology based on beam forming is limited by array element number and Rayleigh limit, the problems of serious space spectrum leakage and wide main lobe are always not solved well, and the multi-target azimuth estimation effect needs to be improved; although the subspace decomposition-based spatial spectrum synthesis technology breaks through the Rayleigh limit and realizes high-resolution spatial spectrum synthesis, the subspace decomposition-based spatial spectrum synthesis technology is greatly influenced by the requirement of an input signal-to-noise ratio, cannot realize effective synthesis on a spatial spectrum under a lower signal-to-noise ratio and has a poor estimation effect on a weak target azimuth.
Compressed sensing has been widely applied to the relevant research field as an emerging theory that changes the "nyquist" sampling theory. In the target orientation estimation, on the basis of compressed sensing space target airspace sparsity, synthesizing a space spectrum by constructing a corresponding sensing matrix and a measured value, and then estimating the target orientation by the space spectrum; in subsequent combination with practical application, researchers have proposed some methods for improving the performance of compressed sensing in target orientation estimation. In general, the existing target orientation estimation methods based on compressed sensing are all realized in a frequency domain, and have certain performance degradation problems under the condition of low signal-to-noise ratio.
Disclosure of Invention
The invention aims to solve the problem of performance degradation of an array signal processing method based on frequency domain compressed sensing under the condition of low signal-to-noise ratio, and provides an array signal processing method and system based on complex domain compressed sensing.
In order to achieve the above object, the present invention provides an array signal processing method based on complex domain compressed sensing, including:
carrying out complex analytic transformation on array signals received by the sonar array to obtain complex analytic data of each array element, carrying out time delay compensation, correlation and accumulation processing on the complex analytic data of each array element in a complex domain according to the estimated orientation, constructing an observation sequence and a complex domain perception matrix, and realizing array signal processing by adopting a complex domain compressed sensing method.
As an improvement of the above method, the method specifically comprises:
step 1) carrying out complex analytic transformation on array signals received by the sonar array by adopting complex analytic wavelet transformation, and constructing complex analytic data of each array element and complex analytic data of the m-th array element in a complex domainComprises the following steps:
wherein, M is 1,2, … M, M is array element number of sonar array, xm(t) andrespectively real part data and imaginary part data, j is an imaginary part symbol, and t represents a time domain;
step 2) at the nth scanning angle thetanThe complex analysis data of the m-th array elementPush buttonPerforming time delay compensation to obtain the complex analysis data of the mth array element after the time delay compensationComprises the following steps:
wherein d is array element spacing of the sonar array, c is sound velocity, N is 1,2, … N, and N is the number of scanning angles;
step 3) constructing a covariance matrix after time delay compensation in a complex domainTo pairZeroing the main diagonal elements and aligning the zeroed covariance matrixPerforming accumulation processing to obtain corresponding spatial spectrum P (theta)n) Comprises the following steps:
Step 4) according to P (theta) ═ P (theta)1),P(θ2),…,P(θN)]T,[·]TFor matrix transposition, let I be the observation sequence, for the corresponding spatial spectrum P (θ)n) Performing transformation processing, expressing in the form of observation sequence and complex domain perception matrix, and constructing complex domain perception matrixWherein, A (theta)n) Satisfies the following formula:
step 5), obtaining a space signal sparse coefficient S (t) by solving the following convex optimization problem:
min||S(t)||1
step 6) processing the space signal sparse coefficient S (t) to obtain a synthetic space spectrum P (theta) which is as follows:
P(θ)=|S(t)|2
the peak position of the synthesized spatial spectrum P (θ) is searched for target detection.
A complex domain compressed sensing based array signal processing system, the system comprising: the method comprises the following steps: the device comprises a complex analysis transformation module, a time delay compensation and correlation accumulation module, a complex domain observation sequence and perception matrix construction module and a complex domain compressed perception processing module; wherein,
the complex analytic transformation module is used for carrying out complex analytic transformation on array signals received by the sonar array to obtain complex analytic data of each array element;
the time delay compensation and correlation accumulation module is used for carrying out time delay compensation, correlation and accumulation processing on the complex analysis data of each array element in the complex domain according to the estimated azimuth;
the complex domain observation sequence and sensing matrix construction module is used for constructing an observation sequence and a complex domain sensing matrix;
the complex domain compressed sensing processing module is used for realizing array signal processing by adopting a complex domain compressed sensing method.
As an improvement of the above system, the specific processing procedure of the complex analytic transformation module includes:
carrying out complex analytic transformation on array signals received by the sonar array by adopting complex analytic wavelet transformation, and constructing complex analytic data of each array element and complex analytic data of the m-th array element in a complex domainComprises the following steps:
wherein, M is 1,2, … M, M is array element number of sonar array, xm(t) andreal and imaginary data, respectively, j is the imaginary symbol, and t represents the time domain.
As an improvement of the above system, the specific processing procedure of the delay compensation and correlation accumulation module includes:
at the nth scan angle thetanThe complex analysis data of the m-th array elementPush buttonPerforming time delay compensation to obtain the complex analysis data of the mth array element after the time delay compensationComprises the following steps:
wherein d is array element spacing of the sonar array, c is sound velocity, N is 1,2, … N, and N is the number of scanning angles;
constructing a time-delay compensated covariance matrix in a complex domainTo pairZeroing the main diagonal elements and aligning the zeroed covariance matrixPerforming accumulation processing to obtain corresponding spatial spectrum P (theta)n) Comprises the following steps:
As an improvement of the above system, the specific processing procedure of the complex-domain observation sequence and perception matrix construction module includes:
according to P (theta) ═ P (theta)1),P(θ2),…,P(θN)]T,[·]TFor matrix transposition, let I be the observation sequence, for the corresponding spatial spectrum P (θ)n) Performing transformation processing, expressing in the form of observation sequence and complex domain perception matrix, and constructing complex domain perception matrixWherein, A (theta)n) Satisfies the following formula:
as an improvement of the above system, the specific processing procedure of the complex domain compressed sensing processing module includes:
obtaining the spatial signal sparsity coefficient s (t) by solving the following convex optimization problem:
min||S(t)||1
processing the space signal sparse coefficient S (t) to obtain a synthetic space spectrum P (theta) as follows:
P(θ)=|S(t)|2
the peak position of the synthesized spatial spectrum P (θ) is searched for target detection.
Compared with the prior art, the invention has the advantages that:
1. according to the method, correlation and accumulation processing are carried out on array data in a complex domain, the signal-to-noise ratio contained in data at each position of a sensing matrix is improved, so that the data at each position of the sensing matrix has certain array gain, and compared with an observation sequence used in an array signal processing method based on frequency domain compressed sensing, the signal-to-noise ratio contained in the data at each position of the sensing matrix is improved, and further the minimum requirement of the frequency domain compressed sensing method on the input signal-to-noise ratio is improved;
2. compared with a frequency domain compressed sensing method, the method provided by the invention has the advantages that the minimum requirement on the input signal-to-noise ratio is reduced by nearly 10lgMdB (M is the number of channels) under the same detection probability, and the detection capability of a weak target is improved.
Drawings
Fig. 1 is a schematic view of a horizontal towed-line array sonar structure adopted in embodiment 1 of the present invention;
FIG. 2 is a flow chart of the array signal processing method based on complex domain compressed sensing of the present invention;
FIG. 3 is a graph showing the probability of detection of a target obtained by 200 independent statistics using the MVDR method, the MUSIC method, the FCS method and the CCS method of the present invention, respectively;
FIG. 4 shows the spatial spectra (SNR ═ 0dB) obtained by the MVDR method, the MUSIC method, the FCS method and the CCS method of the present invention, respectively;
FIG. 5 shows the spatial spectrum (SNR-10 dB) obtained by the MVDR method, the MUSIC method, the FCS method and the CCS method of the present invention, respectively;
FIG. 6 shows the spatial spectrum (SNR-15 dB) obtained by the MVDR method, the MUSIC method, the FCS method and the CCS method of the present invention, respectively;
FIG. 7 is an azimuth history graph output using the MVDR method;
FIG. 8 is a diagram of an azimuth history output using the MUSIC method;
FIG. 9 is a chart of azimuth history output using the FCS method;
FIG. 10 is a chart of azimuth history output using the CCS method of the present invention.
Detailed Description
Aiming at the problem of performance degradation of an array signal processing method based on frequency domain compressed sensing under the condition of low signal-to-noise ratio, the invention provides an array signal processing method based on complex domain compressed sensing according to the relation of phase shift and time delay in array signal processing. The method comprises the steps of carrying out complex analytic transformation on received signals of the linear array, carrying out time delay compensation, correlation and accumulation processing on each array element signal in a complex domain according to an estimated azimuth, constructing a complex domain sensing matrix and a measured value, and realizing array signal processing by adopting a complex domain compression sensing method. The numerical simulation and actual measurement data processing results show that compared with a frequency domain compressed sensing method, the method has the advantages that the minimum requirement on the input signal-to-noise ratio is reduced by nearly 10lgMdB (M is the number of channels) under the same detection probability, and the detection capability of a weak target is improved.
The technical solution of the present invention will be described in detail below with reference to the accompanying drawings and examples.
Example 1
Before describing the method of the present invention in detail, the horizontal drawing line array used in the method of the present invention will be described. Fig. 1 is a schematic structural view of a horizontal towed array sonar, which comprises 6 parts, a display control and signal processor 1, a deck cable 2, a winch 3, a cable guide frame 4, a towing cable 5 and a horizontal towed array 6. The horizontal towing line array 6 is connected with a deck cable 2 on a winch 3 through a towing cable 5, and the towing cable 5 is also arranged on a cable guide frame 4; and the signal received by the horizontal dragging linear array 6 is transmitted to the display control and signal processor 1. It should be noted that the method of the present invention is not limited to horizontal towed linear arrays, and may be sonar arrays with multiple array elements.
The process of the present invention is further illustrated below.
Frequency domain compressed sensing method
Order: the frequency domain expression of the M-element equidistant linear array receiving signals is as follows:
X(f)=A(f,θ)S(f)+V(f) (1)
wherein X (f) ═ X1(f),X2(f),…,XM(f)]TIn the form of a frequency domain vector of the linear array of received signals, S (f) ([ S ]1(f),S2(f),…,SK(f)]TIn the form of a frequency domain vector of K target radiation signals in space, V (f) ═ V1(f),V2(f),…,VM(f)]TIn the form of a frequency domain vector of noise contained in a linear array reception signal, a (f, θ) ═ a (f, θ)1),a(f,θ2),…,a(f,θK)]Is a linear array manifold, thetak(K1, 2, … K) e theta is the orientation of the kth target relative to the linear array,is a direction vector. Wherein [ ·]TAnd c represents the matrix transfer operator, c represents the speed, and d represents the distance between adjacent array elements.
If the entire space is divided into sufficiently small directional intervalsAnd assuming every possible directionAll corresponding to a potential target signalThus, N target signals are constructedComplete array manifold constructed at this timeWherein,
the frequency domain wideband model shown in equation (1) can be expressed as:
where J is 1,2, …, and J is the band number.
It is clear that the following description of the preferred embodiments,Sa(fj) Only the corresponding direction in (theta)1,θ2,…θNThe target energy on is large, while the other direction is a sufficiently small value, i.e. Sa(fj) Is a sparse representation of the frequency domain of the signal space. Comparing the compressed sensing models, if anyTo be seen as a sequence of observations,as a sensing matrix, Sa(fj) For the sparse coefficient component to be solved, Va(fj) To measure noise, the spatial signal sparsity coefficient S can be solved by solving the following convex-down optimization problema(fj):
At this time, the passing type(4) Calculated Sa(fj) I.e. the spatial spectrum estimation value of the j-th sub-band.
P(fj,θ)=|Sa(fj)|2 (4)
Repeating the above process for all sub-bands and summing, so as to obtain a total space spectrum P (theta), and searching the peak position of the space spectrum P (theta) to realize target azimuth estimation.
The method is an array signal processing method based on frequency domain compressed sensing.
Complex domain compressed sensing method
Basic principle
According to the array signal processing method model based on frequency domain compressed sensing, when an observation sequence and a sensing matrix are constructed by using single frequency domain data, the signal-to-noise ratio contained in the array element domain data is enhanced only by using the self information of each array element and not by using the correlation difference characteristic between the signals and the noises of each array element. In contrast, the present invention provides a complex domain processing method, i.e., an array signal processing method based on complex domain compressed sensing, in which correlation and accumulation processing are performed on array data in a complex domain, so as to improve the signal-to-noise ratio of data at each position of a sensing matrix, as shown in fig. 2.
Firstly, complex analytic wavelet transform is adopted to carry out complex analytic transform on array data, and complex analytic data of each array element signal x (t) are constructed in a complex domain Is x (t) the corresponding imaginary data.
Second, at a scan angle θnIn the above, the data is repeatedly analyzed for each array elementPush buttonAnd carrying out time delay compensation to obtain data after time delay compensation.
In the formula, d is the array element spacing, c is the sound velocity, N is 1,2, … N, N is the number of scanning angles, M is 1,2, … M, and M is the number of array elements.
Constructing a time-delay compensated covariance matrix in a complex domainTo pairZeroing the main diagonal elements and aligning the zeroed sum-pair covariance matricesAnd accumulating to obtain a corresponding spatial spectrum as follows:
in the formula,I=[1,1,…,1]implementation ofAnd (5) performing medium element accumulation processing.
Then, equation (7) is transformed and expressed in the form of an observation sequence or a perceptual matrix, that is, I is regarded as an observation sequence, and P (θ) ([ P (θ) ]1),P(θ2),…,P(θN)]T,[·]TFor matrix transposition, for sparse coefficient components to be solved, a complex domain sensing matrix is usedThe design is as follows:
finally, the spatial signal sparsity coefficient s (t) is solved by solving the following convex optimization problem:
and (4) processing the S (t) to obtain a synthesized space spectrum P (theta) of the method, and searching the peak position of the space spectrum P (theta) to realize target detection.
P(θ)=|S(t)|2 (10)
Method capability analysis of the invention
In the processing frequency band, the correlation between the space target signals and the correlation between the background noises are made to be zero, and the power spectrum of the background noises is made to beSignal power ofThenEach location data may be represented as:
in addition, the array signal processing method based on frequency domain compressed sensing is used for forming observation sequences, and the observation sequencesEach position data contains a signal-to-noise ratio of:
also, as can be seen from equation (11), after complex domain processing,each position data contains a signal-to-noise ratio of:
comparing the equations (12) and (13), the array signal processing method based on complex domain compressed sensing makes each position data of the sensing matrix have a certain array gain by correlating and accumulating each array element data, and improves the signal-to-noise ratio of each position data of the sensing matrix compared with the observation sequence used by the array signal processing method based on frequency domain compressed sensing, thereby improving the minimum requirement of the array signal processing method based on frequency domain compressed sensing on the input signal-to-noise ratio.
The effect of the process of the invention is compared below with reference to example 1.
In order to further verify that the method can effectively realize array signal processing under the condition of low signal-to-noise ratio, and realize detection on the target, the following numerical simulation analysis is carried out.
The following table is a numerical simulation parameter set.
TABLE 1
FIG. 3 shows the probability of target detection obtained by 200 independent statistics of the MVDR method, the MUSIC method, the array signal processing method based on frequency domain compressed sensing (referred to as FCS method in the present invention) and the method of the present invention (referred to as CCS method in the present invention) under the condition that the SNR is-24 dB-0 dB
As can be seen from the results in fig. 3, compared with the FCS method, the CCS method performs correlation and accumulation processing on each array element signal, so that effective target detection is achieved under the condition of low signal-to-noise ratio; compared with the FCS method, the CCS method reduces the minimum requirement on the input signal-to-noise ratio by 6 dB.
FIG. 4 shows the spatial spectra (SNR ═ 0dB) obtained by the MVDR method, the MUSIC method, the FCS method and the CCS method;
FIG. 5 shows the spatial spectra (SNR-10 dB) obtained by MVDR method, MUSIC method, FCS method, and CCS method;
fig. 6 shows the spatial spectra (SNR-15 dB) obtained by the MVDR method, the MUSIC method, the FCS method, and the CCS method.
From the simulation results, it can be obtained: compared with the MVDR method and the MUSIC method, the FCS method and the CCS method can realize high-resolution detection and resolution on adjacent targets; however, as the signal-to-noise ratio is reduced, the resolving power of the MVDR method and the MUSIC method is reduced more severely, under the condition that the signal-to-noise ratio is-10 dB, the MVDR method and the MUSIC method cannot achieve the resolution of two targets, and the FCS method and the CCS method also maintain the target resolving power under the condition of high signal-to-noise ratio, but compared with the condition of high signal-to-noise ratio, the spatial spectrum obtained by the FCS method is seriously polluted by noise, and effective detection of the two targets cannot be achieved under the condition of-15 dB.
The effect of the process of the invention is compared below with reference to example 2.
A plurality of targets exist in the data processing time period, wherein the relative array orientations of 30 degrees, 60 degrees, 80 degrees, 100 degrees, 130 degrees and 140 degrees are ships with large tonnage, and the sound level is about 120 dB-130 dB @1 kHz.
The 4 methods were performed as described in example analysis 1. Fig. 7 to 10 are output time history charts of 4 methods. Specifically, the method comprises the following steps:
FIG. 7 is a diagram of an output azimuth history of the MVDR method;
FIG. 8 is a diagram of the output azimuth history of the MUSIC method;
FIG. 9 is a diagram of the output azimuth history of the FCS method;
FIG. 10 is a diagram of the CCS method output bearing history.
As can be seen from the results shown in fig. 7 to 10, the spatial spectrum obtained by the MUSIC method and the FCS method cannot effectively detect the 130 ° target within a time period of 0 to 180s, and the azimuth history chart formed by the spatial spectrum obtained by the CCS method can clearly show the target tracks (carrying platforms) of 30 °, 60 °, 80 °, 100 °, 130 °, 140 ° and 160 °, so that the target azimuth is clearly recognizable; although the MVDR method can effectively detect a plurality of targets in the time interval, the detection effect of the target azimuth at 20 degrees is inferior to that of the CCS method due to the unknown pulse target at 40 degrees, and the background level is much higher than that of the CCS method.
The data processing result verifies that the CCS method can effectively detect the target under the condition of low signal-to-noise ratio.
Example 2
Based on the above method, embodiment 2 of the present invention provides an array signal processing system based on complex domain compressed sensing, where the system includes: the method comprises the following steps: the device comprises a complex analysis transformation module, a time delay compensation and correlation accumulation module, a complex domain observation sequence and perception matrix construction module and a complex domain compressed perception processing module; wherein,
the complex analytic transformation module is used for carrying out complex analytic transformation on array signals received by the sonar array to obtain complex analytic data of each array element;
the time delay compensation and correlation accumulation module is used for carrying out time delay compensation, correlation and accumulation processing on the complex analysis data of each array element in the complex domain according to the estimated azimuth;
the complex domain observation sequence and sensing matrix construction module is used for constructing an observation sequence and a complex domain sensing matrix;
the complex domain compressed sensing processing module is used for realizing array signal processing by adopting a complex domain compressed sensing method.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention and are not limited. Although the present invention has been described in detail with reference to the embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the spirit and scope of the invention as defined in the appended claims.
Claims (2)
1. A method for array signal processing based on complex domain compressed sensing, the method comprising:
carrying out complex analytic transformation on array signals received by the sonar array to obtain complex analytic data of each array element, carrying out time delay compensation, correlation and accumulation processing on the complex analytic data of each array element in a complex domain according to the estimated orientation, constructing an observation sequence and a complex domain sensing matrix, and realizing array signal processing by adopting a complex domain compression sensing method;
the method specifically comprises the following steps:
step 1) carrying out complex analytic transformation on array signals received by the sonar array by adopting complex analytic wavelet transformation, and constructing complex analytic data of each array element and complex analytic data of the m-th array element in a complex domainComprises the following steps:
wherein, M is 1,2, … M, M is array element number of sonar array, xm(t) andrespectively real part data and imaginary part data, j is an imaginary part symbol, and t represents a time domain;
step 2) at the nth scanning angle thetanThe complex analysis data of the m-th array elementPush buttonPerforming time delay compensation to obtain a complex solution of the m-th array element after time delay compensationAnalyzing dataComprises the following steps:
wherein d is array element spacing of the sonar array, c is sound velocity, N is 1,2, … N, and N is the number of scanning angles;
step 3) constructing a covariance matrix after time delay compensation in a complex domainTo pairZeroing the main diagonal elements and aligning the zeroed covariance matrixPerforming accumulation processing to obtain corresponding spatial spectrum P (theta)n) Comprises the following steps:
step 4) according to P (theta) ═ P (theta)1),P(θ2),…,P(θN)]T,[·]TFor matrix transposition, let I be the observation sequence, for the corresponding spatial spectrum P (θ)n) Performing transformation processing, expressing in the form of observation sequence and complex domain perception matrix, and constructing complex domain perceptionMatrix arrayWherein, A (theta)n) Satisfies the following formula:
step 5), obtaining a space signal sparse coefficient S (t) by solving the following convex optimization problem:
min||S(t)||1
step 6) processing the space signal sparse coefficient S (t) to obtain a synthetic space spectrum P (theta) which is as follows:
P(θ)=|S(t)|2
the peak position of the synthesized spatial spectrum P (θ) is searched for target detection.
2. An array signal processing system based on complex domain compressed sensing, the system comprising: the device comprises a complex analysis transformation module, a time delay compensation and correlation accumulation module, a complex domain observation sequence and perception matrix construction module and a complex domain compressed perception processing module; wherein,
the complex analytic transformation module is used for carrying out complex analytic transformation on array signals received by the sonar array to obtain complex analytic data of each array element;
the time delay compensation and correlation accumulation module is used for carrying out time delay compensation, correlation and accumulation processing on the complex analysis data of each array element in the complex domain according to the estimated azimuth;
the complex domain observation sequence and sensing matrix construction module is used for constructing an observation sequence and a complex domain sensing matrix;
the complex domain compressed sensing processing module is used for realizing array signal processing by adopting a complex domain compressed sensing method;
the specific processing procedure of the complex analysis transformation module comprises the following steps:
carrying out complex analytic transformation on array signals received by the sonar array by adopting complex analytic wavelet transformation, and constructing complex analytic data of each array element and complex analytic data of the m-th array element in a complex domainComprises the following steps:
wherein, M is 1,2, … M, M is array element number of sonar array, xm(t) andrespectively real part data and imaginary part data, j is an imaginary part symbol, and t represents a time domain;
the specific processing procedure of the delay compensation and correlation accumulation module comprises the following steps:
at the nth scan angle thetanThe complex analysis data of the m-th array elementPush buttonPerforming time delay compensation to obtain the complex analysis data of the mth array element after the time delay compensationComprises the following steps:
wherein d is array element spacing of the sonar array, c is sound velocity, N is 1,2, … N, and N is the number of scanning angles;
construction of epochs in the complex domainDelay compensated covariance matrixTo pairZeroing the main diagonal elements and aligning the zeroed covariance matrixPerforming accumulation processing to obtain corresponding spatial spectrum P (theta)n) Comprises the following steps:
the specific processing process of the complex domain observation sequence and perception matrix construction module comprises the following steps:
according to P (theta) ═ P (theta)1),P(θ2),…,P(θN)]T,[·]TFor matrix transposition, let I be the observation sequence, for the corresponding spatial spectrum P (θ)n) Performing transformation processing, expressing in the form of observation sequence and complex domain perception matrix, and constructing complex domain perception matrixWherein, A (theta)n) Satisfies the following formula:
the specific processing procedure of the complex domain compressed sensing processing module comprises the following steps:
obtaining the spatial signal sparsity coefficient s (t) by solving the following convex optimization problem:
min||S(t)||1
processing the space signal sparse coefficient S (t) to obtain a synthetic space spectrum P (theta) as follows:
P(θ)=|S(t)|2
the peak position of the synthesized spatial spectrum P (θ) is searched for target detection.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110442884.6A CN113189570B (en) | 2021-04-23 | 2021-04-23 | Array signal processing method and system based on complex domain compressed sensing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110442884.6A CN113189570B (en) | 2021-04-23 | 2021-04-23 | Array signal processing method and system based on complex domain compressed sensing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113189570A CN113189570A (en) | 2021-07-30 |
CN113189570B true CN113189570B (en) | 2022-02-01 |
Family
ID=76978251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110442884.6A Active CN113189570B (en) | 2021-04-23 | 2021-04-23 | Array signal processing method and system based on complex domain compressed sensing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113189570B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114036975B (en) * | 2021-10-19 | 2022-05-17 | 中国科学院声学研究所 | Target signal extraction method based on frequency domain-wavenumber domain deconvolution |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5481503A (en) * | 1984-12-12 | 1996-01-02 | Martin Marietta Corporation | Apparatus for and method of adaptively processing sonar data |
CN101799535A (en) * | 2009-11-27 | 2010-08-11 | 西安电子科技大学 | Method for estimating target direction by multiple input multiple output (MIMO) radar |
CN112684408A (en) * | 2020-11-06 | 2021-04-20 | 西安交通大学 | Near-field signal source positioning method, system and equipment based on sparse Bayesian learning |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ATE549602T1 (en) * | 2002-01-23 | 2012-03-15 | Cidra Corporate Services Inc | DEVICE AND METHOD FOR MEASURING PARAMETERS OF A MIXTURE WITH SOLID PARTICLES SUSPENDED IN A FLUID FLOWING IN A PIPE |
ATE433241T1 (en) * | 2003-06-18 | 2009-06-15 | Panasonic Corp | METHOD AND DEVICE FOR RECEIVING DIGITAL MULTI CARRIER SIGNALS USING A WAVELET TRANSFORMATION |
US9097783B2 (en) * | 2006-04-28 | 2015-08-04 | Telecommunication Systems, Inc. | System and method for positioning using hybrid spectral compression and cross correlation signal processing |
US10310075B2 (en) * | 2016-04-18 | 2019-06-04 | The Boeing Company | Synthetic aperture radar data compression and transmission |
CN106546846B (en) * | 2016-10-18 | 2019-12-10 | 天津大学 | Electric energy quality signal detection device based on compressed sensing blind source signal separation technology |
CN106772225B (en) * | 2017-01-20 | 2019-03-26 | 大连大学 | Compressed sensing based Beam Domain DOA estimation |
US11806155B2 (en) * | 2017-05-11 | 2023-11-07 | Sridhar (sri) KRISHNAN | System and method for QRS complex detection in compressively sensed electrocardiogram data |
CN107480471B (en) * | 2017-07-19 | 2020-09-01 | 福建师范大学 | Sequence similarity analysis method based on wavelet transform characteristics |
CN109444810B (en) * | 2018-12-24 | 2022-11-01 | 哈尔滨工程大学 | Mutual-prime array non-grid DOA estimation method under nonnegative sparse Bayesian learning framework |
CN109995686A (en) * | 2019-03-25 | 2019-07-09 | 哈尔滨工程大学 | A kind of sparse underwater acoustic channel estimation method of complex field |
CN111693960A (en) * | 2020-06-11 | 2020-09-22 | 中山大学 | Signal processing method of variable-frequency compressed sensing radar |
-
2021
- 2021-04-23 CN CN202110442884.6A patent/CN113189570B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5481503A (en) * | 1984-12-12 | 1996-01-02 | Martin Marietta Corporation | Apparatus for and method of adaptively processing sonar data |
CN101799535A (en) * | 2009-11-27 | 2010-08-11 | 西安电子科技大学 | Method for estimating target direction by multiple input multiple output (MIMO) radar |
CN112684408A (en) * | 2020-11-06 | 2021-04-20 | 西安交通大学 | Near-field signal source positioning method, system and equipment based on sparse Bayesian learning |
Also Published As
Publication number | Publication date |
---|---|
CN113189570A (en) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109212500B (en) | High-precision KA-STAP (K-ary-based adaptive-noise) covariance matrix estimation method based on sparse reconstruction | |
US8179300B2 (en) | Method for suppressing clutter in space-time adaptive processing systems | |
CN102323583B (en) | Super-resolution linear array three-dimensional synthetic aperture radar imaging method | |
CN110058193B (en) | Digital multi-beam angle measurement method and system based on single receiving channel | |
CN103383449B (en) | ESPRIT algorithm based short-range clutter suppression method for airborne radar | |
CN113376603B (en) | Sub-band space-time adaptive processing method of broadband airborne phased array radar | |
CN106483516A (en) | Radar clutter space-time adaptive processing method based on priori | |
CN103116162B (en) | High-resolution sonar location method based on sparsity of objective space | |
CN108931766B (en) | Non-uniform STAP interference target filtering method based on sparse reconstruction | |
CN111913157A (en) | Sea clutter suppression method based on radar signal space-time decorrelation model | |
CN109001687A (en) | Airborne radar space-time self-adaptive filtering method based on generalized sidelobe cancellation structure | |
CN111413666A (en) | Design method of array frequency and direction measurement combined receiver | |
CN107886068A (en) | A kind of unknown line spectrum object detection method based on time-space domain line enhancement device | |
CN112462363B (en) | Non-uniform sparse polarization array coherent target parameter estimation method | |
CN105911527B (en) | Airborne radar space-time adaptive processing method based on EFA and MWF | |
CN103197297A (en) | Radar moving target detection method based on cognitive framework | |
CN105572668B (en) | A kind of moving target radial velocity method of estimation converted based on Radon | |
CN113189570B (en) | Array signal processing method and system based on complex domain compressed sensing | |
CN112255595A (en) | Sea clutter data preprocessing method based on simulated airborne measurement | |
CN108896963B (en) | Airborne radar space-time self-adaptive dimension reduction processing method | |
CN106802408A (en) | Airborne non-positive side battle array short range clutter range ambiguity suppressing method based on sparse recovery | |
CN111175727A (en) | Method for estimating orientation of broadband signal based on conditional wave number spectral density | |
CN109061597B (en) | Ionospheric clutter suppression method based on blind source separation and time-frequency ridge domain filtering | |
CN105572642A (en) | Space-time adaptive processing method based on two-stage architecture | |
CN113109760B (en) | Multi-line spectrum combined DOA estimation and clustering method and system based on group sparsity |
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 |