IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS. VOL. 9. NO. 6. AUGUST 1991
`
`895
`
`A Discrete Multitone Transceiver System for HDSL
`Applications
`
`Jacky S. Chow, Student Member, IEEE, Jerry C Tu, Student Member, IEEE,
`and John M. Cioffi, Senior Member, IEEE
`
`Abstract-A discrete multitone (DMT) transceiver design for
`high-bit-rate digital subscriber line (HDSL) access is presented
`and analyzed. The DMT transmitter and receiver structure and
`algorithms are detailed, and the computational requirements
`of DMT for HDSL are estimated. We find that at a sampling
`rate of 640 kHz, using an appropriate combination of a short
`FIR equalizer and a length-512 DMT system, 1.6 Mb/ s data
`transmission is possible within the CSA at an error rate of 10- 1
`on a single twisted pair. We also show significant performance
`margin can be achieved when two coordinated twisted pairs are
`?sed to deliver a total data rate of 1.6 Mb/ s. We conclude that,
`m terms of a performance-per-computation figure of merit, the
`DMT system is an excellent candidate for HDSL implementa(cid:173)
`tion.
`
`I. INTRODUCTION
`
`MUL TITONE modulation methods use an optimized
`
`frequency division allocation of energy and bits to
`maximize the reliably achievable data rates that can be
`transmitted over bandlimited communication channels.
`These systems are easily described because they use fre(cid:173)
`quency division multiplexing to transport a single-input
`data stream on several carriers within the usable fre(cid:173)
`quency band of the channel. The simplicity of multitone
`methods accounts for their use in some of the earliest data
`transmission modems, such as the Collin's Kineplex mo(cid:173)
`dems of 1957 [l). While the simplicity and appeal ofmul(cid:173)
`titone modulation has existed for decades, it is only very
`recently that their optimal performance has been proven
`theoretically (see [2], [3]) and demonstrated in practice
`with recent products, such as Telebit's Trailblazer mo(cid:173)
`dems [4) or NEC's multitone voiceband and groupband
`modems [5]. These highest-performance products are sold
`at only a fraction of the cost of comparable performance
`single-channel or QAM modems.
`The excellent high-performance/cost tradeoff of multi(cid:173)
`tone modulation also makes it a strong candidate for the
`transceiver implementation for high-bit-rate digital sub(cid:173)
`scriber lines (HDSL), a 1.6 Mb/ s digital subscriber ser(cid:173)
`vice on twisted-pair channels of up to two miles. This
`paper presents a multitone modulation structure that is of
`
`Manuscript received November 14, 1990; revised May 17, 1991. This
`work was supported in part by a gift from Bell Communications and by a
`contract with the University Technology Transfer Institute.
`The authors are with the Information Systems Laboratory. Department
`of Electrical Engineering, Stanford University, Stanford, CA 94305.
`IEEE Log Number 9102134.
`
`particularly high performance and low cost for the HDSL
`application. We generally find that, even on worst-case
`channels, a 1.6 Mb/ s (1.536 Mb/ s user data plus 64 kb/ s
`control channel) data rate is achievable at an error rate of
`10- 7 on a single twisted pair (with 6 dB margin). Fur(cid:173)
`thermore, we find that the described multitone system ex(cid:173)
`hibits excellent margins (nearly 18 dB) for the so-called
`dual-duplex option of providing 1.6 Mb/ s on two coor(cid:173)
`dinated twisted pairs. 1
`The intersymbol interference, crosstalk, and attenua(cid:173)
`tion of 24-gauge and/or 26-gauge twisted pair is particu(cid:173)
`larly severe within the two-mile carrier serving area (CSA)
`requirements, as we discuss in Section II. Section II also
`briefly examines and reviews theoretically and practically
`achievable maximum data rates for HDSL. Equalization
`and coding methods have to be carefully applied in the
`HDSL application to render the high transported data rates
`reliable. In the approach presented, a particularly effec(cid:173)
`tive combination of a small amount of time-domain equal(cid:173)
`ization and of a larger amount of multitone modulation is
`used to transmit the data reliably. The equalizer is used
`to shorten the effective response length of the data chan(cid:173)
`nel, while the multitone modulation is used to mitigate
`the remaining intersymbol interference and the effects of
`crosstalk, as is described in Section III. The multicarrier
`system is fully adaptive, employing both an adaptive al(cid:173)
`gorithm for update of the adaptive filters and also a spe(cid:173)
`cial channel identification procedure that is used to ini(cid:173)
`tialize and update
`the
`transmitter bit allocation
`to
`frequency-indexed subchannels. This relatively simple
`procedure is essential for maximizing performance. We
`also estimate the complexity in terms of the number of
`arithmetic operations in Section III and illustrate that the
`multitone transceiver can be implemented with program(cid:173)
`mable signal processors, although those processors may
`need to be custom designed for lowest cost.
`The performance in achievable data rate of the multi(cid:173)
`carrier system is analyzed and characterized in Section IV
`as a function of input power, channel characteristic, num(cid:173)
`ber of tones or carriers, sampling rate, noise margins, and
`crosstalk.
`
`1We note that results in this paper use channel models with less atten(cid:173)
`uation than those of previous studies in [6)-[9), based on measured line
`results of several industrial sources.
`
`0733-8716/91/0800-0895$01.00 © 1991 IEEE
`
`... IU I lijJ L lijJ~
`
`llllll 1mt
`
`!lilt
`
`I l 11
`
`HTC Corp., HTC America, Inc. - Ex. 1018, Page 1
`IPR2018-01555 and IPR2018-01581 (HTC and Apple v. INVT SPE)
`
`

`

