Lin-wei Tao, Ying-min Wang. Target localization based on joint measurement of amplitude and frequency in a LOFAR field[J]. Journal of Zhejiang University Science A, 2014, 15(2): 130-137.
@article{title="Target localization based on joint measurement of amplitude and frequency in a LOFAR field", author="Lin-wei Tao, Ying-min Wang", journal="Journal of Zhejiang University Science A", volume="15", number="2", pages="130-137", year="2014", publisher="Zhejiang University Press & Springer", doi="10.1631/jzus.A1300263" }
%0 Journal Article %T Target localization based on joint measurement of amplitude and frequency in a LOFAR field %A Lin-wei Tao %A Ying-min Wang %J Journal of Zhejiang University SCIENCE A %V 15 %N 2 %P 130-137 %@ 1673-565X %D 2014 %I Zhejiang University Press & Springer %DOI 10.1631/jzus.A1300263
TY - JOUR T1 - Target localization based on joint measurement of amplitude and frequency in a LOFAR field A1 - Lin-wei Tao A1 - Ying-min Wang J0 - Journal of Zhejiang University Science A VL - 15 IS - 2 SP - 130 EP - 137 %@ 1673-565X Y1 - 2014 PB - Zhejiang University Press & Springer ER - DOI - 10.1631/jzus.A1300263
Abstract: To estimate the motion parameters of a moving target before its passing by the closest point of approach (CPA) point in a low frequency analyzing and recording (LOFAR) field, an error-free theoretical method based on the joint measurement of target radiated noise’s amplitude and frequency was presented. First, the error-free theoretical equations for target characteristic frequency, absolute velocity, the CPA, and amplitude of the radiation noise were derived by three equal interval measured values of the target amplitude and frequency. Then, the method to improve the calculation accuracy was given. Finally, the simulation and experiments were conducted in the air and showed the correctness of this method. By using one single piece of LOFAR, this method can calculate four target parameters and the relative error of each estimated parameter is less than 10%.
Darkslateblue:Affiliate; Royal Blue:Author; Turquoise:Article
Article Content
1. Introduction
Low frequency analyzing and recording (LOFAR) is a very important application in the passive measurement field. Generally, LOFAR with no directivity can only passively receive radiated noise of the target. Line spectrum information contained in the noise is recorded and processed to calculate the motion parameters of the target. With characteristics of low cost, small volume, light weight, etc., the LOFAR is applied widely in various fields, for example, underwater target localization (Wu and Sun, 1999; Ling and Wang, 2007).
Methods with single pieces of LOFAR for target parameter estimation are as follows. The classical algorithm for LOFAR, namely algorithm of Doppler-closest point of approach (CPA) was proposed in (Tao 2009). Tao (2009) calculated the velocity of the target and CPA approximately by recording the two symmetrical target frequencies of the CPA point. An improved algorithm of the Doppler-CPA was also proposed in (Tao, 2009), where the calculation of the absolute velocity of the target and the CPA error-free theoretical formula were deduced by introducing time information. All the methods above can only be applied after the target passes by the CPA point, which is limited in use. The approximate calculation method of the target speed and CPA was presented in (Yu, 2011), which used any of the three point of frequency measurements and information of the change rate of the target line spectrum. Actually, the Doppler variation of the target is originally a slow process with noise. Therefore, using the rate of change to calculate the target motion parameters may cause some large errors.
Target localization with joint multiple LOFAR is another active direction, such as the classical LOFAR FIX (LOFIX) principle and analysis of localization accuracy in (Hu et al., 2009), the classical hyperbolic fixing (HYFIX) principle and analysis of localization accuracy in (Sun et al., 2010). The passive range estimation was completed by using a vector sensor array (Li et al., 2012). The passive parameter estimation and its accuracy were analyzed from different aspects (Shang et al., 1985; Mao, 2001; Brian and Kam, 2002; Xue et al., 2005; Ma, 2007; Cockrell and Schmidt, 2010). To improve the accuracy of the target frequency measurement, some methods of image processing were applied in previous studies (Fan et al., 2002; Yang et al., 2010).
This paper focuses on using a single piece of LOFAR. The radiated noise’s amplitude of the target at any time (not restricted to passing by the CPA necessarily) is collected by three types of equal intervals. The error-free calculation formulas for the four types of status information consisted of absolute velocity of the target, characteristic frequency of the target, the closest distance, and the absolute amplitude, are then derived.
2. Measurement model
2.1. Measurements of frequency
As shown in Fig. 1, the target moves with speed v, and the CPA is D. The characteristic frequency of the line-spectrum is fT, ti is measure time, ri is distance from LOFAR to target, and φi is the relative bearing between the moving forward direction of the target and LOFAR in the ith measurement, i=1, 2, 3.
Fig.1 Motion diagram between target and LOFAR
By the Doppler shift equation, the line-spectrum frequency of the detected target by using LOFAR is (Tao and Wang, 2008)
,
where fi is the ith measured target frequency, and c is the speed of sound.
2.2. Measurements of amplitude
LOFAR is mainly used for recording and processing the strong line spectrum in the noise radiated by the target. According to (Liu and Lei, 1993), the line spectrum noise mainly consists of mechanical noise and propeller noise in the background noise of the target. The frequency ranges from 10 Hz to 2000 Hz.
Regarding the target as a point source, according to the spherical wave theory of sound propagation, the transmission loss is
,
where α is the absorption coefficient. The first item of the transmission loss is extended loss, and the second item is sound absorption loss.
The empirical formula of the low frequency band absorption coefficient which is derived from Liu and Lei (1993) is
.
The absorption coefficient is 0.000279–0.12 in the frequency band of 10–2000 Hz. If the distance is 10 km, the extended loss will be 80 dB, and the maximum absorption loss is 1.2 dB. Therefore, the absorption loss can be ignored when the frequency is in the low band.
Assuming that there are only extended loss and expansion of spherical waves, the amplitude data of the target noise which is received by LOFAR are (Liu and Lei, 1993)
,
where p is the pressure of sound, and A is the amplitude of sound pressure which is equivalent to the unit distance from the target.
2.3. Equation of observation
As shown in Fig. 1, the frequency and amplitude of the target are measured and recorded at equal interval times t1, t2, and t3 (t2−t1=t3−t2=Δt). The measurement of frequency is shown in Eq. (1).
The measurement of amplitude is
,
where pi is the signal amplitude of the target in the ith measurement, i=1, 2, 3. AT is the absolute amplitude which is equivalent to the unit distance from the target. According to the geometric principle,
.
3. Calculation of the target’s parameters
3.1. Characteristic frequency
Eq. (1) can be rewritten as
,
and Eq. (5) can be rewritten as
,
then
.
Substituting Eq. (9) into Eq. (6), and substituting the measured values of t1, t2, and t3, we have
.
Because f1, f2, f3 and p1, p2, p3 have been obtained, according to Eq. (10) the characteristic frequency fT of the target can be resolved as follows:
.
3.2. Absolute velocity
According to the principle of trigonometric function and using the measured values of t1, t2, and Eq. (7)2+Eq. (8)2=1, we have
Eq. (12) can also be expressed in a matrix:
.
The velocity v and D/AT can be solved as fT has been known.
3.3. Closest point of approach (CPA)
When fT, v, and D/AT have been known, tgφi can be calculated on the basis of Eq. (9). Then the closest point of approach D can be derived from Eq. (6) by use of the two values tgφ1 and tgφ2. Finally, the absolute amplitude AT at the distance which is equivalent to the unit distance from the target can be resolved.
4. Algorithm improvement
Only three measured data have been used in the method mentioned above. To improve the calculation accuracy in practice, the best results can be obtained by using the statistical method with more data.
4.1. Improvement of target characteristic frequency
For the characteristic frequency of the target, the arithmetic average of multi-point calculation results fT1, fT2, …, fTn that can be obtained from Eq. (11) is used as the final characteristic frequency.
.
4.2. Improvement of target velocity
Using the multi-point data, Eq. (13) can be rewritten as
.
For , , and ,
then in the least squares sense, there is
.
5. Computer simulation
5.1. Error-free simulation
First of all, the correctness of the formula will be verified by use of the computer simulation without consideration of any error.
We set the target speed of 5 m/s, characteristic line spectrum of 1000 Hz, radiation sound source level of 134 dB, and the equivalent voltage of 1 m away from the target AT=501.90 V (which is calculated on the condition that the sensitivity of acoustic receiving transducer is −203 dB, and the gain of the preamplifier circuit is 120 dB).
LOFAR is located in the origin of the coordinates, and the closest distance D to the target is 100 m, 1000 m, and 10000 m. Matlab is used as the simulation software. The sample rate of the frequency and amplitude is 1 Hz, sound speed is 1500 m/s, and time zero is set at the moment that the target passes by the CPA point (the CPA point in Fig. 1). The information of Doppler and amplitude of the targets which contain −300 s, −200 s, and −100 s is used for the calculation.
The calculated results of the simulation model and sample data without errors are shown in Table 1.
Because this algorithm is theoretically accurate and error-free, the calculation error depends entirely on the error of the calculation process. Matlab acted as the numerical simulation software due to its high calculation accuracy. In the actual process, the error of each calculated result can reach 10−14 orders of magnitude. Table 1 does not list all the significant figures as it is limited by the table size. As in the case of 1000.00 m, the actual calculated velocity is 4.99999999999998, and the distance is 1000.0000000000024. The correctness of this algorithm can be proved from the simulation results.
Table 1
Calculated results without error
Actual closest distance (m)
Calculated results
Frequency (Hz)
Velocity (m/s)
Closest distance (m)
Absolute amplitude (V)
100.00
999.99
5.00
100.00
501.90
1000.00
1000.00
4.99
999.99
501.90
10000.00
1000.00
4.99
9999.99
501.90
5.2. Simulation with error
Parameters of simulation with error are consistent with those of error-free simulation in section 5.1. The measurement error of frequency and amplitude is given below. Liu and Lei (1993) showed that in the deep sea the marine environment noise level is about 55 dB near the frequency of 1000 Hz, and the equivalent noise voltage is 0.0563 V (the transducer sensitivity is −203 dB, and the gain of the preamplifier circuit is calculated as 120 dB). In the simulation, the amplitude of the measured signal is given as follows:
,
where p(k) is the magnitude of the signal at time k, r(k) is the distance from the target to the LOFAR at time k, and nA is the Gaussian white noise with mean square error of 0.0536 V and with a mean of zero.
At present, the detection of the frequency is mainly using an adaptive line spectrum enhancer, which can greatly improve the frequency measurement under the colored noise. Its measurement accuracy can reach 0.1 Hz, which is given by
,
where f(k) is the measured value of the frequency at time k, fT(k) is the real value of the target’s frequency at time k, and nF is the Gaussian white noise with a mean square error of 0.1 and mean of zero.
Eqs. (14) and (16) are applied when calculating the target parameters. The calculated results and the relative errors of the calculated results are shown in Table 2.
Table 2
Calculated results with noise
Actual closest distance (m)
Calculated results
Frequency (Hz)
Velocity (m/s)
Closest distance (m)
Absolute amplitude (V)
100.00
1000.97 (0.10%)
5.02 (0.40%)
101.77 (1.77%)
495.89 (1.22%)
1000.00
998.96 (0.10%)
5.20 (4.00%)
970.82 (2.91%)
483.26 (7.12%)
10000.00
1003.32 (0.33%)
6.33 (26.60%)
8505.25 (14.95%)
590.25 (17.58%)
Data in parenthesis are relative errors
As shown in Table 2, the calculated results are very accurate at close range (1000.00 m). The maximum relative error of distance is 2.91%, which can meet the actual requirements. Considering that we are only using a single piece of passive LOFAR, such results are still satisfactory.
The error increases dramatically at far distance, e.g., in the case of 10000.00 m, the calculated error of the velocity is up to 26.60%. The main reason is that the variation of the target’s frequency and amplitude is reduced with the increasing distance, and the variation is almost the same as that of the amplitude of the noise. Therefore, the calculated error is very large, and may even result in a singular matrix which cannot be calculated.
6. Experiment in the air
6.1. Principle of the experiment
A straight road was selected to be the experiment site, and a car was used to simulate the moving target with high-speed. The experimental equipment was divided into two parts: signal generation and acquisition. The car was used as the carrier of the signal generator which can produce a single frequency signal. The signal through a power amplifier would drive the high-power speaker to simulate the target with a single line spectrum. The signal acquisition part consisted of a high-sensitivity microphone, filter amplifier, spectrum analyzer, data logger, and oscilloscope. The microphone picked up the ambient noise in the air and the characteristic frequency of the target. After filtering and amplification, the data was finally collected and stored by the data acquisition instrument. The spectrum analyzer and oscilloscope as the monitoring devices were used to observe the noise signal in the frequency domain and time domain (Fig. 2).
Fig.2 Diagram of experiment plan
In the experiment (Fig. 2), the car traveling at a constant speed on the highway can simulate the targets which passed by the microphone, namely the CPA point at a uniform speed. The signal acquisition part collected the environmental noise and the target line spectrum which were stored by the data logger.
The experiment content was mainly to change the different distances of the CPA (D in Fig. 2) and the car speed v, and then the collected sound data were used to calculate the frequency f of the signal, the absolute amplitude A, the speed v of the car, and the CPA distance D. The correctness of the algorithm and the effects of different factors on the algorithm will be verified in the experiment.
6.2. Processing method
The collected noise spectrum is a time domain signal. To estimate the frequency of each time accurately, the short-time Fourier transform method is used. In each fast Fourier transform (FFT), the length is to be 256 sampling points, the FFT overlap rate is 50%, and the length of the FFT algorithm is 1024. When the FFT of each time period is obtained, the frequency point with the maximum amplitude will be chosen as the measured frequency in this time period. The calculated frequency signal is to be handled by smooth filtering, namely using a 32-order finite impulse response (FIR) low-pass filter whose low-pass cut off frequency is 100 Hz.
6.3. Experiment results and analysis
In the process of the experiment, the microphone sensitivity is −45 dB, and the microphone frequency response band is 50 Hz–18 kHz. The gain of the amplifier is 60 dB, and the pass band frequency range of the filter is 0–20 000 Hz. The sampling rate of the data logger is 50 000 Hz. Four different CPAs are selected in the experiment, i.e., 7, 11.6, 15, and 20 m (due to the environmental limitation of the actual road, the maximum distance is 20 m). The speed of the car is respectively selected as 20, 40, 60, and 80 km/h. The frequency of the signal is 1000 Hz. The battery is chosen to be the energy source to drive the signal whose amplitude remains constant with the root mean square (RMS) of 15.55 VRMS. The speed of sound in the air is taken as c=340 m/s while processing the data.
Some typical processed results are shown in Fig. 3. Fig. 3 is the time-frequency cumulative results, in which the target passes by the CPA point whose distance is 11.60 m, the speed is respectively set as 20, 40, 60, and 80 km/h, and the characteristic frequency is 1000 Hz. The higher the brightness of the figure is, the greater the amplitude of the signal means. When the vehicle approaches, the frequency of the signal is greater than the actual frequency (1000 Hz), and the Doppler is positive. When the car departs, the frequency of the signal is decreasing, and the Doppler is negative.
Fig.3 Processed results (a)-(d) of time-frequency
The turning points at the frequency curve indicate the time when the target passed by the CPA point. It can be seen that the change rate of the CPA point is greater when the speed of the vehicle is higher.
There are totally 13 times valid processes in the entire experiment. The frequency and the amplitude of the signal are handled on the computer. The results are shown in Table 3.
Table 3
Calculated parameters in air experiment
Serial No.
Actual distance (m)
Actual frequency (Hz)
Actual velocity (km/h)
Actual amplitude of the signal (VRMS)
Calculated distance (m)
Calculated frequency (Hz)
Calculated velocity (km/h)
Calculated amplitude (VRMS)
1
7.0
1000
80
15.55
7.03
1002.98
78.04
15.54
2
7.0
1000
60
15.55
6.87
1001.25
62.68
15.20
3
7.0
1000
40
15.55
6.95
999.50
41.91
15.98
4
11.6
1000
20
15.55
11.01
996.25
20.04
16.01
5
11.6
1000
40
15.55
11.46
1002.05
38.13
15.67
6
11.6
1000
60
15.55
10.24
997.43
57.43
14.28
7
11.6
1000
80
15.55
11.50
1001.24
79.29
15.64
8
15.0
1000
80
15.55
15.57
995.63
77.68
16.06
9
15.0
1000
60
15.55
16.05
1001.45
58.36
14.90
10
15.0
1000
40
15.55
15.33
1002.54
38.46
15.90
11
20.0
1000
60
15.55
21.98
1001.20
58.03
15.99
12
20.0
1000
40
15.55
18.74
999.35
39.79
16.58
13
20.0
1000
80
15.55
18.38
998.59
77.22
14.04
The experiment has achieved good results (Table 3). Most of the processes (except the experiment No. 6) estimate the motion parameters of the target well. The relative error of the calculated distance has an average value of 3.95% and the maximum value of 9.90%, while the relative error of the calculated frequency has an average value of 0.19% and the maximum value of 0.38%. The relative error of the calculated velocity has an average value of 2.85% and the maximum value of 4.77%, and the relative error of the calculated amplitude has an average value of 3.19% and the maximum value of 9.71%.
According to the analysis of the experiment results, the closest distance has the greatest impact on the calculation error, which is consistent with the results of the computer simulation. The greater the distance is, the larger the calculation error will be. The most accurate estimated parameter is frequency, due to the frequency being stronger in anti-interference ability than other parameters. At the same time, there is no error propagation in the calculation process, when this method is directly used to measure the frequency. The estimation accuracy of velocity is higher than that of the distance and absolute amplitude, because the velocity is directly due to the Doppler frequency shift. The accuracy of the velocity will decrease, if the accuracy of the frequency is lower. The estimation accuracy of the closest distance and the absolute amplitude is the worst one. The distance calculation depends on the change of the frequency. When the distance is far, the relative velocity is reduced, and the change rate of the frequency decreases, which is equal to the reduction of the signal noise ratio (SNR). Therefore, the accuracy of the calculated results is reduced. Moreover, for the amplitude of the signal, the external interference is added directly to the amplitude information, leading to poor estimation accuracy. For example, in the experiment of No. 6, a heavy truck passed by, so some strong broadband interference was introduced in the vicinity of the low frequency band of 300–500 Hz at the time of 0.075 s–0.18 s. The strong interference, shown in Fig. 3c, resulted in the reduction of the calculation accuracy of the parameter.
7. Conclusions
In this paper, radiated noise’s amplitude information of the target was added to the parameter calculation process of LOFAR firstly. In theory, the error-free formulas were derived by using the amplitude and frequency information of the target, through any three equal interval measurements. Four important parameters, that is the speed of the target, the closest distance, characteristic frequency, and absolute amplitude, were obtained.
The method does not need to overcome the disadvantages of the previous methods. The calculation is completed without the requirement of the CPA point. It improves the calculation speed and efficiency of the LOFAR, and makes the tactics more flexible. Moreover, it is an error-free method in theory, and has important theoretical significance in the algorithm research of the LOFAR fields.
The simulation and experimental results in air show the validity of the method. Relative error of each estimated parameter is less than 10%. Moreover, this method only uses one single piece of LOFAR, with the accurate estimation of the target distance, speed, frequency, and amplitude at close range. Hence, it has a very good application prospects for many engineering fields, such as road speed measurement, namely using a single microphone to measure the speed of a vehicle and amplitude level of the noise. Further research will focus on the case of how to improve the estimation accuracy of the target parameter over a long distance.
[1] Brian, G.F., Kam, W.L., 2002. Passive ranging errors due to multipath distortion of deterministic transient signals with application to the localization of small arms fire. The Journal of the Acoustical Society of America, 111(1):117-128.
[2] Cockrell, K.L., Schmidt, H., 2010. Robust passive range estimation using the waveguide invariant. The Journal of the Acoustical Society of America, 127(5):2780-2789.
[3] Fan, Y.Y., Tao, B.Q., Xiong, K., 2002. Feature extraction of ship radiated noise by 1½-spectrum. Acta Acustica, 02:137-145.
[4] Hu, Z.X., Sun, M.T., Su, W.G., 2009. Simulation analysis of LOFIX fixing accuracy for passive omnidirectional sonobuoy. Electronics Optics & Control, 16(12):26-29.
[5] Li, J., Sun, G.Q., Han, Q.B., 2012. Acoustics vector sensor linear array passive ranging based on waveguide invariant. , Proceedings of the 3rd International Conference on Ocean Acoustics, 576-586. :576-586.
[6] Ling, G.M., Wang, Z.M., 2007. Sonobuoy technology and its development direction. Acoustic and Electronic Engineering, (in Chinese),3:1-5.
[7] Liu, B.S., Lei, J.Y., 1993. Principle of Underwater Acoustics. (in Chinese), Harbin Engineering University Press,Harbin, China :
[8] Ma, J.G., 2007. Passive Localization Technology of Time Reversal. PhD Thesis, (in Chinese), Harbin Engineering University,China :
[9] Mao, W.N., 2001. An overview of passive localization for under-water acoustics. Journal of Southeast University, 31(6):1-4.
[10] Shang, E.C., Clay, C.S., Wang, Y.Y., 1985. Passive harmonic source ranging in waveguides by using mode filter. The Journal of the Acoustical Society of America, 78(1):172-175.
[12] Tao, L.W., 2009. Study on Array Sonobuoy Signal Processing System and Algorithm. PhD Thesis, (in Chinese), Northwestern Polytechnical University,Xian, China :
[13] Tao, L.W., Wang, Y.M., 2008. New algorithm for sonobuoy Doppler-CPA. Journal of System Simulation, 20(23):6353-6355.
[14] Wu, Y.F., Sun, N.H., 1999. The key equipment of submarine: present situation and development of sonobuoy. Technical Acoustics, 18:95-96.
[15] Xue, S.H., Ye, Q.H., Huang, H.N., 2005. Passive-range estimation using near-field MVDR dual focused beamformers. Applied Acoustics, 24(3):177-181.
[16] Yang, P., Yuan, B.C., Zhou, S., 2010. A line spectrum estimation method of underwater target radiated noise base on the 1½D spectrum. , International Conference on Innovative Computing and Communication and Asia-Pacific Conference on Information Technology and Ocean Engineering, 297-299. :297-299.
[17] Yu, T., 2011. Doppler direct location for underwater target. Journal of China Academy of Electronics and Information Technology, (in Chinese),6(3):328-330.
Open peer comments: Debate/Discuss/Question/Opinion
Open peer comments: Debate/Discuss/Question/Opinion
<1>