System and method for determining the time difference of arrival of a frequency shift keyed signal at two separate receivers

• US 6,137,842 A
• Filed: 11/28/1997
• Issued: 10/24/2000
• Est. Priority Date: 11/28/1997
• Status: Expired due to Fees

• Current Assignee: United States Department Of Health And Human Services (Government of the United States of America)
• Original Assignee: the united states of america as represented by the secretary of the navy
• Inventors: Grossnickle, Peter C.
• Primary Examiner(s): Chin, Stephen
• Assistant Examiner(s): Fan, Chieh M.
• Application Number:
  US08/979,633
• Time in Patent Office: 1,061 Days
• Field of Search: 375/224, 375/225, 375/324, 375/334, 375/347, 375/349, 375/354, 375/272, 342/389, 342/442, 342/450, 342/451, 455/415, 455/456
• US Class Current: 375/334
• CPC Class Codes: G01S 5/06 Position of source determined by co-ordinating a plurality of position lines defined by path-difference measurements G01S5/12 takes precedence

Abstract

A system determines the time difference of arrival of a frequency shift kd signal (FSK) source at two separate receiver sites. The system includes a first receiver located at a first site for generating a first series of signal samples in response to detecting an FSK signal from a FSK signal source. A second receiver located at a second site generates a second series of signal samples in response to detecting the FSK signal. A first computer estimates first frequency shift times of the FSK signal in response to receiving the first series of signal samples and a time reference signal, and generates a first output signal representing the first frequency shift times. A second computer estimates second frequency shift times of the FSK signal using the second series of signal samples and time reference signal, and determines the path difference between the FSK signal source and each of the first and second receivers using the first output signal containing the first frequency shift times, and the second frequency shift times. The invention obviates the need for providing a replica of the first signal samples to the second computer.

Claims

1. A system for determining the time difference of arrival of a signal between two separate locations, comprising:

a first receiver located at a first site for generating first signal samples representing an FSK signal in response to receiving the FSK signal at time T_A ;

a second receiver located at a second site for generating second signal samples representing the FSK signal in response to receiving the FSK signal at time T_B, where T_A ≠

T_B ;

a first computer for estimating a first set of frequency shift times of the FSK signal in response to receiving the first signal samples and a time reference signal, and for generating a first output signal representing the first set of frequency shift times; and

a second computer for estimating a second set of frequency shift times of the FSK signal using the second signal samples and time reference signal, and for determining the time difference of arrival of the FSK signal between the first and second receivers using the first output signal, and the second set of frequency shift times in accordance with the relation;

space="preserve" listing-type="equation">TDOA={[t.sub.2 (1)-t.sub.1 (1)]+[t.sub.2 (2)-t.sub.1 (2)]+[t.sub.2 (3)-t.sub.1 (3)] . . . +[t.sub.2 (K)-t.sub.1 (K)]}/(K),

where {t₁ (1), t₁ (2), t₁ (3), . . . t₁ (K)} represent a K number of first frequency shift times within an N number of the first signal samples, {t₂ (1), t₂ (2), t₂ (3), . . . t₂ (K)} represent a K number of second frequency shifts times within an N number of the second signal samples, and where K and N are positive integers.

2. The system of claim 1 wherein said first and second receivers each are digital receivers.

3. A method for determining the time difference of arrival of a signal between two separate locations, comprising the steps of:

detecting an FSK signal with a first receiver at a first location;

generating first signal samples representing the FSK signal;

detecting the FSK signal with a second receiver at a second location;

generating second signal samples representing the FSK signal;

estimating a first set of frequency shift times of the FSK signal from the first signal samples and a time reference signal;

generating a first output signal representing the first set of frequency shift times;

estimating a second set of frequency shift times of the FSK signal from the second series of signal samples and time reference signal; and

determining the time difference of arrival of the FSK signal at the first and second receivers using the first output signal, and the second set of frequency shift times, wherein the time difference of arrival is determined in accordance with the relation;

space="preserve" listing-type="equation">TDOA={[t.sub.2 (1)-t.sub.1 (1)]+[t.sub.2 (2)-t.sub.1 (2)]+[t.sub.2 (3)-t.sub.1 (3)] . . . +[t.sub.2 (K)-t.sub.1 (K)]}/K,

where {t₁ (1), t₁ (2), t₁ (3), . . . t₁ (K)} represent a K number of the first frequency shift times within an N number of the first signal samples, {t₂ (1), t₂ (2), t₂ (3), . . . t₂ (K)} represent a K number of the second frequency shifts times within an N number of the second signal samples, and K and N are positive integers.

4. A system for determining the time difference of arrival of a signal at two separate locations, comprising:

a first receiver located at a first site for generating first signal samples representing a signal that transitions between first and second states;

a second receiver located at a second site for generating second signal samples representing the signal;

a first computer for estimating a first set of transition times within the signal in response to receiving the first signal samples and a time reference signal, and for generating a first output signal representing the first set of transition times, where the first transition times occur when the signal transitions between the first and second states; and

