US20080278837A1 - Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive - Google Patents

Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive Download PDF

Info

Publication number
US20080278837A1
US20080278837A1 US12/178,426 US17842608A US2008278837A1 US 20080278837 A1 US20080278837 A1 US 20080278837A1 US 17842608 A US17842608 A US 17842608A US 2008278837 A1 US2008278837 A1 US 2008278837A1
Authority
US
United States
Prior art keywords
digitized samples
data
field
determining
preamble
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US12/178,426
Inventor
Chengzhi Pan
Andrei Vityaev
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Avago Technologies International Sales Pte Ltd
Original Assignee
Chengzhi Pan
Andrei Vityaev
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chengzhi Pan, Andrei Vityaev filed Critical Chengzhi Pan
Priority to US12/178,426 priority Critical patent/US20080278837A1/en
Publication of US20080278837A1 publication Critical patent/US20080278837A1/en
Assigned to BANK OF AMERICA, N.A., AS COLLATERAL AGENT reassignment BANK OF AMERICA, N.A., AS COLLATERAL AGENT PATENT SECURITY AGREEMENT Assignors: BROADCOM CORPORATION
Assigned to AVAGO TECHNOLOGIES GENERAL IP (SINGAPORE) PTE. LTD. reassignment AVAGO TECHNOLOGIES GENERAL IP (SINGAPORE) PTE. LTD. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: BROADCOM CORPORATION
Assigned to BROADCOM CORPORATION reassignment BROADCOM CORPORATION TERMINATION AND RELEASE OF SECURITY INTEREST IN PATENTS Assignors: BANK OF AMERICA, N.A., AS COLLATERAL AGENT
Abandoned legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G11INFORMATION STORAGE
    • G11BINFORMATION STORAGE BASED ON RELATIVE MOVEMENT BETWEEN RECORD CARRIER AND TRANSDUCER
    • G11B20/00Signal processing not specific to the method of recording or reproducing; Circuits therefor
    • G11B20/10Digital recording or reproducing
    • G11B20/10009Improvement or modification of read or write signals
    • G11B20/10018Improvement or modification of read or write signals analog processing for digital recording or reproduction
    • G11B20/10027Improvement or modification of read or write signals analog processing for digital recording or reproduction adjusting the signal strength during recording or reproduction, e.g. variable gain amplifiers
    • GPHYSICS
    • G11INFORMATION STORAGE
    • G11BINFORMATION STORAGE BASED ON RELATIVE MOVEMENT BETWEEN RECORD CARRIER AND TRANSDUCER
    • G11B20/00Signal processing not specific to the method of recording or reproducing; Circuits therefor
    • G11B20/10Digital recording or reproducing
    • G11B20/10009Improvement or modification of read or write signals
    • GPHYSICS
    • G11INFORMATION STORAGE
    • G11BINFORMATION STORAGE BASED ON RELATIVE MOVEMENT BETWEEN RECORD CARRIER AND TRANSDUCER
    • G11B20/00Signal processing not specific to the method of recording or reproducing; Circuits therefor
    • G11B20/10Digital recording or reproducing
    • G11B20/10009Improvement or modification of read or write signals
    • G11B20/10037A/D conversion, D/A conversion, sampling, slicing and digital quantisation or adjusting parameters thereof
    • GPHYSICS
    • G11INFORMATION STORAGE
    • G11BINFORMATION STORAGE BASED ON RELATIVE MOVEMENT BETWEEN RECORD CARRIER AND TRANSDUCER
    • G11B20/00Signal processing not specific to the method of recording or reproducing; Circuits therefor
    • G11B20/10Digital recording or reproducing
    • G11B20/10009Improvement or modification of read or write signals
    • G11B20/10222Improvement or modification of read or write signals clock-related aspects, e.g. phase or frequency adjustment or bit synchronisation
    • GPHYSICS
    • G11INFORMATION STORAGE
    • G11BINFORMATION STORAGE BASED ON RELATIVE MOVEMENT BETWEEN RECORD CARRIER AND TRANSDUCER
    • G11B27/00Editing; Indexing; Addressing; Timing or synchronising; Monitoring; Measuring tape travel
    • G11B27/10Indexing; Addressing; Timing or synchronising; Measuring tape travel
    • G11B27/19Indexing; Addressing; Timing or synchronising; Measuring tape travel by using information detectable on the record carrier
    • G11B27/28Indexing; Addressing; Timing or synchronising; Measuring tape travel by using information detectable on the record carrier by using information signals recorded by the same method as the main recording
    • G11B27/30Indexing; Addressing; Timing or synchronising; Measuring tape travel by using information detectable on the record carrier by using information signals recorded by the same method as the main recording on the same track as the main recording
    • G11B27/3027Indexing; Addressing; Timing or synchronising; Measuring tape travel by using information detectable on the record carrier by using information signals recorded by the same method as the main recording on the same track as the main recording used signal is digitally coded

