Distance Measurement Errors in Sile
Distance Measurement Errors in Sile
Abstract
The secretiveness of sonar operation can be achieved by using continuous frequency-modulated sounding signals
with reduced power and significantly prolonged repeat time. The application of matched filtration in the sonar
receiver provides optimal conditions for detection against the background of white noise and reverberation, and a
very good resolution of distance measurements of motionless targets. The article shows that target movement
causes large range measurement errors when linear and hyperbolic frequency modulations are used. The
formulas for the calculation of these errors are given. It is shown that for signals with linear frequency
modulation the range resolution and detection conditions deteriorate. The use of hyperbolic frequency
modulation largely eliminates these adverse effects.
Keywords: silent sonar, frequency modulation, continuous wave, matched filtering, distance measurement,
errors, Doppler effect.
© 2012 Polish Academy of Sciences. All rights reserved
1. Introduction
When used in military applications, sonars and radars should be difficult to intercept.
Stealth is a feature offered by silent radars using continuous signals and linear frequency
modulation [1-6]. Because they operate at significantly reduced power, compared to pulse
radars, the distance between the signal emitted by a silent radar and the enemy’s listening
system is shorter. Analogous silent sonars are not manufactured. The most likely reason for
that is the Doppler effect and the errors it causes in measuring the distance to moving targets
[7-9]. This is not an issue with radars, however, because the speed of electromagnetic wave
propagation is about 200 000 greater than the speed of acoustic wave propagation in water.
Despite that, literature on radars covers this aspect extensively, especially the resolution of
target distance measurements using pulse sonars and the negative effect of the Doppler effect
on detection capacity. These problems are analyzed using the narrowband ambiguity function
for selecting signals that will ensure the desired resolution and detection capacity [6, 10].
We will analyze errors of silent sonar distance measurements using the broadband
ambiguity function [11] because emitted acoustic signals feature a high relative bandwidth.
We will look at two types of signals, namely continuous signals with linear and hyperbolic
frequency modulation. Frequency modulated continuous wave (FM-CW) radars use the
Fourier transform of the product of the signal emitted and echo signal for detection purposes.
This article will cover detection with matched filtering. This type of detection is optimal when
known signals are received against the background of Gaussian noise [12]. We will use the
results of the analysis to determine the parameters of the silent sonar that will ensure that
distance measurement errors are acceptable.
________________________________________________________________________________________________________________________________________________________________________________
Article history: received on Feb. 6, 2012; accepted on Apr. 23, 2012; available online on May 18, 2012; DOI: 10.2478/v10178-012-0027-6.
J. Marszal, R. Salamon: DISTANCE MEASUREMENT ERRORS IN SILENT FM-CW SONAR WITH MATCHED FILTERING
The sonar’s transmitter emits a periodical signal x(t) which we will write down as:
x(t ) s (t nT ) , (1)
n 0
f
dfh
fh
dfl
fl
t
T T
0 Td Td
Fig. 1. The frequency of the emitted signal (solid line) and echo signal (dotted line).
The receiver is used for determining the correlation function of the emitted signal and echo
signal:
t
z (t ) s (t ' ) y * (t 't )dt' (6)
0
Metrol. Meas. Syst., Vol. XIX (2012), No. 2, pp. 321-332.
where fl=f0-B/2, fh= f0+B/2 are the upper and lower frequency of both signals respectively.
Fig. 2 shows the spectrum of signal s(t) with linear frequency modulation and Fig. 3 shows
the same with hyperbolic frequency modulation.
Fig. 2. Module of the LFM signal spectrum. Fig. 3. Module of the HFM signal spectrum.
For v=0 at the output of a matched filter the signal we obtain is proportional to the
periodically recurring autocorrelation function delayed by 0=2R/c. Fig. 4 shows a single
period of the LFM signal at the output of a matched filter and Fig. 5 shows a magnified
fragment of the signal. Numerical calculations were made for f0=10 kHz, B=2kHz and T=10 s.
The assumed target distance is R=3 km, and the propagation speed of the acoustic wave where
the time scale has been replaced with the distance scale is c=1500 m/s.
As you can see from the charts, the maximum of the matched filter output signal occurs
exactly at the assumed target distance. The ambiguity of the distance measurement is the
result of the correlation function width as shown in Fig. 5 and is inversely proportional to
signal spectrum width. For the data used in the charts, the ambiguity defined as the distance
between distance R=3 km and the contiguous zero of function z(R) is: R=c/2B=37.5 cm. In
practice, distances are never measured with this theoretical accuracy (which is usually not
necessary at all). This is primarily because we do not know the exact speed of sound in the
sea, acoustic wave propagation is not a straight line and noise or other factors do not occur
[13]. The speed of sound in the sea varies significantly ranging from about 1450 m/s to about
1540 m/s. As a result, when we use the mean value, the relative error in distance measurement
J. Marszal, R. Salamon: DISTANCE MEASUREMENT ERRORS IN SILENT FM-CW SONAR WITH MATCHED FILTERING
can be 3%. The error can be limited by taking frequent measurements of sound velocity
range resolutions in advanced underwater acoustic systems.
Fig. 4. LFM signal at the matched filter. Fig. 5. Magnified fragment of the signal in Fig. 4.
The errors in distance measurements attributed to the above factors are not confined to the
silent sonar only and can be seen in all underwater acoustic systems. Further in the paper we
will focus on analysing the predominant source of errors in silent sonar distance
measurements, i.e. the Doppler effect.
When analyzing the influence of the Doppler effect on errors in distance measurements, we
will use the broadband ambiguity function A(,d) defined as [11]:
A( , d ) | s(t ) s*[d (t ))dt | , (10)
where is the variance between delay 0 and the moment when the matched filter output
signal reaches its maximum. It can be demonstrated that constant delay 0 has no effect on the
results of the analysis [9].
If, for reasons of simplification, we neglect the difference between periods T and Td, then
by using relation (8), we have:
B B B B
A( , d ) | exp[ j 2 ( f 0 t )t ] exp{ j 2 [ f 0 d (t )]d (t )}dt | . (11)
2 2T 2 2T
Following the transformations the above formula reads as follows:
B B B
A( , d ) | exp{ j 2 (1 d )[ f
0
2 2T
(1 d )t )]t}exp( j 2 d 2t )dt | .
T
(12)
Let us note that the expression is the Fourier transform of a signal with linear frequency
modulation where the role of frequency f(,d) is taken on by the function:
B
f ( , d ) d 2 . (13)
T
As you can see in Fig. 2, the spectrum of the LFM signal is almost rectangular. By
differentiating the phase of the signal we obtain:
Metrol. Meas. Syst., Vol. XIX (2012), No. 2, pp. 321-332.
B B
f ( , d ) (1 d )[ f 0 (1 d )t )] , t=[0,T]. (14)
2 T
The boundary frequencies for time t=0 and t=T/d are:
B ,
f l ( , d ) (1 d )( f 0 ) (15)
2
B B .
f h ( , d ) (1 d )( f 0 ) (16)
2 d
By inserting relation (13), and following some simple transformations, we obtain the
boundary values of delay , with the relevant part of the ambiguity function contained
between them. They are respectively:
1 d T B 2v f 1
l 2
( f0 ) T ( 0 ) , (17)
d B 2 c B 2
1 d T B B 2v f 3
h 2
( f0 ) T ( 0 ) . (18)
d B 2 d c B 2
The mean value of the delay is equal to:
1 d T Bd 2v f 1
m 2
( f0 ) T( 0 ) (19)
d B 2 c B 2
and the width of the ambiguity function in cross-section d=const is:
1 d 2 4v
2
T T. (20)
d c
For velocity v=0 the width of the ambiguity function is not equal to zero because by
inserting (12) d=1 into the formula we obtain:
B
A( ,0) | exp( j 2 T t )dt | .
(21)
By integrating the exponential function within limits (0,T) and following some simple
transformation we obtain:
sin B .
A( ,0) | T | (22)
B
The width of the ambiguity function calculated between the first zeros of the above
function is:
2
min . (23)
B
This is the lowest width of the ambiguity function which applies to velocities that meet the
condition v<c/2BT. Let us note that an analogous result was obtained in the case of the
autocorrelation function shown in Fig. 5.
Fig. 6 shows the cross-sections of ambiguity functions for three values of parameter d and
B/f0=0.2, which were numerically determined using formula (10). By analogy, just as in the
previous analysis, for reasons of simplification the period of both functions was assumed to
be T=1 s. For d=0.99 the numerical values of delay and width of the ambiguity function are
given based on formulas (19) and (20).
J. Marszal, R. Salamon: DISTANCE MEASUREMENT ERRORS IN SILENT FM-CW SONAR WITH MATCHED FILTERING
Fig. 6. Cross-sections of the ambiguity function for Fig. 7. LFM signals at matched filter output.
the LFM signal (simplified version).
As you can see, increases in velocity v have three negative consequences, namely:
˗ the shift of the ambiguity function increases,
˗ the function becomes wider,
˗ the function drops.
Delay m produces errors in distance measurements. The mean error at moment in time t=T
when the measurement is taken, based on the formula (19) is equal:
c m vT f
Rm vT 0 . (24)
2 2 B
When the velocity v sign changes, so does the sign of distance error. If the target is moving
closer to the sonar, error Rm has a minus sign. The result is that the distance we have
determined is shorter than the actual distance. It is the reverse as the target moves away.
As the ambiguity function widens, the target distance becomes ambiguous. Using formula
(20), it can be expressed as:
c
R 2vT . (25)
2
The formula above also describes the sonar’s range resolution, i.e. the minimal difference
in the distance between two targets which can be discerned as separate objects.
As the ambiguity function drops, the signal to noise ratio deteriorates which makes target
detection less likely.
These effects provide a good basis for selecting the parameters of silent sonar. This will be
further discussed in Chapter 5.
The analysis so far is based on the simplified assumption of the sounding signal and echo
signal period being identical. As a result, simple formulas could be derived.
Fig. 7 illustrates what happens when we abandon this assumption. It shows matched filter
output signals, shifted in time by delay 0=2R/c. The velocities of the target and the other
parameters of the signals are the same as in Fig. 6. As you can see, the variance between the
periods of the sounding signal and echo signal has a significant impact on the signals. The
delays of the signals, however, which are important from the perspective of distance
measurement errors, are practically unchanged. There is a slight change in signal width but
this is the result of signal increase and the energy conservation law.
To recapitulate, the relations derived using the simplified method can be successfully used
for sonar design.
Metrol. Meas. Syst., Vol. XIX (2012), No. 2, pp. 321-332.
It is the auto-correlation function of the HFM signal contained in formula (9) with a shift
on the time axis equal -d. There is no relation between its maximal value and width and the
Doppler effect deviation which in this sense makes it Doppler effect-invariant [14, 15]. This is
illustrated in Fig. 8 which shows the ambiguity function for three values of the parameter d
and B/f0=0.2. Fig. 9 shows a magnified fragment of the ambiguity function calculated for
spectrum width B=2 kHz. As you can see the width of the cross-section is 1/B.
The error in target distance measurement is the result of delay d. By substituting the value
d=1+2v/c into formula (29) and using the relation (24) we obtain:
c d vT f
Rm vT 0 . (32)
2 2 B
As you can see, the distance measurement error using hyperbolic frequency modulation is
the same as for linear frequency modulation. The Doppler effect does not deteriorate detection
conditions because the ambiguity function’s maximal values are in no relation to target speed.
The sonar’s range resolution is also constant and equal to the resolution of a stationary target.
J. Marszal, R. Salamon: DISTANCE MEASUREMENT ERRORS IN SILENT FM-CW SONAR WITH MATCHED FILTERING
Fig. 8. Cross-sections of the ambiguity function for Fig. 9. Magnified fragment of the ambiguity function.
the HFM signal (simplified version).
Just as in the case of the LFM signal, the conclusions are derived from a simplified
analysis under the assumption that the duration of the HFM signal did not change as the target
was moving. Fig. 10 shows HFM signals at the output of a matched filter which take account
of the change in the function’s duration. A comparison with the signals in Fig. 8 reveals some
additional negative consequences of the Doppler effect with the pulses splitting and
decreasing. The first consequence deteriorates range resolution while the other affects the
conditions of detection. The maximal height of the signals depends on delay 0 (target
distance). It follows the rule that the shorter the delay the bigger the height. In the case of
shorter delays side signals diminish which does not, however, affect the distance between
them and the tall lines. This is illustrated in Fig. 11.
Fig. 10. HFM signal at matched filter Fig.11. HFM signal at matched filter
output (0=2.5 s, B/f 0=0.2). output (0=1 s, B/f 0=0.2).
The delay sign depends on the direction of speed and the changes in the sounding signal
frequency. In the examples above, the target was moving away from the sonar and frequency
was rising with time. Fig. 12 shows a situation in which ascending and descending frequency
signals are emitted consecutively. The first target located at distance R=1 km is moving away
from the sonar at v=3.75 m/s. The second target located at distance R=4 km is moving closer
to the sonar at v=7.5 m/s. As you can see, the signals are symmetrical in relation to the right
target distances. This can be used to achieve a significant reduction in distance measurement
errors as illustrated in Fig. 13. The targets are located at R=1.5 km (the bigger one) and R=3
km (the smaller one). Dotted lines show the signals as frequencies rise and fall and the solid
Metrol. Meas. Syst., Vol. XIX (2012), No. 2, pp. 321-332.
line is the result, i.e. the combination of both signals. As you can see, it occurs at target
locations. False echoes occurring between the actual positions of the targets can be eliminated
using signal selection which goes beyond the scope of this article.
Fig. 12. HFM signal at matched filter for rising Fig. 13. Result of distance measurement error
and falling frequencies (T= 10 s, B/f0=0.2). reduction (T= 10 s, B/f0=0.2).
The analysis above has shown that in the case of the silent sonar distance measurement
errors are primarily and most likely caused by the Doppler effect. The application and
parameters of the sonar depend on whether the errors are acceptable. Errors are caused by
external factors which are unknown and independent of the sonar’s parameters and by those
related to the sonar’s parameters which we can control during the design and operation of the
device. External factors include the speed and distance of targets versus the sonar. The second
group of factors includes the speed of the sonar’s carrier and two sonar parameters, i.e. the
duration of the sounding signal T and its relative spectrum width B/f0.
From the perspective of target speeds, silent sonars can be divided into those detecting
stationary objects (bottom, navigation obstacles, naval mines), objects moving at low speeds
(divers, underwater vehicles) and those designed to detect and follow underwater targets
(submarines, torpedoes).
Stationary target sonars are affected by the Doppler effect when their medium (ship,
underwater vehicle, helicopter) moves. Distance measurement errors in this case can be
reduced to acceptable limits by measuring the medium’s speed vector and using a rescaled
sounding signal for matched filtration, which is the signal s[d()t)], where d() is equal to:
2
d( ) 1 | v | cos . (33)
c
In the above formula | v | is the measured speed of the sonar’s carrier and is the angle
between the speed vector direction and the angle of the axis of a specific sonar beam. Modern
day sonars will usually produce a number of deflected narrow beams covering a wide angular
sector of simultaneous observation, [16]. Matched filtering in these sonars is performed for
each beam signal separately. When compensated, the sonar’s own speed does not increase the
number of numerical operations performed in real time. What is necessary, however, is an
additional generation of signals s[d()t)], an operation that must be performed with each
change of speed. Sets of the signals can be stored in computer memory.
The above method of compensating for the sonar’s own speed should be applied in the
majority of silent sonars, irrespective of what they are used for, because it significantly
J. Marszal, R. Salamon: DISTANCE MEASUREMENT ERRORS IN SILENT FM-CW SONAR WITH MATCHED FILTERING
reduces the negative effects of the Doppler effect. Exceptions include stationary sonars which
can be used as a security measure for harbours, roadstead, inland water structures, etc.
By selecting the right parameters silent sonars designed for diver and underwater vehicle
observation can ensure that distance measurement errors are within acceptable limits. This can
be achieved primarily thanks to slow speeds of the targets and small sonar ranges. The speed
of a diver carrying a load is v 0.5 m/s, and v 1 m/s without one. The speed of underwater
vehicles is hardly ever above v2 m/s. When used in the sonar’s transmitting array, typical
piezoelectric transducers can help emit signals that have the relative spectrum width of
B/f0=0.2. By substituting the value into formula (32), the distance measurement error can be
written down as:
R 5vT . (34)
Today’s multiple beam pulse sonars used in similar applications offer ranges Rm up to
several hundred meters. In a silent sonar duration T must meet the inequality:
T 2 Rm / c (35)
For range Rm=750 m, the minimal duration of T is 1 s. By inserting this value into
formula (34) it ranges from 2.5 m for v=0.5 m/s to 10 m for v=2 m/s. The same accuracy is
perfectly sufficient when the target is a diver or an underwater vehicle.
To ensure the secretiveness of the silent sonar period T should be as long as possible. What
is needed is a compromise between an acceptable error in distance reading and the required
secretiveness. Let us compare the silent sonar with the pulse sonar without matched filtering.
Let us assume that both sonars operate at a frequency of f0= 100 kHz, the silent sonar’s
transmitter emits HFM signals with a spectrum width of B=20 kHz and duration T=1s, and
the pulse sonar’s transmitter emits signals whose duration is i=0.1 ms. As a result, the
required bandwidth of the pulse sonar’s receiver is Bi=1/1=10 kHz.
Let us assume that the signal power at the input to the pulse sonar’s receiver is Pi, and P at
the input to the silent sonar’s receiver. Let us assume further that the required output signal-
to-noise ratio at the output of both receivers is identical and that the probability of detection
and false alarm is identical in both cases. The pulse sonar’s receiver detects the envelope
which means that the output signal to noise ratio is almost equal to the input signal to noise
ratio and amounts to:
Pi
SNRi . (36)
NBi
where N is the power spectral density of sea noise.
The output signal-to-noise ratio in a silent sonar with matched filtering is equal to:
PT
SNR . (37)
N
The above relations show that the signal power ratio at sonar outputs is equal to:
Pi
TBi . (38)
P
For the above values the quotient is 104 which means that the silent sonar’s range is
similar to that of the pulse radar but emits 104 less power. This is a slightly overestimated
figure due to the Doppler effect and its impact on the output signal.
With different powers emitted by the sonars in question, the ranges of their intercept by the
enemy passive system can vary significantly. Let us assume that the passive receiver’s input
signal power needed for intercept is identical for signals emitted by the pulse and silent sonar.
Metrol. Meas. Syst., Vol. XIX (2012), No. 2, pp. 321-332.
Fig. 14. Interception range of the 100kHz Fig. 15. Interception range of the 10 kHz
pulse sonar versus silent sonar pulse sonar versus silent sonar
(a- ocean, b- Baltic Sea, c- fresh water). (a- ocean, b- Baltic Sea, c- fresh water).
The tactical speed of submarines when under water ranges from 5 m/s (older ships) to 13
m/s (most recent ships). Combined with speeds like these, the long duration of T is the source
of serious distance measurement errors. Formula (34) shows that they will range from 500 m
to 1300 m for the above speeds. The error can be reduced by increasing the width of the
sounding signal spectrum, a method made possible by today’s broadband transducer
technology. The range of the silent sonar can also be reduced to ensure that it stays secretive.
As an example, for B/f0=0.5 and T=10, the errors can be reduced to 100 m and 260 m which is
acceptable for distances of ship detection. A radical improvement in distance measurement
J. Marszal, R. Salamon: DISTANCE MEASUREMENT ERRORS IN SILENT FM-CW SONAR WITH MATCHED FILTERING
accuracy can be achieved by applying algorithms using the echo signal symmetry as shown
in Fig. 13.
6. Conclusions
The main reason why silent sonars misread the distances is the Doppler effect. These errors
have acceptable values for high frequency short range sonars designed for detecting divers
and underwater vehicles. The high sound absorption makes the use of such sonar particularly
advantageous in waters with low salinity. The advantage of low frequency long-range silent
sonars over pulse sonars is that they are more secretive. The drawback, however, is that they
make serious mistakes in measuring distance, a feature which can be limited by using
additional post-detection processing of echo signals. This article is limited to the analysis of
the impact of the Doppler effect on the accuracy of distance measurement. This problem can
be considered more broadly, for example, using the method described in article [17].
Acknowledgements
This work was supported by the Ministry of Science and Higher Education, Poland (grant
ON517 624239).
References
[1] Kay, L. (1959). A comparison between pulse and frequency-modulation echo-ranging system. Journal
Brit.I.R.E., 19(2), 105-113.
[2] Kay, L. (1959). An experimental comparison between pulse and frequency-modulation echo-ranging
system. Journal Brit.I.R.E., 1960, 20(10), 785-796.
[3] Fuller, K.L. (1990). To see and not be seen. IEE Proceedings-F, 137(1), 1-10.
[4] Griffiths, H.D. (1990). New ideas in FM Radar. Electron. Commun. Eng. Journal, 2(5), 185-194.
[5] Skolnik, M. (1990). Radar Handbook. Second Edition. McGraw-Hill.
[6] Stove, A.G. (1992). Linear FMCW Radar Techniques. IEE Proceedings-F, 139(5), 343-350.
[7] Kramer, S.A. (1967). Doppler and acceleration tolerance of high-gain wide-band linear FM correlation
SONAR. Proc. IEEE, 33, 627-636.
[8] Marszal, J., Salamon, R., Zachariasz, K., Schmidt, A. (2011). Doppler effect in CW FM sonar.
Hydroacoustics, 14, 157-164.
[9] Salamon, R., Marszal, J., Schmidt, J., Rudnicki, M. (2011). Silent sonar with matched filtration.
Hydroacoustics, 14, 199-208.
[10] Levanon, N., Mozeson, E. (2004). Radar signals. John Wiley&Sons. New Jersey.
[11] Speiser, J.M. (1967). Wide-band ambiguity function. IEEE Trans. Inform. Theory, IP-13, 122-129.
[12] Lathi, B.P., Ding, Z. (2010). Modern digital and analog communication systems. Oxford University
Press. New York.
[13] Hodges R. P. (2010). Underwater acoustics. John Wiley&Sons.
[14] Kroszczyński, J.J. (1969). Pulse compression by means of linear-period modulation. Proc. IEEE 57,
1260-1266.
[15] Yang, J., Sarkar, T.K. (2006). Doppler-invariant property of hyperbolic frequency modulated waveform.
Microwave and optical technology letters, 48(8), 1174-1179.
[16] Salamon, R. (2006). Sonar systems. Wyd. GTN. Gdansk. (in Polish)
[17] Kowalczyk, A., Hanus, R., Szlachta, A. (2011). Investigation of the statistical method of time delay
estimation based on conditional averaging of delayed signal. Metrol. Meas. Syst., 18(2), 335-342.