Abstract:
The invention relates to methods of processing signals subjected to interference. It consists, when this signal is formed of a first wideband channel and of a second narrowband channel modulating one and the same carrier in which these signals are cut off for the duration of the interference, in using a single phase loop and a single code loop to process these two channels simultaneously. It makes it easier to process signals in satellite based navigation systems of the GPS type when they are scrambled by the DME-type distance measuring system signals, as well as more generally making it possible to increase the robustness of the processing by joint processing of the bands broadcast to the users.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
     The present Application is based on International Application No. PCT/FR03/00744, filed on Mar. 7, 2003, entitled “METHOD FOR TREATING A SIGNAL IN THE PRESENCE OF INTERFERENCE”, which in turn corresponds to French Application FR 02/03633filed on Mar. 22, 2002, and priority is hereby claimed under 35 USC §119 based of these applications. Each of these applications are hereby incorporated by reference in their entirety into this application. 
     FIELD OF THE INVENTION 
     The present invention relates to methods of signal processing in the presence of interference, more particularly when this signal is composed of a single carrier modulated according to two frequency bands, one wide and the other narrow, and when the interference consists of discrete spectral lines originating for example from pulses distributed according to a quasi-random mode. More generally, the invention makes it possible to increase the robustness of the processing of satellite based navigation signals composed of a carrier modulated by two frequency bands. 
     DESCRIPTION OF THE RELATED ART 
     The invention applies in particular in this framework to the so-called E5b band modulation of the satellite based navigation systems currently under study, and when these systems are scrambled by the signals of so-called DME distance measuring systems. It applies also to other frequency bands, such as L 1 , L 2 , G 1 , G 2 , E 6 , C and possibly L 5 . These bands are those of the GPS, GLONASS and GALILEO systems. To simplify the description, the latter will be limited in the remainder of this text to the so-called E5 signal provided for the Galileo navigation system, when the latter is scrambled by interference originating from the signals of the DME type distance measuring system, which is generally coupled with VOR beacons. 
     The spectrum of this signal is represented in  FIG. 1 . It is formed by a first signal having a main lobe in a narrow band called NB of 2 MHz for example and by a second signal having a main lobe in a wide band called WB of 20 MHz for example, on carriers at the frequency F 0 . In the example considered the power emitted is distributed equally between the wide and narrow bands, so that for identical power on the two signals, the main lobe of the spectrum in the case of the wide band is not as high as in the case of the narrow band, as may be seen in the figure. The frequencies of the spreading codes are synchronized between the wide band and the narrow band and the data streams are the same between these two bands. These signals are therefore synchronous and in phase (possibly in quadrature). The narrowband signal may be of the BPSK (Binary Phase Shift Keying) or BOC (Binary Offset Carrier) type, while the wideband signal may also be of the BPSK or BOC type. 
     Given the powers involved in satellites on the one hand, and in DME systems on the other hand, the DME unwanted signals cause interference spectral lines whose power is appreciably greater than useful signals. 
     To eliminate the scrambling originating from these DME signals, it is known to perform a processing in which the signal received is cut off when such a scrambling signal is detected, this quite obviously entailing a disturbance to the reception of the useful signals. This system is known by the term “blanking”. 
     According to the number of DME signals received, the proportion of cutoffs, or of “blanking”, is greater or lesser and ultimately one no longer receives enough signal and there is then no means of recovering a useful signal. 
     In practice, such a useful signal is nevertheless obtained and it can be processed and the desired information can be derived from it. 
     SUMMARY OF THE INVENTION 
     To improve the results, the invention proposes a method of signal processing in the presence of interference, this signal comprising a first wideband channel and a second narrowband channel modulating one and the same carrier, in which single phase loop and a single code loop to process the two channels simultaneously. 
     According to another characteristic, in the phase loop a signal is generated at the frequency of the carrier, it is sine and cosine phase-shifted and the first channel and the second channel are demodulated synchronously with these two phase-shifted signals so as to obtain a signal I WB  and a signal Q WB  and a signal I NB  and a signal Q NB  respectively. 
     According to another characteristic, in the code loop a signal is generated at the frequency of the code of the wideband signal, a first punctual code and a first Δ code are generated on the basis of this signal, the first Δ code is multiplied respectively with the signal I NB  and the signal Q NB , the first punctual code is multiplied with the signal I NB  and the signal Q NB , the signals resulting from these multiplications are integrated to obtain signals I PNB , I ΔNB , Q PNB  and Q ΔNB , a second Δ signal and a second punctual signal are generated on the basis of the code frequency signal, the second Δ signal is multiplied respectively with the signal I WB  and the signal Q WB , the second punctual signal is multiplied by the signal I WB  and the signal Q WB , and the signals resulting from these multiplications are integrated to obtain signals I PWB , I ΔWB , Q PWB  and Q ΔWB . 
     According to another characteristic, the signals I PWB , I ΔWB , Q PNB  and Q ΔNB  are used in a first phase discriminator common to the two signals, wideband and narrowband, to obtain the phase of the carrier and to drive, after filtering, the oscillator making it possible to obtain said signal at the frequency of the carrier. 
     According to another characteristic, in the first phase discriminator the input signals are mixed before the discrimination proper according to proportions determined by the formulae:
 
 I   p opt   =αI   p WB +β Ip NB  
 
 Q   p opt   =αQ   p WB   +βQ   p NB  
 
     According to another characteristic, the parameters α and β are determined by the formulae: 
             α   =       σ   NB   2         σ   NB   2     +     σ   WB   2                     β   =       σ   WB   2         σ   NB   2     +     σ   WB   2               