`896
`
`IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS, VOL. 9. NO. 6, AUGUST 1991
`
`II. SUBSCRIBER LOOP CHARACTERIZATION AND
`MODELING
`In any study of HDSL transceivers, the channel models
`are critical to specification of performance. We outline in
`detail the HDSL channel models used throughout the study
`of DMT HDSL transceivers. We note that performance
`projections based on different channel models can lead to
`erroneous claims and comparisons. A detailed tutorial on
`HDSL technologies is presented by Lechleider in this is(cid:173)
`sue [10].
`
`A. CSA Loop Characteristics
`The carrier serving area (CSA) is an identifiable subset
`of the current subscriber loop population in the U.S., and
`our HDSL study is restricted to loops that satisfy minimal
`CSA requirements. The CSA consists of mainly 24- and
`26-gauge twisted pair, with lengths on 24-gauge wire to
`12 kft and on 26-gauge wire to 9 kft. CSA loops can also
`contain bridge taps of limited length; a complete specifi(cid:173)
`cation of the CSA appears in [11). Impairments in the
`CSA data transmission channel include:
`• intersymbol interference (ISi),
`• crosstalk noise coupled from adjacent loops within
`the same cable bundle,
`• impulsive noise (caused by switching transients,
`lightning, and other electrical machinery),
`• echo noise resulting from the combination of imper(cid:173)
`fect hybrids and gauge changes on the line as well as sig(cid:173)
`nals reflected from bridged taps,
`• electronics noise (such as quantization noise in ana(cid:173)
`log-to-digital converters (ADC) and thermal noise in the
`analog portion of the transmitter and receiver), and
`• inductive noise at 60 Hz and its harmonics caused by
`nearby power lines.
`We consider a set of representative CSA channels, pro(cid:173)
`vided by Bell Communications Research [12), that have
`the respective loop cable arrangements illustrated in Fig.
`1. We use a letter designation (A-H) on these channels
`for reference. These channels are generally representative
`of some of the worst, in terms of channel impulse re(cid:173)
`sponse, expected to be found within the CSA limits. In
`general, longer wire length, smaller wire size, and pres(cid:173)
`ence of more and longer bridge taps results in poorer
`achievable data rates. We also note that on-premises wir(cid:173)
`ing is often impossible to characterize statistically, so a 6
`dB margin is imposed upon the performance of any HDSL
`transceiver to ensure that 10- 7 error rate is maintained
`[ 13). In all loop models in this paper, we include this 6
`dB margin simply by increasing loop attenuation for the
`loops in Fig. 1 by 6 dB (see also Footnote 1).
`
`B. Channel Modeling
`A line simulation program, based on standard trans(cid:173)
`mission line modeling and a standard set of twisted-pair
`characteristics, was used to generate an impulse response
`for each of the channels in Fig. 1. The effect of the cou(cid:173)
`pling transformer at each end of the loop termination was
`
`Channel 25
`
`1778'/24 190'/24
`
`A
`
`6451 '/26
`
`305'/26
`
`Channel 29
`B
`
`llll2'/2f
`
`180'/22
`
`Channel 23 196'/24
`C
`1786'/26
`6514'/26
`
`I 80'/26
`
`Channel 211
`140'/21
`E o--11-64-3-'/24--
`
`Channel 2 IO
`8100'/26
`F
`
`800'/22
`
`Channel 27
`G
`
`12000'/24
`
`Channelo,.26 _ _ _ _ -o
`11522'/24
`D
`Fig. l. CSA channels: length (in ft.)/gauge.
`
`Channel 16
`H
`
`9000'/26
`
`included by adding a double-pole double-zero highpass
`digital filter with a cutoff frequency of 300 Hz. These
`characteristics were found to be in near exact agreement
`with those provided by separate programs at various in(cid:173)
`dustrial groups working in the HDSL field [14), [12). One
`such time-domain impulse response is shown in Fig. 2 for
`Channel H.
`Crosstalk noise coupled from adjacent wire pairs in the
`cable bundle is the dominant noise impairment in CSA
`loops. Various characterizations have been proposed [15),
`[16), and we will use [15), which models crosstalk noise
`as near end crosstalk (NEXT) with the following transfer
`function:
`
`(1)
`
`where KNEXT is the coupling coefficient. Since far end
`crosstalk noise is attenuated by the channel, the resulting
`noise power is negligible compared to near end crosstalk.
`Hence, we will only consider NEXT. The coupling coef(cid:173)
`ficient KNEXT for a 50-pair cable was empirically deter(cid:173)
`mined to be approximately 10- 13 [17).
`To model the effects of residual echo after cancellation,
`we note that a typical CSA loop attenuation is on the order
`of 24 dB, and a good echo canceller can presumably re(cid:173)
`duce the echo to a level of 30-40 dB below the received
`signal. Assuming an average input power of 10 mW (10
`dBm), this translates into -44 to -54 dBm of total echo
`noise power across the two-sided bandwidth of 640 kHz,
`which corresponds to a noise power spectral density of
`6.22 x 10- 11 to 6.22 x 10- 12 mW /Hz. This residual
`echo noise is assumed to be additive white Gaussian noise
`(AWGN), and we fix this A WGN to have a power spectral
`density of -110 dBm/Hz for subsequent analysis and
`simulation. The electronics noise, as mentioned earlier,
`is typically on the order of -140 to -170 dBm/Hz [6)
`and is thus, by comparison to echo noise, negligible. Of
`course, crosstalk noise is larger than either echo or elec(cid:173)
`tronics noise.
`While impulsive noise is more difficult to characterize,
`some [ 18) have attempted to model its effect. Impulsive
`noise usually consists of millivolt pulses (say, up to 40
`m V) that last a few µs to I 000 µs.
`Inductive noise can be modeled as a fixed noise level
`at 60 Hz, but the highpass filter simulating the trans(cid:173)
`formers eliminates the inductive noise.
`1) Pole-Zero Modeling: To facilitate simulation of
`CSA loops for performance analysis, we constructed a
`pole-zero (IIR) model for each channel impulse response
`
`- -~ -
`
`HTC Corp., HTC America, Inc. - Ex. 1018, Page 2
`IPR2018-01555 and IPR2018-01581 (HTC and Apple v. INVT SPE)
`
`

`

