# Patent

## Title
Method and apparatus for reduced-state viterbi detection in a read channel of a magnetic recording system

## Abstract

A method and apparatus are disclosed for improving the maximum data rate of reduced-state Viterbi detectors with local feedback in magnetic recording systems. A read channel signal is processed in a magnetic recording device by precomputing branch metrics, intersymbol interference estimates or intersymbol interference-free signal estimates for speculative sequences of one or more channel symbols; selecting one of the precomputed values based on at least one decision from at least one corresponding state; and selecting a path having a best path metric for a given state.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application is a continuation application of U.S. patent application Ser. No. 10/853,090, filed May 25, 2004 now abandoned and a continuation-in-part of pending U.S. application Ser. No. 11/234,446, filed Sep. 26, 2005, which is a divisional of U.S. application Ser. No. 09/834,668, filed Apr. 13, 2001, now U.S. Pat. No. 7,000,175, each incorporated by reference herein.


FIELD OF THE INVENTION
The present invention relates generally to equalization, detection and decoding techniques and, more particularly, to sequence estimation techniques with reduced complexity.
BACKGROUND OF THE INVENTION
A magnetic recording read channel converts an analog read channel into an estimate of the user data recorded on a magnetic medium. Read heads and magnetic media introduce noise and other distortions into the read signal. As the information densities in magnetic recording increase, the intersymbol interference (ISI) becomes more severe as well, (i.e., the channel impulse response becomes longer). In read channel chips, a Viterbi detector is typically used to detect the read data bits in the presence of intersymbol interference and noise. When the channel impulse response is long, however, the hardware complexity associated with the Viterbi detector becomes prohibitively large, as the number of states considered by the Viterbi detector grows exponentially with the length of the channel impulse response. A number of techniques have been proposed or suggested for reducing the complexity of Viterbi detectors.
For example, the hardware complexity of the Viterbi detector can be reduced by using a reduced-state trellis that considers only a shortened impulse response, and canceling intersymbol interference due to the tail of the impulse response for each state by using past survivor symbols as local feedback. See, e.g., J. W. M. Bergmans, “Digital Baseband Transmission and Recording,” Kluwer Academic Publishers, 326 (1996) or U.S. Pat. No. 6,690,754, issued to Haratsch et al., entitled “Method and Apparatus for Reducing the Computational Complexity and Relaxing the Critical Path of Reduced-State Sequence Estimation (RSSE) Techniques,” incorporated by reference herein.
The error rate performance of reduced-state Viterbi detectors with local feedback can approach the performance of full-state Viterbi detectors without local feedback that implement maximum likelihood sequence estimation (MLSE). The maximum achievable data rate of a Viterbi detector implementation with local feedback, however, is considerably lower compared to a Viterbi detector implementation without local feedback, as significantly more operations have to be performed within one clock period. A need therefore exists for a method and apparatus for performing reduced-state Viterbi detection with local feedback at the high data rates that are required by evolving high-end storage applications.
SUMMARY OF THE INVENTION
Generally, a method and apparatus are disclosed for improving the maximum data rate of reduced-state Viterbi detectors with local feedback in magnetic recording systems. The maximum data rate that may be achieved by the disclosed reduced-state Viterbi detectors is improved by precomputing a number of candidate branch metrics and performing pipelined selection of one of the precomputed values. Precomputing the branch metrics for possible symbol combinations in the channel memory makes it possible to shorten the critical path. In an alternative embodiment, intersymbol interference estimates or inter symbol interference-free signal estimates are precomputed, and branch metric are computed based on the selected precomputed values.
A read channel signal is processed in a magnetic recording device by precomputing branch metrics, intersymbol interference estimates or intersymbol interference-free signal estimates for speculative sequences of one or more channel symbols; selecting one of the precomputed values based on at least one decision from at least one corresponding state; and selecting a path having a best path metric for a given state.
A more complete understanding of the present invention, as well as further features and advantages of the present invention, will be obtained by reference to the following detailed description and drawings.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic block diagram of a conventional system model for a baseband communications channel with ISI and additive noise;
FIG. 2 illustrates a trellis diagram for a channel with memory L=1;
FIG. 3 illustrates a trellis diagram for a channel having a memory L=4;
FIG. 4 illustrates a reduced-state trellis diagram corresponding to the full-state trellis of FIG. 3, for a channel having a memory L=4 and a shortened channel memory K=1;
FIG. 5 is a schematic block diagram for an exemplary conventional reduced-state Viterbi detector with local feedback;
FIG. 6 illustrates a detailed state-parallel implementation of a reduced-state Viterbi detector with local feedback corresponding to the trellis of FIG. 4;
FIG. 7 is a schematic block diagram of a reduced-state Viterbi detector that incorporates precomputation of the branch metrics;
FIG. 8 is a schematic block diagram showing the selection of a precomputed branch metric by the reduced-state Viterbi detector of FIG. 7 using survivor symbols;
FIG. 9 is a schematic block diagram of a reduced-state Viterbi detector that incorporates precomputation of the ISI-free signal estimates;
FIG. 10 is a schematic block diagram showing the selection of a precomputed ISI-free signal estimate by the reduced-state Viterbi detector of FIG. 9 using survivor symbols;
FIG. 11 is a schematic block diagram showing the selection of a precomputed intersymbol interference estimate using survivor symbols;
FIG. 12 is a schematic block diagram of a reduced-state Viterbi detector incorporating pipelining of the branch metric selection;
FIG. 13 is a schematic block diagram showing the pipelined selection of a branch metric by the reduced-state Viterbi detector of FIG. 12 using ACS decisions;
FIG. 14 is a schematic block diagram of a reduced-state Viterbi detector that incorporates pipelining of the ISI-free signal estimate selection;
FIG. 15 is a schematic block diagram showing the pipelined selection of an ISI-free signal estimate by the reduced-state Viterbi detector of FIG. 14 using ACS decisions;
FIG. 16 is a schematic block diagram showing the pipelined selection of an inter symbol interference estimate using ACS decisions;
FIG. 17 is a functional block diagram of the ACS operation performed in FIGS. 8, 10, 11, 13, 15 and 16;
FIG. 18 is a functional block diagram of a read channel detector that implements noise-predictive data detection and uses one of the reduced-state Viterbi detectors of FIGS. 7-16 incorporating features of the invention; and
FIG. 19 is a functional block diagram of a read channel detector that implements signal-dependent noise-predictive data detection and uses one of the reduced-state Viterbi detectors of FIGS. 7-16 incorporating features of the invention.