in which the parameters σ NB  and σ WB  are estimated on the basis of the number of samples deleted C and of the total number of samples N 0  by the formulae:
 σ wb   2 =1/( C/N   0 WB ) σ Nb   2 =1/( C/N   0 NB ) 
     According to another characteristic, the signals I PWB , I ΔWB , Q PWB , Q ΔWB , I PNB , I ΔNB , Q PNB  and Q ΔNB  are used in a second phase discriminator common to the two channels, wideband and narrowband, to obtain the setting of the two codes and to drive, after filtering, the oscillator making it possible to obtain said signal at the code frequency. 
     According to another characteristic, in the second phase discriminator the input signals are mixed before the discrimination proper according to proportions determined by the formulae:
 
 I   Δopt =( aI   ΔWB   +bI   ΔNB )
 
 Q   Δopt =( aQ   ΔWB   +bQ   ΔNB )
 
     According to another characteristic, the parameters a and b are determined by the formulae: 
             a   =           ρ   WB   2     ⁢     σ   NB   2         (         ρ   NB   2     ⁢     σ   WB   2       +       ρ   WB   2     ⁢     σ   NB   2         )       ×     1     S   WB                     b   =           ρ   NB   2     ⁢     σ   WB   2         (         ρ   NB   2     ⁢     σ   WB   2       +       ρ   WB   2     ⁢     σ   NB   2         )       ×     1     S   NB               
in which the terms S WB  and S MB  correspond to the slopes of the response curves of the discriminators as measured at the central level of these curves, and the terms ρ WB  and ρ MB  are the normalized slopes of the discriminators as obtained by dividing the slope S by a factor σ Δ  determined by the formula:
 σ Δ =[1 −A (2 d )] 
in which A is the normalized authorization function of the filtered spreading code such that A(0)=1, and 2d is the value of the deviation between the early code and the late code used to generate the Δ code.
 
     Still other objects and advantages of the present invention will become readily apparent to those skilled in the art from the following detailed description, wherein the preferred embodiments of the invention are shown and described, simply by way of illustration of the best mode contemplated of carrying out the invention. As will be realized, the invention is capable of other and different embodiments, and its several details are capable of modifications in various obvious respects, all without departing from the invention. Accordingly, the drawings and description thereof are to be regarded as illustrative in nature, and not as restrictive. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
       Other features and advantages of the invention will become clearly apparent in the following description, presented by way of nonlimiting example and with regard to the appended figures which represent: 
         FIG. 1 , a graph of the spectrum of the signal to be processed; and 
         FIG. 2 , a block diagram of a device allowing the implementation of the method according to the invention. 
     
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     The input signals of the device making it possible to implement the invention are filtered in standard fashion, retaining in the case of the wideband signal the whole of the main lobe of the spectrum, so as to preserve the accuracy of the code, and eliminating the sidelobes of the narrowband signal insofar as they exist. This makes it possible to limit the probability of interference by the DME signals, which are statistically distributed over the whole band. As may be seen in  FIG. 1 , a large part of the DME signals are situated outside the main lobe of the narrowband signal and the temporal occupancy of DME interference is statistically ten times smaller for the narrow band than for the wide band, in the example considered. 
     If the “blanking” technique is then used, the losses, expressed as the number of samples C deleted over the total number of samples N0 are given by:
 
 C/N   0 =1−θ  (1)
 
