Receiver for data transmission system with nonlinearities

Abstract: System for transmitting a data signal at a symbol rate 1/T through a noisy dispersive channel to a data receiver; said channel introducing intersymbol interference and noise into the transmitted data signal; and said receiver estimating the most likely sequence of transmitted data symbols by keeping track of candidate data sequences that are recursively updated on the basis of likehood measures which are determined with the help of means for estimating hypothesized channel outputs in the absence of noise comprising one or more look-up tables, the possibility is provided to adapt each look-up table under the control of digits of said candidate data sequences, in response to an error signal that is representative for the difference of the channel output signal and the output signal of said look-up table.

Description: The invention relates to a system according to the preamble of claim 1. Such a system is know from an article by G.D. Forney, Jr. entitled "Maximum-Likelihood Sequence Estimations of Digital Sequence Estimation of digital Sequences in the Presence of Intersymbol Interference"", IEEE Trans. Inform. Theory, Vol. IT-18, No. 3, pp. 363-378, May 1972. In the forthcoming decades, digital optical storage is expected to find widespread for storage of computer data as well as of digitized audio and video signals. In systems of this type, binary information is stored as a sequence of this type, binary information is stored as a sequence of pits and lands in an optical medium. Imperfections in writing process may cause pits and lands to differ in shape or length. In the replay process, this asymmetry manifests itself in the form of nonlinear intersymbols interference (ISI). This disturbance comes on top of linear ISI that arises due, for example to limitations of optical resolution, and to noise generated in e.g. laser diodes and preamplifiers. Conventional reception techniques are, in general, only able to deal with the latter two imperfections. This even applies for the most powerful variety, whose representatives are commonly referred to as Viterbi detectors. Conventional Viterbi detectors form an estimate of the most likely transmitted data sequence, assuming that only linear ISI and noise are present. To this end, they maintain a list of candidate data sequences that are referred to as survivors. These survivors are recursively extended, and a selection process takes place on the basis of likelihood measures that are calculated for each survivor by comparing the actual channel output signal with a hypothesized output signal that would result if noise were absent and the concerned survivor would have been transmitted. The aforementioned means to form these hypothesized channel output signals conventionally consist of linear weighing networks that operate on a given number of the most recent symbols of concerned survivors. Basic elements of this detection process are described in more detail, for example, in the aforementioned article by Forney. The conventional Viterbi detector developed in this article has the disadvantage that its complexity grows rapidly to unmanageable proportions as channel memory lengths increase. To overcome this problem, various simplified versions of the conventional Viterbi detector have been developed in recent years, as exemplified, for example, by an article by J.W.M. Bergmans, S.A. Rajput and F.A. M. van de Laar entitled "On the Use of Decision Feedback for Simplifying the Viterbi Detector", Philips J. Res., Vol. 42, No. 4, pp. 399-428, 1987. Partly as a consequence of these efforts towards simplification, conventional Viterbi detectors now find application in such areas as voiceband modems and digital magnetic recording. To date, however, nonlinear ISI has discouraged their use in optical storage. An article by M.F. Mesiya, P.J. McLane, and L.L. Campbell entitled "Maximum Likelihood Sequences Transmitted over Nonlinear Channels", IEEE Trans. Commun., Vol COM-25, No. 7, pp. 633-643, July 1977, puts forth a novel type of Viterbi detector that distinguishes itself from conventional Viterbi detectors in that it can handle nonlinear ISI. Unfortunately this ability comes, in general, at the cost of a greatly increased complexity. US 4,015,238 discloses a Viterbi decoder for decoding a signal being encoded using a convolutional code. Said decoder is only arranged for coping with noise, and is not intended for dealing with intersymbol interference. It is the object of the present invention to provide a novel category of Viterbi detectors with the ability to handle a mixture of linear and nonlinear ISI and noise, and with a complexity comparable to that of Viterbi detectors of the abovementioned prior art. To this end the receiver is arranged according to the characterizing part of claim 1. As look-up tables are ideally suited to store nonlinear input-output relations, they can account for any nonlinearity of the channel. Furthermore, because they can account for any nonlinearity of the channel. Furthermore, because they all operate on digital data symbols, these look-up tables are often more conveniently implemented in digital technology than the linear weighing networks in receivers of the above-mentioned prior art. In this way the ability to handle nonlinear ISI is generally accompanied by a decreased receiver complexity. An especially simple version of the receiver according to the invention has only two candidate data sequences, which are recursively updated on the basis of a likelihood measure that is representative for the difference of a function of the likelihoods of both candidate data sequences, and that is determined with the help of means for estimating hypothesized channel outputs in the absence of noise. To facilitate its application at high data rates, this receiver is further characterized in that said likelihood measure is determined by selection among precomputed candidate values of said likelihood measure. A special version of the receiver according to the invention in which all look-up tables have only one entry is characterized in that said look-up tables take the form of registers that store hypothesized channel output symbols in the absence of noise. In practice, the precise characteristics of a storage channel usually vary dynamically around a nominal average as a result of e.g. mechanical vibrations and tracking or defocussing errors. As a rule, the performance of Viterbi detectors degrades rapidly when the actual characteristics of the channel come to differ from the ones reflected in said means for estimating hypothesized channel output signals in the absence of noise. The detrimental effects of this channel-receiver mismatch are illustrated, for example, in an article by K.A. Schouhamer Immink entitled "Coding Methods for High-Density Optical Recording", Philips J. Res., Vol. 41, No. 4, pp. 410-430, 1986. It is a further object of the invention to reduce performance degradation due to said channel-receiver mismatch. As a first means to this end, the receiver is further characterized in that each look-up table is adapted under the control of digits of said candidate data sequences, in response to an error signal that is representative for the difference of the channel output signal and the output signal of said look-up table. As a second means to this end, the receiver according to the invention is further characterized in that each look-up table is adapted in response to an error signal that is representative for the difference of a delayed version of the channel output signal and the output signal of said look-up table when addressed by one or more delayed digits of said candidate data sequences. By making the look-up tables adaptive, they attain the ability to track variations of the channel characteristics. This greatly reduces said channel-receiver mismatch. For the special case that look-up tables take the form of registers that store hypothesized channel output symbols in the absence of noise, an adaptive version of the receiver according to the invention is further characterized in that each register takes the form of a digital counter that is adapted, under the control of one or more delayed digits of said candidate data sequences, in response to an error signal that is representative for the difference of a delayed version of the channel output signal and the contents of said counter. An especially simple version of a receiver according to the invention is further characterized in that each likelihood measure is representative for an accumulated version of a function that essentially equals the modulus of the difference of the actual channel output signal and a hypothesized channel output signal in the absence of noise. To facilitate its application at high data rates, this receiver is further characterized in that said function is determined by selection among precomputed candidate values of said function. Backgrounds and embodiments of the invention and their advantages will now be further explained with reference to the drawing, in which Fig. 1 shows a functional discrete-time model of a system for transmitting data symbols a k at a symbol rate 1/T through a noisy dispersive channel CHN to a data receiver REC. Fig. 2 shows a conceptual model of the computations that are performed for any survivor in a receiver of the above-mentioned prior art to guide the selection process; Fig. 3 shows a conceptual model of the computations that are performed for any survivor in a receiver according to the invention to guide the selection process; Fig. 4 shows an adaptive version of the conceptual model of Fig. 3; Fig. 5 shows an adaptive version of the conceptual model of Fig. 3 in which adaptation is based on delayed digits of the survivor. Fig. 6 shows a model of a 2-state Viterbi detector with linear feedback according to the above-mentioned prior art; Fig. 7 shows a conceptual model of an adaptive 2-state Viterbi detector with nonlinear feedback according to the invention; Fig. 8 shows a conceptual model of an adaptive precomputation unit for a receiver according to the invention; Fig. 9 shows a conceptual model of an adaptive 2-state Viterbi detector with nonlinear feedback according to the invention that uses adaptive pre-computation units according to Fig. 8; Fig. 10 shows bit error characteristics that were obtained by simulation for a conventional receiver according to Fig. 6 and a receiver according to the invention that conforms to Figs. 8 and 9; Fig. 11 illustrates the nonlinearity mechanism in the system that underlies the simulation results of Fig. 10; Fig. 12 shows the transfer characteristics of the linear part of the recording channel that underlies the simulation results of Fig. 10. In all figures, corresponding elements are denoted by the same symbols. In the following description a discrete-time modelling of the transmission system and the arrangement is used, as the general notion of the invention can be presented in the most simple way with reference to such a modelling. This does not lead to a loss of generality as the present modelling can be derived unambiguously from the parameters of the actual continuous-time system as described, for example, in the book entitled "Digital Communications", by J.G. Proakis, McGraw-Hill, New York, 1983, Chapter 6, and especially Section 6.3, pp. 351-357. Fig. 1 shows a functional discrete-time model of a system for transmitting data symbols a k at a symbol rate 1/T through a noisy dispersive channel CHN to a data receiver REC. To simplify the exposition, it will henceforth be assumed that the transmitted data signal a k is binary with a k ∈ {-1, +1}. This assumption is not meant to be restrictive. With self-evident modifications, the invention is equally applicable to multilevel or complex-valued data signals, as encountered in e.g. digital voiceband communication systems. The channel CHN of Fig. 1 models the cascade of the actual continuous-time channel, a possible receiving filter and/or equalizer, and a synchronous sampling operation at the data rate 1/T. The discrete-time output signal r k of channel CHN can be described as r k = f ( a k ) + n k , where n k is a white Gaussian noise signal and f(.) is a deterministic function of a data vector a k = [a k-M , a k-M+1 ,....,a k ] T , where ' T ' denotes transposition. The nonnegative integer M is referred to as the memory length of the channel. In the absence of nonlinear ISI, f( a k ) assumes the linear form f( a k ) = f T a k , where f = [f₀,...,F M ] T is a vector of length M+1 whose components specify the impulse response of the channel. The receiver REC in Fig. 1 operates on r k in order to produce decisions â k-D about a delayed version a k-D of a k , where D is a nonnegative integer that is referred to as the detection delay. Conventional receivers REC are, in general, not able to handle nonlinear ISI. This applies even for the most powerful variety of conventional receivers, whose representatives are generally known as Viterbi detectors. These detectors form an estimate of the most likely transmitted data sequence, assuming that only linear ISI and noise are present. To this end, they maintain at any instant k-1 a list with a given number N of candidate data vectors s i k-1 = â i k-D ,....,â i k-1 T for all i ∈ 0,...,N-1 that are referred to as survivors. In the course of the detection process, this list is recursively updated on the basis of likelihood measures that are calculated for all survivors. Central elements of this process are described, for example, in said papers by Forney and Bergmans et al.. To clarify the process in more detail, Fig. 2 depicts a basic model of the likelihood calculations that are associated to any survivor s i / k-1 in a receiver according to the above-mentioned prior art. In Fig. 2, a measure of accumulated likelihood j i / k-1is associated to the survivor s i / k-1. In accordance with the usual nomenclature, this measure of accumulated likelihood will henceforth be referred to as a metric for the sake of compactness. As a first step in the recursion that is carried out at moment k, the oldest digit â i / k-D of survivor s i / k-1is neglected, and the two possible digits -1 and +1 are appended in order to obtain two new candidate survivors s ij k = â i k-D+1 ,...â i k-1 ,I j for j=0 and 1, respectively, where the index function I(j) is defined according to I(0)=-1 and I(1)=+1. For the sake of completeness it is noted that this index function would assume integer values in 0,..., |A|-1 for a data alphabet that consists of |A| rather than 2 distinct symbols. Along this line it is straightforward to apply the invention to multilevel or complex-valued data alphabets. As explained earlier, the forthcoming exposition is entirely cast in terms of a binary alphabet, for which |A| = 2, solely for the purpose of simplifying the presentation of the invention as much as possible. Associated to the extended survivors s ij / k are metrics J ij k =J i k-1 +G r k - f T â ij k , where â ij k = â i k-M ,....,â i k-1 ,I j is a vector whose components are the M+1 most recent components of s ij / k, while G(.) is a deterministic function for which the choice G(x)=x² for all x ∈ R is common. The components f T â ij / k are generated by linear weighing networks LW ij that operate on the M most recent digits of s ij / k, and can be recognized as hypothesized channel output samples that would result on moment k if noise ware absent and s ij / k were transmitted. In the usual case that G(x)=x² for all x ∈ R, the metrics J ij / k can be interpreted as accumulated Euclidean distances between the actual channel output signal r k and hypothesized channel output signals f T â ij / k. As time proceeds, the detector seeks to minimize this distance across all considered survivors. To update its list of survivors, the detector compares the metrics J ij / k of the extended survivors for all i E {0,...N-1} and j E {0,1}, and makes a selection on this basis. The details of this selection depend on the precise type of Viterbi detector, and will not be described here in further detail as they are immaterial to the invention. Because of their linear form, the outputs of the linear weighing networks LW ij can only serve as hypothesized channel outputs in the absence of noise for channels CHN that do not introduce nonlinear ISI. For this reason Viterbi detectors that conform to Fig. 2 are intrinsically unable to handle nonlinear ISI. Central to the invention is the observation that the linear weighing networks LW ij operate on vectors â ij / k with a finite number M+1 of binary data symbols. This makes it possible to replace said networks by look-up tables in which the possible hypothesized channel output signals in the absence of noise are all stored. This is shown in Fig. 3. Fig. 3 is identical to Fig. 2 except for look-up tables LUT ij that replace the linear weighing networks LW ij . Each table is addressed by a total of M+1 binary data symbols and must therefore contain a total of 2 M+1 entries. Even for values of M as large as e.g. 10 this poses no instrumentational problems when use is made of currently available random access memories. In fact, since possible outputs are directly looked up rather than calculated, the configuration of Fig. 3 will generally be easier to implement than that of Fig. 2. Furthermore, since look-up tables can store fully arbitrary input-output relations, the outputs h( â ij / k) of LUT ij can serve as hypothesized channel outputs in the absence of noise by choosing h( a k ) = f( a k ) for all a k . In this way nonlinear ISI can be fully dealt with without increasing the complexity of the receiver. This compares favourably with a novel Viterbi detector that is described in said article by Mesiya et al.. This detector also distinguishes itself from Viterbi detectors according to the above-described prior art in that it can handle nonlinear ISI. However, for the Viterbi detector of Mesiya et al. the ability to handle nonlinear ISI comes, in general, at the cost of a greatly increased complexity. For example, in Fig. 3 of said article it is illustrated that for M=3 the detector operates on a total of 4 input signals, as opposed to the single input signal r k for the receiver according to the invention. Furthermore, for each of these 4 input signals various summations need to be performed per symbol interval T and survivor (see eq. (26) of said article by Mesiya et al. and the accompanying explanations), as opposted to the single look-up operation in the receiver according to the invention. In Fig. 3, the tables LUT i1 and LUT i0 are addressed by vectors â i,1 / k and â i0 / k whose most recent bits â i1 / k = +1 and â i0 / k=-1 are known a priori. This enables both tables to be halved in size, with LUT i1 and LUT i0 storing only the parts of f( a k ) for which a k =+1 and -1, respectively. In various types of Viterbi detector, including those described in said article by Bergmans et al., several of the most recent bits of â ij / k are known a priori. This enables further reductions in the size of each table. For certain Viterbi detectors, including the one described in said article by Forney, the complete vectors â ij / k are even known a priori. This allows each table to degenerate into a single register that stores the corresponding hypothesized channel output in the absence of noise. An adaptive version of a detector according to the invention that used registers in the form of adaptive digital counters will be described later. Prior knowledge about the channel characteristics may be used to identify the function f(.). The tables LUT ij can then be filled with appropriate values on the basis of condition (8). Unfortunately, the precise characteristics of a storage channel usually vary dynamically as a result of e.g. mechanical vibrations and tracking or defocussing errors. The function f(.) then fluctuates in time, and it becomes impossible for a time-invariant function h(.) as stored in the look-up tables LUT ij of Fig. 3 to fully match f(.). For the particular case of linear functions f(.) and h(.), the detrimental effects of such a channel-receiver mismatch were studied, for example, in said article by Schouhamer Immink. In this article, it is illustrated that even small mismatches may result in serious performance degradations, especially at high information densities. To avoid such degradations, it is desirable for the function h(.) to track any variations of f(.). This is possible by making the look-up tables LUT ij adaptive, as illustrated in Fig. 4. In Fig. 4, the function h(.) as stored in the look-up tables LUT ij is updated recursively according to h'( â ij k ) = h( â ij k ) + µ·d ij k ·e ij k for all i E {0,...,N-1} and j ∈ {0,1}. Compared to the old table entries h( a ij / k), the new entries h'( a ij / k) are ideally improved estimates of f( a ij / k). The error signals e ij k = r k - h( â ij k ) indicate how well the hypothesized channel output signals in the absence of noise resemble the actual channel output signal, and the recursion of (9) seeks to minimize this difference iteratively. A more detailed description of the so-called LMS adaptation algorithm that forms the basis for the recursion of (9) can be found, for example, in an article by P.J. van Gerwen, N.A.M. Verhoecks and T.A.C.M. Claasen entitled "Design Considerations for a 144 kbit/s Digital Transmission Unit for the Local Telephone Network", IEEE J. Selected Areas in Commun., Vol. SAC-2, No. 2, pp. 314-323, 1984. This article discusses application of the LMS algorithm to table look-up filters, and is as such closely related to the application presently being discussed. For proper operation of the LMS algorithm, the data vector a k which by (1) underlies r k should coincide with the data vector â ij / k of the table that is being updated. Among all possible data vectors â ij / k, the one with greatest accumulated likelihood is most likely to satisfy this prerequisite. For this reason, the selector signals d ij / k of expression (9) are chosen according to d ij k = 1 if J ij k < J nl k for all n ∈ 0,...,N-1 , and 1 ∈ {0,1}, 0, Else. Thus at any instant k at most one of the tables is updated. It can be observed that the selector signals d ij / k of (11) are entirely based on information that is generated as an integral part of the detection process. For this reason they can be generated with a minimum of extra hardware. The adaptation constant µ in (9) enables a tradeoff between speed of convergence of the tables and steady state excess mean-square error. To simplify a digital implementation of the algorithm, µ is usually chosen to be of the form 2 -W for some positive integer W, so that the multiplication by µ in (9) amounts to a shift over W bit positions. For the sake of compactness, these and other aspects of the LMS algorithm and its implementation are not described here in further detail, as they are well documented in the literature, see e.g. said article by van Gerwen et at.. A detailed description of an adaptive receiver according to the invention and based on a simplified version of the LMS algorithm will be given later. In practice it is advisable to initialize the table entries h( â ij / k) in accordance with the average characteristics of the channel, provided that these are known a priori. In this way the adaptation algorithm only has to bridge the distance between the average and the actual characteristics, thereby allowing comparatively rapid convergence. For the situation that tables are initialized randomly, simulations reveal that convergence periods can often be reduced significantly by adding an appropriate amount of noise to the receiver input signal during the initial phase of adaptation. A disadvantage of the configuration of Fig. 4 is that very recent estimated data symbols play a role in the adaptation process. By nature of Viterbi detection, these symbols are less reliable estimated of the transmitted data signal than the older digits that also form part of the maintained survivors. More specifically, even if a given survivor s ij / k has greatest current likelihood, its most recent digits (e.g. â ij / k and â ij / k-1) may not coincide with the corresponding transmitted digits. Especially for functions f(.) with a weak dependence on these most recent transmitted digits this may in fact occur quite frequently. By (9) and (10) this would equally often cause erroneous table entries to be updated, a problem that may hamper or even preclude convergence of the table contents to the proper values. To overcome this problem it is necessary to base adaptation on more reliable and hence delayed estimates of the transmitted data signal. A natural possibility to this end is outlined in Fig. 5. When the six switches SW 0 / 0,...SW 1 / 2 of Fig. 5 are in the position "detect", detection proceeds exactly as in Fig. 4, with the look-up tables LUT ij addressed by the estimated data vectors â ij / k. For adaptation the switches are placed in the position "adapt". In this case a delayed data vector b i k = [â i k-M-P ,...,â i k-P ] addresses every look-up table LUT ij . If the delay P is taken sufficiently large, then the digits A i / k-M-P,..., â i / k-P are for any i likely to be reliable estimates of the actually transmitted data symbols a k-M-P ,...,a k-P . Greatest reliability is obtained by selecting the maximum value P=D-M, as shown in Fig. 5, so that the digits â i / k-M-P,...,a i / k-P of b and i / kare the M+1 oldest maintained digits â i / k-D,...,â i / k-D+M of survivor s and i / k-1. To compensate for the data delay of P symbol intervals, the received signal r k is also delayed over P symbol intervals in order to form a delayed error signal which is used to update the look-up table LUT ij according to the LMS logarithm As the data estimates used in this adaptation process are all relatively reliable, it is unnecessary in (14) to use selector signals to condition adaptation on current or past likelihood measures. A disadvantage of the configuration of Fig. 5 is that each table is read out twice per symbol interval for calculation of error signals that play a role in detection and adaptation, respectively. Relative to Fig. 4, where these two functions are combined, this lowers the largest attainable data throughput. To overcome this problem it is possible to base adaptation on delayed versions of the error signals that were calculated for detection P symbol intervals earlier. This also makes it unnecessary to delay the received signal r k . A simplified version of this possibility will be described later. To exemplify the preceding notions, two versions of a two-state Viterbi detector with nonlinear feedback according to the invention will now be developed. To facilitate the explanation, a related Viterbi detector of prior art will be described first. Fig. 6 shows a conceptual model of a two-state Viterbi detector with linear feedback as described in the aforementioned article by Bergmans et al.. This detector has two survivors s 0 / k-1 and s 1 / k-1 conforming to eq. (4), with associated metrics J i / k for i=0,1. The four extended survivors s ij / k for i,j ∈ {0,1} are defined as in (5) and have metrics J ij / k according to (6). Four linear weighing networks LW ij with i,j ∈ {0,1} calculate the four possible weighted sums f T â ij / k of eq. (6). As explained earlier, the vecot f specifies the impulse response of the channel. This may be achieved, for example, with the help of adaptive techniques, as described in the aforementioned book by Proakis, chapter 6, pp. 410-412. Details of these techniques as applied in receivers of prior art are immaterial to the invention and therefore not discussed or shown here. Among the extended survivors s 0j / k and s 1j / k, a Compare.Select unit CS j selects the new survivor s j / k with associated metric J j / k for moment k according to the rule This rule is applied for both j=0 and j=1. It can be noted from (15) that the most recent digits â j / k of the new survivors s j / k are always equal to the index function I(j), i.e. â 0 / k=-1 and â 1 / k=+1 for all k. Because of the recursive nature of the detection process, this means that the most recent digits a i / k-1 of the old survivors s i / k-1 must also equal â 0 / k-1=-1 and â 1 / k-1 irrespective of k. Thus, for the extended survivors s ij / k, the two most recent digits â ij / k-1 and â ij / k are both known a priori to equal I(i) and I(j), respectively. For any Viterbi detector, it is desirable to have a detection delay D much greater than the channel memory length M, as explained, for example, in the aforementioned article by Forney. In this case, the oldest digits â k-D and â 0 / k-D are both comparatively reliable estimates of the transmitted digit a k-D . In Fig. 4, â 1 / k-D is arbitrarily chosen to be the detector output â k-D . This choice is only illustrative and is not meant to be restrictive, for it is clear that a different choice, such as â k-D =â 0 / k-D, might be equally appropriate. Further backgrounds and clarifications of the detection process of Fig. 6 are not expounded here as they are described in detail in said article by Bergmans et al.. A disadvantage of the detector of Fig. 6 is that the metric values j i / k are, by (6) and (15), a non-decreasing function of time in the usual case that the function G(.) is nonnegative definite. This may cause problems of overflow in a digital implementation of the detector. From (15) it can be noted that only differences between metrics play a role in the selection of new survivors. This observation may be used to re-normalize metric values in such a way that they are no longer a non-decreasing function of time. To this end, the modified metrics Q k , Q 0 / k and Q 1 / k are defined as Q k = J 1 k -J 0 k , and Q i k = J i k -J 0 k-1 for i ∈ {0,1} and all k. By making use of (6), it is straightforward to reformulate (15) in terms of these modified metrics as where the error signals e ij / k are defined as e ij k = r k - f T a ij k for i,j E {0,1}. Furthermore, from (16) and (17) it is seen that Q k can be directly calculated from Q 0 / k and Q 1 / k according to Q k = Q 1 k -Q 0 k . Thus the complete selection process can be recast in terms of only three difference metrics Q k , Q 0 / k and Q 1 / k. As these metrics all fluctuate around zero in value, problems of overflow in a digital receiver implementation can be avoided by a proper choice of digital word-lengths. An adaptive receiver according to the invention that incorporates these difference metrics is shown in Fig. 7. For the sake of completeness it is mentioned here that the use of these difference metrics is in itself not new, witness e.g. an article by R.W. Wood and D.A. Petersen entitled "Viterbi Detection of Class IV Partial Response on a Magnetic Recording Channel", IEEE Trans. Commun., Vol. COM-34, No. 5, pp. 454-461, May 1986, and in particular the passage that surrounds expressions (11) and (12) of this article. In the receiver of Fig. 7, two compare-select units CS⁰ and CS¹ are used to control the selection process of eq. (18) for j=0 and j=1, respectively. Based on input signals Q k +G[e 1j / k] and G[e 0j / k], compare-select unit CS j produces an output signal Q j / k according to eq. (18), and a selector signal d j / k according to If Q k-1 +G e 1j k <G e 0j k Then d j k :=1 Else d j k :=0; By comparing (18) and (21), it is seen that this selector signal d j / k can be used to control the selection of survivors according to the rule If d j k =1 Then s j k :=s 1j k Else s j k :=s 0j k ; which applies for both j=0 and j=1. To implement this selection process, two shift registers SR⁰ and SR¹ store the digits [â 0 / k-D,...â 0 / k-2] and [â 1 / k-D,...,â 1 / k-2] of the survivors s 0 / k-1 and s 1 / k-1, respectively. As explained earlier, the most recent digits of these survivors are a priori known to be â 0 / k-1=-1 and â 1 / k-1=1, and these digits are therefore represented in the form of fixed logical levels +1 and -1 that are connected to the inputs of both shift registers. For shift register SR⁰, a selector signal d 0 / k=0 indicates by (22) that the new survivor s 0 / k should be s 00 / k-1. By (7), s 00 / k is just a shifted version of s 0 / k-1, with the oldest digit â 0 / k-D removed and a most recent digit â 0 / k=-1 appended. Thus the selection s 0 / k:= s 00 / k can be realized with a SHIFT LEFT operation on shift register SR⁰. Similarly, a selector signal d 0 / k=1 indicates that the new survivor s 0 / k should be s 10 / k. By (7), s 10 / k is just a shifted version of s 1 / k-1, with the oldest digit â 1 / k-D removed and a most recent digit â 0 / k=-1 appended. Thus the selection s 0 / k:= s 10 / k can be realized with a PARALLEL LOAD operation in which shift register SR⁰ is loaded with a shifted version of the contents of shift register SR¹ and the most recent digit â 1 / k-1=+1. This is symbolically indicated in Fig. 7 by means of the skew arrows that run from SR¹ and â 1 / k-1=1 to SR⁰. For shift register SR¹, the selector signal d 1 / k similarly indicates either a SHIFT LEFT operation for d 1 / k=1, or a LOAD PARALLEL operation from shift register SR⁰ and digit â 0 / k-1=-1 for a 1 / k=-1. For the sake of completeness it is necessary to mention a potential problem that may occur in a direct implementation of the shift register configuration of Fig. 7 as a consequence of poor matching of the propagation delays of both shift registers. If shift register SR¹ happens to have a significantly smaller propagation delay than shift register SR⁰, then a PARALLEL LOAD operation on SR⁰ may cause one or more digits of the new survivor s 1 / k rather than the desired ones of the old survivor SR 0 / k-1 to be loaded into SR¹. Similarly, if the propagation delay of shift register SR⁰ is significantly smaller than that of SR¹, then a LOAD PARALLEL operation on shift register SR¹ may cause one or more digits of the new survivor s 0 / k rather than the desired ones of the old survivor s 0 / k-1 to be loaded into SR¹. Both possibilities are clearly undesirable. To avoid this problem, it is possible in a practical implementation of the receiver of Fig. 7 to choose shift registers SR¹ and SR⁰ with well-matched propagation delays, or to latch the cross-couplings between both shift registers. As the diagram of Fig. 7 is merely meant to provide a conceptual model of a receiver according to the invention, possibilities to avoid this implementation-level problem will not be elucidated here in any further detail. In Fig. 7, the look-up tables LUT ij with i,j ∈ {0,1} are addressed by the digits [â i / k-M,..., â i / k-2] of survivor s i / k-1. By (5) and (7), these digits coincide with the corresponding digits of the address vectors a ij / k. The remaining two digits are known a priori to be â ij / k-1=I(i), as explained for the receiver of Fig. 6. As argued earlier, it is unnecessary to add such a priori known digits to the address vector of look-up tables. This allows each table to be 4 times smaller in size than for an address vector of the "full" length M+1. Mathematically speaking, this four-fold reduction is possible because the four look-up tables LUT ij cover disjunct portions of the domain of h( a k ), viz. those portions for which a k-1 =I(i) and a k =I(j). The adaption mechanism for the look-up tables LUT ij is identical to that of Fig. 4 and is therefore not explained in further detail. At any instant k, the selector signals d ij / k are such that only the table that corresponds to the most likely extended survivor is updated. These selector signals are produced by a selector unit SEL that operates, for example, on the signals Q k , d 0 / k and d 1 / k according to the following truth table: Q k d 0 / k d 1 / k d 00 / k d 01 / k d 10 / k d 11 / k <0 0,1 1 0 0 0 1 <0 0,1 0 0 0 1 0 >0 1 0,1 0 1 0 0 >0 0 0,1 1 0 0 0 To explain this table, it can be noted from (16) and the accompanying explanation that any positive value of Q k indicates that the new survivor s 0 / k is more likely than its counterpart s 1 / k, and vice versa for a negative value of Q k . Thus, for positive Q k it is necessary to distinguish between the two extended survivors s 00 / k and s 01 / k that underly s 0 / k. The signal d 0 / k can be used to this end, as it specifies, by eq. (22), exactly which of these two extended survivors forms s 0 / k. Similarly, for Q k <0 the signal d 1 / kspecifies which of the two selector signals d 01 / k and d 11 / k is to be 1, while the other two selector signals are zero. An attractive choice for the function G(.) in Fig. 7 is G(x)= |x| for all x ∈ R, since this function can be easily realized or approximated with digital circuitry. Simulations for the receiver of Fig. 7 reveal that in many cases this choice yields performances that are essentially equivalent to those for the usual function G(x)=x² for all x ∈ R, which is comparatively difficult to realize or approximate with digital circuitry. With respect to the receiver of Fig. 7, it remains to mention that the new value Q k is determined from the signals Q 1 / k and Q 0 / k according to eq. (20) by means of a summator, while a delay unit stores Q k for use during the next symbol interval. Furthermore, the oldest digit â 1 / k-D serves as the output â k-D of the receiver, as in Fig. 6. For the sake of brevity, further aspects of the receiver are not elaborated here as they are either sufficiently self-evident or sufficiently similar to aspects that were discussed before. The receiver of Fig. 7 is attractive in that it combines a complexity no greater than that of its linear counterpart of Fig. 6 with the ability to handle any form of linear or nonlinear ISI. Together with the receiver of Fig. 6, it shares the disadvantage that implementation may become difficult at very high data rates, as encountered in e.g. digital storage of video signals. One cause of this difficulty is that the formation of the signals G(e ij / k) in Fig. 7 requires a table look-up operation, a subtraction and application of the function G, which together may require more time than is permissible. A technique to overcome this problem will now be described. To explain the technique, it is noted from Fig. 7 that the signals G(e ij / k) depend only on r k and the digits â i / k-M,...â i / k-2 of survivors s i / k-1. Since r k is known, this leaves a total of 2 M+1-2 =2 M-1 possible values for G(e ij / k). By precomputing all these values it becomes possible to make a direct selection of G(e ij / k) on the basis of the actual digits â i / k-M,...,â i / k-2, which requires substantially less time than the sequence of operations just mentioned. For the simplest possible case, viz. M=2 and 2 M-1 =2, an adaptive precomputation unit APU ij that makes use of this possibility is depicted in Fig. 8. In the system of Fig. 8, two digital up/down counters C 0ij and C 1ij replace the w M-1 =2 table entries of table LUT ij in Fig. 7. More specifically, these counters account for extended survivors with data vectors â ij / k =[-1,I(i),I(j)] T and [+1,I(i),I(j)] T by storing nonlinear function values h 0ij =h([-1,I(i), I(j)] T ) and h 1ij =h([1,I(i),I(j)] T ), respectively. Details about the use of digital circuits to represent and process discrete-time signals are not discussed here as they are well described in the literature, see e.g. the book "Theory and Application of Digital Signal Processing" by L.R. Rabiner and B. Gold, Prentice-hall, N.J., 1975. Two adders serve to form error signals e 0ij / k and e 1ij / k by subtracting h 0ij and h 1ij from r k . Subsequent application of the function G(.) yields signals G(e 0ij / k) and G(e 0ij / k), one of which corresponds with the signal G(e ij / k) that is to be formed. To effectuate the proper selection, the digit â i / k-2 of survivor s i / k controls s switch SW g according to If â i k-2 =1 Then G e ij k :=G e 1ij k Else G e ij k :=G e 0ij k . Since G(e 0ij / k) and G(e 1ij / k) are both precomputed, the only delay incurred in generating G(e ij / k) arises from the selection (23) and can be made very small. In the circuit of Fig. 8, the digit â i / k-2 is used to deselect one of the two possible data vectors â i / k=[-1,I(i),I(j)] T and â i / k=[+1,I(i),I(j)] T . This feedback operation is the counterpart of the linear feedback operation that takes place in a more implicit manner in the conventional receiver of Fig. 7, as explained, for example, in the aforementioned article by Bergmans et al.. As a consequence of this feedback operation, only 4 out of the 8 possible vectors â i / k remain to be considered in the total of 4 adaptive precomputation units APU⁰⁰,...,APU¹¹, as opposed to the 8 vectors that would have to be considered in a Viterbi detector without feedback. Conceptually, this simplification may be justified by noting that the digit â i / k-2 that is fed back is the oldest and thus most reliable digit of the data vector â i / k under consideration. Thus the probability that an erroneous selection takes place in the circuit of Fig. 8 is relatively small. Properties and conceptual backgrounds of Viterbi detectors with and without feedback are not elaborated here in further detail as they are well described in the literature, see e.g. the aforementioned article by Bergmans et al.. The configuration of Fig. 8 includes a simplified version of the mechanism of Fig. 5 for adaptation of the counters C 0ij and C 1ij on the basis of delayed digits â i / k-M-P,...,â i / k-P. As explained earlier, the mechanism of Fig. 5 is preferable over the one of Figs. 4 and Fig. 7 in that it lowers convergence problems for functions f( a k ) with a weak dependence on the most recent digits of a k , such as a k and a k-1 . Said simplification stems from a sign operation that is performed on the error signals e 0ij / kand e 1ij / k to obtain one-bit error signals that are conveniently handled with digital circuitry. A switch SW e with a feedback function similar to that of SW g is controlled by â i / k-2 to obtain the one-bit and undelayed counterpart sgn of the error signal of eq. (13). This signal sgn is applied to a binary shift register that introduces a delay of P symbol intervals T. The delayed error signal sgn serves to update the contents of the counters according to the sign algorithm Details about this simplified version of the LMS algorithm can be found, for example, in an article by N. Holte and S. Stueflotten entitled "A New Digital Echo Canceller for Two-Wire Subscriber Lines", IEEE Trans. Commun., Vol. COM-29, pp. 1573-1581, Nov. 1981. This article discusses the application of the Sign-algorithm to adaptive table look-up filters, and is as such closely related to the application presently being discussed. In (24), d 1ij / k is a binary selector signal to be described, not to be confused with the signal d ij / k of Fig. 7. Furthermore, q is the quantization step size that corresponds to an increment or decrement of counter C 1ij by one unit. By an appropriate finite word-length represnetation of the quantities in Fig. 8, it is possible to select a suitably small value for q. Details of this finite word-length representation are not elaborated here as they are well described in the literature, see e.g. the aforementioned book by Rabiner and Gold. Because of the delay over P symbol intervals T, the delayed error signal sgn is by (1) a function of the delayed data vector a k-P . As explained in the aforementioned article by Holte and Stueflotten, for a proper operation of the sign algorithm of (24), only that counter C 1ij should be updated for which [I(1),I(i),I(j)] T = a k-P . When P is taken sufficiently large, the vectors b and 0 / k and b and 1 / k of (11) are both comparatively reliable estimates of a k-P , and either of them can thus be used for the formation of the selector signal d 1ij / k according to the rule IF â n k-P-2 =I 1 AND â n k-P-1 =I i AND â n k-P =I j THEN d 1ij k = 1 ELSE d 1ij k =0 which is applied for either n=0 or n=1 and for all 8 possible combinations (1,i,j). Alternative rules in which both vectors b and 0 / k and b and 1 / k are used for the formation of essentially equivalent selector signal signals are, of course, equally suitable but are not described here for the sake of brevity. In the configuration of Fig. 8, the signals d 1ij / k and sgn are connected to the COUNT ENABLE and UP/DOWN inputs of counters C 1ij in order to realize the iteration of (24). Depending on such implementation details as the type of counters used, slight modifications of the configuration of Fig. 8 may be needed to realize the iteration of eq. (24) in a convenient manner. In this respect Fig. 8 is merely illustrative, and is not meant to restrict in any sense the use of the sign-algorithm as described above. For the sake of completeness it is necessary to mention that the configuration of Fig. 8 can be easily modified for use of the LMS rather than sign-algorithm by omitting the Sign-operations in Fig. 1. The counters of Fig. 8 should then be replaced by digital accumulators that store h 1ij / k and can be updated in steps q. that may assume a multitude of sizes. Intermediate forms of the LMS and sign algorithms arise when such an accumulator is used in combination with a multi-bit quantizer instead of the sign operation in Fig. 8. Also, u or q can be variable rather than fixed. For example, for rapid convergence it is attractive to start adaptation with a relatively large value of u or q. Subsequently, u or q may be decreased gradually or step-wise to a value that is appropriate for small steady-state adaptation errors. Thus, with respect to the precise adaptation algorithm used, Fig. 8 is meant to be illustrative rather than restrictive. It is straightforward to generalize the configuration of Fig. 8 to channel memory lengths M greater than 2. In this case, the M-1 digits â i / k-M,..., â i / k-2 give a total of 2 M-1 possible values of G(e ij / k). To pre-compute these values, 2 M-1 adaptive counters with accompanying circuitry are needed, and the switches SW g and SW e are then used to select the signals G(e ij / k) and Sgn under control of the M-1 digits â i / k-M,..., â i / k-2. Further details of this generalization are not discussed here as they should be sufficiently self-evident after the foregoing explanations. Application of adaptive precomputation units like the one of Fig. 8 facilitates the attainment of high data rates. Fig. 9 depicts a model of a two-state Viterbi detector according to the invention in which the precomputation units of Fig. 8 are applied. The detector of Fig. 9 rather distinguishes itself from the one of Fig. 7 in that it employs a faster method of calculating Q k . In Fig. 7, calculation of Q k can not start before the selection process in compare/select units CS⁰ and CS¹ is completed. In Fig. 9, on the other hand, these actions occur largely in parallel. To explain this parallelism, it is noted from expression (18) that Q j k =Min ( G e 0j k ,Q k-1 +G e 1j k ) for both j=0 and j=1. With (19) this means that Q k =Min ( G e 01 k ,Q k-1 +G e 11 k )-Min( G e 00 k ,Q k-1 +G e 10 k ). Hence Q k assumes one out of four possible values G e 01 k -G e 00 k , G e 01 k - ( Q k-1 +G e 10 k , Q k-1 +G e 11 k -G e 00 k ), and Q k-1 +G e 11 k - ( Q k-1 +G e 10 k )) = G e 11 k -G e 10 k ). In Fig. 9, these 4 possible values are calculated with the help of 4 summators, and concurrently the comparators S⁰ and S¹ produce the logical signals d 0 / kand d 1 / k of eq. (21). Upon completion, the actual value of Q k is merely selected from the 4 possible values in a selection circuit S q under control of d 0 / k and d 1 / k. From eqs. (18), (20) and (22) it may be seen that these two bits provide exactly enough information for this selection. As selection is usually a faster process than addition, the configuration for calculating Q k in Fig. 9 may be applicable at higher data rates than the one of Fig. 7, though at the cost of additional hardware, notably 3 additional adders. The selector signals d 1ij / k for the adaptive precomputation units APU ij in Fig. 9 are generated by a decoder DEC that operates on the digits â 1 / k-D, â 1 / k-D+1 and â 1 / k-D+2 according to the decoding rule of eq. (25) for n=1 and the largest possible delay P=D-2. Further details of the receiver of Fig. 9 are not discussed here as they are entirely similar to those of the receiver of Fig. 7. For illustrative purposes it is mentioned here that the complete receiver according to Figs. 8 and 9 may be implemented with approximately 80 digital integrated circuits from the standard ECL 100K series as described, for example, in the "F100K ECL data book", Fairchild Camera and Instrument Corporation, Mountain View, California, 1982. In this implementation, internal signals of the receiver are represented with a word-length of at most 6 bits. The attainable data rate amounts to approximately 50 Mbit/s. Even for digital video storage applications this may be an appropriate value. To illustrate the merits of receivers according to the invention, Fig. 10 depicts bit error characteristics that were obtained by simulation for a receiver of prior art conforming to Fig. 7 (curve a.) and one according to the invention that conforms to Figs. 8 and 9 (curve b.). Both receivers operate on the output of an equalized optical recording channel with nonlinear ISI and a memory length M=2, to which random (NRZ) data is applied. The nonlinearities arise from a systematic difference in the length of the pits and lands that represent runs of zeros and ones. The curves of Fig. 10 pertain to a situation with severe nonlinear ISI, in which systematic errors in the writing process cause runs of zeros and ones to be T/2 seconds shorter and longer than their nominal value, respectively. This situation is illustrated in Fig. 11. In this Figure, the upper trace A depicts the NRZ waveform that is applied to the channel, while the lower trace B depicts the corresponding pattern of pits and lands that is assumed to be recorded on the optical medium. The systematic difference in the lengths of pits and lands manifests itself in the replayed signal as severe nonlinear ISI. In addition to this nonlinear ISI, the replayed signal is taken to contain linear ISI as a result of the channel bandwidth limitations that are reflected in Fig. 12. The curve that is labeled C in Fig. 12 depicts the transfer characteristic of the linear part of the channel. The loss of around 20 dB at the Nyquist frequency 1/(2T) is characteristic for recording at high information densities, and results in severe linear ISI in the replay signal. As both simulated receivers are only able to handle the comparatively small memory length M=2, an equalizer operating on the replayed signal is used to shorten the memory length of the channel into a memory length M of approximately 2 symbol intervals. Techniques for designing this equalizer are not discussed here as they are well described in the literature, see e.g. an article by D.D. Falconer and F.R. Magee, Jr Entitled "Adaptive Channel Memory Truncation for Maximum Likelihood Sequence Estimation", Bell Syst. Tech. J., Vol. 52, pp. 1541-1562, Nov. 1973. The amplitude-frequency characteristics of the equalizer are depicted in the curve that is labeled D in Fig. 12. Both the equalizer and the linear part of the channel have linear phase characteristics. A third disturbance, white Gaussian noise that models possible noise sources in the system, is added to the output signal of the channel, i.e. just before the input of the equalizer. Fig. 10 confirms the superiority of the receiver according to the invention (curve b.) over its conventional counterpart (curve a.) in dealing with nonlinear ISI. While the former receiver is unable to achieve useful performance levels even at very high signal-to-noise ratio's, the latter one already achieves bit error rates of around 10⁻⁴ for signal-to-noise ratio's of about 16 dB. Additional simulations reveal that this represents a loss of only 3 to 4 dB with respect to a corresponding situation without nonlinearities. Thus a receiver according to the invention may provide an attractive degree of insensitivity to nonlinear ISI, unlike its predecessors of prior art. For the sake of completeness, it is mentioned that the performance of the conventional receiver can be improved with respect to curve a. in Fig. 10 by preceding it by a DC-blocking circuit. This possibility stems from the fact that a major effect of the nonlinearity mechanism of Fig. 11 is a shift of DC-level. This effect is easily dealt with by conventional means, such as a DC-blocking circuit. For nonlinearity mechanisms with a more complicated nature than the one of Fig. 11, however, such a circuit may prove to be largely ineffective, whereas the receiver according to the invention is effective irrespective of the precise nonlinearity mechanism.