DETAILED DESCRIPTION
The present invention increases the maximum data rate that may be achieved by reduced-state Viterbi detectors. According to one aspect of the invention, branch metrics, ISI-free signal estimates or ISI estimates are precomputed, and the correct values are selected based on survivor symbols or ACS decisions. In this manner, the computations of ISI estimates, ISI-free signal estimates or branch metrics are removed from the critical path. According to another aspect of the invention, branch metrics, ISI-free signal estimates or ISI estimates can be selected in a pipelined fashion using a multiplexer network structure that corresponds to the structure of the trellis considered by the detector.
For a detailed discussion of reduced-state Viterbi detection with local feedback, which is also known as Reduced-State Sequence Estimation (RSSE), (Delayed) Decision-Feedback Sequence Estimation (DFSE), and Parallel Decision-Feedback Equalization (PDFE), see, for example, U.S. Pat. No. 6,690,754 to Haratsch et al., entitled “Method and Apparatus for Reducing the Computational Complexity and Relaxing the Critical Path of Reduced-State Sequence Estimation (RSSE) Techniques,” incorporated by reference herein, and the references cited therein. See also, Lee and Messerschmidt, “Digital Communication,” Kluwer Academic Publishers, 2nd ed. (1994).
FIG. 1 is a schematic block diagram of a conventional system model for a baseband communications channel 100 with ISI and additive noise. While the exemplary embodiment is discussed in the context of baseband communications, the techniques discussed herein can also be applied to passband communications systems, as would be apparent to a person of ordinary skill in the art. Further, while it is assumed that trellis-coded modulation (TCM) is not employed for ease of illustration, the disclosed techniques generalize to communication systems using TCM or other modulation schemes.
The modulator 110 maps an information symbol bn into a channel symbol an. For ease of illustration, it is assumed that the number of information bits per information symbol is one. In other words, the information symbol bn is equivalent to a single information bit bn. The modulator 110 maps an information symbol bn to a two-level channel symbol an according to following rule:




  
    
      
        
          a
          n
        
        =
        
          {
          
            
              
                
                  
                    -
                    1
                  
                  ,
                
              
              
                
                  
                    b
                    n
                  
                  =
                  0
                
              
            
            
              
                
                  1
                  ,
                
              
              
                
                  
                    b
                    n
                  
                  =
                  1
                
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




The techniques discussed herein can easily be applied to other modulation schemes and more than two signal levels. For a discussion of reduced-state Viterbi detection for an exemplary modulation scheme with five signal levels, see, U.S. patent application Ser. No. 09/471,920, entitled, “Method and Apparatus for Shortening the Critical Path of Reduced Complexity Sequence Estimation Techniques,” incorporated by reference herein.
The ISI channel 100 is modeled as an FIR filter, and the channel output at time n is given by




  
    
      
        
          
            r
            n
          
          =
          
            
              
                z
                n
              
              +
              
                w
                n
              
            
            =
            
              
                
                  ∑
                  
                    i
                    =
                    0
                  
                  L
                
                ⁢
                
                  
                    f
                    i
                  
                  ·
                  
                    a
                    
                      n
                      -
                      i
                    
                  
                
              
              +
              
                w
                n
              
            
          
        
        ,
      
    
    
      
        (
        2
        )
      
    
  




where zn is the ISI channel output, {fi}, 0≦i≦L are the channel coefficients, L is the channel memory, and wn is noise. The decision of a detector 120 that corresponds to bn is denoted by b′n.

The ISI channel output zn depends on the current channel symbol an and the past L transmitted channel symbols {an−i}, 1≦i≦L. This output can be described as a function of the L past transmitted channel symbols using a finite state machine (FSM) model, where the channel state at time n is defined by

αn=(an−1,an−2, . . . , an−L).  (3)

The channel state is equivalently defined in terms of the L past transmitted information bits:

βn=(bn−1,bn−2, . . . , bn−L).  (4)

It is apparent from equations (3) or (4) that the number of channel states is given by

2L.  (5)

To simplify the notation, the integer value corresponding to the vector (bn−1, . . . , bn−L+1, bn−L) will be used to represent the channel state βn. For example, 0n will stand for βn (0, . . . , 0,0) and 1n will stand for βn=(0, . . . , 0,1).
The FSM process describing the ISI channel 100 can be visualized using a trellis diagram 200, shown in FIG. 2, for a channel with memory L=1. For the considered exemplary uncoded channel model, a trellis state at time n is denoted by σn, and is equal to the channel state, i.e., σn=βn. In FIG. 2, solid lines correspond to survivor paths, dotted lines to discarded transitions, and dashed lines to path extensions. There are two channel states, and two branches corresponding to the information symbols bn=0 and bn=1 leave each state σn to reach respective successor states {σn+1}. It can be seen from equation (5) that the number of channel states grows exponentially with respect to the channel memory
FIG. 2 depicts the operation of the Viterbi algorithm at time step n. At this point, the Viterbi algorithm has already determined the survivor path into state 0n, which corresponds to the surviving state sequence {0n, 1n−1, 0n−2, 1n−3, . . . }. The survivor path into state 1n corresponds in this example to the state sequence {1n, 0n−1, 0n−2, 1n−3, . . . }. Based on these two survivor paths, the Viterbi algorithm decides on the survivor paths into states 0n+1, and 1+1, in the manner described below.
First, the Viterbi algorithm calculates branch metrics for the state transitions from σn to σn+1. For a channel with additive white Gaussian noise, the optimum blanch metric is the Euclidean distance between the received symbol rn and the ideal ISI channel output zn that corresponds to the respective state transition. For a transition from state σn, the branch metric is given by




  
    
      
        
          
            
              λ
              n
            
            ⁡
            
              (
              
                
                  σ
                  n
                
                ,
                
                  a
                  n
                
              
              )
            
          
          =
          
            
              
                (
                
                  
                    r
                    n
                  
                  -
                  
                    z
                    n
                  
                
                )
              
              2
            
            =
            
              
                (
                
                  
                    r
                    n
                  
                  -
                  
                    
                      ∑
                      
                        i
                        =
                        0
                      
                      L
                    
                    ⁢
                    
                      
                        f
                        i
                      
                      ·
                      
                        a
                        
                          n
                          -
                          i
                        
                      
                    
                  
                
                )
              
              2
            
          
        
        ,
      
    
    
      
        (
        6
        )
      
    
  




where an is the channel symbol that is associated with a transition from state σn to a successor state σn+1. The techniques described herein are independent from the way branch metrics are computed, i.e., branch metrics can also by computed by using the absolute value of the difference between the received symbol rn and the ideal ISI channel output zn.

In the trellis 200 of FIG. 2, there ale two path extensions into any state σn+1, e.g., state 0n+1 can be reached from states 0n and 1n. Out of the two path extensions into a particular state σn+1, the Viterbi algorithm keeps only the one extension with the smallest path metric, as it corresponds to the most likely path into this state. The metric for the path that emanates from state σn and leads into σn+1 is calculated by adding the path metric for the preceding state σn, Γn(σn) and the branch metric λn (σn, an) for the transition.
The three operations to determine the best survivor path into a new state σn+1, i.e., adding up corresponding path metrics of predecessor states σn and branch metrics for the extensions into the new state σn+1, comparing the path metrics of these extended sequences, and selecting the extension with the minimum path metric as the survivor sequence for the new state, are referred to as add-compare-select (ACS), which can be described by the following equation:

Γn+1(σn+1)=min (Γn(σn)+λn(σn,an)).  (7)

As previously indicated, the invention can also be applied when branch metrics are computed differently. As known in the art, for certain branch metric definition, the best path into a state is given by the path with the maximum (instead of minimum) path metric. For such cases, the ACS operation described by equation (7) involves a maximum instead of a minimum operation.
In the example of FIG. 2, the two survivor sequences into states 0n and 1n merge into a single path at time step n−2. In general, all survivor paths merge into a single path after some detection delay D with high probability. Thus, information symbols can be uniquely detected from this time step on. Therefore, it is possible to implement the Viterbi algorithm with a fixed detection delay. It is not required to process the whole transmitted sequence before the first information symbols can be detected. Generally, the detection delay D should be approximately five times the memory of the underlying FSM process. For ISI channels, the memory is equal to L. Typically, a good value for D is determined by running error rate simulations for different values of D
FIG. 3 illustrates a trellis 300 describing an ISI channel having a memory L=4. A trellis state at time n is denoted by σn, and for the considered exemplary uncoded channel model, it is equal to the channel state, i.e., σn=βn. There are 16 channel states, and two branches corresponding to the information symbols bn=0 and bn=1 leave each state σn to reach respective successor states {σn+1}.
Reduced-State Viterbi Detection with Local Feedback
As indicated above, the disadvantage of MLSE is that its complexity grows exponentially with the channel memory. Considering fewer states for the detection of the most likely data sequence reduces the required hardware or computational effort. Reduced-state Viterbi Detection with local feedback accomplishes this by merging several separate states into one single reduced state and keeping only one survivor path per reduced state. The ISI that is not considered in the reduced state is cancelled for each reduced-state by using channel symbols from the corresponding survivor path in a local feedback fashion. Reduced-state Viterbi detection with local feedback is also known as “Reduced-State Sequence Estimation (RSSE)”, “(Delayed) Decision-Feedback Sequence Estimation”, “Parallel Decision-Feedback Equalization”, etc.
In the simplest variant of RSSE, a reduced state β′n is obtained by not considering all L information symbols, but only the past K information symbols for the definition of a trellis state:

β′n=(bn−1,bn−2, . . . , bn−K),0≦K≦L,  (8)

where K is referred to as the truncated channel memory. The number of states in the reduced-state trellis is then given by

2K.  (9)

The reduced state β′n does not contain information about the ISI caused by the channel symbols (an−K−1, an−K−2, . . . , an−L). Conceptually, this reduced state is obtained by grouping all original states βn as defined in Equation (4) with the same information symbol sequence (bn−1, bn−2, . . . , bn−K), but different sequences (bn−K−1, bn−K−2, . . . , bn−L) into one single reduced state β′n. Therefore, this reduced state does not make any statement about the ISI associated with the channel coefficients (fK+1, fK+2, . . . , fL). But an estimate for this ISI component can be computed by considering the respective channel symbols from the survivor sequence into this state. The ISI corresponding to a state is not known a-priori as in MLSE, but must be determined at each detection step by using channel symbols from the corresponding survivor path. Let σn denote a state in the reduced-state trellis, i.e., σn=β′n. The ISI estimate un (σn) for a state σn is calculated at time step n as




  
    
      
        
          
            
              u
              n
            
            ⁡
            
              (
              
                σ
                n
              
              )
            
          
          =
          
            
              
                ∑
                
                  i
                  =
                  1
                
                K
              
              ⁢
              
                
                  f
                  i
                
                ·
                
                  a
                  
                    n
                    -
                    i
                  
                
              
            
            +
            
              
                ∑
                
                  i
                  =
                  
                    K
                    +
                    1
                  
                
                L
              
              ⁢
              
                
                  f
                  i
                
                ·
                
                  
                    
                      a
                      ^
                    
                    
                      n
                      -
                      i
                    
                  
                  ⁡
                  
                    (
                    
                      σ
                      n
                    
                    )
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        10
        )
      
    
  




