Commonly assigned U.S. Pat. No. 5,442,315 is incorporated for its showing of phase locked loops used in data detection.
This invention relates to multi-channel data detection, and, more particularly, to sample timing of multiple read channels, wherein the signal clocking of the signals received by the multiple read channels are correlated.
Data storage having multiple channels typically employs movable media in which data is recorded on one pass of the media, and the data is read back and detected at a subsequent time, possibly on a different pass of the media and possibly on a different drive than that which recorded the data. One example of movable media is magnetic tape which has a plurality of parallel tracks for recording. The parallel tracks are typically written simultaneously such that the signal clocking of the signals written on parallel media are correlated.
Symbol timing recovery during readback represents one of the most critical functions in data storage read channels. Sampling an analog readback signal at the right time instant is important for achieving good overall performance. Among the challenges presented are the presence of disturbances such as dropout events, instantaneous speed variations, and signal distortion of various origins, which make timing recovery difficult. As areal recording densities become higher and higher, SNR (signal to noise ratio) margins are decreased, making satisfactory timing recovery an even more challenging task.
Timing recovery typically is based on a PLL (phase locked loop) for each channel whose purpose is to accurately estimate the timing offsets before sampling the analog signal. Problems in this context, especially with degraded channel conditions as explained above, are those of temporary “loss of lock” or “cycle slip”. These terms refer to a phenomenon where the phase adjustments of the timing control loop stabilize around an undesirable phase value that is located one or several symbol interval durations apart from the desired operating point. This phenomenon often results in long bursts of bit and symbol errors which may exceed the error correction capability of the error correction codes, such as Reed-Solomon codes, leading to severe performance degradation or even permanent error situations.
The conventional approaches have been to optimize the operation of the PLL for maximum noise rejection and loop jitter minimization, to increase the loop robustness by using more reliable decisions, etc.
The incorporated '315 patent, inter alia, takes advantage of the multiple channels by providing a global clock which generates a global average frequency signal by averaging the content of the frequency registers employed by the PLLs of the individual channels. The global average frequency is employed by the PLL of each track which adds its own scaled phase error to it.
Common sample timing control for and methods related to sample timing of multiple read channels, wherein the signal clocking of the signals received by the multiple read channels are correlated.
In one embodiment, a common sample timing control comprises multiple phase error inputs, each indicating phase error of one of the read channels. Logic responsive to the multiple phase error inputs is configured to weight and crosscouple the phase error indication of each phase error input with the phase error indication of each other phase error input, and to apply common gain. Feedback logic is responsive to the crosscoupling and is configured to provide a sample timing phase estimate for each read channel.
In a further embodiment, the weighting and crosscoupling are applied with respect to the integral of the phase error indication of each phase error input.
In another embodiment, the weight for each phase error input is related to the noise variance associated with the phase error indications for the phase error input being weighted.
In a further embodiment, the weight for each phase error input is determined in accordance with the formula:
where wk(i) is the weight for channel i at time k, and where ρk(i) is the inverse of the estimate of the noise variance associated with the phase error indications for channel i at time k, and where N is the total number of channels.
In another embodiment, the weight for each phase error input is determined in accordance with the formula:
where wk(i) is the weight for channel i at time k, and where ρk(i) is the inverse of the estimate of the noise variance associated with the phase error indications for channel i at time k, and where N is the total number of channels.
In a further embodiment, the weighting and applied gain are determined in accordance with the formula:
βk(i)=wk(i)√{square root over (ζ)},
where βk(i) is the gain for channel i at time k, where wk(i) is the weight for channel i at time k, and where √ζ is a common gain and is related to the ratio of the noise variance that characterizes the frequency error variations to the noise variance that characterizes the phase error indications for all channels.
In another embodiment, a read detection signal clocking system for a magnetic tape drive comprises multiple phase inputs, each indicating the phase of a sample taken with respect to one of a plurality of read channels of the magnetic tape drive, the read channels configured to provide samples of data from magnetic tape read by the magnetic tape drive; phase error logic responsive to each of the multiple phase inputs configured to determine phase error of the sample of one of the read channels and provide a phase error indication at a phase error input with respect to each of the plurality of read channels; logic responsive to the phase error inputs configured to weight and crosscouple the phase error indication of each phase error input with the phase error indication of each other phase error input, and to apply common gain; and feedback logic responsive to the crosscoupling configured to provide a sample timing phase estimate for each read channel.
In still another embodiment, a common state-space multi-channel digital sample timing phase control of multiple read channels for correlated signals, the signals correlated with respect to the timing of the signals, comprises multiple phase inputs, each indicating the phase of at least one digital sample taken with respect to one of a plurality of read channels at a current state, the read channels configured to provide digital samples of the correlated signals. Phase error logic responsive to each of the multiple phase inputs is configured to determine phase error of the sample of one of the read channels at the current state and provide a phase error indication at a phase error input with respect to each of the plurality of read channels. Logic responsive to the phase error inputs is configured to weight and crosscouple the phase error indication of each phase error input with the phase error indication of each other phase error input, and to apply common gain; and feedback logic responsive to the crosscoupling and configured to provide a sample timing phase estimate for the next state of each read channel.
For a fuller understanding of the present invention, reference should be made to the following detailed description taken in conjunction with the accompanying drawings.
This invention is described in preferred embodiments in the following description with reference to the Figures, in which like numbers represent the same or similar elements. While this invention is described in terms of the best mode for achieving this invention's objectives, it will be appreciated by those skilled in the art that variations may be accomplished in view of these teachings without deviating from the spirit or scope of the invention.
Referring to
Referring to
Timing recovery typically is based on a PLL (phase locked loop) for each channel whose purpose is to accurately estimate the timing offsets before sampling the analog signal. Many versions of PLLs exist, including the PLLs discussed in the incorporated U.S. Pat. No. 5,442,315.
The output signal for each of the channels from the head 17 of
The resultant samples for each channel are applied to a data detector 61, 62 . . . 68. One type of data detector is a maximum likelihood sequence detector which compares the incoming signals to defined particular expected signals, thus generating path metrics, and maintains a path memory of possible data sequences and selects the data sequence that has the maximum likelihood of being correct. The resultant data is output on line 71, 72 . . . 78 and the parallel data may further be combined by subsequent circuitry as needed by the multi-channel detection system of
To achieve interpolating to the right time within the received digital samples involves detecting any phase error still present in the signal after sample interpolation 51, 52 . . . 58 and correcting for it.
State-space logic 80, in one embodiment, comprises multiple phase error inputs 91, 92 . . . 98, each indicating phase error of one of the read channels. The state-space logic is configured to respond to the multiple phase error inputs and weight and crosscouple the phase error indication of each phase error input with the phase error indication of each other phase error input, and to apply common gain. The weighting and crosscoupling, in one embodiment is with respect to the integral of the phase error indication of each phase error input, as will be discussed.
In one embodiment, the weight for each phase error input is related to the noise variance for the phase error input being weighted, and may be determined by weighting logic 101, 102 . . . 108. In one example, the noise variance may comprise the noise variance of the phase error indications for the channel. Alternatively, the weighting may be based on any suitable measure of noise variance for the signal sampling phase.
Feedback logic is responsive to the crosscoupling and is configured to provide a sample timing phase estimate 111, 112 . . . 118 for each read channel.
The sample timing phase estimate adjusts the sample times 51, 52 . . . 58 to tend to correct phase error in the signal.
To model the timing process, the frequency offset is expressed as a “first-order auto-regression”
fk+1=fk+vk. (1)
Here, fk denotes the frequency or frequency offset at time k and vk denotes the realization of a zero-mean additive white Gaussian noise (AWGN) process with autocorrelation function E{vkvl}=qδkl. Physically, in the example of a tape system, {vk} includes the effects of motor-speed and tape-path induced frequency variations as well as various noise phenomena that occur during the write process.
The sampling phase at time k+1, θk+1, is obtained from the sampling phase and frequency at time k as:
θk+1=θk+fk. (2)
For the multichannel case, the frequency and phase evolutions given by (1) and (2) are used for each one of the N channels.
These random processes are represented in vector form as vk=[vk(1) . . . vk(N)]′, where the prime denotes vector transposition. The vector vk therefore denotes a vector zero-mean Gaussian random process, with autocorrelation matrix Q, that is, E{vkv′l}=Qδkl. At the output of the first integration steps on each channel, the frequency offset signals fk(1), . . . , fk(N) are obtained. The second integration steps on each channel provide the sampling phase signals θk(1), . . . θk(N). The observations of the sampling phases are perturbed by the N measurement noise processes nk(1), . . . , nk(N). These noise processes can be represented in vector form as nk=[nk(1) . . . nk(N)]′. The vector nk then represents a vector zero-mean Gaussian noise process, with autocorrelation matrix denoted by R, that is E{nkn′l}=Rδkl. It may be assumed that the noise processes nk(1), . . . , nk(N) are independent Gaussian processes. Therefore, the matrix R is a diagonal matrix with diagonal elements denoted by R=diag(r1, . . . , rN). Finally, the noisy observations of the sampling phases are denoted by yk(1), . . . , yk(N).
The optimum state estimator for the state-space model of
In general, the illustrated multichannel PLL structure involves a weighted sum of the phase errors (the yk(i)−θk/k−1(i)'s) for the proportional as well as for the integral terms of the multi-input multi-output (MIMO) loop filter. The cross-coupling gains αk(ij), βk(ij), i≠j, give the optimum coupling between channels.
Herein, the frequency offset variations in the individual channels are indicated as correlated. As mentioned above, this is captured by stating that the noise processes vk(1), vk(2), . . . , vk(N) are correlated. In accordance with the present invention, the timing process in a multichannel tape system comprises replacing the vector noise process vk by the scalar noise process vk=vk(1)= . . . =vk(N).
The actual loop gains are, in general, time dependent. Herein, a steady-state solution is developed, by employing the steady-state loop gains of the optimum estimator for the above state-space model and also taking into account that in an actual system the process noise vk will have a much lower power than the power of the measurement noise nk(1) and nk(2), which leads to the gain definitions:
where α and β, with proper upper indices, comprise the proportional and integral loop coefficients, respectively; where ri represents the variance (or the power) of the measurement noise on the channel having number i, that is, this is the variance of the noise process nk(i); and where q represents the power of the noise process vk.
The two outputs of the MIMO loop filter are integrated to yield the sampling phase estimates θk/k−1(1) and θk/k−1(2).
In the general case of an N-channel system, the solution is given by:
where ρj=1/rj. The noise variance rj associated with channel j may be regarded as the phase jitter experienced on this channel.
Finally, defining ζj=q/rj, equations (4) and (5) are expressed in the form
Referring to
Phase error inputs 91 and 92 each indicates phase error of one of the read channels, subtracting the sample timing phase estimate 111, 112 from the phase of the input signal. The state-space logic 80 is configured to respond to the multiple phase error inputs 91, 92 and weight 101, 102 and crosscouple the phase error indication of each phase error input with the phase error indication of each other phase error input, and to apply common gain 130. The weighting and crosscoupling, in the illustrated embodiment is with respect to the integral of the phase error indication of each phase error input.
Feedback logic 141, 142 is responsive to the crosscoupling and configured to provide a sample timing phase estimate 111, 112 for each read channel.
The notation βk(i)=wk(i)√{square root over (ζ)} is used for the loop coefficients of the integral terms. The noncrosscoupled loop coefficients for the proportional terms are denoted in this figure by γ1 and γ2, where the proportional parts of the loop filters are realized per channel, i.e., are not global.
In accordance with the present invention, the combination of the phase errors for the proportional terms is omitted without incurring any significant performance penalty. The combination of the phase errors is implemented for the integral terms only.
Furthermore, the set of parameters ζi=q/ri are selected following a traditional, e.g., a loop time-response based, approach which initially does not take the loop couplings into account. This means that the ζi may be chosen to be equal, ζi=ζ∀i, comprising a common gain 130. Then, coupling across channels is introduced for the integral terms only, according to equation (7), as
where ρk(i) is the inverse of the estimate of the measurement noise power for channel i at time k, i.e., an estimate of 1/ri, and N is the total number of channels.
Alternatively, the weighting coefficient is defined to omit the square-root, defining the weighting coefficient as
where wk(i) is the weight for channel i at time k, and where ρk(i) is the inverse of the estimate of the noise variance of said phase error indications for channel i at time k.
There can be various ways to compute ρk(i). One is to derive it from a comparison of the phase of the input signal to the phase locked loop of the considered channel i and the signal phase of the expected bit cell, i.e., by computing the variance of the so-obtained error signal and furthermore taking the inverse of the estimated variance.
Alternatively, ρk(i) is approximated by computing the variance of the signal provided by the timing-error detector and again taking the inverse of it.
Estimates of the common gain 130 may be made by estimating ζ of equation (8) by estimating the parameters q and r, or employ a value that the designer has found appropriate by using any appropriate technique. An example is to employ the loop coefficient value (for the integral term) of a classical PLL design, where the loop coefficients are chosen by considering the time response of the PLL such as the settling time, the amount of overshoot, etc. This value of ζ may be kept fixed and the time variability taken into account through the weighting coefficients w.
Similarly, the γ terms may be fixed: here the notation α is not used since the filters are not crosscoupled for the proportional terms. As mentioned above for ζ, these γ coefficients can be chosen in any manner that the designer sees fit, e.g., using classical time-response based PLL design.
Those of skill in the art will understand that changes may be made with respect to the methods discussed above, including changes to the ordering of the steps. Further, those of skill in the art will understand that differing specific component arrangements may be employed than those illustrated herein.
While the preferred embodiments of the present invention have been illustrated in detail, it should be apparent that modifications and adaptations to those embodiments may occur to one skilled in the art without departing from the scope of the present invention as set forth in the following claims.
Number | Name | Date | Kind |
---|---|---|---|
5986834 | Colineau et al. | Nov 1999 | A |
6549561 | Crawford | Apr 2003 | B2 |
7110485 | Bertschmann et al. | Sep 2006 | B2 |
7245677 | Pare, Jr. | Jul 2007 | B1 |
7274762 | Schenk et al. | Sep 2007 | B2 |
7433142 | Bui et al. | Oct 2008 | B2 |
20050144534 | Hutchins et al. | Jun 2005 | A1 |
20090257540 | Tan | Oct 2009 | A1 |
Number | Date | Country |
---|---|---|
0901126 | Mar 1999 | EP |
1404051 | Mar 2004 | EP |
Number | Date | Country | |
---|---|---|---|
20090097603 A1 | Apr 2009 | US |