CN107728109A - A kind of noncooperative target radiated noise measurement and positioning technology - Google Patents
A kind of noncooperative target radiated noise measurement and positioning technology Download PDFInfo
- Publication number
- CN107728109A CN107728109A CN201710837758.4A CN201710837758A CN107728109A CN 107728109 A CN107728109 A CN 107728109A CN 201710837758 A CN201710837758 A CN 201710837758A CN 107728109 A CN107728109 A CN 107728109A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- munder
- omega
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-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/22—Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of noncooperative target radiated noise measurement and positioning technology, based on vector hydrophone e measurement technology, acoustic pressure, the vibration velocity reception model of vector hydrophone are constructed, utilizes the acoustic pressure vibration velocity joint processing technology of vector hydrophone, using cross-spectrum sound intensity technique, the measurement azimuth of baseline is obtained；Target location coordinate information measured by every group of baseline is resolved by triangle Convergence method.Positioning for more vector hydrophones, the combination of two of vector hydrophone is can be considered, then again by Data Fusion technology, the measurement result of comprehensive all baselines, determine the position of each measurement point of moving target.Postpositive disposal finally is carried out using Kalman filtering algorithm, the track of moving target is further optimized.Data fusion technique combination Kalman filtering algorithm can improve positioning precision, and the movement locus of target is fast, accurately oriented in a small range, solve the problems, such as that two hydrophone positioning precisions are poor, tracking velocity is slower.
Description
Technical field
The present invention relates to noise testing technical field, and in particular to a kind of noncooperative target radiated noise measurement and positioning skill
Art.
Background technology
Noise testing technology (noise measuring technique) noise testing includes various noise sources and noise field
The measurement of fundamental characteristics parameter；The sound absorption used in Noise measarement and sound insulating material, the acoustical behavior measure of damping material；
Technology efficiency evaluation measurement of the control measure such as sound absorption, sound insulation, Amortization ＆＆ Damping, vibration isolation etc..In addition, research noise on human body
The subjective assessment for influenceing and endangering, being carried out to noise, formulate the work such as various Standards of Environmental Noise and noise testing offer is provided
The foundation of science.It is accurately finished these measurement needs of work and uses various technological means.
Noise testing technology can receive the radiated noise of target, obtain its acoustic information, can provide and refer to for vibration and noise reducing
Lead suggestion.And between submarine target and measurement point distance accurate measurement, be target sound source level calculate key.Passive sonar
Target ship-radiated noise is broadband, and the noise intensity on different frequency is different.Reflect sound source radiation noise intensity pair
The amount of frequency dependence, referred to as sound source spectrum level, be at 1 meter away from sound source radiated noise near certain single-frequency in 1HZ bandwidth
The sound intensity relative to the decibels of reference sound intensity, represented with SLs.Due to the complicated mechanism of ship-radiated noise, SLs is difficult to lead to
Theoretical calculation is crossed to obtain, it is necessary to carry out actual measurement.In conventional noise measuring system, it is to utilize cooperation beacon, obtains same pacing
Away from information.But it can not be realized for noncooperative target, this method.Because vector hydrophone simultaneous can obtain the scalar sum arrow of target
Information is measured, measurement capability can be lifted, used in noise measurement system.Therefore herein based on vector hydrophone measurement skill
Art, utilize vector hydrophone orientation estimated result, using data fusion method, the data message for obtaining more vector hydrophones
It is fully utilized, quickly and accurately orients the movement locus of target in a small range, solve domestic noise testing at present
The problem of noncooperative target can not being accurately positioned in system.
Deng Xiu China et al. have studied method (the Deng Xiu China self-conductance canister missdistances measurement side that synchronic distance accurately measures
Method research marine electronic engineerings, 2012, Vol.32No.10), it is provided with synchronous acoustic marker in target, can not realize to non-
The measurement of cooperative moving targets；Wu Yan group et al. have studied vector hydrophone bearings-only target motion analysis method (Wu Yanqun,
Waterborne target motion analysis acoustic techniques of the Hu Yongming based on single vector hydrophone, 2010,1000-3630 (2010) -04-
0361-04), it can carry out orientation estimation using single vector hydrophone to target, a wide range of interior to moving target progress
Tracking measurement, but this method tracking velocity is slower, positioning precision can not meet measurement request.
Based on this, the present invention combines vector hydrophone e measurement technology with data fusion method, obtains to a greater degree
Target and the information content of environment must be tested, the target information that comprehensive more set measuring systems obtain, effectively improves system
Performance, improve positioning precision.On the other hand combine Kalman filtering algorithm is filtered processing to target motion conditions again,
Optimized the movement locus of target, it is met the needs of noise testing uses.Method proposed by the present invention make use of
Data fusion technique has handled the data message that more vector hydrophones measure, and the data for making to measure have played the work of maximum
With being adapted to a small range is quick, the high-precision movement locus for orienting target, improve the efficiency of experiment, in engineering
It is easier to realize.
The content of the invention
The present invention proposes a kind of noncooperative target radiated noise measurement and positioning technology, combines vector data fusion and card
Kalman Filtering algorithm, effectively non-cooperative moving targets can be tracked with measurement, obtain kinematic parameter, preferably improve vector
The positioning precision of hydrophone and the stability of measuring system.
The technical solution adopted for the present invention to solve the technical problems comprises the following steps：
(1) the line measurement battle array model being made up of four array elements, square measurement battle array model, four array elements composition six are established respectively
Bar baseline, they with six triangles of target configuration, form 12 azimuths respectively；
(2) receipt signal model of vector hydrophone is established, obtains the acoustic pressure data P (t) that vector hydrophone receives, X side
To vibration velocity Vx(t), Y-direction vibration velocity Vy(t), by the acoustic pressure vibration velocity joint processing technology of vector hydrophone, using cross-spectrum sound intensity technique
Carry out orientation estimation；
(3) spectrum analysis is carried out to reception signal, extraction envelope obtains effective band, and frequency domain is carried out in narrow bandwidth range and is melted
Conjunction is handled；
(4) 6 baselines and 12 obtained azimuths are combined, is crossed method by triangle, is resolved every baseline and surveyed
The target location obtained；The target location calculated to multigroup vector hydrophone carries out data and merged again；
(5) to the moving object measurement track of generation in (4), it is optimized using Kalman filtering algorithm.
The acoustic pressure vibration velocity joint processing technology of the vector hydrophone of described step (2) specifically includes：
Target is combined with six baselines respectively, obtaining the reception signal at i-th of array element is：
Wherein subscript s represents semaphore, and subscript n represents noisiness；pi(t) the acoustic pressure letter that i-th of array element receives is represented
Number, vxi(t) the vibration velocity signal in the horizontal direction that i-th of array element of expression receives, vyi(t) represent what i-th of array element received
Vibration velocity signal in vertical direction, θ are the horizontal azimuth of incident acoustic wave.
The cross-spectrum sound intensity technique of described step (2) specifically includes：
The acoustic pressure amount p that will be obtainedi(r, t), vibration velocity amount vxi(r,t)、vyi(r, t) makees Fourier transformation, is carried out in frequency domain
Signal transacting can obtain frequency domain sound intensity information：
The horizontal azimuth for estimating each frequency using cross-spectrum sound intensity technique is：
Wherein ω is angular frequency,For the average sound intensity in x directions,For the average sound intensity in y directions,
For the estimate of horizontal azimuth.
Described step (3) specifically includes：
Power Spectral Estimation is carried out to whole frequency domain with FFT first, the line spectrum of signal is found by analysis, to extraction of line spectrum
Envelope, to a narrow bandwidth range (f where envelope1,f2,f3,…,fn) do Frequence zooming analysis；For same echo signal,
A series of azimuth estimation value of targets is obtained using cross-spectrum sound intensity technique；The target location then calculated to multigroup vector hydrophone
Carry out data to merge again, to the weighted comprehensive for being accurately positioned result and should be all baseline positioning results of target location
Weight selected is variance counting backward technique, wherein DiRepresent the variance of i-th group of measurement data：
Described step (4) specifically includes：
The distance between array element i and array element j：
Bi_hydrophone cross bearing schematic diagram as shown in figure 1, solve measurement battle array coordinate system in R andFor：
Described step (5) specifically includes：
The recurrence formula of Kalman filtering algorithm is as follows：
P (k+1 | k)=Φ P (k | k) Φ '+Γ Q (k) Γ '；
K (k+1)=P (k+1 | k) H ' (k+1 | k) S-1(k+1)；
S (k+1)=H (k+1) P (k+1 | k) H ' (k+1)+R (k+1)；
P (k+1 | k+1)=[I-K (k+1) H (k+1)] P (k+1 | k)；
Wherein, Q (k) δkl=E [g (k) .g'(l)], R (k) δkl=E [w (k) .w'(l)].
The beneficial effects of the present invention are：This method sufficiently make use of the data message that vector hydrophone obtains, and will swear
Amount hydrophone acoustic pressure vibration velocity joint processing technology and Data fusion technique are effectively combined together with.And employ Kalman's filter
Ripple algorithm, the track of moving target is optimized so that the positioning precision of moving target is significantly improved, can
To meet the measurement request to noncooperative target, there is stronger engineering practicability.
Brief description of the drawings
Fig. 1 is flow chart of the present invention；
Fig. 2 is double vector hydrophone cross bearing schematic diagrames；
Fig. 3 is square cloth station schematic diagram；
Fig. 4 is target trajectory simulation analysis result；
Fig. 5 is X-direction velocity measuring simulation analysis result；
Fig. 6 is Y-direction velocity measuring simulation analysis result；
Fig. 7 is X-axis position detection error simulation analysis result；
Fig. 8 is Y-axis position detection error simulation analysis result.
Embodiment
The present invention is further described with example below in conjunction with the accompanying drawings.
(1) measurement model of four-vector hydrophone is established, four array element combination of two form six measurement baselines, and they divide
Not with six triangles of target configuration, azimuths are measured so as to obtain 12.
So that square measures battle array as an example, square array is respectively positioned in xoy planes with moving target, four-vector hydrophone
[1,2,3,4] is distributed on square four summit that the length of side is a, is arranged counterclockwise, and coordinate is respectively (xi,yi),
I=1,2,3,4.Array element 1 and array element 2 are located in x-axis, and equidistantly distributed is in origin O both sides.Array element 3 and array element 4 are located at respectively
The surface of array element 2, array element 3.The initial position of target is located at array element 1 and the midpoint of the intersection of array element 4, as shown in Figure 3.
(2) receipt signal model of vector hydrophone is established, obtains the acoustic pressure data P (t) that vector hydrophone receives, X side
To vibration velocity Vx(t), Y-direction vibration velocity Vy(t).Make target be combined respectively with six baselines, then can obtain at i-th of array element
Reception signal is：
Subscript s represents semaphore, and subscript n represents noisiness.Wherein pi(t) the acoustic pressure letter that i-th of array element receives is represented
Number, vxi(t) the vibration velocity signal in the horizontal direction that i-th of array element of expression receives, vyi(t) represent what i-th of array element received
Vibration velocity signal in vertical direction.θ is the horizontal azimuth of incident acoustic wave.
The acoustic pressure amount p that will be obtainedi(r, t), vibration velocity amount vxi(r,t)、vyi(r, t) makees Fourier transformation, is carried out in frequency domain
Signal transacting can obtain frequency domain sound intensity information：
The horizontal azimuth for estimating each frequency using cross-spectrum sound intensity technique is：
ω is angular frequency in formula,For the average sound intensity in x directions,For the average sound intensity in y directions,
For the estimate of horizontal azimuth.
(3) to improve positioning precision, spectrum analysis need to be carried out to reception signal, effective band is obtained, in narrow bandwidth range
Frequency domain fusion treatment is carried out, carrys out positioning result of all baselines of composite measurement battle array to target.
Power Spectral Estimation is carried out to whole frequency domain with FFT first, the line spectrum of signal is found by analysis, to extraction of line spectrum
Envelope, to a narrow bandwidth range (f where envelope1,f2,f3,…,fn) do Frequence zooming analysis.
Vector hydrophone carries out direction finding using cross-spectrum sound intensity technique, corresponding to each frequency f of reception signal1,f2,
f3,…,fn, an azimuth information can be estimated according to formula (2.4)Therefore, for same echo signal, utilize
Cross-spectrum sound intensity technique can obtain a series of azimuth estimation value of targetsFor each group of vector
Hydrophone, corresponding to each frequency of reception signal, one group of target location coordinate data S can be calculated according to formula (8)1
(x,y),S2(x,y),S3(x,y),…,Sn(x, y), utilize the size of each frequency sound intensityIt is each to integrate
The target location result of frequency.
The target location then calculated to multigroup vector hydrophone carries out data and merged again, to the accurate fixed of target location
Position result should be the weighted comprehensive of all baseline positioning results, i.e.,：
Weight selected is variance counting backward technique, wherein DiRepresent the variance of i-th group of measurement data：
(4) six baselines are combined into by four array elements, respectively with six triangles of target configuration, form 12 azimuths.Profit
With the distance between obtained level orientation value θ and baseline, the position coordinates where target is determined by triangle Convergence method.
Specially：
The distance between array element i and array element j：
Bi_hydrophone cross bearing schematic diagram as shown in figure 1, solve measurement battle array coordinate system in R andFor：
(5) processing is filtered to target motion conditions using Kalman filtering.Kalman filtering algorithm (recurrence formula)
It is as follows：
P (k+1 | k)=Φ P (k | k) Φ '+Γ Q (k) Γ '
K (k+1)=P (k+1 | k) H'(k+1 | k) S-1(k+1)
S (k+1)=H (k+1) P (k+1 | k) H'(k+1)+R (k+1)
P (k+1 | k+1)=[I-K (k+1) H (k+1)] P (k+1 | k)
Wherein：Q(k).δkl=E [g (k) .g'(l)], R (k) δkl=E [w (k) .w'(l)].
The embodiment of content of the invention each several part is illustrated above, combines Data fusion technique and karr
More vector hydrophones joint track and localization technology of graceful filtering algorithm, it can effectively improve the positioning precision of system.Below with just
Exemplified by square cloth station, simulation result is analyzed.
Instance parameter sets as follows：Quaternary square array is laid as shown in Figure 2.The horizontal position that No.1 vector hydrophone is laid
It is set to (- 100,0)；The horizontal level that No. two vector hydrophones are laid is (100,0), the horizontal position that No. three vector hydrophones are laid
It is set to (100,200)；The horizontal level that No. four vector hydrophones are laid is (- 100,200).Assuming that target is simple signal, letter
Number frequency is 150Hz, sample rate 4096, and initial time signal to noise ratio is 20dB, the initial position of target in the horizontal plane for (-
100,100), the initial velocity of X-direction is Vx=1m/s, Y-direction initial velocity are Vy=0.
Fig. 3 is target trajectory simulation analysis result, and Fig. 4 is velocity measuring simulation analysis result, and Fig. 5 detects for position
Errors simulation analysis result.
The simulation result of complex chart 3, Fig. 4 and Fig. 5 is understood：
(1) multivariate vector hydrophone Passive Positioning computation, the movement locus of target can be relatively accurately depicted,
Demonstrate the reliability of the algorithm and the validity of localization method.More vector hydrophone combination Kalman filtering algorithms, can enter one
Step improves positioning precision.
(2) the distance between target and hydrophone can be accurately obtained using azimuth information, carries out the calculating of Acoustic Wave Propagation,
So as to accurately obtain the radiated noise level of target.
Claims (6)
1. a kind of noncooperative target radiated noise measurement and positioning technology, it is characterised in that specifically comprise the following steps：
(1) the line measurement battle array model being made up of four array elements, square measurement battle array model are established respectively, and four array elements form six bases
Line, they with six triangles of target configuration, form 12 azimuths respectively；
(2) receipt signal model of vector hydrophone is established, obtains the acoustic pressure data P (t) that vector hydrophone receives, X-direction is shaken
Fast Vx(t), Y-direction vibration velocity Vy(t), by the acoustic pressure vibration velocity joint processing technology of vector hydrophone, carried out using cross-spectrum sound intensity technique
Estimate in orientation；
(3) spectrum analysis is carried out to reception signal, extraction envelope obtains effective band, carried out in narrow bandwidth range at frequency domain fusion
Reason；
(4) 6 baselines and 12 obtained azimuths are combined, is crossed method, is resolved measured by every baseline by triangle
Target location；The target location calculated to multigroup vector hydrophone carries out data and merged again；
(5) to the moving object measurement track of generation in (4), it is optimized using Kalman filtering algorithm.
2. a kind of noncooperative target radiated noise measurement and positioning technology according to claim 1, it is characterised in that described
The acoustic pressure vibration velocity joint processing technology of the vector hydrophone of step (2) specifically includes：
Target is combined with six baselines respectively, obtaining the reception signal at i-th of array element is：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>x</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>y</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
Wherein subscript s represents semaphore, and subscript n represents noisiness；pi(t) sound pressure signal that i-th of array element receives, v are representedxi
(t) the vibration velocity signal in the horizontal direction that i-th of array element of expression receives, vyi(t) Vertical Square that i-th of array element receives is represented
Upward vibration velocity signal, θ are the horizontal azimuth of incident acoustic wave.
3. a kind of noncooperative target radiated noise measurement and positioning technology according to claim 1, it is characterised in that described
The cross-spectrum sound intensity technique of step (2) specifically includes：
The acoustic pressure amount p that will be obtainedi(r, t), vibration velocity amount vxi(r,t)、vyi(r, t) makees Fourier transformation, and signal is carried out in frequency domain
Processing can obtain frequency domain sound intensity information：
<mrow>
<mover>
<msub>
<mi>I</mi>
<mi>x</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mrow>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<mover>
<msub>
<mi>I</mi>
<mi>y</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mrow>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
The horizontal azimuth for estimating each frequency using cross-spectrum sound intensity technique is：
<mrow>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>tan</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>(</mo>
<mfrac>
<mrow>
<mover>
<msub>
<mi>I</mi>
<mi>y</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<msub>
<mi>I</mi>
<mi>x</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>&omega;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
<mo>;</mo>
</mrow>
Wherein ω is angular frequency,For the average sound intensity in x directions,For the average sound intensity in y directions,For water
Equal azimuthal estimate.
4. a kind of noncooperative target radiated noise measurement and positioning technology according to claim 1, it is characterised in that described
Step (3) specifically includes：
Power Spectral Estimation is carried out to whole frequency domain with FFT first, finds the line spectrum of signal by analysis, to extraction of line spectrum bag
Network, to a narrow bandwidth range (f where envelope1,f2,f3,…,fn) do Frequence zooming analysis；For same echo signal, profit
A series of azimuth estimation value of targets is obtained with cross-spectrum sound intensity technique；The target location then calculated to multigroup vector hydrophone is entered
Row data merge again, to the weighted comprehensive for being accurately positioned result and should be all baseline positioning results of target location
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</munder>
<munder>
<mi>&Sigma;</mi>
<mi>j</mi>
</munder>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</munder>
<munder>
<mi>&Sigma;</mi>
<mi>j</mi>
</munder>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</munder>
<munder>
<mi>&Sigma;</mi>
<mi>j</mi>
</munder>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
</mrow>
</munder>
<munder>
<mi>&Sigma;</mi>
<mi>j</mi>
</munder>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
Weight selected is variance counting backward technique, wherein DiRepresent the variance of i-th group of measurement data：
<mrow>
<msub>
<mi>&omega;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>/</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>.</mo>
</mrow>
5. a kind of noncooperative target radiated noise measurement and positioning technology according to claim 1, it is characterised in that described
Step (4) specifically includes：
The distance between array element i and array element j：
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>&NotEqual;</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Bi_hydrophone cross bearing schematic diagram as shown in figure 1, solve measurement battle array coordinate system in R andFor：
6. a kind of noncooperative target radiated noise measurement and positioning technology according to claim 1, it is characterised in that described
Step (5) specifically includes：
The recurrence formula of Kalman filtering algorithm is as follows：
P (k+1 | k)=Φ P (kk) Φ '+Γ Q (k) Γ '；
K (k+1)=P (k+1 | k) H'(k+1k) S-1(k+1)；
S (k+1)=H (k+1) P (k+1 | k) H'(k+1)+R (k+1)；
P (k+1 | k+1)=[I-K (k+1) H (k+1)] P (k+1 | k)；X (k+1 | k+1)=X (k+1 | k)+K (k+1 | k)+K
(k+1)·V(k+1)；
Wherein, Q (k) δkl=E [g (k) .g'(l)], R (k) δkl=E [w (k) .w'(l)].
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710837758.4A CN107728109A (en) | 2017-09-18 | 2017-09-18 | A kind of noncooperative target radiated noise measurement and positioning technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710837758.4A CN107728109A (en) | 2017-09-18 | 2017-09-18 | A kind of noncooperative target radiated noise measurement and positioning technology |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107728109A true CN107728109A (en) | 2018-02-23 |
Family
ID=61206597
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710837758.4A Pending CN107728109A (en) | 2017-09-18 | 2017-09-18 | A kind of noncooperative target radiated noise measurement and positioning technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107728109A (en) |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108614268A (en) * | 2018-04-26 | 2018-10-02 | 中国人民解放军91550部队 | The acoustics tracking of low altitude high speed airbound target |
| CN109283491A (en) * | 2018-08-02 | 2019-01-29 | 哈尔滨工程大学 | A kind of unmanned plane positioning system based on vector probe unit |
| CN110196407A (en) * | 2019-05-06 | 2019-09-03 | 东南大学 | A kind of single vector hydrophone signal arrival bearing's estimation method based on frequency estimation |
| CN110244260A (en) * | 2019-06-17 | 2019-09-17 | 杭州电子科技大学 | Submarine target high-precision DOA estimation method based on acoustic energy flow vector compensation |
| CN110703199A (en) * | 2019-10-22 | 2020-01-17 | 哈尔滨工程大学 | Quaternary cross array high-precision azimuth estimation method based on compass compensation |
| CN111427010A (en) * | 2020-04-20 | 2020-07-17 | 中国电子科技集团公司电子科学研究院 | ROV underwater positioning system and positioning method |
| CN111427011A (en) * | 2020-04-20 | 2020-07-17 | 中国电子科技集团公司电子科学研究院 | Submarine asset position calibration method and system |
| CN111708006A (en) * | 2020-05-28 | 2020-09-25 | 哈尔滨工程大学 | Target line spectrum detection method suitable for unmanned platform detection sonar |
| CN112415467A (en) * | 2020-11-06 | 2021-02-26 | 中国海洋大学 | Single-vector subsurface buoy target positioning implementation method based on neural network |
| CN112816940A (en) * | 2020-12-23 | 2021-05-18 | 中国船舶重工集团有限公司第七一0研究所 | Target distance estimation method and device based on sound pressure and particle vibration velocity |
| CN113204009A (en) * | 2021-05-12 | 2021-08-03 | 深圳康佳电子科技有限公司 | Running step length reminding method and device, terminal and computer readable storage medium |
| CN113671442A (en) * | 2021-07-30 | 2021-11-19 | 青岛海纳水下信息技术有限公司 | Underwater unmanned cluster navigation positioning method based on vector hydrophone technology |
| CN113702960A (en) * | 2021-06-29 | 2021-11-26 | 哈尔滨工程大学 | High-precision speed measurement method for underwater mobile platform based on time delay and Doppler frequency shift |
| CN113777594A (en) * | 2021-08-11 | 2021-12-10 | 上海船舶电子设备研究所(中国船舶重工集团公司第七二六研究所) | Method and system for testing minimum detectable signal-to-noise ratio of nonlinear sonar |
| CN115079088A (en) * | 2022-06-10 | 2022-09-20 | 杭州电子科技大学 | Target DOA estimation method based on frequency domain acoustic energy flow instantaneous phase difference weighting |
| CN115792806A (en) * | 2022-10-24 | 2023-03-14 | 哈尔滨工程大学 | Non-cooperative line spectrum distributed underwater sound positioning method |
| CN117807356A (en) * | 2024-02-29 | 2024-04-02 | 齐鲁工业大学(山东省科学院) | Double-vector hydrophone positioning method based on improved sparrow algorithm optimized particle filtering |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040145968A1 (en) * | 2003-01-29 | 2004-07-29 | John Brittan | Method for processing dual sensor seismic data to attenuate noise |
| CN101650220A (en) * | 2009-09-04 | 2010-02-17 | 合肥工业大学 | Method for correcting finite difference error of cross-spectrum sound intensity |
| CN102226837A (en) * | 2011-04-08 | 2011-10-26 | 哈尔滨工程大学 | Vector circle array acoustic pressure and vibration velocity combined direction finding method on cylindrical form baffle condition |
| CN202329798U (en) * | 2011-11-16 | 2012-07-11 | 中国船舶重工集团公司第七一五研究所 | Two-dimensional vector hydrophone based on piezoelectric ceramic |
| CN103454616A (en) * | 2013-08-27 | 2013-12-18 | 西北工业大学 | Method for estimating orientation of cross type velocity gradient hydrophone |
| CN103605108A (en) * | 2013-07-29 | 2014-02-26 | 哈尔滨工程大学 | High-precision remote direction estimation method of acoustic vector array |
| CN105589066A (en) * | 2015-12-14 | 2016-05-18 | 西北工业大学 | Method for estimating parameters of underwater constant-speed vehicle based on vertical vector array |
| CN106680762A (en) * | 2016-12-15 | 2017-05-17 | 哈尔滨工程大学 | Sound vector array orientation estimation method based on cross covariance sparse reconstruction |
-
2017
- 2017-09-18 CN CN201710837758.4A patent/CN107728109A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040145968A1 (en) * | 2003-01-29 | 2004-07-29 | John Brittan | Method for processing dual sensor seismic data to attenuate noise |
| CN101650220A (en) * | 2009-09-04 | 2010-02-17 | 合肥工业大学 | Method for correcting finite difference error of cross-spectrum sound intensity |
| CN102226837A (en) * | 2011-04-08 | 2011-10-26 | 哈尔滨工程大学 | Vector circle array acoustic pressure and vibration velocity combined direction finding method on cylindrical form baffle condition |
| CN202329798U (en) * | 2011-11-16 | 2012-07-11 | 中国船舶重工集团公司第七一五研究所 | Two-dimensional vector hydrophone based on piezoelectric ceramic |
| CN103605108A (en) * | 2013-07-29 | 2014-02-26 | 哈尔滨工程大学 | High-precision remote direction estimation method of acoustic vector array |
| CN103454616A (en) * | 2013-08-27 | 2013-12-18 | 西北工业大学 | Method for estimating orientation of cross type velocity gradient hydrophone |
| CN105589066A (en) * | 2015-12-14 | 2016-05-18 | 西北工业大学 | Method for estimating parameters of underwater constant-speed vehicle based on vertical vector array |
| CN106680762A (en) * | 2016-12-15 | 2017-05-17 | 哈尔滨工程大学 | Sound vector array orientation estimation method based on cross covariance sparse reconstruction |
Non-Patent Citations (3)
| Title |
|---|
| HU BO ET AL.: "Underwater patch near-field acoustical holography based on particle velocity and vector hydrophone array", 《SCIENCE CHINA PRESS AND SPRINGER-VERLAG BERLIN HEIDELBERG》 * |
| 孙勇 等: "多基地声纳系统定位精度分析与最优布站", 《计算机仿真》 * |
| 莫世奇: "矢量水听器的数据融合研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108614268A (en) * | 2018-04-26 | 2018-10-02 | 中国人民解放军91550部队 | The acoustics tracking of low altitude high speed airbound target |
| CN108614268B (en) * | 2018-04-26 | 2021-12-07 | 中国人民解放军91550部队 | Acoustic tracking method for low-altitude high-speed flying target |
| CN109283491A (en) * | 2018-08-02 | 2019-01-29 | 哈尔滨工程大学 | A kind of unmanned plane positioning system based on vector probe unit |
| CN110196407A (en) * | 2019-05-06 | 2019-09-03 | 东南大学 | A kind of single vector hydrophone signal arrival bearing's estimation method based on frequency estimation |
| CN110196407B (en) * | 2019-05-06 | 2022-06-24 | 东南大学 | Single-vector hydrophone signal incoming wave direction estimation method based on frequency estimation |
| CN110244260B (en) * | 2019-06-17 | 2021-06-29 | 杭州电子科技大学 | Underwater target high-precision DOA estimation method based on acoustic energy flow vector compensation |
| CN110244260A (en) * | 2019-06-17 | 2019-09-17 | 杭州电子科技大学 | Submarine target high-precision DOA estimation method based on acoustic energy flow vector compensation |
| CN110703199A (en) * | 2019-10-22 | 2020-01-17 | 哈尔滨工程大学 | Quaternary cross array high-precision azimuth estimation method based on compass compensation |
| CN110703199B (en) * | 2019-10-22 | 2021-09-24 | 哈尔滨工程大学 | Quaternary cross array high-precision azimuth estimation method based on compass compensation |
| CN111427011A (en) * | 2020-04-20 | 2020-07-17 | 中国电子科技集团公司电子科学研究院 | Submarine asset position calibration method and system |
| CN111427010A (en) * | 2020-04-20 | 2020-07-17 | 中国电子科技集团公司电子科学研究院 | ROV underwater positioning system and positioning method |
| CN111708006A (en) * | 2020-05-28 | 2020-09-25 | 哈尔滨工程大学 | Target line spectrum detection method suitable for unmanned platform detection sonar |
| CN112415467A (en) * | 2020-11-06 | 2021-02-26 | 中国海洋大学 | Single-vector subsurface buoy target positioning implementation method based on neural network |
| CN112816940A (en) * | 2020-12-23 | 2021-05-18 | 中国船舶重工集团有限公司第七一0研究所 | Target distance estimation method and device based on sound pressure and particle vibration velocity |
| CN112816940B (en) * | 2020-12-23 | 2023-06-06 | 中国船舶重工集团有限公司第七一0研究所 | Target distance estimation method and device based on sound pressure and particle vibration velocity |
| CN113204009A (en) * | 2021-05-12 | 2021-08-03 | 深圳康佳电子科技有限公司 | Running step length reminding method and device, terminal and computer readable storage medium |
| CN113204009B (en) * | 2021-05-12 | 2023-12-22 | 深圳康佳电子科技有限公司 | Running step length reminding method, running step length reminding device, terminal and computer readable storage medium |
| CN113702960A (en) * | 2021-06-29 | 2021-11-26 | 哈尔滨工程大学 | High-precision speed measurement method for underwater mobile platform based on time delay and Doppler frequency shift |
| CN113702960B (en) * | 2021-06-29 | 2022-10-21 | 哈尔滨工程大学 | High-precision speed measurement method for underwater maneuvering platform based on time delay and Doppler frequency shift |
| CN113671442A (en) * | 2021-07-30 | 2021-11-19 | 青岛海纳水下信息技术有限公司 | Underwater unmanned cluster navigation positioning method based on vector hydrophone technology |
| CN113777594A (en) * | 2021-08-11 | 2021-12-10 | 上海船舶电子设备研究所(中国船舶重工集团公司第七二六研究所) | Method and system for testing minimum detectable signal-to-noise ratio of nonlinear sonar |
| CN115079088A (en) * | 2022-06-10 | 2022-09-20 | 杭州电子科技大学 | Target DOA estimation method based on frequency domain acoustic energy flow instantaneous phase difference weighting |
| CN115792806A (en) * | 2022-10-24 | 2023-03-14 | 哈尔滨工程大学 | Non-cooperative line spectrum distributed underwater sound positioning method |
| CN115792806B (en) * | 2022-10-24 | 2024-02-20 | 哈尔滨工程大学 | Non-cooperative line spectrum distributed underwater sound positioning method |
| CN117807356A (en) * | 2024-02-29 | 2024-04-02 | 齐鲁工业大学(山东省科学院) | Double-vector hydrophone positioning method based on improved sparrow algorithm optimized particle filtering |
| CN117807356B (en) * | 2024-02-29 | 2024-05-10 | 齐鲁工业大学(山东省科学院) | Double-vector hydrophone positioning method based on improved sparrow algorithm optimized particle filtering |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107728109A (en) | A kind of noncooperative target radiated noise measurement and positioning technology | |
| CN102004244B (en) | Doppler direct distance measurement method | |
| CN104035065A (en) | Sound source orienting device on basis of active rotation and method for applying sound source orienting device | |
| CN104569625B (en) | A kind of large-scale antenna directional diagram measuring method based on rotatable auxiliary antenna | |
| CN107607943A (en) | The high method of survey of delay Doppler radar altimeter based on interferometric phase auxiliary | |
| CN103217211A (en) | Substation noise source distribution measuring method based on synthetic aperture principle | |
| CN108614268A (en) | The acoustics tracking of low altitude high speed airbound target | |
| CN105866636B (en) | Transformer substation positioning method based on time difference positioning | |
| CN101566683A (en) | Method for estimating target azimuth formed based on phase difference beams | |
| CN102540177B (en) | Target positioning method based on 3D ray-tracing algorithm | |
| CN105954746A (en) | Landform correction meter wave radar height measurement method based on broadcast automatic mutual supervisory signals | |
| CN103744052A (en) | Dual-satellite time difference measurement direction-finding method and apparatus for aerial target positioning | |
| CN103558602A (en) | Simulated annealing locating method for multi-base sonar configuration mode | |
| Sun et al. | Array geometry calibration for underwater compact arrays | |
| CN103529428A (en) | Method for passively positioning point sound source based on spatial ten-element array | |
| CN110456304A (en) | Airborne DF and location method | |
| CN108594193A (en) | A kind of radar system bias estimation method based on fixed target and noncooperative target | |
| Chen et al. | TDOA/FDOA mobile target localization and tracking with adaptive extended Kalman filter | |
| CN105824019A (en) | Optimized beam alignment method for large distributed space surveillance radar | |
| CN107202975A (en) | A kind of a burst of first attitude error rectification method of two-dimensional vector | |
| CN105572637A (en) | Far-field sound source positioning system and method | |
| CN111460362A (en) | Sound source positioning data complementation method based on quaternary microphone array group | |
| CN106646413A (en) | Radar networking vertical line crossing integration positioning method and error calculating method thereof | |
| Sun et al. | Underwater asynchronous navigation using single beacon based on the phase difference | |
| Xu et al. | A Simulation method for USBL Localization Performance Analysis using a Stereo Array |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180223 |
|
| RJ01 | Rejection of invention patent application after publication |