where ân−i(σn) is the channel symbol that corresponds to the survivor sequence into state σn and that is associated with trellis step n−i The first term on the right hand side of equation (10) computes the ISI component that is known a-priori due to the definition of the reduced state in equation (8). The second term on the right hand side of equation (10) is the ISI component caused by channel taps that were ignored in the reduced-state definition of equation (8). This ISI term is calculated at each detection step for a given state by using respective survivor symbols as local feedback.

With the ISI estimate un(σn), the branch metric for the transition that emanates from state σn to reach a successor state σn+1 and corresponds to channel symbol an can be computed as:

λn(σn,an)=(rn−f0·an−un(σn))2.  (11)

As in MLSE, the most likely survivor path into the state σn+1 with the path metric Γn+1(σn+1) among the path extensions from all possible predecessor states {σn} is determined with an ACS operation:




  
    
      
        
          
            Γ
            
              n
              +
              1
            
          
          ⁡
          
            (
            
              σ
              
                n
                +
                1
              
            
            )
          
        
        =
        
          
            min
            
              
                {
                
                  σ
                  n
                
                }
              
              →
              
                σ
                
                  n
                  -
                  1
                
              
            
          
          ⁢
          
            
              (
              
                
                  
                    Γ
                    n
                  
                  ⁡
                  
                    (
                    
                      σ
                      n
                    
                    )
                  
                
                +
                
                  
                    λ
                    n
                  
                  ⁡
                  
                    (
                    
                      
                        σ
                        n
                      
                      ,
                      
                        a
                        n
                      
                    
                    )
                  
                
              
              )
            
            .
          
        
      
    
    
      
        (
        12
        )
      
    
  




The version of RSSE where a reduced state is defined by considering just a truncated channel memory as in equation (8) is referred to as (Delayed) Decision-Feedback Sequence Estimation (DFSE), described, for example, in A. Duel-Hallen and C. Heegard, “Delayed Decision-Feedback Sequence Estimation,” IEEE Transaction on Communications, 428-436 (May 1989). A reduced-state trellis can also be constructed by applying set partitioning principles to the channel symbol alphabet, as suggested in M. V. Eyuboglu and S. U. Qureshi, “Reduced-State Sequence Estimation With Set Partitioning and Decision-Feedback,” IEEE Transactions on Communications, 13-20 (January 1988). Recently, even more general rules for the definition of reduced states were given in B. E. Spinnler and J. Huber, “Design of Hyper States for Reduced-State Sequence Estimation,”, AEÜ (Electronics and Communication), 17-26 (1996). The present invention can be applied to such general RSSE methods. In addition, the present invention can be applied to another subclass of RSSE, referred to as Parallel Decision-Feedback Equalization, described in Lee and Messerschmidt, “Digital Communication,” 2nd ed. (1994). These publications are each incorporated by reference herein.
Now, RSSE will be explained for the case that L=4 and K=1. Then, a state in the reduced-state trellis is defined according to equation (8) as:

β′n=(bn−1)  (13)

and the number of states in the reduced-state trellis is equal to 21=2. FIG. 4 illustrates the reduced-state trellis 400 corresponding to the full state trellis 300 of FIG. 3 that describes an ISI channel having a memory L=4. A state at time n in the reduced-state trellis is denoted by σn, i.e., σn=β′n. There are two channel states, and two branches corresponding to the information symbols bn=0 and bn=1 leave each state σn to reach respective successor states {σn+1}.

FIG. 5 is a schematic block diagram for an exemplary conventional reduced-state Viterbi detector 500 with local feedback. As shown in FIG. 5, the reduced-state detector 500 includes a decision-feedback unit (DFU) that computes separate ISI estimates for each trellis state according to equation (10) using local feedback, a branch metric unit (BMW) that computes branch metrics for all transitions, an add-compare-select unit (ACSU) that determines the best survivor path into each state, and a survivor memory unit (SMU) that stores the survivor paths.
As shown in FIG. 5, due to the local feedback the critical path 510 is comprised of a recursive loop that includes each of the processing blocks (i.e., the BMU, ACSU, SMU and DFU). As all operations along this critical path 510 have to be performed within one clock period, this recursive loop limits the maximum achievable data rate. Therefore, the maximum data rate of a reduced-state Viterbi detector with local feedback is significantly lower than the maximum data rate of a Viterbi detector without local feedback, which is only limited by the ACS function.
FIG. 6 illustrates a detailed state-parallel reduced-state Viterbi detector implementation 600 with local feedback corresponding to the trellis 400 of FIG. 4, having a memory L=4 and a shortened channel memory K=1. ân−4(0n) is the survivor symbol for time step n−4 from the survivor path into state 0n sn+1(0n+1) is the ACS decision for the two path extensions into state 0n+1. The part of the SMU that stores the L−K survivor symbols ân−K−1(σn), ân−K−2(σn), . . . , ân−L(σn) for each reduced state is implemented with a register-exchange-architecture, as these decisions are required for the computation of ISI estimates in the DFU without delay. The implementation of the SMU using a register-exchange architecture is described, e.g., in R. Cypher and C. B. Shung, “Generalized Trace-Back Techniques for Survivor Memory Management in the Viterbi Algorithm,” Journal of VLSI Signal Processing, 85-94 (1993). Because the discussed exemplary channel uses two signal levels, the multipliers in the DFU can be implemented with a shift operation. The squaring operation for the Euclidean distance computation in the BMU can be approximated using random logic or a look-up table.
Reduced-state Viterbi detection with local feedback that implements, e.g., RSSE, is associated with less computational complexity than full-state Viterbi detection that implements MLSE for the same channel memory L, as it processes less states. However, this comes at the expense of a significantly longer critical path, which is drawn in FIG. 6 using dotted lines. The critical path comprises one symbol multiplication and L−K additions in the DFU (the first term in the right hand side of equation (10) can be computed outside the loop), one addition, subtraction and squaring operation in the BMU, one add-compare in the ACSU, and a 2-to-1 MUX in the SMU. All the operations along this critical path must be completed within one symbol period and cannot be pipelined. In contrast to this, the critical path in a Viterbi detector just comprises the ACS operation. Therefore, the maximum data rate of a reduced-state Viterbi detector implementation with local feedback is potentially significantly lower compared to a Viterbi detector that per forms MLSE. Furthermore, the maximum throughput of a reduced-state Viterbi detector implementation with local feedback depends on the channel memory such that it decreases for increasing L.
Precomputation Architecture
The present invention employs two techniques to increase the maximum data rate that may be achieved by the reduced-state sequence estimator 500. First, as discussed below in conjunction with FIGS. 7-11, branch metrics, ISI-free signal estimates or ISI estimates are precomputed, and the correct values are selected based on survivor symbols or ACS decisions. In this manner, the computations of ISI estimates, ISI-free signal estimates or branch metrics are removed from the critical path. Second, branch metrics, ISI-free signal estimates or ISI estimates can be selected in a pipelined fashion using a multiplexer network structure that corresponds to the structure of the trellis considered by the detector; as shown in FIGS. 12-16.
FIG. 7 is a schematic block diagram of a reduced-state Viterbi detector 700 with local feedback incorporating features of the present invention. As shown in FIG. 7, the branch metric computation is removed from the critical path by precomputing branch metrics at stage 710. For a detailed discussion of the precomputation of branch metrics for a trellis, see, for example, U.S. patent application Ser. No. 09/471,920, entitled, “Method and Apparatus for Shortening the Critical Path of Reduced Complexity Sequence Estimation Techniques,” incorporated by reference herein.
The correct branch metrics are selected at stage 720, discussed below in conjunction with FIG. 8, based on survivor symbols. As previously indicated, the data rate can be increased by precomputing branch metrics and selecting the appropriate ones based on past survivor symbols. As shown in FIG. 7, a pipeline stage 715 can be inserted between the branch metric precomputation 710 and branch metric selection 720. In the implementation shown in FIG. 7, the critical path comprises the branch metric selection 720, ACSU 730 and SMU 740. However, the computation of ISI estimates and branch metrics is not part of the critical path in contrast to the conventional reduced-state Viterbi detection implementation shown in FIG. 5.
In the exemplary uncoded channel model described above, the input into the reduced-state detector is given by equation (2), and a state in the reduced-state trellis is defined by equation (8). A state in the reduced-state trellis is denoted by σn, i.e., σn=β′n. A branch metric for the transition from state σn to σn+1 that corresponds to the information bit sequence bn, bn−1, . . . bn−L is given by




  
    
      
        
          
            
              λ
              n
            
            ⁡
            
              (
              
                
                  b
                  n
                
                ⁢
                
                  b
                  
                    n
                    -
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                ...
                ⁢
                
                  
                
                ⁢
                
                  b
                  
                    n
                    -
                    L
                  
                
              
              )
            
          
          =
          
            
              (
              
                
                  r
                  n
                
                -
                
                  
                    ∑
                    
                      i
                      =
                      0
                    
                    L
                  
                  ⁢
                  
                    
                      f
                      i
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      a
                      
                        n
                        -
                        i
                      
                    
                  
                
              
              )
            
            2
          
        
        ,
      
    
    
      
        (
        14
        )
      
    
  




