This invention is directed to a continuous phase modulation detector. In particular, this invention is directed to a method for continuous phase modulation detection. More particularly, this invention is directed to a multi-h continuous phase modulation detector.
The Advanced Range Telemetry (ARTM) program is a United States Department of Defense tri-service telemetry modernization project whose goal is to assure that all testing and training ranges are able to use telemetry as necessary to carry out their respective missions. Multi-h Continuous Phase Modulation (CPM) has been selected by the ARTM Joint Programs Office as the Tier II ARTM waveform, because it offers significant improvements over both legacy telemetry waveforms such as pulse width modulation/frequency modulation (“PCM/FM”) and the previous Tier I waveform known as the Feher-patented quadrature-phase-shift keying (“FQPSK”) in terms of spectral containment and detection efficiency, while retaining a constant envelope characteristic.
The ARTM Tier II modulation format is a multi-h continuous phase modulation. Those skilled in the art will appreciate that the multi-h continuous phase modulation format has a constant envelope and narrow bandwidth. Current implementations of receivers for multi-h continuous phase modulation experience several difficulties, including that the branch metrics are solely a function of the data in the multi-symbol observation window. That is, the influence of previous observations is not passed along in the form of a cumulative path metric. The skilled artisan will appreciate that the performance improves as the multi-symbol observation length increases; however, the penalty for this is that trellis complexity increases exponentially with increasing observation length. In addition, the current implementations perform poorly for practical multi-symbol observation lengths with respect to the Advanced Range Telemetry Tier II modulation format. Thus, the existing optimal maximum likelihood sequence estimation receiver for continuous phase modulation may have high complexity, both in trellis size and coherent demodulation requirements.
In view of the aforementioned needs, there is provided in accordance with the present invention an improved, noncoherent receiver capable of allowing multi-symbol observation.
In accordance with the present invention, there is provided a continuous phase modulation detector.
Further, in accordance with the present invention, there is provided a method for continuous phase modulation detection.
Still further, in accordance with the present invention, there is provided a noncoherent receiver capable of allowing multi-symbol observation.
In accordance with the present invention, there is provided a continuous phase modulation detector. The continuous phase modulation detector includes receiver means adapted to receive digitally modulated signals having a generally continuous phase. The detector also includes observation means adapted to perform multi-symbol observations on received digitally modulated signals. Memory means are included in the detector and adapted to store historic observation data corresponding to multi-symbol observations performed by the observation means. The detector further includes adjustment means.
In one embodiment of the present invention, the receiver means is noncoherent and preferably has a trellis structure. The observation means allow for adjusting of a multi-symbol observation length and provide for acquiring cumulative observation data. In a preferred embodiment, controlled use of acquired cumulative observation data is provided, wherein the reliance on past observations is adjusted recursively in accordance with cumulatively acquired observation data. Preferably, the adjustment is based on a “forget factor”. Using the adjusted cumulative metric, the detector of this embodiment is able to perform well while keeping the multi-symbol observation length to a minimum. In one embodiment complex-valued cumulative observation data is evaluated. In another preferred embodiment evaluation of real-valued observation data is performed. These embodiments are equally applicable to both PCM/FM and ARTM Tier II waveforms. In the context of PCM/FM, a two-symbol observation length (4 trellis states) is a few tenths of a dB inferior to the optimal maximum likelihood sequence estimating receiver, and is 3.5 dB superior to conventional FM demodulation. In the context of ARTM Tier II, the same two symbol observation length (64 states) is 2 dB inferior to the maximum likelihood sequence estimating receiver and 4 dB superior to FM demodulation.
Further, in accordance with the present invention, there is provided a method for continuous phase modulation detection. The method begins with the receipt of digitally modulated signals having a generally continuous phase. In a preferred embodiment of the present invention, a noncoherence reception of digitally modulated signals is provided. Multi-symbol observations are then performed on the received digitally modulated signals. In accordance with a predetermined performance, a multi-symbol observation length is adjusted and cumulative observation data resulting from multi-symbol observations is then acquired. Historic observation data corresponding to multi-symbol observations performed on the digitally modulated signals is then stored in a memory. In a preferred embodiment, the amount of acquired cumulative observation data being stored is selectively adjusted according to the stored historic observation data.
In this embodiment of the present invention, the use of a cumulative metric is controlled, wherein the reliance on past observations is adjusted recursively according to the cumulatively acquired observation data. In the preferred embodiment, the adjustment is based on a forget factor. Acquired cumulative observation data is evaluated, wherein in one embodiment complex-valued observation data is evaluated. In another preferred embodiment, evaluation is performed for real-valued observation data.
Still other objects and aspects of the present invention will become readily apparent to those skilled in this art from the following description wherein there is shown and described a preferred embodiment of this invention, simply by way of illustration of one of the best modes suited for to carry out the invention. As it will be realized by those skilled in the art, the invention is capable of other different embodiments and its several details are capable of modifications in various obvious aspects all without from the invention. Accordingly, the drawing and descriptions will be regarded as illustrative in nature and not as restrictive.
The subject invention is described in connection with the attached drawings which are for the purpose of illustrating the preferred embodiment only, and not for the purpose of limiting the same, wherein:
The present invention is directed to a noncoherent receiver capable of allowing multi-symbol observation. In particular, the present invention is directed to a continuous phase modulation detector and method for continuous phase modulation detection.
Continuous phase modulation refers to a general class of digitally modulated signals in which the phase is constrained to be continuous. The complex-baseband signal is expressed as:
where T is the symbol duration, h(i) are the modulation indices, α={αi} are the information symbols in the M-ary alphabet {±1, ±3, . . .±(M−1)}, and q(t) is the phase pulse. The subscript notation on the modulation indices is defined as:
h(i)≡h(i mod N
where Nh is the number of modulation indices (for the special case of single-h continuous phase modulation, Nh=1). The phase pulse q(t) is related to the frequency pulse f(t) by the relationship:
The frequency pulse is time-limited to the interval (0,LT) and is subject to the constraints:
In light of the constraints on f(t) and q(t), Equation (2) is suitably written as:
The term θ(t,αn) is a function of the L symbols being modulated by the phase pulse. For h(i)=2k(i)/p (k(i), p integers), the phase state θn−L takes on p distinct values 0, 2π/p,2·2π/p, . . . , (p−1) 2π/p. The total number of states is pML−1, with M branches at each state. Each branch is defined by the L+1-tuple σn=(θn−L, αn−L+1, αn−L+2, . . . , αn). The Advanced Range Telemetry Tier II modulation is M=4, h={4/16, 5/16} (Nh=2), 3RC (raised cosine frequency pulse of length L=3).
In accordance with the present invention, the model for the received complex-baseband signal is denoted by the equation:
r(t)=s(t,α)ejφ(t)+n(t) (7)
wherein n(t)=x(t)+jy(t) is complex-valued additive white Gaussian noise with zero-mean and single-sided power spectral density N0. The phase shift φ(t) introduced by the channel is unknown in general.
Those skilled in the art will appreciate that there are a plurality of instances wherein this signal model is considered. For example and without limitation, the binary continuous phase frequency shift keying (“CPFSK”)case assumes φ(t) to be uniformly distributed over the interval [−π,π]. It is also assumed to be slowly varying so that it is constant over a multi-symbol observation interval NT. The receiver correlates the received signal against all possible transmitted sequences of length NT and outputs the maximum likelihood decision on the middle bit in the observation.
With respect to the more general continuous phase modulation example, φ(t) is modeled as a slowly varying process with the Tikhonov distribution. The Tikhonov distribution is parameterized by β and has three important special cases: the fully coherent case where β=∞, the noncoherent case where β=0 and φ(t) reduces to a uniformly distributed value over [−π, π], and the partially coherent case where 0<β<∞. A practical receiver is then given for the noncoherent case (β=0), which is a generalization of the CPFSK receiver. This more general receiver has the complex-valued decision variable:
where {tilde over (α)} is a hypothesized data sequence and the observation interval is N1+N2=N symbol times. The term {tilde over (θ)}k−L accumulates the phase of the hypothesized symbols after they have been modulated by the length-LT phase pulse e−jθ(τ,{tilde over (α)}); it is necessary to match the phase of the individual length-T segments of the integral in Equation (8). Equation (9) shows that this metric is suitably computed recursively using the Viterbi algorithm with a trellis of ML+N−2 states. It is important to point out that the recursion does not maintain a cumulative path metric, but rather functions as a sliding window that sums N individual length-T correlations (each rotated by the proper phase). The receiver does not perform a traceback operation to determine the output symbol, but instead outputs the symbol {tilde over (α)}n corresponding to the metric λ{tilde over (α)}(n) with the largest magnitude (the symbol {tilde over (α)}n is the N1-th symbol in the length-N observation, which is not necessarily the middle symbol). Since φ(t) is assumed to be constant over the N-symbol observation interval, the magnitude of the metric λ{tilde over (α)}(n) is statistically independent of the channel pulse.
There are two difficulties with the receiver described by Equation (8). The first difficulty is the number of states grows exponentially with the observation interval N. The second difficulty is that, depending on the particular continuous phase modulation scheme, a large value for N is capable of being required to achieve adequate performance.
According to the present invention, the preceding difficulties are addressed by the receiver described the recursive metric:
wherein the forget factor α is in the range 0≦α≦1. The term {circumflex over (θ)}n−L(i) represents the phase contribution of all previous symbol decisions {circumflex over (α)}k(i) for the i-th state in the trellis. Each state in the trellis stores two values: a cumulative metric λ{tilde over (α)}(n−1), and a cumulative phase {circumflex over (θ)}n−L(i). The receiver uses a traceback matrix of length DD to output the symbol {circumflex over (α)}n−DD(i) corresponding to the state whose metric has the largest magnitude. Here, the branch metric λ{tilde over (α)}(n) is only a function of the L symbols being modulated by the phase pulse q(t), thus the number of states is ML−1. For the special case of α=1 this branch metric reduces to:
This identifies an important tradeoff. As α approaches unity, the branch metric in Equation (11) approaches the one in Equation (15). The metric in Equation (15) is a loose approximation to an infinitely long observation interval because it “remembers” previous observations through the use of a cumulative metric. The optimal maximum likelihood sequence estimating receiver also uses a cumulative metric to recursively compute a correlation from (∞,(n+1)T). The only difference is the non-coherent receiver cannot account for the phase states θn−L (shown in Equation (6)) in the trellis since the magnitude of the metrics (rather than the real part for the maximum likelihood sequence estimating receiver case) is used to determine survivors. However, when the slowly varying channel phase φ(t) is taken into account, the branch metric in Equation (15) will trace a curved path in the complex plane as φ(t) changes. This will reduce the magnitude of the metric and increase the probability that the competing paths through the trellis will have metrics with a magnitude larger than the true path. As α approaches zero, the branch metrics “forget” the infinite past more quickly and allow φ(t) to change more rapidly with less impact on the magnitude of the branch metrics.
Those of ordinary skill in the art will appreciate that the metric, described in Equation (11), is capable of being extended to more closely approximate an infinitely long observation interval. The reason for the inherently loose approximation in Equation (11) is that the trellis only allows for ML−1 states, when the underlying continuous phase modulation signal is described by pML−1 states, where the p-fold increase is due to the phase states θn−L. The extended metric for an observation interval of length N≧1 is given by:
It will be understood by those skilled in the art that an important difference between Equations (11)-(13) and Equations (16)-(18) is that N−1 symbols have been removed from the cumulative phase {circumflex over (θ)}n−L−N+1(i) to form {tilde over (θ)}n−L, which is associated with the branch metric. Thus, as paths merge and survivors are determined, more options are kept open in the trellis. The number of states in this trellis is ML−N−2.
As used hereinafter, the receiver defined in Equations (8)-(10) is denoted as “Receiver-A”, and the receiver defined in Equations (16)-(18) as “Receiver-B”. The skilled artisan will appreciate that Equations (11)-(13) define Receiver-B, wherein N=1. Both receivers have the parameter N, which is the multi-symbol observation length. Receiver-B is also parameterized by the forget factor α.
An alternate embodiment of Receiver-B is given by:
λ{tilde over (α)}(n)=λ{tilde over (α)}(n−1)+Re{e−j{circumflex over (θ)}
Q{tilde over (α)}(n)=αQ{tilde over (α)}(n−1)+(1−α)e−j{circumflex over (θ)}
The receiver defined in Equations (19)-(20) is denoted as “Receiver-C”. The skilled artisan will appreciate that Receiver-C is different from Receiver-B in that the cumulative metric λ{tilde over (α)}(n)is real-valued, and the noncoherent phase is resolved by the phase reference Q{tilde over (α)}(n). Those skilled in the art will understand that Receiver-C is similar to Receiver-B, such that Receiver-C is parameterized by the forget factor α and multi-symbol observation interval N. Receiver-C also uses the same variables, z{tilde over (α)}(n), {circumflex over (θ)}n−L−N+1(i), and {tilde over (θ)}n−L, as are found in Receiver-B. It will be apparent to the skilled artisan that due to the similarities between Receivers-B and -C, the performance results discussed below are given only for Receiver-B, but can be regarded as typical for Receiver-C.
The first continuous phase modulation scheme considered is the PCM/FM waveform, which is M=2,h=7/10,2RC, illustrated as
The next continuous phase modulation scheme in the simulations is the Advanced Range Telemetry Tier II waveform, which is M=4,h=7/10, {4/16, 5/16}, 3RC.
Up to this point, consideration has only been given to the performance with respect to the case of perfect symbol timing and carrier phase. Since the motivation for a noncoherent receiver is the case where the carrier phase is not known and assumed to be varying, a simple model will be introduced for variations in the carrier phase. Let
φn=φ(nT)=φn−1+vnmod2π (19)
where {vn} are independently and identically distributed Gaussian random variables with zero mean and variance δ2. This models the phase noise as a first order Markov process with Gaussian transition probability distribution. For perfect carrier phase tracking, δ=0.
The invention extends to computer programs in the form of source code, object code, partially compiled or otherwise, and code intermediate sources, or in any other form suitable for use in the implementation of the invention. Computer programs are suitably standalone applications, software components, scripts or plug-ins to other applications. Computer programs embedding the invention are advantageously embodied on a carrier, being any entity or device capable of carrying the computer program: for example, a storage medium such as ROM or RAM, optical recording media such as CD-ROM or magnetic recording media such as floppy discs. The carrier is any transmissible carrier such as an electrical, electromagnetic, or optical signal conveyed by electrical or optical cable, or by radio or other means. Computer programs are suitably downloaded across the Internet from a server. Computer programs are also capable of being embedded in an integrated circuit. Any and all such embodiments containing code that will cause a computer to perform substantially the invention principles as described, will fall within the scope of the invention.
The foregoing description of a preferred embodiment of the invention has been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise form disclosed. Obvious modifications or variations are possible in light of the above teachings. The embodiment was chosen and described to provide the best illustration of the principles of the invention and its practical application to thereby enable one of ordinary skill in the art to use the invention in various embodiments and with various modifications as are suited to the particular use contemplated. All such modifications and variations are within the scope of the invention as determined by the appended claims when interpreted in accordance with the breadth to which they are fairly, legally and equitably entitled.
This application is a continuation-in-part of U.S. patent application Ser. No. 11/252,108, filed Oct. 17, 2005, which claims priority to U.S. Provisional Patent Application Ser. No. 60/619,101, filed Oct. 15, 2004.
Number | Date | Country | |
---|---|---|---|
60619101 | Oct 2004 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 11252108 | Oct 2005 | US |
Child | 11369627 | Mar 2006 | US |