Page 524 - Fiber Optic Communications Fund
P. 524

Digital Signal Processing                                                          505


           From Eq. (11.35), we find
                                           {             }
                                                kN
                                     −1        ∑       M       ̂
                                       Arg           (̃y )  ≡ Δ ≅Δ .                   (11.36)
                                                                      l
                                                                k
                                                      l
                                     M
                                             l=(k−1)N+1
                                                                                       ̂
                  ̂
           Here, Δ is the phase estimate of the kth block. The signal sample ̃y is multiplied by exp (iΔ ) to obtain
                                                                                        k
                                                                  l
                  k
           the estimate of the transmitted signal,
                                                          ̂
                                              ̂ x = ̃y exp (iΔ ).                        (11.37)
                                               l   l       k
           The computation in each block can be carried out using separate signal processors. Finally, the signal samples
           in each block are combined using a multiplexer to obtain the serial data. The block size should be chosen
                                                          ′
           carefully. If N is too small, the impact of the noise term n in Eq. (11.32) can not be ignored. If N is too
                                                          l
           large, the laser phase may drift and Δ may not remain constant within each block. The block size should
                                          l
           be optimized based on the laser linewidth.
           11.5.1  Phase Unwrapping
           The function Arg() in Eq. (11.36) can not distinguish between phases that differ by 2 and it returns the results
           in the interval [−, ]. If the phase is  + , > 0, the function Arg() returns a phase of − + . This is known
           as phase wrapping and it could lead to symbol errors. Special techniques have to be used to unwrap phases.
           Consider the following example: suppose Δ in the current block k is ≪ , > 0, and let Δ be roughly
                                               l                                       l
                                                          ̂
           constant over the block. From Eq. (11.36), it follows that Δ = .Now,let Δ of the next block, k + 1jump
                                                          k              l
           by ∕M, i.e., Δ of the (k + 1)th block is  + ∕M. From Eq. (11.36) for the (k + 1)th block, we find
                        l
                             {                             }
                               (k+1)N
                       −1       ∑         [  (           )]    −1
                         Arg        A exp −i M + M +   =    Arg{exp [−i(M + )]}
                                                l
                                     l
                       M                                        M
                               l=kN+1
                                                               M − 
                                                             =
                                                                  M
                                                             =Δ ̂  .                     (11.38)
                                                                  k+1
           Clearly, the estimated phase  − ∕M is different from the actual phase,  + ∕M. This is because of
           the phase wrapping done by the function of Arg(). Phase wrapping in the context of coherent optical
                                                                                             ̂
           communication has been studied in Refs. [14, 15]. Let the carrier phase prior to the unwrapping be Δ .If
                                                                                              k
           we add 2∕M to Δ ̂  , the phase for the (k + 1)th block after the phase unwrapping is
                           k+1
                                        Δ k+1  =Δ ̂ k+1  +  2  =  + ∕M,          (11.39)
                                                       M
           which is actually the phase of the (k + 1)th block. In general, an integral multiple of 2∕M is added to the
           carrier phase. The carrier phase after the phase unwrapping can be written as
                                              ̂
                                      Δ =Δ + m2∕M,   k = 1, 2,...,K,                 (11.40)
                                               k
                                         k
           where                                (                 )
                                                      Δ k−1  −Δ ̂ k
                                       m = Floor 0.5 +              .                      (11.41)
                                                         2∕M
           Here, Floor() returns the nearest integer toward −∞. Suppose the phases of the kth and (k − 1)th block are both
           . In this case, m = 0 and the phase unwrapping block of the phase estimator Eq. (11.40) does not add 2∕M.
   519   520   521   522   523   524   525   526   527   528   529