Definitions

  • the read channel When data is read from a hard disk drive, the read channel must be adequately synchronized so that it may correctly read the data stored in the hard disk drive. In addition, the receiver gain must be adjusted in order to properly read the data. The presence of noise may significantly affect the ability of the data to be properly read when inadequate gain control and timing recovery mechanisms are implemented. In such an instance, the performance of the hard disk drive may be significantly affected.
  • Various aspects of the invention provide for one or more methods and systems for acquiring gain and timing information from one or more preambles of one or more data sectors of a hard disk drive.
  • the gain and timing information is used to optimize performance of the read channel receiver as it reads data stored in the hard disk drive.
  • a method of synchronizing a read channel of a hard disk drive in reference to data that is read by the read channel, comprises first generating two or more preliminary phase estimates based on one or more preambles of one or more data sectors in the hard disk drive, second generating a frequency offset and a phase offset by applying linear regression to the two or more preliminary phase estimates, and using the frequency offset and the phase offset to perform the synchronizing.
  • a method of determining a frequency offset and phase offset associated with reading data from a data sector of a hard disk drive comprises first generating one or more consecutive samples obtained from a preamble of the data sector of the hard disk drive, second generating two or more preliminary phase estimates using the one or more consecutive samples, applying linear regression techniques using the two or more preliminary phase estimates, and computing the frequency offset and the phase offset using one or more registers.
  • a method of determining a frequency offset and a phase offset during reading of a data sector of a hard disk drive comprises determining sine and cosine projections using a number of consecutive samples obtained from a preamble of the data sector of the hard disk drive, using the sine and cosine projections to generate two or more preliminary phase estimates, and using linear regression to generate the frequency offset and the phase offset using the two or more preliminary phase estimates.
  • a method of estimating the gain or amplitude of a data stream read by a read channel receiver of a hard disk drive comprises determining one or more sine and cosine projections of a number of consecutive samples of one or more preambles of one or more data sectors in the hard disk drive.
  • a system for synchronizing a read channel of a hard disk drive comprises a first circuitry and/or software used to compute preliminary phase offset estimates using a sequence of consecutive data samples obtained by sampling a preamble of a data sector of the hard disk drive, and a second circuitry used to compute frequency and phase offsets using two or more of the preliminary phase offset estimates.
  • a system for determining gain or amplitude of a data stream that is read from a data sector of a hard disk drive comprises a first circuitry and/or software used to compute sine and cosine projections of a sequence of consecutive samples.
  • the sequence is obtained from sampling a preamble of the data sector wherein the sampling occurring at fixed time intervals.
  • FIG. 1 is a block diagram of an organizational structure of a sector of data that is read by a read channel receiver in a hard disk drive, in accordance with an embodiment of the invention.
  • FIG. 2 illustrates a block diagram of a system in a read channel receiver of a hard disk drive, used to determine gain, frequency offset, and phase offset, in reference to a data stream, in accordance with an embodiment of the invention.
  • FIG. 3 is a block diagram of a subsystem, in a read channel receiver or a hard disk drive, that computes estimates of the preliminary phase offset and gain, for one or more sets of consecutive preamble samples, in accordance with an embodiment of the invention.
  • FIG. 4 is a block diagram of a subsystem of a read channel receiver that computes frequency and phase offsets used to synchronize the read channel receiver of a magnetic hard disk drive, in accordance with an embodiment of the invention.
  • FIG. 5 is an operational flow diagram illustrating a method of generating frequency offset, phase offset, and gain (or amplitude) of a data stream read by a read channel receiver of a hard disk drive, in accordance with an embodiment of the invention.
  • aspects of the invention provide for one or more methods and systems for acquiring gain and timing information from one or more preambles of one or more data sectors of a hard disk drive.
  • the gain and timing information is used to optimize performance of the read channel receiver as it reads data stored in the hard disk drive.
  • aspects of the one or more methods and systems improves the performance of a hard disk drive. The performance is improved by implementing an improved system and method of acquiring timing and gain information/parameters by the read channel receiver of the hard disk drive.
  • aspects of the invention employ one or more systems and methods of determining the frequency and phase offsets of a hard disk drive read channel, used for timing and gain acquisition in the read channel receiver.
  • the one or more systems and methods utilize the data provided by each preamble of a data sector in a hard disk drive.
  • the frequency and phase offsets are used to more accurately synchronize the read channel to the data stream being read.
  • aspects of the invention provide an estimate of the amplitude or gain of the read data stream.
  • use of various aspects of the invention allow a reduction in the number of bytes used in one or more preambles of one or more sectors of data, thereby improving performance and the storage capacity of the hard disk drive.
  • FIG. 1 is a block diagram of an organizational structure of a sector of data that is read by a read channel receiver in a hard disk drive, in accordance with an embodiment of the invention.
  • the sector comprises one or more sections.
  • the one or more sections comprises a preamble field 104 , a first sync field 108 , a first data stream field 112 , a second sync field 116 , a second data stream field 120 , and a post-amble field 124 .
  • the preamble field 104 may comprise one or more bytes that are used for timing recovery or synchronization of the following one or more data streams.
  • the one or more data streams are read by the read channel receiver.
  • the preamble field 104 comprises 15 bytes.
  • the preamble field 104 comprises less than 15 bytes. In a representative embodiment of the invention, the preamble field 104 , used to synchronize to the data stream, comprises 6 bytes. Aspects of the invention provide a reduction in the number of bytes utilized in each of the preambles of each of the preamble fields 104 of each sector so as to increase the overall data storage capacity of the magnetic hard disk drive.
  • a byte of the preamble comprises the sequence of bits “11001100”, corresponding to a sinewave. The byte represents two periods of the sinewave.
  • the first sync field 108 contains one or more bytes that notifies the read channel receiver or other circuitry that data immediately follows.
  • the first data stream field 112 contains a first stream of data bytes that is read by the magnetic hard disk drive.
  • the first data stream field 112 is followed by a second sync field 116 .
  • the second sync field 116 is used to notify or indicate the start of another section of data.
  • the second data stream field 120 contains a second stream of data bytes.
  • the second data stream field 120 is followed by a post-amble field 124 used to indicate that the end of the one or more data streams has been reached for the sector.
  • FIG. 2 illustrates a block diagram of a system in a read channel receiver of a hard disk drive, used to determine gain, frequency offset, and phase offset, in reference to a data stream, in accordance with an embodiment of the invention.
  • the system comprises a processing/conditioning subsystem 202 , phase and gain computational subsystem 204 , and a frequency and phase offset computational subsystem 208 .
  • One or more preamble samples are provided by a read magnetic head (or transducer head) of the read channel receiver, after being digitized. Typically, the read magnetic head outputs a sinusoidal signal when the magnetic disk media is read.
  • the sinusoidal signal is output to the processing/conditioning subsystem 202 , which may comprise the following sequence of components; a preamplifier, an analog processing and/or conditioning circuitry, an analog to digital converter, and a digital processing and/or conditioning circuitry.
  • the processing/conditioning subsystem 202 may receive the sinusoidal signal output by the read magnetic head so as to process, condition, and sample the sinusoidal signal into digitized samples for input into the phase and gain computational subsystem 204 .
  • the sequence of digitized samples may be represented by the following equation:
  • the variable, a represents amplitude
  • f 0 represents an a priori known center frequency
  • k represents integer variable that indexes the one or more digital samples provided to the phase and gain computational subsystem 204
  • n 0 represents noise.
  • f 0 1/N where N is an integer divisible by 4.
  • aspects of the invention utilize a sequence of digitized samples, represented by equation (1), that are input into the phase and gain computational subsystem 204 .
  • N 4, or the number of samples in one period equals four digitized samples, (e.g., the bits 1100 ) corresponding to one period of the sinewave previously mentioned in reference to FIG. 2 .
  • the number of digitized samples utilized by the phase and gain computational subsystem 204 comprises 8 samples, which is the number of samples contained in 2 periods of the previously mentioned sinusoidal signal.
  • aspects of the invention do not limit the number of bits used to generate the digitized samples, in a representative embodiment, 6 bits are used to generate the digitized samples. As a consequence, the one or more values provided by each of the digitized samples has 6 bits of resolution, allowing its values to vary from ⁇ 32 to +31.
  • FIG. 3 is a block diagram of a subsystem, in a read channel receiver or a hard disk drive, that computes estimates of the preliminary phase offset and gain, for one or more sets of consecutive preamble samples, in accordance with an embodiment of the invention.
  • the phase and gain computational subsystem shown in FIG. 3 , computes the preliminary phase offset and gain estimates by determining the sine and cosine projections using one or more digitized preamble samples.
  • the phase and gain computational subsystem comprises a sine component (or projection) processing hardware/software 304 , a cosine component processing hardware/software 308 , a memory 312 , and a gain computation hardware/software 316 .
  • the sine component processing hardware/software 304 and the cosine component processing hardware/software 308 each provide a projection of the exemplary eight samples, provided by y(k), either in the sine direction or the cosine direction.
  • the sine component processing hardware/software 304 and the cosine component processing hardware/software 308 may each comprise a processor, such as a central processing unit, used to execute commands and/or perform digital computations.
  • the sine component processing hardware/software 304 and the cosine component processing hardware/software 308 may each comprise any logic circuitry and/or set of executable instructions recognized by the processor, which is capable of computing its respective sine or cosine components (or projections).
  • the set of executable instructions may be stored within the memory 312 .
  • the logic circuitry and/or software computes the preliminary phase offset estimate by applying an equation (to be introduced as equation (6), later in this document), that utilizes the outputs provided by the sine component processing hardware/software 304 and cosine component processing hardware/software 308 .
  • the memory 312 serves as a look up table in which the sine and cosine projections, computed by the sine component processing hardware/software 304 and cosine component processing hardware/software 308 , may be appropriately mapped into a preliminary phase offset estimate, to be used by the frequency and phase offset computational subsystem, later shown in FIG. 4 .
  • the memory 312 may comprise a random operating memory or a non-volatile memory.
  • the memory 312 may comprise any circuitry capable of storing data in which one or more values may be appropriately indexed and/or mapped.
  • the gain computation hardware/software 316 may comprise any logic circuitry and/or set of executable instructions (i.e., software) recognized by the processor, which is capable of computing the gain estimate.
  • the logic circuitry and/or software computes the gain estimate by applying an equation (to be introduced as equation (6), later in this document), that utilizes the outputs provided by the sine component processing hardware/software 304 and cosine component processing hardware/software 308 .
  • the sine projection may be represented by the following set of equations.
  • Equation (2) is obtained by substituting equation (1) into P sin (k), for y(k). Equation (1) may be easily represented by the following equation:
  • y ⁇ ( k ) a ⁇ sin ⁇ ⁇ 2 ⁇ ⁇ ⁇ ⁇ ( k N + 1 8 ) ⁇ cos ⁇ ⁇ 2 ⁇ ⁇ ⁇ ⁇ ( ⁇ ⁇ ⁇ fk + ⁇ 0 - 1 8 ) + a ⁇ cos ⁇ ⁇ 2 ⁇ ⁇ ⁇ ⁇ ( k N + 1 8 ) ⁇ sin ⁇ ⁇ 2 ⁇ ⁇ ⁇ ( ⁇ ⁇ ⁇ fk + ⁇ 0 - 1 8 ) + n 0 ⁇ ( k )
  • n 1 , n 2 , n 3 denote various noise components, resulting from white noise or the offset in frequency.
  • the cosine projection may be represented by the following set of equations:
  • equation (3) is obtained by substituting equation (1) into P cos (k), for y(k).
  • the variables n 5 , n 6 , n 7 denote various noise components, resulting from white noise or the offset in frequency.
  • the phase estimate may be determined for the kth sample, where the
  • phase at the kth sample is given by
  • ⁇ ⁇ ( k ) ⁇ ⁇ ⁇ fk + ⁇ 0 - 1 8 ,
  • phase estimate may be found by taking the arctangent of Q(k), and may be represented by the following equation:
  • n Q (k) and n ⁇ (k) represent noise in each of equations (4) and (5).
  • the amplitude or gain is found by calculating the square root of the sum of the squares of the sine and cosine projections.
  • FIG. 4 is a block diagram of a subsystem of a read channel receiver that computes frequency and phase offsets used to synchronize the read channel receiver of a magnetic hard disk drive, in accordance with an embodiment of the invention.
  • the frequency and phase computational subsystem shown in FIG. 4 , comprises a first accumulator 404 , a second accumulator 412 , a counter 424 for the integer variable, k, which indexes the one or more preliminary phase offset estimates that are input, a first multiplier 428 , a second multiplier 432 , a third multiplier 436 , a fourth multiplier 440 , a fifth multiplier 444 , a phase offset adder 448 , and a frequency offset adder 452 .
  • the first accumulator 404 comprises a first accumulator adder 412 and a first register 408 while the second accumulator 412 comprises a second accumulator adder 420 and a second register 416 .
  • the first accumulator 404 accumulates (or provides a cumulative sum of) preliminary phase offset estimates provided by the phase and gain computational subsystem.
  • the second accumulator 412 accumulates successive terms and functionally operates similar to that of the first accumulator 404 .
  • the two factors that generate the successive terms comprise a count of the counter 424 and its associated preliminary phase offset estimate.
  • the first adder 412 within the first accumulator 404 adds the preliminary phase offset estimate to the accumulated value stored in the first register 408 .
  • the second adder 416 within the second accumulator 412 adds the terms (associated with the product of the preliminary phase offset estimate and the count) to the accumulated value stored in the second register 416 .
  • the values stored in the accumulators 404 , 412 are used to compute the frequency and phase offsets.
  • the frequency offset is computed by summing two products, a first product and a second product.
  • the first product is obtained by multiplying the value in the first accumulator 404 with ⁇ (L) while the second product is obtained by multiplying the value in the second accumulator 412 with ⁇ (L).
  • the phase offset is computed by summing two products (a third product and a fourth product) and a constant.
  • the third product is obtained by multiplying the value in the first accumulator 404 with ⁇ (L) while the fourth product is obtained by multiplying the value in the second accumulator 412 with ⁇ (L).
  • the third product, fourth product, and a constant (1 ⁇ 8 in this representative embodiment) are added to generate the phase offset.
  • the aforementioned counter, adders, multipliers, and registers may be implemented using digital logic circuitry.
  • aspects of the invention comprise utilizing a linear regression technique for determining the phase and frequency offsets of a read channel in a read channel receiver.
  • the technique is employed in implementing the embodiment illustrated in FIG. 4 .
  • the function F( ⁇ , ⁇ ) is a quadratic function in ⁇ and ⁇ so to find the minimum it is sufficient to take the partial derivatives and set them to zero.
  • the partial derivatives of F( ⁇ , ⁇ ) with respect to ⁇ and ⁇ are as follows;
  • the vector x comprises sampling times.
  • the sampling times comprises fixed sampling intervals.
  • L corresponds to the number of samples used in the linear regression technique (i.e., number of preliminary phase estimates used to compute a frequency offset and a phase offset).
  • the inverse of the matrix shown in equation (8) only depends on the variable, L, and does not depend on y, and as a result, may be determined in advance.
  • the terms within the invertible 2 ⁇ 2 matrix of equation (8) may be represented by the following equation:
  • Equation (9) the terms in the 2 ⁇ 2 matrix of equation (9), ⁇ (L), ⁇ (L), ⁇ (L), and ⁇ (L), are all functions of the predetermined variable L. These terms may be obtained by inverting the 2 ⁇ 2 matrix shown in equation (8).
  • equation (9) may be implemented using two registers, as shown in equations (10a) and (10b).
  • the two registers, R xy and R y are substituted into the matrix equation (9).
  • the matrix equation (9) may be expanded into the following two equations (11a) and (11b):
  • ⁇ ⁇ ( k ) ⁇ ⁇ ⁇ fk + ⁇ 0 - 1 8 ,
  • equations (10a) and (10b) Using equations (10a) and (10b), and applying the sequence of preliminary phase offset estimates ⁇ circumflex over ( ⁇ ) ⁇ (1), ⁇ circumflex over ( ⁇ ) ⁇ (2), . . . , ⁇ circumflex over ( ⁇ ) ⁇ (L) provided by the phase and gain computational circuitry, the following equations are obtained, in which ⁇ circumflex over (f) ⁇ (L) corresponds to the frequency offset while ⁇ circumflex over ( ⁇ ) ⁇ 0 (L) corresponds to the phase offset.
  • R xy ( k ) R xy ( k ⁇ 1)+ k ⁇ circumflex over ( ⁇ ) ⁇ ( k ) (12a)
  • R y ( k ) R y ( k ⁇ 1)+ ⁇ circumflex over ( ⁇ ) ⁇ ( k ) (12b)
  • the embodiment provides an implementation that adjusts to a maximum number of digital samples, L, as one or more preambles of one or more disk sectors are input into the read channel of a magnetic hard disk drive.
  • Equations (12a) and (12b) represent the implementation of the first portion of FIG. 4 while equations (13a) and (13b) represent the implementation of the second portion of FIG. 4 .
  • FIG. 5 is an operational flow diagram illustrating a method of generating frequency offset, phase offset, and gain (or amplitude) of a data stream read from data stored in a data sector, by a read channel receiver of a hard disk drive, in accordance with an embodiment of the invention.
  • a preamble is read from the data sector of one or more data sectors in a hard disk drive. The one or more bytes contained within the preamble is used to generate the frequency offset, phase offset, and gain of the data stream.
  • the number of bytes used in the preamble is less than the typical number of bytes used—15 bytes
  • a sequence of digitized samples is provided to the phase and gain computational subsystem.
  • a processing/conditioning circuitry performs the processing and conditioning.
  • the sequence of digitized samples comprises the latest fixed number of digitized samples.
  • the last fixed number of digitized samples is provided to the phase and gain computational subsystem In a representative embodiment, this fixed number, x, comprises the number of samples used to determine the sine and cosine projections using the phase and gain computational subsystem.
  • the value for the fixed number, x is 8. This fixed number may be programmed by a user, prior to performing the computations provided by the phase and gain computational subsystem.
  • the phase and gain computational subsystem determines the sine and cosine projections (or sine and cosine components) based on the sequence of digitized samples provided Equations (2) and (3) are used to compute these projections.
  • Q(k) the quotient is computed for P cos(k)/P sin(k).
  • a look-up table (implemented in memory, for example) is used to determine the arctangent of Q(k), ⁇ ( ⁇ ).
  • the gain or amplitude is output by the phase and gain computational subsystem, for use by the hard disk drive.
  • the phase and gain computational subsystem outputs ⁇ ( ⁇ ) (preliminary phase offset estimates) into the frequency and phase computational subsystem, in which, the processing described in reference to FIG. 4 , is accomplished.
  • the frequency and phase computational subsystem outputs the phase and frequency offsets associated with the sequence of samples provided to the frequency and phase computational subsystem.
  • an assessment is made whether the last sample has been generated by the preamble. If the end of the preamble of the data sector has been sampled, then the process ends for this particular data sector.
  • a new data sector may be read, if necessary, which contains another preamble. If a new data sector is read, then the process commences with another cycle, starting at step 504 . Otherwise, the process resumes with step 532 , at which, the next digitized sample is introduced at the next sample clock, using the remaining un-sampled preamble data. Then, the process repeats itself at step 504 , when a new sequence of digitized samples is introduced to the phase and gain computational subsystem. In this exemplary embodiment, 7 of the previous 8 samples are used in the next sequence of samples. The 7 previously used samples are shifted by one sample clock, so as to accommodate the latest new sample.
  • additional performance improvements may be obtained by utilizing knowledge of noise statistics and applying this to the noise variables previously used in the preceding equations.
  • the performance improvements rely on an analysis of the noise variables n 1 , n 2 , n 3 , n 5 , n 6 , and n 7 used in equations (2) and (3).
  • Bayesian Linear Regression may be employed when a priori knowledge of the frequency offset distribution is used.
  • Colored Noise Linear Regression may be employed when noise correlation is taken into account while performing a linear regression.
  • d 1 (k) and d 2 (k) may be calculated as:
  • n 1 ⁇ ( k ) + n 2 ⁇ ( k ) - ⁇ ⁇ ⁇ fa ⁇ ⁇ sin ⁇ ⁇ 2 ⁇ ⁇ ⁇ ( ⁇ ⁇ ⁇ fk + ⁇ 0 - 1 8 ) + 2 ⁇ ⁇ ⁇ ⁇ fab ⁇ ⁇ sin ⁇ ⁇ 2 ⁇ ⁇ ⁇ ⁇ ( ( 2 N + ⁇ ⁇ ⁇ f ) ⁇ k + ⁇ d + ⁇ 0 - 1 8 ) ,
  • n 5 ⁇ ( k ) + n 6 ⁇ ( k ) - ⁇ ⁇ ⁇ fa ⁇ ⁇ cos ⁇ ⁇ 2 ⁇ ⁇ ⁇ ( ⁇ ⁇ ⁇ fk + ⁇ 0 - 1 8 ) + 2 ⁇ ⁇ ⁇ ⁇ fab ⁇ ⁇ cos ⁇ ⁇ 2 ⁇ ⁇ ⁇ ( ( 2 N + ⁇ ⁇ ⁇ f ) ⁇ k + ⁇ d + ⁇ 0 - 1 8 )
  • n 0 is AWGN (additive white Gaussian noise)
  • n 3 is zero-mean Gaussian.
  • the correlation matrix of n 3 is as follows:
  • r n 3 (k,k ⁇ 1) is constant w.r.t. k.
  • r n 3 ⁇ ( k , k - l ) ⁇ ⁇ n 0 2 ⁇ 2 MN ⁇ ( 1 - ⁇ l ⁇ MN ) , ⁇ l ⁇ ⁇ MN - 1 0 , ⁇ l ⁇ ⁇ MN
  • the correlation function takes a similar form.
  • the variable, L comprises the total number of samples provided by the preamble of the data sector.
  • the noise power density spectrum may be obtained by taking the fast Fourier transform of the correlation function as follows:
  • N 3 ( k ⁇ ) FFT ( Fr n 3 ( l ))
  • aspects of the invention may incorporate values for N and M that provide embodiments having significantly reduced projection noise.
  • N 3 ⁇ ( 0 ) ⁇ l ⁇ Fr n 3 ⁇ ( l ) ⁇ ⁇ n 0 2
  • the correlation matrix of n 7 in formula (3) is similar to that of n 3 :
  • the noise properties may be derived as follows:
  • r n ⁇ ⁇ ( k , k - l ) ⁇ 2 MN ⁇ ⁇ 0 2 ⁇ ( 1 - ⁇ l ⁇ MN ) ⁇ cos ⁇ ⁇ 2 ⁇ ⁇ ⁇ ⁇ ⁇ fl ⁇ ⁇ 2 MN ⁇ ⁇ 0 2 ⁇ ( 1 - ⁇ l ⁇ MN ) ( ⁇ ⁇ ⁇ fl ⁇ ⁇ 1 )
  • ⁇ 0 ⁇ n 0 2 ⁇ ⁇ ⁇ ⁇ a
  • An estimate of the noise level may be based on the following, such that autocorrelation of ⁇ tilde over ( ⁇ ) ⁇ (k) is zero when the delay is no larger than M*N
  • two alternative methods may be used for linear regression: linear regression with colored noise, and Bayesian linear regression.
  • k 0, . . . , L ⁇ 1 ⁇
  • [ ⁇ ⁇ ⁇ f ⁇ 0 ] [ k T ⁇ R n ⁇ - 1 ⁇ k k T ⁇ R n ⁇ - 1 ⁇ i n ⁇ 1 i n ⁇ 1 T ⁇ R n ⁇ - 1 ⁇ k i n ⁇ 1 T ⁇ R n ⁇ - 1 ⁇ i n ⁇ 1 ] - 1 ⁇ [ k 7 ⁇ R n ⁇ - 1 ⁇ ⁇ ⁇ i n ⁇ 1 T ⁇ R n ⁇ - 1 ⁇ ⁇ ⁇ ] + [ 0 1 / 8 ]
  • This formula replaces formulas (10), (11), (12), (13), and (14), previously described.
  • Bayesian linear regression When prior information of ⁇ f and ⁇ ′ 0 is available, or the information has been adaptively learned real time, Bayesian linear regression may be employed.
  • the likelihood function of the preliminary phase offset estimate is given by:
  • the preceding equation uses a Gaussian distribution for ⁇ f and an even distribution results for ⁇ ′ 0 ;
  • Bayesian estimation of ⁇ f is defined as:
  • ⁇ 0 2 In the event that ⁇ 0 2 is stable in different sections, one or more old a ⁇ 0 2 values may be used to calculate f in advance. In the event f varies in different sections, ⁇ 0 2 may be estimated in real-time and the formula may require a small adjustment yielding:
  • f 1 , f 2 , f 3 , f 4 , f 5 are constants.

