CN112748376B - Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method - Google Patents

Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method Download PDF

Info

Publication number
CN112748376B
CN112748376B CN202011601964.3A CN202011601964A CN112748376B CN 112748376 B CN112748376 B CN 112748376B CN 202011601964 A CN202011601964 A CN 202011601964A CN 112748376 B CN112748376 B CN 112748376B
Authority
CN
China
Prior art keywords
positioning
underwater
matrix
magnetic
target
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
CN202011601964.3A
Other languages
Chinese (zh)
Other versions
CN112748376A (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 CN202011601964.3A priority Critical patent/CN112748376B/en
Publication of CN112748376A publication Critical patent/CN112748376A/en
Application granted granted Critical
Publication of CN112748376B publication Critical patent/CN112748376B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/0094Sensor arrays
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/02Measuring direction or magnitude of magnetic fields or magnetic flux
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices

Abstract

The invention discloses an underwater three-dimensional space passive magnetic field positioning method based on an MP-MUSIC method, and belongs to the field of underwater passive magnetic field positioning. The underwater magnetic field source target positioning method mainly positions the underwater magnetic field source target, can overcome underwater interference, and ensures that accurate positioning can be realized in a complex marine environment. The positioning method adopts the magnetic sensor receiving array, has the advantages of small volume, portability and arrangement, reduces the calculation complexity by optimizing the MP-MUSIC algorithm, realizes the quick optimization approach of the positioning algorithm path, realizes the accurate positioning in the complex large-scene marine environment by two processes of fuzzy positioning and refined positioning, and has the positioning accuracy reaching the meter order of magnitude in a detection area of one hundred meters. The technology can be applied to positioning of the leakage position of the submarine power supply cable, detection and tracking of the moving position of an underwater vehicle, the position of a frogman and the like, and has a great application value.

Description

Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method
Technical Field
The invention relates to an underwater three-dimensional space passive magnetic field positioning method, in particular to an underwater three-dimensional space passive magnetic field positioning method based on an MP-MUSIC method, and belongs to the field of underwater magnetic field positioning.
Background
The traditional underwater detection positioning technology is mainly sonar detection, but the underwater sound propagation environment in a shallow sea water area is extremely deteriorated and is strongly interfered by echo reverberation of sonar, so that effective identification and tracking of a target cannot be realized. The underwater electromagnetic field is not easily influenced by environmental factors, and the target detection in shallow sea water areas is extremely advantageous. Underwater detection and communication based on electromagnetic field are used as an effective supplement of sound detection, and research surge is raised in scientific research institutions of various countries in recent years. The underwater target electromagnetic field positioning technology is a technology that a magnetic sensor equipped on a positioning device is used as an electromagnetic field signal receiving array, and an excitation magnetic field signal generated by a positioning target is measured and processed to estimate a target space position parameter.
In the last 70 th century, kuipers, raab proposed locating objects using electromagnetic tracking, and Kuipers et al used this magnetic field to locate objects by controlling the current in electrical coils to create a swinging magnetic field around the pointer vector. The northwest university of 2012 tenses, and li bin adopts an electromagnetic vector sensor to estimate DOA and polarization parameters of a plurality of signal sources in seawater, and based on a TLS-ESPRIT algorithm, successfully distinguishes 5 uncorrelated monochromatic underwater electromagnetic field sources. In 2014, the method for positioning the magnetic dipole by the Gong balance of the reworking group 710 of the Chinese ship by utilizing the magnetic field vector and the magnetic field tensor data improves the original single-point magnetic tensor, and the new magnetic tensor measurement system has higher underwater target positioning accuracy. In 2018, xulongfei, zheng Qusong and the like of the university of northwest industry research the problem of positioning underwater non-transparent objects by an electromagnetic field by adopting a total field-scattering field source adding technology in a time domain finite difference method of a shift operator, so that the calculation efficiency is improved, and the calculation time is shortened.
The underwater target electromagnetic field positioning technology utilizes measurement electromagnetic field signals to realize target positioning, and the equipment sensors only need a certain number of magnetic sensors, are easier to carry and set, and have lower cost compared with optical detection equipment and underwater acoustic detection equipment. The underwater electromagnetic field positioning technology has a good application scene in the detection and tracking of the shallow sea underwater target, so that the innovative exploration and the deep research on the related technical problems of the underwater electromagnetic field positioning still have important scientific research value and practical technical application significance.
Disclosure of Invention
The invention provides an underwater three-dimensional space passive magnetic field positioning method based on an MP-MUSIC optimization algorithm, which has strong magnetic resistance and interference resistance, can overcome underwater interference, and can accurately position a magnetic dipole source target in an underwater environment through algorithm inversion.
The purpose of the invention is realized by the following steps:
an underwater three-dimensional space passive magnetic field positioning method based on an MP-MUSIC method is characterized by comprising the following steps:
(1) In the length of L meterIn the underwater environment within the range of the positioning area with the width of W meters and the height of H meters, the positioning target is n magnetic dipole sources, and the magnetic dipole moment direction is
Figure GDA0003876741690000021
The magnetic dipole source is positioned
Figure GDA0003876741690000022
Arranging a circular magnetic sensor array underwater, wherein the number of the sensors is k, and the positions of sensor array elements are
Figure GDA0003876741690000023
I is more than or equal to 1 and less than or equal to k, the distance vector from the magnetic dipole source to the sensor array element is
Figure GDA0003876741690000024
(2) The magnetic induction intensity of the magnetic dipole source in the three-dimensional space measured by each magnetic sensor is as follows:
Figure GDA0003876741690000025
μ r for relative permeability, mu 0 Permeability of 4 pi x 10 in vacuum -7 H/m, N' is the number of turns of the coil in the closed current loop, I is the closed current for generating the magnetic dipole source, R is the radius of the coil, N i For underwater interference, T represents matrix transposition, and gain matrix is defined as a i (r i ) Let B i =a i (r i )m+n i The magnetic induction including three directions of intensity
Figure GDA0003876741690000026
Wherein
Figure GDA0003876741690000027
(3) Storing the magnetic induction intensity data measured by the k magnetic sensor array into an array
Figure GDA0003876741690000028
Where B is i Is the ith array element testThe magnetic induction of the quantity, D is a 3k × 1 matrix, and the array constituted by the gain matrix is A = [ a ] 1 T (r 1 ),a 2 T (r 2 ),…,a k T (r k )] T Where a is i (r i ) Is the gain matrix corresponding to the ith array element, A is 3k × 3 matrix, the corresponding matrix form operation relationship is D = Am + N,
Figure GDA0003876741690000029
is an underwater disturbance;
(4) Solving the covariance matrix R DD =E{DD H Where R is DD Is a 3k x 3k order matrix, H represents the conjugate transpose, E {. Cndot. } represents the desired operation, and R is obtained by matrix transformation DD =AE{mm H }A H2 I, where σ 2 I is a 3k multiplied by 3k order unit matrix for interference noise power;
(5) For covariance matrix R DD Decomposing the characteristic value to obtain R DD =U∑U H Where U is a 3k × 3k order eigenvector matrix and Σ is a diagonal matrix λ consisting of eigenvalues 1 ≥λ 2 ≥…λ 3k ≧ 0, defined by λ 12 ,…λ n Form a one-dimensional matrix sigma S From λ n+1n+2 ,…,λ 3k Form a one-dimensional matrix sigma N N magnetic dipole sources;
(6)∑ S the characteristic vector group corresponding to the medium characteristic value is defined as U S ,∑ N The feature vector group corresponding to the medium feature value is defined as U N Then U = [ U ] may be represented S ,U N ]Wherein U is S Is a signal subspace, U N Is a noise subspace, an orthogonal projection matrix P of the signal subspace is obtained 3k×3k =I-U S U S H
(7) Calculating the target function of the MP-MUSIC algorithm as the minimum lambda of the generalized characteristic decomposition min (A H P 3k×3k A,A H A);
(8) The process of traversing and solving the minimum characteristic value by the optimized MP-MUSIC algorithm is divided into a fuzzy traversal process and a refining traversal process, the traversal process traverses all positions in a positioning area and searches for the globally minimum generalized characteristic value, and the method comprises the following steps:
(8-1) fuzzy traversal process: dividing the target location area into N × M × K volumes
Figure GDA0003876741690000031
The center point position of the region block
Figure GDA0003876741690000032
Representing the region, j is more than or equal to 1 and less than or equal to NxMxK, r tj Including position information of fuzzy traversal process, distance vector from magnetic dipole source to sensor array element
Figure GDA0003876741690000033
Traverse r t1 ,r t2 ,…,r tN×M×K Calculating a gain matrix A, substituting the gain matrix A into the generalized minimum eigenvalue lambda for solving different position points min (A H P 3k×3k A,A H A) The central point corresponding to the minimum characteristic value is the target position obtained by fuzzy traversal, and the area block where the position is located is the area block where the target is likely to appear in the whole situation;
(8-2) refining the traversal process: the area block which is possibly generated by the target and is solved by the fuzzy traversal process is taken as a new positioning area of the target, and the area is subdivided into NN multiplied by MM multiplied by KK volumes
Figure GDA0003876741690000034
Traversing the position of the central point of the area block according to the step (8-1) to find the minimum characteristic value;
(8-3) when the refinement traversal process is finished, the minimum characteristic value is the global minimum generalized characteristic value;
(9) And the characteristic value corresponding point is the position of the target magnetic dipole source estimated by the MP-MUSIC algorithm, the algorithm is finished, and the positioning is completed.
The invention also includes such features:
in the step (1), k magnetic sensor arrays;
in step (7), the objective function is λ min (A H P 3k×3k A,A H A) Therefore, the position of the target source can be located by finding the minimum characteristic value;
and (8) positioning through a fuzzy traversal process and a refining traversal process.
Compared with the prior art, the invention has the beneficial effects that:
the receiving array of the multi-polarization circular magnetic sensor is adopted, so that the sensitivity is high, the magnetic field signal is accurately received, the effective receiving and estimation of the magnetic field signal in any direction polarization mode of a radiation source are ensured, the volume of the magnetic sensor array is small, the power consumption is low, the chip packaging is basically realized, the carrying is convenient, and the setting is easy; the optimized MP-MUSIC algorithm has clear steps, stable performance, small calculation complexity and calculation amount, higher calculation speed than that of a common algorithm, can realize quick positioning, overcomes underwater interference and ensures high precision. The technology can be applied to positioning of the leakage position of the submarine cable, the moving position of an underwater vehicle, detection of the position of a frogman and the like, and has a great application value.
Drawings
FIG. 1 is a schematic view of a positioning device within a positioning area;
FIG. 2 is a schematic diagram of a receiving array arrangement for k magnetic induction sensors;
FIG. 3 is an overall flowchart of the method for positioning the passive magnetic field in the underwater three-dimensional space based on the MP-MUSIC algorithm;
FIG. 4 is a flow chart for fuzzy traversal;
FIG. 5 is a flow diagram for a refinement traversal.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
Example 1: an underwater detection region shown in the figure 1 is constructed by using simulation software, and a three-dimensional space coordinate axis is established, wherein the length unit is m. At coordinates (7, 0, 2) a magnetic dipole source is arranged, with the direction vector of the magnetic dipole moment set to (0, 1). The receiving magnetic sensors are arranged according to a k magnetic sensor receiving array in fig. 2, the radius of the receiving array is taken as 20m, and the height difference of each annular array is taken as 9m. And (3) acquiring the magnetic induction intensity data of the receiving sensor in a simulation manner, and adding Gaussian white noise to simulate underwater interference. Supposing that the position of a target magnetic dipole source is unknown, scrambling data is obtained through a magnetic sensor receiving array, the data is substituted into an optimized MP-MUSIC positioning algorithm program to obtain an optimal estimated position, the optimal estimated position is compared with a target source position (7, 0, 2) set by simulation, and finally whether the positioning algorithm can meet the index requirement is judged. The method comprises the following specific steps:
s100: in a relative magnetic permeability of
Figure GDA0003876741690000041
The seawater in which a magnetic dipole source to be positioned is arranged consists of a closed electrified coil, the number of turns of the coil is 1, the current of a closed circuit is 1A, and the radius of the coil is 1m.
S200:24 magnetic sensor array as shown in fig. 2, each of the magnetic sensor arrays has 8 magnetic sensors, the radius of the receiving array is 20m, the height difference of each of the annular arrays is 9m, and the coordinate positions of the receiving magnetic sensors 1,2, \ 8230;, 24 are set to (0, 20, 19), (14.1421, 19), (20, 0, 19), (14.1421, -14.1421, 19), (0, -20, 19), (-14.1421, -14.1421, 19), (0, 20, 19), (-14.1421, 19), (0, 20, 10), (14.1421, 10), (20, 0, 10), (14.1421, -14.1421, 10), (0, -20, 10), (-14.1421, 10), (-20, 0, 10), (-14.1421, 10), (0, 20, 1), (14.1421, 1), (20, 0, 1), (14.1421, -14.142, 1), (20, 0, 1), (-14.1421, 1), (-20, 0, 1), (-14.1421, 1). The magnetic sensor array measures magnetic induction intensity, so the array measurement data is written as D = [ B = 1 T ,B 2 T ,…,B k T ] T ,B i The magnetic induction measured for the magnetic sensor array element i. In order to simulate underwater interference, gaussian white noise n is added into the measured magnetic induction intensity data i ,B i And n i The signal-to-noise ratio of (A) is 20dB, so the measured data is B i =B i +n i The magnetic induction in the x direction in this example was measured to be [ -3.767 × 10 [ ] -5 ,3.328×10 -5 ,8.219×10 -5 ,3.032×10 -5 ,1.482×10 -5 ,-1.734×10 -5 ,-6.580×10 -5 ,-2.584×10 -5 ,1.366×10 -5 ,2.928×10 -5 ,1.147×10 -4 ,4.208×10 -5 ,2.009×10 -6 ,-2.928×10 -5 ,-9.329×10 -5 ,9.495×10 -6 ,1.185×10 -4 ,-3.741×10 -5 ,-7.039×10 -5 ,-5.529×10 -5 ,-3.327×10 -5 ,7.069×10 -5 ,-7.591×10 -6 ,-4.295×10 -6 ]。
S310: the positioning process is shown in FIG. 3, and the existing magnetic induction data D is measured w And substituting the data into an MP-MUSIC optimization algorithm. Setting the position of a magnetic dipole source to be positioned as r t =[x t ,y t ,z t ]Magnetic dipole moment direction of m = [0, 1 =] T
S320: for measured data D w To calculate the covariance matrix R DwDw Decomposition of eigenvalues to yield R DwDw =U∑U H Splitting the eigenvector component into signal subspaces U S Sum noise subspace U N Determining the signal subspace orthographic projection P k×k ⊥=I-U S U S H And I is a k multiplied by k order identity matrix.
S330: solving for the minimum eigenvalue λ at the updated estimated position min (A w H P k×k A w ,A w H A w ) In order to accurately position the target, the operation of updating the estimated position solution characteristic value is divided into a fuzzy traversal process and a detailed traversal process, and the specific operations are as follows:
s331: the fuzzy traversal process is as shown in fig. 4, the length, the width and the height of the underwater detection region as shown in fig. 1 are respectively divided into 20, 20 and 4 blocks to form 1600 5 × 5 × 5 cube region blocks, the central point position of each region block is used as the possible occurrence position of the target source, and the search of the global fuzzy minimum characteristic value is carried out. Calculating the vector r from each magnetic sensor array element to the central point position i Obtaining a gain matrix A, selecting the measurement data D from the x direction w Calculated P k×k And a gain matrix A corresponding to the direction w Calculating a generalized eigenvalue λ (A) w H P k×k A w ,A w H A w ) Traversing all central points, sequencing feature values from small to large to obtain the minimum feature value of 0.596, finding a position point corresponding to the minimum feature value, and finding the optimal estimation position r after the fuzzy traversal process t1 =[7.5,-2.5,2.5]。
S332: the thinning traversal process is as shown in fig. 5, the area block where the target source is located obtained in the fuzzy traversal process is used as the target positioning area in the thinning traversal process, and on the basis, the target source is accurately positioned. Dividing a 5 multiplied by 5 cubic region block into 255 cubic region blocks of 1 multiplied by 1, continuously taking the central point of each region block as a point which can appear in a target, traversing all points according to a fuzzy traversal process to solve a minimum characteristic value d to obtain a minimum characteristic value of 0.565, finding out an optimal estimation position r obtained in a thinning traversal process t2 =[6.5,-0.5,2.5]。
S400: after fuzzy traversal and refinement traversal, the position r of the final positioning target magnetic dipole source is obtained t =[6.5,-0.5,2.5]And the positioning precision reaches 1m.
S500: comparing the actual magnetic dipole source position (7, 0, 2) with the magnetic dipole source position r estimated by the positioning algorithm t =[6.5,-0.5,2.5]The positioning precision within the range of the positioning area in the figure 1 reaches 1m, and the high-precision positioning requirement is met; the minimum characteristic value is solved twice through fuzzy traversal and detailed traversal, the optimal position of the target source is searched, and the fast optimization approach of the positioning algorithm path is realized; and magnetic induction intensity data with the signal-to-noise ratio of 20dB is received, the positioning error can be controlled within 1%, and the advantage of overcoming underwater interference is met.
In conclusion: the invention discloses an underwater three-dimensional space passive magnetic field positioning method based on an MP-MUSIC optimization algorithm, and belongs to the field of underwater passive magnetic field positioning. The invention mainly positions the underwater magnetic field source target, can overcome underwater interference and ensure that accurate positioning can be realized in a complex marine environment. The positioning method adopts the magnetic sensor receiving array, has the advantages of small volume, portability and arrangement, reduces the calculation complexity by optimizing the MP-MUSIC algorithm, realizes the quick optimization approach of the positioning algorithm path, realizes the accurate positioning in the complex large-scene marine environment by two processes of fuzzy positioning and refined positioning, and has the positioning accuracy reaching the meter order of magnitude in a detection area of one hundred meters. The technology can be applied to positioning of the electric leakage position of the submarine power supply cable, detection and tracking of the moving position of an underwater vehicle, the position of a frogman and the like, and has a great application value.

