GB2191649A - Radio direction-finding - Google Patents

Abstract

A method and a system for radio direction-finding by measuring the Time of Arrival (ToA) of the leading edge of signals from a distant source at two relatively closely spaced receiving elements. The times at which the instantaneous detected amplitudes of the received signals first exceed a minimal threshold value just above the noise level are measured, and this means that the measured time is not affected by multipath which involves more than a few metres additional path length for the indirect, delayed signal. A suitable timing circuit is disclosed (Fig. 2). By making ToA measurements on three coplanar, non-collinear receivers, directions of incidence in three dimensions can be determined (Figs. 4,5). <IMAGE>

Description

SPECIFICATION Radio direction-finding The invention relates to a method and to a system for determining the direction of incidence of electromagnetic wave signals from a distant source by determining the difference between the time of arrival (ToA) of the leading edge of wave signals received respectively from the source at the two elements of at least one pair of mutually-spaced wave-receiving elements.

The invention further relates to a timing circuit for measuring time with a high resolution, for example one nanosecond.

ToA (also known as Time Difference of Arrival, TDOA) direction-finding (DF) with a long baseline, i.e.

wherein the wave-receiving elements are spaced miles apart, is known from, for example, the article "Passive Direction Finding and Signal Location" by A.R. Baron et al, Microwave Journal, September 1982, p'ages 59-76: see particularly page 59 and pages 6670. A major disadvantage of ToA DF using a long baseline in many practicai situations is that if, as is normally the case, the direction of incidence is to be determined over a substantial range of directions, there is a substantial interval of time over which signals from the same source may arrive at one wave-receiving element of the pair relative to the other, the difference between the ToAs depending on the position of the source relative to the pair of elements.If there is a plurality of sources from which signals may be received, for example sources emitting pulsed signals with a substantial pulse repetition frequency (PRF), then there is a significantly probability that the signals whose ToAs at the two elements are compared come from different sources rather than the same source; the greater the spacing between the elements, and consequently the longerthe above-mentioned period, the greater is the probability. It is then necessary to compare one or more characteristic parameters of the signals received at the two elements, for example frequency, PRF or pulse length, to ascertain whether they come from the same or different sources. Not only does this require substantial further equipment, but it substantially increases the time taken to ascertain the direction of incidence of the signals.

ToA DF using a short baseline, for example 24 feet, is known from US patent 3 936 831. The use of a short baseline has the advantage (although there is no mention of it in the US patent) that the above-mentioned interval within which signals from the same source can arrive at the two elements is so small that there is a high probability in practical situations that signals from a different source will not arrive in that period.

However, the use of a short baseline imposes the difficulty of determining time differences of the order of tens of nanoseconds or less. The above-mentioned US patent proposes a system wherein a capacitor is charged at a fast, linear rate from the constant-current source, charging being started by the arrival of a pulse signal at one element and stopped by the arrival of a pulse signal at the other element; the time difference is then effectively multiplied by transferring the voltage on the capacitor to a further capacitor which is discharged at a much slower constant rate. However, the circuitry disclosed for performing these functions would not in practice be suitable for the very short time differences involved.For example, the current from the constant-current source could not be switched between zero and its full value in a time which is short compared with the time difference involved. Moreover, Figure 5 of the US patent, which is a graph of a count representative of the measured time difference against time delay (actual time difference), shows a predominantly substantially rectilinear variation from about 2000 nanoseconds down to about 150 nanoseconds; at this point, there is an abrupt change of slope, with what appears to be a hypothetical extrapolation to the origin of the graph. This indicates that the circuitry would not in fact operate as intended for the time differences of 0-50 nanoseconds that would actually need to be measured.

It is an object of the invention to provide an improved method and system for short-baseline ToA DF. It is a further object of the invention to provide an improved timing circuit.

According to a first aspect of the invention, a method of determining the direction of incidence of electromagnetic wave signals from a distant source from the time of arrival of the leading edge of the wave signals comprises: receiving said signals at a plurality of mutually spaced wave-receiving elements, detecting the respective instantaneous amplitude of the signal received at each element, measuring the times at which the detected amplitudes of wave signals received respectively at at least two of said elements first exceed a minimal threshold value such that signals can be satisfactorily distinguished from noise and which is substantially less than the minimum peak value of signals whose direction of incidence is to be determined by said method, the time being measured in such a manner that the measured time is generally unaffected by multipath propagation, determining the difference between the measured times in respect of one pair or of a plurality of pairs of said elements, wherein the two elements of said one pair of each of at least two of said plurality of pairs are sufficiently close together that the length of the interval of time within which signals from the same source must arrive at the two elements is so short that there is a high probability in operation that no signals from another source will arrive in that interval, and deriving a representation of the direction of incidence from the time difference(s) utilising the relationship cos oz=cAt/d where a is the angle between the direction of incidence of the signals and the line joining the two elements of a said pair, d is the distance between those two elements, 3t is the time difference in respect of that pair of elements, and c is the free-space velocity of electromagnetic waves.