where an−i is the channel symbol that corresponds to the information bit bn−i. This is the same branch metric that was referred to as λn(σn, an) in the context of equation (11). To account for all possible bit sequences, 2L+1 branch metrics have to be precomputed. For a transition from state σn to σn+1, there are 2L−K branch metric candidates that correspond to the same bit sequence bn−1, bn−2, . . . , bn−K, which is determined by σn, but different speculative bit sequences bn−K−1, bn−K−2, . . . , bn−L. For each state and transition in the reduced-state trellis, the appropriate branch metric is selected based on the L−K survivor symbols ân−K−1, ân−K−2, . . . , ân−L that correspond to this state. The selection of the branch metric associated with the transition from state σn=0n to σn+1=0n+1 is shown in FIG. 8, where L=4 and K=1.

FIG. 8 is a schematic block diagram showing the selection of a branch metric 720 using survivor symbols, the ACS operation 730 and the survivor memory operation 740, as performed by the reduced-state Viterbi detector of FIG. 7. Only the selection of one branch metric, and the ACS operation and survivor memory operation for one state are shown. The detector of FIG. 7 would implement the circuits of FIG. 8 for all required branch metrics and states. In FIG. 8, λn(00xxx) is the selected branch metric for a transition from state 0n that corresponds to the bit sequence bn=0, bn−1=0, bn−2={circumflex over (b)}n−2(0n), bn−3={circumflex over (b)}n−3(0n) and bn−4={circumflex over (b)}n−4(0n). A 2L−K-to-1 multiplexer 810 is required to select the correct branch metric among the precomputed ones. The critical path just comprises the multiplexer 810 for the branch metric selection and an add-compare-select 820. This is significantly shorter compared to a conventional RSSE implementation, as the computation of branch metrics is outside the critical path. Except for the multiplexer 810, the critical path in this RSSE architecture with precomputed branch metrics is the same as in a Viterbi detector that implements MLSE without any decision-feedback.
To reduce the hardware for the precomputation, the branch metrics do not have to be fully precomputed. Instead, ISI-free signal estimates can be precomputed, and branch metrics are then calculated using the correct ISI-free signal estimates that are selected based on past survivor symbols as shown in FIG. 9. An ISI-free signal estimate that corresponds to the bit sequence bn−1, bn−2, . . . , bn−L is given by




  
    
      
        
          
            q
            n
          
          ⁡
          
            (
            
              
                b
                
                  n
                  -
                  1
                
              
              ⁢
              
                b
                
                  n
                  -
                  2
                
              
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              
                b
                
                  n
                  -
                  L
                
              
            
            )
          
        
        =
        
          
            r
            n
          
          -
          
            
              ∑
              
                i
                =
                1
              
              L
            
            ⁢
            
              
                f
                i
              
              ·
              
                
                  a
                  
                    n
                    -
                    i
                  
                
                .
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




In total, 2L ISI-free signal estimates have to be precomputed, but only 2K+1 branch metrics are calculated in this architecture. However, the Euclidean distance metric computation,

λn=(qn−f0an)2  (16)

is in the critical path, while the computation of ISI-free signal estimates is outside the critical path.

FIG. 9 is a schematic block diagram of a reduced-state Viterbi detector 900 incorporating pipelining of the ISI-free signal estimate computation in accordance with the present invention. As shown in FIG. 9, the ISI-free signal estimate computation is removed from the critical path by precomputing ISI-free signal estimates at stage 910.
The correct ISI-free signal estimates are selected at stage 920, discussed below in conjunction with FIG. 10, based on survivor symbols. As previously indicated, the data rate can be increased by precomputing ISI-free signal estimates and selecting the appropriate ones based on past survivor symbols. As shown in FIG. 9, a pipeline stage 915 can be inserted between the ISI-free signal estimate precomputation 910 and ISI-free signal estimate selection 920. In the implementation shown in FIG. 9, the critical path comprises the ISI-free signal estimate selection 920, branch metrics computation 925, ACSU 930 and SMU 940. However, the computation of ISI-free signal estimates is not part of the critical path in contrast to the conventional reduced-state Viterbi detection implementation shown in FIG. 5.
FIG. 10 is a schematic block diagram showing the selection of an ISI-free signal estimate 920 using survivor symbols, the branch metric computation 925, the ACS operation 930 and the survivor memory operation 940, as performed by the reduced-state Viterbi detector of FIG. 9. Only the selection of one ISI-free signal estimate, the computation of one branch metric, and the ACS operation and survivor memory operation for one state are shown. The detector of FIG. 9 would implement the circuits of FIG. 10 for all required ISI-free signal estimates, branch metrics and states. In FIG. 10, qn(0xxx) is the selected ISI-free signal estimate for a transition from state 0n that corresponds to the bit sequence bn−1=0, bn−2={circumflex over (b)}n−2(0n), bn−3={circumflex over (b)}n−3(0n) and bn−4={circumflex over (b)}n−3(0n). A 2L−K-to-1 multiplexer 1010 is required to select the correct ISI-free signal estimate among the precomputed ones. The critical path just comprises the multiplexer 1010 for the ISI-free signal estimate selection, the branch metric computation 1015 and an add-compare-select 1020. This is significantly shorter compared to a conventional RSSE implementation, as the computation of the ISI-free signal estimates is outside the critical path.
To reduce the hardware for the precomputation even further, ISI estimates instead of ISI-free signal estimates can be precomputed, and ISI-free signal estimates and branch metrics are then calculated using the correct ISI estimates that are selected based on past survivor symbols as shown in FIG. 11. The inverse of an ISI estimate that corresponds to the bit sequence bn−1, bn−2, . . . , bn−L is given by




  
    
      
        
          
            q
            n
            ′
          
          ⁡
          
            (
            
              
                b
                
                  n
                  -
                  1
                
              
              ⁢
              
                b
                
                  n
                  -
                  2
                
              
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              
                b
                
                  n
                  -
                  L
                
              
            
            )
          
        
        =
        
          -
          
            
              ∑
              
                i
                =
                1
              
              L
            
            ⁢
            
              
                f
                i
              
              ·
              
                
                  a
                  
                    n
                    -
                    i
                  
                
                .
              
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




In total, 2L ISI estimates have to be precomputed, but only 2K ISI-free signal estimates and 2K+1 branch metrics are calculated in this architecture. An ISI-free signal estimate is then computed according to:

qn=rn+q′n,  (18)

and the corresponding branch metric is given by

λn=(qn−f0·an)2.  (19)