Claims (3)

1. An underwater three-dimensional space passive magnetic field positioning method based on an MP-MUSIC method is characterized by comprising the following steps:
(1) In the underwater environment in the positioning area range with the length of L meters, the width of W meters and the height of meter, the positioning targets are n magnetic dipole sources, and the magnetic dipole moment direction is
Figure FDA0003876741680000011
The magnetic dipole source is positioned
Figure FDA0003876741680000012
Arranging a circular magnetic sensor array underwater, wherein the number of the sensors is k, and the positions of sensor array elements are
Figure FDA0003876741680000013
The distance vector from the magnetic dipole source to the sensor array element is
Figure FDA0003876741680000014
(2) The magnetic induction intensity of the magnetic dipole source in the three-dimensional space measured by each magnetic sensor is as follows:
Figure FDA0003876741680000015
μ r for relative permeability, mu 0 Permeability of 4 pi x 10 in vacuum -7 H/m, N' is the number of turns of the coil of the closed current loop, I 0 For generating a closed current for a magnetic dipole source, R is the coil radius, n i Is waterLower interference, T represents matrix transposition, and gain matrix is defined as a i (r i ) Let B i =a i (r i )m+n i The magnetic induction includes intensities in three directions, i.e.
Figure FDA0003876741680000016
Wherein
Figure FDA0003876741680000017
(3) Storing the magnetic induction intensity data measured by the k magnetic sensor arrays into an array D = [ B = 1 T ,B 2 T ,…,B k T ] T Where B is i Is the magnetic induction intensity measured by the ith array element, D is a 3k multiplied by 1 matrix, and the array formed by the gain matrix is A = [ a ] 1 T (r 1 ),a 2 T (r 2 ),…,a k T (r k )] T Where a is i (r i ) Is the gain matrix corresponding to the ith array element, A is 3k × 3 matrix, the corresponding matrix form operation relationship is D = Am + N,
Figure FDA0003876741680000018
is an underwater disturbance;
(4) Solving the covariance matrix R DD =E{DD H Here R DD Is a 3k x 3k order matrix, H represents the conjugate transpose, E {. Cndot. } represents the desired operation, and R is obtained by matrix transformation DD =AE{mm H }A H2 I, wherein σ 2 I is a 3k multiplied by 3k order identity matrix for interference noise power;
(5) For covariance matrix R DD Decomposing the characteristic value to obtain R DD =U∑U H Where U is a 3k × 3k order eigenvector matrix and Σ is a diagonal matrix λ consisting of eigenvalues 1 ≥λ 2 ≥…λ 3k Is ≥ 0, defined by λ 12 ,…λ n Form a one-dimensional matrix sigma S From λ n+1n+2 ,…,λ 3k Form a one-dimensional matrix∑ N N magnetic dipole sources;
(6)∑ S the feature vector group corresponding to the medium feature value is defined as U S ,∑ N The characteristic vector group corresponding to the medium characteristic value is defined as U N Then, it can represent U = [ ] S ,U N ]Wherein U is S Is a signal subspace, U N Is a noise subspace, an orthogonal projection matrix P of the signal subspace is obtained 3k×3k ⊥=I-U S U S H
(7) Calculating the target function of the MP-MUSIC algorithm as the minimum lambda of the generalized characteristic decomposition min (A H P 3k×3k A,A H A);
(8) The process of traversing and solving the minimum characteristic value by the optimized MP-MUSIC algorithm is divided into a fuzzy traversal process and a refining traversal process, the traversal process traverses all positions in a positioning area and searches for the global minimum generalized characteristic value, and the steps are as follows:
(8-1) fuzzy traversal process: dividing the target location area into N × M × K volumes
Figure FDA0003876741680000021
The center point position of the region block
Figure FDA0003876741680000022
Representing the region, j is more than or equal to 1 and less than or equal to NxMxK, r tj Including position information of fuzzy traversal process, distance vector from magnetic dipole source to sensor array element
Figure FDA0003876741680000023
Traverse r t1 ,r t2 ,…,r tN×M×K Calculating a gain matrix A, substituting the gain matrix A into the generalized minimum eigenvalue lambda for solving different position points min (A H P 3k×3k A,A H A) The central point corresponding to the minimum characteristic value is the target position obtained by fuzzy traversal, and the area block where the position is located is the area block where the target is likely to appear in the whole situation;
(8-2) refining the traversal process: will blurThe target possible area block solved by the traversal process is used as a new target positioning area, and the area is subdivided into NN multiplied by MM multiplied by KK volumes
Figure FDA0003876741680000024
Traversing the position of the central point of the area block according to the step (8-1) to find the minimum characteristic value;
(8-3) when the refinement traversal process is finished, the minimum characteristic value is the global minimum generalized characteristic value;
(9) And the characteristic value corresponding point is the position of the target magnetic dipole source estimated by the MP-MUSIC algorithm, the algorithm is finished, and the positioning is completed.
2. The method for underwater three-dimensional passive magnetic field positioning based on MP-MUSIC method as claimed in claim 1, wherein in step (7), the objective function is λ min (A H P 3k×3k A,A H A) Therefore, the position of the target source can be located by finding the minimum feature value.
3. The MP-MUSIC method based underwater three-dimensional passive magnetic field positioning method of claim 1, wherein in step (8), the positioning is performed by fuzzy traversal process and refinement traversal process.
CN202011601964.3A 2020-12-30 2020-12-30 Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method Active CN112748376B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011601964.3A CN112748376B (en) 2020-12-30 2020-12-30 Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011601964.3A CN112748376B (en) 2020-12-30 2020-12-30 Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method