a second computer for estimating a second set of transition times within the signal using the second series of signal samples and time reference signal, and for determining the time difference of arrival of the signal at the first and second receivers using the first output signal, and the second set of transition times, wherein the time difference of arrival is determined in accordance with the relation;

space="preserve" listing-type="equation">TDOA={[t.sub.2 (1)-t.sub.1 (1)]+[t.sub.2 (2)-t.sub.1 (2)]+[t.sub.2 (3)-t.sub.1 (3)] . . . +[t.sub.2 (K)-t.sub.1 (K)]}/K,

where {t₁ (1), t₁ (2), t₁ (3), . . . t₁ (K)} represent a K number of first frequency shift times within an N number of the first signal samples, {t₂ (1), t₂ (2), t₂ (3), . . . t₂ (K)} represent a K number of second frequency shifts times within an N number of the second signal samples, and K and N are positive integers.

5. The system of claim 4 wherein the communications signal is a phase shift keyed signal.

6. The system of claim 5 wherein the communications signal is an amplitude shift keyed signal.

7. A method for determining the time difference of arrival of a signal at two separate locations, comprising the steps of:

detecting a signal that transitions between first and second states with a first receiver at a first location;

generating first signal samples representing the signal;

detecting the signal with a second receiver at a second location;

generating second signal samples representing the signal;

estimating a first set of transition times of the signal from the first signal samples and a time reference signal, where the first transition times occur when said signal transitions between the first and second states;

generating a first output signal representing the first transition times;

estimating a second set of transition times of the signal from the second signal samples and time reference signal, where the second transition times occur when said signal transitions between the first and second states; and

determining the time difference of arrival of the signal at the first and second receivers using the first output signal and the second transition times, wherein said time difference of arrival is determined in accordance with the relation;

space="preserve" listing-type="equation">TDOA={[t.sub.2 (1)-t.sub.1 (1)]+[t.sub.2 (2)-t.sub.1 (2)]+[t.sub.2 (3)-t.sub.1 (3)] . . . +[t.sub.2 (K)-t.sub.1 (K)]}/K,

where {t₁ (1), t₁ (2), t₁ (3), . . . t₁ (K)} represent a K number of first frequency shift times within an N number of the first signal samples, {t₂ (1), t₂ (2), t₂ (3), . . . t₂ (K)} represent a K number of second frequency shifts times within an N number of the second signal samples, and K and N are positive integers.

8. A system for determining the time difference of arrival of a signal between two separate locations, comprising:

a first receiver for receiving a signal that transitions between two different states and for transforming said signal into an n number of pairs of first in-phase and quadrature signal samples, where n is a positive integer index;

a first computer for determining;

a) phase values P₁ (n) for each of said first in-phase and quadrature signal samples; and

b) transition times t₁ (1), t₁ (2), t₁ (3), . . . t₁ (k₁) when said signal transitions between said first and second states, where k₁ represents the number of transitions detected from said n number of first in-phase and quadrature signal samples, and said transition times are determined from said phase values P₁ (n);

a second receiver for transforming said signal into an n number of pairs of second in-phase and quadrature signal samples; and

a second computer for determining;

a) phase values P₂ (n) for each of said second in-phase and quadrature signal samples;

b) transition times t₂ (1), t₂ (2), t₂ (3), . . . t₂ (k₂) when said signal transitions between said first and second states, where k₂ represents the number of transitions detected from said n number of second in-phase and quadrature signal samples, k₂ =k₁, and said transition times are determined from said phase values P₂ (n); and

c) a time difference of arrival TDOA of said signal between said first and second receivers, where

space="preserve" listing-type="equation">TDOA={[t.sub.2 (1)-t.sub.1 (1)]+[t.sub.2 (2)-t.sub.1 (2)]+[t.sub.2 (3)-t.sub.1 (3)] . . . +[t.sub.2 (k.sub.1)-t.sub.1 (k.sub.1)]}/k.sub.1.

9. The system of claim 8 wherein:

said first computer is programmed to;

transform said phase values P₁ (n) into first unwrapped phase samples;

group said first unwrapped phase samples into first group sets;

perform a least squared fit on each said first group set to define a set of first linear functions;

define a sequence of first adjacent linear functions from said first linear functions; and

determine a set of first intersections of said first adjacent linear functions, where said first intersections represent said transition times t₁ (1), t₁ (2), t₁ (3) . . . and t₁ (k₁); and

said second computer is programmed to;

transform said phase values P₂ (n) of said into second unwrapped phase samples;

group said second unwrapped phase samples into second group sets;

perform a least squared fit on each said second group set to define a set of second linear functions;

define a sequence of second adjacent linear functions from said second linear functions; and

determine a set of second intersections of said second adjacent linear functions, where said second intersections represent said transition times t₂ (1), t₂ (2), t₂ (3) . . . and t₂ (k₂).

10. A method for determining the time difference of arrival of a signal between two separate locations, comprising the steps of:

receiving a signal at a first location where said signal transitions between two different states;

transforming said signal into an n number of pairs of first in-phase and quadrature signal samples, where n is a positive integer index;

determining phase values P₁ (n) for each of said first in-phase and quadrature signal samples;

determining transition times t₁ (1), t₁ (2), t₁ (3), . . . t₁ (k₁) when said signal transitions between said first and second states, where k₁ represents the number of transitions detected from said n number of first in-phase and quadrature signal samples, and said transition times are determined from said phase values P₁ (n);

receiving said signal at a second location;

transforming said signal into an n number of second in-phase and quadrature signal samples;

determining phase values P₂ (n) for each of said second in-phase and quadrature signal samples;

determining transition times t₂ (1), t₂ (2), t₂ (3), . . . t₂ (k₂) when said signal transitions between said first and second states, where k₂ represents the number of transitions detected from said n number of second in-phase and quadrature signal samples, k₂ =k₁, and said transition times are determined from said phase values P₂ (n); and

determining a time difference of arrival TDOA of said signal between said first and second receivers, where;

space="preserve" listing-type="equation">TDOA={[t.sub.2 (1)-t.sub.1 (1)]+[t.sub.2 (2)-t.sub.1 (2)]+[t.sub.2 (3)-t.sub.1 (3)] . . . +[t.sub.2 (k.sub.1)-t.sub.1 (k.sub.1)]}/k.sub.1.

11. The method of claim 10 further including the steps of:

transforming said phase values P₁ (n) into first unwrapped phase samples;

grouping said first unwrapped phase samples into first group sets;

performing a least squared fit on each said first group set to define a set of first linear functions;

defining a sequence of first adjacent linear functions from said first linear functions;

determining a set of first intersections of said first adjacent linear functions, where said first intersections represent said transition times t₁ (1), t₁ (2), t₁ (3) . . . and t₁ (k₁);

transforming said phase values P₂ (n) of said into second unwrapped phase samples;

grouping said second unwrapped phase samples into second group sets;

performing a least squared fit on each said second group set to define a set of second linear functions;

defining a sequence of second adjacent linear functions from said second linear functions; and

determining a set of second intersections of said second adjacent linear functions, where said second intersections represent said transition times t₂ (1), t₂ (2), t₂ (3) . . . and t₂ (k₂).

12. A method for determining the time difference of arrival of a signal between two separate locations, comprising the steps of:

receiving a signal at a first location where said signal transitions between two different states;

transforming said signal into an n number of pairs of first in-phase and quadrature signal samples, where n is a positive integer index;

determining phase values P₁ (n) for each of said first in-phase and quadrature signal samples;

determining first transition times t₁ (1), t₁ (2), t₁ (3), . . . t₁ (k₁) when said signal transitions between said first and second states, where k₁ represents the number of transitions detected from said n number of first in-phase and quadrature signal samples, and said transition times are determined from said phase values P₁ (n);

receiving said signal at a second location;

transforming said signal into an n number of second in-phase and quadrature signal samples;

determining phase values P₂ (n) for each of said second in-phase and quadrature signal samples;

determining second transition times t₂ (1), t₂ (2), t₂ (3), . . . t₂ (k₂) when said signal transitions between said first and second states, where k₂ represents the number of transitions detected from said n number of second in-phase and quadrature signal samples, k₂ =k₁, and said transition times are determined from said phase values P₂ (n); and

determining a time difference of arrival TDOA of said signal between said first and second receivers from said first transition times transition times t₁ (1), t₁ (2), t₁ (3), . . . t₁ (k₁), and second transition times t₂ (1), t₂ (2), t₂ (3), . . . t₂ (k₂).

13. The method of claim 12 wherein said time difference of arrival is determined in accordance with the following algorithm:

space="preserve" listing-type="equation">TDOA={[t.sub.2 (1)-t.sub.1 (1)]+[t.sub.2 (2)-t.sub.1 (2)]+[t.sub.2 (3)-t.sub.1 (3)] . . . +[t.sub.2 (k.sub.1)-t.sub.1 (k.sub.1)]}/k.sub.1.

14.

14. The method of claim 12 further including the steps of:

transforming said phase values P₁ (n) into first unwrapped phase samples;

grouping said first unwrapped phase samples into first group sets;

performing a least squared fit on each said first group set to define a set of first linear functions;

defining a sequence of first adjacent linear functions from said first linear functions;

determining a set of first intersections of said first adjacent linear functions, where said first intersections represent said first transition times t₁ (1), t₁ (2), t₁ (3) . . . and t₁ (k₁);

transforming said phase values P₂ (n) of said into second unwrapped phase samples;

grouping said second unwrapped phase samples into second group sets;

performing a least squared fit on each said second group set to define a set of second linear functions;

defining a sequence of second adjacent linear functions from said second linear functions; and

determining a set of second intersections of said second adjacent linear functions, wherein said second intersections represent said second transition times t₂ (1), t₂ (2), t₂ (3) . . . and t₂ (k₂).