FIG. 11 is a schematic block diagram of an architecture that is a derivation of the architecture shown in FIG. 10. In FIG. 11, ISI estimates are precomputed, whereas in FIG. 10 ISI-free signal estimates are precomputed. In FIG. 11, an ISI-free signal estimate and branch metric are computed based on a selected ISI estimate. Survivor symbols are used to select the correct ISI estimate. While the architecture of FIG. 11 is associated with less hardware complexity than the architecture of FIG. 10, the critical path in FIG. 11 not only includes the branch metric computation 1115, but also the computation of an ISI-free signal estimate 1112. However; the critical path of this architecture is still shorter than the critical path of a conventional reduced-state Viterbi implementation shown in FIG. 6.
It is noted that the improved data rate achieved by the present invention comes at the expense of increased hardware complexity, as shown, for example, in FIGS. 8, 10 and 11, as several branch metric, ISI-free signal estimate or ISI estimate candidates are precomputed pet state transition in the trellis, while only one of these precomputed values is selected for the ACS operation. In contrast, the architecture of FIG. 6, that is associated with a significantly longer critical path, computes only one branch metric, ISI-free signal estimate and ISI estimate pet state transition.
Pipelined Selection
In the precomputation architectures of FIGS. 8 (branch metrics), 10 (ISI-free signal estimates) and 11 (intersymbol interference estimates), a 2L−K-to-1 multiplexer 810, 1010, 1110 lies in the critical path. Although the computation of the branch metrics, ISI-free signal estimates or ISI estimates is not part of the critical path anymore, the delay due to the multiplexer 810, 1010, 1110 still depends on the channel memory. In a straightforward tree-wise implementation of the 2L−K-to-1 multiplexer 810, 1010, 1110 using 2L−K-1 2-to-1 multiplexers, the delay is equal to the delay of L−K 2-to-1 multiplexers, potentially mitigating the speed-up achieved by precomputing branch metrics, ISI-free signal estimates or ISI estimates.
The present invention recognizes that when branch metrics, ISI-free signal estimates or ISI estimates are precomputed L−K time steps in advance, they can be selected using L−K levels of 2-to-1 multiplexers that are driven by ACS decisions and where each level is associated with a pipeline stage. However, only a single 2-to-1 multiplexer is part of the critical path, and the delay associated with the selection of collect values becomes independent of the channel memory.
FIG. 12 is a schematic block diagram of a reduced-state Viterbi detector 1200 incorporating pipelining of the branch metric selection. As shown in FIG. 12, the branch metric computation is removed from the critical path by precomputing branch metrics at stage 1210. The correct blanch metrics are selected at stage 1220, discussed below in conjunction with FIG. 13, based on ACS decisions. As shown in FIG. 12, a pipeline stage 1215 can be inserted between the branch metric precomputation 1210 and branch metric selection 1220.
FIG. 13 is a schematic block diagram showing the pipelined selection of a blanch metric 1220 by the reduced-state Viterbi detector of FIG. 12 using ACS decisions. The exemplary embodiment shown in FIG. 13 shows the pipelined selection of branch metrics with three pipeline stages for L=4 and K=1. The critical path in this architecture just includes a 2-to-1 multiplexer, such as a multiplexer in stage 1310, 1320 or 1330, and an add-compare 1340 in the ACSU. The critical path is significantly shorter compared to FIG. 8, where an 8-to-1 multiplexer 810 lies in the critical path. In fact, it has the same length as in a Viterbi detector that implements MLSE without any decision-feedback.
In FIG. 13, branch metrics are precomputed L−K=3 time units in advance and then selected over three clock periods based on corresponding ACS decisions. As the ACS decision sn+1(σn+1) determines the survivor symbol ân−K(σn+1), the appropriate branch metric can be selected among candidates that correspond to the same symbol sequence an+L−K, an+L−K−1, . . . , an−K+1, but a different past symbol an−K. In FIG. 13, sn+1(0n+1) determines the surviving branch metric λn+3 (0000x) in the top multiplexer in stage 1310 among the candidates λn+3 (00000) and λn+3 (00001), where x is a dummy variable for the symbol an−1. Similarly, sn+1(0n+1) determines the surviving branch metric λn+2 (000xx) in the top multiplexer in stage 1320 among the candidates λn+2 (0000x) and λn+2 (0001x), where xx is a placeholder for the symbols an−1 and an−2. Finally, after selection by the multiplexer in stage 1330, just one surviving branch metric per transition remains, which is used in the ACSU 1340. The branch metric selection resembles the selection of survivor symbols in a SMU that is implemented according to the register exchange architecture. The total number of 2-to-1 multiplexers is 2L−K-1 and thus equal to the number of 2-to-1 multiplexers required for a straightforward tree-wise implementation of the 8-to-1 multiplexer in FIG. 8. Therefore, except for pipeline registers, the branch metric selection in FIG. 13 is associated with about the same complexity as the branch metric selection in FIG. 8.
Generally, branch metrics are precomputed two or more time units in advance and then selected over two or more clock periods. Each pipeline stage has one or more functional units, each comprised of a multiplexer and an associated pipeline register. An input to each functional unit comprises a plurality of precomputed branch metrics that form a subgroup and an output of each subgroup comprises a selected precomputed branch metric. As shown in FIG. 13, at each stage, such as stage 1310, the selected precomputed branch metrics are combined into one or more new subgroups for a following stage, such as stage 1320. In the exemplary embodiment shown in FIG. 13, the subgroups are combined based on the bit pattern of the selected precomputed branch metrics. For example, the selected precomputed branch metrics ln+2(0000x) and ln+2(0001x) have similar bit patterns, other than the second-to-least significant bit and are combined for selection in the second stage 1320. Both branch metrics ln+2(0000x) and ln+2(0001x) correspond to a bit pattern that agrees in the first three bit positions.
FIG. 14 is a schematic block diagram of a reduced-state Viterbi detector 1400 that incorporates pipelining of the ISI-free signal estimate selection. As shown in FIG. 14, the ISI-free signal estimate computation is removed from the critical path by precomputing ISI-free signal estimates at stage 1410. The correct ISI-free signal estimates are selected at stage 1420, discussed below in conjunction with FIG. 15, based on ACS decisions. As shown in FIG. 14, a pipeline stage 1415 can be inserted between the ISI-free signal estimate precomputation 1410 and ISI-free signal estimate selection 1420.
FIG. 15 is a schematic block diagram showing the pipelined selection of an ISI-free signal estimate 1420 by the reduced-state Viterbi detector of FIG. 14 using ACS decisions. The selection of ISI-free signal estimates resembles the selection of branch metrics, discussed above in conjunction with FIG. 13. The exemplary embodiment shown in FIG. 15 shows the pipelined selection of an ISI-free signal estimate with three pipeline stages for L=4 and K=1. The critical path in this architecture just includes a 2-to-1 multiplexer, such as a multiplexer in stage 1510, 1520 or 1530, an add-compare 1540 in the ACSU and a branch metric computation 1535. The critical path is significantly shorter compared to FIG. 10, where an 8-to-1 multiplexer 1010 lies in the critical path.
In FIG. 15, ISI-free signal estimates are precomputed L−K=3 time units in advance and then selected over three clock periods based on corresponding ACS decisions. As the ACS decision sn+1(σn+1) determines the survivor symbol ân−K (σn+1), the appropriate ISI-free signal estimate can be selected among candidates that correspond to the same symbol sequence an+L−K, an+L−K−1, . . . , an−K+1, but a different past symbol an−K. In FIG. 15, sn+1 (0n+1) determines the surviving intersymbol interference estimate qn+3 (000x) in the top multiplexer in stage 1510 among the candidates qn+3 (0000) and qn+3 (0001), where x is a dummy variable for the symbol an−1. Similarly, sn+1 (0n+1) determines the surviving branch metric qn+2 (00xx) in the top multiplexer in stage 1520 among the candidates qn+2 (000x) and qn+2 (001x), where xx is a placeholder for the symbols an−1 and an−2. Finally, after selection by the multiplexer in stage 1530, just one surviving ISI-free signal estimate per transition remains, which is used for the branch metric computation 1535.
Generally, intersymbol interference-free signal estimates are precomputed two or more time units in advance and then selected over two or more clock periods. Each pipeline stage has one or more functional units, each comprised of a multiplexer and an associated pipeline register. An input to each functional unit comprises a plurality of precomputed intersymbol interference-free signal estimates that form a subgroup and an output of each subgroup comprises a selected precomputed intersymbol interference-free signal estimate. As shown in FIG. 15, at each pipeline stage, such as stage 1510, the selected precomputed intersymbol interference-free signal estimates are combined into one or more new subgroups for a following pipeline stage, such as stage 1520. In the exemplary embodiment shown in FIG. 15, the subgroups are combined based on the bit pattern of the selected precomputed intersymbol interference-free signal estimates. For example, the selected precomputed intersymbol interference-free signal estimates qn+2(000x) and qn+2(001x) have similar bit patterns, other than the second-to-least significant bit and are combined for selection in the second stage 1520. Both intersymbol interference-free signal estimates qn+2(000x) and qn+2(001x) correspond to a bit pattern that agrees in the first two bit positions.
FIG. 16 is a schematic block diagram showing the pipelined selection of an intersymbol interference estimate using ACS decisions according to the invention. In FIG. 16, ISI estimates instead of ISI-free signal estimates are selected in a pipelined fashion. The ISI estimates are precomputed according to equation (17). The intersymbol interference estimate selection resembles the selection of the ISI-free signal estimate, discussed above in conjunction with FIG. 15. The exemplary embodiment shown in FIG. 16 shows the pipelined selection of intersymbol interference estimates with three pipeline stages for L=4 and K=1. The critical path in this architecture just includes a 2-to-1 multiplexer; such as a multiplexer in stage 1610, 1620 or 1630, an ISI-free signal estimate computation 1632, a branch metric computation 1635, and an add-compare 1640 in the ACSU. The critical path is significantly shorter compared to FIG. 11, where an 8-to-1 multiplexer 1110 lies in the critical path.
In FIG. 16, intersymbol interference estimates are precomputed L−K=3 time units in advance and then selected over three clock periods based on corresponding ACS decisions. As the ACS decision sn+1 (σn+1) determines the survivor symbol ân−K (σn+1), the appropriate intersymbol interference estimate can be selected among candidates that correspond to the same symbol sequence an+L−K, an+L−K−1, . . . , an−K+1, but a different past symbol an−K. In FIG. 16, sn+1 (0n+1) determines the surviving intersymbol interference estimate q′n+3 (000x) in the top multiplexer in stage 1610 among the candidates q′n+3 (0000) and q′n+3 (0001), where x is a dummy variable for the symbol an−1. Similarly, sn+1 (0n+1) determines the surviving branch metric q′n+2 (00xx) in the top multiplexer in stage 1620 among the candidates q′n+2 (000x) and q′n+2 (001x), where xx is a placeholder for the symbols an−1 and an−2. Finally, after selection by the multiplexer in stage 1630, just one surviving intersymbol interference estimate per transition remains, which is used for the ISI-free signal estimate computation 1632.
FIG. 17 is a functional block diagram 1700 of the ACS operation performed in FIGS. 8, 10, 11, 13, 15, and 16. As shown in FIG. 17, the ACS block 1700 includes an add function 1710, compare function 1720 and select function 1730. The exemplary add function 1710 includes two adders. The exemplary compare function 1720 could be implemented using a subtractor, where the sign bit of the subtractor output controls the selector 1730. The select function 1730 comprises a multiplexer, controlled by the output of the compare function 1720. FIG. 17 shows an exemplary two-way ACS implementation for a trellis with two transitions per state A 4-way ACS structure for a trellis with four transitions per state is shown in United States Patent Application entitled “Method and Apparatus for Multiple Step Viterbi Detection with Local Feedback,” filed simultaneously herewith, assigned to the assignee of the present invention and incorporated by reference herein.
Among other benefits, the present invention allows for a VLSI implementation of reduced-state Viterbi detectors with local feedback for data rates that are significantly increased relative to conventional designs. Even larger data rate increases can be achieved when two or more trellis steps are processed within once clock period using a multi-step trellis. To achieve this additional speed advantage, the invention disclosed here can be combined with the multi-step detection method disclosed in United States Patent Application entitled “Method and Apparatus for Multiple Step Viterbi Detection with Local Feedback,” filed contemporaneously herewith and incorporated by reference herein. The invention uses an architecture that is very regular making it suitable for high-speed implementation. Viterbi detectors with local feedback can achieve better error rate performance than postprocessor-based structures in the magnetic recording application. Therefore, reduced-state Viterbi detection with local feedback is an attractive detector structure for future read channel chips. The use of reduced-state Viterbi detection with local feedback in the magnetic recording application is described in E. F. Haratsch, “Viterbi Detector Architectures for Magnetic Recording,” 2003 International Symposium on VLSI Technology, Systems, and Applications, 243-46, Oct. 6-8, 2003. Post-processor based detector structures are discussed in Z. A. Keirn et al., “On the Use of Redundant Bits for Magnetic Recording: Single Parity Codes and Reed-Solomon ECC,” IEEE Transactions on Magnetics, 225-30 (January 2004), and the references therein.
Trace-Back Survivor Memory
Another benefit of the invention is that ACS decisions can be used to select precomputed branch metrics, ISI-free signal estimates or ISI estimates as shown in FIGS. 12-16. When ACS decisions are used to select precomputed values, the SMU 1240 and 1440 can be implemented using a trace-back structure, as survivor symbols are not used for local feedback anymore. The details of a trace-back survivor memory architecture can be read, e.g. in R. Cypher and C. B. Shung, “Generalized Trace-Back Techniques for Survivor Memory Management in the Viterbi Algorithm,” Journal of VLSI Signal Processing, 85-94 (1993), or in H.-L. Lou, “Implementing the Viterbi algorithm”, IEEE Signal Processing Magazine, 42-52 (September 1995), or in O. J. Joeressen and H. Meyr, “Viterbi Decoding with Dual Timescale Traceback Processing,” IEEE International Symposium on Personal, Indoor and Mobile Radio Communications, 213-217 (September 1995), each incorporated by reference herein.
In a register-exchange survivor memory implementation, survivor symbols for each state are stored and updated at each detection step. In a trace-back implementation, however, ACS decisions are stored as pointers in a memory, and the detected symbols ale obtained by tracing back the pointers that correspond to a survivor path. As the trace-back architecture does not require the updating of all survivor symbols at each detection step, it is associated with less power consumption than the register-exchange architecture. However, the trace-back architecture is associated with larger detection latency and therefore not suitable for the reduced-state Viterbi detector shown in FIG. 6, where zero delay survivor symbols are required for the local feedback to compute ISI estimates and branch metrics. However, the disclosed architectures shown in FIGS. 12-16 use ACS decision to select precomputed branch metrics, ISI-free signal estimates or ISI estimates, therefore making it possible to implement the survivor memory SMU using a trace-back architecture. In this case, the trace-back SMU will be associated with significantly less power consumption than a corresponding register-exchange SMU implementation.
Magnetic Recording Read Channels
The techniques described herein can be employed, e.g., to detect data in the presence of intersymbol interference and noise in magnetic recording read channels. The disclosed reduced-state Viterbi detectors with local feedback improve the detection of read data bits compared to post-processor based structures. In particular, the invention can be used to implement a read channel that performs noise-predictive data detection and achieves the ever increasing high data rates that are required by evolving storage applications. For a discussion of noise-predictive detection in magnetic recording, see, e.g., R. D. Cideciyan et al., “Noise Predictive Maximum Likelihood Detection Combined With Parity-Based Post-Processing,” IEEE Trans on Magnetics, 714-20 (March 2001), and E. F. Haratsch, “Viterbi Detector Architectures for Magnetic Recording,” International Symposium on VLSI Technology, Systems, and Applications, 243-46 (October 2003).
The simplified block diagram for a read channel incorporating noise-predictive reduced-state Viterbi detection is shown in FIG. 18, where signals received at the input of the finite response (FIR) equalizer are in fact signals that have been processed by the analog front-end, which typically includes a variable gain amplifier; continuous time filter and A/D converter. The FIR equalizer 1810 shapes the channel impulse response such that the signals at the output of the FIR equalizer yn can be described by the equation:




  
    
      
        
          
            y
            n
          
          =
          
            
              
                ∑
                
                  i
                  =
                  0
                
                M
              
              ⁢
              
                
                  h
                  i
                
                ·
                
                  a
                  
                    n
                    -
                    i
                  
                
              
            
            +
            
              v
              n
            
          
        
        ,
      
    
    
      
        (
        20
        )
      
    
  