Publications (2)

Publication Number Publication Date
CN112748376A CN112748376A (en) 2021-05-04
CN112748376B true CN112748376B (en) 2022-12-13

Family

ID=75647120

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011601964.3A Active CN112748376B (en) 2020-12-30 2020-12-30 Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method

Country Status (1)

Country Link
CN (1) CN112748376B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107064862A (en) * 2017-02-22 2017-08-18 西安电子科技大学 Electromagnetic vector sensor linear array decorrelation LMS MUSIC method for parameter estimation

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702000A (en) * 2009-11-26 2010-05-05 华北电力大学(保定) Phased ultrasound locating method and system of partial discharge of transformer
JP5741115B2 (en) * 2011-03-25 2015-07-01 ソニー株式会社 POSITIONING DEVICE, POSITIONING METHOD, PROGRAM, AND RECORDING MEDIUM
JP5865822B2 (en) * 2012-04-17 2016-02-17 日東電工株式会社 Method for forming magnetic field space
US9812875B2 (en) * 2014-09-05 2017-11-07 Qualcomm Incorporated Systems and methods for adjusting magnetic field distribution using ferromagnetic material
US20200225655A1 (en) * 2016-05-09 2020-07-16 Strong Force Iot Portfolio 2016, Llc Methods, systems, kits and apparatuses for monitoring and managing industrial settings in an industrial internet of things data collection environment
CN108983145A (en) * 2018-08-27 2018-12-11 西安电子科技大学 Electromagnetic Vector Sensor Array Arrival Estimation of Wide-Band Coherent Source localization method
CN109633555A (en) * 2019-02-28 2019-04-16 哈尔滨理工大学 Submarine target localization method based on hereditary MUSIC algorithm

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107064862A (en) * 2017-02-22 2017-08-18 西安电子科技大学 Electromagnetic vector sensor linear array decorrelation LMS MUSIC method for parameter estimation

