CN113189543B - Interference suppression method based on motion compensation robust principal component analysis - Google Patents

Interference suppression method based on motion compensation robust principal component analysis Download PDF

Info

Publication number
CN113189543B
CN113189543B CN202110459376.9A CN202110459376A CN113189543B CN 113189543 B CN113189543 B CN 113189543B CN 202110459376 A CN202110459376 A CN 202110459376A CN 113189543 B CN113189543 B CN 113189543B
Authority
CN
China
Prior art keywords
signal
interference
matrix
signals
principal component
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
Application number
CN202110459376.9A
Other languages
Chinese (zh)
Other versions
CN113189543A (en
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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202110459376.9A priority Critical patent/CN113189543B/en
Publication of CN113189543A publication Critical patent/CN113189543A/en
Application granted granted Critical
Publication of CN113189543B publication Critical patent/CN113189543B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
    • G01S5/22Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

An interference suppression method based on motion compensation robust principal component analysis relates to the technical field of underwater acoustic array signals, and aims at the problem that in the prior art, signal detection estimation performance is low due to the influence of interference signals in sea water, and the method comprises the following steps: receiving underwater broadband signals and motion information of a sound source to be detected by using an array, wherein the broadband signals comprise expected signals and interference signals; step two: performing motion compensation on the underwater broadband signal acquired by the array according to the motion information of the sound source to be detected; step three: vectorizing the underwater broadband signal sequence compensated in the second step; step four: and separating the expected signal and the interference signal in the underwater broadband signal after column vectorization by using a robust principal component analysis method. The method has the advantages that the method can separate the expected signal from the interference signal without affecting the signal phase information, thereby inhibiting the influence of interference on the signal detection performance and improving the underwater sound signal detection performance under the condition of low signal-to-interference ratio.

Description

Interference suppression method based on motion compensation robust principal component analysis
Technical Field
The invention relates to the technical field of underwater acoustic array signals, in particular to an interference suppression method based on motion compensation robust principal component analysis.
Background
The ocean is rich in rich animal and plant resources and mineral resources, and has extremely high national defense strategic value, so that the development and exploration of the ocean are increasingly emphasized by various countries. Because of the specificity of the seawater medium, electromagnetic waves cannot be transmitted in water in a long distance, so that sound waves can only be used as a carrier for underwater information transmission, and the underwater acoustic array signal processing is an important foundation for carrying out communication, detection, positioning, navigation and other works underwater by utilizing the sound waves, and is an important research direction of underwater acoustic technology.
In the process of developing and detecting the ocean through the underwater acoustic technology, an underwater acoustic transducer array is usually required to receive acoustic signals, then the acoustic signals are converted into electric signals through an analog-digital converter, and then the electric signals are processed to perform the detection, positioning and other works of the underwater target through solid lines. The ocean is a replicated channel environment and the underwater acoustic transducer array receives a large number of interfering signals, including mechanical noise generated by the surface vessel, ocean environmental noise generated by the marine organisms and the sea water, and electrical noise generated by the signal processing equipment, while the sea receives the desired signal, which can seriously affect the detection estimation performance of the signal.
Disclosure of Invention
The purpose of the invention is that: aiming at the problem that the signal detection estimation performance is low due to the influence of interference signals in seawater in the prior art, an interference suppression method based on motion compensation robust principal component analysis is provided.
The technical scheme adopted by the invention for solving the technical problems is as follows:
an interference suppression method based on motion compensation robust principal component analysis comprises the following steps:
step one: receiving underwater broadband signals and motion information of a sound source to be detected by using an array, wherein the broadband signals comprise expected signals and interference signals;
step two: performing motion compensation on the underwater broadband signal acquired by the array according to the motion information of the sound source to be detected;
step three: vectorizing the underwater broadband signal sequence compensated in the second step;
step four: and separating the expected signal and the interference signal in the underwater broadband signal after column vectorization by using a robust principal component analysis method.
Further, the motion information of the sound source to be measured in the first step specifically includes velocity and acceleration information when the sound source moves.
Further, the motion compensation in the second step specifically includes the following steps:
step two,: calculating the motion distance of the sound source to be measured between two adjacent emission moments by using the speed and acceleration information of the sound source to be measured;
step two: estimating the time delay difference of two adjacent expected signals reaching the matrix according to the motion distance of the sound source to be measured between two adjacent emission moments and the depth of the sound source to be measured;
step two, three: and performing time delay compensation on the underwater broadband signal in the time domain by using the time delay difference.
Further, the depth of the sound source to be measured is obtained by a pressure sensor arranged on the sound source to be measured.
Further, the specific steps of the second step are as follows:
first, let sound source be at t 10 The moment is positioned under the matrix, at t 11 Transmitting acoustic signals at time t 12 Arrive at the matrix at a moment and then at t 21 Transmitting acoustic signals at time t 22 The time delay difference of the two pulse signals reaching the receiving matrix is that
Figure BDA0003041680100000021
Wherein r is 1 And r 2 Respectively at t 11 Time sum t 21 The distance between the target and the matrix, c is the propagation speed of the sound wave under water;
assuming that the target always makes uniform variable-speed linear motion and the depth is unchanged, according to the geometric relationship:
Figure BDA0003041680100000022
wherein v is 0 At t as the target 10 The speed at the moment, a is the movement acceleration of the target, h is the depth of the target,
and obtaining the time delay difference between the received signals at each moment and the starting moment according to the above formula, and obtaining the motion-compensated observation matrix D.
Further, the specific steps of the third step are as follows:
step three: connecting the columns of the matrix with L multiplied by N in one dimension to form a column vector with L multiplied by N multiplied by 1 in one dimension;
step three, two: m times of sampling the column vector obtained in the third step to form an observation matrix with the dimension of (L multiplied by N) multiplied by M;
wherein, N is the number of array elements forming an array, L is the snapshot number of each sampling, and M is the sampling times.
Further, the specific steps of the fourth step are as follows:
firstly, the observation matrix D after motion compensation is regarded as being composed of a low-rank expected signal matrix L and a sparse interference signal matrix S, and the solving process is expressed as follows:
Figure BDA0003041680100000023
wherein I * Represents an infinite norm, i.e., the maximum of the absolute values of the individual elements in the vector, |·|| 1 Representing the 1-norm, i.e. the sum of the absolute values of the individual elements in the vector,
Figure BDA0003041680100000031
and->
Figure BDA0003041680100000032
And respectively estimating values of L and S, solving the optimization problem, and separating the expected signal from the interference signal.
Further, the interference signal in the first step includes ship radiation noise, marine environment noise and circuit noise.
Further, the motion information of the sound source to be measured in the first step is obtained by a motion sensor arranged on the target to be measured and is transmitted through a cable.
Further, the motion sensor includes a compass and an accelerometer.
The beneficial effects of the invention are as follows:
the motion data of the object to be measured in the process of transmitting signals are obtained through the motion sensor, then the array receiving signals are compensated according to the motion data, and finally the expected signals and the interference signals are separated through the robust principal component analysis method. The method has the advantages that the method can separate the expected signal from the interference signal without affecting the signal phase information, thereby inhibiting the influence of interference on the signal detection performance and improving the underwater sound signal detection performance under the condition of low signal-to-interference ratio.
Drawings
FIG. 1 is a schematic diagram of a receive matrix;
FIG. 2 is a primitive number 1 receive signal;
FIG. 3 is a graph of signal contrast before and after robust principal component analysis;
FIG. 4 is a graph showing correlation peaks before and after robust principal component analysis;
fig. 5 is a schematic overall flow chart of the present invention.
Detailed Description
The first embodiment is as follows: referring to fig. 5, a specific description is given of an interference suppression method based on motion compensation robust principal component analysis according to the present embodiment, including the following steps:
step one: receiving underwater broadband signals and motion information of a sound source to be detected by using an array, wherein the broadband signals comprise expected signals and interference signals;
step two: performing motion compensation on the underwater broadband signal acquired by the array according to the motion information of the sound source to be detected;
step three: vectorizing the underwater broadband signal sequence compensated in the second step;
step four: and separating the expected signal and the interference signal in the underwater broadband signal after column vectorization by using a robust principal component analysis method.
The principle of the fourth step is as follows: the broadband signal received by the array in the first step is formed by superposing the expected signal and the interference signal, the signal in the first step is subjected to motion compensation, the compensated signal is subjected to column vectorization, the column vectorized signal is processed by a robust principal component analysis method, and the expected signal and the interference signal are superposed and separated in the received signal.
The array is fixedly arranged on a working ship on the water surface, can transmit or receive underwater acoustic signals, and converts the acoustic signals into electric signals for processing.
A signal source emitting a desired signal is fixedly installed on an underwater object to be measured, and the received signal includes the desired signal and a plurality of interference signals, wherein the desired signal is emitted by a sound source installed on the underwater object to be measured, and the interference signals include mechanical interference generated by a surface ship, marine environmental interference generated by marine organisms and seawater, and electrical interference generated by a signal processing apparatus. Sources of interference signals include ship radiation noise, marine environment noise, circuit noise and the like, and belong to non-Gaussian interference signals.
The motion information of the sound source to be measured is obtained by testing motion sensors such as a compass and an accelerometer which are arranged on the target to be measured, and the motion state of the sound source at the moment when the target emits the expected signal can be reflected. The motion information of the sound source to be measured can be transmitted to the surface operation ship through a cable, and also can be transmitted to the surface operation ship through an acoustic signal.
In the second step, motion compensation is to calculate the motion distance of the target to be measured between two adjacent transmitting moments according to the speed and acceleration information of the target to be measured, then estimate the time delay difference of two adjacent expected signals reaching the matrix by combining the depth of the target to be measured, and then perform time delay compensation on the sampling signals in the time domain, so that the moments of the expected signals reaching the matrix at different sampling moments are the same, and thus the expected signal sequence satisfies the low rank characteristic in the time domain.
The robust principal component analysis method in the fourth step is a data dimension reduction algorithm, and when a matrix is formed by superposing a low-rank matrix and a sparse matrix, the two matrices can be separated through the robust principal component algorithm, and the method has good robustness to non-Gaussian signals.
The target sound source transmits pulse signals at certain time intervals, the water surface matrix receives signals at the same time intervals, each time the signals are received into a sampling period, and data received in each sampling period form a signal matrix with L multiplied by N dimensions.
Each sampling time can obtain a data matrix, n sampling times can obtain n data matrices, to carry out step four, the matrix must be converted into a column vector, n times correspond to n column vectors, then the n column vectors are combined into a new matrix, and the step is to process the new matrix.
The column vectorization is to convert an L×N signal matrix in one sampling period into a column vector with a dimension of (L×N) ×1, and obtain an observation matrix with a dimension of (L×N) ×M through M sampling periods. Since the object to be measured is moving, the propagation times of the desired signals to the matrix are different, resulting in a matrix of desired signals being full-rank.
The motion compensation is performed on observed data, namely, the motion distance of a target to be detected between two adjacent transmitting moments is calculated through the speed and acceleration information of the target to be detected, then the time delay difference of two adjacent expected signals reaching the matrix is estimated by combining the depth of the target to be detected, and then the time delay compensation is performed on sampling signals in a time domain, so that the time when the expected signals reach the matrix at different sampling moments is the same, and thus, the expected signal sequence meets the low rank characteristic in the time domain. Let sound source be t 10 The moment is positioned under the matrix, at t 11 Transmitting acoustic signals at time t 12 The moment reaches the matrix. Then at t 21 Transmitting acoustic signals at time t 22 The time delay difference of the two pulse signals reaching the receiving matrix is that
Figure BDA0003041680100000051
Wherein r is 1 And r 2 Respectively at t 11 Time sum t 21 The distance between the target and the matrix, c, is the propagation velocity of the sound wave under water.
Wherein, the "broadband signal" in the first step is emitted by the sound source.
Assuming that the target always makes uniform variable-speed linear motion and the depth is unchanged, the method can obtain according to the geometric relationship
Figure BDA0003041680100000052
Wherein v is 0 At t as the target 10 The speed at the moment, a is the movement acceleration of the target, and h is the depth of the target.
And obtaining the time delay difference between the received signal at each moment and the initial moment according to the formula, and obtaining the motion-compensated observation matrix D.
The robust principal component analysis method is a data dimension reduction algorithm, and when a matrix is formed by superposing a low-rank matrix and a sparse matrix, the two matrices can be separated through the robust principal component algorithm, and the method has good robustness to non-Gaussian signals. After motion compensation, the observation matrix D may be regarded as being composed of a low-rank desired signal matrix L and a sparse interference signal matrix S, and the solving process may be expressed as the following optimization problem.
Figure BDA0003041680100000053
Solving the optimization problem can accurately separate the expected signal from the interference, thereby achieving the effect of interference suppression.
Examples: the receiving matrix is a uniform linear array composed of 8 array elements, and the array element spacing is 5cm, as shown in figure 1. The sound source to be measured is located 1km below the water surface, the transmitted signal is a linear frequency modulation signal with the bandwidth of 9-15kHz and the pulse width of 15ms, and the interference signal is a complex interference signal acquired from the marine environment. Data is collected for 30ms every sampling period, and the received signal at the first sampling period is shown in fig. 2.
The movement speed of the target at the initial moment is set to be 1m/s, and the acceleration is set to be 0.02m/s 2 When the speed reaches 5m/s, the linear motion starts at a uniform speed, and the sound source emits a pulse signal every 2 s. The time delay difference between each sampling time and the starting time can be calculated according to the formula (2), and is shown in the following table.
Figure BDA0003041680100000054
Figure BDA0003041680100000061
The signals of 10 sampling periods are subjected to motion compensation according to the table, and then an observation matrix is formed through column vectorization. And (3) solving an optimization problem shown in a formula (3), namely separating a desired signal from interference, wherein the separated signal is shown in fig. 3. As can be seen from fig. 3, the desired signal is almost completely separated from the received signal, and noise is greatly suppressed.
To verify whether the robust principal component analysis affects the phase of the desired signal very much, the positions of occurrence of the maximum values of correlation peaks before and after the robust principal component analysis are compared by a cross-correlation method (as shown in fig. 4). It can be seen from fig. 4 that the correlation peaks appear at the same positions and that the side lobes of the correlation peaks are greatly reduced, which also proves that the present invention can effectively suppress noise.
It should be noted that the signals described in the present invention are all acoustic signals, and the acoustic signals may be classified into wideband signals and narrowband signals according to frequencies. The wideband signal is the signal in step one.
It should be noted that the detailed description is merely for explaining and describing the technical solution of the present invention, and the scope of protection of the claims should not be limited thereto. All changes which come within the meaning and range of equivalency of the claims and the specification are to be embraced within their scope.

Claims (8)

1. An interference suppression method based on motion compensation robust principal component analysis is characterized by comprising the following steps:
step one: acquiring an underwater broadband signal by using an array, wherein the broadband signal comprises a desired signal and an interference signal, and then acquiring motion information of a sound source to be detected;
step two: performing motion compensation on underwater sound signals received by the array according to the motion information of the sound source to be detected;
step three: vectorizing the broadband signal column compensated in the second step;
step four: separating the desired signal from the interfering signal using a robust principal component analysis;
the motion compensation in the second step comprises the following specific steps:
step two,: acquiring speed and acceleration information of a target to be detected;
step two: calculating the movement distance of the target to be measured between two adjacent emission moments;
step two, three: estimating the time delay difference of two adjacent expected signals reaching the matrix by combining the depth of the target to be detected;
step two, four: performing time delay compensation on the sampling signal in the time domain;
the specific steps of the second step are as follows:
first, let sound source be at t 10 The moment is positioned under the matrix, at t 11 Transmitting acoustic signals at time t 12 Arrive at the matrix at a moment and then at t 21 Transmitting acoustic signals at time t 22 The time delay difference of the two pulse signals reaching the receiving matrix is that
Figure FDA0004218998420000011
Wherein r is 1 And r 2 Respectively at t 11 Time sum t 21 The distance between the target and the matrix, c is the propagation speed of the sound wave under water;
assuming that the target always makes uniform variable-speed linear motion and the depth is unchanged, according to the geometric relationship:
Figure FDA0004218998420000012
wherein v is 0 At t as the target 10 The speed at the moment, a is the movement acceleration of the target, h is the depth of the target,
and obtaining the time delay difference between the received signals at each moment and the starting moment according to the above formula, and obtaining the motion-compensated observation matrix D.
2. The interference suppression method based on motion compensation robust principal component analysis according to claim 1, wherein the depth of the target is obtained by a pressure sensor provided on the target.
3. The interference suppression method based on motion compensation robust principal component analysis according to claim 2, wherein: the specific steps of the third step are as follows:
step three: connecting the columns of the matrix with L multiplied by N in one dimension to form a column vector with L multiplied by N multiplied by 1 in one dimension;
step three, two: m times of sampling the column vector obtained in the third step to form an observation matrix with the dimension of (L multiplied by N) multiplied by M;
wherein, N is the number of array elements forming an array, L is the snapshot number of each sampling, and M is the sampling times.
4. The interference suppression method based on motion compensation robust principal component analysis according to claim 3, wherein the specific steps of the fourth step are as follows:
firstly, the observation matrix D after motion compensation is regarded as being composed of a low-rank expected signal matrix L and a sparse interference signal matrix S, and the solving process is expressed as follows:
Figure FDA0004218998420000021
wherein, the liquid crystal display device comprises a liquid crystal display device, representation of |· | | | | an infinite norm of the sum of the norms, i.e. the maximum of the absolute values of the individual elements in the vector, I.I 1 Representing the 1-norm, i.e. the sum of the absolute values of the individual elements in the vector,
Figure FDA0004218998420000022
and->
Figure FDA0004218998420000023
And respectively estimating values of L and S, solving the optimization problem, and separating the expected signal from the interference signal.
5. The interference suppression method based on motion compensation robust principal component analysis according to claim 1, wherein: the interference signal in the first step is a non-Gaussian interference signal, and comprises ship radiation noise, ocean environment noise and circuit noise.
6. The interference suppression method based on motion compensation robust principal component analysis according to claim 1, wherein: and in the first step, the motion information of the sound source to be detected is transmitted through a cable.
7. The interference suppression method based on motion compensation robust principal component analysis according to claim 1, wherein: in the first step, the motion information of the sound source to be measured is obtained by a motion sensor arranged on the target to be measured.
8. The method for interference suppression based on motion compensated robust principal component analysis of claim 7, wherein: the motion sensor includes a compass and an accelerometer.
CN202110459376.9A 2021-04-27 2021-04-27 Interference suppression method based on motion compensation robust principal component analysis Active CN113189543B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110459376.9A CN113189543B (en) 2021-04-27 2021-04-27 Interference suppression method based on motion compensation robust principal component analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110459376.9A CN113189543B (en) 2021-04-27 2021-04-27 Interference suppression method based on motion compensation robust principal component analysis

Publications (2)

Publication Number Publication Date
CN113189543A CN113189543A (en) 2021-07-30
CN113189543B true CN113189543B (en) 2023-07-14

Family

ID=76979403

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110459376.9A Active CN113189543B (en) 2021-04-27 2021-04-27 Interference suppression method based on motion compensation robust principal component analysis

Country Status (1)

Country Link
CN (1) CN113189543B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6018317A (en) * 1995-06-02 2000-01-25 Trw Inc. Cochannel signal processing system
JP2005017309A (en) * 2004-09-29 2005-01-20 Francois Bernard System including acoustic array equipped with acoustic velocity meter for measuring position of underwater equipment
CN104101871A (en) * 2013-04-15 2014-10-15 中国科学院声学研究所 Narrowband interference suppression method and narrowband interference suppression system used for passive synthetic aperture
CN107015201A (en) * 2017-06-09 2017-08-04 国家深海基地管理中心 A kind of ultra short baseline locating system is emergent to clock method and system
CN109669160A (en) * 2019-02-21 2019-04-23 哈尔滨工程大学 A kind of detection method of underwater transient state acoustical signal
CN109870695A (en) * 2019-02-21 2019-06-11 哈尔滨工程大学 Noncooperative target multinode hydrolocation method based on deep seafloor reflected sound
CN110703259A (en) * 2019-10-22 2020-01-17 哈尔滨工程大学 Underwater acoustic array channel phase consistency calibration method based on moving sound source
WO2021037252A1 (en) * 2019-08-28 2021-03-04 华南理工大学 Fractional fourier transform echo imaging method and system based on probe compensation

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2415113B (en) * 2004-06-12 2006-10-18 Sonardyne Internat Ltd Robust underwater communication system
US10228449B2 (en) * 2012-03-09 2019-03-12 The United States Of America As Represented By The Secretary Of The Army Method and system for jointly separating noise from signals
CN104181815B (en) * 2014-08-19 2017-02-22 哈尔滨工程大学 Ship movement compensation control method based on environmental estimation
CN106595834B (en) * 2016-11-10 2019-01-04 西北工业大学 A method of obtaining the horizontal longitudinal correlation of the big depth sound field in deep-sea
CN107705295B (en) * 2017-09-14 2021-11-30 西安电子科技大学 Image difference detection method based on robust principal component analysis method
CN108880586B (en) * 2018-06-28 2019-11-08 中国人民解放军战略支援部队信息工程大学 A kind of broadband weak signal enhancement method and apparatus
CN111273238B (en) * 2020-01-06 2022-09-27 中国航天科工集团八五一一研究所 SAR (synthetic aperture radar) wide-band and narrow-band interference simultaneous inhibition method based on low-rank recovery
CN111275669B (en) * 2020-01-13 2022-04-22 西安交通大学 Priori information guided four-dimensional cone beam CT image reconstruction algorithm
CN111239697B (en) * 2020-02-11 2023-02-10 西北工业大学 Multidimensional domain combined SAR broadband interference suppression method based on low-rank matrix decomposition
CN111381212A (en) * 2020-03-04 2020-07-07 哈尔滨工程大学 Virtual ultrashort baseline positioning method based on subarray division
CN111812648B (en) * 2020-07-22 2022-01-04 东南大学 Multichannel synthetic aperture radar RPCA amplitude-phase combined target detection method and device
CN112014791B (en) * 2020-08-28 2024-01-19 安徽合擎智能机器人有限公司 Near-field source positioning method of array PCA-BP algorithm with array error

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6018317A (en) * 1995-06-02 2000-01-25 Trw Inc. Cochannel signal processing system
JP2005017309A (en) * 2004-09-29 2005-01-20 Francois Bernard System including acoustic array equipped with acoustic velocity meter for measuring position of underwater equipment
CN104101871A (en) * 2013-04-15 2014-10-15 中国科学院声学研究所 Narrowband interference suppression method and narrowband interference suppression system used for passive synthetic aperture
CN107015201A (en) * 2017-06-09 2017-08-04 国家深海基地管理中心 A kind of ultra short baseline locating system is emergent to clock method and system
CN109669160A (en) * 2019-02-21 2019-04-23 哈尔滨工程大学 A kind of detection method of underwater transient state acoustical signal
CN109870695A (en) * 2019-02-21 2019-06-11 哈尔滨工程大学 Noncooperative target multinode hydrolocation method based on deep seafloor reflected sound
WO2021037252A1 (en) * 2019-08-28 2021-03-04 华南理工大学 Fractional fourier transform echo imaging method and system based on probe compensation
CN110703259A (en) * 2019-10-22 2020-01-17 哈尔滨工程大学 Underwater acoustic array channel phase consistency calibration method based on moving sound source

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Novel Fast Coherent Detection Algorithm for Radar Maneuvering Target With Jerk Motion;Zhang, Jiancheng等;IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING;第10卷(第5期);1792-1803 *
基于单水听器的移动声源运动参数估计;孟路稳;罗夏云;程广利;尚建华;张明敏;;船舶力学(第11期);全文 *
复杂条件下多通道SAR运动目标检测与参数估计方法研究;张文鹏;中国博士学位论文全文数据库信息科技辑(第01期);I136-274 *
第二讲 合成孔径声纳成像及其研究进展;张春华;刘纪元;;物理(第05期);全文 *

Also Published As

Publication number Publication date
CN113189543A (en) 2021-07-30

Similar Documents

Publication Publication Date Title
CN108226933B (en) Deep sea broadband target depth estimation method based on fringe interference structure
CN111580048B (en) Broadband sound source depth estimation method using single-vector hydrophone
CN103076594B (en) Method for positioning underwater sound pulse signal by double array elements on basis of cross-correlation
CN104678384B (en) Method for estimating underwater target speed by using sound pressure difference cross-correlation spectrum analysis of beam fields
CN107505596B (en) MIMO active detection signal design and detection system and method based on dual-extension underwater acoustic channel environment
CN102645265A (en) Ship radiated noise level measuring method based on virtual time reversal mirror
CN103267966A (en) Orientation estimation method for single two-dimensional differential pressure type vector hydrophone
CN114167423A (en) Radar sea wave parameter measuring method based on depth regression network
Hursky et al. High-frequency (8–16 kHz) model-based source localization
CN113189543B (en) Interference suppression method based on motion compensation robust principal component analysis
CN111679248B (en) Target azimuth and distance combined sparse reconstruction positioning method based on seabed horizontal L-shaped array
CN111458678B (en) Passive distance measurement method based on time-frequency interference spectrum and radiation noise sound intensity measurement
CN115220026B (en) Matching field ice-crossing positioning method based on ice surface horizontal array
US8130592B2 (en) Correlation method for SONAR or RADAR
CN107888372B (en) Underwater sonar communication system based on chaotic oscillator array element
CN103513249A (en) Broadband coherent mold base signal processing method and system
CN114909610A (en) Method for detecting and positioning leakage of underwater oil and gas pipeline and control system
Yu et al. Multiple moving targets detection and parameters estimation in strong reverberation environments
CN113126029A (en) Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment
Talebpour et al. Time Delay Estimation in Underwater Environment using Cross-Correlation Based Techniques
Ostrowski et al. Underwater navigation system based on Doppler shift–measurements and error estimations
CN112684437A (en) Passive distance measurement method based on time domain warping transformation
Magalhães et al. ALMA 2015: Sea trial of an underwater target localization technique using Hausdorff distance
CN115236592B (en) Ice sound positioning method based on single-array-element time-frequency curve matching
CN110208745B (en) Underwater sound positioning method based on self-adaptive matched filter

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