where an are the data symbols defined as in equation (1), hi are the equalization target coefficients, M is the equalization target order, and vn is the noise at the output of the FIR equalizer. The equalization target is chosen such that its frequency spectrum matches the characteristics of the read channel well. The impulse response associated with the equalization target can be described by the equation:

H(D)=h0+h1·D+h2·D2+ . . . +hM·DM.  (21)

The error rate performance of a lead channel can be improved by employing a noise-predictive FIR (NP-FIR) filter 1820 after the FIR equalizer 1810 that whitens the noise. The impulse response associated with the NP-FIR can be characterized with the polynomial:

P(D)=p0+p1·D+p2·D2+ . . . +pN·DN,  (22)

where pi, 0≦i≦N are the coefficients and N is the order of the NP-FIR filter.

The subsequent reduced-state Viterbi detector considers a channel response with the polynomial:

F(D)=f0+f1·D+f2·D2+ . . . +fM+N·DM+N=H(D)·P(D),  (23)

and the signals at the input of the reduced-state Viterbi detector are given by:





  
    
      
        
          
            r
            n
          
          =
          
            
              
                ∑
                
                  i
                  =
                  0
                
                L
              
              ⁢
              
                
                  f
                  i
                
                ·
                
                  a
                  
                    n
                    -
                    i
                  
                
              
            
            +
            
              w
              n
            
          
        
        ,
      
    
    
      
        (
        24
        )
      
    
  




