# Patent

## Title
Adaptive read and write systems and methods for memory cells

## Abstract

Adaptive memory read and write systems and methods are described herein that adapts to changes to threshold voltage distributions of memory cells as of result of, for example, the detrimental affects of repeated cycling operations of the memory cells. The novel systems may include at least multi-level memory cells, which may be multi-level flash memory cells, and a computation block operatively coupled to the multi-level memory cells. The computation block may be configured to compute optimal or near optimal mean and detection threshold values based, at least in part, on estimated mean and standard deviation values of level distributions of the multi-level memory cells. The optimal or near optimal mean and detection threshold values computed by the computation block may be subsequently used to facilitate writing and reading, respectively, of data to and from the multi-level memory cells.

## Background
CROSS-REFERENCES TO RELATED APPLICATIONS
The present application claims priority to U.S. Patent Application No. 60/864,468, filed Nov. 6, 2006, entitled “Adaptive Read and Write Systems and Methods for Flash Memory,” and U.S. Patent Application No. 60/910,325, filed Apr. 5, 2007, entitled “Parameter Estimation for NV Memory,” the entire disclosures of which are hereby incorporated by reference in its entirety for all purposes.


TECHNICAL FIELD
Embodiments of the present invention relate to the field of data memory devices, and more particularly, to storage and retrieval of data to and from memory cells.
BACKGROUND
Memory cells, such as flash memory cells, may store data by trapping granulized amounts of charge in, for example, an isolated region of a transistor. In such devices, data retrieval from a memory cell is typically made possible by applying a read voltage to the transistor and subsequently estimating the readout current which is determined by the amount of charge trapped in the cell.
An example of a basic type of memory cell is one that may store 1-bit of information. In such a memory cell, the memory cell may hold or not hold a charge to indicate, for example, logic 1 when a charge is stored, and to indicate logic 0, when no charge is stored.
In contrast, “multi-level memory cells” may be able to store more than 1-bit of information by taking advantage of the ability of such memory cells to hold varying amounts of charge or charge levels. For example, suppose the maximum number of trapped charge allowed in a multi-level memory cell is Q. It may then be possible to store more than 1 bit of information in such a cell by storing a granulized amount of charge between 0 and Q, and subsequently estimating the amount of charge stored during readout of the cell. Thus, for example, 2 bits of information may be stored in one multi-level memory cell by trapping any one of, for example, four levels of charges: 0, Q/3, 2Q/3, Q. This process of trapping charges may be referred to as programming.
In practice, it is often difficult to precisely program a multi-level memory cell with a desired amount of charges. Indeed, the actual programmed amount of charges approximately follows a Gaussian distribution centered on a desired charge level. The variance of the distribution may be determined by the programming method as well as the physical properties of the memory cell. Consequently, the threshold voltage distributions of flash memory cells are also Gaussian.
FIG. 1 illustrates four threshold voltage distributions (herein “level distributions”) for a 2-bit memory cell. The four level distributions depicted are associated with four different levels of charge that may be stored in a memory cell, each level distribution having its own mean and variance. As depicted in FIG. 1, the intersections of the four charge levels (level 0, level 1, level 2, and level 3) define three detection thresholds (t1, t2, and t3) That is, the three detection thresholds (t1, t2, and t3) are located where curves of two adjacent level distributions intersect.
In order to properly write and read data to and from a multi-level memory cell, two things should be known: the detection thresholds and the means of the level distributions of the multi-level memory cell. In particular, the detection thresholds (e.g., t1, t2, and t3) may be needed in order to read data from the memory cells, and the means (e.g., m1, m2, m3, and m4) of the level distributions may be needed in order to write data to the memory cell. That is, the detection thresholds are needed during a read operation of a multi-level memory cell in order to determine whether the charge stored in the memory cell is at level 0, level 1, level 2, or level 4. In contrast, the means of the level distributions are needed during a write operation of a multi-level memory cell in order to more accurately target the amount of charge to be programmed into the memory cell.
For example, in order to determine whether the total charge stored in a multi-level memory cell is in level 0 during a read operation, the value of the first detection threshold (t1) should be known. By knowing the value of t1, one would simply determine whether the charge stored (or not stored since level 0 could be zero charge) in the memory cell is less than t1 in order to determine whether the stored charge is at level 0. Similarly, in order to determine whether the charge stored in the memory cell is at level 1, you would determine whether the charge stored in the memory cell is between t1 and t2.
In contrast, in order to target the right amount of charge to program into a multi-level memory cell during a write operation, the means (herein “mean values”) of the level distribution should be known. For example, referring back to FIG. 1, if one wanted to store level 2 amount of charge in the memory cell, one would need to know the second mean value (m1) in order to properly program the memory cell. By targeting m1 amount of charge to be stored in the memory cell, error may be minimized since m1 is located at the top of the Gaussian curve.
Unfortunately, memory cells, such as the multi-level flash memory cells described above, may be subject to retention loss after undergoing read and/or write cycling. As a result, the mean and variance of the level distributions change after cycling (e.g., read and write operations) as illustrated FIG. 2. In order to account for the degradation of such memory cells and to minimize error during read and write operations of such memory cells, memory read/write systems need to track not only the changes to the level distributions, but also to adaptively adjust the read and write processes to mitigate the detrimental effects of repeated cycling operations.
SUMMARY OF INVENTION
According to various embodiments of the present invention, adaptive memory read and write systems and methods are provided that may adjust to level distributions changes of memory cells. For the embodiments, the systems may include at least multi-level memory cells, which may be multi-level flash memory cells, and a computation block operatively coupled to the multi-level memory cells. The computation block may be configured to compute optimal or near optimal mean and detection threshold values based, at least in part, on estimated mean and standard deviation values of level distributions of the multi-level memory cells. The optimal or near optimal mean and detection threshold values computed by the computation block may be subsequently used to facilitate writing and reading, respectively, of data to and from the multi-level memory cells.
The computation block included in the systems may be operatively coupled to the multi-level memory cells via an estimation block, the estimation block being configured to determine the estimated mean and standard deviation values used by the computation block to compute the optimal or near optimal mean and detection threshold values. In some embodiments of the present invention, the multi-level memory cells may include one or more pilot cells to store therein predetermined data. For these embodiments, the estimation block may be further configured to determine the estimated mean and standard deviation values using the one or more pilot cells.
The systems may further include a look-up table that is operatively coupled to the computation block. The look-up table may be configured to store the estimated mean and standard deviation values determined by the estimation block and/or the optimal or near optimal mean and detection threshold values computed by the computation block.
The systems may also include a read block and a write block that are each operatively coupled to the look-up table and configured to read and write, respectively, data from and to the multi-level memory cells. The read block may read data from the multi-level memory cells based, at least in part, on optimal or near optimal detection threshold values that may be stored in the look-up table. Alternatively, the optimal or near optimal detection threshold values may be provided directly from the computation block. In contrast, the write block may be configured to write data to the multi-level memory cells based, at least in part, on optimal or near optimal mean values that may be stored in the look-up table. Alternatively, the optimal or near optimal mean values may be provided directly from the computation block.
In some embodiments of the present invention, the multi-level memory cells may include at least one M-level memory cell having four level distributions. For these embodiments, the computation block may be further configured to compute near optimal detection threshold values of the M-level memory cell according to the equation:




  
    t
    i
  
  =
  
    
      
        
          m
          
            i
            -
            1
          
        
        ⁢
        
          σ
          i
        
      
      +
      
        
          m
          
            i
            +
            1
          
        
        ⁢
        
          σ
          
            i
            -
            1
          
        
      
    
    
      
        σ
        i
      
      +
      
        σ
        
          i
          -
          1
        
      
    
  




