Abstract:
In a phase discriminator device for receiving, as input, a complex signal whose argument represents a phase error, and for producing, as output, an estimate of the phase error for each signal sample Z n+1  received, the device includes a frequency discriminator and a computation part for determining the phase estimate obtained at an instant (N+1)T 1 . The computation part determines the phase estimate according to predetermined relations, and T 1  is the time interval between two samples Z n  and Z n+1 , received consecutively.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application claims priority to foreign French patent application No. FR 10 02152, filed on May 21, 2010, the disclosure of which is incorporated by reference in its entirety. 
       FIELD OF THE INVENTION 
       [0002]    The technical field of the invention is that of devices for demodulating a signal in a reception system for which an accurate estimation of the phase is required. The invention applies, for example, to satellite radionavigation systems for which the information contained in the phase of the signal helps in accurately determining the position of a receiver. 
         [0003]    The invention relates to a device for discriminating or estimating the phase or phase error during the demodulation of a signal subject to significant phase variations, for example due to a rapid trend in the frequency due to the Doppler effect. In the subsequent description the use of the term Doppler frequency is made with reference to a frequency bias due to the velocity of the sender or receiver, with respect to the nominal frequency in the absence of motion of the carrier. 
       BACKGROUND OF THE INVENTION 
       [0004]    The aim of a satellite radionavigation system is to determine the position of a receiver by using the principle of triangulation. The pseudo-distances measured on the basis of the signals received by three satellites are utilized to determine the position of the receiver. 
         [0005]    In addition to the accuracy of positioning, two other parameters characterize the performance of a radio-navigation receiver: the capacity to operate with a low signal-to-noise ratio and the resilience to phase dynamics due to the motions of the carrier. In the case where the carrier is an aircraft, the accelerations of motion that it is liable to cause result in a significant Doppler frequency variation in the signals received by the receiver. 
         [0006]    The means for demodulating the signals received by the receiver generally use the known phase-locked loop principle to accurately estimate the phase of the signals received so as to determine an accurate positioning. The performance of a phase-locked loop (PLL) is characterized notably by two criteria, its resilience to rapid variations in Doppler frequency and its capacity to operate at a low signal-to-noise ratio. A compromise is required between these two parameters, thus, the increase in resilience to Doppler frequency variations is achieved at the price of lower resistance to noise. Moreover, it requires an increase in the passband of the loop, thus giving rise to an overhead in terms of computational load for the receiver. 
         [0007]    The resilience to dynamics of the known phase-locked loops is restricted by the phase discriminator which exhibits a limited linear span. The resilience to dynamics of a phase-locked loop represents its capacity to tolerate a significant change between two stationary states. In the case of slackness of the said loop, also called the carrier loop, during the transient phase between the two stationary states, the phase error may depart from the discriminator capture zone and this may give rise to the divergence of the loop and its dropout. 
         [0008]    In order to increase the extent of the phase loop capture zone, it is known to implement a so-called “winder” discriminator which records the number of revolutions over the phase error, the consequence of which is that the capture zone becomes infinite in terms of phase error. It remains, however, limited in terms of frequency error or phase velocity error. This limitation is directly related to the sampling frequency used by the digital device which executes the loop. In a satellite radionavigation receiver, the limitation of the phase loop sampling frequency is related to the duration of coherent integration applied to demodulate the received signal. 
         [0009]    This limit in terms of frequency may be penalizing in the case of applications, for example airborne, inducing large variations in phase velocity. 
         [0010]    Another known solution is described in the applicant&#39;s French patent application No. 07 01931 which pertains to a device for receiving satellite signals comprising a phase loop with delay compensation. This solution makes it possible to increase the passband of the loop without rendering it unstable, the effect of which is to increase the resilience to swing of the loop without increasing the system pre-detection band, therefore without overburdening the computational load of the loops. 
         [0011]    However, this solution, although making it possible to reduce the slackness of the phase loop and thus to improve the resilience to dynamics does not change the operating limit in terms of frequency of the discriminator which still constitutes the limiting factor. 
         [0012]    Another known solution consists in making direct use of a frequency-locked loop but in this case, the phase measurement is much less accurate, this being penalizing to performance in terms of accuracy and integrity of the positioning measurement. 
       SUMMARY OF THE INVENTION 
       [0013]    The invention substantially increases the resilience to phase dynamics without increasing the operating point in terms of signal-to-noise ratio of the system and while limiting the computational load. 
         [0014]    The present invention alleviates the limitations of the aforementioned known solutions by using a frequency-extended winder discriminator. 
         [0015]    The invention includes a phase discriminator device receiving as input a complex signal whose argument represents a phase error Δφ and producing as output an estimate θ n+1  of the phase error Δφ for each signal sample Z n+1  received, comprising computation means suitable for determining the phase estimate obtained at an instant (N+1)T 1  according to the following relations: 
         [0000]      θ n+1 =θ n +Δθ n+1 +Δθ′ n+1  
 
         [0000]      Δθ′ n+1   =Discri   Frequency ( n+ 1)·2 π·T   1  
 
         [0000]      Δθ n+1 =arctan( I′   n+1   ,Q′   n+1 )
 
         [0000]    
       
      
       I′ 
       n+1 
       +jQ′ 
       n+1 
       =Z 
       n+1 
       ·Z′ 
       * 
       n+1  
      
     
         [0000]        Z′   n+1 =exp( jΔθ′   n+1 )· Z   n  
 
         [0000]    where T 1  is the time interval between two samples Z n  and Z n+1  received consecutively, Discri frequency (n+1) is an estimate of the frequency error ∂Δφ/∂t of the complex signal produced by a frequency discriminator that the device comprises, and Z′* n+1  is the complex conjugate of Z′ n+1 . 
         [0016]    In embodiments of the invention, the discriminator device furthermore comprises first means of coherent integration of the signal producing from M complex signal samples Z n,k  at the frequency F′ e =1/T 0  a complex sample Z n  at the frequency F e =1/T 1 , each time interval of duration T 1  being subdivided into M time intervals of duration T 0 , the frequency discriminator determining the frequency error estimate Discri frequency (n) on the basis of the samples Z n,k . 
         [0017]    In another embodiment of the invention, the frequency discriminator comprises at least first means for computing at least two frequency error hypotheses at frequencies −F 0  and +F 0  with F 0 =F e /2: 
         [0000]    
       
         
           
             
               
                 Z 
                 
                   + 
                   1 
                 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   k 
                   = 
                   0 
                 
                 
                   M 
                   - 
                   1 
                 
               
                
               
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
                  
                 
                    
                   
                     
                       + 
                       j 
                     
                      
                     
                         
                     
                      
                     π 
                      
                     
                       k 
                       M 
                     
                   
                 
               
             
           
         
       
       
         
           
             
               
                 
                   Z 
                   
                     - 
                     1 
                   
                 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   ∑ 
                   
                     k 
                     = 
                     0 
                   
                   
                     M 
                     - 
                     1 
                   
                 
                  
                 
                   
                     Z 
                      
                     
                       ( 
                       
                         n 
                         , 
                         k 
                       
                       ) 
                     
                   
                    
                   
                      
                     
                       
                         - 
                         j 
                       
                        
                       
                           
                       
                        
                       π 
                        
                       
                         k 
                         M 
                       
                     
                   
                 
               
             
             , 
           
         
       
     
         [0000]    and second means for computing the frequency error estimate according to the following two relations: 
         [0000]    
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   { 
                   
                     
                       
                          
                         
                           
                             Z 
                             
                               + 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                     - 
                     
                       
                          
                         
                           
                             Z 
                             
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                   
                   } 
                 
                 
                   ( 
                   
                     32 
                      
                     
                         
                     
                      
                     
                       
                         CT 
                         1 
                       
                       / 
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
               
                
               
                   
               
                
               or 
             
           
         
       
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 Re 
                  
                 
                   [ 
                   
                     
                       ( 
                       
                         
                           
                             Z 
                             
                               + 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                         - 
                         
                           
                             Z 
                             
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                       
                       ) 
                     
                     · 
                     
                       
                         
                           Z 
                           0 
                         
                          
                         
                           ( 
                           n 
                           ) 
                         
                       
                       * 
                     
                   
                   ] 
                 
               
               
                 ( 
                 
                   16 
                    
                   
                       
                   
                    
                   
                     
                       CT 
                       1 
                     
                     / 
                     
                       π 
                       2 
                     
                   
                 
                 ) 
               
             
           
         
       
     
         [0000]    with C the estimated power of the received signal and Z 0 (n)=Z n    
         [0018]    In an embodiment of the invention, the frequency error estimate is replaced with the average of the frequency error estimates obtained at the instants n and n+1. 
         [0019]    In another embodiment of the invention, the frequency error hypotheses are corrected with the frequency error estimate obtained at the previous instant: 
         [0000]    
       
         
           
             
               
                 Z 
                 
                   + 
                   1 
                 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   k 
                   = 
                   0 
                 
                 
                   M 
                   - 
                   1 
                 
               
                
               
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
                  
                 
                    
                   
                     
                       + 
                       j 
                     
                      
                     
                         
                     
                      
                     π 
                      
                     
                       k 
                       M 
                     
                   
                 
                  
                 
                    
                   
                     
                       - 
                       j 
                     
                      
                     
                         
                     
                      
                     2 
                      
                     
                         
                     
                      
                     
                       π 
                       · 
                       
                         
                           Discri 
                           Frequency 
                         
                          
                         
                           ( 
                           
                             n 
                             - 
                             1 
                           
                           ) 
                         
                       
                     
                   
                 
               
             
           
         
       
       
         
           
             
               
                 Z 
                 
                   - 
                   1 
                 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   k 
                   = 
                   0 
                 
                 
                   M 
                   - 
                   1 
                 
               
                
               
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
                  
                 
                    
                   
                     
                       - 
                       j 
                     
                      
                     
                         
                     
                      
                     π 
                      
                     
                       k 
                       M 
                     
                   
                 
                  
                 
                    
                   
                     
                       - 
                       j 
                     
                      
                     
                         
                     
                      
                     2 
                      
                     
                         
                     
                      
                     
                       π 
                       · 
                       
                         
                           Discri 
                           Frequency 
                         
                          
                         
                           ( 
                           
                             n 
                             - 
                             1 
                           
                           ) 
                         
                       
                     
                   
                 
               
             
           
         
       
     
         [0000]    and that the frequency error estimate determined at the instant n−1 is added to that obtained at the instant n. 
         [0020]    In a variant embodiment of the invention, the said first computation means carry out a number 2S of frequency error hypotheses at the frequency values ranging from −SF 0  to +SF 0  in steps of F 0  according to the relations 
         [0000]    
       
         
           
             
               
                 Z 
                 
                   + 
                   p 
                 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   k 
                   = 
                   0 
                 
                 
                   M 
                   - 
                   1 
                 
               
                
               
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
                  
                 
                    
                   
                     
                       + 
                       j 
                     
                      
                     
                         
                     
                      
                     π 
                      
                     
                         
                     
                      
                     p 
                      
                     
                       k 
                       M 
                     
                   
                 
                  
                 
                     
                 
                  
                 and 
               
             
           
         
       
       
         
           
             
               
                 
                   Z 
                   
                     - 
                     p 
                   
                 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   ∑ 
                   
                     k 
                     = 
                     0 
                   
                   
                     M 
                     - 
                     1 
                   
                 
                  
                 
                   
                     Z 
                      
                     
                       ( 
                       
                         n 
                         , 
                         k 
                       
                       ) 
                     
                   
                    
                   
                      
                     
                       
                         - 
                         j 
                       
                        
                       
                           
                       
                        
                       π 
                        
                       
                           
                       
                        
                       p 
                        
                       
                         k 
                         M 
                       
                     
                   
                 
               
             
             , 
           
         
       
     
         [0000]    the said frequency discriminator furthermore comprises means for determining the index p 0  associated with the maximum value ∥Z p     0   (n)∥ 2  in the set [−(S−1) . . . , − 2 , − 1 ,  0 ,  1 ,  2 , . . . (S−1)], with Z 0 (n)=Z n  and the frequency error estimate is determined according to one of the following two relations: 
         [0000]    
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   { 
                   
                     
                       
                          
                         
                           
                             Z 
                             
                               
                                 p 
                                 0 
                               
                               + 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                     - 
                     
                       
                          
                         
                           
                             Z 
                             
                               
                                 p 
                                 0 
                               
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                   
                   } 
                 
                 
                   ( 
                   
                     32 
                      
                     
                         
                     
                      
                     
                       
                         CT 
                         1 
                       
                       / 
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
               
               + 
               
                 
                   p 
                   0 
                 
                  
                 
                   F 
                   0 
                 
                  
                 
                     
                 
                  
                 or 
               
             
           
         
       
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   Re 
                    
                   
                     [ 
                     
                       
                         ( 
                         
                           
                             
                               Z 
                               
                                 
                                   p 
                                   0 
                                 
                                 + 
                                 1 
                               
                             
                              
                             
                               ( 
                               n 
                               ) 
                             
                           
                           - 
                           
                             
                               Z 
                               
                                 
                                   p 
                                   0 
                                 
                                 - 
                                 1 
                               
                             
                              
                             
                               ( 
                               n 
                               ) 
                             
                           
                         
                         ) 
                       
                       · 
                       
                         
                           
                             Z 
                             
                               p 
                               0 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                         * 
                       
                     
                     ] 
                   
                 
                 
                   ( 
                   
                     16 
                      
                     
                         
                     
                      
                     
                       
                         CT 
                         1 
                       
                       / 
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
               
               + 
               
                 
                   p 
                   0 
                 
                  
                 
                   F 
                   0 
                 
               
             
           
         
       
     
         [0021]    In another embodiment of the invention, the phase discriminator device furthermore comprises means for estimating the signal-to-noise ratio RŜB, means for comparing the signal-to-noise ratio with a threshold predetermined as a function of a compromise between the extent of the frequency range of the device and its resistance to noise and means for deactivating and activating the frequency discriminator as a function of the result of the comparison. 
         [0022]    In another embodiment of the invention, the phase discriminator device further comprises means for estimating the noise power σ b   2 , for example according to the relation 
         [0000]    
       
         
           
             
               
                 σ 
                 b 
                 2 
               
               = 
               
                 
                   ∑ 
                   
                     k 
                     = 
                     0 
                   
                   
                     M 
                     - 
                     1 
                   
                 
                  
                 
                   
                      
                     
                       
                         Z 
                         b 
                       
                        
                       
                         ( 
                         
                           n 
                           , 
                           k 
                         
                         ) 
                       
                     
                      
                   
                   2 
                 
               
             
             , 
           
         
       
     
         [0000]    with Z b (n,k) being measured noise samples and the signal-to-noise ratio being estimated according to one of the following relations: 
         [0000]    
       
         
           
             
               R 
                
               
                 S 
                 ⋀ 
               
                
               B 
             
             = 
             
               
                 
                   
                     
                        
                       
                         
                           Z 
                           
                             p 
                             0 
                           
                         
                          
                         
                           ( 
                           n 
                           ) 
                         
                       
                        
                     
                     2 
                   
                   
                     σ 
                     b 
                     2 
                   
                 
                  
                 
                     
                 
                  
                 or 
                  
                 
                     
                 
                  
                 R 
                  
                 
                   S 
                   ⋀ 
                 
                  
                 B 
               
               = 
               
                 
                   
                     
                        
                       
                         
                           Z 
                           0 
                         
                          
                         
                           ( 
                           n 
                           ) 
                         
                       
                        
                     
                     2 
                   
                   
                     σ 
                     b 
                     2 
                   
                 
                 . 
               
             
           
         
       
     
         [0023]    The invention also includes a device for receiving satellite radionavigation signals receiving a continuous signal modulated by a carrier, of phase φ received , including means making it possible to generate a local carrier, of phase φ local , and means for demodulating and filtering the signal, furthermore including second means of coherent integration of the said signal which deliver as output M samples Z n,k  at a sampling frequency F′ e  and a phase-locked loop comprising a one phase discriminator device according to the invention which receives as input a complex signal whose argument represents the phase error Δφ=φ local −φ received  and which makes it possible to lock the phase φ local  of the local carrier onto that φ received  of the received carrier. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0024]    The above and other aspects will become more readily apparent on reading the following detailed description given by way of nonlimiting example and taken in view of the accompanying drawings, in which: 
           [0025]      FIG. 1 , a basic diagram of the demodulation means implemented by a satellite radionavigation receiver; 
           [0026]      FIG. 2 , an illustration of the operating principle of a phase discriminator of arctangent type according to the prior art applied to a pilot channel; 
           [0027]      FIG. 3 , an illustration of the operating principle of a phase discriminator of arctangent type according to the prior art applied to a data channel; 
           [0028]      FIG. 4 , an illustration of the operating principle of a “winder” phase discriminator according to the prior art applied to a pilot channel, in the case where the phase loop does not drop out; 
           [0029]      FIG. 5 , an illustration of the operating principle of a “winder” phase discriminator according to the prior art applied to a pilot channel, in the case where the phase loop drops out; 
           [0030]      FIG. 6 , an illustration of the operating principle of a “winder” phase discriminator according to the prior art applied to a data channel; 
           [0031]      FIG. 7 , an illustration of the operating principle of a phase discriminator according to the invention applied to a pilot channel; 
           [0032]      FIG. 8 , a comparative block diagram of a “winder” phase discriminator and of a phase discriminator according to the invention; 
           [0033]      FIG. 9 , a basic diagram of the trend of the phase error over coherent integration sub-intervals; 
           [0034]      FIG. 10 , an illustration of the principle used to construct the frequency discriminator according to the invention; 
           [0035]      FIG. 11 , an illustration of a variant of a frequency discriminator according to the invention; 
           [0036]      FIG. 12 , an illustration of a third variant of a frequency discriminator according to the invention; and 
           [0037]      FIG. 13 , an illustration of the operating principle of a phase discriminator according to the invention applied to a data channel. 
       
    
    
     DETAILED DESCRIPTION 
       [0038]      FIG. 1  shows diagrammatically the principle of the demodulation of the signal received by a satellite radionavigation receiver  10 . The processing actions situated to the left of the dashed line constitute a continuous or discrete system with a high sampling frequency F c , for example equal to 100 MHz. They are carried out preferably by integrated circuits of ASIC or FPGA type which are compatible with a high frequency such as this. The processing actions situated to the right of the dashed line constitute a system sampled at a low frequency, for example equal to 50 Hz, and are carried out by software, for example DSP signal processing software. 
         [0039]    The received signal  101  is demodulated by a local carrier  109  of phase φ local (t) by way of a complex multiplier  102 . The term “local carrier” designates a signal generated locally at a second frequency different from that of the received signal, this second frequency being controlled by virtue of a digitally-controlled oscillator  108 . A second demodulation  103  is also applied with the aim of demodulating the received signal with the local spreading code. Means  105  of coherent integration allow notably the filtering of the demodulated received signal  104 . The complex signal Z(n)=I(n)+jQ(n) thus obtained is transmitted to a phase discriminator  106  which produces as output an estimation θ n  of the phase error Δφ n =φ local (n·T)−φ received (n·T) of the local carrier with respect to the phase of the carrier of the received signal. This estimation feeds a loop corrector  110 . Such a corrector  110  can consist for example of an acceleration integrator which executes the following relation: A(n+1)=A(n)+T·K3·θ(n), associated with a velocity integrator which executes the following relation: V(n+1)=V(n)+T·A(n)+T·K2·θ(n). The digitally-controlled oscillator  108  can also constitute a phase integrator executing the following relation: φ local (n+1)·T)=φ local (n·T)+T·V(n)+T·K1·θ(n) if it receives as input a velocity control equal to: V(n)+K1·θ(n). In this example, the phase loop thus produced is a 3 rd -order loop since it contains three integrators. The acceleration integrator makes it possible to support a constant acceleration with no slackness error in the steady state. On the other hand in the transient state, subsequent to a change of acceleration, the phase slackness error may be significant and the discriminator  106  must be capable of following the maximum error without saturating so that the loop can return to a stationary state with no cycle skipping which would be prejudicial for the accuracy of the measurements 
         [0040]    The loop corrector  110  delivers a velocity control  107  which is produced as input to a digitally-controlled oscillator  108  which delivers at its output the phase of the local carrier φ local (t). Thus the phase-locked loop implemented makes it possible to clamp the phase φ local (t) of the local carrier onto the phase φ received (t) of m the carrier of the received signal. 
         [0041]    The known phase discriminators are of “winder” type. They exhibit a limitation because of the 2π ambiguity on the phase of the signal. When the phase error has altered by π radians, i.e. half a revolution of a circle, between two measurement instants, an ambiguity in the sign of this error persists, namely π or −π. By construction the discriminator  106  knows how to measure a phase error variation included in the interval [+π−π[. In the presence of a large variation in the phase, the discriminator measures a phase error equal to (+π+ε) as being equal to (−π+ε) thereby leading to an error of 2π, liable to cause the loop to drop out. 
         [0042]    As long as the variation in the phase error is strictly limited to π in absolute value between two measurement instants the discriminator operates linearly. Beyond this the loop is liable to drop out. The limit of such a “winder” discriminator device in terms of phase velocity error (or phase error rate) is π/T in radians per second or λ/T in metres per second, λ being the wavelength and T being the duration between two correlation samples, or else πF e , with F e  the sampling frequency of the device which executes the phase discriminator  106 . 
         [0043]    With such a device, resilience to a high phase velocity error can only be achieved by increasing the sampling frequency F e  thereby giving rise to an increase in the computational load. 
         [0044]    The received signal  101  can be obtained from a data channel or from a pilot channel. In the first case, the signal  101  is composed of a carrier modulated by a spreading code and by data. In the second case, the signal  101  is composed solely of a carrier and of a spreading code. 
         [0045]    The spreading code is known to the receiver, contrary to the data which contain information about the position of the satellites. 
         [0046]      FIG. 2  illustrates the operating principle of a phase discriminator of arctangent type according to the prior art applied to a pilot channel. 
         [0047]    At the top of the figure are represented two complex symbols Z n  and Z n+1  in the complex plane (I,Q) and their respective phases θ n  and θ n+1 . The two symbols correspond to two successive samples of the complex signal obtained as output from the means  105  of coherent integration which allow the filtering of the received signal after demodulation by the local carrier and the local code. The estimate of the phase error θ n  is obtained by computing the arctangent between the real part and the imaginary part of the symbol Z n . At the bottom of  FIG. 2  is represented the trend of the phase estimate θ n  as a function of time t=nT=n/F e . It can be seen that the output of the discriminator is limited between the values π and −π in radians, thereby reducing the robustness of the loop to phase velocity variations. When the slackness error exceeds π radians, the loop may diverge and therefore drop out because of the feedback action at the output of the corrector which may change sign. 
         [0048]      FIG. 3  illustrates the operating principle of a phase discriminator of arctangent type according to the prior art applied to a data channel. 
         [0049]    In this second example, the symbols Z n , Z n+1  obtained after coherent integration  105  over the interval [nT, (n+1)T] coincide with binary data [−1, +1] which modulate the signal. These data are not known to the receiver, thus a change of sign of the data at the instant (n+1)T will result in a change of sign between Z n  and Z n+1  that must not be confused with a phase variation. The consequence is that the output of the discriminator is now limited to the interval [−π/2, π/2] as illustrated at the bottom of  FIG. 3 . 
         [0050]    The principle of the “winder” phase discriminator according to the prior art consists in using the measured phase error θ n  as reference to resolve the ambiguity in the phase error θ n+1  and remove the limitations of the discriminator of arctangent type. 
         [0051]    We have θ n+1 =θ n +Δθ n+1  with Δθ n+1 =arctangent(I n+1 ′, Q n+1 ′) ε]−π, π] and Z′ n+1 =I n+1 ′+j Q n+1 ′=Z n+1 . Z n * with Z n+1 =(I n+1 +j Q n+1 ), Z n =(I n +j Q n ) and Z n *=(I n −j Q n ) the conjugate value of Z n . 
         [0052]      FIG. 4  illustrates two examples for which the “winder” phase discriminator according to the prior art, applied to a pilot channel, operates without dropout. In these two examples, the phase jump Δθ n+1  between two consecutive symbols is less than π in absolute value. The trend of the phase θ n  as a function of time is continuous and does not cause any dropout. 
         [0053]      FIG. 5  presents an example for which the phase jump between the two consecutive symbols Z n  and Z n+1  is greater than π in absolute value. In this case, and as illustrated in the bottom part of  FIG. 5 , the output of the phase discriminator is no longer limited in itself. On the other hand, a limit persists in the variation of the phase error Δφ n  between two measurements. As long as this variation remains lower in absolute value than π then operation is linear as illustrated in the bottom part of  FIG. 4  and the phase loop does not drop out. On the other hand if this variation exceeds π in absolute value then θ n+1  no longer reflects the real phase error Δφ n+1 =φ local ((n+1)·T)−φ received ((n+1)·T) and the phase loop may drop out. 
         [0054]      FIG. 6  illustrates the operation of the “winder” phase discriminator applied to a data channel for which the phase jump limitation is now equal to π/2 in absolute value. Because of the ambiguity in the sign of Z n+1  due to the modulation by the data, we refer to the right half-plane limited by the perpendicular to the axis I′ which contains the symbol Z n . This amounts to limiting Δθ n+1  in the interval ]−π/2, π/2], thus reducing by two the tolerance limit in terms of Doppler frequency error. 
         [0055]    We have θ n+1 =θ n +Δθ n+1 , Δθ n+1 =arctangent(Q n+1 ′ I n+1 ′) ε[−π/2, π/2] 
         [0000]    with I n+1 ′+j Q n+1 ′=Z n+1 . Z n *, Z n+1 =(I n+1 +j Q n+1 ), Z n =(I n +j Q n ) and Z n *=(I n −j Q n ) the conjugate value of Z n . 
         [0056]    The symbol Z n+1  is estimated by such a discriminator as being the symbol Z′ n+1 =−Z n+1  for which the phase error Δθ n+1  which separates it from the previous symbol Z n  is lower in absolute value than π/2. 
         [0057]      FIG. 7  illustrates the principle of a phase discriminator according to the invention. 
         [0058]    In order to remove the limitation on the variation in phase error inherent to known phase discriminators, the present invention is based on the addition of a frequency discriminator making it possible to measure, upstream of the phase discriminator, the velocity of the phase error or Doppler frequency error. 
         [0059]    Aspects of the invention may be explained with reference to the following relations and in accordance with the example of  FIG. 7 . Two consecutive symbols Z n  and Z n+1  are represented in the complex plane (I,Q). The phase jump between these two symbols is greater than π thus corresponding to a typical case for which the known phase discriminators do not enable correct convergence of the phase loop to be ensured. 
         [0060]    The invention includes using an intermediate symbol Z′ n+1 =exp j(Δθ n+1 ′). Z n  constructed by virtue of a frequency discriminator which makes it possible to measure the Doppler frequency error ∂Δφ/∂t or local carrier frequency error with respect to the frequency of the received carrier. The phase error Δθ n+1 ′ is obtained on the basis of the following relation: 
         [0000]      Δθ n+1   ′=Discri   Frequency ( n+ 1)·2 π·T=Discri   Frequency ( n+ 1)·2 π/F   e  
 
         [0000]    where Discri Frequency (n+1) is an estimation of the mean Doppler frequency error (or variation of the phase error) between the instants nT and (n+1)T. 
         [0061]    The phase error between the two consecutive symbols can thus be decomposed into a first part Δθ n+1 ′ that can be predicted by virtue of the estimation of the Doppler frequency error and a second part Δθ n+1  which corresponds to the non-predictable phase error. This results in the following relations: 
         [0000]      θ n+1 =θ n +Δθ n+1 ′+Δθ n+1   (1)
 
         [0000]    with Δθ n+1 =arctangent(I′ n+1 , Q′ n+1 )ε]−π, π] and 
         [0000]        I   n+1   ′+jQ   n+1   ′=Z   n+1   ·Z   n+1 ′*  (2)
 
         [0062]    The use of a frequency discriminator makes it possible to pre-correct the signal so as to refer the residual phase error Δθ n+1  to an interval ]−π, π] compatible with the operation of a “winder” phase discriminator. 
         [0063]      FIG. 8  shows diagrammatically for comparative purposes a block diagram of a “winder” phase discriminator  91  and an overview of a phase discriminator  92  according to the invention. The latter uses a frequency discriminator  921  applied to the complex signal Z(n,k) received before coherent integration at a higher sampling frequency. Applying the frequency discriminator to the input signal of the means  105  of coherent integration makes it possible to estimate a higher Doppler frequency than if it were applied to the output signal of the said means  105  which exhibit a low sampling frequency due to the sub-sampling performed by the means of coherent integration  105 . The means  105  can also be separated into two successive stages as will be explained further on in the description. A first stage performs a coherent integration over a duration T 0  for example a duration of 1 ms, the result of which is used by the frequency discriminator  921  and then a second stage finalizes the coherent integration over a duration T 1 , for example a duration of 20 ms which corresponds to the sampling period of the signal processing software which executes the phase discriminator according to the invention. 
         [0064]    The frequency discriminator  921  produces as output an estimate of predictable phase error variation Δ′θ n  which is utilized in accordance with relations (1) and (2) to produce the signals I′ n  and Q′ n  on the basis of which the variation in non-predictable phase error Δθ n  is computed. The latter is added to the predictable phase error variation Δ′θ n  delivered as output from the frequency discriminator  921  so as to produce an accurate global phase error variation estimate which is thereafter produced at the input of an integrator  922 . 
         [0065]    The principle used for the construction of the frequency discriminator according to the invention is now explained. 
         [0066]    The means  105  of coherent integration carry out notably a coherent integration of the demodulated received signal over a time interval of duration T 1 . This interval is cut up into a number M of elementary integration sub-intervals of duration T 0 . The integration means  105  provide the phase discriminator according to the invention with the results of the coherent integrations over a duration T 1 , which is performed by the software part of the demodulation system described in  FIG. 1 . In parallel, two filters, left and right, which carry out two symmetric Doppler frequency error hypotheses, are applied to the results of the coherent integrations over a duration T 0  so as to construct a frequency discriminator making it possible to extrapolate the phase error variation. 
         [0067]      FIG. 9  shows diagrammatically the trend of the phase of the received signal φ received (t) on the one hand and of the local phase φ local (t) on the other hand. F e =1/T 1  is the sampling frequency, or processing frequency, of the phase discriminator device according to the invention. The symbol Z(n) of phase θ(n) is spread over the whole of the interval of duration T 1 . Its phase θ(n) corresponds to the mean phase over this time interval. The interval [(n−1)T 1  nT 1 ] is cut up into a number M of sub-intervals I(n,k) of duration T 0  for which a symbol Z(n,k) is produced by the integration means  105  which thus carry out a piecewise integration. 
         [0068]    A first Doppler error hypothesis corresponding to the frequency F 0 =F e /2=1/(2T 1 ) is carried out on the basis of the symbols Z(n,k): 
         [0000]    
       
         
           
             
               
                 Z 
                 
                   + 
                   1 
                 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   k 
                   = 
                   0 
                 
                 
                   M 
                   - 
                   1 
                 
               
                
               
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
                  
                 
                    
                   
                     
                       + 
                       j 
                     
                      
                     
                         
                     
                      
                     π 
                      
                     
                       k 
                       M 
                     
                   
                 
               
             
           
         
       
     
         [0069]    A second Doppler error hypothesis corresponding to the frequency −F 0  is also carried out: 
         [0000]    
       
         
           
             
               
                 Z 
                 
                   - 
                   1 
                 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   k 
                   = 
                   0 
                 
                 
                   M 
                   - 
                   1 
                 
               
                
               
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
                  
                 
                    
                   
                     
                       - 
                       j 
                     
                      
                     
                         
                     
                      
                     π 
                      
                     
                       k 
                       M 
                     
                   
                 
               
             
           
         
       
     
         [0070]    Computation of the quantity ∥Z +1 (n)∥ 2 −∥Z −1 (n)∥ 2  makes it possible to obtain a frequency error estimate proportional to the Doppler frequency error if the latter remains limited to the linear zone between −F 0  and +F 0  such as is illustrated in  FIG. 10 . 
         [0071]    The frequency discriminator according to the invention produces as output the following estimate expressed in Hertz: 
         [0000]    
       
         
           
             
               
                 
                   Discri 
                   Frequency 
                 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   { 
                   
                     
                       
                          
                         
                           
                             Z 
                             
                               + 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                     - 
                     
                       
                          
                         
                           
                             Z 
                             
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                   
                   } 
                 
                 
                   ( 
                   
                     32 
                      
                     
                         
                     
                      
                     
                       
                         CT 
                         1 
                       
                       / 
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
               
             
              
             
                 
             
           
         
       
     
         [0000]    with C the power of the signal received, as estimated for example by way of the following relation C=∥Z(n)∥ 2 . 
         [0072]    In another embodiment of the invention, a third Doppler error hypothesis, termed centred, is utilized. This hypothesis is determined by the following relation: 
         [0000]    
       
         
           
             
               
                 Z 
                 0 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   ∑ 
                   
                     k 
                     = 
                     0 
                   
                   
                     M 
                     - 
                     1 
                   
                 
                  
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
               
               = 
               
                 Z 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
             
           
         
       
     
         [0073]    In this case the frequency discriminator according to the invention produces as output the following estimate expressed in Hertz: 
         [0000]    
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 Re 
                  
                 
                   [ 
                   
                     
                       ( 
                       
                         
                           
                             Z 
                             
                               + 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                         - 
                         
                           
                             Z 
                             
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                       
                       ) 
                     
                     · 
                     
                       
                         
                           Z 
                           0 
                         
                          
                         
                           ( 
                           n 
                           ) 
                         
                       
                       * 
                     
                   
                   ] 
                 
               
               
                 ( 
                 
                   16 
                    
                   
                       
                   
                    
                   
                     
                       CT 
                       1 
                     
                     / 
                     
                       π 
                       2 
                     
                   
                 
                 ) 
               
             
           
         
       
     
         [0000]    with Z 0 (n)*, the complex conjugate of Z 0 (n). 
         [0074]      FIG. 11  illustrates another embodiment of the frequency discriminator according to the invention. At the instant (n+1)T 1 , the estimates Discri frequency (n) and Discri frequency (n+1) of the phase variations at the instants n and n+1 are available. Discri frequency (n) is representative of the variation of the phase error over the time interval [(n−1)T 1  nT 1 ] and Discri frequency (n+1) is representative of the variation of the phase error over the time interval [nT 1  (n+1)T 1 ]. These two estimates are represented in  FIG. 11  respectively by the slopes of the straight lines followed by the phase error δφ as a function of time over each interval of duration T 1 . In order to improve the accuracy of the frequency discriminator, the average of two successive estimations is used so as to produce as output the following result: 
         [0000]    
       
         
           
             
               
                 
                   Discri 
                   Frequency 
                 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               + 
               
                 
                   Discri 
                   Frequency 
                 
                  
                 
                   ( 
                   
                     n 
                     + 
                     1 
                   
                   ) 
                 
               
             
             2 
           
         
       
     
         [0075]    In another embodiment of the invention, the frequency discriminator output computed at the instant (n−1)T 1  is used to correct the left and right Doppler error hypotheses performed initially at the frequencies F 0  and −F 0 . The following relations are used: 
         [0000]    
       
         
           
             
                 
             
              
             
               
                 
                   Z 
                   
                     + 
                     1 
                   
                 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   ∑ 
                   
                     k 
                     = 
                     0 
                   
                   
                     M 
                     - 
                     1 
                   
                 
                  
                 
                   
                     Z 
                      
                     
                       ( 
                       
                         n 
                         , 
                         k 
                       
                       ) 
                     
                   
                    
                   
                      
                     
                       
                         + 
                         j 
                       
                        
                       
                           
                       
                        
                       π 
                        
                       
                         k 
                         M 
                       
                     
                   
                    
                   
                      
                     
                       
                         - 
                         j 
                       
                        
                       
                           
                       
                        
                       2 
                        
                       
                           
                       
                        
                       
                         π 
                         · 
                         
                           
                             Discri 
                             Frequency 
                           
                            
                           
                             ( 
                             
                               n 
                               - 
                               1 
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
             
           
         
       
       
         
           
             
                 
             
              
             
               
                 
                   Z 
                   
                     - 
                     1 
                   
                 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   ∑ 
                   
                     k 
                     = 
                     0 
                   
                   
                     M 
                     - 
                     1 
                   
                 
                  
                 
                   
                     Z 
                      
                     
                       ( 
                       
                         n 
                         , 
                         k 
                       
                       ) 
                     
                   
                    
                   
                      
                     
                       
                         - 
                         j 
                       
                        
                       
                           
                       
                        
                       π 
                        
                       
                         k 
                         M 
                       
                     
                   
                    
                   
                      
                     
                       
                         - 
                         j 
                       
                        
                       
                           
                       
                        
                       2 
                        
                       
                           
                       
                        
                       
                         π 
                         · 
                         
                           
                             Discri 
                             Frequency 
                           
                            
                           
                             ( 
                             
                               n 
                               - 
                               1 
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
             
           
         
       
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   Re 
                    
                   
                     [ 
                     
                       
                         ( 
                         
                           
                             
                               Z 
                               
                                 + 
                                 1 
                               
                             
                              
                             
                               ( 
                               n 
                               ) 
                             
                           
                           - 
                           
                             
                               Z 
                               
                                 - 
                                 1 
                               
                             
                              
                             
                               ( 
                               n 
                               ) 
                             
                           
                         
                         ) 
                       
                       · 
                       
                         
                           
                             Z 
                             0 
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                         * 
                       
                     
                     ] 
                   
                 
                 
                   ( 
                   
                     16 
                      
                     
                       
                         CT 
                         1 
                       
                       / 
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
               
               + 
               
                 
                   Discri 
                   Frequency 
                 
                  
                 
                   ( 
                   
                     n 
                     - 
                     1 
                   
                   ) 
                 
               
             
           
         
       
     
         [0076]    In another embodiment of the invention, such as illustrated in  FIG. 12 , in order to further widen the operating domain in terms of velocity slackness error of the phase discriminator, a number S of Doppler error hypotheses are performed at the following frequency values: −S F 0 , . . . , −2 F 0 , −F 0 , 0, F 0 , 2F 0 , . . . SF 0 . 
         [0077]    These N hypotheses may be written in a general manner with the aid of the following relation: 
         [0000]    
       
         
           
             
               
                 Z 
                 p 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   k 
                   = 
                   0 
                 
                 
                   M 
                   - 
                   1 
                 
               
                
               
                 
                   Z 
                    
                   
                     ( 
                     
                       n 
                       , 
                       k 
                     
                     ) 
                   
                 
                  
                 
                    
                   
                     
                       - 
                       j 
                     
                      
                     
                         
                     
                      
                     
                       π 
                       · 
                       p 
                     
                      
                     
                       k 
                       m 
                     
                   
                 
               
             
           
         
       
     
         [0000]    for p varying −S to S. 
         [0078]    The index p 0  corresponding to the maximum value ∥Z p     0   (n)∥ 2  among all the computed hypotheses is retained so as to construct the frequency discriminator around this maximum with the aid of one of the following two relations: 
         [0000]    
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   { 
                   
                     
                       
                          
                         
                           
                             Z 
                             
                               
                                 p 
                                 0 
                               
                               + 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                     - 
                     
                       
                          
                         
                           
                             Z 
                             
                               
                                 p 
                                 0 
                               
                               - 
                               1 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                          
                       
                       2 
                     
                   
                   } 
                 
                 
                   ( 
                   
                     32 
                      
                     
                         
                     
                      
                     
                       
                         CT 
                         1 
                       
                       / 
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
               
               + 
               
                 
                   p 
                   0 
                 
                  
                 
                   F 
                   0 
                 
               
             
           
         
       
       
         
           
             
               
                 Discri 
                 Frequency 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   Re 
                    
                   
                     [ 
                     
                       
                         ( 
                         
                           
                             
                               Z 
                               
                                 
                                   p 
                                   0 
                                 
                                 + 
                                 1 
                               
                             
                              
                             
                               ( 
                               n 
                               ) 
                             
                           
                           - 
                           
                             
                               Z 
                               
                                 
                                   p 
                                   0 
                                 
                                 - 
                                 1 
                               
                             
                              
                             
                               ( 
                               n 
                               ) 
                             
                           
                         
                         ) 
                       
                       · 
                       
                         
                           
                             Z 
                             
                               p 
                               0 
                             
                           
                            
                           
                             ( 
                             n 
                             ) 
                           
                         
                         * 
                       
                     
                     ] 
                   
                 
                 
                   ( 
                   
                     16 
                      
                     
                         
                     
                      
                     
                       
                         CT 
                         1 
                       
                       / 
                       
                         π 
                         2 
                       
                     
                   
                   ) 
                 
               
               + 
               
                 
                   p 
                   0 
                 
                  
                 
                   F 
                   0 
                 
               
             
           
         
       
     
         [0079]    In another embodiment of the invention, the use of the frequency discriminator is conditioned on the estimated signal-to-noise ratio. The frequency-extended phase discriminator according to the invention makes it possible to substantially increase the phase-wise operating zone of the phase-locked loop. However the use of the frequency discriminator is liable to render the phase loop less robust to thermal noise or to noise due to interference. To alleviate this drawback, the output of the frequency discriminator is forced to zero when the estimated signal-to-noise ratio falls below a given threshold, for example a threshold of 10 dB. 
         [0080]    The estimation of the signal-to-noise ratio is performed on the basis of the Doppler error hypotheses previously computed on the basis of one of the following relations: 
         [0000]    
       
         
           
             
               R 
                
               
                 S 
                 ⋀ 
               
                
               B 
             
             = 
             
               
                 
                   
                     
                        
                       
                         
                           Z 
                           
                             p 
                             0 
                           
                         
                          
                         
                           ( 
                           n 
                           ) 
                         
                       
                        
                     
                     2 
                   
                   
                     σ 
                     b 
                     2 
                   
                 
                  
                 
                     
                 
                  
                 or 
                  
                 
                     
                 
                  
                 R 
                  
                 
                   S 
                   ⋀ 
                 
                  
                 B 
               
               = 
               
                 
                   
                      
                     
                       
                         Z 
                         0 
                       
                        
                       
                         ( 
                         n 
                         ) 
                       
                     
                      
                   
                   2 
                 
                 
                   σ 
                   b 
                   2 
                 
               
             
           
         
       
     
         [0000]    where σ b   2  is a noise estimation carried out by the receiver over the integration interval of duration T 1 , for example, with the aid of the relation 
         [0000]    
       
         
           
             
               
                 σ 
                 b 
                 2 
               
               = 
               
                 
                   ∑ 
                   
                     k 
                     = 
                     0 
                   
                   
                     M 
                     - 
                     1 
                   
                 
                  
                 
                   
                      
                     
                       
                         Z 
                         b 
                       
                        
                       
                         ( 
                         
                           n 
                           , 
                           k 
                         
                         ) 
                       
                     
                      
                   
                   2 
                 
               
             
             , 
           
         
       
     
         [0000]    where Z b (n,k) are noise samples computed, over the time interval of duration T 1 , on a dedicated demodulation channel containing only noise as output. 
         [0081]      FIG. 13  illustrates the operation of the phase discriminator according to the invention, applied to a data channel. The principle used is identical, except for the difference that the phase error Δθ n+1  exhibits an ambiguity of π/2 instead of π in the case of a pilot channel. The tolerance to dynamics in phase error is therefore lower and the maximum phase error velocity for which the phase discriminator according to the invention will operate correctly is substantially half that in the case of a pilot channel. 
         [0082]    In another embodiment of the invention, the second means of coherent integration  105  are removed and the frequency discriminator device  921  carries out an estimation of the frequency error directly on the basis of the signal received at the frequency F e . The frequency error estimate can also be carried out by a device external to the invention and is in this case provided as input parameter to the phase discriminator according to the invention. 
         [0083]    The present invention applies to all radionavigation receivers, and more broadly to any receiver including means for demodulating a carrier for which there is a desire to accurately estimate the phase or phase error. It makes it possible either to increase the resilience to phase dynamics while limiting the computational load, or to reduce the phase loop sampling frequency for equivalent performance. 
         [0084]    In the case of a receiver suitable for demodulating a signal whose spectrum has not been spread, the means  105  of coherent integration are merged with filtering and sub-sampling means. In the case of a radionavigation receiver, the means  105  of coherent integration implicitly carry out sub-sampling of the signal.