# Patent

## Title
System and method for obtaining accurate symbol rate and carrier phase, frequency, and timing acquisition for minimum shift keyed waveform

## Abstract

The present invention relates generally to communication systems, both wired and wireless, employing a continuous phase modulation (“CPM”) waveform with a minimum shift keying (“MSK”) preamble. The present inventive system and method uses information from contiguous Fourier Transforms taken on contiguous data blocks to determine baud rate, phase, frequency offset, and bit timing of the CPM waveform or can be used to determine the frequency of continuous wave waveform. More particularly, the inventive system and method is applicable to the military satellite communications UHF frequency band for deciding whether a signal of interest is.

## Background
RELATED APPLICATIONS
The instant application claims the priority benefit of U.S. application Ser. No. 10/158,923, filed Jun. 3, 2002 now U.S. Pat. No. 7,158,588, entitled “SYSTEM AND METHOD FOR OBTAINING ACCURATE SYMBOL RATE AND CARRIER PHASE, FREQUENCY, AND TIMING ACQUISITION FOR MINIMUM SHIFT KEYED WAVEFORM” by the same inventors, the entirety of which is incorporated herein by reference.
BACKGROUND OF THE INVENTION
The present invention relates generally to communication systems, both wired and wireless, employing a continuous phase modulation (“CPM”) waveform. One such CPM waveform is known as minimum shift keying (“MSK”) modulation. The present inventive system and method is applicable to all communication systems and radio frequency bands which utilize an MSK preamble, defined below, to determine baud rate, phase, frequency offset, and bit timing. More particularly, the inventive system and method is applicable to the military satellite communications UHF frequency band for deciding whether a signal of interest is present.
Many communication systems or networks, both wired (e.g., Ethernet) and wireless (e.g., HF, VHF, UHF radio), utilize a preamble to determine the modulation carrier frequency and phase. A MSK waveform with an alternating sequence, e.g., 1,1,0,0,1,1,0,0, . . . , has a characteristic frequency spectrum, sometimes referred to as the “MSK Tones” which also, in addition to carrier frequency and phase, provides modulation symbol rate and accurate baud timing of the MSK waveform. The preceding MSK alternating sequence may be written as [(1·2),(0·2)]m which may be generalized in the following form: [(1·n),(0·n)]m where the variable “n” may be referred to as the “symbol repetition factor” and the variable “m” may be referred to as the “symbol pair repetition factor”. Other MSK waveforms that fit this general pattern, e.g., 1,1,1,0,0,0,1,1,1,0,0,0, . . . , which can be written as [(1·3),(0·3)]m and 1,1,1,1,0,0,0,0,1,1,1,1,0,0,0,0, . . . , which can be written as [(1·40,(0·4)]m, are all contemplated by the present inventive system and method.
There are several problems with using an MSK waveform preamble for waveform detection and parameter estimation which must be overcome in order to allow for accurate measurement of symbol rate, frequency, phase, timing, and signal strength (measured, for example, as signal-to-noise ratio, carrier-to-interference ratio, etc.). In most communication systems it is highly desirable to limit the preamble time and begin to transmit information-carrying data as soon as possible since the preamble essentially displaces data-carrying capability. However, limiting the preamble time has the effect of limiting the signal energy of the preamble which increases the difficulty in finding the characteristic MSK tones in the frequency spectrum due to noise. The transition in energy at the start of the preamble exacerbates this problem because the unit step in energy at the start of the preamble induces tones in the frequency spectrum, making it harder to distinguish the MSK tones from the noise. This forces the acquisition algorithm to pay close attention to gain control because amplitude changes, which are usually implemented in gain “steps”, also add tones to the frequency spectrum. Additionally, tones are introduced to the frequency spectrum due to amplitude changes which are typically implemented in predetermined step changes by the gain control algorithm. Therefore, the gain control algorithm must attempt to minimize the number of step changes in order to minimize the introduction of tones into the frequency spectrum and thereby maximize the ability to detect the MSK tones.
The implementation of the Fourier Transform (“FT”) or the Laplace Transform (“LT”) are common methods to convert time domain data to frequency domain data for analysis. The FT has discrete bins which contain the energy as correlated with a frequency offset for each bin. The Fourier Transform may be implemented as a Discrete Fourier Transform (“DFT”) or as a Fast Fourier Transform (“FFT”) in those devices that are computationally-limited. It is to be understood that any method for converting time domain data to frequency domain data, such as Fourier Transforms, Laplace Transforms, discrete cosine, etc., are contemplated by the present invention and any method for converting time domain data to frequency domain data may be referred to herein, individually and/or collectively, as a Fourier Transform, or “FT” as would be consistent in the context as used.
The use of a FFT to find MSK tones at a specific spacing equal to one-quarter of the symbol rate of the received signal is described in the paper “An Innovative Synchronization Preamble for UHF MILSATCOM”, authored by Mark Miller, Mark A. Harris, and Donald R. Stephens (the “Miller Paper”), which is hereby incorporated herein by reference. This paper only describes the use of the FFT to find the MSK Tones and implies the use of a correlation function to determine the characteristic spacing of the MSK Tones. It describes the use of the FFT Bin number to find the Carrier Frequency and the Phase value of that center bin to find the Carrier Phase. It describes the use of the Phase difference between the carrier signal and the +/−1 MSK Tones to determine symbol phase (timing). However, use of the method disclosed in the Miller Paper results in a lot of false alarms, e.g., reporting a detection on noise, a foreign signal, or on an impulse signal. Part of the problem with the procedure used in the Miller paper is that the use of the center frequency of the bin with the largest amplitude is at best a rough estimate of the frequency of the carrier. The invention herein described includes the functionality of signal detection as well as a much more accurate method of signal parameter estimation. The present inventive system and method overcomes the problems of the prior art as detailed below.
Generally, there are two major issues with the FT approach to finding the characteristic MSK tones. The first problem deals with computational horsepower required to perform an FT. For example, a FT or DFT may be computationally burdensome, typically requiring N2 operations where N is the number of frequency bins. A FFT, by contrast, typically only requires N*log(N) operations and is therefore less computationally burdensome. However, with either an FT, DFT, or FFT, the fewer the number of bins the less frequency resolution is attainable. Therefore, an undesirable tradeoff is required between computational intensity and frequency resolution. Typical prior art solutions strike a compromise between computational intensity and frequency resolution by merely using the center frequency bin of the FT. The present inventive system and method overcomes the compromise problem by taking two contiguous FTs and, generally, determining the difference in phase for the bin of the center (carrier) frequency of the contiguous FTs to accurately determine the waveform frequency. The actual implementation used will be discussed in detail further below.
The second major issue with the FT implementation is that any frequency which is not an exact integer multiple of the data sample rate divided by the number of FT bins ends up with energy split between two adjacent bins of the FT. This effect tends to hide the characteristic MSK tones in the surrounding noise. The present inventive system and method solves this problem by multiplying the input data by a constant tone which may correspond to exactly 1⁄2 bin frequency, or some other fraction of bin frequency, and then performing another FT on the input data that has been multiplied by the constant tone and comparing the results with the results of the FT performed on the non-multiplied input data. The procedure creates two FT's on the same time-domain data. In one case, the carrier frequency will be located more closely to the center of an FT bin. In the other case, the carrier frequency will be located closer to the edge of the bin which may cause energy spillover into the next, adjacent FT bin. This energy spillover is undesired because it reduces the apparent signal strength in relation to the noise energy. The worst case situation occurs when tone appears exactly on a bin edge. In this case, the tone energy will be equally distributed between the two adjacent FT bins. For the single bin of interest (the transmitted carrier frequency), the apparent signal to noise ratio is 3 dB (a power ratio of 2) less than the signal to noise ratio which would be apparent if the FT was modified to locate the transmit carrier tone at the center of an FT bin. With two candidate FT's, choosing the transmit carrier frequency which has the highest magnitude will provide the best possible representation of the transmitted signal (in the frequency domain). This will improve the apparent signal to noise ratio which will improve the signal estimation performance in noisy environments.
The MSK preamble has an additional problem because a Fourier Transform of time domain data which is not an exact multiple of MSK symbols of the [(1·n),(0·n)]m pattern has an effect of varying the amplitude of the characteristic tones with respect to each other in a “walking” type pattern, i.e., the tone energy is constant but the tones repeated throughout the spectrum will not have a symmetric pattern in amplitude. This effect increases the difficulty in finding the characteristic (MSK) tones in the noise. To resolve the problem, the inventive system and method utilizes a polyphase resampler (interpolation and decimation method) to exactly place the sample rate of the time domain signal at an exact integer multiple of the symbol rate of the MSK waveform.
Thus there is a need for a system and method which can detect a continuous phase modulation waveform with a shortened MSK preamble and overcome the limitations of prior art systems/methods.
The present inventive system and method separates the detection and estimation functions. During detection, the baud rate is calculated as a first estimated parameter which is then utilized by the receiver for tailoring the signal sample rate and bandwidth to better match the incoming signal before calculating the transmitted carrier frequency, phase, and bit timing. The signal to noise ratio improvement which results from the tailoring process (of resampling and filtering) consequently improves the accuracy of the calculated values. The method described in the Miller paper does not separate the detection and estimation functions and therefore cannot accommodate the tailoring of the sample rate and bandwidth to match the incoming signal.
Additionally, the present inventive system and method performs a baud rate calculation by sorting the FT tones in order of amplitude and measuring the bin distance between the closest two tones. The expected value of the two closest tones is equal to the symbol rate divided by 2n where “n” is the number of bit repeats in the preamble from the [(1·n),(0·n)]m form of the MSK preamble. This baud rate calculation greatly improves the baud rate estimation algorithm as compared to the correlation described in the Miller paper because the inventive method takes advantage of the large signal to noise ratio apparent in a FT bin that contains a large signal level. The correlation method from the Miller paper integrates the noise in the entire range of frequencies which cover the transmit carrier frequency (e.g., +/−1500 Hz) and the MSK Tones (e.g., +/−28,000 Hz), which results in a sum of +/−29,500 Hz added into the correlation function. Additive White Gaussian Noise is, by definition, equally distributed in each frequency bin and the correlation process proposed in the Miller paper would integrate the noise in all of those bins. The sorting method used by the present invention excludes energy in those bins which do not contain enough signal energy to cause a signal detection. Therefore, the present invention excludes the noise from the baud rate estimation calculation as will become apparent in the detailed description below.
Furthermore, the detection calculation of the present invention determines the amplitude of largest tone (typically the carrier frequency) and adds the amplitudes of the +/−1 MSK Tones and the +/−2 MSK Tones in the signal detection calculation. The Miller paper does not address detection of the incoming signal at all.
The present invention also uses two adjacent FT windows to measure the phase difference of the carrier frequency between the two windows to thereby accurately determine the carrier frequency. Moreover, the present invention uses 1⁄2 tone spacing in a third FT window in order to place the amplitude of the carrier frequency in the most advantageous location so that carrier energy is not dispersed between two FT bins, as it would be if the carrier frequency were located at (or near) the edge of a bin. It is to be understood that the present invention contemplates offsetting the third FT by any amount, and is not limited in any way to just 1⁄2 tone spacing, so as to place the carrier frequency in the middle of the bin. The Miller paper is silent on the use of adjacent FT bins as well as on the use of a third FT spaced apart by a half tone.
One embodiment of the present invention avoids the problems of the prior art by using two or more contiguous Fourier Transforms for detecting a continuous phase modulation waveform and determine the characteristics of the waveform such as frequency, phase, timing, and signal strength.
Accordingly, it is an object of the present invention to obviate many of the above problems in the prior art and to provide a novel system and method for detecting a continuous phase modulation waveform.
It is another object of the present invention to provide a novel system and method for determining the waveform characteristics of a continuous phase modulation waveform.
It is yet another object of the present invention to provide a novel system and method for using two or more Fourier Transforms in the detection of a continuous phase modulation waveform with an MSK preamble.
It is still another object of the present invention to provide a novel system and method for using two or more Fourier Transforms to determine the characteristics of a received continuous phase modulation waveform with an MSK preamble.
It is a further object of the present invention to provide a novel system and method for using two or more Fourier Transforms for acquiring a continuous phase modulation waveform with an MSK preamble.
It is yet a further object of the present invention to provide a novel system and method of determining the frequency of a received and detected continuous wave waveform.
These and many other objects and advantages of the present invention will be readily apparent to one skilled in the art to which the invention pertains from a perusal of the claims, the appended drawings, and the following detailed description of the preferred embodiments.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a diagram showing an incoming signal being divided into blocks, contiguous Fourier Transforms (FT1 and FT2) being taken, resampling the blocks and another set of contiguous Fourier Transforms (FT3 and FT4) being taken, multiplication of a block of data by a tone and another Fourier Transform (FT5) being taken.
FIGS. 2A through 2E are a flow chart indicating the major steps for detecting a signal according to the present invention (FIGS. 2A and 2B) and the major steps for estimating characteristics of a signal according to the present invention (FIGS. 2C through 2D).
FIG. 3 is a flow chart indicating the major steps in determining the frequency of a received continuous wave waveform according to the present invention.



DESCRIPTION OF PREFERRED EMBODIMENTS
With reference to the drawings, like numerals represent like components/steps throughout the several drawings.
Many communication systems, both wired and wireless, use Continuous Phase Modulation (“CPM”) waveforms. One particular type of CPM waveform is the Minimum Shift Keyed (“MSK”) waveform. As is known in the art, an MSK waveform with an alternating sequence, such as 1,1,0,0,1,1,0,0, . . . , generalized to the form: [(1·n),(0·n)]m, also referred to herein as an “MSK preamble”, has a characteristic frequency spectrum from which can be used to determine parameters of the received signal, such as carrier frequency, carrier phase, modulation symbol rate and baud rate. The present inventive system and method makes use of two or more contiguous Fourier Transforms, Discrete Fourier Transforms, and/or Fast Fourier Transforms to determine the aforementioned parameters. It is to be understood that any MSK waveform of the generalized form, i.e., where n=1, 2, 3, . . . , is contemplated by the present inventive system and method.
An MSK waveform of the above generalized form has a frequency spectrum which has a tone at the carrier frequency and characteristic tones which are spaced at the waveform symbol rate divided by 2n. For example, for the alternating sequence 1,1,0,0,1,1,0,0, . . . , the characteristic tones are spaced apart from the carrier frequency by the waveform symbol rate divided by 4. Although this property of the MSK waveform is well documented, the present inventive system and method describes an approach to the accurate detection and estimation of the MSK preamble that is both innovative and easily implemented.
If the input to a FT is a time domain signal, the output of an FT is a frequency domain signal in a predetermined number of bins. Each bin, by itself, can give an estimate of the frequency and phase of the incoming signal. However, the estimate is typically too inaccurate to allow for accurate parameter estimation. The FT frequency domain output counts the number of cycles of the input time domain waveform. For example, bin 5 of a FT will typically represent frequencies between 41⁄2 and 51⁄2 Hz, or cycles. Bin zero (sometimes referred to herein as bin 0), typically the first bin of the FT, is the bin that includes a zero offset between the transmitter and the receiver, which implies that the transmitter and receiver are synchronized. For example, for a FT with 256 bins, numbered 0 to 255, bin 0 is the zero offset bin with the offset increasing from bin 1 through bin 255 such that bin 255 has the same magnitude, but opposite direction, offset of bin 1. Therefore, the bins can be ordered in a series with bin 0 in the middle of the series and the same magnitude offsets being the same number of bins away from bin 0. For example, the ordering of the FT bins in a series for a maximum offset of 6 bins in either direction from the zero offset bin (bin 0) would be:
























BIN
BIN
BIN
BIN
BIN
BIN
BIN 0
BIN 1
BIN 2
BIN 3
BIN 4
BIN 5
BIN 6


250
251
252
253
254
255









It is to be understood that in the above explanation, the use of 256 bins is exemplary only and the present invention contemplates the use of any number of bins. Additionally, the placement of Bin 0 is arbitrary and can be placed anywhere in the series of bins.
The present inventive system and method makes use of contiguous, non-overlapping FTs, from contiguous, non-overlapping samples, in order to accurately determine the parameters of the received signal. In doing so, the magnitude and phase of the various FT bins are utilized. For example, the phases of the contiguous FTs are used to determine frequency of the received signal. The more contiguous FTs that are used, the higher the resolution and the more accurate the estimate of the frequency. It is to be understood that although the present description may recite only two contiguous FTs, the present inventive system and method contemplates the use of two or more contiguous FTs.
With attention now to FIG. 1, an incoming signal 10 is received by the receiver. What is meant herein by “received” is that the raw incoming data, at the predetermined sampling rate, enters the processing circuitry at the receiver. The incoming signal may be any CPM signal with an MSK preamble transmitted over any medium. Specifically, the signal may be a signal on the Military Satellite Communications UHF Frequency band. In order to detect the incoming signal 10 and determine the characteristics of the signal, the signal is divided into N equal blocks of a predetermined number of samples, as is known in the art. One embodiment of the present invention utilizes a block size of 256 samples.
The block of data designated Block 1 may generally be thought of as the block of samples that contain the initial received energy of the incoming CPM signal. The subsequent blocks are each contiguous with their immediate neighbors. The first Fourier Transform, designated as FT1, is taken on the data Block 1. FT2 then begins with the sample immediately after FT1 ends, so that FT1 and FT2 are contiguous without any overlap. The information derived from either FT1 or FT2 as well as other known information about the incoming signal, such as maximum allowable carrier offset and maximum symbol rate, are sufficient for detecting the incoming signal 10. A FT may be performed on the other data blocks, but are not shown in FIG. 1 for clarity purposes.
The incoming signal 10, or parts thereof, are resampled to thereby create the resampled incoming signal 20. Typically, data Block 1 and data Block 2, which were used for detection purposes, will not be resampled, although resampling those two data blocks is contemplated by an embodiment of the present invention. The resampled rate is typically lower than the sampling rate in order to ensure adequate sampling of the incoming signal. However, it is contemplated that the resampled rate may also be the same as or faster than the original sampling rate. The FT3 is taken on the resampled data Block 3 and the FT5 is taken on the resampled data Block 4 such that the FT3 and the FT5 are contiguous. A FT may also be performed on the other resampled data blocks, but those FTs are not shown in FIG. 1 for the sake of clarity. The resampled data Block 3 is modified by multiplying by a tone, exp(jω/s) where “ω” is the frequency and “s” is the new sample rate, then FT4 is taken on the modified resampled data Block 3. Preferably, the tone shifts the data in resampled Block 3 by 1⁄2 the bin spacing. The 1⁄2 bin spacing shift is exemplary only and is not meant to limit the application in any way. Other bin spacings are contemplated, but shifting by 1⁄2 bin spacing allows for the greatest amount of resolution as will become apparent in the discussion below.
As discussed in more detail below, the information from the FT3, FT4, and FT5 operations, in addition to the information provided above for detection, allow the characteristics of the incoming signal 10 to be accurately estimated.
With reference now to FIGS. 2A and 2B, a flow chart is depicted showing the major steps for detection of a CPM waveform with an MSK preamble. At the step 202, a sample rate is chosen. Since the symbol rate and maximum allowable carrier offset for the incoming signal is unknown at the time of detection, the sample rate that is chosen must be high enough to allow for (a) the maximum allowable carrier frequency excursion as well as (b) the MSK characteristic tones from the MSK preamble. The MSK tones are typically spaced apart by the baud rate, which for an MSK preamble of the form [(1·n),(0·n)]m, where n=2, i.e., a preamble of the form 1,1,0,0,1,1,0,0, . . . , is the symbol rate divided by 4. Generally, the baud rate equals the symbol rate divided by 2n. For example, for an MSK preamble with a 28,000 symbols/sec symbol rate and a ±1500 Hz maximum allowable carrier offset, a sample rate of 64,000 complex samples/sec (e.g. I and Q samples) will provide a good initial frequency estimate for the incoming signal. It is to be understood that the above example is not intended to limit the present invention in any manner but rather is included to assist in the understanding of the invention.
During the step 204 partitions the incoming signal into data blocks each containing a predetermined number of samples. One useful partitioning of samples is 256 samples per block, since 256 is a power of 2, although other block sizes are also contemplated by the present invention. As shown in FIG. 1, the incoming signal 10 is partitioned into N blocks of S samples.
At the step 206, the FT bins are normalized by any means known in the art, such as: [I2+Q2]1/2; or adding the larger of I or Q to 1⁄2 of Q or 1⁄2 of I, respectively, so that the sum of the amplitudes of all the samples divided by the bin size is unity. There can be any number of bins for the FT, consistent with the discussion above for computation intensity and resolution. For a block size of 256 samples/block, a preferable number of FT bins is 256. Alternative numbers of bins are contemplated by the present invention. Therefore, for a sampling rate of 64,000 samples/sec and 256 bins, the bin size, or frequencies (or samples) per bin, is 250. It is preferable, though not required, to choose the sample rate such that the allowable symbol rates for the incoming signal all appear as an exact integer multiple of the bin size, Hz/bin. As would be known to those of skill in the art, the FT used would typically be a complex FT, although other types of FT are also contemplated by the present inventive system and method. The RMS amplitude difference between the samples in a first portion of data Block 1 and the samples in a second portion of data Block 1 is also determined. The first portion of data Block 1 comprises those samples in data Block 1 that are not included in FT1. The second portion of data Block 1 comprises those samples in data Block 1 that are included in FT1. The RMS amplitude difference will be used later as one indication of detection of the incoming signal.
At the step 208, the Tone Bin Distance (“TBD”) is determined. The TBD is the number of bins between the MSK tones as determined by the baud rate. The TBD is approximately equal to the baud rate (samples/sec, or Hz) divided by the bin size (Hz/bin). The TBD will be used, as discussed below, to determine a cumulative amplitude which is used as one indication of detection of the incoming signal. A first FT, FT1 in FIG. 1, is performed on the first data block, Block 1 in FIG. 1, that typically contains the initial energy of the incoming signal. The FT1 is performed by any means known in the art. If the signal does not appear in the Block 1, the FT1 is discarded, along with the samples in the Block 1, and the FT2 performs the same tasks as would have FT1. Therefore, the procedure can be regarded as Block 1 is the block of data within which the incoming signal 10 appears.
At the step 210, the baud rate of the incoming signal is determined. The baud rate may be determined by sorting the bins of FT1 in amplitude order, designating a set of the largest amplitude bins as containing the MSK tones, and determining the closest bin distance (by the difference in bin numbers) between any two of the adjacent MSK tone bins. Typically, the three or five bins with the largest amplitudes are designated as the MSK tones, although other numbers of bins are contemplated. For example, if the five bins with the largest amplitudes, sorted in amplitude order, are bin 0, bin 253, bin 4, bin 9, and bin 248, the closest bin distance between any two of the bins is 3 bins between bin 0 and bin 253. The baud rate can then be determined by the closest bin distance (number of bins) and the bin size (samples/bin). If the determined baud rate is approximately equal to an expected symbol rate divided by 2n, then there is possible signal detection, as will be discussed further below. Once the baud rate is determined, the bandwidth of the incoming signal can be determined by techniques known in the art, and a variable bandwidth filter can be adjusted so as to take advantage of the knowledge of the bandwidth of the incoming signal.
With reference now to the step 216, a Carrier Frequency Window (“CFW”) is determined. The CFW contains a number of bins as a function of the maximum allowable carrier offset. For example, if the maximum allowable carrier offset is 1500 Hz, and the bin size (as determined above) is 250 Hz/bin, then a Carrier Bin Offset (“CBO”) is 6 bins. The CFW is comprised of the zero bin offset, bin 0, plus those bins that are one CBO greater than and one CBO less than bin 0. For this example, the CFW is shown below:
























BIN
BIN
BIN
BIN
BIN
BIN
BIN 0
BIN 1
BIN 2
BIN 3
BIN 4
BIN 5
BIN 6


250
251
252
253
254
255









The above example is for illustrative purposes only and is not to be construed as limiting the invention in any way.
With attention now drawn to the steps 218 and 220, the cumulative amplitude of each bin in the CFW, including the amplitudes of the MSK tone for each respective bin, is determined. The use of the amplitudes of just the bins in the CFW does not result in a sufficiently low error rate to be useful. Therefore, the amplitudes of the respective ±MSK tones are added in to allow for a more robust detection algorithm. The present invention contemplates determining the cumulative amplitude of the CFW bin and the ±1 MSK tones as well as determining the cumulative amplitude of the CFW bin and the ±1 and ±2 MSK tones.
As an example, assume that the CFW is as shown above (where the maximum allowable carrier offset is 1500 Hz), for a FT with 256 bins, and a TBD of 28 bins (as determined from a baud rate of 7000 Hz divided by a bin size of 250 Hz/bin). The TBD of 28 indicates that an MSK tone should appear every 28 bins from bin 0, i.e., bins 28 and 228, 56 and 200, 84 and 172, etc. Therefore, for this example, in determining the cumulate amplitude for each of the bins in the CFW, including the ±1 MSK tones, the amplitudes of the following bins (row wise) would be added together:












Bin # of
Bin # of −1
Bin # of


Bin in
MSK
+1 MSK


CFW
Tone
Tone

















250
222
22


251
223
23


252
224
24


253
225
25


254
226
26


255
227
27


0
228
28


1
229
29


2
230
30


3
231
31


4
232
32


5
233
33


6
234
34









It shall be understood by those of skill in the art that the above example is illustrative only and that the procedure exemplified can be applied to other assumed values.
Alternatively, the determination of the cumulative amplitude can include the bins immediately adjacent to the bin in the CFW and/or the bins adjacent to the ±1 MSK tones and/or the bins adjacent to the ±2 MSK tones (if the bins for the ±2 MSK tones are used in calculating the cumulative amplitude. The use of the adjacent bins is useful in the event the carrier frequency is located at the edge of a bin. As an example, the cumulative amplitude can be determined for each of the bins in the CFW, including the ±1 MSK tones, and the adjacent bins for both the bins in the CFW and the ±1 MSK tones by adding the amplitudes of the following bins (row wise):















Bin # of

Bin # of −1

Bin # of



Bin in
Adjacent
MSK
Adjacent
+1 MSK
Adjacent


CFW
Bins
Tone
Bins
Tone
Bins




















250
249, 251
222
221, 223
22
21, 23


251
250, 252
223
222, 224
23
22, 24


252
251, 253
224
223, 225
24
23, 25


253
252, 254
225
224, 226
25
24, 26


254
253, 255
226
225, 227
26
25, 27


255
254, 0 
227
226, 228
27
26, 28


0
255, 1 
228
227, 229
28
27, 29


1
0, 2
229
228, 230
29
28, 30


2
1, 3
230
229, 231
30
29, 31


3
2, 4
231
230, 232
31
30, 32


4
3, 5
232
231, 233
32
31, 33


5
4, 6
233
232, 234
33
32, 34


6
5, 7
234
233, 235
34
33, 35









With reference now to FIG. 2B, specifically to the steps 222, 224, and 226, a determination of whether the incoming signal is detected is made. In order to show detection of the signal:

    
    
        1. at the step 222, the baud rate as determined in the step 210 must be equal to approximately a known symbol rate for the incoming signal divided by 2n where n defines the MSK preamble [(1·n),(0·n)]m;
        2. at the step 224, the FT window does not experience a large change of amplitude as determined by the difference in RMS amplitude from the step 212; and
        3. at the step 226, the largest cumulative amplitude from the steps 218 and 220 is less than a predetermined threshold.
    
    


At the step 228, if all three of the above tests are passed, the incoming signal has been detected. If any of the above tests are not passed, there is no detection, as indicated at the step 230.
With attention now to FIGS. 2C and 2D, the flow chart from FIGS. 2A and 2B is continued. FIGS. 2C and 2D depict the major steps for determining certain characteristics of a CPM waveform with an MSK preamble. Alternatively, the procedure depicted in FIGS. 2C and 2D can be performed independently from the detection procedures in FIGS. 2A and 2B, so long as the incoming waveform has been detected.
At the step 250 in FIG. 2C, the data blocks of the incoming signal, e.g., the incoming signal 10 in FIG. 1, are resampled at a predetermined resampling rate that is preferentially greater than the sampling rate chosen in the step 202. In one embodiment, the blocks of data that are resampled are the last half of the block in which the incoming signal 10 appears and enough of the remaining data blocks after that point for completing the determination of the waveform characteristics of the incoming signal. Typically, this is at least 2 data blocks, as shown in FIG. 1. Preferentially, the resampling rate is 2n where n, the symbol repetition factor, is from the aforementioned general equation for the MSK preamble [(1·n),(0·n)]m. However, the resampling rate may be any reasonable rate as contemplated by the invention. The resampling changes the time base of the data in the data blocks.
The resample rate is preferentially a multiple of the MSK signal symbol rate. For example, if there are four samples per symbol in the radio and using the previous example of an MSK preamble with a 28,000 symbols/sec rate, the resample rate would preferentially be 112,000 samples/sec. By requiring that the resample rate be an integer multiple of the symbol rate removes the MSK tone “walking” phenomenon, which is caused by boundary effects of the FT bins. Prior art solutions to this problem incorporate the use of FT “windows” to overcome the MSK tone “walking” phenomenon. However, FT “windowing” degrades the signal energy. The tying of the resample rate to an integer multiple of the symbol rate has the additional benefit of filtering out out-of-band noise and thus improving the signal-to-noise ratio of the received signal.
With reference now to the steps 252 through 268, loops in the procedure occur for taking FT3, FT4, and FT5, which are shown in FIG. 1. The operations performed in the steps 252, 254, and 258 correspond to the operations performed on the steps 206, 208, and 216 and will not be discussed further here. The step 260 corresponds to the step 218 with the insertion of the substep of determining the largest cumulative amplitude for the FT for which the loop is currently being run. The step 262 is a decision step to determine when FT4 will be run. For example, during the first pass through the loop, FT3 is taken on the resampled data Block 3 of FIG. 1. As previously discussed for FIG. 1, another FT, FT4, will be taken on resampled data Block 3 after the samples have been multiplied by a tone in the step 264: exp(jω/s) where “ω” is the frequency and “s” is the resampled rate, then FT4 is taken on the modified resampled data Block 3. Preferably, the tone shifts the data in resampled Block 3 by 1⁄2 the bin spacing, although other shifts of the data may be used as appropriate. After the loop of the steps 252 through 260 is repeated for FT4, the step 262 is answered in the affirmative thereby opening the path to the step 266. If an FT is to be performed on the next data block, which it will for FT5, the step 268 increments the next block of data and FT5 is performed. Note that FT5 is not performed on samples that have been multiplied by the tone exp(jω/s). The loop can be repeated as many times as necessary in order to obtain the data required for determining the characteristics of the incoming signal, as discussed further below.
Alternatively, since ideally FT5 does not need to pass through all the steps in the loop, the only steps in the loop that need to be performed for FT5 are the normalizing, step 252, and the performing the FT, step 254. This is indicated in FIG. 2C by the dotted line between the step 254 and the step 270.
With reference now to the step 270, the bin containing the carrier frequency is determined from the bins of FT3 and FT4. The bin with the largest cumulative amplitude from FT3 and FT4 is the bin that contains the carrier frequency. The center frequency of the bin that contains the carrier frequency is the carrier frequency estimate.
With attention now directed at FIG. 2D, the carrier frequency estimate from the step 270 must be refined by determining a phase difference in the step 272 and a frequency offset in the step 274.
At the step 272, the phase difference is determined between the phase value of the carrier frequency bin and the phase value of the bin from FT5 that has the same bin number as the carrier frequency bin. This phase difference corrects for the FT frequency which is not located at the exact center of the FT bin. At the step 274, the frequency offset is determined by dividing the phase difference in the step 272 by the time duration of one of the FTs. At the step 276, the carrier frequency is determined by adding the frequency offset from the step 274 to the carrier frequency estimate from the step 270.
With reference now to FIG. 2E, the phase and timing of the incoming signal will now be determined from the steps 278 through 286. At the step 278, the phase of the bins containing the ±1 MSK tones for the carrier frequency bin is determined. One method of determining the phase is by taking arctan(Q/I) for these bins, although any method for determining phase can be used. At the step 280, a delay offset is determined by calculating the difference in phase between the ±1 MSK bins and dividing the difference in phase by the bin distance between the ±1 MSK tones, which may be the same as twice the TBD. At the step 281, the baud rate offset is determined, by any method known in the art. The baud rate offset is a timing offset which is used to determine where the FT started with respect to the start of the preamble. At the step 282, the phase of the carrier frequency bin is determined by known methods and at the step 284, a modified phase for the carrier frequency bin is determined by subtracting the delay offset determined in the step 280 from the phase of the carrier frequency bin determined in the step 282. At the step 285, the frequency offset calculated in the step 274 is subtracted from the modified phase for the carrier frequency bin determined in the step 284 to thereby determine the phase of the incoming signal at the step 286. With reference now to the step 288, the signal level of the incoming signal is determined from the gain values from one of the FT bin normalizing steps. Accordingly, the frequency, phase, and signal level characteristics of the incoming signal are determined.
With attention now drawn to FIG. 3, the present inventive system and method can also be used to determine, by the use of contiguous Fourier Transforms, the frequency of a received and detected continuous wave (“CW”) waveform. At the step 302, the CW waveform is received and detected by any means known in the art. At the step 304, a first FT is performed on a first set of samples. At the step 306, a second FT is performed on a second set of samples where the first and second set of samples are contiguous, as are the first and second FTs. At the step 308, the bin number of the bin with the largest amplitude, from both the first and second FTs, is determined. At the step 310, a phase difference is determined between the phase values for the bin with the largest amplitude, as determined in the step 308, and the bin with the same bin number as the largest amplitude bin but is from the FT that does not contain the bin with the largest amplitude. At the step 312, a frequency offset is determined by dividing the phase difference determined in the step 310 by the time duration of one FT. At the step 314, the frequency of the CW waveform is determined by adding the frequency offset determined in the step 312 to the center frequency for the bin with the largest amplitude determined in the step 308.
It can be seen by those of skill in the art that the carrier frequency of a CPM waveform with a MSK preamble can be determined by information derived from contiguous FT blocks.
While preferred embodiments of the present invention have been described, it is to be understood that the embodiments described are illustrative only and that the scope of the invention is to be defined solely by the appended claims when accorded a full range of equivalence, many variations and modifications naturally occurring to those of skill in the art from a perusal hereof.

## Claims

1. In a method for determining a carrier frequency of a continuous phase modulation (“CPM”) waveform containing a binary data sequence comprising a pattern of a predetermined number (N) of first binary tones followed by N second binary tones, the pattern repeating a predetermined number of times, said method using a Fourier Transform (“FT”), the improvement comprising the steps of:
performing a first FT of said CPM waveform to generate a first set of FT bins and a second FT of said CPM waveform to generate a second set of FT bins, wherein said second FT is contiguous with said first FT;
determining a phase difference between phase values from a first FT bin of said first set of FT bins and a second FT bin of said second set of FT bins wherein said first FT bin and said second FT bin have the largest amplitudes in the first set and second set of FT bins, respectively;
dividing said phase difference by the time duration of either said first FT or said second FT to determine a frequency offset; and
determining the carrier frequency of said CPM waveform at least in part from said first FT and said second FT by adding said frequency offset to a frequency of the first or second FT bin having the largest amplitude.

2. The method of claim 1 wherein the CPM waveform is a minimum shift keyed waveform.

3. The method of claim 1 wherein the Fourier Transform is a Discrete Fourier Transform (“DFT”).

4. The method of claim 1 wherein the Fourier Transform is a Fast Fourier Transform (“FFT”).

5. The method of claim 1 wherein said predetermined number of tones is two.

6. The method of claim 1 wherein said predetermined number of tones is greater than two.