CN108344987A - Multiple submatrixes synthetic aperture sonar frequency domain system function based on numerical computation method - Google Patents
Multiple submatrixes synthetic aperture sonar frequency domain system function based on numerical computation method Download PDFInfo
- Publication number
- CN108344987A CN108344987A CN201810008390.5A CN201810008390A CN108344987A CN 108344987 A CN108344987 A CN 108344987A CN 201810008390 A CN201810008390 A CN 201810008390A CN 108344987 A CN108344987 A CN 108344987A
- Authority
- CN
- China
- Prior art keywords
- computation method
- numerical computation
- system function
- synthetic aperture
- calculated
- 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
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
Abstract
The invention discloses a kind of multiple submatrixes synthetic aperture sonar frequency domain system function based on numerical computation method, includes the following steps：S1, phase dwell point is calculated using numerical computation method；S2, two-dimensional frequency system function is calculated using the phase dwell point that numerical computation method obtains；S3, the orientation pulse pressure dwell point for azimuth match filtering is calculated using numerical computation method；S4, according to orientation pulse pressure dwell point computer azimuth to the phase of matched filtering function；S5, two-dimentional coupled phase is calculated according to two-dimensional frequency system function and azimuth match filter function.The present invention effectively overcomes the Approximation Problem of traditional multiple submatrixes synthetic aperture sonar frequency domain system function, improves the precision and synthetic aperture imaging quality of multiple submatrixes synthetic aperture sonar system function.The present invention not only contributes to the raising of picture quality, while can also replace the performance of other imaging algorithms of time domain imaging algorithms rapid evaluation.
Description
Technical field
The invention belongs to field of signal processing, more particularly to the multiple submatrixes synthetic aperture sonar frequency based on numerical computation method
Domain system function.
Background technology
Multiple submatrixes synthetic aperture sonar imaging is regarded to the inverse process of linear system as, then the structure pair of frequency domain system function
It is most important for fast imaging.Multiple submatrixes synthetic aperture sonar system one way oblique distance course with radical sign containing there are two, because
This is difficult to be utilized algebraic method and calculates parsing and accurate frequency domain system function.Traditional fast imaging algorithm is all built upon accurately
On the approximate basis of system function, this causes what imaging performance was limited to approximation system function pair accurate system function to approach journey
Degree.And high-precision system function is not only advantageous to the raising of picture quality, while also help rapid evaluation other imaging calculate
The performance of method.
Invention content
It is an object of the invention to avoid the Approximation Problem of existing method frequency domain system function, provide one kind can be improved it is more
The multiple submatrixes synthetic aperture sonar frequency domain system based on numerical computation method of submatrix synthetic aperture sonar frequency domain system function precision
System function.
The purpose of the present invention is achieved through the following technical solutions：
A kind of multiple submatrixes synthetic aperture sonar frequency domain system function based on numerical computation method, including following step are provided
Suddenly：
S1, phase dwell point is calculated using numerical computation method；
S2, two-dimensional frequency system function is calculated using the phase dwell point that numerical computation method obtains；
S3, the orientation pulse pressure dwell point for azimuth match filtering is calculated using numerical computation method；
S4, according to orientation pulse pressure dwell point computer azimuth to the phase of matched filtering function；
S5, two-dimentional coupled phase is calculated according to two-dimensional frequency system function and azimuth match filter function.
As a further improvement, in the step S1, phase dwell point is calculated using numerical computation methodIts table
It is up to formula：
Wherein subscript i (i ∈ [1, M]) indicates that the index of reception array element, M indicate to receive the number of array element in receiving array；Indicate distance in two-dimensional space, the specified point mesh that azimuthal coordinates are (r, 0)
Round trip oblique distance course between mark and each sending and receiving array element；Table
Show first-order partial derivative of the round trip oblique distance course about the orientation slow time；diIt indicates to receive between array element and transmitting array element for i-th
Distance；R indicates distance；V indicates sonar bearer rate；T indicates the orientation slow time.Expression is obtained based on numerical computation method
The phase dwell point arrived.fτ、ftIndicate distance to instantaneous frequency and Azimuth Doppler Frequency respectively；fcIndicate carrier frequency；C tables
Show the underwater sound velocity of sound.
As a further improvement, in the step S2, two-dimensional frequency system function Ψ is calculatedi(fτ,ft), expression
Formula is：
As a further improvement, in the step S3, it is calculated using numerical computation method and is filtered for azimuth match
The orientation pulse pressure dwell point of waveSolution, expression formula is：
As a further improvement, in the step S4, phase of the computer azimuth to matched filtering function
Its expression formula is：
As a further improvement, in the step S5, two-dimentional coupled phase is calculatedIts expression formula
For：
Multiple submatrixes synthetic aperture sonar frequency domain system function provided by the invention based on numerical computation method, including it is following
Step：S1, phase dwell point is calculated using numerical computation method；S2, the phase dwell point meter obtained using numerical computation method
Calculate two-dimensional frequency system function；S3, the orientation pulse pressure dwell point for azimuth match filtering is calculated using numerical computation method；
S4, according to orientation pulse pressure dwell point computer azimuth to the phase of matched filtering function；S5, according to two-dimensional frequency system function and side
Position calculates two-dimentional coupled phase to matched filtering function.The present invention effectively overcomes traditional multiple submatrixes synthetic aperture sonar frequency domain system
The Approximation Problem of system function, improves the precision and image quality of multiple submatrixes synthetic aperture sonar system function.The present invention is not only
Be conducive to the raising of picture quality, while the performance of other imaging algorithms of time domain imaging algorithms rapid evaluation can also be replaced.
Description of the drawings
Using attached drawing, the invention will be further described, but the embodiment in attached drawing does not constitute any limit to the present invention
System, for those of ordinary skill in the art, without creative efforts, can also obtain according to the following drawings
Other attached drawings.
Fig. 1 is the calculation process signal of the multiple submatrixes synthetic aperture sonar frequency domain system function based on numerical computation method
Figure.
Fig. 2 is the two-dimensional imaging geometry of multiple submatrixes synthetic aperture sonar.
Fig. 3 is the imaging results based on numerical computations frequency domain system function.
Specific implementation mode
It is below in conjunction with the accompanying drawings and specific real in order to make those skilled in the art more fully understand technical scheme of the present invention
Applying example, the present invention is described in further detail, it should be noted that in the absence of conflict, embodiments herein and
Feature in embodiment can be combined with each other.
As shown in Figure 1, the multiple submatrixes synthetic aperture sonar frequency domain system letter based on numerical computation method described in the present embodiment
Several implementing procedures mainly consist of the following steps：Phase dwell point is calculated using numerical computation method；Utilize numerical computations side
The phase dwell point that method obtains calculates two-dimensional frequency system function；It is calculated using numerical computation method and is filtered for azimuth match
Orientation pulse pressure dwell point；According to orientation pulse pressure dwell point computer azimuth to the phase of matched filtering function；According to two-dimensional frequency
It is that function and azimuth match filter function calculate two-dimentional coupled phase.The present invention effectively overcomes traditional multiple submatrixes synthetic aperture
The Approximation Problem of sonar frequency domain system function improves the precision and image quality of multiple submatrixes synthetic aperture sonar system function.
The present invention not only contributes to the raising of picture quality, while can also replace other imaging algorithms of time domain imaging algorithms rapid evaluation
Performance.
Fig. 2 show multiple submatrixes synthetic aperture sonar system two-dimensional imaging geometry, without loss of generality, it is assumed that in two-dimensional space
There are an ideal point target, azimuthal coordinates 0 are r apart from coordinate.Platform with speed v during at the uniform velocity being advanced, hair
It penetrates array element and emits the broadband signal unrelated with position to positive side-looking direction with fixed pulse recurrence frequency simultaneously.According to two dimension
Imaging geometry, when by time t, the azimuthal coordinates for emitting array element are vt, receive array element, transmitting array element i-th at this time
With the round trip oblique distance course R of point targeti(t；R) it is：
Wherein subscript i indicates to receive the index of array element；R indicates distance；T indicates the orientation slow time；diIt indicates to connect for i-th
Receive the distance between array element and transmitting array element；C indicates the underwater sound velocity of sound.
The derivative about the slow time t of orientation of round trip oblique distance course is：
According to numerical computation method, phase dwell point can be obtainedIts calculation formula is：
Wherein fτ、ftIndicate distance to instantaneous frequency and Azimuth Doppler Frequency respectively；fcIndicate carrier frequency.
Based on the phase dwell point that numerical computations obtain, two-dimensional frequency system function can be expressed as：
Based on numerical computation method, the orientation pulse pressure dwell point filtered for azimuth match can be calculatedIts
Calculation formula is：
According to orientation pulse pressure dwell point, the phase function for being available for azimuth match filtering is：
According to two-dimensional frequency system function and orientation matched filtering function, the two dimension coupling between azimuth-range can be obtained
Phase, i.e.,：
The frequency domain system function that multiple submatrixes synthetic aperture sonar can be obtained after being handled according to above-mentioned steps, with distance to point
The imaging algorithm of block, and the frequency domain system function based on numerical computations can carry out target imaging, and imaging results are as schemed
Shown in 3.
The above described is only a preferred embodiment of the present invention, be not intended to limit the present invention in any form, though
So the present invention has been disclosed as a preferred embodiment, and however, it is not intended to limit the invention, any technology people for being familiar with this profession
Member, without departing from the scope of the present invention, when the technology contents using the disclosure above make a little change or modification
For the equivalent embodiment of equivalent variations, as long as being the content without departing from technical solution of the present invention, according to the technical essence of the invention
To any simple modification, equivalent change and modification made by above example, in the range of still falling within technical solution of the present invention.
Claims (6)
1. the multiple submatrixes synthetic aperture sonar frequency domain system function based on numerical computation method, which is characterized in that including following step
Suddenly：
S1, phase dwell point is calculated using numerical computation method；
S2, two-dimensional frequency system function is calculated using the phase dwell point that numerical computation method obtains；
S3, the orientation pulse pressure dwell point for azimuth match filtering is calculated using numerical computation method；
S4, according to orientation pulse pressure dwell point computer azimuth to the phase of matched filtering function；
S5, two-dimentional coupled phase is calculated according to two-dimensional frequency system function and azimuth match filter function.
2. the multiple submatrixes synthetic aperture sonar frequency domain system function according to claim 1 based on numerical computation method,
It is characterized in that：In the step S1, phase dwell point is calculated using numerical computation methodIts expression formula is：
Wherein subscript i (i ∈ [1, M]) indicates that the index of reception array element, M indicate to receive the number of array element in receiving array；Indicate distance in two-dimensional space, the specified point that azimuthal coordinates are (r, 0)
Round trip oblique distance course between target and each sending and receiving array element；
Indicate first-order partial derivative of the round trip oblique distance course about the orientation slow time；diIndicate i-th receive array element and transmitting array element it
Between distance；R indicates distance；V indicates sonar bearer rate；T indicates the orientation slow time；It indicates to be based on numerical computation method
Obtained phase dwell point；fτ、ftIndicate distance to instantaneous frequency and Azimuth Doppler Frequency respectively；fcIndicate carrier frequency；c
Indicate the underwater sound velocity of sound.
3. the multiple submatrixes synthetic aperture sonar frequency domain system function according to claim 2 based on numerical computation method,
It is characterized in that：In the step S2, two-dimensional frequency system function Ψ is calculatedi(fτ,ft), expression formula is：
4. the multiple submatrixes synthetic aperture sonar frequency domain system function according to claim 3 based on numerical computation method,
It is characterized in that：In the step S3, it is calculated using numerical computation method resident for the orientation pulse pressure of azimuth match filtering
PointSolution, expression formula is：
5. the multiple submatrixes synthetic aperture sonar frequency domain system function according to claim 4 based on numerical computation method,
It is characterized in that：In the step S4, phase of the computer azimuth to matched filtering functionIts expression formula is：
6. the multiple submatrixes synthetic aperture sonar frequency domain system function according to claim 5 based on numerical computation method,
It is characterized in that：In the step S5, two-dimentional coupled phase is calculatedIts expression formula is：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810008390.5A CN108344987B (en) | 2018-01-04 | 2018-01-04 | Numerical value calculation-based method for multi-subarray synthetic aperture sonar frequency domain function |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810008390.5A CN108344987B (en) | 2018-01-04 | 2018-01-04 | Numerical value calculation-based method for multi-subarray synthetic aperture sonar frequency domain function |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108344987A true CN108344987A (en) | 2018-07-31 |
| CN108344987B CN108344987B (en) | 2022-06-14 |
Family
ID=62960371
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810008390.5A Active CN108344987B (en) | 2018-01-04 | 2018-01-04 | Numerical value calculation-based method for multi-subarray synthetic aperture sonar frequency domain function |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108344987B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109116364A (en) * | 2018-08-03 | 2019-01-01 | 中国人民解放军91388部队 | Based on the equivalent single gust of conversion method of the tactic multiple submatrixes SAS echo-signal of time domain |
| CN113640807A (en) * | 2021-06-23 | 2021-11-12 | 中国人民解放军海军工程大学 | Multi-subarray synthetic aperture sonar intra-pulse Doppler frequency shift compensation line-by-line imaging method |
| CN116500626A (en) * | 2023-06-30 | 2023-07-28 | 海底鹰深海科技股份有限公司 | Multi-receiving array element data receiving-transmitting combination conversion method |
| CN116990791A (en) * | 2023-09-25 | 2023-11-03 | 海底鹰深海科技股份有限公司 | Multi-receiving-array synthetic aperture sonar echo simulation method |
| CN113640807B (en) * | 2021-06-23 | 2024-04-30 | 中国人民解放军海军工程大学 | Multi-subarray synthetic aperture sonar intra-pulse Doppler frequency shift compensation progressive imaging method |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1227681A (en) * | 1996-08-22 | 1999-09-01 | 金桥技术有限公司 | Symbol-matched filter having a low silicon and power requirement |
| US8437431B1 (en) * | 2007-09-20 | 2013-05-07 | Gregory Hubert Piesinger | Sequential decoder fast incorrect path elimination method and apparatus for pseudo-orthogonal coding |
| CN203965527U (en) * | 2014-07-11 | 2014-11-26 | 沈阳华兴防爆器材有限公司 | A kind of electrostatic admeasuring apparatus |
| WO2015126505A2 (en) * | 2013-12-06 | 2015-08-27 | Lynch Jonathan J | Methods and apparatus for reducing noise in a coded aperture radar |
| CN105676220A (en) * | 2014-11-21 | 2016-06-15 | 中国航空工业集团公司雷华电子技术研究所 | High-resolution bunching SAR auto-focusing method based on two-dimensional inverse filtering |
| CN106794009A (en) * | 2014-10-07 | 2017-05-31 | 蝴蝶网络有限公司 | Ultrasonic signal processing circuit and relevant device and method |
| CN106772326A (en) * | 2016-12-05 | 2017-05-31 | 中国人民解放军91388部队 | A kind of multiple submatrixes synthetic aperture sonar phase error analysis method |
| CN106842210A (en) * | 2016-12-05 | 2017-06-13 | 中国人民解放军91388部队 | A kind of new multiple submatrixes synthetic aperture sonar fast imaging algorithm |
| CN107219512A (en) * | 2017-03-29 | 2017-09-29 | 北京大学 | A kind of sound localization method based on acoustic transfer function |
-
2018
- 2018-01-04 CN CN201810008390.5A patent/CN108344987B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1227681A (en) * | 1996-08-22 | 1999-09-01 | 金桥技术有限公司 | Symbol-matched filter having a low silicon and power requirement |
| US8437431B1 (en) * | 2007-09-20 | 2013-05-07 | Gregory Hubert Piesinger | Sequential decoder fast incorrect path elimination method and apparatus for pseudo-orthogonal coding |
| WO2015126505A2 (en) * | 2013-12-06 | 2015-08-27 | Lynch Jonathan J | Methods and apparatus for reducing noise in a coded aperture radar |
| CN203965527U (en) * | 2014-07-11 | 2014-11-26 | 沈阳华兴防爆器材有限公司 | A kind of electrostatic admeasuring apparatus |
| CN106794009A (en) * | 2014-10-07 | 2017-05-31 | 蝴蝶网络有限公司 | Ultrasonic signal processing circuit and relevant device and method |
| CN105676220A (en) * | 2014-11-21 | 2016-06-15 | 中国航空工业集团公司雷华电子技术研究所 | High-resolution bunching SAR auto-focusing method based on two-dimensional inverse filtering |
| CN106772326A (en) * | 2016-12-05 | 2017-05-31 | 中国人民解放军91388部队 | A kind of multiple submatrixes synthetic aperture sonar phase error analysis method |
| CN106842210A (en) * | 2016-12-05 | 2017-06-13 | 中国人民解放军91388部队 | A kind of new multiple submatrixes synthetic aperture sonar fast imaging algorithm |
| CN107219512A (en) * | 2017-03-29 | 2017-09-29 | 北京大学 | A kind of sound localization method based on acoustic transfer function |
Non-Patent Citations (4)
| Title |
|---|
| STANDARD: ""Air Interface for Fixed Broadband"", 《IEEE》 * |
| 杨博 等: ""基元间隔不等线列阵指向性仿真研究"", 《声学技术》 * |
| 殷钊 等: ""一种多子阵合成孔径声纳成像算法研究"", 《舰船电子工程》 * |
| 王阳阳 等: ""弹载前斜视成像技术研究"", 《工程科技Ⅱ辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109116364A (en) * | 2018-08-03 | 2019-01-01 | 中国人民解放军91388部队 | Based on the equivalent single gust of conversion method of the tactic multiple submatrixes SAS echo-signal of time domain |
| CN109116364B (en) * | 2018-08-03 | 2023-01-03 | 中国人民解放军91388部队 | Multi-subarray SAS echo signal equivalent single-array conversion method based on time domain sequential arrangement |
| CN113640807A (en) * | 2021-06-23 | 2021-11-12 | 中国人民解放军海军工程大学 | Multi-subarray synthetic aperture sonar intra-pulse Doppler frequency shift compensation line-by-line imaging method |
| CN113640807B (en) * | 2021-06-23 | 2024-04-30 | 中国人民解放军海军工程大学 | Multi-subarray synthetic aperture sonar intra-pulse Doppler frequency shift compensation progressive imaging method |
| CN116500626A (en) * | 2023-06-30 | 2023-07-28 | 海底鹰深海科技股份有限公司 | Multi-receiving array element data receiving-transmitting combination conversion method |
| CN116500626B (en) * | 2023-06-30 | 2023-09-19 | 海底鹰深海科技股份有限公司 | Multi-receiving array element data receiving-transmitting combination conversion method |
| CN116990791A (en) * | 2023-09-25 | 2023-11-03 | 海底鹰深海科技股份有限公司 | Multi-receiving-array synthetic aperture sonar echo simulation method |
| CN116990791B (en) * | 2023-09-25 | 2023-11-28 | 海底鹰深海科技股份有限公司 | Multi-receiving-array synthetic aperture sonar echo simulation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108344987B (en) | 2022-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Reis et al. | Design and experimental validation of a USBL underwater acoustic positioning system | |
| CN108344987A (en) | Multiple submatrixes synthetic aperture sonar frequency domain system function based on numerical computation method | |
| CN107942284B (en) | Underwater direction of arrival estimation method and device based on two-dimensional orthogonal non-uniform linear array | |
| CN106324597B (en) | The translational compensation and imaging method of big corner ISAR radar based on PFA | |
| Chitre et al. | Ambient noise imaging in warm shallow waters; robust statistical algorithms and range estimation | |
| WO2020072041A1 (en) | Systems and methods for stereo radar tracking | |
| CN106842210A (en) | A kind of new multiple submatrixes synthetic aperture sonar fast imaging algorithm | |
| CN108535682B (en) | Underwater two-dimensional DOA estimation method and device based on rotation non-uniform double-L array | |
| CN106896360A (en) | A kind of FPGA implementation method of SAR signal processing algorithms | |
| CN108445444A (en) | Correct the direct localization method of volume Kalman filtering | |
| CN108519576B (en) | Underwater direction of arrival estimation method and device based on included angle adjustable non-uniform linear array | |
| CN104569917A (en) | Self-adaption positioning and navigation system and method for AUV (automatic underwater vehicle) platform | |
| CN109581275B (en) | Two-dimensional underwater DOA estimation method and device based on non-circular signal and three-dimensional orthogonal array | |
| CN104407340A (en) | Device and method for calibrating lineup of dragging linear array | |
| CN109581388B (en) | Near-field wide-view-angle beam forming method of real-time three-dimensional imaging sonar | |
| CN109444901A (en) | A kind of multiple submatrixes SAS sub-aperture image method under isomerous environment | |
| Dubrovinskaya et al. | Underwater localization via wideband direction-of-arrival estimation using acoustic arrays of arbitrary shape | |
| WO2017063032A1 (en) | Determining elevation and bearing information of a remote point | |
| CN106569180B (en) | Prony method-based orientation estimation algorithm | |
| CN103513238A (en) | A regularization least square subspace crossing target direction finding method | |
| CN103487796B (en) | A kind of method utilizing underwater acoustic channel Statistically invariant feature to realize passive ranging | |
| CN107092009A (en) | A kind of indoor orientation method and device | |
| CN112710984B (en) | Passive positioning method and system based on homotopy continuation | |
| CN109541573A (en) | A kind of element position calibration method being bent hydrophone array | |
| AU2015324511A1 (en) | Phase center alignment for fixed repetition rate synthetic aperture systems |
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 |