EP1566033B1 - Method and system for applying viterbi type psk demodulation for optimum correlation of gps signals - Google Patents
Method and system for applying viterbi type psk demodulation for optimum correlation of gps signals Download PDFInfo
- Publication number
- EP1566033B1 EP1566033B1 EP02779087A EP02779087A EP1566033B1 EP 1566033 B1 EP1566033 B1 EP 1566033B1 EP 02779087 A EP02779087 A EP 02779087A EP 02779087 A EP02779087 A EP 02779087A EP 1566033 B1 EP1566033 B1 EP 1566033B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- signal
- phase
- receiver
- gps
- creating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7097—Interference-related aspects
- H04B1/7103—Interference-related aspects the interference being multiple access interference
- H04B1/7105—Joint detection techniques, e.g. linear detectors
- H04B1/71057—Joint detection techniques, e.g. linear detectors using maximum-likelihood sequence estimation [MLSE]
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03178—Arrangements involving sequence estimation techniques
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/18—Phase-modulated carrier systems, i.e. using phase-shift keying
- H04L27/22—Demodulator circuits; Receiver circuits
- H04L27/233—Demodulator circuits; Receiver circuits using non-coherent demodulation
- H04L27/2332—Demodulator circuits; Receiver circuits using non-coherent demodulation using a non-coherent carrier
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
- H04L2027/003—Correction of carrier offset at baseband only
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0024—Carrier regulation at the receiver end
- H04L2027/0026—Correction of carrier offset
- H04L2027/0038—Correction of carrier offset using an equaliser
- H04L2027/0042—Correction of carrier offset using an equaliser the equaliser providing the offset correction per se
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0044—Control loops for carrier regulation
- H04L2027/0046—Open loops
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0044—Control loops for carrier regulation
- H04L2027/0063—Elements of loops
- H04L2027/0067—Phase error detectors
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0083—Signalling arrangements
- H04L2027/0085—Signalling arrangements with no special signals for synchronisation
Definitions
- the present invention relates to a method and system for increasing the sensitivity of the demodulation of direct-sequence spread spectrum modulated-signals and more particularly to the application of a Viterbi demodulation type algorithm to the output of a coherent correlation of a GPS C/A LI signal to improve the demodulation sensitivity.
- GPS global positioning system
- the GPS receiver in the handset unit In order for mobile location to be successfully based on AGPS, it is necessary that the GPS receiver in the handset unit have significantly higher sensitivity than is nominally required in a typical stand alone GPS unit. The reason is that the handset unit is usually oriented in positions or in locations that are not favorable for GPS signal reception. Hence a signal to noise ratio (SNR) sensitivity target is typically around 17 dB-Hz.
- SNR signal to noise ratio
- C/A L1 band which is typically used for mobile position location has a bandwidth of about 2 MHz, this implies that the input raw GPS, prior to processing, is more than 40 dB below the thermal noise floor.
- the C/A signal from a GPS SV space vehicle
- a DS-SS direct sequence spread spectrum
- BPSK binary phase shift keying
- the GPS receiver correlates the received signal with a locally generated DS-SS code signal.
- AGPS scheme traditional DS-SS correlation is also done.
- BS host base station
- the mobile GPS receiver is still required to determine the code delay and doppler to a finer resolution than that available from the host BS such that mobile location is possible.
- the optimum receiver would correlate the signal in a coherent fashion over an integration time period that is sufficiently long to provide about 11 dB SNR at the correlator output. This will typically provide an adequate probability of detection with a reasonable false alarm rate.
- the coherent integration epoch needs to be limited. Also valid GPS readings are still required even if the user does not hold the receiver steady. Hence, typically the coherent integration time is limited to 10 msec or less. As the available coherent integration epoch is not sufficient to obtain the sensitivity required, non-coherent summations of sequential coherent correlation integration outputs are used. However, non-coherent processing is a very inefficient means of further enhancing the SNR of a signal as it discards certain known statistical aspects of the signal.
- the method and system of the present invention overcome the above by providing a way to demodulate the RF phase trajectory of a DS-SS correlation, using quasi-coherent integration for a measurement epoch of arbitrary duration.
- the receiver in the present invention downshifts an incoming signal using a local oscillator.
- the receiver then despreads the signal.
- the output from the despreader would be a DC signal.
- the phase trajectory is generally a phasor.
- a trellis is created with a fixed number of phase states evenly distributed between zero and 360°.
- the state transitions in the trellis are limited based on known phase slew rate limitations. Probabilities for these state transitions are similarly based on these same factors.
- a Viterbi algorithm is then applied to the phase of the despread signal. Likelihood metrics are created for each path based on the phase of the despread signal as well as path probabilities. The system and method then uses the Viterbi algorithm to find the highest likelihood path for the phase trajectory.
- the present invention therefore provides a of improving sensitivity in the demodulation of a received signal by a receiver over an arbitrary measurement time epoch, said method comprising the steps of: correlating said received signal with a local replica of pseudo noise code in a coherent fashion over time intervals in said time epoch creating a correlation signal; creating a trellis of evenly distributed phase state nodes at each time interval, said creating step comprising: defining a plurality of phase states representing the phases evenly quantized over 0 to 360 degrees; defining possible state transitions from and to each phase state node; creating paths between phase state nodes in one time interval and phase state nodes in another time interval according to said possible state transitions; assigning a transition probability to each path; and creating a likelihood metric for each path based on a measured phase of the correlation signal and the transition probability for the path, said measured phase of the correlation signal having a random process which is for creation of the likelihood metric approximated utilizing a first order Markov process; and utilizing a Viterbi algorithm on said
- the present invention still further provides a receiver for receiving a direct sequence spread spectrum signal, said receiver comprising: an antenna for receiving the direct sequence spread spectrum signal; a downconverter for downconverting the received signal, producing a downconverted signal; an analog to digital converter to convert the downconverted signal to a digital signal; a despreader for despreading and coherently correlating the digital signal to a known signal, creating a despread signal; and a processor for applying a Viterbi algorithm to a trellis created for the despread signal, the processor being arranged to: break the despread signal into time intervals; create the trellis of evenly distributed phase state nodes at each time interval by: defining a plurality of phase states representing the phases evenly quantized over 0 to 360 degrees; defining possible state transitions from and to each phase state node; creating paths between phase states node in one time interval and phase state nodes in another time interval according to said possible state transitions; assigning a transition probability to each path; and creating a likelihood metric for
- the present invention still further provides a wireless mobile handset comprising a receiver in accordance with the present invention.
- a system and method are described that provide a practical way of performing an optimal correlation of the DS-SS GPS signal using the known statistical properties of the signal.
- the method makes use of the generally known slew rate limitation of the LO instability and SV doppler.
- the main aspect of the method is a Maximum Likelihood Viterbi procedure similar to the processing used in trellis decoding of standard PSK communication signals.
- the method is applied exclusively to GPS signals in the following disclosure, one skilled in the art will realize that the method is also applicable for precision tracking of CDMA pilot signals for applications of high sensitivity mobile position location.
- the antenna 10 of a receiver 1 receives a direct sequence spread spectrum signal.
- the signal is then passed through band-pass filter 20 to isolate the signal.
- band-pass filter 20 is centred at 1.574 Ghz, which is the frequency of the GPS C/A L1 carrier.
- band-pass filter 20 The output from band-pass filter 20 is then downconverted by multiplying it with the signal from the radio frequency local oscillator 30 in multiplier 40, producing a baseband signal. This baseband signal is then passed through low-pass filter 50 and an analog to digital converter 60 to create a spread spectrum digitized signal.
- the spread spectrum digitized signal is then passed through a digital signal processor 70.
- Digital signal processor 70 is comprised of a despreader 80 and a processor 90.
- Despreader 80 is comprised of multiplier 84 and signal generator 83 which generates coarse acquisition LI signal 82 which is multiplied with the spread spectrum digitized signal.
- despreader 80 In a perfect system, the output from despreader 80 would be a DC signal. However, due to frequency and phase instabilities in local oscillator 30, and due to uncertainty of the GPS SV doppler, the phase trajectory is generally a time varying phasor.
- the output of despreader 80 is thus processed by processor 90, which applies a Viterbi algorithm as described below.
- the trellis of Figure 2 is only one example.
- One skilled in the art will realize that other trellis state transition probabilities are possible, and that a probability trellis can be created to more accurately model these other probabilities.
- Figure 3 a system is illustrated in which the slew rate can be higher but a lower probability is assigned to higher slew rates.
- M The number of phase states, "M" is a function of the phase resolution required. If M is too small then there will be a small demodulation loss due to the discrete phase quantization.
- the standard Viterbi algorithm is used to eliminate the trellis paths of lower likelihood.
- a finite number of paths converge onto each of the M trellis phase states.
- the path selected for each path is the one associated with the highest likelihood. Consequently, at each interval step, P paths converge to each trellis state and subsequently P new path segments are generated, as shown in the above examples for the (k+1) th trellis interval.
- conditional probability metric can be approximated as (p k -S m ) 2 where S m represents the assigned phase of the m th state.
- S m represents the assigned phase of the m th state.
- S m ) is irrelevant as it is added to each candidate trellis path and therefore does not result in a path differentiating factor.
- Practical simulation experience has shown that this simple path likelihood metric is effective even though the approximation to the conditional probability of P(p k
- Another very important factor is that the differentiating part of the likelihood metric when the above approximation is used becomes independent of the variance of the conditional probability of p k . This is a key factor in the robustness and practicality of the method.
- the Viterbi demodulator was implemented as described above and tracked the phase excursion of the correlation samples x k .
- the accuracy of the demodulated phase was compared to the actual phase to determine a signal loss factor that was dependent on phase errors only.
- Figure 5 shows a plot of this loss factor as a function of the SNR of the correlation samples. As expected, the loss generally decreases as a function of increased SNR. At very low SNR levels, the phase tracking becomes so erratic that the Viterbi demodulation is not effective.
- Figure 5 can be used to determine the number of samples required for a given SNR.
- the samples at 0 dB SNR As stated earlier, to satisfy the requirements of adequate detection and false alarm statistics, approximately 11 dB of SNR is required at the output of the overall correlation process.
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Error Detection And Correction (AREA)
Abstract
Description
- The present invention relates to a method and system for increasing the sensitivity of the demodulation of direct-sequence spread spectrum modulated-signals and more particularly to the application of a Viterbi demodulation type algorithm to the output of a coherent correlation of a GPS C/A LI signal to improve the demodulation sensitivity.
- Due to Federal Communications Commission (FCC) requirements for identifying mobile phone location, vendors are beginning to incorporate global positioning system (GPS) receivers into their handsets. This is an assisted form of GPS where the host base station provides key search parameter information to the handset such that the GPS satellite search and subsequent delay measurement becomes significantly less computationally extensive for the handset. In the code-division multiple access (CDMA) context, the messaging between the base station and the handset to support the assisted GPS (AGPS) is documented in "Position Determination Service Standards for Dual Mode Spread Spectrum Systems," TIA/EIA/IS-801-1, published by the Telecommunications Industry Association (TIA).
- In order for mobile location to be successfully based on AGPS, it is necessary that the GPS receiver in the handset unit have significantly higher sensitivity than is nominally required in a typical stand alone GPS unit. The reason is that the handset unit is usually oriented in positions or in locations that are not favorable for GPS signal reception. Hence a signal to noise ratio (SNR) sensitivity target is typically around 17 dB-Hz. Given that the GPS Coarse/Acquisition (C/A) L1 band which is typically used for mobile position location has a bandwidth of about 2 MHz, this implies that the input raw GPS, prior to processing, is more than 40 dB below the thermal noise floor.
- There is also a significant push to drive the cost of the overall handset down which limits the quantity of the processing which can be done on the received GPS signals to extract them from the noise. This also limits the quality of the various radio frequency (RF) oscillator components used in the basket's receiver. Consequently, it can be expected that there will be some additional local oscillator (LO) noise that needs to be compensated for. Notably, there is no closed automatic frequency correction (AFC) loop possible during the GPS measurement as the input SNR is too low. The lack of AFC, and the use of low cost, low power, RF LO components, implies that there will inevitably be a significant finite offset, drift, and some instability associated with the LO down conversion frequency.
- The C/A signal from a GPS SV (space vehicle) is a DS-SS (direct sequence spread spectrum) with a chip rate of 1.023 Mcps. It is modulated as binary phase shift keying (BPSK) on a 1.574 GHz carrier. The GPS receiver correlates the received signal with a locally generated DS-SS code signal. In the AGPS scheme, traditional DS-SS correlation is also done. However, detailed information regarding the doppler shift of the SV GPS signal and code offset is available from the host base station (BS) which significantly reduces the search and detection effort. Nevertheless, the mobile GPS receiver is still required to determine the code delay and doppler to a finer resolution than that available from the host BS such that mobile location is possible.
- Based on standard assumptions regarding the noise in the GPS signal channel, the optimum receiver would correlate the signal in a coherent fashion over an integration time period that is sufficiently long to provide about 11 dB SNR at the correlator output. This will typically provide an adequate probability of detection with a reasonable false alarm rate. However, in the case of a mobile GPS receiver, due to the instability of the RF LO and the uncertainly of the SV doppler, the coherent integration epoch needs to be limited. Also valid GPS readings are still required even if the user does not hold the receiver steady. Hence, typically the coherent integration time is limited to 10 msec or less. As the available coherent integration epoch is not sufficient to obtain the sensitivity required, non-coherent summations of sequential coherent correlation integration outputs are used. However, non-coherent processing is a very inefficient means of further enhancing the SNR of a signal as it discards certain known statistical aspects of the signal.
- Various examples of the use of Viterbi algorithms for phase trajectory determinations are known. For example US Patent No. 6,477,208 to Huff teaches a method and apparatus for processing a received digitally-modulated carrier signal to coherently demodulate a signal utilizing a composite trellis diagram. However, the method of Huff is only applicable to cases of discrete jumps in phase based on specific modulations. When using a continuous, free running oscillator, the phase steps are continuous in time. Huff does not teach a method for approximating these phase steps.
- The method and system of the present invention overcome the above by providing a way to demodulate the RF phase trajectory of a DS-SS correlation, using quasi-coherent integration for a measurement epoch of arbitrary duration. The receiver in the present invention downshifts an incoming signal using a local oscillator. The receiver then despreads the signal.
- In a perfect system the output from the despreader would be a DC signal. However, due to frequency and phase instabilities in the radio frequency local oscillator, and due to uncertainty of the GPS SV doppler, the phase trajectory is generally a phasor.
- In order to provide increased sensitivity, a trellis is created with a fixed number of phase states evenly distributed between zero and 360°. The state transitions in the trellis are limited based on known phase slew rate limitations. Probabilities for these state transitions are similarly based on these same factors.
- A Viterbi algorithm is then applied to the phase of the despread signal. Likelihood metrics are created for each path based on the phase of the despread signal as well as path probabilities. The system and method then uses the Viterbi algorithm to find the highest likelihood path for the phase trajectory.
- The present invention therefore provides a of improving sensitivity in the demodulation of a received signal by a receiver over an arbitrary measurement time epoch, said method comprising the steps of: correlating said received signal with a local replica of pseudo noise code in a coherent fashion over time intervals in said time epoch creating a correlation signal; creating a trellis of evenly distributed phase state nodes at each time interval, said creating step comprising: defining a plurality of phase states representing the phases evenly quantized over 0 to 360 degrees; defining possible state transitions from and to each phase state node; creating paths between phase state nodes in one time interval and phase state nodes in another time interval according to said possible state transitions; assigning a transition probability to each path; and creating a likelihood metric for each path based on a measured phase of the correlation signal and the transition probability for the path, said measured phase of the correlation signal having a random process which is for creation of the likelihood metric approximated utilizing a first order Markov process; and utilizing a Viterbi algorithm on said trellis to perform a maximum likelihood estimation of the phase trajectory of said correlation signal with said quantized resolution of phase states over 0 to 360° throughout the measurement time epoch.
- The present invention still further provides a receiver for receiving a direct sequence spread spectrum signal, said receiver comprising: an antenna for receiving the direct sequence spread spectrum signal; a downconverter for downconverting the received signal, producing a downconverted signal; an analog to digital converter to convert the downconverted signal to a digital signal; a despreader for despreading and coherently correlating the digital signal to a known signal, creating a despread signal; and a processor for applying a Viterbi algorithm to a trellis created for the despread signal, the processor being arranged to: break the despread signal into time intervals; create the trellis of evenly distributed phase state nodes at each time interval by: defining a plurality of phase states representing the phases evenly quantized over 0 to 360 degrees; defining possible state transitions from and to each phase state node; creating paths between phase states node in one time interval and phase state nodes in another time interval according to said possible state transitions; assigning a transition probability to each path; and creating a likelihood metric for each path based on a measured phase of said despread signal and the transition probability for the path, the measured phase of the despread signal having a random process which is for creation of the likelihood metric approximated utilizing a first order Markov process; and utilize the Viterbi algorithm on said trellis to perform a maximum likelihood estimation of the phase trajectory of said despread signal with said quantized resolution of phase states over 0 to 360°throughout the time interval.
- The present invention still further provides a wireless mobile handset comprising a receiver in accordance with the present invention.
- Preferred embodiments of the present invention will now be described in greater detail and will be better understood when read in conjunction with the following drawings in which:
- Figure 1 is a schematic diagram of the receiver of the present invention;
- Figure 2 is an example trellis diagram for use with the method of the present invention;
- Figure 3 is an alternative trellis diagram for use with the method of the present invention;
- Figure 4 is a graph illustrating an approximation of the probability distribution function of a phase sample as a periodic gaussian pulse on top of a constant function; and
- Figure 5 is a graph illustrating the losses in decibels verses the signal to noise ratio of the samples received in decibels.
- In this disclosure, a system and method are described that provide a practical way of performing an optimal correlation of the DS-SS GPS signal using the known statistical properties of the signal. In particular, the method makes use of the generally known slew rate limitation of the LO instability and SV doppler.
- The main aspect of the method is a Maximum Likelihood Viterbi procedure similar to the processing used in trellis decoding of standard PSK communication signals. Although the method is applied exclusively to GPS signals in the following disclosure, one skilled in the art will realize that the method is also applicable for precision tracking of CDMA pilot signals for applications of high sensitivity mobile position location.
- Reference is now made to Figure 1. The
antenna 10 of areceiver 1 receives a direct sequence spread spectrum signal. The signal is then passed through band-pass filter 20 to isolate the signal. In a preferred embodiment, band-pass filter 20 is centred at 1.574 Ghz, which is the frequency of the GPS C/A L1 carrier. - The output from band-
pass filter 20 is then downconverted by multiplying it with the signal from the radio frequencylocal oscillator 30 inmultiplier 40, producing a baseband signal. This baseband signal is then passed through low-pass filter 50 and an analog todigital converter 60 to create a spread spectrum digitized signal. - The spread spectrum digitized signal is then passed through a digital signal processor 70. Digital signal processor 70 is comprised of a
despreader 80 and a processor 90.Despreader 80 is comprised ofmultiplier 84 andsignal generator 83 which generates coarseacquisition LI signal 82 which is multiplied with the spread spectrum digitized signal. - In a perfect system, the output from
despreader 80 would be a DC signal. However, due to frequency and phase instabilities inlocal oscillator 30, and due to uncertainty of the GPS SV doppler, the phase trajectory is generally a time varying phasor. The output ofdespreader 80 is thus processed by processor 90, which applies a Viterbi algorithm as described below. - Consider the total GPS signal measurement time epoch, denoted by Ttot as being divided into K equal length segments of Tc each. It is assumed that the GPS correlator will coherently integrate for each of the segment lengths Tc resulting in K correlation samples X=[x1,x2,...,xk]T. It is also assumed that Tc is short enough such that uncertainties in GPS doppler, RF LO offset or frequency/phase instabilities of the RFLO are such that the resulting signal loss in the coherent correlation output samples in X is not significant. Typically Tc will be somewhere between 1 to 10 msec.
- X=[x1,x2,...,xK]T is mapped into corresponding sample phases denoted by P=[p1,p2,...,pK]T. These phases will be uniformly distributed between 0 and 360°. In order to implement a practical Viterbi algorithm, discrete phases are defined as being uniformly spaced such that:
- M = number of discrete phases
- s0,s1,...sM-1 are the discrete phase states given by sn=(m/M) 360°
- Let k denote the kth state interval corresponding to the kth sample time or xk. Reference is now made to Figure 2, which shows an example trellis diagram for better illustration of the following description.
- Consider the mth state of the kth interval. If there are no phase dynamics due to the RF LO or doppler, then it will map into state m of the (k+1)th interval with a probability of 1. However, due to the uncertainty in the phase dynamics, the mth state can map into state (m-1), m, (m+1) of the (k+1)th interval. Clearly, with these trellis state transition probabilities, the modeled slew rate of the phase dynamics can be accurately limited.
- The trellis of Figure 2 is only one example. One skilled in the art will realize that other trellis state transition probabilities are possible, and that a probability trellis can be created to more accurately model these other probabilities. For example, in Figure 3, a system is illustrated in which the slew rate can be higher but a lower probability is assigned to higher slew rates.
- Intuitively such slew rate limiting lattice structures as described above effectively "low pass filter" the phase samples to reduce the influence of noise in correlation samples X.
- The number of phase states, "M", is a function of the phase resolution required. If M is too small then there will be a small demodulation loss due to the discrete phase quantization.
- At each update, the standard Viterbi algorithm is used to eliminate the trellis paths of lower likelihood. (A good description of the Viterbi algorithm is found in "Digital Communications" by J. Proakis, McGraw Hill , 3rd edition 1995.) Briefly, at the kth update a finite number of paths converge onto each of the M trellis phase states. For example, in the trellis structure of Figure 2, there are three candidate paths and in the trellis structure of Figure 3, there are five such paths. Each of these paths will have a likelihood metric associated with it. The path selected for each path is the one associated with the highest likelihood. Consequently, at each interval step, P paths converge to each trellis state and subsequently P new path segments are generated, as shown in the above examples for the (k+1)th trellis interval.
-
- This is an example only as there are a variety of numerical methods of extracting the phase from xk. In order to update the likelihood metric for each surviving path at the kth trellis interval, an estimate of probability of the received phase being pk conditioned on the phase state and the trellis path corresponding to the particular candidate path being considered is required. As the SNR of the correlator outputs, xk, is typically low (usually less than 0 dB), the probability distribution function (PDF) of the phase sample pk is not gaussian. Also the phase trajectory is dependent in a complicated statistical way on the candidate trellis path being considered. Consequently, finding the actual conditional probability is impractical. However, there are several assumptions and approximations that can be reasonably made that significantly simplify the likelihood probability metric. These are:
- 1. Model the phase random process as first order Markov
- 2. Approximate the PDF of pk as a periodic gaussian pulse on top of a constant function as illustrated in Figure 4.
- Given this approach, the conditional probability metric can be approximated as (pk-Sm)2 where Sm represents the assigned phase of the mth state. The influence of the constant portion of P(pk|Sm) is irrelevant as it is added to each candidate trellis path and therefore does not result in a path differentiating factor. Practical simulation experience has shown that this simple path likelihood metric is effective even though the approximation to the conditional probability of P(pk|Sm) is crude.
- Another very important factor is that the differentiating part of the likelihood metric when the above approximation is used becomes independent of the variance of the conditional probability of pk. This is a key factor in the robustness and practicality of the method.
- Below there is provided a simple MatLab simulation to demonstrate the effectiveness of the present method.
- In the simulation, an array of 2000 complex samples of xk were generated based on the assumption that the signal component was a constant amplitude exponential with a phase increment of 0.03 radians. Random noise was added to the samples of xk that was zero mean gaussian to both in-phase and quadrature phase components. The noise components were independent from sample to sample and also the noise in the quadrature and in-phase components were independent of each other. The variance of the noise components was selected to satisfy a given SNR of the correlation samples xn for the simulation.
- The Viterbi demodulator was implemented as described above and tracked the phase excursion of the correlation samples xk. The accuracy of the demodulated phase was compared to the actual phase to determine a signal loss factor that was dependent on phase errors only.
- Figure 5 shows a plot of this loss factor as a function of the SNR of the correlation samples. As expected, the loss generally decreases as a function of increased SNR. At very low SNR levels, the phase tracking becomes so erratic that the Viterbi demodulation is not effective.
- Figure 5 can be used to determine the number of samples required for a given SNR. Consider the samples at 0 dB SNR. As stated earlier, to satisfy the requirements of adequate detection and false alarm statistics, approximately 11 dB of SNR is required at the output of the overall correlation process. Hence if the samples xk of 0 dB SNR are generated after coherent correlation of Tc, the number of samples required, K, is given by the formula:
or, using Figure 5, - Solving this, K=16.
- As indicated, the additional 1 dB added to the right hand side of the above formula is to offset the coherence integration losses represented in Figure 5. Comparing this with the traditional scheme which uses non-coherent summation of the correlator output samples, a minimum of K=26 is required. Hence, for a given GPS hypothesis dwell time, sensitivity can be significantly better with the Viterbi style demodulator as represented in this disclosure.
- The above-described embodiments of the present invention are meant to be illustrative of preferred embodiments and are not intended to limit the scope of the present invention. Also, various modifications, which would be readily apparent to one skilled in the art, are intended to be within the scope of the present invention. The only limitations to the scope of the present invention are set forth in the following claims appended hereto.
Claims (26)
- A method of improving sensitivity in the demodulation of a received signal by a receiver (1) over an arbitrary measurement time epoch, said method comprising the steps of:correlating (84) said received signal with a local replica of pseudo noise code (82) in a coherent fashion over time intervals in said time epoch creating a correlation signal;creating a trellis of evenly distributed phase state nodes at each time interval, said creating step comprising:defining a plurality of phase states representing the phases evenly quantized over 0 to 360 degrees;defining possible state transitions from and to each phase state node;creating paths between phase state nodes in one time interval and phase state nodes in another time interval according to said possible state transitions;assigning a transition probability to each path; andcreating a likelihood metric for each path based on a measured phase of the correlation signal and the transition probability for the path, said measured phase of the correlation signal having a random process which is for creation of the likelihood metric approximated utilizing a first order Markov process; andutilizing a Viterbi algorithm (90) on said trellis to perform a maximum likelihood estimation of the phase trajectory of said correlation signal with said quantized resolution of phase states over 0 to 360°throughout the measurement time epoch.
- The method of claim 1, wherein the possible state transitions and the probability of the paths are assigned to reflect properties of said receiver.
- The method of claim 2 wherein the step of creating possible state transitions for each node is performed based on a known phase slew rate limitation of said receiver (1).
- The method of claim 3, wherein the known phase slew rate limitation is calculated from the instability of a radio frequency local oscillator in said receiver (1).
- The method of any one of claims 1 to 4, wherein the received signal is a direct sequence spread spectrum signal.
- The method of any one of claims 1 to 4, wherein the received signal is a global positioning system (GPS) coarse/acquisition LI signal generated by a space vehicle (SV).
- The method of any one of claims 1 to 4 wherein the received signal is a code-division multiple access (CDMA) pilot signal.
- The method of claim 6, wherein the GPS SV creates a doppler shift in the phase trajectory, and the known phase slew rate limitation is calculated from the uncertainty of the GPS SV doppler shift.
- The method of claim 6 or claim 8, wherein the known phase slew rate limitation is calculated from both the instability of a radio frequency local oscillator (30) in said receiver (1) and the uncertainty of the GPS SV doppler.
- The method of any one of claims 1 to 9, wherein the likelihood metric is created based on an approximation of a probability distribution function of the phase of said correlation signal.
- The method of claim 10, wherein the approximation is to model the probability distribution function of the phase as a periodic gaussian pulse on top of a constant function.
- The method of any one of claims 1 to 11 wherein said receiver (1) is a mobile receiver.
- A receiver (1) for receiving a direct sequence spread spectrum signal, said receiver comprising:an antenna (10) for receiving the direct sequence spread spectrum signal;a downconverter (40) for downconverting the received signal, producing a downconverted signal;an analog to digital converter (60) to convert the downconverted signal to a digital signal;a despreader (80) for despreading and coherently correlating the digital signal to a known signal (82), creating a despread signal; anda processor (90) for applying a Viterbi algorithm to a trellis created for the despread signal, the processor (90) being arranged to:break the despread signal into time intervals;create the trellis of evenly distributed phase state nodes at each time interval by:defining a plurality of phase states representing the phases evenly quantized over 0 to 360 degrees;defining possible state transitions from and to each phase state node;creating paths between phase states node in one time interval and phase state nodes in another time interval according to said possible state transitions;assigning a transition probability to each path; andcreating a likelihood metric for each path based on a measured phase of said despread signal and the transition probability for the path, the measured phase of the despread signal having a random process which is for creation of the likelihood metric approximated utilizing a first order Markov process; andutilize the Viterbi algorithm on said trellis to perform a maximum likelihood estimation of the phase trajectory of said despread signal with said quantized resolution of phase states over 0 to 360° throughout the time interval.
- The receiver (1) of claim 13, wherein the received signal is a global positioning system coarse/acquisition (C/A) LI signal generated by a space vehicle (SV).
- The receiver (1) of claim 13 or claim 14, wherein the received signal is a CDMA pilot signal.
- The receiver (1) of any one of claims 13 to 16, wherein the receiver (1) is a mobile receiver.
- The receiver (1) of any one of claims 13 or 15, wherein the received signal is a GPS C/A L1 signal.
- The receiver (1) of claim 17, wherein the possible state transitions and the probability of the paths are assigned to reflect properties of said receiver.
- The receiver (1) of claim 18, wherein the possible state transitions for each node are based on a known phase slew rate limitation of said receiver.
- The receiver (1) of claim 19, wherein the known phase slew rate limitation is calculated from the instability of a radio frequency local oscillator in said receiver.
- The receiver (1) of claim 14, wherein the GPS SV creates a doppler shift in a phase trajectory of said received signal, and the known phase slew rate limitation is calculated from the uncertainty of said GPS SV doppler shift.
- The receiver (1) of claim 14 or claim 21, wherein the possible state transitions for each node are based on a known phase slew rate limitation of said receiver.
- The receiver (1) of claim 22, wherein the known phase slew rate limitation is calculated from both the instability of a radio frequency local oscillator (30) in the receiver and the uncertainty of the GPS SV doppler.
- The receiver (1) of any one of claims 13 to 23, wherein the likelihood metric is created based on an approximation of a probability distribution function of the phase of said despread signal.
- The receiver (1) of claim 24, wherein said approximation is to model the probability distribution function of the phase as a periodic gaussian pulse on top of a constant function.
- A wireless mobile handset comprising the receiver (1) of any one of claims 13 to 25.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AT02779087T ATE343888T1 (en) | 2002-11-27 | 2002-11-27 | METHOD AND SYSTEM FOR APPLYING VITERBI-TYPE PSK DEMODULATION FOR OPTIMUM CORRELATION OF GPS SIGNALS |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CA2002/001843 WO2004049663A1 (en) | 2002-11-27 | 2002-11-27 | Method and system for applying viterbi type psk demodulation for optimum correlation of gps signals |
Publications (2)
Publication Number | Publication Date |
---|---|
EP1566033A1 EP1566033A1 (en) | 2005-08-24 |
EP1566033B1 true EP1566033B1 (en) | 2006-10-25 |
Family
ID=32331600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP02779087A Expired - Lifetime EP1566033B1 (en) | 2002-11-27 | 2002-11-27 | Method and system for applying viterbi type psk demodulation for optimum correlation of gps signals |
Country Status (7)
Country | Link |
---|---|
US (2) | US7599452B2 (en) |
EP (1) | EP1566033B1 (en) |
AU (1) | AU2002342489A1 (en) |
CA (1) | CA2504648C (en) |
DE (1) | DE60215707T2 (en) |
HK (1) | HK1080236A1 (en) |
WO (1) | WO2004049663A1 (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2002342489A1 (en) * | 2002-11-27 | 2004-06-18 | Research In Motion Limited | Method and system for applying viterbi type psk demodulation for optimum correlation of gps signals |
US7642959B2 (en) * | 2003-03-21 | 2010-01-05 | Global Locate, Inc. | Method and apparatus for processing satellite signals at a satellite positioning system receiver |
US8010124B2 (en) * | 2003-03-24 | 2011-08-30 | Sony Ericsson Mobile Communications Ab | Methods, systems and computer program products for providing location determination information to an assisted location service |
EP2279484B1 (en) * | 2008-05-19 | 2011-10-12 | Nxp B.V. | Rfid transponder |
US20110050496A1 (en) * | 2009-09-02 | 2011-03-03 | Qualcomm Incorporated | Energy Domain Based Peak Reconstruction Methods And Apparatuses |
US8630362B1 (en) | 2011-05-02 | 2014-01-14 | Urbain A. von der Embse | QLM co-state MAP trellis |
KR101343531B1 (en) | 2012-11-12 | 2013-12-20 | 국방과학연구소 | Method of recognizing pri modulation type based on hidden markov model |
US8917786B1 (en) * | 2013-05-09 | 2014-12-23 | Urbain Alfred von der Embse | QLM communications faster than Shannon rate |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5353278A (en) * | 1990-05-08 | 1994-10-04 | Sharp Kabushiki Kaisha | Magneto-optical disk and method of manufacturing the same |
US5375129A (en) * | 1990-07-19 | 1994-12-20 | Technophone Limited | Maximum likelihood sequence detector |
US5237587A (en) * | 1992-11-20 | 1993-08-17 | Magnavox Electronic Systems Company | Pseudo-noise modem and related digital correlation method |
FR2709028B1 (en) * | 1993-08-13 | 1995-10-20 | Matra Communication | Method for selecting the propagation paths used to receive messages transmitted by CDMA radiocommunication. |
DE69532577T2 (en) * | 1994-08-08 | 2004-12-16 | Ntt Mobile Communications Network Inc. | METHOD FOR DELAY DEMODULATION OF DPSK SHAFTS WITH LINEAR PREDICTION |
US5832041A (en) * | 1994-10-21 | 1998-11-03 | Philips Electronics North America Corporation | 64 QAM signal constellation which is robust in the presence of phase noise and has decoding complexity |
US5729558A (en) * | 1995-03-08 | 1998-03-17 | Lucent Technologies Inc. | Method of compensating for Doppler error in a wireless communications system, such as for GSM and IS54 |
US5825887A (en) * | 1995-12-28 | 1998-10-20 | Trimble Navigation Limited | Transmitting and receiving apparatus for full code correlation operation under encryption for satellite positioning system |
US6356595B1 (en) * | 1997-10-14 | 2002-03-12 | Sony Corporation | Method and apparatus for decoding continuously coded convolutionally encoded messages |
US6477208B1 (en) * | 1997-10-30 | 2002-11-05 | Comtier | Composite trellis system and method |
US5974091A (en) * | 1997-10-30 | 1999-10-26 | Communication Network Systems | Composite trellis system and method |
US6370201B1 (en) * | 1997-11-04 | 2002-04-09 | L-3 Communications Corp. | Simplified branch metric calculation in pragmatic trellis decoders |
US6577271B1 (en) * | 1999-03-30 | 2003-06-10 | Sirf Technology, Inc | Signal detector employing coherent integration |
US6633615B1 (en) * | 2000-01-31 | 2003-10-14 | Agere Systems Inc. | Trellis transition-probability calculation with threshold normalization |
US6721366B1 (en) * | 2000-04-19 | 2004-04-13 | Lockheed Martin Corporation | Phase tracking apparatus and method for continuous phase modulated signals |
JP4076202B2 (en) * | 2000-08-07 | 2008-04-16 | 富士通株式会社 | Spread spectrum signal receiver and receiving method |
US6940927B2 (en) * | 2001-07-17 | 2005-09-06 | Lockheed Martin Corporation | Simplified symbol timing tracking circuit for a CPM modulated signal |
US7120207B2 (en) * | 2001-12-31 | 2006-10-10 | Nokia Corporation | Transmission method and radio receiver |
US7154965B2 (en) * | 2002-10-08 | 2006-12-26 | President And Fellows Of Harvard College | Soft detection of data symbols in the presence of intersymbol interference and timing error |
AU2002342489A1 (en) * | 2002-11-27 | 2004-06-18 | Research In Motion Limited | Method and system for applying viterbi type psk demodulation for optimum correlation of gps signals |
-
2002
- 2002-11-27 AU AU2002342489A patent/AU2002342489A1/en not_active Abandoned
- 2002-11-27 WO PCT/CA2002/001843 patent/WO2004049663A1/en not_active Application Discontinuation
- 2002-11-27 CA CA002504648A patent/CA2504648C/en not_active Expired - Lifetime
- 2002-11-27 US US10/534,735 patent/US7599452B2/en active Active
- 2002-11-27 EP EP02779087A patent/EP1566033B1/en not_active Expired - Lifetime
- 2002-11-27 DE DE60215707T patent/DE60215707T2/en not_active Expired - Lifetime
-
2006
- 2006-02-21 HK HK06102296A patent/HK1080236A1/en not_active IP Right Cessation
-
2009
- 2009-09-03 US US12/553,285 patent/US8457255B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
HK1080236A1 (en) | 2006-04-21 |
DE60215707T2 (en) | 2007-02-08 |
AU2002342489A1 (en) | 2004-06-18 |
US8457255B2 (en) | 2013-06-04 |
CA2504648A1 (en) | 2004-06-10 |
EP1566033A1 (en) | 2005-08-24 |
CA2504648C (en) | 2008-08-12 |
US20060034350A1 (en) | 2006-02-16 |
WO2004049663A1 (en) | 2004-06-10 |
DE60215707D1 (en) | 2006-12-07 |
US20100027594A1 (en) | 2010-02-04 |
US7599452B2 (en) | 2009-10-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8457255B2 (en) | Method and system for applying Viterbi type PSK demodulation for optimum correlation of GPS signals | |
US7623562B2 (en) | Initial synchronization acquiring device and method for parallel processed DS-CDMA UWB system and DS-CDMA system's receiver using the same | |
US6201954B1 (en) | Method and system for providing an estimate of the signal strength of a received signal | |
US8208513B2 (en) | Spread-spectrum receiver and reception method | |
JP4425467B2 (en) | Pilot strength measurement and multipath delay searcher for CDMA receiver | |
US8798214B2 (en) | Minimum finger low-power demodulator for wireless communication | |
US7020180B2 (en) | Method and apparatus for acquiring pilots over code space and frequency errors in a CDMA communication system | |
US6134260A (en) | Method and apparatus for frequency acquisition and tracking for DS-SS CDMA receivers | |
US7852905B2 (en) | System and method for despreading in a spread spectrum matched filter | |
EP0892528A2 (en) | Carrier recovery for DSSS signals | |
US6950458B2 (en) | Method and apparatus for acquiring slot timing and frequency offset correction and storage medium storing control program therefor | |
US20060072650A1 (en) | Method and system for managing, controlling, and combining signals in a frequency selective multipath fading channel | |
US20120033768A1 (en) | Non-coherent detection apparatus and method for ieee 802.15.4 lr-wpan bpsk receiver | |
WO2000014893A1 (en) | User terminal parallel searcher | |
US7106783B2 (en) | Method and apparatus for searching multipaths of mobile communication system | |
US20070253471A1 (en) | Gps Receiver Using Differential Correlation | |
EP1260032B1 (en) | Method and apparatus for fast detection of locations of path rays in a rake receiver | |
Persson | CDMA code synchronization using segmented matched filter with accumulation and best match selection | |
US7711038B1 (en) | System and method for despreading in a spread spectrum matched filter | |
Casandra et al. | Performance Evaluation of a Tracking Algorithm for Galileo E1 Signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 20050530 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LI LU MC NL PT SE SK TR |
|
AX | Request for extension of the european patent |
Extension state: AL LT LV MK RO SI |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1080236 Country of ref document: HK |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR IE IT LI LU MC NL PT SE SK TR |
|
AX | Request for extension of the european patent |
Extension state: AL LT LV MK RO SI |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: BE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: AT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: LI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: CH Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 Ref country code: IT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT;WARNING: LAPSES OF ITALIAN PATENTS WITH EFFECTIVE DATE BEFORE 2007 MAY HAVE OCCURRED AT ANY TIME BEFORE 2007. THE CORRECT EFFECTIVE DATE MAY BE DIFFERENT FROM THE ONE RECORDED. Effective date: 20061025 |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20061127 |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20061130 |
|
REF | Corresponds to: |
Ref document number: 60215707 Country of ref document: DE Date of ref document: 20061207 Kind code of ref document: P |
|
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1080236 Country of ref document: HK |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070125 Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070125 Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070125 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: ES Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070205 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: PT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070326 |
|
ET | Fr: translation filed | ||
NLV1 | Nl: lapsed or annulled due to failure to fulfill the requirements of art. 29p and 29m of the patents act | ||
LTIE | Lt: invalidation of european patent or patent extension |
Effective date: 20061025 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
26N | No opposition filed |
Effective date: 20070726 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20070126 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20061127 Ref country code: TR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CY Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20061025 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R082 Ref document number: 60215707 Country of ref document: DE Representative=s name: MERH-IP MATIAS ERNY REICHL HOFFMANN, DE |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R081 Ref document number: 60215707 Country of ref document: DE Owner name: BLACKBERRY LIMITED, WATERLOO, CA Free format text: FORMER OWNER: RESEARCH IN MOTION LTD., WATERLOO, ONTARIO, CA Effective date: 20140925 Ref country code: DE Ref legal event code: R082 Ref document number: 60215707 Country of ref document: DE Representative=s name: MERH-IP MATIAS ERNY REICHL HOFFMANN, DE Effective date: 20140925 Ref country code: DE Ref legal event code: R082 Ref document number: 60215707 Country of ref document: DE Representative=s name: MERH-IP MATIAS ERNY REICHL HOFFMANN PATENTANWA, DE Effective date: 20140925 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 14 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 15 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: PLFP Year of fee payment: 16 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 20211124 Year of fee payment: 20 Ref country code: GB Payment date: 20211129 Year of fee payment: 20 Ref country code: DE Payment date: 20211126 Year of fee payment: 20 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R071 Ref document number: 60215707 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: PE20 Expiry date: 20221126 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Free format text: LAPSE BECAUSE OF EXPIRATION OF PROTECTION Effective date: 20221126 |