for 1≦i≦M−1, where ti is the detection threshold between level i−1 and level i; mi, 0≦i≦M−1 are the estimated means of the M-level distributions; σi, 0≦i≦M−1 are the estimated standard deviations of the M-level distributions.
In some embodiments of the present invention, the multi-level memory cells may include at least one M level memory cell having M levels. For these embodiments, the computation block may be configured to compute an ith level near optimal mean value ({tilde over (m)}) of the M level memory cell according to the equation:




  
    
      m
      ~
    
    i
  
  =
  
    
      m
      0
    
    +
    
      
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  i
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            i
          
        
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  M
                  -
                  2
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            
              M
              -
              1
            
          
        
      
      ⁢
      L
    
  




where mi is the estimated mean value of the ith level of the M level memory cell, σi is the estimated standard deviation value of the ith level of the M level memory cell, and L is equal to mM-1−m0. These and other aspects of various embodiments of the present will be described in greater detail in the following description.



BRIEF DESCRIPTION OF THE DRAWINGS
The present invention will be described by way of exemplary embodiments, but not limitations, illustrated in the accompanying drawings in which like references denote similar elements, and in which:
FIG. 1 illustrates four threshold voltage distributions of an exemplary 2-bit memory cell;
FIG. 2 illustrates the four threshold voltage distributions of the exemplary 2-bit memory cell of FIG. 1 after cycling;
FIG. 3 illustrates an adaptive read/write memory system, in accordance with various embodiments of the present invention;
FIG. 4 illustrates an exemplary table structure of a look-up table, in accordance with various embodiments of the present invention;
FIG. 5 illustrates an adaptive signal detection flow process for reading multi-level memory cells using computed optimal or near optimal detection threshold values, in accordance with various embodiments of the present invention;
FIG. 6 illustrates an adaptive flow process for a write operation of multi-level memory cells, in accordance with various embodiments of the present invention;
FIG. 7 illustrates a 3-level threshold voltage distribution (“3-level distribution”) of an exemplary multi-level memory cell, in accordance with various embodiments of the present invention; and
FIG. 8 illustrates a 4-level threshold voltage distribution (“4-level distribution) of an exemplary multi-level memory cell, in accordance with various embodiments of the present invention.



