Abstract:
The present invention provides a method and circuits for digital demodulation of FM and PM modulated signals. In both cases a direct digital phase digitizer is used to obtain the instantaneous phases of the input signal. Digital signal processing circuits comprising only of registers, adders and subtractors, is used to extract the modulating signals from the instantaneous phase information.

Description:
FIELD OF THE INVENTION  
       [0001]     The present invention relates generally to signal demodulators, and particularly to decoders of Phase and Frequency modulations.  
       BACKGROUND OF THE INVENTION  
       [0002]     Phase and frequency modulations are widely used in communication, and the decoding of such modulated signals has been the subject of many articles and inventions.  
         [0003]     In phase modulation, the instantaneous phase deviation of the modulated signal from its unmodulated value is proportional to the instantaneous amplitude of the modulating signal. For a general modulating signal v m (t), the instantaneous phase deviation is, Θ(t)=k Θ v m (t), wherein k Θ  is the phase deviation constant in radians per volts. For V m  defined as the maximum value of |v m (t)|, it is convenient to define a “normalized” v(t)=[v m (t)]/V m , and in this notation, Θ(t)=k Θ V m v(t), and the maximum phase shift, k Θ V m =ΔΘ=m p  is called the modulation index for phase modulation.  
         [0004]     In terms of m p  the phase modulated signal is written as F pm (t)=A cos[ω c t+m p v(t)], and the instantaneous phase deviation is Θ(t)=m p v(t), radians.  
         [0005]     The instantaneous frequency of the modulated signal is  
         ω   ⁡     (   t   )       =         ω   c     +       ⅆ     Θ   ⁡     (   t   )           ⅆ   t         =       ω   c     +       m   p     ⁢         ⅆ     v   ⁡     (   t   )           ⅆ   t       .               
 
         [0006]     Frequency modulation results when the deviation δω of the instantaneous frequency ω(t) from the carrier frequency ω c  is directly proportional to the instantaneous amplitude of the modulating voltage.  
         [0007]     Since  
         ω   ⁡     (   t   )       =         ⅆ   ϕ       ⅆ   t       =       ω   r     +       ⅆ     Θ   ⁡     (   t   )           ⅆ   t               
 
 the frequency deviation δω of ω(t) from ω c  is given by  
         δω   ⁡     (   t   )       =         ω   ⁡     (   t   )       -     ω   c       =         ⅆ     Θ   ⁡     (   t   )           ⅆ   t       .           
 
         [0009]     In frequency modulation δω(t) is proportional to the modulating voltage v m (t), as δω(t)=k m v m (t), in which k m  is the sensitivity of the modulator in rad/s/V. Since Θ(t), and δω(t) are related, as shown above, then Θ(t)=∫ 0   1 k m v m (t)dt+Θ(0), and assuming Θ(0)=0,  
         [0010]     And F FM =A cos└ωω c t+k m ∫ 0   1 v m (t)dt┘.  
         [0011]     Demodulators of Phase and Frequency modulated signals, are used to extract the modulating signal v m (t) from the modulated signal. Such demodulators typically comprise of tuned circuits, or phase locked loops, in which a phase or frequency deviation causes a change in the output voltage, which is directly related to the magnitude of the phase or frequency deviation.  
         [0012]     A simple FM demodulator is based on an LC tank resonator circuit. In resonance the amplitude versus the frequency response of the tuned circuit has the shape of a bell, as shown in  FIG. 1 . The LC tank is tuned to resonate, such that the center frequency (carrier frequency) of the FM modulated signal, and therefore any change in the signal frequency causes a change in the voltage on the resonator.  
         [0013]      FIG. 2 , shows the most commonly used FM demodulator. To analyze this circuit, consider the voltage relationship in the coupled circuit L 1 C 1 -L 2 C 2 , in  FIG. 2 , and  FIG. 3 . If the impedance coupled into the primary circuit is negligible in comparison with the primary self impedance, then the primary current I 1  will be  
         I   1     =       V   1       j   ⁢           ⁢   ω   ⁢           ⁢     L   1               
         [0014]     The voltage induced in the secondary by I 1  is jωMI 1 ; at the secondary resonance frequency ω c  the secondary current will be:  
         I   2     =         j   ⁢           ⁢     ω   c     ⁢     MI   1         R   2       =       M     L   1       ⁢       V   1       R   2               
 
         [0015]     The voltage across the capacitor is: 2V 2 =jω c C 2 I 2 , and therefore  
         V   2     =         j   ⁢           ⁢     ω   c     ⁢     C   2     ⁢   M       2   ⁢     L   1     ⁢     R   2         ⁢     V   1           
 
         [0016]     At resonance V 2  leads V 1  by 90°. At frequency ω slightly different from ω c ,  
               V   2     =       ⁢           j   ⁢           ⁢   ω   ⁢           ⁢     C   2     ⁢   M       2   ⁢     L   1     ⁢     R   2         ×       V   1       1   +       (     1     R   2       )     ⁡     [       j   ⁢           ⁢   ω   ⁢           ⁢     L   2       +     (       1   /   j     ⁢           ⁢   ω   ⁢           ⁢     C   2       )       ]             ≈                 ≈       ⁢         j   ⁢           ⁢   ω   ⁢           ⁢     C   2     ⁢     MV   1         2   ⁢     L   1     ⁢     R   2         ×     {     1   -       1     R   2       ⁡     [       j   ⁢           ⁢   ω   ⁢           ⁢     L   2       +     1     j   ⁢           ⁢     ωC   2           ]         }                       ⁢   □             
 
 and the phase angle between V 1  and V 2  is given by:  
         Arg   ⁡     (       V   2       V   1       )       =       Arg   ⁡     (     j   +       Δω   ⁢           ⁢     L   2         R   2       -     1     j   ⁢           ⁢   ω   ⁢           ⁢     C   2     ⁢     R   2           )       =       cot     -   1       ⁡     (         ω   ⁢           ⁢     L   2         R   c       -     1     ω   ⁢           ⁢     C   2     ⁢     R   2           )             
 
         [0018]     Then, if ω=ω c +Δω, the phase angle formula is reduced to  
         Arg   ⁡     (       V   2       V   1       )       =         cot     -   1       ⁡     (       Δω   ⁢           ⁢     L   2         R   2       )       =             cot     -   1       ⁡     (       Δω   ⁢           ⁢   Q       ω   c       )       wherein     ⁢           ⁢   Q     =       ω   ⁢           ⁢     L   2         R   2               
 
 which is the Q of the secondary resonant circuit. 
 
         [0020]     For Δω&gt;0, if ω&gt;ω c  then the angle is &lt;90°, and if ω&lt;ω c  then the angle is &gt;90°, as shown in  FIG. 4 .  
         [0021]     Referring to  FIG. 2 , assuming that C o  and C c  are RF shorts, and that RFC is an RF open, the point “A” in  FIG. 2 , may be considered RF ground, and the RF voltage V a ′ applied to the top diode, and the R o C o  network is V a ′=V 1 +V 2 . The RF voltage V b ′ applied to the bottom diode is V b ′=V 1 −V 2 . The phase relationship between V 1  and V 2  is shown in  FIG. 4 .  
         [0022]      FIG. 5 , illustrates the variation of V a ′, and V b ′ as a function of variations in frequency. The basic relationships are: 
    V a ′=V b ′when ω=ω c ,     |V a ′|&lt;|V b ′| when ω&lt;ω c ,     V a ′/&gt;|V b ′| when ω&gt;ω c .      
         [0026]     The demodulated output from the FM demodulator is V d =|V a ′|−|V b ′|.  
         [0027]     Since the phase deviation is the derivative of the frequency modulation, the phase demodulation is obtained by differentiating the frequency demodulator output.  
         [0028]     An alternative method of FM demodulation, uses Phase Locked Loops (PLL), as the means to track frequency and phase deviations, and extract the modulating signal.  
         [0029]     Referring to  FIG. 6 , showing a typical phase locked loop, comprised of a phase/frequency detector ( 61 ), a loop filter ( 62 ), and a voltage controlled oscillator ( 63 ). The control mechanism of the loop tries to maintain the loop in the locked state, wherein the signal at the output of the VCO ( 65 ) matches the phase and frequency of the reference input ( 66 ). As the reference input ( 66 ), which is the modulated signal, changes its phase or frequency, the loop follows those changes in order to stay locked. The voltage output of the loop filter ( 64 ) is directly proportional to the modulating signal.  
         [0030]     In this invention, some new, completely digital, decoders for phase and frequency modulated signals are described. These decoders do not contain any tuned circuits, nor do they employ any phase locked loops. These decoders are inherently wideband, and the bandwidth of their operation is determined mainly by the frequency of the clock signal used.  
         [0031]     This invention uses a direct phase sampler (DPS) to provide numerical information identifying the instantaneous phase of the input signal.  FIG. 7 , shows a block diagram of a direct phase sampler. The DPS receives the modulated RF signal at its input, and on each clock transition it produces a digital number, indicating the instantaneous phase of the input signal at the time of the clock transition. Such phase samplers have been described before in other patens, and is not the subject of this invention.  
         [0032]     Referring to  FIG. 8 , which shows a block diagram of an embodiment of a PM demodulator. The output of the DPS ( 100 ) is applied to a differencing circuit ( 10 ), which subtracts the phase of the input signal at the time of the last clock transition, from the phase of the input signal at the time of the previous clock transition, thus providing the phase difference between any two clocks. The output of the differencing circuit is applied both to a “p” deep running averager ( 20 ), and to the “A” input of a subtractor ( 30 ). The averager calculates the average of phase differences over the last consecutive “p” clock periods. The requirement on the averager length is that if t c  is the sampling clock period, and t m(max)  is the period of the lowest frequency in the PM modulating bandwidth, then  
       p   &gt;     10   ⁢       t     m   ⁡     (   max   )           t   c             
 
         [0033]     The subtractor ( 30 ), which follows the averager ( 20 ), subtracts the average phase difference calculated by the averager, from the instantaneous phase difference calculated by the differencing circuit. The resulting output is the variation in phase difference from clock transition to clock transition. A sine lookup table ( 40 ), which follows the subtractor ( 30 ), converts the phase variations information generated by the subtractor ( 30 ), into an amplitude voltage output ( 45 ), which is essentially the demodulation of the PM modulated input signal.  
         [0034]      FIG. 10 , shows a block diagram of an FM demodulator. This demodulator is very similar to the PM demodulator, except that the output of the differencing circuit ( 10 ) is not applied directly to the “A” input of the subtractor ( 30 ), but instead is connected to a “q” deep running averager ( 50 ), who&#39;s output connects to the “A” input of the subtractor ( 30 ). The requirements for the size of “q” depend on the application. In general, “q” is much smaller than “p”. However, the size of “q” effects the decoder sensitivity, and its fidelity or linearity. The smaller “q” is, the higher is the demodulator&#39;s sensitivity, but the lower is its fidelity. 
     
    
     DESCRIPTION OF THE DRAWINGS  
       [0035]      FIG. 1 , shows the amplitude versus frequency response of an LC tank resonator circuit.  
         [0036]      FIG. 2 , shows the schematic of a typical “analog” frequency discriminator.  
         [0037]      FIG. 3 , shows a circuit to illustrate the phase relationship in the discriminator. The resistors represent the internal resistance of the inductors.  
         [0038]      FIG. 4 , shows the phase relationship in the discriminator.  
         [0039]      FIG. 5  shows the voltages on the diodes in the discriminator.  
         [0040]      FIG. 6 , shows a Phase Locked Loop FM/PM discriminator.  
         [0041]      FIG. 7 , shows a block diagram of direct digital phase sampler.  
         [0042]      FIG. 8 , shows a block diagram of a direct digital PM demodulator.  
         [0043]      FIG. 9 , shows an embodiment of a PM demodulator.  
         [0044]      FIG. 10 , shows a block diagram of an FM demodulator.  
         [0045]      FIG. 11 , shows an embodiment of an FM demodulator.  
         [0046]      FIG. 12 , shows an embodiment of a receiver comprising a direct digital demodulator.  
         [0047]      FIG. 13 , shows the waveforms involved in binary to Grey code conversion.  
         [0048]      FIG. 14 , shows an embodiment of a phase to amplitude converter.  
         [0049]      FIG. 15 , shows waveforms involved with phase to amplitude conversion.  
     
    
     DESCRIPTION OF THE INVENTION  
       [0050]      FIG. 9 , shows the details of an embodiment of a PM demodulator. The register ( 12 ), the subtractor ( 13 ), and the register ( 14 ), comprise the phase differencing circuit ( 10 ). The output of the direct digital phase sampler ( 101 ) Θ k (t) is applied simultaneously to the input of the register ( 12 ), and the “A” input of the subtractor ( 13 ). The output of the register ( 13 ) lags behind the input to that register by one clock period, and thus the input “B” to the subtractor ( 13 ) Θ k−1 (t) lags one clock period behind the input “A” to the subtractor ( 13 ) Θ k (t). As a result, the output of the subtractor ( 13 ), which is the difference between inputs “A” and “B” ΔΘ k =Θ k (t)−Θ k−1 (t), is actually the change in the phase of the input signal ( 102 ), over one clock period, which is the instantaneous frequency of the input signal ( 102 ).  
         [0051]     The adder ( 21 ), the “p” deep shift register ( 22 ), the register ( 23 ), the subtractor ( 24 ), and the register ( 25 ), comprise the averager ( 20 ). Assuming that initially all registers and shift registers outputs are “0”. The output of the shift register ( 22 ) will remain “0” for at least “p” clock cycles, as any non “0” data at the input to the shift register ( 22 ) propagates through the shift register in “p” clock periods. The adder ( 21 ) adds new data ΔΘ k  coming from the differencing circuit ( 10 ) [register ( 14 )], with data Φ k  coming out of the subtractor ( 24 ) via the register ( 25 ). While the output of the shift register ( 22 ) ΔΘ k−p  is “0” for “p” clock cycles, the output of the subtractor ( 24 ) Φ k =Λ k −ΔΘ k−p  is the same as the data at its “B” input Λ k . As a result, for the first “p” clock cycles, the adder ( 21 ) accumulates all the phase differences generated by the subtractor ( 13 ) Λ k+1 =Φ k +ΛΘ k =Λ k −ΔΘ k−p +ΔΘ k . The divider ( 26 ), which follows the register ( 25 ) divides the output from the subtractor Φ k  by p, to yield the running average  
         Δ   ⁢           ⁢     Θ   ACG       =         Σ   ⁡     (       Δ   ⁢           ⁢     Θ   k       -     Δ   ⁢           ⁢     Θ     k   -   p           )       p     .         
 
 If p is selected p such that p=2″, then the division can be accomplished by simply discarding the n least significant bits at the output of the averager. The output of the averager is the average phase difference for any clock period. Dividing the average phase difference by the clock period yields the average, or center frequency of the input signal  
         F   c     =         Δ   ⁢           ⁢     Θ   AVG         t   c       .         
 
         [0053]     The output of the averager ( 20 ) ΔΘ AVG  is subtracted by the subtractor ( 31 ) from the instantaneous phase difference ΔΘ k  to yield the phase deviation θ k =ΔΘ AVG −ΔΘ k .  
         [0054]     A sine lookup table ( 41 ) followed by a digital to analog converter ( 42 ), is a convenient way to convert phase information to amplitude information for the demodulator output.  
         [0055]      FIG. 11 , shows an embodiment of an FM demodulator. This demodulator is very similar to the PM demodulator. It uses the same phase sampler ( 101 ), the same differencing circuit ( 10 ), the same averager ( 20 ), the same subtractor ( 30 ) and the same sine lookup table ( 40 ). The only difference is that the input “A” of the subtractor ( 31 ) is not connected to ΔΘ k , the output of the differencing circuit ( 10 ), but instead, an averager ( 50 ), which is comprised of the subtractor ( 51 ), q deep shift register ( 52 ), register ( 53 ), subtractor ( 54 ), register ( 55 ), and divider ( 56 ). The operation of this averager ( 50 ), is similar to the operation of the averager ( 20 ), with the only difference in the length of the shift register which is q&lt;&lt;p.  
         [0056]     In the FM demodulator, the averager ( 50 ), having a shift register much shorter than that of the other averager ( 20 ), produces the instantaneous deviated frequency  
           F   d     =       Δ   ⁢           ⁢     Ξ   AVG         t   c         ,       
 
 wherein ΔΞ AVG  is the average phase difference per clock period output of the averager ( 50 ). The subtractor ( 31 ) subtracts the instantaneous deviated frequency from the center frequency, resulting in the frequency deviation ΔF=F d −F c . The sine lookup table ( 40 ) converts the phase information into amplitude information, to complete the demodulation process. 
 
         [0058]     An alternative method and circuit for converting phase information into voltage amplitude is shown in  FIGS. 13, 14 , and  15 .  
         [0059]     In digital presentation of numbers the bits are assigned values which are power of 2 wherein the least significant bit is assigned the value of 2 0 , the next bit is 2 1 , etc. In binary code presentation the order of values in 4 bits is: 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15. In a Grey code on the other hand, the order of values is: 0, 1, 3, 2, 6, 7, 5, 4, 12, 13, 15, 14, 10, 11, 9, 8. When these values are presented in 4 bits waveforms, the resulting Grey code waveforms are symmetrical, unlike the binary code, which is non-symmetrical. This symmetry feature of Grey code waveforms enables their use in phase to amplitude conversion.  
         [0060]     Obtaining Grey code out of binary code is a straight forward process of EXORing pairs of bits in the form G n =B n ⊕B n+1 . The most significant bit (MSB) in the Grey code is the same as the MSB in binary code, as shown in  FIG. 13 .  
         [0061]      FIG. 14 , shows an embodiment of a phase to amplitude converter. Following the conversion of binary to Grey code ( 1 ), the Grey code bits are further EXORed ( 2 ) and applied to amplifiers with output spanning between a positive supply rail (+V) to a negative supply rail (−V). The resulting voltage waveforms are applied to a resistive network ( 3 ). In the resistive network, currents are summed together on the output resistor (R(out)) to generate a sinewave approximation waveform output ( 4 ).  FIG. 15 , shows the waveforms in converting phase information presented in Grey code into a sinewave approximation waveform.  
         [0062]      FIG. 12 , shows an embodiment of a receiver utilizing a digital demodulator. The input signal ( 212 ) is split in the power splitter ( 201 ) into two equal but lower power version of the input signal. These signals are applied to two RF mixers ( 203 , and  204 ). A local oscillator ( 211 ) generates a signal, at a frequency which when added to (or subtracted from) the input signal ( 212 ) yield a frequency which is in the center of the band of the bandpass filters ( 204 ). The output of the oscillator is passed through a Hybrid Coupler ( 202 ), which splits the oscillator&#39;s output into two signals with equal amplitudes but with a 90° phase relationship. The two signals generated by the hybrid coupler ( 202 ) are applied to the mixers ( 203 , and  204 ) at their LO ports. As a result, the mixers ( 203 , and  204 ) outputs are two signals at a frequency of the center of the bandpass filters, and with a 90° phase relationship between them. These signals pass through the bandpass filters, and are applied to the direct phase digitizer as “I” and “Q” ( 209 , and  210  respectively). The digital processing is operating with a clock (not shown). On every clock cycle, the digitizer ( 206 ) generates a data output representing the phase of the input signal at the time of the clock transition. The digital demodulator ( 207 ) receives the data output from the digitizer ( 206 ) and extracts the modulating signal from the modulated input signal ( 212 ).