According to a second aspect of the invention, a system for determining the direction of incidence of electromagnetic wave signals from a distant source from the time of arrival of the leading edge of the wave signals, comprises: a plurality of mutually spaced wave-receiving elements, means for detecting the respective instantaneous amplitude of the signal received at each element, means for measuring the times at which the detected amplitudes of wave signals received respectively at at least two of said elements first exceed a minimal threshold value such that signals can be satisfactorily distinguished from noise and which is substantially less than the minimum peak value of signals whose direction of incidence is to be determined by said method, the time being measured in such a manner that the measured time is generally unaffected by multipath propagation, means for determining the difference between the measured times in respect of one pair or of a plurality of pairs of said elements, wherein the two elements of said one pair or of each of at least two of said plurality of pairs are sufficiently close together that the length of the interval of time within which signals from the same source must arrive at the two elements is so shortthatthere is a high probability in operation that no signals from another source will arrive in that interval, and means for deriving a representation of the direction of incidence from the time difference(s) utilising the relationship cos (x=côt/d where a is the angle between the direction of incidence of the signals and the line joining the two elements of a said pair, d is the distance between those two elements, 3t is the time difference in respect of that pair of elements, and c is the free-space velocity of electromagnetic waves.

These aspects of the invention involve the recognition that in contrast to other methods and systems for direction-finding such as interferometry, ToA DF can provide a good degree of immunity to multipath propagation, provided that the time of arrival is measured early on the leading edge of the signal and in such a manner that the measured time is not substantially affected by multipath propagation which involves more than short additional path length for the indirect, delayed signal. If the ToA is not measured early on the leading edge of a signal, multipath can significantly distort the shape of the leading edge as the signal approaches its peak value, leading to a substantial error in a measured time difference and hence in the derived direction of incidence.As will be described below, the ToA may be measured in respect of a threshold which is well below the minimum peak amplitude and in such a manner that the measured time is unaffected by variations in signal amplitude within a period of, for example, not more than about 10 nanoseconds after the threshold is first exceeded, so that the system is immune to multipath propagation which involves more than about three metres additional path length, as will generally be the case for a distant source.

The above-mentioned US patent pays particular attention to attempting to eliminate the effects of pulse amplitude on ToA measurement by using a so-called Normalizer, but makes no reference to the possible effects of multipath. In the first of two Normalizer processes described with reference to Figure 7 of the patent, the value of the signal amplitude at which the signal is timed is dependent on the rate of increase of the amplitude, and since this value would necessarily have to be at a minimal threshold for satisfactory distinction of received signals from noise when the rate of increase is at its lowest acceptable value, the amplitude value at which the signal is timed will generally be above the minimal threshold value and hence not as early on the leading edge as it might be.In the second Normalizer process described with reference to Figure 8 of the patent, a second threshold well above a first threshold is used. The Normalizer processes are therefore inherently more susceptible to multipath which affects the shape of the leading edge of the signal. Moreover, in view ofthe relatively slow-acting circuitry described in the patent, the Normalizers would necessarily require the timed point on a pulse to be well after the start of the pulse. The present invention involves the recognition that in practice, it is more desirable to obviate the effect of most multipath propagation: timing errors due to multipath will generally be worse than errors due to different signal amplitudes. Moreover, by making the ToA measurement early on the leading edge, inaccuracies due to different signal amplitudes may in any case be reduced.

According to the third aspect of the invention, a timing circuit comprising a clock pulse generator, a tapped delay device having a plurality of n mutually spaced taps, a latch coupled to the delay device for latching any signal on each of the n taps, and a decoding device coupled to the latch for producing a time representation from the signal(s) latched from the n taps, is characterised in that an input signal to be timed is coupled to the input of the delay device, in that the clock pulse generator is normally operable to clock the latch, in that the circuit comprises inhibiting means responsive to the presence of a signal on at least one of then taps when the latch is clocked to inhibit further clocking of the latch, and in that the decoding device is operable to produce a representation of the interval between the time that said input signal reaches the tap nearest the input of the delay device and the preceding clock pulse.

Said interval may be represented as zero for the case in which the input signal has reached the tap furthest from the input when the latch is clocked, and other intervals represented accordingly.

Atiming circuit embodying the third aspect of the invention may be contrasted with the timing circuit disclosed in GB 2 132 043 A and EP 113 935 A, in which the clock pulse generator is coupled to the input of the delay line, and the input signal to be timed is used to latch the latch.

Suitably, the circuit further comprises a counter for counting the pulses of the clock pulse generator, wherein the inhibiting means are further operable to inhibit further counting of the clock pulses, the outputs of the decoding device and the counter being concatenated.

In order to be able to produce representations of intervals over the majority of the period of the clock pulse generator, the period of the clock pulse generator may be not substantially less than the time delay between the tap nearest to and the tap furthest from the input of the delay device. Suitably, said period is substantially equal to said time delay.

To make good use of the delay device and to provide representations of integral multiples of a fraction of the period of the clock pulse generator, the time delay between each adjacent pair of the n taps may be the same, being equal to T, and the period of the clock pulse generator be nT.

Where the circuit is to be used to time the beginning of signals which persist for at least the period of the clock pulse generator, the inhibiting means may be responsive to the presence of a signal on the tap nearest the input of the delay device when the latch is clocked. This helps to distinguish true signals on a tap beyond the tap nearest the input of the delay device when the latch is clocked, since such a signal might be due to noise. As a further safeguard against false measurements due for example to noise, the decoding device may be operable not to produce said time representation unless when the latch is clocked a signal is present on each of the n taps between the input of the delay device and the tap furthest from the input of the delay device on which a signal is present.

It has been found that attempting to operate a timing circuit of the kind disclosed in the above-mentioned GB and EP published applications to measure time to a resolution of about 1 nanosecond produces difficulties in synchronising the fine count produced by the decoding device and the coarse count produced by the counter. A timing circuit embodying the third aspect of the present invention can be both simpler and more reliable. Furthermore, it has been found advantageous to use the threshold crossing merely to feed an input signal to the delay device, rather than to use it to control gates: the latter is liable in practice to produce undesired distortion of the signal.

A method embodying the first aspect of the invention may involve using three substantially coplanar but substantially non-collinear elements to form at least two said pairs, and deriving a representation of the angle (3 and/or a representation of the angle ss utilising the relationship sin (90-ex)=sin ocos ss in respect of each of said at least two pairs, wherein 0 is the angle between the direction of incidence projected into the plane of the three elements and the normal in said plane to the line joining the two elements of a said pair, and ss is the angle between the direction of incidence and said plane. Directions of incidence in three dimensions may thus be determined from ToA measurements at three elements.

The method may further cornprise determining a parameter representative of the rate of increase of the detected amplitude of the respective signal received at at least one of the elements in the region of said threshold value, and determining the direction of incidence of received signals only if said parameter satifies a criterion representing a minimum rate of increase in said region. Signals whose direction of incidence cannot be determined with satisfactory accuracy can thereby be discarded. Such a method may comprise measuring the time at which the detected amplitude first exceeds an adjacent further threshold value, wherein said parameter is the difference between the measured times in respect of the two threshold values, and wherein said criterion is that said parameter does not exceed a maximum value.Alternatively, such a method may comprise differentiating the increasing detected amplitude at least in said region, wherein said parameter is the rate of increase in detected amplitude derived by differentiation, and wherein said criterion is that said parameter exceeds a minimum value.

A system embodying the second aspect of the invention may comprise means for performing optional features of a method embodying the first aspect of the invention as set forth in Claims 7 to 10.

Embodiments of the invention will now be described, by way of example, with reference to the diagrammatic drawings, in which: Figure lisa block diagram of a ToA DF system embodying the invention and comprising one pair of wave-receiving elements; Figure 2 is a block diagram of a timing circuit embodying the invention suitable for use in a ToA DF system embodying the invention; Figure 3 is a block diagram of a ToA DF system embodying the invention and comprising three collinear wave-receiving elements; Figure 4 illustrates the disposition of three non-collinear wave-receiving elements for an omnidirectional ToA DF system embodying the invention, and Figure 5 illustrates schematically processing to calculate an angle using the elements of Figure 4.

Figure 1 shows a simple system comprising one pair of wave-receiving elements. The system comprises two similar channels A and B respectively. Each of the channels comprises in succession an antenna ANT, an RF amplifier RFA, a detector D, a video amplifier VA, and a timing circuit TC. The antennae may be omnidirectional, or may be directional with their axes substantially parallel. The spacing d between the antennae is chosen to be sufficiently small that the length of the interval of time within which signals from the same source must arrive at the two elements is so short that there is a high probability in operation that no signals from another source will arrive in that interval.If signals may be received from any direction, the length of the interval is twice the time taken for electromagnetic waves to travel the distance d (in free space): the limits of this interval are set by the possibility of signals being incident along the line joining the antennae in one sense or the other, i.e. from left or right in Figure 1. The probability of no signals arriving from another source in that interval will depend on the number of sources from which signals can be received, the frequency with which they emit fresh signals, and the duration of the signals.What probability is sufficiently high will depend on what proportion of incorrect representations of direction, due to measurements having been made in the two channels on signals from different sources, is considered acceptabie. For typical operational situations of the number of sources emitting pulsed signals, their PRF and pulse length, a separation d of the order of 50 feet, giving an interval length of about 100 nanoseconds, is considered to give an acceptably high probability. When an RF signal is incident on the antenna of one of the channels, the RF signal is amplified and detected, and the time at which the instantaneous detected amplitude of the signal, after further amplification, first exceeds a threshold value is measured.This threshold value is chosen to be substantially less than the minimum peak value of signals whose direction of incidence is to be determined, as will be explained in greater detail below; the threshold will usually be much less than the typical peak value.

The timing circuits operate with a common clock (CLOCK). A caiculating unit CALC determines the difference 8t between the measured times and provides therefrom a representation of the direction of incidence, for example of the angle a between the incident signals and the line joining the antennae, utilising the relationship cos a=c8Vd where c is the free-space velocity of electromagnetic waves.

The time at which the threshold value is first exceeded is measured in each timing circuit TC in such a manner that the measured time is unaffected by multipath propagation where the delayed signal has been reflected from a surface not virtually coincident with the direct path between a distant source and the antenna.

An ECL (emitter-coupled logic) circuit arrangement, suitable for the timing circuit TC of Figure 1 and capable of measuring to an accuracy of a nanosecond, is shown in Figure 2. The detected and amplified video signal from the video amplifier VA is applied to a very fast comparator COMP whose output changes from a logic "0" to "1" when the input signal exceeds a threshold voltage VT. The comparator output signal is fed to the input of a tapped delay line TDL having 8 consecutive taps separated by 1 nanosecond intervals. The taps are connected to respective inputs of an 8-bit latch LATCH. The latch is clocked at 8 nanosecond intervals by a 125 MHz clock via a gate G1,the clock signal also being supplied via a further gate G2 to a synchronous counter CNTR which provides a coarse measurement of time.The outputs, labelled 0-7, of the latch are fed to a decoding circuit DCDR; the presence of a signal on at least one output, in this case the first output, 0, is also used to control the gates G1 and G2, the output being connected thereto by a fast feedback loop. The outputs of the counter CNTR and the decoder DCDR are concatenated to give a representation of the time at which the output of the comparator COMP changed from 0 to 1.

In operation, the gates G1 and G2 are normally open. The counter CNTR measures time in 8-nanosecond units, up to a maximum time at least as long as the above-mentioned interval, determined by the antenna spacing d, within which signals from a distant source must be received by both antennae. The latch is similarly clocked at 8-nanosecond intervals, but while the comparator output is 0, there are no signals from the tapped delay line, and the latch outputs remain at zero. When the comparator output changes to 1 (which in the case of a true received signal as opposed to noise will normally persist for longer than the period of the clock), the signal travels along the delay line changing successive tap outputs from 0 to 1. When the latch is next clocked, the series of ones and remaining zeroes is held in the latch. The presence of a "1" on the first output, 0, of the latch closes the gates G1 and G2, preventing further clocking of the latch and the counter. The output of the latch will be one of the following codes: 10000000 11000000 11100000 11110000 11111000 11111100 11111110 11111111 The first of these codes represents the most recent ToA and the last the earliest ToA since the latch was last clocked; the first indicates that 7 nanoseconds should be added to the time represented by the counter CNTR while the last requires zero addition. The decoder DCDR transforms that latch output to binary digits which are concatenated with the counter output.

The threshold the crossing of which is timed is, as previously explained, set at a low value in order largely to avoid errors due to multipath. Setting the threshold at a level substantially below the minimum peak level of signals whose direction of incidence is to be determined, for example 10 dB below the minimum peak level, also provides the advantage of tending to alleviate timing errors which would occur if the signal amplitude crossed the threshold at a slow rate because the amplitude were reaching its peak level. The lowest level at which the threshold can be set will depend on the noise level in the system: if the threshold is set too close to the noise level, the accuracy of timing will be degraded by the random fluctuation in the amplitude of desired signal plus noise, and a positive-going threshold crossing may even be caused by noise alone.

The decoder is in this embodiment arranged to accept only the above-mentioned codes. It consequently accepts only signals whose amplitude remains above threshold long enough to produce a continuous succession of ones in the latch, and rejects any other pattern of zeroes and ones which might result from triggering of the comparator by noise spikes or from a received signal with a slow rate of increase of amplitude.

It will be seen that provided the amplitude remains above threshold long enough for this state to be latched, the measurement of ToA will be unaffected by subsequent variations in amplitude, in particular variations due to a delayed multipath signal which arrives in phase opposition to the original direct-path signal and causes the amplitude to fall below threshold. The maximum period taken to latch the above-threshold state is in this embodiment the length of the delay line, i.e. 8 nanoseconds. This time may be much shorter than the time taken for the amplitude to reach typical peak value.

A delay not less than and substantially equal to the time taken for electromagnetic waves to travel the distanced may be included in one channel before the timing circuit so that time differences are measured with respect to the time of arrival of a signal in the other channel. Suitably, a time "window" is used to prevent unnecessary computation on time differences which are too large for the signals to have come from the same source. Where the above-mentioned delay is included in one channel, this window may be defined as beginning with a time difference of zero and ending with a time difference not less than and substantially equal to twice the time taken for electromagnetic waves to travel the distance d. The use of a time window also provides some protection against random noise signals which cause the detected amplitude to exceed the threshold from causing false measurements.

At least one, and preferably each, of the channels in the system of Figure 1 may comprise a signal validating circuit to ascertain the rate of increase of the amplitude of the signal in that channel in the region of the threshold value, and to cause the system not to determine the direction of incidence unless the rate satisfies a criterion of minimum slope. For this purpose, the output of the video amplifier VA may be supplied to a further timing circuit which measures the time at which the signal amplitude first exceeds an adjacent further threshold value. The difference between the times measured by the two circuits in a channel may be determined, and the direction of incidence determined only if the difference does not exceed a maximum value.

Alternatively, the amplitude increase may be differentiated and the direction of incidence determined only if the rate of increase of amplitude derived by differentiation exceeds a minimum value.

As a further way of distinguishing signals coming from a distant source from noise, at least one of the channels may compries means for determining the peak amplitude of a signal which causes the threshold to be exceeded, and for inhibiting the determination of the direction of incidence unless the peak amplitude is substantially greater than the threshold.

It is considered that a suitable criterion of minimum slope may be that the steepness of the rising edge should be predominantly controlled by the video bandwidth of the system. This inter alia has the effect of reducing the dependence of the measured time on the rate of increase of the RF amplitude and hence tends to achieve the same object as the Normalizers in the above-mentioned US patent. It may be desirable for the video bandwidth to be switchable between a broadband value and a narrowband value. The broadband value may allow more accurate timing of steeply-rising leading edges, but the narrowband value may enable acceptable results to be obtained with more slowly-rising edges, since it may reduce the noise level in the video circuit and allow the threshold to be set to a lower value and hence to a relatively steeper part of the leading edge.

Figure 3 depicts a modification of the system of Figure 1 comprising three coplanar and collinearly disposed antennae, L, M, N respectively in respective channels each as each of the channels shown in Figure 1. The spacing of each of the two pairs of elements LM and MN are equal (each being d) and each satisfies the above-mentioned criterion that the spacing is sufficiently small that the length of the interval of time within which signals from the same source must arrive at the two elements of a pair is so short that there is a high probability in operation that no signals will arrive from another source in that interval; the spacing 2d between antennae Land N may however be too large to satisfy this criterion. Nevertheless, the difference between the times of arrival of signals at antennae Land N may be used to provide a more accurate representation of the direction of incidence than could be provided by the system of Figure 1 if one or more steps are taken to reduce the possibility that the time difference measured on one of the pairs of antennae LM, MN does not relate to the same source as the time difference measured on the other pair.For example, as indicated in Figure 3, the time differences measured in relation of each pair of antennae, tLM and tMN respectively, may be compared, and only if their values are equal to within a small tolerance is the direction of incidence determined from the difference between the times of arrival at antennae Land N; the probability that signals from different sources should result in substantially equal time differences being measured between the elements of the pairs LM and MN is small, and even if the signals should have come from different sources, the resultant error in the indicated direction of incidence in relation to the source from which signals were first received at one pair of antennae will be small.

The systems so far described provide only a representation of the direction of incidence that defines an angle to the line joining a pair of wave-receiving elements, and hence the surface of a cone whose axis is said line. Whose sources are known to lie substantially in a single plane including said line and where the wave-receiving elements are directional, this may be sufficient (although it should be borne in mind that the accuracy with which a can be determined decreases as a decreases from 90 degrees to 0); however, when signals may be received from each side of said line, and particularly when sources are not restricted to a single plane, it is desirable to perform measurements on at least one further wave-receiving element which is not collinearwith the one pair of elements.Figure 4 depicts the deposition of three substantially coplanar but not collinear elements A, B, C respectively, forming an arbitrarily-shaped triangle. The spacings AB, BC, CA each satisfy the above-mentioned criterion of being sufficiently small. By measuring the times of arrival of signals at each element, the direction of incidence may be determined for the general case of distant sources in 3-dimensional space as follows.

Let the spacings AB, BC, CA be d, d2, d3 respectively. Let the angles CAB and ABC be m and n respectively.

Let the length of the perpendicular from C onto AB be a, and the distance from A to the intersection of said perpendicular with AB be b, so that a=d3 sin m and b=d3 cos m. (Thus b is negative if m > 90 degrees.) Let the angle between the direction of incidence and the normal to AB in the plane of AB and the direction of incidence be oy (so that pry=(90 degrees) and cos a=sin y), the angle between the normal to AB in the plane ABC and the direction of incidence projected into that plane be (3 (typically the azimuth angle), and the angle between the direction of incidence and the plane of ABC be ss (typically the elevation angle).Let the times of arrival at A, B, C be tA, ts, tc respectively.

Then sin If = c(tA-t3)/dl Now sin ag = sin a cos ft Writing x = c(tA-tB)/d1 y = c(tA-tC)/d3 z = c(tB-tC)/d2, we may put x = sin 8 cos A and analogously y = sin (B-m) cos z = sin 8-(180-n)7 cosp = -sin (8+n) cos A Eliminating ss from either of two pairs of these equations, one obtains tan 8 = x sin miEx cos m and tan 8 = -x sin nIEx cos n + z.

These expressions have two-foid ambiguity. To distinguish between -90 degrees < o < 90 degrees and 90 degrees < 0 < 270 degrees, one may note that the denominator of for example the first expression for tan 0 may be expanded as sin 8 cos A cos m - sin(8-m) cos or cos A cos 8 sin m.

Thus the denominator is positive for -90 degreees < 0 < 90 degrees and negative for 90 degrees < 0 < 270 degrees.

Re-writing the first expression for tan o in terms of the times of arrival and multiplying the numerator and denominator by dad3, tan g = (tA-tB)d3 sin m/C(tA-tg)cf3 cos m -(tA-tC)d13 = - a (tA-tB)/Edl (tA-tc) - b(tA-tB)3 = - a (tA-tg)/E(d1-b)(tA-tC)+b(tB-tc)] Writing P = a(tA-tB) Q = (d1-b) (tA-tc) R = b (tB-tC) one obtains tan B = - P/(Q+R) or 9 = - Arc tan CP/(Q+R)].

Figure 5 illustrates schematically the processing to calculate 0 according to this expression. The differences tA-tB, ts-tc, tA-tC are formed from the measured times tA, ts, tc and then scaled to produce the quantities P, Q, R. The quantities Ci and R are summed and divided into P; the angle whose tangent is equal to minus the quotient is then determined, for example from a look-up table in a PROM, to obtain an ambiguous value Oamb of 6. The sign of (Q+Rì is also checked; if (Q+R) > 6, a quantity 6o= 180 degrees is produced, otherwise 6=00.

6amb and 00 are summed to produce an unambiguous value of 6. Having calculated 6, the angle ss may be calculated by substituting 6 in, for example, the expression x = sin a cos A Alternatively, p may be calculated without needing to calculate 6 by eliminating 6 from a pair of the expressions for x, y, z.

The calculations and processing may be simplified for particular cases. For an equilateral triangle of side d, a = J3d/2 and b = d/2, so that the time differences need only be scaled by factors which are independent of d. Alternatively, if m=90 degrees, a = d3 and b = 0, so that the quantity R is zero. If d1 =d3, the scaling factors are again independent of the actual value of the spacing.

An omnidirectional direction-finding system may comprise four receiving elements disposed at the corners of a parallelogram, or more especially a rectangle and more particularly still a square.

The direction of incidence may be calculated from the times of arrival of the earliest-received three signals which are of acceptable quality. This allows for the possibility that signals received at one of the four elements may have been degraded by, for example, an obstruction in the region of the elements.

Comparison of the equations which can be derived from the above two expressions for tan 6 in terms of x andy and ofx and z respectively to relate the error in 6 to errors in x, y and z show that the error is not dependent on which expression is used. The choice of which baselines are considered as primary and secondary baselines for determining 6 is therefore not significant.

As alternative to a clock common to the timing circuits of all the channels as depicted in Figures 1 and 3, each channel may have a respective accurate clock and the clocks be kept in synchronism via a low-bandwidth link. The time measurements and any other other data may be passed to a central processor and control unit via, for example, an optical fibre link.

Claims (20)

1. A method of determining the direction of incidence of electromagnetic wave signals frqm a distance source from the time of arrival of the leading edge of the wave signals, wherein the method comprises: receiving said signals at a plurality of mutually spaced wave-receiving elements, detecting the respective instantaneous amplitude of the signal received at each element, measuring the times at which the detected amplitudes of wave signals received respectively at at least two of said elements first exceed a minimal threshold value such that signals can be satisfactorily distinguished from noise and which is substantially less than the minimum peak value of signals whose direction of incidence is to be determined by said method, the time being measured in such a manner that the measured time is generally unaffected by multipath propagation, determining the difference between the measured times in respect of one pair or of a plurality of pairs of said elements, wherein the two elements of said one pair or of each of at least two of said plurality of pairs are sufficiently close together that the length of the interval of time within which signals from the same source must arrive at the two elements is so short that there is a high probability in operation that no signals from another source will arrive in that interval, and deriving a representation of the direction of incidence from the time difference(s) utilising the relationship cos a=cbt/d where a is the angle between the direction of incidence of the signals and the line joining the two elements of a said pair, d is the distance between those two elements, 8t is the time difference in respect of that pair of elements, and c is the free-space velocity of electromagnetic waves.

2. A method as claimed in Claim 1 using three substantially coplanar but substantially non-collinear elements to form at least two said pairs, the method comprising deriving a representation of the angle 0 and/or a representation of the angle ss utilising the relationship sin (90-a)=sin n cos ss in respect of each of said at least two pairs, wherein o is the angle between the direction of incidence projected into the plane of the three elements and the normal in said plane to the line joining the two elements of a said pair, and ss is the angle between the direction of incidence and said plane.

3. A method as claimed in Claim 1 or 2 which further comprises determining a parameter representative of the rate of increase of the detected amplitude of the respective signal received at at least one of the elements in the region of said threshold value, and determining the direction of incidence of received signals only if said parameter satisfies a criterion representing a minimum rate of increase in said region.

4. A method as claimed in Claim 3 which comprises measuring the time at which the detected amplitude first exceeds an adjacent further threshold value, wherein said parameter is the difference between the measured times in respect of the two threshold values, and wherein said criterion is that said parameter does not exceed a maximum value.

5. A method as claimed in Claim 3 which comprises differentiating the increasing detected amplitude at least in said region, wherein said parameteris the rate of increase in detected amplitude derived by differentiation, and wherein said criterion is that said parameter exceeds a minimum value.

6. A system for determining the direction of incidence of electromagnetic wave signals from a distant source from the time of arrival of the leading edge of the wave signals, wherein the system comprises: a plurality of mutually spaced wave-receiving elements, means for detecting the respective instantaneous amplitude of the signal received at each element, means for measuring the times at which the detected amplitudes of wave signals received respectively at at least two of said elements first exceed a minimal threshold value such that signals can be satisfactorily distinguished from noise and which is substantially less than the minimum peak value of signals whose direction of incidence is to be determined by said method, the time being measured in such a manner that the measured time is generally unaffected by multipath propagation, means for determining the difference between the measured times in respect of one pair or of a plurality of pairs of said elements, wherein the two elements of said one pair or of each of at least two of said plurality of pairs are sufficiently close together that the length of the interval of time within which signals from the same source must arrive at the two elements is so short that there is a high probability in operation that no signals from another source will arrive in that interval, and means for deriving a representation of the direction of incidence from the time difference(s) utilising the relationship cos a=c3t/d where a is the angle between the direction of incidence of the signals and the line joining the two elements of a said pair, d is the distance between those two elements, 3t is the time difference in respect of that pair of elements, and c is the free-space velocity of electromagnetic waves.

7. A system as claimed in Claim 6 comprising three substantially coplanar but substantially non-collinear elements disposed to form at least two said pairs, wherein said means for deriving a representation of the direction of incidence comprises means for deriving a representation of the angle 0 and/or a representation of the angle ss utilising the relationship sin (90-a)=sin 0cos ss in respect of each of said at least two pairs, wherein 0 is the angle between the direction of incidence projected into the plane of the three elements and the normal in said plane to the line joining the two elements of a said pair, and p is the angle between the direction of incidence and said plane.

8. A system as claimed in Claim 6 or 7 comprising means for determining a parameter representative of the rate of increase of the detected amplitude of the respective signal received at at least one of the elements in the region of said threshold value, and means for inhibiting the determination of the direction of incidence if said parameter does not satisfy a criterion representing a minimum rate of increase in said region.

9. A system as claimed in Claim 8 wherein the parameter-determining means comprise means for measuring the time at which the detected amplitude first exceeds an adjacent further threshold value, said parameter being the difference between the measured times in respect of the two threshold values, and said criterion being that said parameter does not exceed a maximum value.

10. A system as claimed in Claim 8 wherein the parameter-determining means comprise means for differentiating the increasing detected amplitude at least in said region, said parameter being the rate of increase in detected amplitude derived by differentiation, and said criterion being that said parameter exceeds a minimum value.

11. A timing circuit comprising a clock pulse generator, a tapped delay device having a plurality of n mutually spaced taps, a latch coupled to the delay device for latching any signal on each of the n taps, and a decoding device coupled to the latch for producing a time representation from the signal(s) latched from the n taps, characterised in that an input signal to be timed is coupled to the input of the delay device, in that the clock pulse generator is normally operable to clock the latch, in that the circuit comprises inhibiting means responsive to the presence of a signal on at least one of the n taps when the latch is clocked to inhibit further clocking of the latch, and in that the decoding device is operable to produce a representation of the interval between the time that said input signal reaches the tap nearest the input of the delay device and the preceding clock pulse.

12. A circuit as claimed in Claim 11 further comprising a counter for counting the pulses of the clock pulse generator, wherein the inhibiting means are further operable to inhibit further counting of the clock pulses, the outputs of the decoding device and the counter being concatenated.

13. A circuit as claimed in Claim 11 or 12 wherein the period of the clock pulse generator is not substantially less than the time delay between the tap nearest to and the tap furthest from the input of the delay device.

14. A circuit as claimed in Claim 13 wherein said period is substantially equal to said time delay.

15. A circuit as claimed in Claim 14 wherein the time delay between each adjacent pair of then taps is the same, being equal to T, and the period of the clock pulse generator is nT.

16. A circuit as claimed in any of Claims 11 to 15 wherein the inhibiting means are responsive to the presence of a signal on the tap nearest the input of the delay device when the latch is clocked.

17. A circuit as claimed in any of Claims 11 to 16 wherein the decoding device is operable not to produce said time representation unless when the latch is clocked a signal is present on each of the n taps between the input of the delay device and the tap furthest from the input of the delay device on which a signal is present.

18. A method of radio direction-finding, substantially as herein described with reference to Figures 1 and 3, to Figures 1,4 and 5, to Figures 1 and 2, to Figures 1,2 and 3, or to all the Figures of the drawings.

19. A system for radio direction-finding, substantially as herein described with reference to Figure 1, to Figures 1 and 3, to Figures 1,4 and 5, to Figures 1 and 2, to Figures 1,2 and 3, or to all the Figures of the drawings.

20. Atiming circuit substantially as herein described with reference to Figure 2.