DETAILED DESCRIPTION
In the following detailed description, reference is made to the accompanying drawings which form a part hereof wherein like numerals designate like parts throughout, and in which is shown by way of illustration embodiments in which the invention may be practiced. It is to be understood that other embodiments may be utilized and structural or logical changes may be made without departing from the scope of the present invention. Therefore, the following detailed description is not to be taken in a limiting sense, and the scope of embodiments in accordance with the present invention is defined by the appended claims and their equivalents.
Various operations may be described as multiple discrete operations in turn, in a manner that may be helpful in understanding embodiments of the present invention; however, the order of description should not be construed to imply that these operations are order dependent.
For the purposes of the instant description, the phrase “A/B” means A or B. For the purposes of the instant description, the phrase “A and/or B” means “(A), (B), or (A and B).” For the purposes of the instant description, the phrase “at least one of A, B and C” means “(A), (B), (C), (A and B), (A and C), (B and C) or (A, B and C).” For the purposes of the instant description, the phrase “(A)B” means “(B) or (AB),” that is, A is an optional element.
The description may use the phrases “in various embodiments,” or “in some embodiments,” which may each refer to one or more of the same or different embodiments. Furthermore, the terms “comprising,” “including,” “having,” and the like, as used with respect to embodiments of the present invention, are synonymous.
According to various embodiments of the present invention, adaptive read and write memory systems and methods are provided that may adapt to changes to level distributions of multi-level memory cells. In particular, the novel systems may be configured to compute new mean values of level distributions and/or new detection threshold values of multi-level memory cells (herein “memory cells”) after the memory cells have degraded as a result of, for example, repeated cycling. For purposes of the following description, the computed new mean and detection threshold values will be referred to as “optimal or near optimal” values. That is, the phrase “optimal or near optimal” as used herein are in reference to the mew mean and detection threshold values that may be computed using either an optimal solution, which may be a more complex solution requiring more computational power, or a simpler near optimal (approximate) solution. As will be described herein, these values may be calculated whenever a read operation is performed on multi-level memory cells.
FIG. 3 depicts an adaptive read/write memory system in accordance with various embodiments of the present invention. The read/write memory system (herein “system”) 10, as in conventional read/write memory systems, may include memory cells 12, which may be multi-level flash memory cells, a read block 14, a signal processing and decoding block 18, and write block 20. However, unlike conventional systems, the system 10 may further include a mean and standard deviation estimation block 22, a computation block 24, and a look-up table 26, operationally coupled together. One or more of the components depicted, such as the mean and standard deviation estimation block 22 and the computation block 24, may be implemented using hardware components, such as, for example, application specific integrated circuit (ASIC), and/or software.
In brief, and as will be described in greater detail herein, the mean and standard deviation estimation block (herein “estimation block”) 22 may be configured to calculate estimated mean and standard deviation values of level distributions of the memory cells 12 during, for example, a read operation of the memory cells 12. The computation block 24 may be configured to compute optimal or near optimal mean and detection threshold values based on the estimated mean and standard deviation values provided by the estimation block 22. As will be described herein, the optimal or near optimal mean values computed may be used in order to adaptively write data into the memory cells 12 while the detection threshold values computed may be used in order to adaptively read data stored in the memory cells 12.
The optimal or near optimal mean and detection threshold values computed by the computation block 24 and the estimated mean and standard deviation values calculated by the estimation block 22, in various embodiments of the present invention, may be stored in the look-up table 26. The read block 14, the signal processing and decoding block 18, and the write block 20 may use selected values stored in and provided by the look-up table 26 to perform various operations. Alternatively, such values may be directly provided by the estimation block 22 and the computation block 24 as indicated by reference 28.
As briefly described above, the estimation block 22 may calculate estimated mean and standard deviation values of level distributions of the memory cells 12 during a read operation. The estimated mean and standard deviation values may be calculated as intermediate steps in order for the computation block 24 to eventually compute the optimal or near optimal mean and detection threshold values based, at least in part, on the estimated mean and standard deviation values. Estimated mean and standard deviation values may be calculated for each level distribution of a memory cell. In various embodiments and as will be described herein, the estimation of the mean and standard deviation values may be achieved via training or online adaptation.
For example, the estimated mean and standard deviation values may be calculated by using pilot memory cells having known or predetermined data. That is, certain memory cells 12 may be designated as pilot memory cells, where the data stored in these memory cells are predefined and known. The read block 14 may then exploit these pilot cells for estimating the mean and standard deviations as described in, for example, co-pending U.S. patent application Ser. No. 11/738,263, filed Apr. 20, 2007, entitled “Channel Estimation for Multi-Level Flash Memories Using Pilot Signals,” which is hereby incorporated by reference in its entirety for all purposes. Such a method for estimating the mean and standard deviations (i.e., variances) of level distributions of a memory cell is referred to herein as a “training” technique. Alternatively, online adaptation techniques may be employed for estimating the mean and standard deviations. For instance, the LMS (least-mean-squares) algorithm may be used to estimate the mean and standard deviations based on the data recovered from the memory cells.
Based on the estimated mean and standard deviation values calculated by the estimation block 22, the computation block 24 may compute optimal or near optimal mean and detection threshold values for a memory cell or a group of memory cells. Specifics on how the optimal or near optimal mean and detection threshold values may be computed will be described in detail herein. The calculated optimal or near optimal mean and detection threshold values may then be stored in look-up table 26. An example of look-up table 26 is depicted in FIG. 4. In particular, FIG. 4 depicts an exemplary table structure 40 of look-up table 26 of FIG. 3, in accordance with various embodiments of the present invention.
In table structure 40, the “block index” column 41 on the far left is in reference to a block of memory cells. The second and third columns 42 and 43 from the left are for estimated mean and standard deviation values calculated by the estimation block 22. The two columns 44 and 45 on the right are for the optimal or near optimal mean and detection threshold values as computed by the computation block 24. Thus, in this example, a block or a group of memory cells may be associated with common estimated mean and standard deviation values, as well as common optimal or near optimal mean and detection threshold values.
Since a multi-level memory cell may have multiple level distributions, multiple estimated mean and standard deviation values, as well as multiple optimal or near optimal mean and detection threshold values may be calculated and stored in the table structure 40 for a memory cell (in this example, for each group of memory cells). Thus, for the second and third columns, the “estimated mean values” column and the “estimated standard deviations” column, there are multiple mean (m0, m1, . . . ) and standard deviation (σ0, σ1, . . . ) values for each of the multiple levels of a multi-level cell (see FIG. 1). Similarly, there may be multiple values included in the fourth and fifth columns 44 and 45 for the “near optimal mean values” and the “near optimal detection thresholds.”
Although the computation block 24, in some embodiments of the present invention, may compute the optimal or near optimal mean and detection threshold values as soon as the estimated mean and standard deviations are calculated by the estimation block 22, as appears to be case in the above described embodiment, in some alternative embodiments, the optimal or near optimal mean and detection threshold values may be computed at some later time after the estimated mean values and estimated standard deviation values have already been stored in the table 26.
As described previously, the optimal or near optimal detection threshold values computed may be used during a read operation of one or more multi-level memory cells. In using the optimal or near optimal detection threshold values, error resulting from reading a degraded memory cell as a result of repeated cycling may be minimized.
Referring now to FIG. 5, an adaptive signal detection flow process for reading multi-level memory cells based, at least in part, on computed optimal or near optimal detection threshold values in accordance with various embodiments of the present invention is illustrated. The process 50 may begin when the latest estimated mean values and standard deviation values for level distributions of the multi-level memory cells are calculated at 52. The latest estimated mean values and standard deviation values may be calculated based on, for example, readback signals of pilot memory cells or online adaptation techniques as described previously.
Based at least on the estimated mean values and standard deviation values, optimal or near optimal detection threshold values may be computed at 54. The optimal or near optimal detection threshold values may be computed using a linear solution, which may be an approximate or near optimal solution, or Newton's method, which may be an optimal solution), both of which will be described in greater detail herein. Signal detection (i.e., reading) of the readback signal may then be performed using the estimated mean values and standard deviation values, and the computed near optimal detection threshold values at 56.
In contrast to the computed optimal or near optimal detection threshold values, the computed optimal or near optimal mean values may be used during a write operation to program a memory cell. That is, although optimal or near optimal mean values may be calculated (along with the near optimal detection threshold values) for a multi-level memory cell during or after a read operation of the multi-level memory cell, the optimal or near optimal mean values may not be used until a subsequent write (i.e., programming) operation of the memory cell. The computed optimal or near optimal mean values may be used to more reliably program multi-level memory cells during a write operation, particularly for example, those multi-level memory cells that have been repeatedly cycled.
FIG. 6 illustrates an adaptive flow process for a write operation of multi-level memory cells, according to various embodiments of the present invention. The process 60 may begin when a determination is made as to whether new optimal or near optimal mean values for the memory cells are available at 62. In some instances, new optimal or near optimal mean values for the memory cells may have already been calculated and stored, such as those that may be stored in look-up table 26. If so, the memory cells may be programmed (i.e., written) in accordance with the new optimal or near optimal mean values at 68. If not, then the most recent estimates of mean and standard deviation values for the level distributions of the memory cells are obtained either from the look-up table 26, or are calculated at 64. Based on the estimated mean and standard deviation values, new optimal or near optimal mean values are computed at 66 using a linear solution (i.e., near optimal solution) or Newton's method (i.e., optimal solution). After computing the optimal or near optimal mean values, the memory cells may be programmed (i.e., written) according to the optimal or near optimal mean values at 68.
In order to calculate optimal or near optimal mean and detection threshold values of multi-level memory cells, it is recognized that many parameters associated with multi-level memory cells including, for example, the means and standard deviations of the lowest and highest level distributions (e.g., the level distributions of level 0 and level 4 of FIG. 1) are functions of the memory cells and are not easily controllable. Given these values, however, it may be possible to optimize a read/write memory system by adjusting the means of level distributions (except for the means associated with the lowest and highest level distributions) during programming of the memory cells. In addition, for hard decision detection in a read operation, the detection threshold values may be optimized according to the current or most recent level distributions for minimum probability of error. In other words, optimal or near optimal mean and detection threshold values may be determined for write and read operations of multi-level memory cells in order to reduce error if certain parameters such as the means and standard deviations of the lowest and highest level distributions are assumed to be determined by the device characteristics.
In order to obtain an optimal solution for computing optimal mean and detection threshold values, reference is now made to FIG. 7, which depicts an exemplary 3-level threshold voltage distribution (“3-level distribution”). For this example, the respective standard deviations for the 3 level distributions are σ0, σ1, and σ2. For purposes of illustration, assume that the level distribution means fall in the range between 0 and L, where L is the distance or range between the mean of the lowest level distribution (i.e., the left-most level distribution in FIG. 7) and the mean of the highest level distribution (i.e., the right-most level distribution in FIG. 7). Then the mean of the middle level distribution may be denoted by d. Once d is given, the crossing points of the probability density functions (pdf's) may also be determined. In this example the distances of the crossing points from d are denoted by x0(d) and x1(d), respectively.
From signal detection theory, it is known that the optimal detection thresholds for multi-level memory cells are the crossing points of the pdf's. In the following, the pdf's are first shown to maintain the same value at the crossing points (i.e., detection thresholds) when minimum probability of error is achieved. The shaded areas in FIG. 7 correspond to the error regions. It is not difficult to see that the probability of making an error in detection is given by




  
    P
    e
  
  =
  
    
      
        ∫
        
          
            d
            -
            
              
                x
                0
              
              ⁡
              
                (
                d
                )
              
            
          
          
            σ
            0
          
        
        ∞
      
      ⁢
      
        N
        ⁡
        
          (
          
            0
            ,
            1
          
          )
        
      
    
    +
    
      
        ∫
        
          
            
              x
              0
            
            ⁡
            
              (
              d
              )
            
          
          
            σ
            1
          
        
        ∞
      
      ⁢
      
        N
        ⁡
        
          (
          
            0
            ,
            1
          
          )
        
      
    
    +
    
      
        ∫
        
          
            
              x
              1
            
            ⁡
            
              (
              d
              )
            
          
          
            σ
            1
          
        
        ∞
      
      ⁢
      
        N
        ⁡
        
          (
          
            0
            ,
            1
          
          )
        
      
    
    +
    
      
        ∫
        
          
            L
            -
            d
            -
            
              
                x
                1
              
              ⁡
              
                (
                d
                )
              
            
          
          
            σ
            2
          
        
        ∞
      
      ⁢
      
        N
        ⁡
        
          (
          
            0
            ,
            1
          
          )
        
      
    
  




where N(0, 1) denotes the standard Gaussian distribution function with zero mean and variance 1. Taking the derivative of Pe with respect to d, it follows that




  
    
      
        
          
            ∂
            
              P
              e
            
          
          
            ∂
            d
          
        
        =
        
          0
          =
          
            
              
                g
                ⁡
                
                  (
                  
                    
                      d
                      -
                      
                        
                          x
                          0
                        
                        ⁡
                        
                          (
                          d
                          )
                        
                      
                    
                    
                      σ
                      0
                    
                  
                  )
                
              
              ·
              
                1
                
                  σ
                  0
                
              
              ·
              
                
                  ∂
                  
                    (
                    
                      d
                      -
                      
                        
                          x
                          0
                        
                        ⁡
                        
                          (
                          d
                          )
                        
                      
                    
                    )
                  
                
                
                  ∂
                  d
                
              
            
            +
            
              
                g
                ⁡
                
                  (
                  
                    
                      
                        x
                        0
                      
                      ⁡
                      
                        (
                        d
                        )
                      
                    
                    
                      σ
                      1
                    
                  
                  )
                
              
              ·
              
                1
                
                  σ
                  1
                
              
              ·
              
                
                  ∂
                  
                    (
                    
                      
                        x
                        0
                      
                      ⁡
                      
                        (
                        d
                        )
                      
                    
                    )
                  
                
                
                  ∂
                  d
                
              
            
            +
            
              
                g
                ⁡
                
                  (
                  
                    
                      
                        x
                        1
                      
                      ⁡
                      
                        (
                        d
                        )
                      
                    
                    
                      σ
                      1
                    
                  
                  )
                
              
              ·
              
                1
                
                  σ
                  1
                
              
              ·
              
                
                  ∂
                  
                    (
                    
                      
                        x
                        1
                      
                      ⁡
                      
                        (
                        d
                        )
                      
                    
                    )
                  
                
                
                  ∂
                  d
                
              
            
            +
            
              
                g
                ⁡
                
                  (
                  
                    
                      L
                      -
                      d
                      -
                      
                        
                          x
                          1
                        
                        ⁡
                        
                          (
                          d
                          )
                        
                      
                    
                    
                      σ
                      2
                    
                  
                  )
                
              
              ·
              
                1
                
                  σ
                  2
                
              
              ·
              
                
                  ∂
                  
                    (
                    
                      L
                      -
                      d
                      -
                      
                        
                          x
                          1
                        
                        ⁡
                        
                          (
                          d
                          )
                        
                      
                    
                    )
                  
                
                
                  ∂
                  d
                
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          1
          )
        
      
    
  