Also Published As

Publication number Publication date
CN112748376A (en) 2021-05-04

Similar Documents

Publication Publication Date Title
Xu et al. Optimal sensor placement for 3-D angle-of-arrival target localization
CN111504318B (en) Ocean navigation auxiliary method based on multi-magnetic dipole inversion
Yu et al. Multi-objectivization-based localization of underwater sensors using magnetometers
CN110132281B (en) Underwater high-speed target high-precision autonomous acoustic navigation method based on inquiry response mode
Zhao et al. Open-lake experimental investigation of azimuth angle estimation using a single acoustic vector sensor
Bereketli et al. Experimental results for direction of arrival estimation with a single acoustic vector sensor in shallow water
Pyo et al. Beam slice-based recognition method for acoustic landmark with multi-beam forward looking sonar
Xin et al. A TOA/AOA underwater acoustic positioning system based on the equivalent sound speed
Yang et al. A triangle matching algorithm for gravity-aided navigation for underwater vehicles
Sun et al. High-precision underwater acoustical localization of the black box based on an improved TDOA algorithm
CN111323752A (en) Far and near field transition interval sound source positioning method
CN112748376B (en) Underwater three-dimensional space passive magnetic field positioning method based on MP-MUSIC method
Deng et al. Underwater map-matching aided inertial navigation system based on multi-geophysical information
Hildebrandt et al. Combining cameras, magnetometers and machine-learning into a close-range localization system for docking and homing
CN111679248A (en) Target azimuth and distance combined sparse reconstruction positioning method based on seabed horizontal L-shaped array
CN111220146A (en) Underwater terrain matching and positioning method based on Gaussian process regression learning
CN110082611A (en) A kind of localization method of field measurement device
CN108629357B (en) Data fusion method and system for underwater vehicle
Huang An Underwater Target Tracking Algorithm Based on Extended Kalman Filter
Tian et al. Underwater target tracking based on the feature-aided GM-PHD method
Khalil et al. Multiple object localization in underwater wireless communication systems using the theory of gravitation
RU2559565C2 (en) Method of determining spatial position of extended objects located at depth, primarily under water, and electromagnetic line locator, primarily ship electromagnetic line locator for carrying out said method
Alexandri et al. A reverse bearings only target motion analysis (BO-TMA) for improving AUV navigation accuracy
Zhai et al. Research on an accurate passive location algorithm based on the underwater target acoustic-magnetic fusion in the near field
Xu et al. Placement of surface sensors for cooperative localization of submerged autonomous underwater vehicles

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