Claims: System for transmitting a data signal at a symbol rate 1/T through a noisy dispersive channel (CHN) to a data receiver (REC); said channel (CHN) introducing intersymbol interference and noise into the transmitted data signal, said receiver (REC) comprising means for estimating the most likely sequence of transmitted data symbols by keeping track of candidate sequences ( S ) and associated likelihood measures (I), means for determining said likelihood measures by comparison of the actual channel output signal (r k ) with a hypothetical channel output signal that would have been received if noise was absent and said candidate data sequence would have been transmitted, and means for estimating said hypothetical channel output signal from said candidate sequence, characterized in that said means for estimating said hypothetical channel output signals comprise one or more look-up tables storing said possible hypothetical channel output signals corresponding to said candidate data sequences. System according to claim 1, characterized in that the means for estimating the most likely sequence of transmitted data symbols is arranged for keeping track of two candidate sequences (S¹, S°) and a difference likelihood measure (Q k ) that is representative for the difference of a function of the likelihoods of both candidate data sequences. System according to claim 2, further characterized in that said difference likelihood measure (Q k ) is determined by selection among precomputed candidate values of said likelihood measure. System according to claim 1, 2 or 3, further characterized in that said look-up tables (LUT) take the form of registers that store hypothesized channel output symbols in the absence of noise. System according to claim 1, 2 3, or 4 further characterized in that each look-up table (LUT) is adapted under the control of digits of said candidate data sequences, in response to an error signal (e k ) that is representative for the difference of the channel output signal and the output signal of said look-up table. System according to claim 1, 2 or 3, further characterized in that each look-up table (LUT) is adapted in response to an error signal (e k' ) that is representative for the difference of a delayed version of the channel output signal and the output signal of said look-up table (LUT) when addressed by one of more delayed digits of said candidate data sequences. System according to claim 4, further characterized in that each register takes the form of a digital counter (C) that is adapted under the control of one or more delayed digits of said candidate data sequences, in response to an error signal that is representative for the difference of a delayed version of the channel output signal and the contents of said counter (C). System according to claims 1, 2, 3, 4, 5, 6 or 7, further characterized in that each likelihood measure is representative for an accumulated version of a function that essentially equals the modulus of the difference of the actual channel output signal and a hypothesized channel output signal in the absence of noise. System according to claim 8, further characterized in that said function is determined by selection among precomputed candidate values G(l k ) of said function.

Inventors: BERGMANS JOHANNES WILHELMUS MA, IZUMITA MORISHI, MITA SEIICHI, BERGMANS, JOHANNES WILHELMUS MARIA, IZUMITA, MORISHI, MITA, SEIICHI
Applicants: HITACHI LTD, PHILIPS ELECTRONICS NV, HITACHI, LTD., PHILIPS ELECTRONICS N.V.
Application year: 1990
Main classification: H04L1
IPC codes: H04L25, H04L1, G11B20