where




  
    g
    ⁡
    
      (
      x
      )
    
  
  =
  
    
      1
      
        
          2
          ⁢
          π
        
      
    
    ⁢
    
      
        e
        
          
            -
            
              x
              2
            
          
          /
          2
        
      
      .
    
  




Noting that



  
    
      
        g
        ⁡
        
          (
          
            
              d
              -
              
                
                  x
                  0
                
                ⁡
                
                  (
                  d
                  )
                
              
            
            
              σ
              0
            
          
          )
        
      
      ·
      
        1
        
          σ
          0
        
      
    
    =
    
      
        
          g
          ⁡
          
            (
            
              
                
                  x
                  0
                
                ⁡
                
                  (
                  d
                  )
                
              
              
                σ
                1
              
            
            )
          
        
        ·
        
          1
          
            σ
            1
          
        
      
      ⁢
      
        
      
      ⁢
      and
    
  




  
    
      
        g
        ⁡
        
          (
          
            
              
                x
                1
              
              ⁡
              
                (
                d
                )
              
            
            
              σ
              1
            
          
          )
        
      
      ·
      
        1
        
          σ
          1
        
      
    
    =
    
      
        g
        ⁡
        
          (
          
            
              L
              -
              d
              -
              
                
                  x
                  1
                
                ⁡
                
                  (
                  d
                  )
                
              
            
            
              σ
              2
            
          
          )
        
      
      ·
      
        1
        
          σ
          2
        
      
    
  



