The present invention relates, in general, to a method and apparatus for an instantaneous frequency measurement (IFM) receiver system. More specifically, the present invention relates to a method and apparatus for measuring instantaneous frequency of an input signal by directly digitizing the input signal, without using analog delay lines or video detectors.
Electronic surveillance measurement (ESM) receivers commonly require that frequency calculations be performed on a received signal from targets of interest. The frequency of an input signal is often measured using an instantaneous frequency measurement (IFM) receiver. The IFM receiver, generally, uses a difference in phase between a delayed and a non-delayed version of the input signal to calculate the frequency of the input signal.
Two conventional IFM receivers are shown in
As also shown in
Another IFM receiver, generally designated as 17, is shown in
Conventional IFM receivers use analog components, such as hybrids, power dividers and crystal video detectors, as illustrated in
System 10 and system 17 provide different approaches to determining the input frequency of a received signal. In system 10, three 90° hybrids and one 180° hybrid are used. In system 17, on the other hand, two 90° hybrids are shown replaced by two power dividers. In both system 10 and system 17, the input signal is divided into two paths, one path is delayed by a known delay time τ, through an analog delay line, as shown. In both systems, the video signals are digitized and correlated. In building an IFM receiver, multiple correlators with different delay line lengths are typically needed. For example, some receivers may use four correlators and other receivers may use up to seven correlators.
It would be advantageous, if the number of delay lines required could be reduced. It would also be advantageous if the number of hybrids could be reduced. Furthermore, it would be advantageous if the number of crystal video detectors could be reduced.
The present invention provides such advantages by having a reduced number of components. As will be explained, analog delay lines are not necessary for the present invention. As will also be explained, the present invention only requires one 90° hybrid, and does not require any crystal video detectors.
To meet this and other needs, and in view of its purposes, the present invention provides an instantaneous frequency measurement (IFM) receiver including a hybrid for outputting in-phase (I) and quadrature (Q) signals from a received input signal; a first analog-to-digital converter (ADC) for digitizing the I signal to produce a first digital signal at a sampling rate of 1/τ, where τ is a sample time; a second ADC for digitizing the Q signal to produce a second digital signal at the same sampling rate of 1/τ; and a processor configured to delay the first and second digital signals by one sample time of τ, and count number of samples produced having predetermined phase shifts. The processor is also configured to determine a frequency of the received input signal, based on the number of samples having the predetermined phase shifts. The first and second ADCs are each 1-bit analog-to-digital converters. A sampling clock is coupled to each ADC for sampling each ADC at a rate of 1/τ. The processor is configured to count the number of samples having the predetermined phase shifts of 0°, 90°, −90° and 180° over a predetermined time interval.
In another aspect of the invention, the processor is configured to delay the first and second digital signals by at least one of 2τ, 4τ, 8τ, 16τ, 32τ, or 64τ, in addition to the 1τ sample time. The processor is also configured to count the number of samples produced having predetermined phase shifts at the at least one of 2τ, 4τ, 8τ, 16τ, 32τ, or 64τ, in addition to the 1τ sample time.
Yet another aspect of the invention includes a look up table (LUT) for access by the processor, where the LUT includes values of numbers of samples at predetermined phase shifts of 0°, 90°, −90° and 180° as a function of input frequency for determining the frequency of the received input signal. The LUT may also include values of numbers of samples of n23 and n14 as a function of frequency having the following relationship:
n23=n2−n3
n14=n1−n4
where
n1 is the predetermined phase shift of 0°,
n2 is the predetermined phase shift of 90°,
n3 is the predetermined phase shift of −90°, and
n4 is the predetermined phase shift of 180°.
Another embodiment of the present invention includes a method of determining the frequency of a signal having the following steps: (a) receiving a signal S; (b) sampling the signal S to form a digital signal C; (c) delaying C from itself by one sample time τ; (d) counting numbers of samples having predetermined phase shifts, after delaying C; and (e) determining the frequency of the signal S, based on the number of samples counted in step (d).
An aspect of the method includes the steps of: (f) forming in-phase (I) and quadrature (Q) signals from the signal S in step (a); (g) sampling the I and Q signals to form first and second digital signals, respectively; and (h) combining the first and second digital signals to form the combined signal C in step (b). Step (g) of the method includes (1) sampling the I and Q signals using respective one-bit analog-to-digital converters; and (2) sampling the I and Q signals at a sampling rate of 1/τ.
Another aspect of the method includes counting the number of samples having the predetermined phase shifts of 0°, 90°, −90° and 180° over a predetermined interval of time. Step (c) of the method may include delaying the signal C from itself by at least one of 2τ, 4τ, 8τ, 16τ, 32τ, or 64τ, in addition to the 1τ sample time, and step (d) may include counting the number of samples produced having predetermined phase shifts at the at least one of 2τ, 4τ, 8τ, 16τ, 32τ, or 64τ, in addition to the 1τ sample time.
Still another embodiment of the invention includes a method of determining a frequency of a signal having the steps of: (a) receiving a signal S; (b) forming in-phase (I) and quadrature (Q) signals from the signal S; (c) sampling the I and Q signals to form first and second digital signals; (d) combining the first and second digital signals to form a combined signal C; (e) delaying C from itself by one sample time τ; (f) counting numbers of samples having predetermined phase shifts of 0°, 90°, −90° and 180°, after delaying C; and (g) determining the frequency of the signal S, based on the numbers of samples counted in step (f) at each of the predetermined phase shifts. An aspect of step (c) includes sampling the I and Q signals using one bit ADCs during a predetermined time interval.
It is understood that the foregoing general description and the following detailed description are exemplary, but are not restrictive, of the invention.
The invention is best understood from the following detailed description when read in connection with the accompanying drawing. Included in the drawing are the following figures:
a is a diagram of a prior art IFM receiver;
b is a diagram of another prior art IFM receiver;
a–4d are graphical plots of values of n1, n2, n3 and n4 as a function of input frequency, where n1 is a 0° phase shift, n2 is a 90° phase shift, n3 is a −90° phase shift and n4 is 180° phase shift, occurring when the two digital output signals depicted in
a and 5b are graphical plots of values of n23 and n14 as a function of input frequency, where n23 is equal to n2−n3 and n14 is equal to n1–n4, occurring when the two digital output signals depicted in
c and 5d are graphical plots of values of n23 and n14 as a function of input frequency, where n23 is equal to n2−n3 and n14 is equal to n1−n4, occurring when the two digital output signals depicted in
a–6d are graphical plots of values of n23 and n14 as a function of input frequency, with delay times of 1τ and 4τ that are similar, respectively, to
a and 8b are graphical plots of values of n14 and n23 as a function of input frequency for a 1τ delay in the two digital output signals depicted in
c and 8d are graphical plots of values of n14 and n23 as a function of input frequency for a 64τ delay in the two digital output signals depicted in
a and 9b are graphical plots of values of n23 and n14 as a function of input frequency for 1τ delay in the two digital output signals depicted in
c and 9d are graphical plots of the values of n23 and n14 as a function of input frequency for 4τ delay in the two digital output signals depicted in
Referring to
As shown, the received signal passes through a 90° hybrid, generally designated as 21, to form a first analog signal I and a second analog signal Q. The second analog signal is delayed from the first analog signal by 90°. One bit analog-to-digital converter (ADC) 22 converts the first analog signal I into a first digital signal. Another one bit ADC, designated as 23, converts the second analog signal Q into a second digital signal. The first and second digital signals are provided to digital processor 25, which determines the input frequency of the received signal.
Generally, an IFM receiver, such as receiver 20, covers a wide bandwidth, for example 2 GHz and more. In order to achieve such wide bandwidth, the sampling frequency of receiver 20 must be high. The inventors realized that a high speed ADC having a low number of bits may be used to achieve speed for sampling the input signal. Accordingly, the inventors use a one bit ADC in the exemplary embodiment of the invention. In the exemplary embodiment of
In order to cover an even wider bandwidth, an in-phase (I) and a quadrature (Q) approach is used by the exemplary embodiment shown in
Receiver 20 is configured, in accordance with the present invention, to provide the phase angle of the input signal at every sampling time, as sampled by clock 24. ADC 22 and ADC 23 are each sampled by clock 24 in a manner that is similar to a D type flip flop providing an output Q signal every sample period. This advantageously simplifies the signal processing required by receiver 20.
In order to cover the desired 2 GHz bandwidth, for example, a sampling frequency of 2.56 GHz is used by clock 24. A corresponding sampling time of 1τ is about 0.39 ns ( 1/2.56 GHz), which is the inverse of the sampling frequency. The sampling time of 1τ may be used as the shortest delay time in the IFM receiver. It will be appreciated that one bit ADCs are available at sampling speeds of up to 50 giga-samples per second (GSPS), allowing the IFM receiver to process very wide bandwidth input signals. Other sampling times may also be used, such as 2τ, 4τ, 8τ, 16τ, 32τ, and/or 64τ, and more. Since all these sampling times are multiples of one sampling time interval (1τ), there is no special timing adjustment that must be made by receiver 20. All of these sampling times, if necessary, may be executed by digital processor 25.
Digital processor 25 provides the delay times to find the phase difference between the received input signal and its delayed version. This phase difference is used to find the frequency of the received input signal. The delay time provided by digital processor 25 is the interval between various sampling times. According to the above example, the shortest delay time equals one sampling time interval 1τ, which is 0.39 ns. Longer delay times may equal multiples of the one sampling time interval 1τ, namely 2τ, 4τ, 8τ, 16τ, 32τ, 64τ, etc.
Since ADC 22 and ADC 23 are advantageously only one bit each, the formed outputs of the first digital signal and the second digital signal are either a+1 level or a−1 level. These two output signals from the two ADCs may be considered as one real component and one imaginary component, as exemplified in
When the input frequency of received signal S is below 2.56 GHz, as for example, having an input frequency from 10 to 2540 MHz, the maximum phase shift between two adjacent sampling times is plus (+) or minus (−) 180°. Accordingly, between two sampling times, there are four possible phase shifts, namely, 0°, 90°, −90° and 180°. In an embodiment of the present invention, plus (+) or minus (−) 180° are considered as one angle equal to 180°. The phase shift may be determined from the angle difference between two adjacent sampling times of tn+1 and tn.
When the input frequency changes from 10 to 2540 MHz, for example, the number of phase shifts measured, or counted, by digital processor 25 over a predetermined interval of time (for example, 100 ns) is shown in
It will be appreciated that in the example shown in
Digital processor 25 combines the first digital signal and the second digital signal to form a combined signal, referred to herein as signal C. The results shown in
Values of n1, n2, n3 and n4, shown in
n14=n1−n4 and
n23=n2−n3.
The results of the number of samples of n23 versus input frequency are shown in
c depicts the number of samples of n23 versus input frequency and
It will be appreciated that the triangular forms of the curves, shown in
The inventors discovered that the output curves shown in
As previously described, the delay times of the combined signal C, when delayed from itself, may be chosen as 1τ, 4τ, 16τ, 64τ, etc. The delays of 1τ, 4τ, and 16τ, for example, may be used to resolve any frequency ambiguity and to generate coarse frequency resolution. The longest delay of 64τ may be used to generate the fine frequency resolution. It will be appreciated that there may be some advantage in ambiguity resolution by using delay times that are multiples of smaller intervals such as 2τ (for example, τ, 2τ, 4τ, etc.). This may require, however, more hardware resources within an FPGA (a fixed programmable gated array used as processor 25). Since the angle at every sampling time is calculated by processor 25, the phase differences between the different delay times may simply be found by the angle differences.
An exemplary frequency encoding method and noise effect used by the present invention is described in the following paragraphs. Referring to the top two figures of
for frequency region 0, n14>0 and n23>0;
for frequency region 1, n14<0 and n23>0;
for frequency region 2 n14<0 and n23<0; and
for frequency region 3 n14>0 and n23<0.
These four frequency regions may be clearly determined when the input signal-to-noise (S/N) is high. When the S/N is low, indeterminate frequency regions near the edges of the four determinate frequency regions may cause erroneous results. These indeterminate frequency regions are marked as A, B, C and D in
From experimental results, the inventors discovered that the widths of regions A, B, C and D are wider for shorter delay times and narrower for longer delay times. This phenomenon may be explained by referring to
Using the above frequency division, a method of the invention divides the 40 MHz range into 4 frequency regions with a fine frequency resolution of 10 MHz. In this frequency encoding method, the present invention assumes that the fine frequency (10 MHz) measured from the longest delay time is encoded correctly. In other words, the A, B, C and D regions in
An approach used by the present invention to resolve the ambiguous regions is to use both frequency ranges, going from the finer frequency resolution to the next higher coarse frequency resolution. In
Using the values shown in
Table 1. Conditions used to determine coarse frequency from finer frequency.
This table may be generated from
In the exemplary embodiment above, the fine frequency resolution is 10 MHz. Usually, finer frequency resolution may be desirable. In order to achieve this finer resolution, more than 2 bits may be generated from the longest delay time. This may be achieved by comparing n14 and n23 and dividing the frequency in finer ranges. If 4 bits are generated from the longest delay time, the corresponding frequency resolution is 2.5 MHz ( 10/4), and if 5 bits are generated from the longest delay time, the corresponding frequency resolution is 1.25 MHz ( 10/8).
It will be appreciated that the 1τ delay, by itself, may provide unambiguous frequency determination. The 2τ delay, on the other hand, has 2 ambiguities, whereas the 4τ delay has 4 ambiguities, etc.
Referring now to
Using step 1006 of the method, processor 25 delays the combined signal C from itself by one sample time 1τ. Using step 1007 next, processor 25 counts the number of samples at 0° phase, +90° phase, −90° phase and 180° phase during a predetermined interval of time. Such predetermined interval of time may be, for example, 100 ns.
Referring to step 1008 of the method, processor 25 counts the number of samples at 0° phase and lets that number equal n1. The processor counts the number of samples at +90° phase and lets that number equal n2. The processor counts the number of samples at −90° and lets that number equal n3. Finally, the processor counts the number of samples at 180° phase and lets that number equal n4. In step 1009 of the method, processor 25 computes n23, which is equal to n2−n3. Similarly, the processor computes n14, which is equal to n1−n4. Using the counted samples over the predetermined time interval, a LUT is consulted in step 1010, which contains results of the number of samples versus input frequency, as shown in
Referring next to
Step 1018 of the method determines the number of samples at 0° phase for one or more of the delays to combined signal C performed in step 1016. The number of samples determined at 0° phase is set to be n1. Similarly, the number of samples determined to be at +90° phase is set to be n2. The number of samples determined at −90° phase is set to n3. Finally, the number of samples determined at 180° phase is set to n4. Again, the number of samples at n1, n2, n3 and n4 are determined for each of one or more of the delays to combined signal C performed in step 1016. In step 1019 of the method, for each sample time delay performed in step 1016, the processor computes n23, which equal n2−n3 and computes n14, which equals n1−n4. Using the determined sample numbers, the method in step 1020 determines the frequency of the received signal by using another LUT, which may contain values similar to the values plotted in
Although the invention is illustrated and described herein with reference to specific embodiments, the invention is not intended to be limited to the details shown. Rather, various modifications may be made in the details within the scope and range of equivalents of the claims and without departing from the invention.
Number | Name | Date | Kind |
---|---|---|---|
4336541 | Tsui et al. | Jun 1982 | A |
4426648 | Tsui et al. | Jan 1984 | A |
4547727 | Tsui et al. | Oct 1985 | A |
4633516 | Tsui | Dec 1986 | A |
4928105 | Langner | May 1990 | A |
4963816 | Tsui et al. | Oct 1990 | A |
5099243 | Tsui et al. | Mar 1992 | A |
5109188 | Sanderson et al. | Apr 1992 | A |
5291125 | Tsui et al. | Mar 1994 | A |
5499391 | Tsui | Mar 1996 | A |
5793323 | Tsui | Aug 1998 | A |
5917737 | Tsui et al. | Jun 1999 | A |
6411076 | Rudish | Jun 2002 | B1 |