and in decibels by:
 
 C/N   0dB =10 log(1−θ)   (2)
 
     In this formula, θ is the proportion of cutoff time, which is statistically 10 times smaller for the narrowband signal than for the wideband signal, if the filtering described above is performed properly. 
     The output of the phase loop, which will be described later, is applied to a digitally controlled oscillator  201  which delivers a signal at the frequency of the carrier. This signal is then sine and cosine phase-shifted in two phase-shifters  202 ,  203 . 
     The cosine signal is then applied to the wideband signal in a multiplier  204  and to the narrowband signal in a multiplier  205 . 
     The sine signal is applied to the wideband signal in a multiplier  206  and to the narrowband signal in a multiplier  207 . 
     The carrier is thus demodulated on each channel so as to obtain on each of these channels the signals known conventionally as I and Q. 
     The output of the code loop, which will be described later, is applied to a digitally controlled oscillator  208  which delivers a signal at the frequency of the code of the wideband signal. The latter is applied to a first spreading code generator  209  which delivers on the one hand the so-called punctual code and on the other hand the so-called Δ code, corresponding to the difference between a late code and an early code, these codes being adapted to that received on the narrowband channel. 
     The Δ code is multiplied with the signal I of the narrowband channel in a multiplier  210  and with the signal Q of the same channel in a multiplier  211 . 
     The punctual code is multiplied with the signal I of this narrowband channel in a multiplier  212  and with the signal Q of the same channel in a multiplier  213 . 
     The code frequency signal output by the oscillator  208  is also applied to a second code signal generator  214 , similar to the generator  209  but adapted to the code present on the wideband channel. It too delivers a punctual code and a Δ code. 
     The Δ code is multiplied with the signal I of the wideband channel in a multiplier  215  and with the signal Q of the same channel in a multiplier  216 . 
     The punctual signal is multiplied with the channel I of this wideband channel in a multiplier  217  and with the signal Q of this channel in a multiplier  218 , 
     The signals output by all these multipliers are integrated in an assembly of integrators  219 , one for each signal resulting from the multiplications. 
     In the wideband channel, the signals discharged by the integrators are respectively I pWB , I ΔWB , Q pWB , Q ΔWB , this notation signifying that one is dealing with the I or Q signal for the punctual code or the Δ code in the wideband channel. 
     In the same manner, the integrators of the narrowband channel discharge four signals I pNB , I ΔNB , Q pNB , Q ΔNB . 
     According to the invention, these signals are used in two different phase discriminators, which are however common to the two signals, wideband and narrowband, one  220  to obtain the phase of the carrier and the other  222  to obtain the setting of the two codes of the channels. 
     According to the invention, the phase discriminator  220  can use a technique known by the name extended discriminator, more adapted in aeronautical applications so as to be able to support a very wide dynamic range. 
     Conventionally, such a discriminator operates on the basis of a I signal and of a Q signal, and according to the invention these signals are optimized signals called Ip opt and Qp opt obtained by mixing the signals I pWB , Q pWB , I pNB , and Q pNB  using the formulae:
 
 I   p opt   =α I   p WB +β IpNB    (3)
 
 Q   p opt   =αQ   p WB   +βQ   p NB    (4)
 
     In these formulae α and β are given by: 
     
       
         
           
             
               
                 
                   α 
                   = 
                   
                     
                       σ 
                       NB 
                       2 
                     
                     
                       
                         σ 
                         NB 
                         2 
                       
                       + 
                       
                         σ 
                         WB 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
             
               
                 
                   β 
                   = 
                   
                     
                       σ 
                       WB 
                       2 
                     
                     
                       
                         σ 
                         NB 
                         2 
                       
                       + 
                       
                         σ 
                         WB 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
     The parameters σ NB  and σ WB  correspond respectively to the noise levels on the wideband channel and on the narrowband channel. They are estimated on the basis of the parameters C/N0 through the formulae:
 
σ wb   2 =1/( C/N   0 WB )   (7)
 