where fi, 0≦i≦L, are the channel coefficients of the channel seen by the reduced-state Viterbi detector, L is the channel memory, and wn is the remaining noise at the output of the NP-FIR filter 1820. The channel memory L would be typically equal to L=M+N, but the reduced-state detector could also consider a channel with shorter channel memory, i.e. L<M+N. The channel coefficients fi are given by the convolution of the equalization target and the impulse response of the NP-FIR filter 1820 (see equation (23)). Equation (24) is equivalent to equation (2). Therefore, the features of the invention can be applied to the read channel application, i.e., a reduced-state can be defined as in equation (8), branch metrics can be precomputed as in equation (14) using the channel coefficients defined by (23), and an ACS operation can be performed as in equation (12). A correct branch metric for a transition in the trellis can be selected using the architectures shown in FIG. 8 or FIG. 13. Instead of branch metrics, ISI-free signal estimates or ISI estimates can be precomputed according to equations (15) and (17) respectively using the channel coefficients defined by (23), and a collect ISI-free signal estimate or ISI estimate can be selected as shown in FIGS. 10, 11, 15 and 16 respectively. Therefore, the reduced-state Viterbi detector 1830 can be implemented using the architectures shown in FIGS. 7-16.

The invention can also be applied when a parity check code is used to achieve coding gain. For example, when a one-bit parity check code is used, a state σn in the full-state trellis is given by:

σn=(sn−1;bn−1,bn−2, . . . , bn−L),  (25)

where sn is the running parity syndrome given by the XOR-sum:

sn=bn⊕sn−1.  (26)

The total number of states in the reduced-state trellis that accounts for the parity check code is given by:

2×2L  (27)

Analogous to equation (8), a state σn in the reduced-state trellis can be defined by considering only the past K information bits or symbols:

σn=(sn−1;bn−1,bn−2, . . . , bn−K),  (28)

and the number of states in the reduced-state trellis is equal to:

2×2K.  (29)

A conventional implementation of a reduced-state Viterbi detector that considers the reduced-state trellis defined according to equation (28) would use equations (10)-(12) to compute ISI estimates, blanch and path metrics. However, due to the local feedback required for the computation of ISI estimates, it is challenging to achieve very high data rates. However, the maximum achievable data rate can again be increased significantly by precomputing branch metrics and selecting correct ones as described above. Alternatively, ISI-free signal estimates or ISI estimates can be precomputed and correct ones can be selected as described above.
Branch metrics can be precomputed for possible bit sequences according to equation (14). As described above for L=4 and K=1, the required branch metric λn (bnbn−1xxx) for a transition from state σn=(sn−1;bn−1) is selected among the precomputed branch metrics λn (bnbn−1000), λn(bnbn−1001), λn(bnbn−1010), λn(bnbn−1011), λn(bnbn−1100), λn(bnbn−1101), λn(bnbn−1110), and λn(bnbn−1111) based on survivor symbols or ACS decisions.
The invention can also be applied to signal-dependent detection, which is sometimes referred to as data-dependent detection and explained in detail in the co-pending United States Patent Application entitled “Method and Apparatus for Generating Filter Tap Weights and Biases for Signal Dependent Branch Metric Computation,” incorporated by reference herein. In signal-dependent detection, more than one signal-dependent (SD) NP-FIR filters operate in parallel to whiten the noise. FIG. 19 illustrates this for the case that two SD NP-FIR filters 1920-1 and 1920-2 are used. The invention can easily be used when there are more than two NP-FIR filters 1920. The reduced-state Viterbi detector 1930 can be implemented using the architectures shown in FIGS. 7-16.
In FIG. 19, the output of the FIR equalizer 1910 is supplied to two SD NP-FIR filters 1920-1 and 1920-2 to produced two signals rn(1) and rn(2). Each SD NP-FIR filter 1920 implements the impulse response defined by equation (22) with a different set of filter coefficients. For example, the first SD NP-FIR filter 1920-1 that produces rn(1) uses a first set of coefficients pi(1), 0≦i≦N, whereas the second SD NP-FIR filter 1920-2 that produces rn(2) uses a second set of coefficients pi(2), 0≦i≦N that can differ from the first set of NP-FIR filter coefficients. The corresponding polynomials that describe the SD-NP FIR filters 1920 are denoted P(D;1) and P(D;2), e.g.,

P(D;1)=p0(1)+p1(1)·D+pN(1)·DN.  (30)

The filter coefficients of the different SD NP-FIR filters 1920 can differ; as in a signal-dependent channel the noise statistics depend on the transmitted data or bit sequence. The generation of filter coefficients for the SD NP-FIR filters is described in co-pending United States Patent Application entitled “Method and Apparatus for Generating Filter Tap Weights and Biases for Signal Dependent Branch Metric Computation,” incorporated by reference herein.

For the considered channel with two SD NP-FIR filters 1920, the reduced-state Viterbi detector 1930 would compute branch metrics considering two different channel impulse responses with the polynomials F(D;1) and F(D;2) that are given by:

F(D;1)=f0(1)+f1(1)·D+ . . . fM+N(1)·DM+N=H(D)·P(D;1), and  (31)

F(D;2)=f0(2)+f1(2)·D+ . . . fM+N(2)·DM+N=H(D)·P(D;2).  (32)

In a signal-dependent channel, the filter coefficients fi that are used to compute a branch metric depend on the transmitted data or bit sequence. For two SD NP-FIR filters 1920, branch metrics are precomputed according to equation (14) for a first group of bit sequences (bn−1bn−2 . . . bn−L) using filter coefficients fi(1), and for a second group of bit sequences (bn−1bn−2 . . . bn−L) the filter coefficients fi(2) are used.
For example, signal-dependent branch metrics for transitions from states σn to σn+1 that correspond to all bit sequences starting with (bnbn−1)=(00) or (bnbn−1)=(11) are computed using channel coefficients fi(1) and the sample rn(1):




  
    
      
        
          
            
              λ
              n
            
            ⁡
            
              (
              
                00
                ⁢
                
                  
                
                ⁢
                
                  b
                  
                    n
                    -
                    2
                  
                
                ⁢
                
                  
                    b
                    
                      n
                      -
                      3
                    
                  
                  ·
                  
                    b
                    
                      n
                      -
                      L
                    
                  
                
              
              )
            
          
          =
          
            
              (
              
                
                  
                    r
                    n
                  
                  ⁡
                  
                    (
                    1
                    )
                  
                
                -
                
                  
                    ∑
                    
                      i
                      =
                      0
                    
                    L
                  
                  ⁢
                  
                    
                      
                        f
                        i
                      
                      ⁡
                      
                        (
                        1
                        )
                      
                    
                    ·
                    
                      a
                      
                        n
                        -
                        i
                      
                    
                  
                
              
              )
            
            2
          
        
        ,
      
    
    
      
        (
        33
        )
      
    
  
  
    
      
        
          
            λ
            n
          
          ⁡
          
            (
            
              11
              ⁢
              
                
              
              ⁢
              
                b
                
                  n
                  -
                  2
                
              
              ⁢
              
                
                  b
                  
                    n
                    -
                    3
                  
                
                ·
                
                  b
                  
                    n
                    -
                    L
                  
                
              
            
            )
          
        
        =
        
          
            
              (
              
                
                  
                    r
                    n
                  
                  ⁡
                  
                    (
                    1
                    )
                  
                
                -
                
                  
                    ∑
                    
                      i
                      =
                      0
                    
                    L
                  
                  ⁢
                  
                    
                      
                        f
                        i
                      
                      ⁡
                      
                        (
                        1
                        )
                      
                    
                    ·
                    
                      a
                      
                        n
                        -
                        i
                      
                    
                  
                
              
              )
            
            2
          
          .
        
      
    
    
      
        (
        34
        )
      
    
  