Abstract

One or more methods and systems are presented for performing gain and timing acquisition of data read from one or more data sectors of a hard disk drive. The gain and timing information may be used to synchronize a read channel receiver of the hard disk drive to the data read from the data sector(s). In a representative embodiment, the one or more methods comprises determining one or more frequency and phase offsets using one or more preambles of one or more data sectors residing in a hard disk drive. In a representative embodiment, the gain of the one or more preambles may be determined. In another representative embodiment, the one or more systems used to determine one or more phase and frequency offsets, and amplitudes, comprise hardware and/or software capable of utilizing a sequence of consecutive samples provided by one or more preambles of one or more data sectors.

Description

    RELATED APPLICATIONS/INCORPORATION BY REFERENCE
  • [Not Applicable]
  • FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT
  • [Not Applicable]
  • MICROFICHE/COPYRIGHT REFERENCE
  • [Not Applicable]
  • BACKGROUND OF THE INVENTION
  • When data is read from a hard disk drive, the read channel must be adequately synchronized so that it may correctly read the data stored in the hard disk drive. In addition, the receiver gain must be adjusted in order to properly read the data. The presence of noise may significantly affect the ability of the data to be properly read when inadequate gain control and timing recovery mechanisms are implemented. In such an instance, the performance of the hard disk drive may be significantly affected.
  • Further limitations and disadvantages of conventional and traditional approaches will become apparent to one of skill in the art, through comparison of such systems with some aspects of the present invention as set forth in the remainder of the present application with reference to the drawings.
  • BRIEF SUMMARY OF THE INVENTION
  • Various aspects of the invention provide for one or more methods and systems for acquiring gain and timing information from one or more preambles of one or more data sectors of a hard disk drive. The gain and timing information is used to optimize performance of the read channel receiver as it reads data stored in the hard disk drive.
  • In a representative embodiment, a method of synchronizing a read channel of a hard disk drive, in reference to data that is read by the read channel, comprises first generating two or more preliminary phase estimates based on one or more preambles of one or more data sectors in the hard disk drive, second generating a frequency offset and a phase offset by applying linear regression to the two or more preliminary phase estimates, and using the frequency offset and the phase offset to perform the synchronizing.
  • In a representative embodiment, a method of determining a frequency offset and phase offset associated with reading data from a data sector of a hard disk drive comprises first generating one or more consecutive samples obtained from a preamble of the data sector of the hard disk drive, second generating two or more preliminary phase estimates using the one or more consecutive samples, applying linear regression techniques using the two or more preliminary phase estimates, and computing the frequency offset and the phase offset using one or more registers.
  • In a representative embodiment, a method of determining a frequency offset and a phase offset during reading of a data sector of a hard disk drive comprises determining sine and cosine projections using a number of consecutive samples obtained from a preamble of the data sector of the hard disk drive, using the sine and cosine projections to generate two or more preliminary phase estimates, and using linear regression to generate the frequency offset and the phase offset using the two or more preliminary phase estimates.
  • In a representative embodiment, a method of estimating the gain or amplitude of a data stream read by a read channel receiver of a hard disk drive comprises determining one or more sine and cosine projections of a number of consecutive samples of one or more preambles of one or more data sectors in the hard disk drive.
  • In a representative embodiment, a system for synchronizing a read channel of a hard disk drive comprises a first circuitry and/or software used to compute preliminary phase offset estimates using a sequence of consecutive data samples obtained by sampling a preamble of a data sector of the hard disk drive, and a second circuitry used to compute frequency and phase offsets using two or more of the preliminary phase offset estimates.
  • In a representative embodiment, a system for determining gain or amplitude of a data stream that is read from a data sector of a hard disk drive comprises a first circuitry and/or software used to compute sine and cosine projections of a sequence of consecutive samples. In this representative embodiment, the sequence is obtained from sampling a preamble of the data sector wherein the sampling occurring at fixed time intervals.
  • These and other advantages, aspects, and novel features of the present invention, as well as details of illustrated embodiments thereof will be more fully understood from the following description and drawings.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 is a block diagram of an organizational structure of a sector of data that is read by a read channel receiver in a hard disk drive, in accordance with an embodiment of the invention.
  • FIG. 2 illustrates a block diagram of a system in a read channel receiver of a hard disk drive, used to determine gain, frequency offset, and phase offset, in reference to a data stream, in accordance with an embodiment of the invention.
  • FIG. 3 is a block diagram of a subsystem, in a read channel receiver or a hard disk drive, that computes estimates of the preliminary phase offset and gain, for one or more sets of consecutive preamble samples, in accordance with an embodiment of the invention.
  • FIG. 4 is a block diagram of a subsystem of a read channel receiver that computes frequency and phase offsets used to synchronize the read channel receiver of a magnetic hard disk drive, in accordance with an embodiment of the invention.
  • FIG. 5 is an operational flow diagram illustrating a method of generating frequency offset, phase offset, and gain (or amplitude) of a data stream read by a read channel receiver of a hard disk drive, in accordance with an embodiment of the invention.
  • DETAILED DESCRIPTION OF THE INVENTION
  • Aspects of the invention provide for one or more methods and systems for acquiring gain and timing information from one or more preambles of one or more data sectors of a hard disk drive. The gain and timing information is used to optimize performance of the read channel receiver as it reads data stored in the hard disk drive. Aspects of the one or more methods and systems improves the performance of a hard disk drive. The performance is improved by implementing an improved system and method of acquiring timing and gain information/parameters by the read channel receiver of the hard disk drive. Aspects of the invention employ one or more systems and methods of determining the frequency and phase offsets of a hard disk drive read channel, used for timing and gain acquisition in the read channel receiver. The one or more systems and methods utilize the data provided by each preamble of a data sector in a hard disk drive. The frequency and phase offsets are used to more accurately synchronize the read channel to the data stream being read. In addition, aspects of the invention provide an estimate of the amplitude or gain of the read data stream. In a representative embodiment, use of various aspects of the invention allow a reduction in the number of bytes used in one or more preambles of one or more sectors of data, thereby improving performance and the storage capacity of the hard disk drive.
  • FIG. 1 is a block diagram of an organizational structure of a sector of data that is read by a read channel receiver in a hard disk drive, in accordance with an embodiment of the invention. As illustrated, the sector comprises one or more sections. The one or more sections comprises a preamble field 104, a first sync field 108, a first data stream field 112, a second sync field 116, a second data stream field 120, and a post-amble field 124. The preamble field 104 may comprise one or more bytes that are used for timing recovery or synchronization of the following one or more data streams. The one or more data streams are read by the read channel receiver. Typically, the preamble field 104 comprises 15 bytes. In a representative embodiment of the invention, the preamble field 104 comprises less than 15 bytes. In a representative embodiment of the invention, the preamble field 104, used to synchronize to the data stream, comprises 6 bytes. Aspects of the invention provide a reduction in the number of bytes utilized in each of the preambles of each of the preamble fields 104 of each sector so as to increase the overall data storage capacity of the magnetic hard disk drive. In a typical implementation, a byte of the preamble comprises the sequence of bits “11001100”, corresponding to a sinewave. The byte represents two periods of the sinewave. The first sync field 108 contains one or more bytes that notifies the read channel receiver or other circuitry that data immediately follows. The first data stream field 112 contains a first stream of data bytes that is read by the magnetic hard disk drive. In this representative embodiment, the first data stream field 112 is followed by a second sync field 116. The second sync field 116 is used to notify or indicate the start of another section of data. The second data stream field 120 contains a second stream of data bytes. The second data stream field 120 is followed by a post-amble field 124 used to indicate that the end of the one or more data streams has been reached for the sector.
  • FIG. 2 illustrates a block diagram of a system in a read channel receiver of a hard disk drive, used to determine gain, frequency offset, and phase offset, in reference to a data stream, in accordance with an embodiment of the invention. The system comprises a processing/conditioning subsystem 202, phase and gain computational subsystem 204, and a frequency and phase offset computational subsystem 208. One or more preamble samples are provided by a read magnetic head (or transducer head) of the read channel receiver, after being digitized. Typically, the read magnetic head outputs a sinusoidal signal when the magnetic disk media is read. The sinusoidal signal is output to the processing/conditioning subsystem 202, which may comprise the following sequence of components; a preamplifier, an analog processing and/or conditioning circuitry, an analog to digital converter, and a digital processing and/or conditioning circuitry. The processing/conditioning subsystem 202 may receive the sinusoidal signal output by the read magnetic head so as to process, condition, and sample the sinusoidal signal into digitized samples for input into the phase and gain computational subsystem 204. The sequence of digitized samples may be represented by the following equation:

  • y(k)=a·sin 2π[(f 0 +Δf)k+θ 0 ]+n 0(k),|Δf 0|>>  (1)
  • In equation (1), the variable, a, represents amplitude, f0 represents an a priori known center frequency, k represents integer variable that indexes the one or more digital samples provided to the phase and gain computational subsystem 204, while n0 represents noise. We assume that f0=1/N where N is an integer divisible by 4. Aspects of the invention provide for at least a system and a method of estimating the phase offset, θ0 and the frequency offset Δf, by way of determining {circumflex over (θ)}0, the estimate of the phase offset, and Δ{acute over (f)}, the estimate of the frequency offset.
  • Aspects of the invention utilize a sequence of digitized samples, represented by equation (1), that are input into the phase and gain computational subsystem 204. In a representative embodiment, N=4, or the number of samples in one period equals four digitized samples, (e.g., the bits 1100) corresponding to one period of the sinewave previously mentioned in reference to FIG. 2. In a representative embodiment, the number of digitized samples utilized by the phase and gain computational subsystem 204 comprises 8 samples, which is the number of samples contained in 2 periods of the previously mentioned sinusoidal signal. Hence, in this representative embodiment, the number of digitized samples utilized in generating a phase and/or gain estimate comprises M*N=8 digitized samples. Although aspects of the invention do not limit the number of bits used to generate the digitized samples, in a representative embodiment, 6 bits are used to generate the digitized samples. As a consequence, the one or more values provided by each of the digitized samples has 6 bits of resolution, allowing its values to vary from −32 to +31.
  • FIG. 3 is a block diagram of a subsystem, in a read channel receiver or a hard disk drive, that computes estimates of the preliminary phase offset and gain, for one or more sets of consecutive preamble samples, in accordance with an embodiment of the invention. The phase and gain computational subsystem, shown in FIG. 3, computes the preliminary phase offset and gain estimates by determining the sine and cosine projections using one or more digitized preamble samples. The phase and gain computational subsystem comprises a sine component (or projection) processing hardware/software 304, a cosine component processing hardware/software 308, a memory 312, and a gain computation hardware/software 316. The sine component processing hardware/software 304 and the cosine component processing hardware/software 308 each provide a projection of the exemplary eight samples, provided by y(k), either in the sine direction or the cosine direction. The sine component processing hardware/software 304 and the cosine component processing hardware/software 308 may each comprise a processor, such as a central processing unit, used to execute commands and/or perform digital computations. The sine component processing hardware/software 304 and the cosine component processing hardware/software 308 may each comprise any logic circuitry and/or set of executable instructions recognized by the processor, which is capable of computing its respective sine or cosine components (or projections). The set of executable instructions may be stored within the memory 312. In a representative embodiment, the logic circuitry and/or software computes the preliminary phase offset estimate by applying an equation (to be introduced as equation (6), later in this document), that utilizes the outputs provided by the sine component processing hardware/software 304 and cosine component processing hardware/software 308. In a representative embodiment, the memory 312 serves as a look up table in which the sine and cosine projections, computed by the sine component processing hardware/software 304 and cosine component processing hardware/software 308, may be appropriately mapped into a preliminary phase offset estimate, to be used by the frequency and phase offset computational subsystem, later shown in FIG. 4. The memory 312 may comprise a random operating memory or a non-volatile memory. It is contemplated that the memory 312 may comprise any circuitry capable of storing data in which one or more values may be appropriately indexed and/or mapped. The gain computation hardware/software 316 may comprise any logic circuitry and/or set of executable instructions (i.e., software) recognized by the processor, which is capable of computing the gain estimate. In a representative embodiment, the logic circuitry and/or software computes the gain estimate by applying an equation (to be introduced as equation (6), later in this document), that utilizes the outputs provided by the sine component processing hardware/software 304 and cosine component processing hardware/software 308. In a representative embodiment, the sine projection may be represented by the following set of equations.
  • P sin ( k ) 2 MN i = - MN / 2 + 1 MN / 2 y ( k + i ) · sin 2 π ( k + i N + 1 8 ) = 2 MN i = - MN / 2 + 1 MN / 2 a · sin 2 2 π ( k + i N + 1 8 ) cos 2 π ( Δ f ( k + i ) + θ 0 - 1 8 ) + 2 MN i = - MN / 2 + 1 MN / 2 a · cos 2 π ( k + i N + 1 8 ) sin 2 π ( k + i N + 1 8 ) sin 2 π ( Δ f ( k + i ) + θ 0 - 1 8 ) + 2 MN i = - MN / 2 + 1 MN / 2 n 0 ( k + i ) · sin 2 π ( k + i N + 1 8 ) = { 2 a MN [ i = - MN / 2 + 1 MN / 2 sin 2 2 π ( k + i N + 1 8 ) ] cos 2 π ( Δ fk + θ 0 - 1 8 ) + n 1 ( k ) } + { 2 a MN [ i = - MN / 2 + 1 MN / 2 cos 2 π ( k + i N + 1 8 ) sin 2 π ( k + i N + 1 8 ) ] sin 2 π ( Δ fk + θ 0 - 1 8 ) + n 2 ( k ) } + 2 MN i = - MN / 2 + 1 MN / 2 n 0 ( k + i ) · sin 2 π ( k + i N + 1 8 ) = a cos 2 π ( Δ fk + θ 0 - 1 8 ) + n 1 ( k ) + n 2 ( k ) + n 3 ( k ) ( 2 )
  • Equation (2) is obtained by substituting equation (1) into Psin(k), for y(k). Equation (1) may be easily represented by the following equation:
  • y ( k ) = a · sin 2 π ( k N + 1 8 ) cos 2 π ( Δ fk + θ 0 - 1 8 ) + a · cos 2 π ( k N + 1 8 ) sin 2 π ( Δ fk + θ 0 - 1 8 ) + n 0 ( k )
  • For each k≧MN/2-1, M*N consecutive samples may be used to generate Psin(k), using {y(k−MN/2+1), . . . , y(k+MN/2)}. The variables n1, n2, n3 denote various noise components, resulting from white noise or the offset in frequency. Similarly, the cosine projection may be represented by the following set of equations:
  • P cos ( k ) 2 MN i = - MN / 2 + 1 MN / 2 y ( k + i ) · cos 2 π ( k + i N + 1 8 ) = 2 MN i = - MN / 2 + 1 MN / 2 a · sin 2 π ( k + i N + 1 8 ) cos 2 π ( k + i N + 1 8 ) cos 2 π ( Δ f ( k + i ) + θ 0 - 1 8 ) + 2 MN i = - MN / 2 + 1 MN / 2 a · cos 2 2 π ( k + i N + 1 8 ) sin 2 π ( Δ f ( k + i ) + θ 0 - 1 8 ) + 2 MN i = - MN / 2 + 1 MN / 2 n 0 ( k + i ) · cos 2 π ( k + i N + 1 8 ) = { 2 a MN [ i = - MN / 2 + 1 MN / 2 sin 2 π ( k + i N + 1 8 ) cos 2 π ( k + i N + 1 8 ) ] cos 2 π ( Δ fk + θ 0 - 1 8 ) + n 5 ( k ) } + { 2 a MN [ i = - MN / 2 + 1 MN / 2 cos 2 2 π ( k + i N + 1 8 ) ] sin 2 π ( Δ fk + θ 0 - 1 8 ) + n 6 ( k ) } + 2 MN i = - MN / 2 + 1 MN / 2 n 0 ( k + i ) · cos 2 π ( k + i N + 1 8 ) = a sin 2 π ( Δ fk + θ 0 - 1 8 ) + n 5 ( k ) + n 6 ( k ) + n 7 ( k ) ( 3 )
  • Again, equation (3) is obtained by substituting equation (1) into Pcos(k), for y(k). The variables n5, n6, n7 denote various noise components, resulting from white noise or the offset in frequency.
  • The phase estimate may be determined for the kth sample, where the
  • phase at the kth sample is given by
  • φ ( k ) = Δ fk + θ 0 - 1 8 ,
  • using the following equation:
  • Q ( k ) P cos ( k ) P sin ( k ) a sin 2 π φ ( k ) + n 7 ( k ) a cos 2 π φ ( k ) + n 3 ( k ) tan 2 πφ ( k ) + n Q ( k ) ( 4 )
  • The phase estimate may be found by taking the arctangent of Q(k), and may be represented by the following equation:
  • φ ^ ( k ) 1 2 π arc tan Q ( k ) = φ ( k ) + n φ ( k ) ( 5 )
  • In practice, one may use a look-up table to calculate the arctangent function, using the memory, as shown in FIG. 3. The variables nQ(k) and nφ(k) represent noise in each of equations (4) and (5).
  • In order to properly equalize the read channel, it may be necessary to determine one or more gain estimates or amplitude estimates using the preamble of the incoming digital samples. Aspects of the invention allow computation of a gain estimate by taking the square root of the sum of squares of the sine and cosine projections. The sine and cosine projections (or components) are referenced using equations (2) and (3). This may be represented by the following equation:

  • ã(k)≡√{square root over (P sin(k)2 +P cos(k)2)}{square root over (P sin(k)2 +P cos(k)2)}a+cos 2πφ(k)n 7(k)+sin 2πφ(k)n 3(k)  (6)
  • As shown in equation (6), the amplitude or gain is found by calculating the square root of the sum of the squares of the sine and cosine projections.
  • FIG. 4 is a block diagram of a subsystem of a read channel receiver that computes frequency and phase offsets used to synchronize the read channel receiver of a magnetic hard disk drive, in accordance with an embodiment of the invention. The frequency and phase computational subsystem, shown in FIG. 4, comprises a first accumulator 404, a second accumulator 412, a counter 424 for the integer variable, k, which indexes the one or more preliminary phase offset estimates that are input, a first multiplier 428, a second multiplier 432, a third multiplier 436, a fourth multiplier 440, a fifth multiplier 444, a phase offset adder 448, and a frequency offset adder 452. The first accumulator 404 comprises a first accumulator adder 412 and a first register 408 while the second accumulator 412 comprises a second accumulator adder 420 and a second register 416. The first accumulator 404 accumulates (or provides a cumulative sum of) preliminary phase offset estimates provided by the phase and gain computational subsystem. The second accumulator 412 accumulates successive terms and functionally operates similar to that of the first accumulator 404. The two factors that generate the successive terms (product of two factors) comprise a count of the counter 424 and its associated preliminary phase offset estimate. Of course, the first adder 412 within the first accumulator 404 adds the preliminary phase offset estimate to the accumulated value stored in the first register 408. Likewise, the second adder 416 within the second accumulator 412 adds the terms (associated with the product of the preliminary phase offset estimate and the count) to the accumulated value stored in the second register 416. The values stored in the accumulators 404, 412 are used to compute the frequency and phase offsets. The frequency offset is computed by summing two products, a first product and a second product. The first product is obtained by multiplying the value in the first accumulator 404 with β(L) while the second product is obtained by multiplying the value in the second accumulator 412 with α(L). The phase offset is computed by summing two products (a third product and a fourth product) and a constant. The third product is obtained by multiplying the value in the first accumulator 404 with δ(L) while the fourth product is obtained by multiplying the value in the second accumulator 412 with γ(L). The third product, fourth product, and a constant (⅛ in this representative embodiment) are added to generate the phase offset. The aforementioned counter, adders, multipliers, and registers may be implemented using digital logic circuitry.
  • Aspects of the invention comprise utilizing a linear regression technique for determining the phase and frequency offsets of a read channel in a read channel receiver. The technique is employed in implementing the embodiment illustrated in FIG. 4. The objective of the linear regression technique is to find ξ and ζ using a linear regression equation, such as yk≈ξxk+ζ for all k, given the set of points with coordinates x=(x1, x2, . . . , xL)7 and y=(y1, y2, . . . , yL)T. More precisely, the variables, ξ and ζ are found that minimizes the L-2 norm of the error. The error may be written in vector form as F(ξ,ζ)=∥y−ξx−ζ∥2=(y−ξx−ζ)T(y−ξx−ζ). The function F(ξ,ζ) is a quadratic function in ξ and ζ so to find the minimum it is sufficient to take the partial derivatives and set them to zero. The partial derivatives of F(ξ,ζ) with respect to ξ and ζ are as follows;
  • F ξ = - 2 x 7 ( y - ξ x - ζ ) F ζ = - 2 I 1 × L ( y - ξ x - ζ )
  • The resulting equations after setting the partial derivatives to zero are as follows:
  • F ξ = - 2 x 7 ( y - ξ x - ζ ) = 0 , F ζ = - 2 I 1 × L ( y - ξ x - ζ ) = 0
  • The equations may be re-written as follows:
  • { x T y - ( x T x ) ξ - ( x T I L × 1 ) ζ = 0 I 1 × L y - ( I 1 × L x ) ξ - ( I 1 × L I L × 1 ) ζ = 0 where I 1 × L denotes a 1 by L ( 7 )
  • vector of all ones. The equations, after rearranging, may be re-written in matrix form as follows:
  • ( ξ ζ ) = ( x T x x T I L × 1 I 1 × L x L ) - 1 ( x T y I 1 × L y ) ( 8 )
  • In a representative embodiment, the vector x comprises sampling times. In a representative embodiment, the sampling times comprises fixed sampling intervals. Hence, for example, the vector x may comprise x=(1, . . . , L)T, where L corresponds to the number of samples used in the linear regression technique (i.e., number of preliminary phase estimates used to compute a frequency offset and a phase offset). The inverse of the matrix shown in equation (8) only depends on the variable, L, and does not depend on y, and as a result, may be determined in advance. The terms within the invertible 2×2 matrix of equation (8) may be represented by the following equation:
  • ( ξ ζ ) = ( α ( L ) β ( L ) γ ( L ) δ ( L ) ) ( x T y I 1 × L y ) ( 9 )
  • As shown in equation (9), the terms in the 2×2 matrix of equation (9), α(L), β(L), γ(L), and ξ(L), are all functions of the predetermined variable L. These terms may be obtained by inverting the 2×2 matrix shown in equation (8). The 2×1 matrix,
  • ( x T y I 1 × L y ) ,
  • of equation (9) may be implemented using two registers, as shown in equations (10a) and (10b).

  • R xy(k)=R xy(k−1)+ky k,  (10a)

  • R y(k)=R y(k−1)+y k  (10b)
  • This leads to a simple implementation for determining ξ and ζ, when Rxy and Ry are substituted into
  • ( x T y I 1 × L y ) .
  • In this representative embodiment, the two registers, Rxy and Ry are substituted into the matrix equation (9). The matrix equation (9) may be expanded into the following two equations (11a) and (11b):

  • ξ=αR xy +βR y,  (11a)

  • ζ=γR xy +δR y  (11b)
  • Recalling that
  • φ ( k ) = Δ fk + θ 0 - 1 8 ,
  • one may apply the aforementioned linear regression equation, yk≈ξxk+ζ for all k, to solve for Δf and θ0. Using equations (10a) and (10b), and applying the sequence of preliminary phase offset estimates {circumflex over (φ)}(1),{circumflex over (φ)}(2), . . . , {circumflex over (φ)}(L) provided by the phase and gain computational circuitry, the following equations are obtained, in which Δ{circumflex over (f)}(L) corresponds to the frequency offset while {circumflex over (θ)}0(L) corresponds to the phase offset.

  • R xy(k)=R xy(k−1)+k{circumflex over (φ)}(k)  (12a)

  • R y(k)=R y(k−1)+{circumflex over (φ)}(k)  (12b)
  • Δ f ^ ( k ) = α ( k ) R xy + β ( k ) R y ( 13 a ) θ ^ 0 ( k ) = 1 8 + γ ( k ) R xy + δ ( k ) R y ( 13 b )
    R xy(0)=0,R y,(0)=0  (14)
  • In this representative embodiment, the equations may be used to implement a system that allows one to determine the values for Δ{circumflex over (f)}(k) and {circumflex over (θ)}0(k) (mathematical terminology representing the estimates for the frequency and phase offsets, respectively) for any value of k={1, 2, 3, . . . , L}. As a consequence, the embodiment provides an implementation that adjusts to a maximum number of digital samples, L, as one or more preambles of one or more disk sectors are input into the read channel of a magnetic hard disk drive.
  • Equations (12a) and (12b) represent the implementation of the first portion of FIG. 4 while equations (13a) and (13b) represent the implementation of the second portion of FIG. 4. Equations (14a) and (14b) indicate that the registers, Rxy and Ry used to implement the embodiment shown in FIG. 4, are initialized with value 0 when k=0.
  • FIG. 5 is an operational flow diagram illustrating a method of generating frequency offset, phase offset, and gain (or amplitude) of a data stream read from data stored in a data sector, by a read channel receiver of a hard disk drive, in accordance with an embodiment of the invention. At step 504, a preamble is read from the data sector of one or more data sectors in a hard disk drive. The one or more bytes contained within the preamble is used to generate the frequency offset, phase offset, and gain of the data stream. In a representative embodiment, the number of bytes used in the preamble is less than the typical number of bytes used—15 bytes After processing and conditioning is performed on the sinusoidal output generated by the read channel receiver in the hard disk drive, a sequence of digitized samples is provided to the phase and gain computational subsystem. As was referenced in the embodiment of FIG. 2, a processing/conditioning circuitry performs the processing and conditioning. The sequence of digitized samples comprises the latest fixed number of digitized samples. At step 506, the last fixed number of digitized samples is provided to the phase and gain computational subsystem In a representative embodiment, this fixed number, x, comprises the number of samples used to determine the sine and cosine projections using the phase and gain computational subsystem. In a representative embodiment, the value for the fixed number, x, is 8. This fixed number may be programmed by a user, prior to performing the computations provided by the phase and gain computational subsystem. At step 508, the phase and gain computational subsystem determines the sine and cosine projections (or sine and cosine components) based on the sequence of digitized samples provided Equations (2) and (3) are used to compute these projections. At step 512, Q(k) the quotient is computed for P cos(k)/P sin(k). At step 516, a look-up table (implemented in memory, for example) is used to determine the arctangent of Q(k), Φ(κ). At step 520, the gain or amplitude is output by the phase and gain computational subsystem, for use by the hard disk drive. At step 524, the phase and gain computational subsystem outputs Φ(κ) (preliminary phase offset estimates) into the frequency and phase computational subsystem, in which, the processing described in reference to FIG. 4, is accomplished. Finally, at step 528, the frequency and phase computational subsystem outputs the phase and frequency offsets associated with the sequence of samples provided to the frequency and phase computational subsystem. Next, at step 530, an assessment is made whether the last sample has been generated by the preamble. If the end of the preamble of the data sector has been sampled, then the process ends for this particular data sector. Thereafter, a new data sector may be read, if necessary, which contains another preamble. If a new data sector is read, then the process commences with another cycle, starting at step 504. Otherwise, the process resumes with step 532, at which, the next digitized sample is introduced at the next sample clock, using the remaining un-sampled preamble data. Then, the process repeats itself at step 504, when a new sequence of digitized samples is introduced to the phase and gain computational subsystem. In this exemplary embodiment, 7 of the previous 8 samples are used in the next sequence of samples. The 7 previously used samples are shifted by one sample clock, so as to accommodate the latest new sample.
  • In one or more embodiments, additional performance improvements may be obtained by utilizing knowledge of noise statistics and applying this to the noise variables previously used in the preceding equations. The performance improvements rely on an analysis of the noise variables n1, n2, n3, n5, n6, and n7 used in equations (2) and (3). For example, Bayesian Linear Regression may be employed when a priori knowledge of the frequency offset distribution is used. In another example, Colored Noise Linear Regression may be employed when noise correlation is taken into account while performing a linear regression.
  • The following set of equations provides a description of the noise variables, n1 and n2, in reference to equation (2) previously described.
  • n 1 ( k ) 2 MN i = - MN / 2 + 1 MN / 2 a · sin 2 2 π ( k + i N + 1 8 ) [ cos 2 π ( Δ fk + θ 0 - 1 8 ) - sin 2 π ( Δ fk + θ 0 - 1 8 ) · 2 πΔ fi ] - 2 MN i = - MN / 2 + 1 MN / 2 a · sin 2 2 π ( k + i N + 1 8 ) cos 2 π ( Δ fk + θ 0 - 1 8 ) = 2 MN ( - 2 πΔ fa ) [ i = - MN / 2 + 1 MN / 2 i · sin 2 2 π ( k + i N + 1 8 ) ] sin 2 π ( Δ fk + θ 0 - 1 8 ) = 2 MN ( - 2 πΔ fa ) d 1 ( k ) sin 2 π ( Δ fk + θ 0 - 1 8 ) n 2 ( k ) 2 MN ( 2 π Δ fa ) [ i = - MN / 2 + 1 MN / 2 i · cos 2 π ( k + i N + 1 8 ) sin 2 π ( k + i N + 1 8 ) ] cos 2 π ( Δ fk + θ 0 - 1 8 ) = 2 MN ( 2 πΔ fa ) d 2 ( k ) cos 2 π ( Δ fk + θ 0 - 1 8 )
  • d1(k) and d2(k) may be calculated as:
  • d 1 ( k ) = 2 MN ( - 1 ) i = - MN / 2 + 1 MN i · 1 2 [ 1 - cos 4 π ( k + i 8 + 1 8 ) ] = b cos 2 π ( 2 k N + θ d ) - 1 2 d 2 ( k ) = 2 MN i = - MN / 2 + 1 MN / 2 I · 1 2 sin 4 π ( k + i N + 1 8 ) = b sin 2 π ( 2 k N + θ d )
  • where b and θd are constant w.r.t. k and almost constant within a range of M and N:
  • Further, for b(M, N): b(1,4)=b(1,8)=b(2,4)=b(2,8)=½;b(1,16)=b(2,16)≈1.2 Hence,
  • n 1 ( k ) + n 2 ( k ) = - πΔ fa sin 2 π ( Δ fk + θ 0 - 1 8 ) + 2 πΔ fab sin 2 π ( ( 2 N + Δ f ) k + θ d + θ 0 - 1 8 ) ,
  • which is negligible since Δf<<θn 0 ·σn 0 is the covariance of n0(k).
  • Similarly, the following set of equations provides a description of the noise variables, n3 and n4, in relation to equation (3) previously described
  • n 5 ( k ) 2 MN 2 π Δ fad 3 ( k ) sin 2 π ( Δ fk + θ 0 - 1 8 ) n 6 ( k ) = 2 MN 2 πΔ fad 4 ( k ) cos 2 π ( Δ fk + θ 0 - 1 8 ) d 3 ( k ) = - b sin 2 π ( 2 k N + θ d ) d 4 ( k ) = b cos 2 π ( 2 k N + θ d ) + 1 2
  • Hence,
  • n 5 ( k ) + n 6 ( k ) = - πΔ fa cos 2 π ( Δ fk + θ 0 - 1 8 ) + 2 πΔ fab cos 2 π ( ( 2 N + Δ f ) k + θ d + θ 0 - 1 8 )
  • is also negligible when Δf<<σn 0.
  • The following set of equations provides a description of the noise variable n3:
  • As previously discussed,
  • n 3 ( k ) = 2 MN i = - MN / 2 + 1 MN / 2 n 0 ( k + i ) · sin 2 π ( k + i N + 1 8 )
  • Assuming n0 is AWGN (additive white Gaussian noise), then n3 is zero-mean Gaussian. The correlation matrix of n3 is as follows:
  • r n 3 ( k , k - l ) E [ n 3 ( k ) n 3 ( k - l ) ] = E [ 2 MN i = - MN / 2 + 1 MN / 2 n 0 ( k + i ) · sin 2 π ( k + i N + 1 8 ) · 2 MN j = - MN / 2 + 1 MN / 2 n 0 ( k + j ) · sin 2 π ( k + j N + 1 8 ) ] = ( 2 MN ) 2 σ n 0 2 [ - MN / 2 + 1 i MN / 2 - MN / 2 + 1 j = i + l MN / 2 sin 2 2 π ( k + i N + 1 8 ) ] = ( 2 MN ) 2 σ n 0 2 [ i = - MN / 2 + 1 MN / 2 sin 2 2 π ( k + i N + 1 8 ) - l 0 : MN / 2 - l + 1 i MN / 2 l < 0 : - MN / 2 + 1 i - MN / 2 + 1 sin 2 2 π ( k + i N + 1 8 ) ] ( l MN - 1 ) = σ n 0 2 [ 2 MN - ( 2 MN ) 2 l 0 : MN / 2 - l + 1 i MN / 2 l < 0 : - MN / 2 + 1 i - MN / 2 + 1 sin 2 2 π ( k + i N + 1 8 ) ] ( l MN - 1 ) r n 3 ( k , k - l ) = 0 ( otherwise )
  • When N=4 as may be found used in one or more embodiments of the present invention, rn 3 (k,k−1) is constant w.r.t. k.
  • r n 3 ( k , k - l ) = { σ n 0 2 2 MN ( 1 - l MN ) , l MN - 1 0 , l MN
  • The correlation function takes a similar form.
  • Fr n 3 ( l ; L ) E [ 1 L k = 0 L - 1 n 3 ( k ) n 3 ( k - l ) ] = 1 L k = 0 L - 1 r n 3 ( k , l ) = { σ n 0 2 [ 2 MN - ( 2 MN ) 2 l 0 : MN / 2 - l + 1 i MN / 2 l < 0 : - MN / 2 + 1 i - MN / 2 + l ( 1 L k = 0 L - 1 sin 2 2 π ( k + i N + 1 8 ) ) ] , l MN - 1 0 , l MN { σ n 0 2 [ 2 MN - ( 2 MN ) 2 l 0 : MN / 2 - l + 1 i MN / 2 l < 0 : - MN / 2 + 1 i - MN / 2 + l ( 1 2 ) ] , l MN - 1 ( iff L >> N / 2 ; or L = iN / 2 , i I ) 0 , l MN = { σ n 0 2 2 MN ( 1 - l MN ) , l MN - 1 0 , l MN
  • The variable, L, comprises the total number of samples provided by the preamble of the data sector. The noise power density spectrum may be obtained by taking the fast Fourier transform of the correlation function as follows:

  • N 3(k ω)=FFT(Fr n 3 (l))
  • It may be shown that
  • k ω N 2 ( k ω ) = Fr n 3 ( 0 ) = σ n 0 2 2 MN 1 MN ,
  • which indicates that a reduction in the sine projection noise is proportional to the product M*N. Hence, aspects of the invention may incorporate values for N and M that provide embodiments having significantly reduced projection noise.
  • Another interesting property is that the dc noise is invariant w.r.t. M*N:
  • N 3 ( 0 ) = l Fr n 3 ( l ) σ n 0 2
  • As may be expected, the correlation matrix of n7 in formula (3) is similar to that of n3:
  • r n 7 ( k , k - l ) = { σ n 0 2 [ 2 MN - ( 2 MN ) 2 l 0 : MN / 2 - l + 1 i MN / 2 l < 0 : MN / 2 + 1 i - MN / 2 + l cos 2 2 π ( k + l N + 1 8 ) ] , l MN - 1 0 , l MN When N = 4 , r n 7 ( k , k - l ) is constant w . r . t . k 0 n r n 7 ( k , k - l ) = { σ n 0 2 2 MN ( 1 - l MN ) , l MN - 1 0 , l MN = r n 3 ( k , k - l )
  • Hence, Frn 3 (0)=Frn 7 (0)
  • When performing a preliminary phase offset estimate as indicated by way of equation (4), the noise properties may be derived as follows:
  • Q ( k ) n 3 ( k ) | P sin ( k ) = a cos 2 π φ ( k ) P cos ( k ) = a sin 2 π φ ( k ) = - sin 2 πφ ( k ) a cos 2 2 πφ ( k ) Q ( k ) n 7 ( k ) | P sin ( k ) = a cos 2 π φ ( k ) P cos ( k ) = a sin 2 π φ ( k ) = - 1 a cos 2 2 πφ ( k ) n Q ( k ) = - sin 2 πφ ( k ) a cos 2 2 πφ ( k ) n 3 ( k ) + 1 a cos 2 πφ ( k ) n 7 ( k )
  • Recall that the preliminary phase offset estimate is given by
  • φ ^ ( k ) 1 2 π arctan Q ( k ) = φ ( k ) + n φ ( k ) .
  • Therefore,
  • n φ ( k ) 1 2 π cos 2 2 πφ ( k ) n Q ( k ) 1 2 π a ( - sin 2 πφ ( k ) n 3 ( k ) + cos 2 πφ ( k ) n 7 ( k ) )
  • And its correlation matrix is given by:
  • r n φ ( k , k - l ) E n φ ( k ) n φ ( k - l ) ( 1 2 π a ) 2 ( sin 2 πφ ( k ) sin 2 πφ ( k - l ) r n 3 ( k , k - l ) + cos 2 πφ ( k ) cos 2 πφ ( k - l ) r n 7 ( k , k - l ) )
  • In above equation, the near-zero cross-correlation property of n3 and n7 is utilized, which may be proven by way of the following equations:
  • xr n 3 , n 7 ( k , k - l ) E [ n 3 ( k ) n 7 ( k - l ) ] = E [ 2 MN i = - MN / 2 + 1 MN / 2 n 0 ( k + i ) · sin 2 π ( k + i N + 1 8 ) · 2 MN j = - MN / 2 + 1 MN / 2 n 0 ( k + j ) · cos 2 π ( k + j N + 1 8 ) ] = ( 2 MN ) 2 σ n 0 2 [ - MN / 2 + 1 i MN / 2 - MN / 2 + 1 j = i + l MN / 2 1 2 sin 4 π ( k + i N + 1 8 ) ] = { 2 MN σ n 0 2 [ 1 MN l 0 : 1 i MN - l l < 0 : 1 - l i MN sin 4 π ( k + i N + 1 8 ) ] , MN - N / 2 l MN - 1 0 , l MN - 1 - N / 2 = { 2 MN σ n 0 2 ( - 1 MN cos k π ) , l = MN - 1 0 , l < MN - 1 ( if N = 4 ) 0
  • Thus, using the cross-correlation of n3 and n7 analysis (assuming that N=4, and rn 3 =rn 7 the correlation of the noise associated with the preliminary phase offset estimate is given by:
  • r n φ ( k , k - l ) = 2 MN σ 0 2 ( 1 - l MN ) cos 2 πΔ fl 2 MN σ 0 2 ( 1 - l MN ) ( Λ fl << 1 ) where , σ 0 = σ n 0 2 π a
  • An estimate of the noise level may be based on the following, such that autocorrelation of {tilde over (α)}(k) is zero when the delay is no larger than M*N
  • σ ^ n 0 2 MN 4 1 L - mn k = nm L - 1 ( a ~ ( k ) - a ~ ( k - mn ) ) 2 MN 4 E [ ( a ~ ( k ) - a ~ ( k - mn ) ) 2 ] MN 4 [ cos 2 2 πφ ( k ) σ n 7 2 - cos 2 2 πφ ( k - mn ) σ n 7 2 + sin 2 2 πφ ( k ) σ n 3 2 - sin 2 2 πφ ( k - mn ) σ n 3 2 ] = MN 4 2 σ n 3 2 ( N = 4 σ n 3 2 = σ n 7 2 , or , σ n 3 2 σ n 7 2 ) = σ n 0 2
  • In one or more embodiments, two alternative methods may be used for linear regression: linear regression with colored noise, and Bayesian linear regression. In the embodiment using linear regression with colored noise, when the correlation of phase estimation noise is considered, the log-likelihood of {nφ(k)} becomes F(Δf,θ′0)=({circumflex over (φ)}−Δfk−θ′0)TRn φ −1({circumflex over (φ)}−Δfk−θ′0), where bold variables comprise matrices or vectors:

  • R n φ −1=(r n φ (k,k−1))k,k-1,{circumflex over (φ)}={φ(k)|k=0, . . . , L−1},k={k|k=0, . . . , L−1}
  • and θ′00−⅛
  • From ∂F(Δf,θ′0)/θΔf=0 and ∂F(Δf,θ′0)/∂θ′0=0, one obtains the matrix equations:
  • [ Δ f θ 0 ] = [ k T R n ϕ - 1 k k T R n ϕ - 1 i n × 1 i n × 1 T R n ϕ - 1 k i n × 1 T R n ϕ - 1 i n × 1 ] - 1 [ k 7 R n ϕ - 1 ϕ ^ i n × 1 T R n ϕ - 1 ϕ ^ ] + [ 0 1 / 8 ]
  • This formula replaces formulas (10), (11), (12), (13), and (14), previously described.
  • When prior information of Δf and θ′0 is available, or the information has been adaptively learned real time, Bayesian linear regression may be employed. The likelihood function of the preliminary phase offset estimate is given by:

  • l(f,θ′ 0;{circumflex over (φ)})=p({circumflex over (φ)}|Δf,θ′ 0)=(2π)−1/2(detR n φ )−1exp(−({circumflex over (φ)}−Δfk−θ′ 0)T R n φ −1({circumflex over (φ)}−Δfk−θ′D)/2)
  • The preceding equation uses a Gaussian distribution for Δf and an even distribution results for θ′0;
  • p Δ f ^ , θ ^ 0 ( Δ f , θ 0 ) = 1 2 π σ Δ f exp ( - ( Δ f ) 2 σ Δ f 2 ) · p θ ^ 0 ( θ 0 ) ,
  • where
  • p θ ~ 0 = { 1 θ 0 < 1 2 0 θ 0 >= 1 2
  • is the probability density function of the initial phase offset.
  • Thus, Bayesian estimation of Δf is defined as:
  • Δ f ^ = θ 0 Δ f l ( f , θ 0 ; ϕ ^ ) p Δ f ^ , θ ^ 0 ( Δ f , θ 0 ) · Δ f Δ f θ 0 θ 0 Δ f l ( f , θ 0 ; ϕ ^ ) p Δ f ^ , θ ^ 0 ( Δ f , θ 0 ) · Δ f θ 0 = θ 0 [ Δ f exp ( a + b Δ f + ( c + h ) ( Δ f ) 2 ) Δ f · Δ f ] · θ 0 θ 0 Δ f exp ( a + b Δ f + ( c + h ) ( Δ f ) 2 ) · Δ f θ 0 where a = - ( ϕ - θ 0 ) T R n ϕ - 1 ( ϕ ^ - θ 0 ) / 2 b = k T R n ϕ - 1 ( ϕ ^ - θ 0 ) c = - k T R n ϕ - 1 k / 2 h = - 1 2 σ Δ f 2
  • In a representative embodiment, mathematical approximation (σ0 2<<1) yields the following concise formulations:

  • Δ{circumflex over (f)}=f{circumflex over (φ)}
  • f = - k T R n ϕ - 1 2 ( c + h ) ( l n × n + R i 1 × n Ri n × 1 ) R = R n ϕ - 1 + R n ϕ - 1 kk T R n ϕ - 1 2 ( c + h )
  • In the event that σ0 2 is stable in different sections, one or more old a σ0 2 values may be used to calculate f in advance. In the event f varies in different sections, σ0 2 may be estimated in real-time and the formula may require a small adjustment yielding:
  • Δ f ^ = f 1 ϕ ^ + γ f 2 ϕ ^ f 3 + γ f 4 + γ 2 f 5 , γ = σ 0 2 σ Δ f 2 ,
  • wherein f1, f2, f3, f4, f5, are constants.
  • While the invention has been described with reference to certain embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teachings of the invention without departing from its scope. Therefore, it is intended that the invention not be limited to the particular embodiment disclosed, but that the invention will include all embodiments falling within the scope of the appended claims.

Claims (31)

1-37. (canceled)
38. A method of estimating gain or amplitude of a data stream that is read from a hard disk drive comprising:
receiving a plurality of digitized samples from reading a preamble field of a sector of a hard drive;
determining a sine component for each of said plurality of digitized samples;
determining a cosine component for each of said plurality of digitized samples;
computing a sum of squares of said sine and cosine components; and
computing a square root of said sum to generate said gain or amplitude.
39. The method of claim 38 wherein said preamble field contains less than 15 bytes of data.
40. The method of claim 38 wherein said preamble field contains 6 bytes of data.
41. The method of claim 38 wherein said preamble field precedes a sync field, said sync field used to notify or indicate that a data field immediately follows said sync field.
42. The method of claim 38 wherein said plurality of digitized samples equals 8 digitized samples.
43. The method of claim 42 wherein said plurality of digitized samples represents two periods of a sinusoidal waveform.
44. A system for estimating gain or amplitude of data that is read from a hard disk drive, said system comprising:
at least one circuitry operative for, at least:
receiving a plurality of digitized samples from reading a preamble field of a sector of a hard drive;
determining a sine component for each of said plurality of digitized samples;
determining a cosine component for each of said plurality of digitized samples;
computing a sum of squares of said sine and cosine components; and
computing a square root of said sum to generate said gain or amplitude.
45. The system of claim 44 wherein said preamble field contains less than 15 bytes of data.
46. The system of claim 44 wherein said preamble field contains 6 bytes of data.
47. The system of claim 44 wherein said preamble field precedes a sync field, said sync field used to notify or indicate that a data field immediately follows said sync field.
48. The system of claim 45 wherein said plurality of digitized samples equals 8 digitized samples.
49. The system of claim 48 wherein said plurality of digitized samples represents two periods of a sinusoidal waveform.
50. A system for estimating gain or amplitude of data that is read from a hard disk drive, said system comprising:
a memory for storing a set of executable instructions;
a processor for executing said executable instructions, wherein said executing performs:
receiving a plurality of digitized samples from reading a preamble field of a sector of a hard drive;
determining a sine component for each of said plurality of digitized samples;
determining a cosine component for each of said plurality of digitized samples;
computing a sum of squares of said sine and cosine components; and
computing a square root of said sum to generate said gain or amplitude.
51. The system of claim 50 wherein said preamble field contains less than 15 bytes of data.
52. The system of claim 50 wherein said preamble field contains 6 bytes of data.
53. The method of claim 50 wherein said preamble field precedes a sync field, said sync field used to notify or indicate that a data field immediately follows said sync field.
54. The method of claim 50 wherein said plurality of digitized samples equals 8 digitized samples.
55. The method of claim 54 wherein said plurality of digitized samples represents two periods of a sinusoidal waveform.
56. A method for optimizing performance of a hard drive by synchronizing a read channel receiver of said hard drive to a data stream that is read, said method comprising:
receiving a plurality of digitized samples from reading a preamble field of a sector of a hard drive;
first determining a sine component for each of said plurality of digitized samples;
second determining a cosine component for each of said plurality of digitized samples;
dividing said cosine component with said sine component to generate a quotient;
computing an arctangent of said quotient to generate a preliminary phase estimate of each of said plurality of digitized samples;
using at least two of said plurality of digitized samples to generate at least two preliminary phase estimates based on said first determining, second determining, dividing, and computing;
generating a frequency offset and a phase offset by applying linear regression to said at least two preliminary phase estimates and an integer variable, said integer variable used for indexing each of said plurality of digitized samples; and
using said frequency offset and said phase offset to perform said synchronizing.
57. The method of claim 56 wherein said preamble field contains less than 15 bytes of data.
58. The method of claim 56 wherein said preamble field contains 6 bytes of data.
59. The method of claim 56 wherein said preamble field precedes a sync field, said sync field used to notify or indicate that a data field immediately follows said sync field.
60. The method of claim 56 wherein said plurality of digitized samples equals 8 digitized samples.
61. The method of claim 60 wherein said plurality of digitized samples represents two periods of a sinusoidal waveform.
62. A system for optimizing performance of a hard drive by synchronizing a read channel receiver of said hard drive to a data stream that is read, said system comprising:
at least one circuitry operative for, at least:
receiving a plurality of digitized samples from reading a preamble field of a sector of a hard drive;
first determining a sine component for each of said plurality of digitized samples;
second determining a cosine component for each of said plurality of digitized samples;
dividing said cosine component with said sine component to generate a quotient;
computing an arctangent of said quotient to generate a preliminary phase estimate of each of said plurality of digitized samples;
using at least two of said plurality of digitized samples to generate at least two preliminary phase estimates based on said first determining, second determining, dividing, and computing;
generating a frequency offset and a phase offset by applying linear regression to said at least two preliminary phase estimates and an integer variable, said integer variable used for indexing each of said plurality of digitized samples; and
using said frequency offset and said phase offset to perform said synchronizing.
63. The system of claim 62 wherein said preamble field contains less than 15 bytes of data.
64. The system of claim 62 wherein said preamble field contains 6 bytes of data.
65. The system of claim 62 wherein said preamble field precedes a sync field, said sync field used to notify or indicate that a data field immediately follows said sync field.
66. The system of claim 62 wherein said plurality of digitized samples equals 8 digitized samples.
67. The system of claim 66 wherein said plurality of digitized samples represents two periods of a sinusoidal waveform.
US12/178,426 2004-07-08 2008-07-23 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive Abandoned US20080278837A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US12/178,426 US20080278837A1 (en) 2004-07-08 2008-07-23 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US10/887,544 US7405894B2 (en) 2004-07-08 2004-07-08 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive
US12/178,426 US20080278837A1 (en) 2004-07-08 2008-07-23 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US10/887,544 Continuation US7405894B2 (en) 2004-07-08 2004-07-08 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive

Publications (1)

Publication Number Publication Date
US20080278837A1 true US20080278837A1 (en) 2008-11-13

Family

ID=35541080

Family Applications (2)

Application Number Title Priority Date Filing Date
US10/887,544 Expired - Fee Related US7405894B2 (en) 2004-07-08 2004-07-08 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive
US12/178,426 Abandoned US20080278837A1 (en) 2004-07-08 2008-07-23 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive

Family Applications Before (1)

Application Number Title Priority Date Filing Date
US10/887,544 Expired - Fee Related US7405894B2 (en) 2004-07-08 2004-07-08 Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive

Country Status (1)

Country Link
US (2) US7405894B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8947811B1 (en) * 2013-03-12 2015-02-03 Western Digital Technologies, Inc. Disk drive using preamble and postamble to adjust asynchronous signal samples of payload

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100660839B1 (en) * 2004-10-07 2006-12-26 삼성전자주식회사 Joint symbol timing and carrier phase recovery circuit for ATSC receiver
US8199422B2 (en) 2010-07-30 2012-06-12 Lsi Corporation Methods and apparatus for gain estimation using servo data with improved bias correction
US8824087B2 (en) * 2012-08-02 2014-09-02 Lsi Corporation Zero gain start bias estimation
US8780470B1 (en) 2013-03-11 2014-07-15 Western Digital Technologies, Inc. Disk drive adjusting digital phase locked loop over sector data with frequency induced phase error measured over preamble

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5001724A (en) * 1989-01-13 1991-03-19 Hewlett-Packard Company Method and apparatus for measuring phase accuracy and amplitude profile of a continuous-phase-modulated signal
US5170415A (en) * 1989-06-14 1992-12-08 Nec Corporation Burst demodulator for establishing carrier and clock timing from a sequence of alternating symbols
US5202842A (en) * 1989-11-29 1993-04-13 Kabushiki Kaisha Okuma Tekkosho Rotational position detecting device which compensates for eccentricity of rotating object
US5297185A (en) * 1992-02-13 1994-03-22 North American Philips Corporation Pattern detection and synchronization circuit
US5341249A (en) * 1992-08-27 1994-08-23 Quantum Corporation Disk drive using PRML class IV sampling data detection with digital adaptive equalization
US5583705A (en) * 1994-10-04 1996-12-10 Quantum Corporation Write precompensation optimization in a PRML channel using a selected PRML signal level
US5835295A (en) * 1996-11-18 1998-11-10 Cirrus Logice, Inc. Zero phase restart interpolated timing recovery in a sampled amplitude read channel
US5862005A (en) * 1994-10-11 1999-01-19 Quantum Corporation Synchronous detection of wide bi-phase coded servo information for disk drive
US6023386A (en) * 1997-10-31 2000-02-08 Cirrus Logic, Inc. Fault tolerant sync mark detector for synchronizing a time varying sequence detector in a sampled amplitude read channel
US6122117A (en) * 1998-02-04 2000-09-19 Fujitsu Limited Method for digitally demodulating position signals recorded on recording media
US6204988B1 (en) * 1998-08-24 2001-03-20 Western Digital Corporation Disk drive capable of autonomously evaluating and adapting the frequency response of its servo control system
US6239934B1 (en) * 1994-03-11 2001-05-29 Kai C. K. Sun Disk drive system and control circuitry adapted for downward compatibility
US6704156B1 (en) * 2000-01-31 2004-03-09 Maxtor Corporation Self-writing of servo patterns in a disk drive using a printed reference pattern
US6754025B1 (en) * 2000-01-12 2004-06-22 Maxtor Corporation Disk drive spindle motor speed and timing control
US6980389B1 (en) * 2004-02-27 2005-12-27 Western Digital Technologies, Inc. Disk drive having a disk including a servo burst pattern having a normal servo burst field and a quadrature servo burst field
US6989956B2 (en) * 2002-07-03 2006-01-24 Fujitsu Limited Head position control method, disk device, and servo track write method
US7054088B2 (en) * 2002-01-21 2006-05-30 Fujitsu Limited Information recording and reproducing apparatus and method, and signal decoding circuit for performing timing recovery
US7400561B2 (en) * 2004-12-14 2008-07-15 Lsi Corporation Single PLL demodulation of pre-formatted information embedded in optical recording medium

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5001724A (en) * 1989-01-13 1991-03-19 Hewlett-Packard Company Method and apparatus for measuring phase accuracy and amplitude profile of a continuous-phase-modulated signal
US5170415A (en) * 1989-06-14 1992-12-08 Nec Corporation Burst demodulator for establishing carrier and clock timing from a sequence of alternating symbols
US5202842A (en) * 1989-11-29 1993-04-13 Kabushiki Kaisha Okuma Tekkosho Rotational position detecting device which compensates for eccentricity of rotating object
US5297185A (en) * 1992-02-13 1994-03-22 North American Philips Corporation Pattern detection and synchronization circuit
US5341249A (en) * 1992-08-27 1994-08-23 Quantum Corporation Disk drive using PRML class IV sampling data detection with digital adaptive equalization
US6239934B1 (en) * 1994-03-11 2001-05-29 Kai C. K. Sun Disk drive system and control circuitry adapted for downward compatibility
US5583705A (en) * 1994-10-04 1996-12-10 Quantum Corporation Write precompensation optimization in a PRML channel using a selected PRML signal level
US5862005A (en) * 1994-10-11 1999-01-19 Quantum Corporation Synchronous detection of wide bi-phase coded servo information for disk drive
US5835295A (en) * 1996-11-18 1998-11-10 Cirrus Logice, Inc. Zero phase restart interpolated timing recovery in a sampled amplitude read channel
US6023386A (en) * 1997-10-31 2000-02-08 Cirrus Logic, Inc. Fault tolerant sync mark detector for synchronizing a time varying sequence detector in a sampled amplitude read channel
US6122117A (en) * 1998-02-04 2000-09-19 Fujitsu Limited Method for digitally demodulating position signals recorded on recording media
US6204988B1 (en) * 1998-08-24 2001-03-20 Western Digital Corporation Disk drive capable of autonomously evaluating and adapting the frequency response of its servo control system
US6754025B1 (en) * 2000-01-12 2004-06-22 Maxtor Corporation Disk drive spindle motor speed and timing control
US6704156B1 (en) * 2000-01-31 2004-03-09 Maxtor Corporation Self-writing of servo patterns in a disk drive using a printed reference pattern
US7054088B2 (en) * 2002-01-21 2006-05-30 Fujitsu Limited Information recording and reproducing apparatus and method, and signal decoding circuit for performing timing recovery
US6989956B2 (en) * 2002-07-03 2006-01-24 Fujitsu Limited Head position control method, disk device, and servo track write method
US6980389B1 (en) * 2004-02-27 2005-12-27 Western Digital Technologies, Inc. Disk drive having a disk including a servo burst pattern having a normal servo burst field and a quadrature servo burst field
US7400561B2 (en) * 2004-12-14 2008-07-15 Lsi Corporation Single PLL demodulation of pre-formatted information embedded in optical recording medium

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8947811B1 (en) * 2013-03-12 2015-02-03 Western Digital Technologies, Inc. Disk drive using preamble and postamble to adjust asynchronous signal samples of payload

Also Published As

Publication number Publication date
US7405894B2 (en) 2008-07-29
US20060007571A1 (en) 2006-01-12

Similar Documents

Publication Publication Date Title
US6628606B1 (en) Coarse frequency offset estimator in orthogonal frequency division multiplexing receiver and method thereof
US9305582B2 (en) Systems and methods for fly-height control using servo data
JP5684187B2 (en) Method for obtaining frequency difference between input signal and reference frequency, discriminating apparatus, GPS receiver and computer program for executing this method
JP4583374B2 (en) Receiver
US7027486B2 (en) System and method for fast code phase and carrier frequency acquisition in GPS receiver
US20080278837A1 (en) Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive
JP4783481B1 (en) Ultrasonic measurement method and ultrasonic measurement apparatus
US8064552B2 (en) Adaptive correlation
US20020061118A1 (en) Electronic watermarking method and system
JP2004282759A (en) Synchronization method and apparatus for initial frequency in ofdm system
WO2000044120A1 (en) Synchronization of ofdm signals
JP2000068893A (en) Method and architecture to correct carrier offset and timing offset of spread code in direct sequence spread spectrum communication system
US9516616B2 (en) Method and apparatus for estimating frequency errors
CN110708267B (en) Method for determining frequency offset information estimation value
JP5101910B2 (en) Frequency estimation method, broadband frequency discriminator, and receiver for wireless position measurement
US9672833B2 (en) Sinusoidal interpolation across missing data
US5073907A (en) Digital phase-lock loop
US7590184B2 (en) Blind preamble detection for an orthogonal frequency division multiplexed sample stream
CN113237492B (en) Initial phase stability determination method, storage medium and electronic device
Davis et al. A new approach to single-tone frequency estimation via linear least squares curve fitting
US7830971B2 (en) Apparatus and method for detecting advanced GCL sequence in wireless communication system
CN114978828B (en) Frequency offset estimation method and device, electronic equipment and storage medium
EP1296493A2 (en) Symbol synchronisation in a multicarrier receiver
KR102133126B1 (en) Frequency estimation apparatus and method
Xu et al. Time Delay Estimation of Ultrasonic Forward and Reverse Signals Based on Phase Difference

Legal Events

Date Code Title Description
STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION

AS Assignment

Owner name: BANK OF AMERICA, N.A., AS COLLATERAL AGENT, NORTH CAROLINA

Free format text: PATENT SECURITY AGREEMENT;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:037806/0001

Effective date: 20160201

Owner name: BANK OF AMERICA, N.A., AS COLLATERAL AGENT, NORTH

Free format text: PATENT SECURITY AGREEMENT;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:037806/0001

Effective date: 20160201

AS Assignment

Owner name: AVAGO TECHNOLOGIES GENERAL IP (SINGAPORE) PTE. LTD., SINGAPORE

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:041706/0001

Effective date: 20170120

Owner name: AVAGO TECHNOLOGIES GENERAL IP (SINGAPORE) PTE. LTD

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BROADCOM CORPORATION;REEL/FRAME:041706/0001

Effective date: 20170120

AS Assignment

Owner name: BROADCOM CORPORATION, CALIFORNIA

Free format text: TERMINATION AND RELEASE OF SECURITY INTEREST IN PATENTS;ASSIGNOR:BANK OF AMERICA, N.A., AS COLLATERAL AGENT;REEL/FRAME:041712/0001

Effective date: 20170119