σ Nb   2 =1/( C/N   0 NB )   (8)
 
     The weightings which are thus applied make it possible to optimize the ratio C/N 0  in the loop. Specifically, the level of the signal is the same in the two channels by assumption, but the scrambling levels are different because the DME interference effect depends on the bandwidth of the channel as was shown above. 
     To obtain these parameters σ, a statistical estimate is made, with the aid of a device (not represented), over samples of the signals Ip and Iq tapped off at a relatively low frequency, 1 kHz for example, carrying this estimate over several seconds, this being possible since the DME signals are constant over distances of several kilometers on account of the intrinsic characteristics of the DME system. 
     The signal thus obtained at the output of the discriminator  220  is filtered in a conventional manner by a filter  221 , then applied to the input of the oscillator  201 , this closing the phase loop. 
     The code discriminator uses signals IΔ and QΔ optimized on the basis of the signals I ΔWB , Q ΔWB , I ΔNB  and Q ΔNB  by applying the formulae:
 
 I   Δopt =( aI   ΔWB   +bI   ΔNB )   (9)
 
 Q   66 opt =( aQ   ΔWB   +bQ   ΔNB )   (10)
 
     In these formulae a and b are given by: 
     
       
         
           
             
               
                 
                   a 
                   = 
                   
                     
                       
                         
                           ρ 
                           WB 
                           2 
                         
                         ⁢ 
                         
                           σ 
                           NB 
                           2 
                         
                       
                       
                         ( 
                         
                           
                             
                               ρ 
                               NB 
                               2 
                             
                             ⁢ 
                             
                               σ 
                               WB 
                               2 
                             
                           
                           + 
                           
                             
                               ρ 
                               WB 
                               2 
                             
                             ⁢ 
                             
                               σ 
                               NB 
                               2 
                             
                           
                         
                         ) 
                       
                     
                     × 
                     
                       1 
                       
                         S 
                         WB 
                       
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
             
               
                 
                   b 
                   = 
                   
                     
                       
                         
                           ρ 
                           NB 
                           2 
                         
                         ⁢ 
                         
                           σ 
                           WB 
                           2 
                         
                       
                       
                         ( 
                         
                           
                             
                               ρ 
                               NB 
                               2 
                             
                             ⁢ 
                             
                               σ 
                               WB 
                               2 
                             
                           
                           + 
                           
                             
                               ρ 
                               WB 
                               2 
                             
                             ⁢ 
                             
                               σ 
                               NB 
                               2 
                             
                           
                         
                         ) 
                       
                     
                     × 
                     
                       1 
                       
                         S 
                         NB 
                       
                     
                   
                 
               
               
                 
                   ( 
                   12 
                   ) 
                 
               
             
           
         
       
     
     In these formulae the terms σ WB  and σ NB  were defined above by formulae (7) and (8). The terms S WB  and S MB  correspond for their part to the slopes of the response curves of the discriminators  220  and  222 , these slopes being measured at the central level of this curve, which corresponds in principle to the response around the carrier frequency. 
     ρ WB  and ρ MB  are the normalized slopes of the discriminator, respectively on the wideband and narrowband channels. To obtain this normalized slope, the corresponding slope S is divided by the factor σ 66 , representative of the standard deviation of the noise on the signals IΔ and QΔ, given by the formula:
 
σ 66 =[1− A (2 d )]  (13)
 
     In this formula, A designates the normalized autocorrelation function of the filtered spreading code, such that A(0)=1, and 2d is the value of the deviation between the early code and the late code used to generate the Δ code. 
     The signal at the output of the discriminator  222  is then filtered very conventionally by a filter  223  and applied to the input of the code oscillator  208 , this closing the code loop. 
     The robustness of the carrier loop is the same as that of the narrowband channel taken in isolation. When the ratio C/N 0  on the wideband channel becomes too low on account of DME interference, the overall ratio C/N 0  of the device tends toward that of the narrowband channel taken in isolation. The carrier loop tracking threshold is therefore actually the threshold of the narrowband channel. 
     The accuracy of the code loop is almost that of the wideband channel taken in isolation when the percentage of time occupancy by the interference is lower than 95%. When this percentage approaches 100%, the accuracy of the code then becomes that of the narrowband channel taken in isolation. 
     Between these two extremities, a synergy is noted stemming from the fact that the robustness on the carrier of one of the channels makes it possible to improve the detection of the code on the other. Stated otherwise the system according to the invention benefits from the robustness of the narrowband channel and the accuracy of the wideband channel when the level of interference remains relatively reasonable. 
     It will also be emphasized that the invention is particularly simple since it uses only a single carrier loop and only a single phase loop, and that it needs no logic detection system to know which of the two channels is to be used as a function of the circumstances. 
     It will be readily seen by one of ordinary skill in the art that the present invention fulfills all of the objects set forth above. After reading the, foregoing specification, one of ordinary skill will be able to affect various changes, substitutions of equivalents and various other aspects of the invention as broadly disclosed herein. It is therefore intended that the protection granted hereon be limited only by the definition contained in the appended claims and equivalents thereof.