Continuing this example, signal-dependent branch metrics for all bit sequences that start with (bnbn−1)=(01) or (bnbn−1)=(10) are computed using the second of channel coefficients fi(2) and the second sample rn(2):




  
    
      
        
          
            
              λ
              n
            
            ⁡
            
              (
              
                01
                ⁢
                
                  
                
                ⁢
                
                  b
                  
                    n
                    -
                    2
                  
                
                ⁢
                
                  b
                  
                    n
                    -
                    3
                  
                
                ⁢
                
                  
                
                ⁢
                ...
                ⁢
                
                  
                
                ⁢
                
                  b
                  
                    n
                    -
                    L
                  
                
              
              )
            
          
          =
          
            
              (
              
                
                  
                    r
                    n
                  
                  ⁡
                  
                    (
                    2
                    )
                  
                
                -
                
                  
                    ∑
                    
                      i
                      =
                      0
                    
                    L
                  
                  ⁢
                  
                    
                      
                        f
                        i
                      
                      ⁡
                      
                        (
                        2
                        )
                      
                    
                    ·
                    
                      a
                      
                        n
                        -
                        i
                      
                    
                  
                
              
              )
            
            2
          
        
        ,
      
    
    
      
        (
        35
        )
      
    
  
  
    
      
        
          
            λ
            n
          
          ⁡
          
            (
            
              10
              ⁢
              
                
              
              ⁢
              
                b
                
                  n
                  -
                  2
                
              
              ⁢
              
                b
                
                  n
                  -
                  3
                
              
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              
                b
                
                  n
                  -
                  L
                
              
            
            )
          
        
        =
        
          
            
              (
              
                
                  
                    r
                    n
                  
                  ⁡
                  
                    (
                    2
                    )
                  
                
                -
                
                  
                    ∑
                    
                      i
                      =
                      0
                    
                    L
                  
                  ⁢
                  
                    
                      
                        f
                        i
                      
                      ⁡
                      
                        (
                        2
                        )
                      
                    
                    ·
                    
                      a
                      
                        n
                        -
                        i
                      
                    
                  
                
              
              )
            
            2
          
          .
        
      
    
    
      
        (
        36
        )
      
    
  




Comparing the signal-dependent branch metric equations (33)-(36) with the non signal-dependent branch metric equation (14), signal-dependent blanch metrics are computed using the same underlying function, but the choice of the sample rn and channel coefficients fi depends on the bit sequence for which the branch metric is computed. The precomputation of signal-dependent branch metrics was illustrated here using two signal-dependent NP-FIR filters and a particular grouping of bit sequences, but it is apparent how signal-dependent blanch metrics are precomputed for more than two signal-dependent NP-FIR filters and other groupings. E.g., all possible bit sequences of length L can be divided into more than two groups, for which separate samples rn and separate sets of channel coefficients fi would be used to precompute branch metrics.
The selection of the correct branch metric, the ACS operation and the SMU are implemented as described above for the non-signal dependent detector. Therefore all the benefits of the invention apply to signal-dependent detection as well. In a further variation of the present invention, signal-dependent ISI-free signal estimates or intersymbol interference estimates can be precomputed, instead of branch metric, as would be apparent to a person of ordinary skill in the art based on the present disclosure. Then, signal-dependent blanch metrics would be computed based on selected signal-dependent ISI-free signal estimates of ISI estimates as described above.
It is to be understood that the embodiments and variations shown and described herein are merely illustrative of the principles of this invention and that various modifications may be implemented by those skilled in the art without departing from the scope and spirit of the invention.

## Claims

1. A method for processing a received signal, said method comprising the steps of:
processing said received signal using an FIR equalizer to generate an equalized received signal;
processing said equalized received signal using a plurality of noise-predictive FIR filters to generate a plurality of signals;
precomputing branch metrics for speculative sequences of one or more channel symbols, wherein at least one of said branch metrics is precomputed using an output of one of said plurality of noise-predictive FIR filters based on a bit pattern;
storing said precomputed branch metrics in at least one pipeline register;
selecting one of said precomputed branch metrics from one of said at least one pipeline register based on one or more of at least one survivor symbol and at least one ACS decision from at least one corresponding state; and
selecting a path having a best path metric for a given state.

2. The method of claim 1, wherein said plurality of signals are noise whitened signals.

3. The method of claim 1, wherein said plurality of noise-predictive FIR filters each use a corresponding set of filter coefficients.

4. The method of claim 1, wherein a set of channel coefficients used to precompute at least one of said branch metrics is dependent on a bit pattern, and wherein said set of channel coefficients is given by a convolution of impulse responses of said FIR equalizer and one of said plurality of noise-predictive FIR filters.

5. The method of claim 1, wherein said processing of said received signal employs a trellis that includes syndrome information for a parity check code.

6. The method of claim 1, wherein said bit pattern is of length L, and all possible bit patterns of length L are divided in two or more subgroups, for which separate signals and separate sets of channel coefficients are used to precompute branch metrics.

7. A signal processor, comprising:
an FIR equalizer to process a received signal to generate an equalized received signal;
a plurality of noise-predictive FIR filters to process said equalized received signal to generate a plurality of signals;
a branch metric unit for precomputing branch metrics for speculative sequences of one or more channel symbols, wherein at least one of said branch metrics is precomputed using an output of one of said plurality of noise-predictive FIR filters based on a bit pattern;
at least one pipeline register for storing said precomputed branch metrics;
at least one multiplexer for selecting one of said precomputed branch metrics from one of said at least one pipeline register based on one or more of at least one survivor symbol and at least one ACS decision from at least one corresponding state; and
an add-compare-select unit for selecting a path having a best path metric for a given state.

8. The signal processor of claim 7, wherein said bit pattern is of length L, and all possible bit patterns of length L are divided in two or more subgroups, for which separate signals and separate sets of channel coefficients are used to precompute branch metrics.

9. A method for processing a received signal, said method comprising the steps of:
precomputing intersymbol interference-free signal estimates for speculative sequences of one or more channel symbols;
storing said precomputed intersymbol interference-free signal estimates in at least one pipeline register;
selecting one of said precomputed intersymbol interference-free signal estimates from one of said at least one pipeline registers based on one or more of at least one survivor symbol and at least one ACS decision from at least one corresponding state; and
selecting a path having a best path metric for a given state.

10. The method of claim 9, further comprising the steps of:
processing said received signal using an FIR equalizer to generate an equalized received signal; and
processing said equalized received signal using a noise-predictive FIR filter to generate a signal.

11. The method of claim 10, wherein said signal is a noise whitened signal.

12. The method of claim 10, wherein channel coefficients used to precompute said intersymbol interference-free signal estimates are a convolution of impulse responses of said FIR equalizer and said noise-predictive FIR filter.

13. The method of claim 9, further comprising the steps of:
processing said received signal using an FIR equalizer to generate an equalized received signal; and
processing said equalized received signal using a plurality of noise-predictive FIR filters to generate a plurality of signals.

14. The method of claim 13, wherein said plurality of noise-predictive FIR filters each use a corresponding set of filter coefficients.

15. The method of claim 13, wherein said step of precomputing intersymbol interference-free signal estimates precomputes at least one intersymbol interference-free signal estimates by using an output of one of said plurality of noise-predictive FIR filters based on a bit pattern.

16. The method of claim 15, wherein said bit pattern is of length L, and alt possible bit patterns of length L are divided in two or more subgroups, for which separate signals and separate sets of channel coefficients are used to precompute intersymbol interference-free signal estimates.

17. The method of claim 13, wherein a set of channel coefficients used to precompute at least one of said intersymbol interference-free signal estimates is dependent on a bit pattern, and wherein said set of channel coefficients is given by a convolution of impulse responses of said FIR equalizer and one of said plurality of noise-predictive FIR filters.

18. The method of claim 9, wherein said processing of said received signal employs a trellis that includes syndrome information for a parity check code.

19. A signal processor, comprising:
an intersymbol interference unit for precomputing intersymbol interference-free signal estimates for speculative sequences of one or more channel symbols;
at least one pipeline register for storing said precomputed intersymbol interference-free signal estimates;
at least one multiplexer for selecting one of said precomputed intersymbol interference-free signal estimates from one of said at least one pipeline register based on one or more of at least one survivor symbol and at least one ACS decision from at least one corresponding state; and
an add-compare-select unit for selecting a path having a best path metric for a given state.

20. The signal processor of claim 19, wherein said received signal is processed using an FIR equalizer to generate an equalized received signal and said equalized received signal is processed using a plurality of noise-predictive FIR filters to generate a plurality of signals, wherein at least one of said intersymbol interference-free signal estimates is precomputed using an output of one of said plurality of noise-predictive FIR filters based on a bit pattern and wherein said bit pattern is of length L, and all possible bit patterns of length L are divided in two or more subgroups, for which separate signals and separate sets of channel coefficients are used to precompute intersymbol interference-free signal estimates.