CN114236516A - Deconvolution beam forming acceleration method based on R-L iterative algorithm - Google Patents

Deconvolution beam forming acceleration method based on R-L iterative algorithm Download PDF

Info

Publication number
CN114236516A
CN114236516A CN202111190318.7A CN202111190318A CN114236516A CN 114236516 A CN114236516 A CN 114236516A CN 202111190318 A CN202111190318 A CN 202111190318A CN 114236516 A CN114236516 A CN 114236516A
Authority
CN
China
Prior art keywords
psf
deconvolution
natural directivity
beam forming
beams
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.)
Pending
Application number
CN202111190318.7A
Other languages
Chinese (zh)
Inventor
张铮
徐雅南
孔强
王伟
李恒光
范泽亚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
715th Research Institute of CSIC
Original Assignee
715th Research Institute of CSIC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 715th Research Institute of CSIC filed Critical 715th Research Institute of CSIC
Priority to CN202111190318.7A priority Critical patent/CN114236516A/en
Publication of CN114236516A publication Critical patent/CN114236516A/en
Pending legal-status Critical Current

Links

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a deconvolution beam forming acceleration method based on an R-L iterative algorithm, which utilizes array type related parameters to calculate a natural directivity function PSF of a uniform linear arrayR＝{PSF1，PSF2，...，PSFrR represents the total number of frequency point numbers; dividing sub-bands with unequal intervals according to the bandwidth, and extracting natural directivity functions of central frequency points of the corresponding sub-bands; performing conventional beam forming calculation of m beams by using array metadata of the uniform line array to obtain power spectrums P of the m beamsM＝{P1，P2，...，Pm}; natural directivity function PSF after subband central frequency point extractionO＝{PSF1，PSF2，...，PSFoPower spectrum P of M wave beamsM＝{P1，P2，...，PmR-L deconvolution calculation is carried out, where o denotes the total number of divided subbands, oR is less than r; and obtaining m wave beam results as a final direction estimation result by using the R-L deconvolution calculation result. The invention overcomes the defects of wide main lobe, high side lobe, poor weak target detection capability and the like of the conventional beam forming, and simultaneously solves the problem of overlarge calculated amount when the traditional deconvolution beam forming is used for processing a broadband random signal.

Description

Deconvolution beam forming acceleration method based on R-L iterative algorithm
Technical Field
The invention relates to the field of underwater acoustic array signal processing, in particular to a deconvolution beam forming acceleration method based on an R-L iterative algorithm.
Background
The horizontal line array is an important component of the sonar system as a basic array. Array gain is obtained through Conventional Beam Forming (CBF) processing, and target detection capability is improved, but CBF main lobe beam width is wide, side lobe level is high, spatial resolution capability is reduced, aliasing is easily generated during multi-target processing, and weak target detection capability under strong signal interference is weak. Deconvolution beam forming is a post-processing method of CBF, which can give consideration to both high resolution capability and stability, but when the bandwidth of the processed signal is wide, the calculation amount is large, and the operation complexity is high. In addition, the deconvolution algorithm based on the Richardson-Lucy iterative method (R-L) has the defect that the boundary blurring appears at the two ends of a beam pattern due to signal truncation, so that the practical engineering application is influenced. On the basis of the traditional deconvolution beam forming, the method of multi-band processing and energy spectrum boundary expansion is adopted, and the higher resolution is obtained, and the forming speed is improved by 50% compared with the traditional deconvolution beam forming.
Disclosure of Invention
One of the purposes of the invention is to provide a deconvolution beam forming acceleration method based on an R-L iterative algorithm, so as to solve the defects of the prior conventional beam forming in the background technology, such as wide main lobe, high side lobe, poor weak target detection capability and the like.
The invention also aims to provide a deconvolution beam forming acceleration method based on the R-L iterative algorithm, and simultaneously solves the problem that the traditional deconvolution beam forming has overlarge calculation amount when processing broadband random signals.
In order to achieve the purpose, the invention provides the following technical scheme:
a deconvolution beam forming acceleration method based on an R-L iterative algorithm comprises the following steps:
natural directivity function PSF (pseudo-random function) for calculating uniform linear array by using array type related parametersR＝{PSF1，PSF2，...，PSFrR represents the total number of frequency point numbers;
dividing sub-bands with unequal intervals according to the bandwidth, and extracting natural directivity functions of central frequency points of the corresponding sub-bands;
performing conventional beam forming calculation of m beams by using array metadata of the uniform line array to obtain power spectrums P of the m beamsM＝{P1，P2，…，Pm}；
Natural directivity function PSF after subband central frequency point extractionO＝{PSF1，PSF2，…，PSFoPower spectrum P of M wave beamsM＝{P1，P2，…，PmPerforming R-L deconvolution calculation, wherein o represents the total number of divided sub-bands, and o is less than R;
and obtaining m wave beam results as a final direction estimation result by using the R-L deconvolution calculation result.
Preferably, the power spectrum P (cos theta) is calculated by
Figure RE-GDA0003511540100000021
Wherein the content of the first and second substances,
Figure RE-GDA0003511540100000022
Figure RE-GDA0003511540100000023
s () is the source distribution function, Bp() As a beam directivity function, K is the number of sound signals, i is the ith sound signal, AiIs the intensity of the ith sound signal, θiIs the ithThe azimuth of incidence of the individual sound signals,
Figure RE-GDA0003511540100000024
and delta () is an impulse function, N is the number of array elements, d is the array element spacing, and lambda is the wavelength.
Preferably, the natural directivity function is calculated by
Figure RE-GDA0003511540100000031
Wherein T is signal frequency, N is the number of array elements, d is the spacing of the array elements, c is the sound velocity in water, theta is the incident azimuth angle of the signal,
Figure RE-GDA0003511540100000032
is the differential value of the beam angle.
Preferably, the power spectrums P of m wave beams are obtained after natural directivity functions of corresponding sub-band central frequency points are obtainedMThen, power spectra P of m wave beams are respectivelyMAnd a natural directivity function PSF after sub-band central frequency point extractionOPerforming period expansion to obtain an expanded power spectrum PM′And an extended natural directivity function PSFO′(ii) a For a periodically extended power spectrum PM′And a periodically extended natural directivity function PSFO′Performing R-L deconvolution calculation; and taking the middle m beam results in the R-L deconvolution results as final azimuth estimation results.
Preferably, said power spectrum P for m beamsMThe cycle extension comprises the following steps:
by a power spectrum PMThe left end point and the right end point of the beam are used as starting points, and the preset number l of beams is respectively expanded outwards;
spreading the beam to power spectrum PMTaking the left end point and the right end point as reference values to obtain a power spectrum P with the period expansionM′＝{P1，P2，…，Pm，…，Pm+2l}，PmRepresenting the power value of the mth beam.
Preferably, the extracted subband central frequency point is the subband central frequency pointDirectivity function PSFOThe cycle extension comprises the following steps:
natural directivity function PSF with subband central frequency pointβ＝{PSF′1，PSF′2，…，PSF′m}βThe left and right endpoints of (b) are starting points, and a preset number l of beams are respectively extended outwards, wherein β represents the number of the subband, and β is {1, 2, 3, …, o }, PSF'mThe natural directivity of the mth wave beam under a certain fixed sub-band central frequency point is shown;
natural directivity function PSF of extended wave beam with sub-band center frequency pointβTaking the left and right endpoints as reference values to obtain the expanded natural directivity function PSFβ′＝{PSF′1，PSF′2，...，PSF′m，...，PSF′m+2l}β；
Obtaining a periodically extended natural directivity function PSFO＝{PSF1′，PSF2′，...，PSFβ′，...，PSFo′}。
Preferably, the preset number l is 0.25 m.
Preferably, the power of the expanded beam is the power spectrum PMA power value corresponding to the endpoint; the natural directivity of the expanded beam is the natural directivity of the corresponding natural directivity function corresponding to the end point.
Preferably, the power of the beam is spread to a power spectrum PMTaking the reference axis of the corresponding end point as a symmetry axis, and carrying out symmetry value taking; and the natural directivity of the expanded wave beam takes the reference axis of the corresponding endpoint of the corresponding natural directivity function as a symmetry axis to carry out symmetrical value taking.
Compared with the prior art, the invention has the beneficial effects that:
the invention overcomes the defects of wide main lobe, high side lobe, poor weak target detection capability and the like of the conventional beam forming, and simultaneously solves the problem of overlarge calculated amount when the traditional deconvolution beam forming is used for processing a broadband random signal. The operation amount of deconvolution is reduced through sub-band division processing; by the periodic expansion of the power spectrum, the problem of boundary ambiguity caused by the R-L iterative algorithm is solved, and meanwhile, the calculation speed is further improved.
Drawings
Fig. 1 is a flowchart of an R-L iterative algorithm-based deconvolution beam forming acceleration method of embodiment 2.
FIG. 2 is a schematic diagram of the arrangement of a horizontal uniform array.
Fig. 3 is a graph of the PSF of a horizontal uniform linear array at each frequency point.
Fig. 4 is a graph of deconvolved beam output after molecular band decimation.
Fig. 5 is a graph of deconvolved beam output after beam spectrum boundary extension.
Fig. 6 is a table for increasing the speed of deconvolution beam forming based on the R-L iterative algorithm in comparison with the conventional deconvolution beam forming speed in embodiment 2.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without any inventive step, are within the scope of the present invention.
Example 1:
referring to fig. 1, the deconvolution beam forming acceleration method based on the R-L iterative algorithm is mainly applied to a uniform linear array, taking a horizontal linear array with N array elements as an example, the horizontal uniform linear array is shown in fig. 2, a plurality of array elements are arranged at equal levels and intervals, the distance between every two adjacent array elements is d, a sound source is in a far field, an acoustic signal reaches each array element in a plane wave form, and the incident direction of the signal is theta. The deconvolution beam forming acceleration method based on the R-L iterative algorithm comprises the following steps.
PSF calculation step: natural directivity function PSF (pseudo-random function) for calculating uniform linear array by using array type related parametersR＝{PSF1，PSF2，...，PSFrR represents the total number of frequency point numbers.
Natural directivity function (i.e. in deconvolution algorithm)The point spread function PSF) is a parameter related to the signal frequency, and the specific calculation formula is:
Figure RE-GDA0003511540100000061
wherein f is the signal frequency, N is the number of array elements, d is the array element spacing, c is the sound velocity in water, theta is the signal incidence azimuth angle,
Figure RE-GDA0003511540100000062
is the differential value of the beam angle theta.
In the invention, the frequency points are numbers for fixed frequencies. If the frequency intervals are all 1Hz, the frequency intervals are divided into 6001 wireless frequency bands from 0Hz, 1MHz, 2Hz, 3Hz, 4Hz and 5Hz … … 6000Hz according to the frequency intervals of 1Hz, and each frequency band is numbered from 1, 2, 3 and 4 … … 6001; these numbers for fixed frequencies are frequency bins. The method calculates the natural directivity function of the uniform horizontal line array at each frequency point through a calculation formula of the natural directivity function, and finally collects the natural directivity functions of all the frequency points to obtain the natural directivity function PSF of the uniform horizontal line arrayR＝{PSF1，PSF2，…，PSFr}，PSFrThe meaning of (1) is m natural directivities corresponding to the m wave beams respectively at the r frequency point.
Molecular band extraction: and dividing sub-bands with unequal intervals according to the bandwidth, and extracting the natural directivity function of the central frequency point of the corresponding sub-band.
The bandwidth division in the invention can be performed according to the actual situation, for example, the total analysis bandwidth is 0-6kHz, and the division is (0-2000) · (2000-.
In the present invention, bandwidth generally refers to the frequency bandwidth occupied by a signal; when used to describe a channel, bandwidth refers to the maximum frequency bandwidth of a signal that can effectively pass through the channel. For analog signals, bandwidth is also known as bandwidth, in hertz (Hz). A subband, i.e. a sub-band, is a portion of a frequency band, which has the unit of hertz (Hz), and refers to the portion of the spectrum that is located between two specific frequency limits. For a signal, the frequency band is the frequency range between the highest frequency and the lowest frequency contained in the signal.
The invention divides the broadband signal into a plurality of narrow sub-bands with different widths, and a natural directivity function of a central frequency point is taken in each sub-band, and finally the natural directivity function PSF of the central frequency points of all the sub-bands is obtainedO＝{PSF1，PSF2，…，PSFO}。
A beam spectrum acquisition step: performing conventional beam forming calculation of m beams by using array metadata of the uniform line array to obtain power spectrums P of the m beamsM＝{P1，P2，...，Pm}，PmRepresenting the power value of the mth beam.
In the invention, because the number of the array elements of the uniform line array is N, when K (K is more than or equal to 1) mutually irrelevant random incident sound signals exist, the incident signal received by each array element can be expressed as follows:
Figure RE-GDA0003511540100000071
wherein i is the target number, AiFor the ith target signal strength, NoisenIs noise, n is array element serial number;
Figure RE-GDA0003511540100000072
a phase shift vector at each array element for the ith target, wherein p is defined as
Figure RE-GDA0003511540100000073
T denotes transposition, j is an imaginary number, and k is a wave number.
According to conventional beamforming theory, a steering vector S ═ S is defined1，s2，…，sN]TWherein s isn＝(1/N)exp(-jk(n-1)dcosθi) Then power spectrumP (cos θ) can be expressed as
Figure RE-GDA0003511540100000074
Wherein the content of the first and second substances,
Figure RE-GDA0003511540100000081
Figure RE-GDA0003511540100000082
s () is the source distribution function, Bp() As a beam directivity function, K is the number of sound signals, i is the ith sound signal, AiIs the intensity of the ith sound signal, θiIs the incident azimuth angle of the ith sound signal,
Figure RE-GDA0003511540100000083
and d is the array element spacing, and lambda is the wavelength and is the convolution.
R-L iterative calculation step: natural directivity function PSF after subband central frequency point extractionO＝{PSF1，PSF2，…，PSFoPower spectrum P of M wave beamsM＝{P1，P2，…，PmR-L deconvolution calculation is carried out, wherein o represents the total number of divided sub-bands, and o < R.
In the invention, R-L deconvolution is calculated as Richardson-Lucy algorithm, which is an iterative algorithm.
The traditional deconvolution beam forming algorithm needs the deconvolution operation of the beam data of each frequency point of the CBF and the PSF of the corresponding frequency point to obtain the source distribution function SS (f, cos theta) of each frequency point, and finally the source distribution functions of each frequency point are added in phase to obtain the source distribution function SS (f, cos theta)
Figure RE-GDA0003511540100000084
From the above formula, the PSF is located in a narrow frequency bandThe difference is not large, and the higher the frequency is, the wider the frequency band with approximate PSF is, as shown in FIG. 3, therefore the invention greatly reduces the calculation amount of the PSF and deconvolution by dividing the broadband signal into a plurality of narrow sub-bands with unequal widths, and taking the PSF of a central frequency point in each sub-band to perform approximate deconvolution calculation.
A value taking step: and taking m wave beam results obtained by the R-L deconvolution calculation result as a final azimuth estimation result.
In the present invention, the power spectra P of the m beamsMThe natural directivity function PSF after being respectively extracted from the central frequency points of each sub-bandO＝{PSF1，PSF2，…，PSFoPerforming R-L deconvolution calculation to obtain o R-L deconvolution calculation results, and adding all the deconvolution calculation results, wherein the result obtained by the addition is the m wave beam results, namely the final output results of the m wave beams.
In the present invention, the power output due to conventional beamforming can be viewed as the convolution of the point source scattering function of the signal with the natural directivity function of the array, i.e., PM＝PSFM*SMRepresents convolution calculation, and under the condition of known power output and natural directivity, the point source scattering function corresponding to the target azimuth can be recovered by means of deconvolution, and the point source scattering function is used for PMAnd PSF of each frequency pointMAll utilize R-L iterative algorithm to make deconvolution calculation to obtain estimated SMValue of, will each SMAnd taking the m beam results obtained by adding the values as a final azimuth estimation result, wherein the final azimuth estimation result is also an output result formed by deconvolution beams.
Fig. 4 shows a beam output diagram of deconvolution beam forming after the extraction of the molecular band, in which the solid line is the CBF result, the dotted line is the deconvolution result of the present invention, and the dot-dash line is the conventional deconvolution result. The detection capability and the resolution of the method are similar to those of the traditional deconvolution method, and are superior to CBF, and the method has smaller calculated amount.
Example 2:
furthermore, the invention discloses a deconvolution beam forming acceleration method based on an R-L iterative algorithm,
step 1: PSF calculation, namely calculating the natural directivity function PSF of the uniform linear array by utilizing the array type related parametersR＝{PSF1，PSF2，…，PSFrR represents the total number of frequency point numbers.
Step 2: and (4) extracting the molecular bands, namely dividing the sub-bands with unequal intervals according to the bandwidth, and extracting the natural directivity functions of the central frequency points of the corresponding sub-bands.
And step 3: obtaining a beam spectrum, performing conventional beam forming calculation on m beams by using array metadata of the uniform line array, and obtaining a power spectrum P of the m beamsM＝{P1，P2，...，Pm}， PmRepresenting the power value of the mth beam.
The above steps in embodiment 2 are the same as those in embodiment 1, but embodiment 2 of the present invention further includes a period extension step, specifically:
and 4, step 4: power spectra P for m beams respectivelyMAnd a natural directivity function PSF after sub-band central frequency point extractionOPerforming period expansion to obtain an expanded power spectrum PM′And an extended natural directivity function PSFo′(ii) a For a periodically extended power spectrum PM′And a periodically extended natural directivity function PSFo′R-L deconvolution calculation is performed.
In step 4, the power spectrum P for m wave beamsMThe cycle extension comprises the following 2 steps:
by a power spectrum PMThe left and right end points of the beam are starting points, and 0.25m wave beams are respectively expanded outwards;
spreading the beam to power spectrum PMThe values of the left end point and the right end point are taken as reference values to obtain a power spectrum P with the period expansionM′＝{P1，P2，…，Pm，…，P1.5m}，PmRepresenting the power value of the mth beam.
As one embodiment of the power spectrum period expansion, the power spectrum P is usedMThe beams are respectively extended outwards by 0.25m at cos theta of +/-1 as left and right end points, and extended at cos theta of-1The power values of the 0.25m beams are all the power values at cos θ -1, and the power values of the 0.25m beams spread at cos θ -1 are all the power values at cos θ -1.
As another specific embodiment of the power spectrum period expansion, the power spectrum PMThe power values of the 0.25m beams spread at cos θ ═ 1 are symmetric to the curves formed from cos θ ═ 1 to cos θ ═ 0.5, the symmetry axis is cos θ ═ 1, the power values of the 0.25m beams spread at cos θ ═ 1 are symmetric to the curves formed from cos θ ═ 1 to cos θ ═ 0.5, and the symmetry axis is cos θ ═ 1.
R-L iteration of a one-dimensional CBF power spectrum results in blurring at the power spectrum boundary, i.e., around cos θ ± 1 for a uniform linear array. The traditional deconvolution beam forming utilizes the periodicity of a beam directivity pattern to provide a method for expanding the integral range of a power spectrum, and the specific formula is
Figure RE-GDA0003511540100000111
Figure RE-GDA0003511540100000112
By the transformation of the above formula, the boundary ambiguity problem can be moved to be about cos theta ≦ 1.5, so that the boundary ambiguity problem is not influenced in the range of | cos theta | ≦ 1 with practical physical significance, however, the boundary extension mode increases the calculation amount of 50% of beam forming while solving the boundary ambiguity problem of the R-L algorithm, and particularly when the number of beams is large, the application of the algorithm in engineering is greatly influenced by the increase of the calculation amount; the invention optimizes the expansion, and periodically expands M/4 wave beams to two sides by the power spectrum value P (cos theta) near the boundary cos theta of +/-1, thereby omitting the wave beam forming calculation process of M/2 wave beams and improving the speed.
The natural directivity function PSF after the extraction of the central frequency point of the sub-bandOThe cycle extension comprises the following 2 steps:
natural directivity function PSF with subband central frequency pointβ＝{PSF′1，PSF′2，...，PSF′m}βThe left and right endpoints of (2) are starting points, and each beam is extended by 0.25m, β represents the subband number, β ═ 1, 2, 3, …, o }, PSF'mThe natural directivity of the mth wave beam under a certain fixed sub-band central frequency point is shown;
natural directivity function PSF of extended wave beam with sub-band center frequency pointβThe values of the left and right end points are taken as reference values to obtain the expanded natural directivity function PSFβ′＝{PSF′1，PSF′2，...，PSF′m，...，PSF′1.5m}β；
Obtaining a periodically extended natural directivity function PSFO′＝{PSF1′，PSF2′，...，PSFβ′，...，PSFo′}。
Natural directivity function PSFOIs two-dimensional, and has to be directed to the natural directivity function PSF of each sub-band central frequency pointβAnd performing period extension, wherein the method of the period extension is the same as that of the period extension of the power spectrum.
And 5: and taking the middle m beam results in the R-L deconvolution results as final azimuth estimation results.
In step 5 of the invention, the power spectrum P after the period expansion is processedM′Natural directivity function PSF extended from each periodβ′And performing R-L deconvolution calculation to obtain o calculation results, and adding all the calculation results to obtain m intermediate beam results as final azimuth estimation results.
The embodiment 2 of the invention solves the problem of boundary ambiguity caused by a deconvolution R-L algorithm on the premise of not increasing the beam forming calculation amount.
Fig. 5 shows a deconvolution beam output map after the beam spectrum boundary is expanded, the solid line is the CBF result, the dotted line is the deconvolution result of the conventional method for performing periodic expansion on the beam domain, and the dotted line is the deconvolution result of the invention after the CBF energy spectrum boundary is expanded. The method is slightly inconsistent with the traditional deconvolution beam forming at the boundary, and the detection performance and the resolution are consistent with those of the traditional deconvolution method and are superior to CBF.
Fig. 6 shows the comparison result of the calculation speed of the method of the present invention with that of the conventional deconvolution beam forming method, which indicates that the overall calculation speed of the method of the present invention is increased by about 50% at a single frequency point, which is the lower limit of the speed increase of the method of the present invention, and the speed increase of the method of the present invention is increased with the increase of the analysis frequency points.

Claims (9)

1. A deconvolution beam forming acceleration method based on an R-L iterative algorithm is characterized by comprising the following steps:
natural directivity function PSF (pseudo-random function) for calculating uniform linear array by using array type related parametersR＝{PSF1，PSF2，...，PSFrR represents the total number of frequency point numbers;
dividing sub-bands with unequal intervals according to the bandwidth, and extracting natural directivity functions of central frequency points of the corresponding sub-bands;
performing conventional beam forming calculation of m beams by using array metadata of the uniform line array to obtain power spectrums P of the m beamsM＝{P1，P2，...，Pm}；
Natural directivity function PSF after subband central frequency point extractionO＝{PSF1，PSF2，...，PSFoPower spectrum P of M wave beamsM＝{P1，P2，...，PmPerforming R-L deconvolution calculation, wherein o represents the total number of divided sub-bands, and o is less than R;
and taking m wave beam results obtained by the R-L deconvolution calculation result as a final azimuth estimation result.
2. The deconvolution beam forming acceleration method based on R-L iterative algorithm of claim 1, characterized in that the power spectrum P (cos θ) is calculated by
Figure FDA0003299319320000011
Figure FDA0003299319320000012
Figure FDA0003299319320000013
Wherein S () is the source distribution function, Bp() As a beam directivity function, K is the number of sound signals, i is the ith sound signal, AiIs the intensity of the ith sound signal, θiIs the incident azimuth angle of the ith sound signal,
Figure FDA0003299319320000021
and delta () is an impulse function, N is the number of array elements, d is the array element spacing, and lambda is the wavelength.
3. The deconvolution beam forming acceleration method based on R-L iterative algorithm of claim 1, characterized in that the calculation method of the natural directivity function is
Figure FDA0003299319320000022
Wherein f is the signal frequency, N is the number of array elements, d is the array element spacing, c is the sound velocity in water, theta is the signal incidence azimuth angle,
Figure FDA0003299319320000023
is the differential value of the beam angle.
4. The deconvolution beam forming acceleration method based on R-L iterative algorithm as claimed in claim 1, characterized in that the power spectrum P of m beams, the natural directivity function of the center frequency point of the corresponding sub-band is obtainedMThen, power spectra P of m wave beams are respectivelyMAnd a natural directivity function PSF after sub-band central frequency point extractionOPerforming period expansion to obtain an expanded power spectrum PM′And extended natural directivity functionPSFO′(ii) a For a periodically extended power spectrum PM′And a periodically extended natural directivity function PSFO′Performing R-L deconvolution calculation; and taking the middle m beam results in the R-L deconvolution results as final azimuth estimation results.
5. The deconvolution beam forming acceleration method based on R-L iterative algorithm of claim 4, characterized in that, the power spectrum P of m beamsMThe cycle extension comprises the following steps:
by a power spectrum PMThe left end point and the right end point of the beam are used as starting points, and the preset number l of beams is respectively expanded outwards; spreading the beam to power spectrum PMThe values of the left end point and the right end point are taken as reference values to obtain a power spectrum P with the period expansionM′＝{P1，P2，...，Pm，...，Pm+2l}，PmRepresenting the power value of the mth beam.
6. The deconvolution beam forming acceleration method based on R-L iterative algorithm as claimed in claim 5, characterized in that said natural directivity function PSF after extracting the center frequency point of the subbandOThe cycle extension comprises the following steps:
natural directivity function PSF with subband central frequency pointβ＝{PSF′1，PSF′2，...，PSF′m}βThe left and right endpoints of (b) are starting points, a preset number l of beams are respectively extended outwards, beta represents the number of a subband, and beta is {1, 2, 3.mThe natural directivity of the mth wave beam under a certain fixed sub-band central frequency point is shown;
natural directivity function PSF of extended wave beam with sub-band center frequency pointβThe values of the left and right end points are taken as reference values to obtain the expanded natural directivity function PSFβ′＝{PSF′1，PSF′2，...，PSF′m，...，PSF′m+2l}β；
Obtaining a periodically extended natural directivity function PSFO′＝{PSF1′，PSF2′，...，PSFβ′，...，PSFo′}。
7. The deconvolution beam forming acceleration method based on R-L iterative algorithm of claim 6, characterized in that the preset number L is 0.25 m.
8. The deconvolution beam forming acceleration method based on R-L iterative algorithm of claim 6, characterized in that the power of the expanded beam is the power spectrum PMA power value corresponding to the endpoint; the natural directivity of the expanded beam is the natural directivity of the corresponding natural directivity function corresponding to the end point.
9. The deconvolution beam forming acceleration method based on the R-L iterative algorithm of claim 6, characterized in that the power of the beam is spread to power spectrum PMTaking the reference axis of the corresponding end point as a symmetry axis, and carrying out symmetry value taking;
and the natural directivity of the expanded wave beam takes the reference axis of the corresponding endpoint of the corresponding natural directivity function as a symmetry axis to carry out symmetrical value taking.
CN202111190318.7A 2021-10-12 2021-10-12 Deconvolution beam forming acceleration method based on R-L iterative algorithm Pending CN114236516A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111190318.7A CN114236516A (en) 2021-10-12 2021-10-12 Deconvolution beam forming acceleration method based on R-L iterative algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111190318.7A CN114236516A (en) 2021-10-12 2021-10-12 Deconvolution beam forming acceleration method based on R-L iterative algorithm

Publications (1)

Publication Number Publication Date
CN114236516A true CN114236516A (en) 2022-03-25

Family

ID=80743125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111190318.7A Pending CN114236516A (en) 2021-10-12 2021-10-12 Deconvolution beam forming acceleration method based on R-L iterative algorithm

Country Status (1)

Country Link
CN (1) CN114236516A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023148455A1 (en) * 2022-02-02 2023-08-10 Safran Data Systems Device, method and program for recording radiofrequency activity of artificial satellites

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023148455A1 (en) * 2022-02-02 2023-08-10 Safran Data Systems Device, method and program for recording radiofrequency activity of artificial satellites

Similar Documents

Publication Publication Date Title
CN109490850B (en) Broadband array self-adaptive beam forming method under main lobe interference
CN102621527B (en) Broad band coherent source azimuth estimating method based on data reconstruction
CN108872970B (en) Grating lobe discrimination method suitable for general equidistant sparse array single-frequency signal beam forming
CN105334508B (en) A kind of grating lobe suppression method of thinned array broad-band EDFA
CN109597041B (en) Segmented linear frequency modulation waveform design method based on coherent FDA
CN110687207B (en) Sub-wavelength level power-discrimination ultrasonic imaging method based on frequency domain processing
CN114236516A (en) Deconvolution beam forming acceleration method based on R-L iterative algorithm
Yu et al. Clutter suppression for high-resolution atmospheric observations using multiple receivers and multiple frequencies
CN109507643A (en) The method of null-broadening in sidelobe cancellation
CN111538005A (en) SAR front-side-looking imaging method based on FPGA and multiple multi-core DSPs
CN105044680B (en) The phase-coded signal design method of the low Doppler sidelobne of multi-peak
Chen et al. Wave measurements with multi-frequency HF radar in the East China Sea
CN109491009B (en) Optical fiber combined array and grating lobe suppression method based on optical fiber combined array
CN105974409B (en) Satellite-borne sliding spotlight MIMO-SAR imaging method based on multi-frequency sub-band concurrence
CN111175727B (en) Method for estimating orientation of broadband signal based on conditional wave number spectral density
Ma et al. Spatiotemporal two-dimensional deconvolution beam imaging technology
CN109725313B (en) SAR (synthetic aperture radar) sea wave imaging method and system, electronic equipment and medium
CN115236589B (en) Polar region under-ice DOA estimation method based on covariance matrix correction
CN109061597B (en) Ionospheric clutter suppression method based on blind source separation and time-frequency ridge domain filtering
CN114563760B (en) Second-order super-beam forming method, equipment and medium based on SCA array
CN113625222B (en) Broadband signal high-resolution direction estimation method based on time delay difference norm weighting
CN113238194B (en) Broadband phased array radar anti-decoy interference method based on fractional domain-frequency domain processing
CN115902791A (en) S-band wave-measuring radar time Doppler spectrum-based sea wave inversion method and system
CN112924926B (en) DOA estimation method under multi-path channel
CN111257891B (en) Deconvolution-based MIMO sonar distance sidelobe suppression method

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