`CHOW et al.: A DISCRETE MULTITONE TRANSCEIVER SYSTEM
`
`897
`
`0 . 0 3 5 , - - -~~ - -~ - - -~ - - -~ - - -~
`
`TABLE I
`CAPACITY AND FLAT INPUT CAPACITY FOR CSA LOOPS, 800 KHZ RATE AND
`WITH 6 DB MARGIN (IN MB/S)
`
`0.03
`
`0.025
`
`0.02
`
`0.015
`
`O.Ql
`
`0.005
`
`o-.
`
`Bellcore line simulation
`-
`---- Pole-Zero Model
`
`Capacity with 6 dB
`Margin
`
`Flat Input Capacity with 6 dB
`Margin
`
`Channel
`
`Crosstalk
`Bound
`
`II dBm
`Power
`
`OdBm
`Power
`
`10 dBm
`Power
`
`17 dBm
`Power
`
`A
`B
`C
`D
`E
`F
`G
`H
`
`2.42
`2.39
`2.32
`2.31
`2.26
`2.29
`2.19
`2.14
`
`2.30
`2.27
`2.20
`2.19
`2.14
`2.17
`2.07
`2.02
`
`1.87
`1.85
`1.78
`1.77
`1.72
`1.75
`1.65
`1.60
`
`2.29
`2.27
`2.19
`2.19
`2.13
`2.16
`2.06
`2.01
`
`2.40
`2.37
`2.30
`2.29
`2.24
`2.27
`2.16
`2.11
`
`(10 dBm or above), the penalty in achievable data rate
`without optimizing the transmit power distribution be(cid:173)
`comes negligible. Thus, we will assume a flat input power
`spectrum is used, and this will simplify analysis and sim(cid:173)
`ulation. Also, this allays any practical concern about loop(cid:173)
`to-loop energy spectrum variation.
`We can infer from this table that 1. 6 Mb/ s can be re(cid:173)
`liably transmitted over most CSA channels. However, an
`excellent transceiver design is required to accomplish this
`goal on either a single pair or on two coordinated pairs
`with a large performance margin such as 12 dB. Even with
`only a 6 dB performance margin on two pairs, 800 kb/ s
`is 35-50% of capacity and more than a trivial receiver is
`required. On the other hand, a data rate of 160 kb/ s on
`CSA channels (as in the case of basic rate ISDN) is well
`below the fundamental limits (although it may be a little
`more difficult on longer (18 kft) channels). Thus, one
`would not expect that a simple receiver used for basic rate
`ISDN would necessarily be good for use in HDSL.
`
`Ill. THE DISCRETE MULTITONE TRANSCEIVER
`This investigation focuses on a discrete multitone
`(DMT) transceiver for HDSL. We briefly describe the
`modulation and demodulation processes and the major
`system components. We also accurately compute com(cid:173)
`putational requirements for the DMT HDSL transceiver
`using Motorola's 56000 DSP as a representative host for
`implementation. A system using multiple Motorola 56000
`DSP's has been used to implement the described DMT
`transceiver in real-time at Stanford University.
`
`A. System Description
`A simplified illustration of the DMT transmitter, re(cid:173)
`ceiver, and channel appears in Fig. 3. The basic concept
`is to use frequency division modulation to divide the
`transmission system into a set of frequency-indexed sub(cid:173)
`channels that appear to be modulated and demodulated
`independently. With a careful allocation of bits and trans(cid:173)
`mit power to the subchannels, it is well known [2], [3]
`that such a system is capable of performing at the highest
`theoretical limits and that no other system can exceed its
`
`-o.oo5a:--------::50=-----1~00----1~50----2~00 _ _ _ _J
`250
`
`Fig. 2. Impulse response of a CSA loop.
`
`at various sampling rates, using standard system identifi(cid:173)
`cation tools. An impulse response obtained for Channel
`H is compared to the impulse response generated by a
`pole-zero model in Fig. 2. Since the pole-zero model
`closely approximates the true response, the use of pole(cid:173)
`zero models for subsequent studies is justified.
`The number of poles and zeroes varies with different
`CSA channels, but they all fall into the same range: 5-8
`zeroes, and 6-12 poles. For future reference, assume the
`number of zeroes for a particular channel is equal to v,
`and the number of poles is equal to L - 1. We also found
`that most of the CSA channels are not minimum phase,
`i.e., some of the zeroes from the pole-zero models lie out(cid:173)
`side the unit circle in the z-plane.
`
`C. Capacity
`The capacity of these loops with the aforementioned
`models is computed in [6], at the sampling rates of 400
`and 800 kHz. As explained later in Section IV-A and [6],
`a sampling rate of 400 kHz is not sufficient to maximize
`the HDSL transceiver performance on most CSA chan(cid:173)
`nels shown in Fig. 1 as well as the capacity with a flat
`input spectrum (called "flat input capacity" in Table I)
`at 800 kHz sampling rate. The two columns under "Ca(cid:173)
`pacity" list the capacity using the optimal transmit power
`distribution [6], and represent the best data rates that can
`be achieved with an 800 kHz sampling rate. "Crosstalk
`Bound'' assumes crosstalk noise is the only noise impair(cid:173)
`ment (and, thus, ignores the contribution of A WGN),
`while "11 dBm Power" includes the effect of A WGN
`with 11 dBm transmit input power. The three columns
`under "Flat Input Capacity" do not optimize transmit
`power distribution; instead, a flat input power spectral
`density is assumed. The effects of near end crosstalk, 6
`dB margin, and A WGN are included for the calculation,
`and the maximum achievable data rate based on the flat
`input power spectrum are shown for three different input
`powers: 1, 10, and 50 mW (0, 10, and 17 dBm, respec(cid:173)
`tively). Notice that, with a sufficiently high input power
`(10 dBm or above), the penalty in achievable data rate
`
`... IU l l~J L l~J~
`
`IIJlll
`
`llllt
`
`1m1.
`
`11
`
`111
`
`HTC Corp., HTC America, Inc. - Ex. 1018, Page 3
`IPR2018-01555 and IPR2018-01581 (HTC and Apple v. INVT SPE)
`
`

`

`898
`
`IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS. VOL. 9. NO. 6. AUGUST 1991
`
`Transmitter
`
`Receiver
`
`Fig. 3. Simplified illustration of DMT transceiver and channel.
`
`performance. There are now at least three commercial
`19 .2 kb/ s voiceband modems, marketed by three differ(cid:173)
`ent companies, that use DMT to attain excellent perfor(cid:173)
`mance. This section describes a DMT system for HDSL.
`The system in Fig. 3 is presumed to be sampled (digi(cid:173)
`tized) at a rate that is at least twice as high as any signal
`frequency, so that sampling theory arguments apply and
`the discrete-time system is equivalent to the actual con(cid:173)
`tinuous-time system. A block of input bits is divided into
`N subsymbols, xn, n = 0, · · · , N -
`I, that are indepen(cid:173)
`dently modulated by N-dimensional sampled-sinusoid
`modulating vectors Pm and then summed to form an
`N-dimensional block of modulated transmit signal sam(cid:173)
`ples ib k = 0, · · · , N -
`I. The nth modulating vector
`is given by:
`
`Pn
`
`[Pn.O P11.I ·'' Pn.N-il
`
`(2)
`
`where
`
`~ _I_ e+j(h/N)kn
`n.k - JN
`P
`
`k, n = 0, · · · , N -
`
`I
`
`(3)
`
`the inverse discrete Fourier transform (IDFT) vector. To
`ensure a real-valued transmit signal, ik = lR(ik), it is nec(cid:173)
`essary that we impose the constraint on the subsymbols
`that (* denotes complex conjugate, lR denotes real part)
`
`n=l,···,N-1.
`(4)
`It is convenient to think of DMT as consisting of N /2
`quadrature amplitude modulation (QAM) channels. In
`practice, of course, the modulation is done digitally and
`no explicit use of any carrier modulation frequencies or
`modulators is required. The channel output samples Yk are
`found by the sum of the convolution of the modulated
`transmit samples, ib with the sampled channel impulse
`response hk and the additive noise nb
`
`(5)
`
`The impulse response hk is assumed to be of finite length,
`and can be described by the v + I samples h0 ,
`, hv.
`• • •
`Its D-transform is h(D) = Ei'.=o hkDk. In the receiver, the
`corresponding demodulation is accomplished by taking the
`inner product of the corresponding N channel outputs, y0 ,
`
`· · · , YN _ 1, with p:. The noise at the channel output pres(cid:173)
`
`ently includes only additive white Gaussian noise; we will
`generalize to the crosstalking case as this development
`proceeds.
`Because of the finite number of subchannels, a cyclic
`prefix (CP) of the transmitted signal samples ik is needed
`once for every transmitted block of bits. The CP prefixes
`the block of transmitted samples i 0 ,
`, iN _ 1 by the
`• • •
`(repeated) samples, iN- "' · · · , iN- 1• We denote the cor(cid:173)
`responding time indexes at the beginning of the block by
`k = - v, - v + I , · · · , - 1 . When v + I is the length
`of the channel impulse response in sample periods, then
`the successive blocks of transmitted data do not interfere
`with each other, and one can write (see [3]) that
`
`(6)
`
`where
`Hn =ho+ h1e-j(21r/N)n + h2e-j(21r/N)n + ...
`+ hve-j(h/N)vn.
`
`(7)
`In practice, the modulation and demodulation process(cid:173)
`ing is implemented with fast Fourier transform (FFT) al(cid:173)
`gorithms that require significantly less computations than
`would occur in the direct implementation shown in Fig.
`3. We discuss specific computational requirements in Sec(cid:173)
`tion III-C.
`A trellis code [ 19] of coding gain 'Y can be applied to
`each subchannel so that the energy allocated to this chan(cid:173)
`nel will sustain a greater number of bits than in the ab(cid:173)
`sence of the applied trellis code. The number of bits can
`be approximated by
`
`SNRn)
`(
`bn = log2 1 + -r-
`
`(8)
`
`where r = 10 dB -
`'Y for a bit error rate of 10- 7 [20]. If
`no code is used, 'Y = 0 dB and r = 10 dB. The signal(cid:173)
`to-noise ratio on the nth subchannel is denoted by SNRn
`and is computed by
`
`n
`
`(9)
`
`SNR = Ex 1Hnl2
`a2
`where Ex = E lxnl2 is the signal energy of the nth subsym(cid:173)
`bol and a 2 is the corresponding white noise energy. In
`practice with HDSL, the dominant noise is crosstalk
`noise, which is not white, so that a 2 -+ a~, a subchannel(cid:173)
`dependent noise energy. When N is chosen sufficiently
`large, the noise on each of the subchannels will be ap(cid:173)
`proximately white within the occupied frequency band and
`will be uncorrelated with the noise on any of the other
`subchannels. In order to compute the bit allocation out(cid:173)
`lined by (8), the channel characteristic must be known at
`the transmitter, which is achieved as described in Section
`III-A-3. The energy allocation is approximated by a flat
`
`HTC Corp., HTC America, Inc. - Ex. 1018, Page 4
`IPR2018-01555 and IPR2018-01581 (HTC and Apple v. INVT SPE)
`
`