to obtain




  
    
      g
      ⁡
      
        (
        
          
            d
            -
            
              
                x
                0
              
              ⁡
              
                (
                d
                )
              
            
          
          
            σ
            0
          
        
        )
      
    
    ·
    
      1
      
        σ
        0
      
    
  
  =
  
    
      g
      ⁡
      
        (
        
          
            L
            -
            d
            -
            
              
                x
                1
              
              ⁡
              
                (
                d
                )
              
            
          
          
            σ
            2
          
        
        )
      
    
    ·
    
      1
      
        σ
        2
      
    
  




The above equation is the result after canceling out the equal terms in Eq. (1). This completes the proof. Extensions to more than 3 levels are similar. For ease of reference, the above property will be referred to, herein, as “Equal Value Property.”
As an example, the “Equal Value Property” may be used to solve the optimization problem for a 4-level distribution, which is depicted in FIG. 8. For the 4-level distribution example depicted in FIG. 8, the mean range is between 0 and L, the means of the level distributions are denoted as m0, m1, m2, and m3, the detection thresholds for the level distributions are denoted as t1, t2, and t3, and the respective standard deviations of the level distributions are denoted as σ0, σ1, σ2, and σ3. The distances between the means and the detection thresholds, as depicted, are x, y, z, and w. In order to determine the optimal solutions for x, y, z, and w, the following set of equations may be used:




  {
  
    
      
        
          
            
              
                1
                
                  σ
                  0
                
              
              ⁢
              
                e
                
                  -
                  
                    
                      x
                      2
                    
                    
                      2
                      ⁢
                      
                        σ
                        0
                        2
                      
                    
                  
                
              
            
            =
            
              
                
                  1
                  
                    σ
                    1
                  
                
                ⁢
                
                  e
                  
                    -
                    
                      
                        y
                        2
                      
                      
                        2
                        ⁢
                        
                          σ
                          1
                          2
                        
                      
                    
                  
                
              
              =
              
                
                  
                    1
                    
                      σ
                      2
                    
                  
                  ⁢
                  
                    e
                    
                      -
                      
                        
                          z
                          2
                        
                        
                          2
                          ⁢
                          
                            σ
                            2
                            2
                          
                        
                      
                    
                  
                
                =
                
                  
                    1
                    
                      σ
                      3
                    
                  
                  ⁢
                  
                    e
                    
                      -
                      
                        
                          w
                          2
                        
                        
                          2
                          ⁢
                          
                            σ
                            3
                            2
                          
                        
                      
                    
                  
                
              
            
          
        
      
      
        
          
            
              x
              +
              
                2
                ⁢
                y
              
              +
              
                2
                ⁢
                z
              
              +
              w
            
            =
            L
          
        
      
    
    .
  




Numerical approaches may be used to solve the above equations. For example, the constraint in Eq. (2) may be integrated into Eq. (3) by defining a constant C viz:




  
    
      
        {
        
          
            
              
                
                  
                    1
                    
                      σ
                      0
                    
                  
                  ⁢
                  
                    e
                    
                      -
                      
                        
                          x
                          2
                        
                        
                          2
                          ⁢
                          
                            σ
                            0
                            2
                          
                        
                      
                    
                  
                
                =
                
                  
                    
                      1
                      
                        σ
                        1
                      
                    
                    ⁢
                    
                      e
                      
                        -
                        
                          
                            y
                            2
                          
                          
                            2
                            ⁢
                            
                              σ
                              1
                              2
                            
                          
                        
                      
                    
                  
                  =
                  
                    
                      
                        1
                        
                          σ
                          2
                        
                      
                      ⁢
                      
                        e
                        
                          -
                          
                            
                              z
                              2
                            
                            
                              2
                              ⁢
                              
                                σ
                                2
                                2
                              
                            
                          
                        
                      
                    
                    =
                    
                      
                        
                          1
                          
                            σ
                            3
                          
                        
                        ⁢
                        
                          e
                          
                            -
                            
                              
                                w
                                2
                              
                              
                                2
                                ⁢
                                
                                  σ
                                  3
                                  2
                                
                              
                            
                          
                        
                      
                      =
                      C
                    
                  
                
              
            
          
          
            
              
                
                  x
                  +
                  
                    2
                    ⁢
                    y
                  
                  +
                  
                    2
                    ⁢
                    z
                  
                  +
                  w
                
                =
                L
              
            
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          4
          )
        
      
    
  




Now, C may be found by using Newton's method to solve



  
    
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          =
          
          ⁢
          
            x
            +
            
              2
              ⁢
              y
            
            +
            
              2
              ⁢
              z
            
            +
            w
            -
            L
          
        
      
    
    
      
        
          =
          
          ⁢
          
            
              g
              ⁡
              
                (
                
                  
                    σ
                    0
                  
                  ,
                  C
                
                )
              
            
            +
            
              2
              ⁢
              
                g
                ⁡
                
                  (
                  
                    
                      σ
                      1
                    
                    ,
                    C
                  
                  )
                
              
            
            +
            
              2
              ⁢
              
                g
                ⁡
                
                  (
                  
                    
                      σ
                      2
                    
                    ,
                    C
                  
                  )
                
              
            
            +
            
              g
              ⁡
              
                (
                
                  
                    σ
                    3
                  
                  ,
                  C
                
                )
              
            
            -
            
              L
              .
            
          
        
      
    
  




  
    
      where
      ⁢
      
        
      
      ⁢
      g
      ⁢
      
        (
        
          σ
          ,
          C
        
        )
      
    
    =
    
      σ
      ⁢
      
        
          
            -
            2
          
          ⁢
          
            ln
            ⁡
            
              (
              
                σ
                ⁢
                
                  
                
                ⁢
                C
              
              )
            
          
        
      
    
  



After taking the derivative of f(x) with respect to C, the following is obtained:




  
    
      
        
          
            d
            
              f
              ⁡
              
                (
                x
                )
              
            
          
          /
          
            d
            C
          
        
        =
        
          
            d
            
              (
              
                x
                +
                
                  2
                  ⁢
                  y
                
                +
                
                  2
                  ⁢
                  z
                
                +
                w
                -
                L
              
              )
            
          
          /
          
            d
            C
          
        
      
    
  
  
    
      
        =
        
          
            
              -
              
                σ
                0
              
            
            
              g
              ⁡
              
                (
                
                  
                    σ
                    0
                  
                  ,
                  C
                
                )
              
            
          
          +
          
            
              
                -
                2
              
              ⁢
              
                σ
                1
              
            
            
              g
              ⁡
              
                (
                
                  
                    σ
                    1
                  
                  ,
                  C
                
                )
              
            
          
          +
          
            
              
                -
                2
              
              ⁢
              
                σ
                2
              
            
            
              g
              ⁡
              
                (
                
                  
                    σ
                    2
                  
                  ,
                  C
                
                )
              
            
          
          +
          
            
              
                -
                
                  σ
                  3
                
              
              
                g
                ⁡
                
                  (
                  
                    
                      σ
                      3
                    
                    ,
                    C
                  
                  )
                
              
            
            .
          
        
      
    
  




Now, C can be found through the following iteration




  
    C
    
      n
      +
      1
    
  
  =
  
    
      C
      n
    
    -
    
      
        
          f
          ⁡
          
            (
            
              C
              n
            
            )
          
        
        
          
            f
            ′
          
          ⁡
          
            (
            
              C
              n
            
            )
          
        
      
      .
    
  




The solution converges very fast for good initial values of C. Once C is solved, the x, y, z, and w values may easily follow. Thus, equation (2) and (3) may be used in order to obtain what has been referred to previously as, the optimal solution.
Although, the above solution may be very accurate, approximations may be used to simplify the computation and obtain a near optimal solution. For example, taking the logarithm of the first equation of Eq. (2), the following may be obtained:




  
    
      
        
          
            x
            2
          
          
            2
            ⁢
            
              σ
              0
              2
            
          
        
        =
        
          
            
              y
              2
            
            
              2
              ⁢
              
                σ
                1
                2
              
            
          
          +
          
            
              ln
              ⁡
              
                (
                
                  
                    σ
                    1
                  
                  
                    σ
                    0
                  
                
                )
              
            
            .
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          5
          )
        
      
    
  