`

`CHOW el al.: A DISCRETE MULTITONE TRANSCEIVER SYSTEM
`
`899
`
`N/2
`subchannels
`
`N/2
`QJ\M signals
`
`---- .-----~ ~---~
`
`R=l.6 Mbps
`
`Bit Parsing
`(M-bit frames)
`and Buffering
`
`Coset-Encoder
`and
`Modulation
`plus Tone Shuffle
`
`Inverse
`Length-N
`FFr
`
`u-Tap
`Cyclic Prefix
`(parallel/serial)
`
`DAC
`and
`Transmit Filter
`(S = 640 kHz)
`
`Unit
`
`Included for
`Trellis Coding
`
`N=512
`
`Fig. 4. Block diagram of a specific DMT HDSL transmitter.
`
`transmit energy spectrum, so that Ex is not a function of
`the frequency index n. Our studies have found that very
`little is lost in using this flat energy distribution in the
`HDSL application. The poorest subchannels carry only
`dummy data and are simply ignored by the receiver. This
`flat transmit spectrum has the effect of rendering the
`crosstalk power spectrum from any DMT HDSL trans(cid:173)
`mitter constant in practice.
`1) The DMT HDSL Transmitter: A specific embodi(cid:173)
`ment of the HDSL transmitter is illustrated in block dia(cid:173)
`gram form in Fig. 4. In the specific implementation of the
`system that we study, " is always limited (through equal(cid:173)
`ization, as discussed later) to a maximum of 8, and later
`in Section IV-A-3) and 4), we fix the sampling rate S to
`be equal to 640 kHz and N = 512.
`For a single twisted pair to carry R = 1.6 Mb/ s (the
`single-duplex option), the input data bit stream is parsed
`into M-bit blocks, where M = (N + P)R/S = 1300. These
`bits are then transformed into (a maximum of) N /2 =
`256 QAM subsymbols that are then applied to a trellis
`encoder. The trellis coder sequentially processes the fre(cid:173)
`quency-indexed subsymbols to avoid the large latency and
`memory requirement that would occur with multiple trel(cid:173)
`lis encoders. This process requires that the subchannel in(cid:173)
`dexes be "shuffled" to avoid any minor correlation be(cid:173)
`tween the noise on adjacent frequency subchannels. The
`output of the trellis encoder is then modulated, as de(cid:173)
`scribed earlier, through the use of an N = 512 inverse
`FFT. Finally, an eight-sample cyclic prefix is placed in
`the beginning of the corresponding block of modulated
`transmit samples, and the extended block of N + "sam(cid:173)
`ples is then applied to the channel through the digital-to(cid:173)
`analog converter (DAC) and line-interface unit. The sam(cid:173)
`pling rate is 640 kHz, which leaves the 256 subchannels
`effectively separated by 1.25 kHz. This corresponds to a
`block symbol period of 512 + 8 = 520 samples or 812. 5
`µs. Eight samples is considerably shorter than most HDSL
`channel impulse responses at 640 kHz. An unusual and
`simple equalizer is used in the receiver to contain inter(cid:173)
`symbol interference to eight samples, or impulse-re(cid:173)
`sponse length to nine samples.
`For the so-called dual-duplex option with two twisted
`pairs, the 1300-bit blocks are parsed into the 512 QAM
`subchannels that exist on both of the channels. There are
`then two IFFT's, two DAC's, and two line-interface units.
`
`However, only one trellis code is necessary with this op(cid:173)
`tion with respect to the single-duplex option. This im(cid:173)
`proves the performance margin but increases the cost of
`the dual-duplex option. If a sufficient margin can be ob(cid:173)
`tained using the single-pair solution (called "single-du(cid:173)
`plex"), then it is preferable from a cost standpoint.
`2) The DMT HDSL Receiver: A block diagram of the
`general DMT HDSL receiver appears in Fig. 5, and a
`specific receiver that corresponds to the specific transmit(cid:173)
`ter of Fig. 4 is shown in Fig. 6. This receiver structure is
`a slight variation of the original receiver proposed in [7].
`Both receiver structures achieve the same performance but
`the one in Fig. 6 results in lower complexity.
`An actual HDSL channel response is of significant du(cid:173)
`ration in practice, and we can model its polynomial-frac(cid:173)
`tion D-transform
`
`h(D)
`
`a(D)
`p(D)
`
`(10)
`
`In practice, we identify this infinite-impulse-response
`model, h(D), by choosing a(D) and p(D) to minimize the
`mean square error between this model and the channel
`output. Using this method, the crosstalk noise also affects
`the settings for a(D) and p(D). Due to the nonminimum
`phase nature of a(D), the use of an IIR equalizer 1 /h(D)
`to eliminate ISi results in unstable filter. A short, fixed
`L-tap feedforward equalizer is thus used to process the
`sampled channel output Yk sequentially. This equalizer
`does not remove intersymbol interference completely, but
`rather confines the impulse response so that its length is
`approximately " + 1 sample periods or less. By setting
`the equalizer to p(D), we will have a minimum mean
`square error approximation of an additive white noise
`channel with impulse response characterized by a(D),
`which has a response length of " + 1 sample periods or
`less. The last N samples of the N + " samples that cor(cid:173)
`respond to the transmit block are extracted from the
`equalizer output. An N-point FFT is performed on the
`equalizer output. The FFT outputs, Zn, (n = 0, 1, · · · ,
`N - 1), are multiplied by N complex, 1-tap adaptive fil(cid:173)
`ters, wm (n = 0, 1, · · · , N -
`1) so that a common
`decision device can be used to estimate the subsymbols
`on each of the subchannels. The initial tap setting for Wn
`
`HTC Corp., HTC America, Inc. - Ex. 1018, Page 5
`IPR2018-01555 and IPR2018-01581 (HTC and Apple v. INVT SPE)
`
`

`

`900
`
`IEEE JOURNAL ON SELECTED AREAS IN COMMUNICATIONS. VOL. 9. NO. 6. AUGUST 1991
`
`Fig. 5. Block diagram of the receiver.
`
`Fig. 6. Specific block diagram of a DMT HDSL receiver.
`
`is:
`
`where
`
`n = 0, 1, · · · , N - 1
`
`(11)
`
`N-1
`An= I; ake-j(27r/N)kn,
`k=O
`
`n = 0, 1, · · · , N - 1
`
`the FFT of a(D). The resulting output data
`
`n = 0, 1, · · · , N - 1
`
`(13)
`
`(12)
`
`ceiver in practice, but are still important. (We note that
`decision feedback approaches that employ Tomlinson
`precoding [21], [22] also need to know the channel at the
`transmit side and require an even greater accuracy than
`DMT [8], even when an adaptive feedforward section is
`used in the DFE receiver.) Rather than pursue channel
`identification in detail (see, for example, [23]-[25] for
`details), we will briefly summarize the basic steps in the
`initialization protocol. The complete startup time (includ(cid:173)
`ing channel identification, SNR calculation, and bit and
`energy allocation) is less than one second.
`1) After receiving a request for transmission signal,
`some loopback and test patterns are exchanged, followed
`by timing and synchronization patterns.
`2) A predefined pseudorandom sequence is sent from
`the transmitter to the receiver for spectral estimation. A
`pole-zero model for the channel is derived from the re(cid:173)
`ceived training pattern using a least-squares fit.
`3) Based on the result of the above channel spectral
`estimation, the receiver sets the tap values for the equal(cid:173)
`izer at the pole settings of the pole-zero model, then com(cid:173)
`putes the SNR for each frequency-indexed subchannel,
`and finally computes the transmit bit and energy alloca(cid:173)
`tion across the frequency band.
`4) The bit assignments are sent to the transmitter in a
`reliable fashion.
`5) For bidirectional transmission, the above steps are
`repeated for the other direction.
`6) After frame synchronization and start commands are
`exchanged, the DMT transceiver system is ready for ac(cid:173)
`tual data transmission.
`We have not considered the implementation of echo
`cancellation in this paper. We note that low-computation/
`high-performance echo cancellation is possible with DMT
`(see [26]).
`
`B. Transmit Power Level
`
`As shown in Section IV-A, we can improve perfor(cid:173)
`mance by increasing the average input power. However,
`due to physical constraints, only a certain range of aver(cid:173)
`age input power level is feasible. As a result, maximum
`average input power or energy is a common constraint
`that is usually imposed at the input of the channel ( or
`equivalently, at the output of the transmitter).
`I) Transmit Power Adjustment: We have assumed an
`average energy constraint of Ex per transmitted subsym(cid:173)
`bol, before adding the cyclic prefix. The IFFT and FFT
`do not change the energy of a block of data, so the energy
`constraint specified at the output of the IFFT in the trans(cid:173)
`mitter is the same as at the input of the IFFT, i.e., the
`total average energy per block at the input of the IFFT is
`limited to NEx. However, the total average energy of a
`v-sample cyclic prefix is not necessarily equal to vEx. We
`will thus briefly investigate the impact of the cyclic ex(cid:173)
`tension samples to the average energy constraint.
`The real and imaginary parts of the noise variance after
`the FFT in the receiver are, in general, unequal for any
`
`can then be decoded, which in the case of an applied trel(cid:173)
`lis code requires a Viterbi decoder as shown. The symbol
`decisions are used (in the coded case) only to derive an
`error signal, so that the adaptive updating mechanism can
`be used to allow for slight channel variation. The corre(cid:173)
`sponding error is
`e,, = d11
`
`n = 0, 1, .. · , N - 1.
`
`(14)
`
`v,,,
`
`-
`
`(15)
`
`Standard LMS algorithm can be used for updating,
`n = 0, 1, · · · , N - 1
`where /1n is the stepsize for each subchannel, and can be
`adjusted independently to optimize the rate of conver(cid:173)
`gence.
`Since ak and Uk are real, An = At-n and Zn = zt-n, so
`that (13)-(15) are computed only through n = N /2.
`3) Channel Identification and Initialization Protocol:
`As the DMT transmitter requires knowledge (either di(cid:173)
`rectly, or indirectly through knowing the bit allocation b11 )
`of the channel, the channel impulse response and crosstalk
`noise spectrum must be estimated during initialization.
`The more accurate the initial estimate, the more reliable
`the transmission with DMT. Accuracy requirements are
`significantly mitigated by the adaptive filters in the re-
`
`HTC Corp., HTC America, Inc. - Ex. 1018, Page 6
`IPR2018-01555 and IPR2018-01581 (HTC and Apple v. INVT SPE)
`
`

`

`CHOW er al.: A DISCRETE MULTITONE TRANSCEIVER SYSTEM
`
`901
`
`specific frequency band due to the correlated, non-i.i.d.
`crosstalk noise at the input to the FFT. This inequality of
`variance is more profound at smaller blocklength (larger
`than 4), and the difference between the real and imaginary
`noise variance gradually diminishes as longer blocklength
`is used. As blocklength approaches infinity, the difference
`becomes zero. To optimize bit error rate throughout the
`subchannels, we would then need to assign unequal signal
`energy on the real and imaginary part of each frequency
`band. A fairly tight bound for the average energy of a
`block of N + v samples (after cyclic extension) is derived
`in the Appendix.
`However, as the blocklength increases, the difference
`between