Note that for practical applications, the ratio between σ1 and σ0 is close to 1, which after taking the logarithm are negligible comparing to other terms in Eq. (5). Ignoring the term ln




  
    (
    
      
        σ
        1
      
      
        σ
        0
      
    
    )
  
  ,




it follows that





  
    x
    
      σ
      0
    
  
  =
  
    
      y
      
        σ
        1
      
    
    .
  




Similar approximations also hold true for z and w. Thus, the near optimal solutions may be obtained as




  {
  
    
      
        
          
            x
            =
            
              
                
                  σ
                  0
                
                
                  
                    σ
                    0
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      1
                    
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      2
                    
                  
                  +
                  
                    σ
                    3
                  
                
              
              ⁢
              L
            
          
        
      
      
        
          
            y
            =
            
              
                
                  σ
                  1
                
                
                  
                    σ
                    0
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      1
                    
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      2
                    
                  
                  +
                  
                    σ
                    3
                  
                
              
              ⁢
              L
            
          
        
      
      
        
          
            z
            =
            
              
                
                  σ
                  2
                
                
                  
                    σ
                    0
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      1
                    
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      2
                    
                  
                  +
                  
                    σ
                    3
                  
                
              
              ⁢
              L
            
          
        
      
      
        
          
            w
            =
            
              
                
                  σ
                  3
                
                
                  
                    σ
                    0
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      1
                    
                  
                  +
                  
                    2
                    ⁢
                    
                      σ
                      2
                    
                  
                  +
                  
                    σ
                    3
                  
                
              
              ⁢
              L
            
          
        
      
    
    .
  




These solutions may be referred to, for purposes of this description as linearized solutions and they are near optimal under practical operating conditions. That is, once x, y, z, and w are solved, approximate (i.e., near optimal) solutions may be determined for the mean and detection threshold values of 4-level memory cells. It should be noted, however, that the above formulation may be easily extended to any M number of levels for M≧2.
Based on the above near optimal solution, and once the estimated mean and standard deviations for all the distributions have been determined using, for example, pilot cells, it may be straightforward to find the near optimal detection threshold values. For example, if the estimated mean and standard deviation values of the level distributions for a 4-level (2 bit/cell) flash memory device are found to be {mi, i=0, 1, 2, 3} and {σi, i=0, 1, 2, 3}, respectively. Utilizing the linearized solution (i.e., near optimal solution), one may obtain:




   
  
    {
    
      
        
          
            
              t
              1
            
            =
            
              
                
                  
                    m
                    1
                  
                  ⁢
                  
                    σ
                    0
                  
                
                +
                
                  
                    m
                    0
                  
                  ⁢
                  
                    σ
                    1
                  
                
              
              
                
                  σ
                  0
                
                +
                
                  σ
                  1
                
              
            
          
        
      
      
        
          
            
              t
              2
            
            =
            
              
                
                  
                    m
                    2
                  
                  ⁢
                  
                    σ
                    1
                  
                
                +
                
                  
                    m
                    1
                  
                  ⁢
                  
                    σ
                    2
                  
                
              
              
                
                  σ
                  1
                
                +
                
                  σ
                  2
                
              
            
          
        
      
      
        
          
            
              t
              3
            
            =
            
              
                
                  
                    m
                    3
                  
                  ⁢
                  
                    σ
                    2
                  
                
                +
                
                  
                    m
                    2
                  
                  ⁢
                  
                    σ
                    3
                  
                
              
              
                
                  σ
                  2
                
                +
                
                  σ
                  3
                
              
            
          
        
      
    
  




where t1, t2 and t3 are the near optimal signal detection thresholds. Thus, the near optimal detection thresholds, t1, t2 and t3, of a 4-level multi-level memory may be solved using the above equations. If more accuracy is desired, exact solutions may always be found by utilizing the Equal Value Property through the Newton's method illustrated previously.
The near optimal mean values for multi-level memory cells may be obtained as follows. For a M level memory cell, assume that the estimated voltage means are denoted as {mi, i=0, 1, . . . M−1} and the corresponding standard deviations as {σi, i=0, 1, . . . M−1}. Due to physical reasons and as previously alluded to, the mean values corresponding to the lowest (m0) and highest levels (mM-1) and the standard deviations are not easily controllable, thus such values are assumed to be predefined and set. However, the values of mi, i=1, 2, . . . , M−2 may be adjusted for optimal performance. By denoting L=mM-1−m0 and utilizing the linearized solutions, the near optimal mean value ({tilde over (m)}) for the i-th level is given by




  
    
      
        
          
            m
            ~
          
          i
        
        =
        
          
            m
            0
          
          +
          
            
              
                
                  σ
                  0
                
                +
                
                  2
                  ⁢
                  
                    
                      ∑
                      
                        k
                        =
                        1
                      
                      
                        i
                        -
                        1
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      σ
                      k
                    
                  
                
                +
                
                  σ
                  i
                
              
              
                
                  σ
                  0
                
                +
                
                  2
                  ⁢
                  
                    
                      ∑
                      
                        k
                        =
                        1
                      
                      
                        M
                        -
                        2
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      σ
                      k
                    
                  
                
                +
                
                  σ
                  
                    M
                    -
                    1
                  
                
              
            
            ⁢
            L
          
        
      
    
    
      
        Eq
        .
        
          
        
        ⁢
        
          (
          6
          )
        
      
    
  




Thus, the above equation may be used to solve for near optimal mean values of M-level distribution memory cells in accordance with various embodiments of the present invention. Alternatively, the Newton's method may be used for more accuracy.
Although specific embodiments have been illustrated and described herein, it will be appreciated by those of ordinary skill in the art and others, that a wide variety of alternate and/or equivalent implementations may be substituted for the specific embodiments illustrated and described without departing from the scope of the present invention. This application is intended to cover any adaptations or variations of the embodiments discussed herein. Therefore, it is manifested and intended that various embodiments of the invention be limited only by the claims and the equivalents thereof.

## Claims

1. An apparatus, comprising:
multi-level memory cells;
an estimation block configured to determine estimated mean and standard deviation values of level distributions of the multi-level memory cells; and
a computation block operatively coupled to the estimation block and configured to compute at least optimal or near optimal detection threshold values based, at least in part, on the estimated mean and standard deviation values, the optimal or near optimal detection threshold values to be used in order to facilitate reading of data stored in the multi-level memory cells,
wherein the multi-level memory cells include at least one M level memory cell having M levels, and the computation block is further configured to compute an ith level near optimal mean value ({tilde over (m)}i) of the M level memory cell according to the equation:





  
    
      m
      ~
    
    i
  
  =
  
    
      m
      0
    
    +
    
      
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  i
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            i
          
        
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  M
                  -
                  2
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            
              M
              -
              1
            
          
        
      
      ⁢
      L
    
  



where mi is the estimated mean value of the ith level of the M level memory cell, σi is the estimated standard deviation value of the ith level of the M level memory cell, L is equal to mM-1-mo, and wherein at least one of the estimated standard deviation values σ0, . . . , σM-1 is non-zero.

2. The apparatus of claim 1, wherein the multi-level memory cells include one or more pilot cells for storing therein predetermined data, and the estimation block is further configured to determine the estimated mean and standard deviation values using the one or more pilot cells.

3. The apparatus of claim 1, wherein the the computation block is further configured to compute near optimal detection threshold values of the M-level memory cell according to the equation:




  
    t
    i
  
  =
  
    
      
        
          m
          
            i
            -
            1
          
        
        ⁢
        
          σ
          i
        
      
      +
      
        
          m
          
            i
            -
            1
          
        
        ⁢
        
          σ
          
            i
            -
            1
          
        
      
    
    
      
        σ
        i
      
      +
      
        σ
        
          i
          -
          1
        
      
    
  



for 1≦i≦M−1, where ti is the detection threshold between level i−1 and level i; mi, 0≦i≦M−1 are the estimated means of the M-level distributions; and σi, 0≦i≦M≦1 are the estimated standard deviations of the M-level distributions.

4. The apparatus of claim 1, further comprising a read block configured to read data in the multi-level memory cells based, at least in part, on the optimal or near optimal detection threshold values.

5. The apparatus of claim 4, further comprising a look-up table configured to provide the optimal or near optimal detection threshold values to the read block.

6. The apparatus of claim 4, further comprising a signal processing and decoding block configured to receive the optimal or near optimal detection threshold values from the look-up table, the optimal or near optimal detection threshold values to be used by the signal processing and decoding block to facilitate processing and decoding of data read from the multi-level memory cells.

7. The apparatus of claim 1, wherein the multi-level memory cells are M-level memory cells, where M is greater than or equal to 2.

8. The apparatus of claim 1, wherein the multi-level memory cells are multi-level flash memory cells.

9. An apparatus, comprising:
multi-level memory cells; and
a computation block configured to compute optimal or near optimal mean and detection threshold values based, at least in part, on estimated mean and standard deviation values of level distributions of the multi-level memory cells, the optimal or near optimal mean and detection threshold values to be used to facilitate writing and reading, respectively, of data to and from the multi-level memory cells,
wherein the multi-level memory cells include at least one M level memory cell having M levels, and the computation block is further configured to compute an ith level near optimal mean value ({tilde over (m)}i) of the M level memory cell according to the equation:





  
    
      m
      ~
    
    i
  
  =
  
    
      m
      0
    
    +
    
      
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  i
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            i
          
        
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  M
                  -
                  2
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            
              M
              -
              1
            
          
        
      
      ⁢
      L
    
  



where mi is the estimated mean value of the ith level of the M level memory cell, σi is the estimated standard deviation value of the ith level of the M level memory cell, and L is equal to mM-1-m0, and wherein at least one of the estimated standard deviation values σ0, . . . , σM-1 is non-zero.

10. The apparatus of claim 9, further comprising a look-up table configured to store the estimated mean and standard deviation values and the optimal or near optimal mean and detection threshold values.

11. The apparatus of claim 10, further comprising a read block configured to read data from the multi-level memory cells, the reading of the data being based, at least in part, on the stored optimal or near optimal detection threshold values.

12. The apparatus of claim 10, further comprising a write block configured to write data to the multi-level memory cells, the writing of the data based, at least in part, on the stored optimal or near optimal mean values.

13. The apparatus of claim 9, wherein the computation block is further configured to compute near optimal detection threshold values of the M-level memory cell according to the equation:




  
    t
    i
  
  =
  
    
      
        
          m
          
            i
            -
            1
          
        
        ⁢
        
          σ
          i
        
      
      +
      
        
          m
          
            i
            -
            1
          
        
        ⁢
        
          σ
          
            i
            -
            1
          
        
      
    
    
      
        σ
        i
      
      +
      
        σ
        
          i
          -
          1
        
      
    
  



for 1≦i≦M−1, where ti is the detection threshold between level i−1 and level i; mi, 0≦i≦M−1 are the estimated means of the M-level distributions; and σi, 0≦i≦M−1 are the estimated standard deviations of the M-level distributions.

14. A method, comprising:
determining estimated mean and standard deviation values of level distributions of multi-level memory cells; and
computing optimal or near optimal mean values based, at least in part, on the estimated mean and standard deviation values, the optimal or near optimal mean values to be used to facilitate writing of data to the multi-level memory cells,
wherein the multi-level memory cells include at least one M level memory cell having M levels, and said computing comprises computing an ith level near optimal mean value ({tilde over (m)}) of the M level memory cell according to the equation:





  
    
      m
      ~
    
    i
  
  =
  
    
      m
      0
    
    +
    
      
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  i
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            i
          
        
        
          
            σ
            0
          
          +
          
            2
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                
                  M
                  -
                  2
                
              
              ⁢
              
                
              
              ⁢
              
                σ
                k
              
            
          
          +
          
            σ
            
              M
              -
              1
            
          
        
      
      ⁢
      L
    
  



where mi is the estimated mean value of the ith level of the M level memory cell, σi is the estimated standard deviation value of the ith level of the M level memory cell, and L is equal to mM-1-m0, and wherein at least one of the estimated standard deviation values σ0, . . . , σM-1 is non-zero.

15. The method of claim 14, wherein the multi-level memory cells include one or more pilot cells for storing predetermined data, and said determining comprises determining the estimated mean and standard deviation values by using the one or more pilot cells.

16. The method of claim 14, further comprising computing optimal or near optimal detection threshold values based, at least in part, on the estimated mean and standard deviation values, the optimal or near optimal detection threshold values to be used to facilitate reading of data from the multi-level memory cells.

17. The method of claim 16, wherein said computing optimal or near optimal detection threshold values further comprises computing optimal or near optimal detection threshold values of the M-level memory cell according to the equation:




  
    t
    i
  
  =
  
    
      
        
          m
          
            i
            -
            1
          
        
        ⁢
        
          σ
          i
        
      
      +
      
        
          m
          
            i
            -
            1
          
        
        ⁢
        
          σ
          
            i
            -
            1
          
        
      
    
    
      
        σ
        i
      
      +
      
        σ
        
          i
          -
          1
        
      
    
  



for 1≦i≦M−1, where ti is the detection threshold between level i−1 and level i; mi, 0≦i≦M−1 are the estimated means of the M-level distributions; and σi,0≦i≦M−1 are the estimated standard deviations of the M-level distributions.

18. The method of claim 16, further comprising reading data from the multi-level memory cells based, at least in part, on the optimal or near optimal detection threshold values.

19. The method of claim 18, further comprising storing the optimal or near optimal detection threshold values to a look-up table, and said reading of data from the multi-level memory cells comprises using the optimal or near optimal detection threshold values stored in the lookup table to facilitate said reading of the data in the multi-level memory cells.

20. The method of claim 19, further comprising signal processing and decoding of data read from the multi-level memory cells, the signal processing and decoding being based, at least in part, on the optimal or near optimal detection threshold values stored in the look-up table.

21. The method of claim 14, further comprising writing data into the multi-level memory cells based, at least in part, on the optimal or near optimal mean values.

22. The method of claim 21, further comprising storing the optimal or near optimal mean values to a look-up table, and said writing is based, at least in part, on the stored optimal or near optimal mean values.

23. The method of claim 14, further comprising storing the determined estimated mean and standard deviation values in a look-up table.