M A IEEE CircuitsCircuitsCircuits G A andand Z and I N SystemsSystemsSystems E Volume 1, Number 2, Second Quarter 2001 ISSN 1531-636X

Filter Banks in Digital Communications 4

A Microphone Array for Hearing Aids 26 CallsCalls forfor PapersPapers andand ParticipationParticipation

8th IEEE International Conference on CALL FOR PAPERS st Electronics, Circuits and Systems 1 IEEE International Conference on Circuits and Systems for Communications ICECS’01 St.Petersburg, Russia June 26–28, 2002 September 2–5, 2001 “Circuits and Systems in Broadband Communication Technologies” Author’s Schedule: The Westin Dragonara Resort, Malta Deadline submission of extended abstract or full paper: December 17, 2001 Deadline for notification of acceptance: February 25, 2002 http://www.eng.um.edu.mt/microelectronics/icecs2001 Deadline for final version: March 29, 2002 ICECS is a major international conference which includes regular, special and poster The conference is sponsored by the IEEE Circuits and Systems Society. The conference sessions on topics covering analogue circuits and signal processing, general circuits and topics include questions and problems that are around the theory and design of circuits systems, digital signal processing, VLSI, multimedia and communication, computational and systems for communications applications. Signal processing, RF design and micro- methods and optimization, neural systems, control systems, industrial and biomedical electronic implementations of such types of circuits and systems are of interest. A cul- applications, and electronic education. tural program including the Hermitage, museums, and beautiful sceneries around St. General Chair Technical Program Chair Petersburg will be available as well. Last but not least, visitors are expected to capture a Dr. Joseph Micallef Prof. Franco Maloberti lifetime impression of the famous season of June in St. Petersburg called “White Nights”. Department of Microelectronics Department of Further Information: Secretariat: [email protected] University of Malta Texas A&M University Technical Program Committee (Prof. A.S.Korotkov): [email protected] Msida MSD06, Malta College Station, Texas 77843 Conference Web site: http://iccsc.spbstu.ru/ Email: [email protected] Email: [email protected] ICECS 2001 Secretariat The 2001 IEEE Workshop on Department of Microelectronics, Faculty of Engineering SiGNAL PROCESSING SYSTEMS (SiPS) University of Malta, Msida MSD 06, Malta, Europe Design and Implementation E-mail: [email protected] Tel: (+356) 346760; Fax: (+356) 343577 Antwerp, Belgium September 26–28, 2001 14th Annual IEEE International ASIC/SOC Conference General Chair: For further information, visit our web site: Joos Vandewalle Dept. of Electrical Engineering, ESAT http://www.imec.be/sips/ “SOC in a Networked World” Univ. of Leuven (K.U.Leuven) September 12–15, 2001 Crystal City Hyatt Hotel Technical Program Co-chairs: http://asic.union.edu Washington, DC Francky Catthoor Marc Moonen IMEC Univ. of Leuven (K.U.Leuven) Leuven, Belgium Leuven, Belgium Driven by the rapid growth of the Internet, communication tech- [email protected] [email protected] nologies, pervasive computing, and consumer electronics, Sys- tems-on-Chip (SoC) have started to become a key issue in today’s electronic industry. The transition from the traditional Application-Specific Integrated Circuit (ASIC) to SoC has led MEMSconference to new challenges in design methods, automation, manufactur- 2001 ing, technology, and test. Conference General Chair Technical Program Chair August 24–26, 2001 Berkeley, California, USA P. R. Mukund John Chickanosky http://www.memsconference.com Rochester Institute of Technology IBM Microelectronics MEMS (Microelectromechanical Systems) Conference 2001 provides a medium for [email protected] [email protected] the exchange of ideas between theoreticians, practitioners, investors, administrators and students to address important issues in the development, manufacturing and commer- For updated conference information, cialization of MEMS. please visit our web site at Further Information http://asic.union.edu For further information either contact [email protected] or see or contact the ASIC/SOC Conference office at 301–527–0900 x104 the conference homepage at: http://www.memsconference.com. IEEE CAS/Notre Dame Workshop European Conference on Circuit Theory and Design, 2001 on Wireless Communications ECCTD’01 and Networking “Circuit Paradigm in the 21st Century” August 9–10, 2001 August 28–31, 2001 Espoo, Finland University of Notre Dame Notre Dame, Indiana Veikko Porra, Conference chair, Helsinki University of Technology (HUT) This workshop is the fourth in the series of IEEE CAS Workshops on Emerging Tech- Olli Simula, Conference co-chair, HUT nologies, which began in Lucerne, Switzerland in June 1998.The workshop also features two keynote speakers, Dr. Jerry Foschini of Bell Laboratories, Lucent Technologies, Inc., FOR MORE INFORMATION PLEASE SEE: http://ecctd01.hut.fi and Dr. Steve Weinstein of NEC-USA. Conference Secretariat: Secretary General: Conference Web site: http://www.nd.edu/~ee/casworkshop Eventra Ltd., [email protected].fi Timo Veijola, ecctd01@hut.fi GENERAL CHAIR TECHNICAL PROGRAM CO-CHAIRS Pietarinkatu 11 B 40 Helsinki University of Technology Yih-Fang Huang, University of Notre Dame Steve Weinstein, NEC-USA FIN-00140 Helsinki, Finland P.O. Box 3000, FIN-02015 HUT, Finland E-mail: [email protected] E-mail: [email protected] Tel. +358 9 6926 949 Tel. +358 9 451 2293 Lang Tong. Cornell University Fax. +358 9 6926 970 Fax. +358 9 451 4818 E-mail [email protected] . . . continued on Page 44 M A IEEE CircuitsCircuits G A and Z and I N SystemsSystems E Volume 1, Number 2, Second Quarter 2001 Filter Banks 4 in Digital Communications by P. P. Vaidyanathan

Efficient and successful communication of messages via imperfect channels is one of the major triumphs of information technology today. With more and more users desiring to share communication channels, the importance of clever exploitation of the bandwidth becomes paramount.

A Microphone Array 26 for Hearing Aids by Bernard Widrow

. . . This method enables the design of highly-directive-hearing in- struments which are comfortable, inconspicuous, and convenient to use. The array provides the user with a dramatic improvement in speech perception over existing hearing aid designs, particularly in the presence of background noise, reverberation, and feedback.

IEEE Circuits and Systems Magazine (ISSN 1531–636X) is published quarterly by the Institute of Electrical and Elec- tronics Engineers, Inc. Headquarters: 3 Park Avenue, 17th Floor, New York, NY, 10016–5997. Responsibility for the con- tents rests upon the authors and not upon the IEEE, the Society, or its members. IEEE Service Center (for orders, subscrip- tions, address changes): 445 Hoes Lane, P.O. Box 1331, Piscataway, NJ 08855–1331. Telephone: +1 732 981 0060, +1 800 678 4333. Individual copies: IEEE members $10.00 (first copy only), nonmembers $20.00 per copy; $7.00 per member per year (included in Society fee) for each member of the IEEE Circuits and Systems Society. Subscription rates available upon request. Copyright and Reprint Permissions: Abstracting is permitted with credit to the source. Libraries are per- mitted to photocopy beyond the limits of the U.S. Copyright law for private use of patrons: 1) those post-1977 articles that carry a code at the bottom of the first page, provided the per-copy fee indicated in the code is paid through the Copyright A Clearance Center, 222 Rosewood Drive, Danvers, MA 01923; and 2) pre-1978 articles without fee. For other copying, reprint, or republication permission, write to: Copyrights and Permissions Department, IEEE Service Center, 445 Hoes Lane, Piscataway, NJ 08855–1331. Copyright © 2001 by the Institute of Electrical and Electronics Engineers, Inc. All rights reserved. Periodicals postage paid at New York, NY, and at additional mailing offices. Postmaster: Send address changes to IEEE Circuits and Systems Magazine, IEEE Operations Center, 445 Hoes Lane, Piscataway, NJ, 08855. Printed in U.S.A.

1 M A IEEE CircuitsCircuits G A and Z and I Columns and Departments SystemsSystems N E 3 From the Editor by Michael K. Sain Editor-in-Chief Michael K. Sain Electrical Engineering Department 33 Transactions University of Notre Dame Guillemin-Cauer Award Notre Dame, IN, USA 46556–5637 Phone: (219) 631–6538 Darlington Award Fax: (219) 631–4393 Outstanding Young Author Coming next issue: E-mail: [email protected] CSVT Best Paper Chaos in Cryptography. Publication Coordinator CAD Best Paper Eric L. Kuehner E-mail:[email protected] VLSI Best Paper VLSI SLIP Special Issue—Call for Papers Editorial Advisory Board Rui J. P. de Figueiredo–Chair 36 Technical Committees University of California, Irvine Nanoelectronics and Gigascale Systems Guanrong Chen City University of Hong Kong by Chung-Yu Wu Wai-Kai Chen University of Illinois, Chicago Leon Chua 39 Meetings University of California, Berkeley Giovanni De Micheli ISCAS 2001 Report Stanford University by Chung-Yu Wu and Michael K. Sain George S. Moschytz Swiss Federal Institute of Technology Mona Zaghloul 40 Society George Washington University Swamy Assumes Publications Post Editorial Board Free Student Membership Promotion Trajkovic Joins Board Orla Feely University College Dublin CASS Awards S.Y. (Ron) Hui Van Valkenburg City University of Hong Kong Ravi P. Ramachandran Technical Achievement Rowan University Education Mehmet Ali Tan Applied Micro Circuits Corporation Meritorious Service Ljiljana Trajkovic Industrial Pioneer Simon Fraser University Chapter-of-the-Year President Names Divisions IEEE Publishing Services Robert Smrek Production Director 41 The ’Umble Ohm IEEE Service Center by Shlomo Karni 445 Hoes Lane P.O. Box 1331 Piscataway, NJ 08855–1331, USA Phone: (732) 562–3944 42 Recognitions 2001 CASS Fellows Frequency of Publication Quarterly 44 Calls for Papers and Participation See also inside front cover and back cover. Magazine Deadlines Final materials for the CAS Magazine must be received Scope: A publication of the IEEE Circuits and Systems Society, IEEE Circuits and Systems Magazine publishes articles presenting novel by the Editor on the following dates: results, special insights, or survey and tutorial treatments on all topics having to do with the theory, analysis, design (computer aided design), and practical implementation of circuits, and the application of circuit theoretic techniques to systems and to signal processing. The cover- age of this field includes the spectrum of activities from, and including, basic scientific theory to industrial applications. Submission of Issue Due Date Manuscripts: Manuscripts should be in English and may be submitted electronically in pdf format to the publications coordinator at First Quarter February 1 [email protected]. They should be double-spaced, 12 point font, and not exceed 20 pages including figures. An IEEE copyright form should be included with the submission. Style Considerations: 1) articles should be readable by the entire CAS membership; Second Quarter May 1 2) average articles will be about ten published pages in length; 3) articles should include efforts to communicate by graphs, diagrams, and Third Quarter August 1 pictures—many authors have begun to make effective use of color, as may be seen in back issues of the Newsletter, available at www.nd.edu/ Fourth Quarter November 1 ~stjoseph/newscas/; 4) equations should be used sparingly.

2 From the Editor Michael K. Sain Editor-in-Chief, IEEE Circuits and Systems Magazine

A Circuits and Systems Carrousel available, and so we will have to rotate through the ow that IEEE Circuits and Systems Magazine topics of interest. So if your topic has not yet ap- Nhas completed its second issue, we are able to peared, please be patient and do not hesitate to send perceive the skeleton begin to form, and we want us your suggestions. to say a few words about the flesh which in due Beyond the articles, which are certainly the most course will fill out the body. obvious of the contents, we have the goal of factor- According to the dictionary, the word carrousel ing in the developments in our Divisions—described has a first meaning of a friendly contest. Here one elsewhere in this issue. These would include our may imagine several groups which are vying for a technical committees, our transactions, our meetings prize. The prize is envisioned as the center of activ- and workshops, and our regional subsocieties. We ity, with the groups arranged around the center. It is have begun this process in our Columns and Depart- a unifying image, with the groups becoming closer ments pages as part of issue one, and we continue it to each other as they move toward the center. The as part of this issue two. Additional plans are being second meaning of carrousel is that of a merry-go- laid for issue three. round, which has multiple interpretations. One Fleshing out the body of the carrousel will of pleasant associated image is that of children riding course require horses and riders. If all those who can, horses mounted on a turning platform. One less do contribute, it will grow into a principal attraction pleasant meaning is that of a humdrum, repetitive on the midway. existence, not necessarily desirable at all. Anyone want a ticket? Even though the second meaning is the more commonly known, we choose the first meaning, because it is a good fit TRANSACTIONS for what the Circuits and Systems Soci- TCAS– I TCAS–VT TCAS– II T–Multimedia ety has in mind for this magazine. Indeed, TCAD T–VLSI EXECUTIVE in the figure on this page, one may see a AWARDS COMMITTEE Circuits and Systems Society Carrousel, with some of the constituent groups rep- resented around the magazine as a focus. IEEE In this simplified diagram, it was not pos- CONFERENCES sible to show every one of the Society en- CIRCUITS BOARD OF ISCAS APCCAS tities; but these are indicative. GOVERNORS AND SYSTEMS ICCAD MWSCAS ICCD ICECS A major goal of the IEEE Circuits and MAGAZINE DAC etc. Systems Magazine is to be a sort of lens at the center of the Society, a lens through which the activities of the various con- STANDING stituencies can be magnified and focussed DISTINGUISHED COMMITTEES for the general membership. Clearly, this LECTURERS TECHNICAL involves the publication of up-to-date ar- COMMITTEES ASP MSA ticles to reach our broad audience. We are CNN&NA VSP CANDE VLSI–S&A already working hard on this aspect, but DSP etc. there is a limit on the amount of space

3 Filter Banks in Digital Communications

by P. P. Vaidyanathan

Abstract—Digital signal processing has played a key role in the de- velopment of telecommunication systems over the last two decades. In recent years digital filter banks have been occupying an increas- ingly important role in both wireless and wireline communication systems. In this paper we review some of these applications of filter banks with special emphasis on discrete multitone modulation which has had an impact on high speed data communication over the twisted pair telephone line. We also review filter bank precoders which have been shown to be important for channel equalization applications.1

1531-636X/06/$10.00©2001IEEE

4 paper we emphasize some of these, Introduction especially the role of filter banks. The fficient and successful communi- aim here is to give an overview of fil- Ecation of messages via imperfect ter banks as applied to digital commu- channels is one of the major triumphs nication systems. Filter banks were of information technology today. With originally proposed for application in

. more and more users desiring to share speech compression more than 25 )

n

( communication channels, the impor- years ago (see references in [4]). To- e tance of clever exploitation of the day they are used for the compression bandwidth becomes paramount. For of image, video, and audio signals, and example telephone lines (twisted-pair the story of their success can be found channels) which were originally in- in many references. More recently fil-

and additive noise

)

is transmitted through a tended to carry speech signals (about ter banks have been used in digital z

)

(

n

( C

4 kHz bandwidth) are today used to communication systems in many x carry several megabits of data per sec- forms. Some of these include perfect ond. This has been possible because of digital transmultiplexing [5, 6], filter- efficient use of high frequency regions bank precoding for channel equaliza- which suffer from a great deal of line tion [7, 8], equalization with fraction- attenuation and noise. As a result of ally spaced sampling [9], and discrete these developments the twisted pair multitone modulation [2, 10, 11]. Fil-

telephone line, which reaches nearly ter banks have been used in high speed Here a sequence of symbols

channel with transfer function every home and office in the western DSL services for internet traffic [3]. Figure 1. A simple model for a digital communication system. world, can today handle high speed They have also been considered for internet traffic as exemplified by popu- blind equalization in wireless channels lar services such as the DSL (digital [12, 13]. x(n) subscriber loop), ADSL, and so forth. As someone pointed out, twisted pair The Noisy Channel copper lines are buried but not dead, In this paper we model the com- thanks to advanced signal processing munication channel as a linear time technology! invariant system with transfer function

detector Communication channels can be C(z) followed by an additive Gaussian wireless or wireline channels, or a noise source e(n) as shown in Fig. 1.

combination of both. In any case they In a digital communication system y(n) introduce linear and nonlinear distor- each sample x(n) comes from a fixed

+

e(n) tions, random noise components, and finite set of values. This set of values noise deterministic interference. The trans- is called a constellation, some standard mission of information with high rate examples being PAM and QAM con- and reliability under such unfavorable stellations explained in Fig. 2. If each

C(z)

conditions has been possible because symbol x(n) is a b-bit number and there channel of fundamental contributions from are fs symbols per second, then the bit many disciplines such as information rate is R = bfs bits per second. For ex- theory, signal processing, linear sys- ample if fs = 1 MHz and b = 4 then tem theory, and mathematics. R = 4 megabits per second (Mbps). x(n) The role of digital signal process- The received signal y(n) is a noisy ing in communication systems has and distorted version of x(n). The de- been quite significant [1, 2, 3]. In this tector at the receiver has to guess the 1 Work supported in parts by the NSF grant MIP value of x(n) based on y(n). The esti- 0703755 and ONR grant N00014-99-1-1002. mated value ˆx(n) belongs to the same

5 (a) 3-bit PAM or 8-PAM –7s –5s –3s –s s3s5s7s

(b) 4-bit QAM or 16-QAM constellation that x(n) came from. If Then the power is the integral of the jω x(n) = ˆx(n) there is no error; otherwise power spectrum Sxx(e ), that is, P = 2π jω the detection is erroneous. In practice ∫0 Sxx(e )dω / 2π. For a given channel there is a nonzero probability of er- the average probability of error Pe at ror Pe in this detection because of the the receiver depends on the transmit- noise e(n) and the intersymbol interfer- ted power P, and the bit rate R . We can ence caused by the channel C(z). The decrease the error probability by trans-

acceptable value of Pe depends on ap- mitting more power. For fixed power, plication. For example it is in the range the error probability increases with bit 10–7 to 10–9 for DSL applications; for rate R . digital speech with mobile phone qual- Note that the power spectrum of

bits ity, larger Pe is acceptable whereas for x(n) tells us how its power is distrib-

= 4 b deep space communication, Pe has to uted in frequency. By carefully shap-

is a quantized

) be quite a bit smaller. ing it, we can increase the achievable

n

( x The model shown in Fig. 1 does rate (for fixed error probability and not work in all situations. For example, transmitted power). The idea is to mobile phone channels are time vary- “pour” more power in the regions ing because of vehicular movement, where the channel gain is large and and a single C(z) cannot be used to rep- noise spectrum is small. This can be resent them successfully. However a pursued in a mathematically rigorous number of practical channels (e.g., the way using fundamentals from infor- wireline telephone channel) can be mation theory [14].

can be regarded as a compex number, taking one

)

n approximated by this. A second remark If the transfer function of the chan-

( x is that the figure implicitly uses dis- nel is known, then the detector can crete time notations. In practice, the equalize it by using the filter 1/C(z). sequence x(n) is converted into a con- This is called the ideal equalizer. It is

tinuous time pulse train ∑nx(n)p(t – nT) also known as the zero-forcing equal-

cient constellations exist [14]. The distance between the constellation fi before imposing it on the channel. The izer for reasons explained in [14]. In output of the channel is sampled and practice 1/C(z) can be approximated digitized to obtain y(n). These details with a stable (possibly FIR) filter. As are not shown in the figure. It should seen from Fig. 3, the effective noise

words (e.g., s in part (a)) can be controlled to control the transmitted power. be understood that C(z) and e(n) are q(n) seen by the receiver is e(n) filtered

possible values from a rectangular constellation as demonstrated in part (b) for the effective discrete time equivalents through 1/C(z). This has the effective

b 2 of the actual channel parameters. noise power spectrum

of

real number as demonstrated in part (a) for 3-bit PAM (also called 8-PAM). (b) For the case of

Figure 2. (a) For the case of Pulse Amplitude Modulation (PAM), the sample

(called 16-QAM). More ef

Quadrature Amplitude Modulation (QAM)

Power Allocation and jω ω ∆ S (e ) j = ee Water-Filling Strategy Sqq (e ) ω 2 . (1) C(e j ) The transmitted signal power P is proportional to the mean square value This ratio summarizes the channel of x(n). Assume that x(n) is a wide completely for the purpose under dis- sense stationary random process [14]. cussion. If C(z) has zeros close to the

6 y(n) x(n) C(z) + 1/C(z) detector x(n)

channel equalizer e(n) Figure 3. A channel, equalized with an ideal equalizer 1/C(z) (also called the zero-forcing equalizer). The effective noise q(n) as seen by noise the detector is nothing but e(n) filtered through 1/C(z).

jω unit circle, then 1/C(z) has poles near Sxx(e ) (total power) increases with the unit circle and the noise gain can choice of λ. The choice of λ therefore be large. In frequency regions where depends on the total available power P. jω the ratio Sqq(e ) is small, we should For fixed total power P and fixed allocate more power. In fact the opti- channel, the capacity C is the maxi- jω Filter Banks mal power distribution Sxx(e ) can be mum rate at which information can be described precisely with the help of transmitted with arbitrarily small error in Digital Fig. 4. This figure says that probability. This should not be con- Communications fused with the actual bit rate R with  λ−S (e jω ) when ≥ 0 jω = qq nonzero error probality Pe. Evidently Sxx (e )  (2)  0 otherwise if Pe is allowed to be large enough then R can be made larger than the capac- Filter Banks where λ is a constant. That is, the trans- ity C. In situations where Pe is required in Digital mitted power at a frequency should be to be reasonably small, C can be regarded jω Communications equal to the gap between λ and Sqq(e ). as a useful bound on acheivable rate R . Notice that in regions where the chan- How do we shape the power spec- jω jω nel is “too bad” (Sqq(e ) > λ) we trans- trum Sxx(e ) to satisfy the water-filling mit no power at all. If we imagine a type of power allocation? This is tricky bowl with its bottom shaped like because we do not have a great deal of ω ω Filter Banks S (ej ), then S (ej ) is the height of freedom to shape things, especially if qq xx in Digital water filling the bowl, with λ denoting x(n) is user generated data! For ex- Communications the uniform water level everywhere. ample x(n) could be binary data or data This is a classic result on power allo- from a PAM or QAM constellation, cation, and is called the water filling and might behave like an iid (indepen- rule [14]. It is clear that the area under dent identically distributed) sequence

ω jω S (e j ) is zero here S xx (e ) xx

λ water level

jω S qq (e ) defines bottom of bowl

ω 0 2π

jω jω Figure 4. Relation between the input power spectrum Sxx(e ) and the effective noise power spectrum Sqq(e ) to achieve maximum rate for fixed total power. This is called the water filling rule.

7 Filter Banks in Digital Communications

x(n) M y(n) x(n) M y(n)

M-fold expander M-fold decimator

x(n) x(n)

n n –1 0123 –3 0123 6 9

example with M=3 example with M=3 y(n) y(n)

n n –1 0123 6 9 –1 0123

Figure 5. The M-fold expander merely inserts M – 1 zeros between adjacent samples, as demonstrated in the figure for M = 3. The M-fold decimator has input-output relation y(n) = x(Mn). Thus, only a subset of input samples are retained. This is demonstrated in the figure for M = 3. Note that the samples automatically get renumbered so that y(1) = x(M), y(2) = x(2M) and so forth.

v0 (n) M v(n) v(n) M v0 (n) z–1 z v (n) 1 M M v1 (n) z–1 z interleaving deinterleaving or unblocking or blocking z–1 z (a) vM – 1(n) M (b) M vM – 1(n)

v (1) v (0) 1 1 v (1) v2 (0) 2 v2 (2) v0 (0) signal v(n) v0 (2) v0 (1) v1 (2) (c) n –1 0123456 789

Figure 6. (a), (b) The operations of interleaving and deinterleaving using multirate building block notation. (c) Demonstration for M = 3. The interleaving operation is also called unblocking. Similarly deinterleaving is also called blocking.

8 with a practically flat power spectrum. M = 3. It is clear that we can regard ω A clever way to approximate optimal v(n) as a time-domain multiplexed or 0 π F power allocation would be to divide TDM version of the individual signals 2 the channel bandwidth into several vk(n). The components vk(n) are also subbands and transmit in each subband called the polyphase components of 3 channel separately [15, 16]. This al- v(n) with respect to M [17]. F

ready suggests a resemblance to fre- 2 quency division multiplexing but the The Digital Transmultiplexer F main difference now is that the differ- Our next step would be to describe 1 ent subband channels carry different the operation of a filter bank structure F effective noise psd parts of a single input stream. How this called the transmultiplexer [18, 5, 17]. 0 is actually accomplished will be ex- It was originally intended to convert F plained in the section Discrete data between time division multi- 0

Multitone Modulation (DMT). plexed (TDM) format and frequency (c) division multiplexed (FDM) format. Decimators, Expanders, Figure 7(a) shows a schematic of this and Multiplexers in all-discrete language. Here Fk(z) are Before proceeding further we re- called transmitting filters or interpola- view a few standard building blocks tion filters. The kth transmitting filter M), scaled – 0 and terminology used in multirate sig- has output f (n nal processing. Most of the details can ∞ 0 f (n)

= − M be found in [17]. The building blocks uk (n) ∑ xk (i) f k(n iM) . (3) ↑ M and ↓ M, shown in Fig. 5, are i=−∞ 0 x (n)

called the expander and decimator re- Figure 7(b) demonstrates how this 0 spectively. Their operation is ex- construction is done for the 0th filter plained in the figure caption. F0(z), assumed to be lowpass. Essen- Two standard operations called tially we draw one copy of the impulse . “blocking” and “interleaving” often response sequence f0( ) around every 1) (b) – arise in communication systems that sample of x0(n) (separated by M) and 0 x (n use filter banks. It is sometimes con- add them up. Thus uk(n) is an interpo- venient to explain these using multirate lated version of xk(n) and has M times building blocks. These are shown in higher rate. The outputs of the filters

Figs. 6(a) and 6(b). The connection F (z), F (z) and so forth are more com- lters). 1 2 fi between the signal v(n) and the plicated waveforms because they are jω “blocked components” vk(n) is indi- bandpass. The filters {Fk(e )} tradi-

nite order

cated in Fig. 6(c) using the example of tionally cover different uniform re- fi

are shown in (b).

)

n

(

0

f

y (n) x (n) u (n) y(n) 0 0 0 H (z) M F0 (z) 0 M

, and (c) frequency responses of the y1 (n) ) x (n) u 1(n) z

1 ( M F (z) H (z) M 0

1 1 F

lter

(a) lters (assumed to be ideal in

channel fi transfer func. u (n) y (n) xM – 1(n) M–1 x(n) M–1 H (z) M F M – 1 (z) C(z) + M–1 M expanders Transmitting Receiving decimators

Only the envelope of the samples of filters filters interpolation

transmitting

channel Figure 7. (a) The digital transmultiplexer, (b) operation of the modulation noise e(n) detected symbol symbol

9 gions of frequency as shown in Fig. The receiving filter bank {Hk(z)}

7(c). The signals uk(n) are analogous separates the signal y(n) into the com-

to modulated versions of the ponents yk(n) which are distorted and

“baseband” sequence xk(n) be- noisy versions of the symbols xk(n). cause the bandwidth is shifted The task at this point is to detect the

to the passband of Fk(z). symbols xk(n) from yk(n) with accept- These are packed into M ad- able error probability. Thus, even

jacent frequency bands though xk(n) is interpolated to get uk(n),

(passbands of the filters) it is not necessary to ensure uk(Mn) =

and added to obtain the xk(n); the crucial issue is to make yk(n)

composite signal x(n). resemble xk(n).

With the filters Fk(z) cho- sen as good bandpass Discrete Multitone Modulation filters, we can regard x(n) (DMT) as a frequency division The transmultiplexer configuration multiplexed or FDM ver- is used in another scheme called dis- sion of the separate signals crete multitone modulation or DMT.

xk(n). By contrast, if Fk(z) are The main difference is in the interpre- –k just delay elements z , then tation of the signals x(n) and xk(n). To the transmitter part is similar to explain this consider Fig. 8(a) which Fig. 6(a) and x(n) is a time-multi- shows the first stage of multitone

plexed version of the M signals xk(n). modulation [15] called the parsing Letting T denote the spacing between stage. Here s(n) represents binary samples of x(n), we see that the data to be transmitted over a channel.

samples of xk(n) for any given k are This data is divided into nonover- separated by a longer duration of MT lapping b-bit blocks. The b bits in each seconds. block are partitioned into M groups,

b 0 b1 b2 b 0 b1 b2 b 0 b1 b2 binary data, s(n)

b-bit block b-bit block b-bit block

(a)

x (n) 0 b0 bits

serial x 1(n) b1 bits s(n) to parallel

x (n) 2 b 2 bits (b)

Figure 8. Parsing stage in discrete multitone (DMT) modulation. A binary stream s(n) is subdivided

into b-bit blocks and each block is subdivided into M subgroups. The kth subgroup defines a bk-bit

symbol xk(n). Note that xk(n) and xk(n + 1) are separated by bT seconds if s(n) and s(n + 1) are separated by T seconds.

10 the kth group being a collection of Biorthogonality and bk bits (demonstrated in the figure for Perfect DMT Systems M = 3). Thus the total number of bits Consider a linear time b per block can be expressed as invariant system with im-

M −1 pulse response h(n) and = transfer function H(z) b ∑bk . (4) k= 0 sandwiched between an expander and decimator The bk bits in the kth group constitute (Fig. 9(a)). It can be shown the kth symbol xk which can therefore that this is equivalent to a be regarded as a bk-bit number. For the linear time invariant sys- nth block, this symbol is denoted as tem with decimated im- xk(n). This is the modulation symbol pulse response g(n) = for the kth band. The collection of h(Mn). In z-transform notation we de- symbols {x0(n), x1(n), … , xM–1(n)} is note this as G(z) = [H(z)]↓M. together referred to as the DMT sym- Using this simple idea, we can bol. The sample xk(n) is typically a understand the operation of the DMT PAM or a QAM symbol (Fig. 2). The system in a very effective way. Thus transmitting filters fk(n) create the the transfer function Dkm(z) from M-fold higher rate signals uk(n) as be- . . xm( ) to yk( ) in Fig. 7(a) is the deci- fore, which are then added to produce mated version of the product-filter the composite signal x(n). In this way, H (z)C(z)F (z). If this is nonzero for various parts of the original binary k m m ≠ k then the symbol yk(n) is affected message s(n) are packed into different by x (i) resulting in interband inter- frequency regions allowed by the m ference. Similarly if Dkk(z) is not a con- channel [15, 16]. stant then y (n) is affected by x (i), Notice that for a given constella- k k tion, the power can be increased or decreased by scaling the distance be- tween the codewords (e.g., by adjust- ing s for the PAM constellation in Fig. (a) x(n) M H(z) M y(n) 2(a)). We therefore have the freedom to allocate different powers for differ- ent subband channels. In this way the classical water-filling rule can be ap- (b) x(n) G(z) y(n) proximated. For a given transmitted power and probability of error, multitone g(n) = h(Mn) orG(z) = [H(z)] modulation yields better bit rate than M single tone modulation (M = 1 case), assuming no channel coding. Figure 9. A transfer function H(z) sandwiched between an Background material on the DMT expander and a decimator is equivalent to another transfer system and more generally on the use function G(z) with impulse response g(n) = h(Mn). of digital filter banks in communica- tions can be found in [6, 19, 16, 10]. The DMT idea is similar in principle i ≠ n, due to the filtering effect of to subband coding [4, 20, 17] where Dkk(z). This is called intraband inter- a signal x(n) to be quantized is first de- ference. If interband and intraband in- composed into subbands. terferences are eliminated, the DMT

11 u (n) y0 (n) x 0(n) 0 y(n) M F0 (z) H0 (z) M

u (n) y1 (n) x1 (n) 1 M F1 (z) H1 (z) M

x (n) y (n) M–1 u M – 1(n) x(n) M–1 M F M – 1(z) C(z) + 1/C(z) HM – 1(z) M expanders Transmitting channel equalizer Receiving decimators filters filters e(n) Figure 10. The DMT system with ideal channel equalizer 1/C(z). This system has the perfect symbol recovery (PR) property

in absence of noise if the filter bank {Hk, Fm} is biorthogonal. See text.

system is said to be free from ideal filters having overlapping re- intersymbol interference (ISI). If we sponses. This idea goes back to early assume that the filters are ideal work on transmultiplexers [5] and is nonoverlapping bandpass filters related to the notion of a biorthogonal stacked as in Fig. 7(c), then there is no filter bank. To explain this, consider interband interference. Furthermore, what happens when the channel intro- suppose that C(z) is completely equal- duces no distortion (C(z) =1). Under ized with the inverse filter 1/C(z) as this condition we have perfect symbol shown in Fig. 10. Then the system is recovery if and only if the transmitting

ISI free, and yk(n) = xk(n) for all k (in and receiving filters satisfy the absence of noise), and we have the per- biorthogonality property defined as fect symbol recovery or PR property. δ Ideal nonoverlapping filters are of Hk(z)Fm(z) ↓M = (k – m) . (5) course unrealizable, and good approxi- This means that the impulse re- mations of such filters are expensive. sponse gkm(n) of the product filter It turns out that perfect symbol recov- ∆ G (z)= H (z)F (z) has the Nyquist(M) ery can be obtained even with non- km k m or zero-crossing property g (Mn) = 0 (6) g (n) km km for k ≠ m and g (Mn) = δ(n) as dem- 1 2 n kk onstrated in Fig. 11 for M =3. This –6 –360 3 condition is readily achieved with careful design of filters. For example, it is possile to design FIR biorthogonal g (n) filter banks with almost any filter length. kk In this paper we shall make the

simplifying assumption that {Fm, Hk} n is biorthogonal (i.e., Eq. (5) holds) and –6 –360213 that the channel transfer function C(z) Figure 11. In a biorthogonal transmultiplexer, we have the perfect symbol recovery is equalized by using the inverse filter property in absence of channel distortion C(z) and channel noise e(n). This is or zero-forcing equalizer 1/C(z) just achieved by constraining the filters such that the products Gkm(z) have the Nyquist before entering the bank of filters (M) property. That is, their impulse responses have the zero crossing property at integer multiples of M as demonstrated above for M = 3. {Hk(z)}. In a biorthogonal DMT sys-

12 tem with zero-forcing equalizer, the expressed in terms of the signal power σ 2 only remaining distortion is due to the Pk, noise variance qk, and number of channel noise. The received symbol bits bk. The exact expression can be can be written as found in many references (e.g., see [14, 11, 21]). The main point is that this y (n) = x (n) + q (n) (7) k k k expression can be inverted to obtain where qk(n) is the channel noise the total power in the symbols xk(n).

filtered through Hk(z)/C(z) and deci- The result takes the form [11, 21] mated (Fig. 12(a)). Thus the variance M −1 M −1 of qk(n) can be calculated using the P = ∑ P =∑β()P (k),b ×σ2 equivalent circuit shown in Fig. 12(b) k e k qk (8) k= 0 k= 0 jω jω j 2 where Sqq(e ) = See(e )/|C(e )| is the equivalent noise spectrum defined in where the exact nature of the function the section, Power Allocation and β(., .) is not of immediate interest. This Water-Filling Strategy. This figure expression says that if the acceptable shows an analysis bank {Hk(z)} probabilities of error at the bit rates whose input is a noise source with ef- {bk} are {Pe(k)}, then the total power fective power spectrum given by (1). P has to be at least as large as the right hand side of (8). If we try to decrease Optimization of Pe(k) for a given bit rate, we need more DMT Filter Banks power. In this section we discuss the opti- The crucial point to note here is mization of filter banks used in DMT that the power P can be minimized by σ 2 systems. The variance of the symbol carefully controlling the variances qk xk(n) in Fig. 7(a) represents its aver- of the noise components qk(n) at the age power Pk. For simplicity assume detector inputs. Given the channel C(z) that xk(n) comes from a bk-bit PAM and the channel noise spectrum See(z), constellation (Fig. 2) with equal prob- the only freedom we have in order to σ 2 ability for all codewords. Assume fur- control qk is the choice of the filters ther that the noise qk(n) is Gaussian Hk(z) (see Fig. 12(b)). But we have to σ 2 control these filters under the con- with variance qk. Then the probabil- ity of error in detecting xk(n) can be straint that {Hk, Fm} is biorthogonal.

q (n) q(n) H0 (z) M 0 spectrum q (n) S (z) k qq H (z) q (n) e(n) 1 M 1 H k(z)/C(z) M channel noise at noise at noise detector (a) detector (b) q (n) HM – 1(z) M M–1 analysis bank Figure 12. (a) The effect of the channel noise at the detector inputs can be modelled as the filtered version of the channel noise. (b) The noise variances at the detector inputs can therefore be calculated as the subband variances in a maximally decimated analysis filter bank, jω jω jω 2 whose input has the power spectrum Sqq(e ) = See(e )/ |C(e )| .

13 1/ M f 0 (n) (a) n

01 M –1 Figure 13. An example of the uniform DFT filter bank. (a) Impulse response of F0(z) and (b) magnitude jω responses of the digital filters Fk(e ). Each filter response is a frequency-shifted version of the preceding F F F filter. The shifts are in uniform increments of 2π/M 0 1 2 where M is the number of filters. If the peak passband response is normalized to 0 dB, the minimum stop band attenuation is about –13 dB.

13 dB (b)

ω

π π 0 2 /M 4 /M 2 π

Since the scaled system {αkHk, Fm /αm} … fk(n + M), fk(n),

is also biorthogonal (as we can show fk(n – M), fk(n – 2M) … (9) using (5)), it appears that the variances σ 2 qk can be made arbitrarily small by covering the kth frequency band. The making αk small. The catch is that the basis has an infinite number of ele-

transmitting filters Fm(z)/αm will have ments, each element being a filter ob- correspondingly larger energy which tained from the preceding element by means an increase in the power actu- a time-shift of M samples. The com- ally fed into the channel. One correct posite signal x(n) which enters the approach to do the optimization would channel is therefore a linear combina- be to impose a power constraint. tion of the basis functions from all the Mathematically this is trickier than channels. We say that a set of M filters

constraining the powers Pk in the sym- {Fk(z)} is orthonormal if these basis

bols xk(n). In the next section we functions are orthogonal to each other, confine the optimization to a class of and each of them is normalized to have filter banks called orthonormal filter unit energy. For perfect symbol recov- banks. In this case the optimization ery (or biorthogonality), the transmit- problem is especially easy to formu- ting and receiving filters in any or- late, and elegant solutions can be found thonormal filter bank are related by

as well. * hk(n) = ƒk(– n) (10) Orthonormal DMT Systems which is called time reversed-conju- Recall from Fig. 7(a) that the gation. This condition means in par-

subband channel signals uk(n) are the ticular that the transmitting and receiv- outputs of interpolation filters, and can ing filters have identical frequency re- be expressed as in Eq. (3). We can re- sponse magnitudes. Orthonormal filter

gard the subchannel signal uk(n) as banks have been extensively studied belonging to a subspace spanned by and documented, see for example [17], the basis functions [20] and references therein. It is pos-

14 Filter Banks in Digital Communications sible to have orthonormal filter banks DFT matrix (with proper normaliza- where the filters are FIR. An example tion) is unitary. The DFT filter bank is is the filter bank where f0(n) is chosen used widely in DMT systems [19] for as a rectangular pulse of length M and certain types of DSL services. fk(n) are the modulated versions The popularity of the DFT based ω DMT filter bank arises from the fact f (n) = f (n)ej kn (11) k 0 that if M is chosen as a power of two with ωk = 2πk / M representing the kth (e.g., M = 512 which is typical) the center-frequency. See Fig. 13. The fre- DFT can be implemented very effi- quency responses are uniformly ciently using the fast Fourier transform jω shifted versions of F0(e ) as shown in (FFT) algorithm. By using bit alloca- the figure. tion in the transform domain, the DFT This is called the DFT filter bank based DMT system can take advantage because it can be implemented with a of the shape of the effective noise spec- DFT matrix and an inverse DFT trum (1) and obtain a performance (IDFT) matrix as shown in Fig. 14. At close to the water-filling ideal (section each instant of time n, the DMT sym- Power Allocation and Water-Filling bol {x0(n), x1(n), … xM–1(n)} is trans- Strategy). In the Examples section we formed into the IDFT domain. The shall provide numerical examples in components vk(n) of the resulting sym- terms of bit rates and transmitted power. bol are interleaved to obtain the chan- nel signal x(n). At the receiver the sig- Optimal Orthonormal nal is de-interleaved and a DFT is per- DMT Systems formed. The results yk(n) are noisy ver- In an orthonormal DMT filter bank sions of the transmitted symbols xk(n). the transmitting and receiving filters Orthonormality of the basis functions have unit energy. So we cannot insert in (9) follows from the fact that the arbitrary scale factors in front of Hk(z)

v (n) y (n) x 0(n) 0 x(n) 0 M C(z) + 1/C(z) M v (n) z –1 x 1(n) 1 z y 1(n) M channel equalizer M e(n) z –1 IDFT z DFT

v (n) z –1 y (n) x M – 1 (n) M–1 z M –1 M M

interleaving deinterleaving Figure 14. DMT system based on the uniform DFT filter bank. The channel equalizer 1/C(z) can be approximated well in many ways. An indirect but effective way to perform equalization is to introduce redundancies such as the cyclic prefix [22]. 15 q(n) M q 0 (n) effective z q (n) Figure 15. Noise model for the DMT filter bank, redrawn in noise M 1 terms of the transform matrix T. If T is the DFT matrix, then z this is the noise model for the DFT based DMT. If T is the KLT matrix for the effective noise q(n) then this corresponds T to an optimal DMT system. See text. z

M q M – 1(n)

to reduce the noise at the detector in- put. Moreover, orthonormality also im- KLT Based DMT Systems plies that the average variance of the Consider for example the class of composite signal x(n) is the average of FIR orthonormal DMT systems, that

the variances of the symbols xk(n). is, systems where Fk(z) and Hk(z) are That is, the actual power entering the FIR. Assume further that the filter channel is proportional to the sum of lengths are constrained to be no larger

powers Pk in the symbols xk(n). Refer than M. An example is the DFT based again to Fig. 12(b) now. For a given DMT system described in the section, channel the effective noise spectrum Orthonormal DMT Systems, with jω Sqq(e ) is fixed (ratio defined in (1)). transmitting filters as in (11). In view Assume further that the integer M of the time-reversed conjugation prop-

(number of subchannels) is fixed. For erty (10), the receiving filters hk(n) are jω n — a given set of error probabilities and bit given by e k /√M for –(M – 1) ≤ n ≤ 0, rates, the required transmitted power and by 0 otherwise. depends only on the noise variances Since the channel noise is filtered σ 2 qk as shown by Eq. (8).We have to by the DFT matrix, the noise model for find an orthonormal filter bank such this system can be drawn as in Fig. 15 that this power is minimized. This is where T represents the DFT matrix. the problem of designing an optimal This is precisely the structure of the orthonormal DMT system. It turns out noise model of Fig. 12(b), drawn in a

that the solution {Hk} depends only on different way. In general any pair of

the effective noise spectrum and not on receiver noise components qk(n) and

the desired values of error probabili- qm(n) have some statistical correlation ties and bit rates. between them. But it is possible to re-

P. P. Vaidyanathan received the B.Sc. (Hons.) degree in physics and the B.Tech. and M.Tech. degrees in radiophysics and electronics, all from the University of Calcutta, India, in 1974, 1977 and 1979, respectively, and the Ph.D degree in electrical and computer engineering from the University of California at Santa Barbara in 1982. In March 1983 he joined the Department of Electrical Engineering at the Calfornia Institute of Technology as assistant professor, and since 1993 has been professor of electrical engineering there. His main research interests are in digital signal processing, multirate systems, wavelet transforms and signal processing for digital communications. Dr. Vaidyanathan served as vice-chairman of the technical program committee for the 1983 IEEE ISCAS, and as technical program chair for the 1992 IEEE ISCAS. He was associate editor for the IEEE Transactions on Circuits and Systems from 1985-1987, and is currently associate editor for the journal IEEE Signal Processing Letters, and consulting editor for the journal Applied and Computational Harmonic Analysis. He was guest editor in 1998 for special issues of the IEEE Transactions on Signal Processing and the IEEE Transactions on Circuits and Systems–II. Dr. Vaidyanathan is the author of the book Multirate Systems and Filter Banks. He was three times a recipient of the Award for Excellence in Teaching at the California Institute of Technology. He also received the NSF’s Presidential Young Investigator award in 1986. In 1989 he received the IEEE ASSP Senior Award for his paper on multirate perfect-reconstruction filter banks. In 1990 he was recipi- ent of the S. K. Mitra Memorial Award from the Institute of Electronics and Telecommunications Engineers, India, for his joint paper in the IETE journal. Dr. Vaidyanathan was elected Fellow of the IEEE in 1991. He received the 1995 F. E. Terman Award of the American Society for Engineering Education, sponsored by Hewlett Packard Co. In 1999 he was chosen to receive the IEEE CAS Society’s Golden Jubilee Medal.

16 pcfb Figure 16. This figure schematically explains what a principal component filter bank (PCFB) does. The dark other blue columns represent the partial sums of subband variances of a PCFB in a class C of orthonormal filter banks, and the light blue columns represent the same sums for an arbitrary filter bank in C. By definition, the PCFB partial sums always dominate. Remember here that the sum of all subband variances is the same for all orthonormal filter banks, and is equal to M times subband 0 subbands subbands all subbands the input variance. 0 and 1 0, 1 and 2 place the DFT matrix with another bank or PCFB for the power spectrum jω unitary matrix T such that qk(n) and Sqq(e ), as shown in [11, 21]. qm(n) are uncorrelated for all n when To explain what a PCFB is, assume k ≠ m. Such a decorrelating matrix T that we are given a class C of M chan- depends only on the power spectrum nel orthonormal analysis filter banks as of the effective noise q(n). It is called in Fig. 12(b). Given an input power jω the M x M KLT matrix for q(n). Essen- spectrum Sqq(e ), a PCFB in C is a filter tially it is a unitary matrix which di- bank such that the output variances σ 2 agonalizes the M x M autocorrelation qk, arranged in nonincreasing order σ 2 ≥ σ 2 matrix of q(n). It can be shown that if qk qk + 1 have a very special prop- T is chosen as the KLT matrix (and its erty. Namely, inverse used in the transmitter) then the partial sums the required power P is minimized. of these vari- If the transmitting and re- The KLT offers the optimal DMT so- ances, lution if we compare all orthonormal σ 2 σ 2 σ 2 ceiving filters are allowed to q0, q0 + q1, DMT systems with FIR filters of 2 2 2 σ q + σ q + σ q , have infinite length with no length ≤ M. The proof follows as a spe- 0 1 2 … (12) cial case of the results in [11, 21]. causality restrictions, then are larger than nonoverlapping brickwall Optimal Orthonormal DMT Systems the correspond- Using Unconstrained Filters ing partial sums filters are allowed. If the transmitting and receiving for any other filters are allowed to have infinite filter bank in length with no causality restrictions, this class.2 This idea is demonstrated then nonoverlapping brickwall filters in Fig. 16. If C is the class of orthonor- are allowed. In fact ideal filters with mal filter banks with filter length ≤ M multiple passbands are allowed as well then the KLT of the input is a PCFB. and could be useful as we shall see If C represents the class of ideal or- later. Assuming orthonormality, the thonormal filter banks (with filters al- transmitting filters are constrained as lowed to be noncausal with unre- jω * jω Fk(e ) = Hk(e ). What is the best stricted lengths) then there is a system- choice of the frequency responses of atic way to construct the PCFB [21]. the receiving filters {Hk(z)} if we wish Unlike the brickwall filter bank of Fig. to minimize transmitted power? The 7(c), each filter can have multiple pass- answer again depends only on the ef- bands. Thus, the PCFB partitions the jω fective noise spectrum Sqq(e ). In fact 2 Readers familiar with singular value decomposi- the optimal choice of {Hk(z)} is the so- tion and principal component reconstruction will called principal component filter notice an analogy.

17 frequency domain in a different way ued. Assume M = 2 (two band DMT). according to the nature of the input The two filters of the PCFB for the spectrum. above power spectrum are shown in For an arbitrary class C such as the Fig. 17(b), and the traditional class of FIR orthonormal filter banks brickwall filter bank response is shown with length constrained by some inte- in Fig. 17(d). For purpose of calcula- ger, the PCFB may not in general ex- tion assume that the desired probabil- ist. The detailed theory of the PCFB is ity of error is 10–9 in each band and that available in [23], and a tutorial review the number of bits per symbol in the

can be found in [21]. two PAM constellations are b0 = 6 and

Assume that the error probabilities b1 = 2. If the sampling rate is 2 MHz and total allowed power are fixed. It this implies a bit rate of 8 Mbits/sec. can then be shown that the bit rate, The average power needed can be

which is proportional to ∑kbk, is maxi- found from (8). It turns out that the mized by the PCFB. Similarly, with power required by the PCFB is nearly appropriate theoretical modelling the 10 times smaller than the power re- information capacity (for fixed total quired by the brickwall filter bank. For power) is also maximized by the PCFB. example if the brickwall system re- Details can be found in [24, 11, 25]. quires 56 milliwatts, then the PCFB uses only about 5.67 milliwatts for the Examples same performance! For DMT systems Assume that the effective noise with larger number of bands, the dif- power spectrum has the hypothetical ference is less dramatic. In fact for form shown in Fig. 17(a). We have very large M the performances are shown only the region 0 ≤ ω ≤ π be- nearly identical [25]. cause we assume in this example that It turns out that for a monotone all time domain quantities are real val- decreasing (or increasing) power spec-

(a) 1000 (c) 1000 500 500 effective effective noise power 1 noise power 1 spectrum ω spectrum ω 0 π/4π/2 3π/4 π 0 π/4π/2 3π/4 π

H 0 (d) ω H0 H1 (b) 0 π/4 3π/4 π brickwall PCFB FB ω filters 0 π/4π/2 3π/4 π H1 ω 0 π/4 3π/4 π

Figure 17. An example showing the difference between brickwall stacking and PCFB stacking. Here M = 2. The PCFB tends to distribute the variances in such a way that one subband variance is maximized and the other is minimized. Each filter in a PCFB can have more than one passband. See text.

18 Figure 18. The channel gain in a twisted pair copper line is a decaying function, with severe attenuation at high frequencies. Moreover, there are sharp dips at various frequencies due to the presence of bridged taps.

Channel gain

bridged tap 2

bridged tap 1 f

2 MHz trum the ideal PCFB is precisely the transmission of baseband speech brickwall filter bank. For a power spec- (about 4 kHz bandwidth) more than trum with many variations, especially 100 years ago, the twisted pair copper bumps and dips as in Fig. 17(a) the wire has therefore come a long way PCFB has filters with many passbands in terms of and its performance differs significantly bandwidth from the brickwall filter bank as dem- utilization Originally intended for trans- onstrated in the preceding example. and commer- A practical example is the effective cial applica- mission of baseband speech noise spectrum in a twisted pair chan- tion. This has (about 4 kHz bandwidth) more nel used for ADSL (asymmetric DSL) given rise to service. The twisted pair is a pair of the popular than 100 years ago, the twisted insulated copper wires that are twisted saying that pair copper wire has therefore at periodic intervals.3 It reaches every the DSL tech- home in the world that has a wireline nology turns come a long way in terms of telephone. The channel gain of the copper into bandwidth utilization and com- twisted pair decays rapidly with fre- gold. quency and wirelength. Figure 18 It is typi- mercial application. This has shows a typical qualitative example. cal to have 50 given rise to the popular say- The two dips in the figure are created twisted pairs by bridged taps which are used in the bundled into ing that the DSL technology to provide additional ser- one cable for turns copper into gold. vice flexibility [2]. In spite of the large several attenuation, it is still possible to use the kilofeet. As a twisted pair over a large frequency result, the most dominant noise is the range (up to a few MHz) and achieve interference created by services from high rates. Several types of DSL ser- other cables. Two kinds of such inter- vices are based on exploiting its band- ference can be distinguished, namely width like this. Originally intended for near end cross talk abbreviated as NEXT, and far end cross talk abbre- 3 The idea of twisting originated from Alexander viated as FEXT. Essentially, NEXT is Graham Bell who invented it around 1880, with a view to cancelling the effect of electromagnetic the cross-talk from a transmitter at the interference. same side of the cable whereas FEXT

19 Noise power –35 dBm/Hz spectrum –80 HAM NEXT –80 –120 FEXT AM f

25 kHz 138 kHz 2 MHz

600 kHz 1.9 MHz Figure 19. Various components contributing to the noise spectrum in the ADSL service offered on a twisted pair line. The figure

shows the noise power spectra in milliwatts per hertz on a dB scale. The unit dBm stands for 10 log10(mW).

jω jω 2 is the cross-talk from a transmitter at trum See(e )/| C(e )| has several the other end of the cable. The statis- bumps and dips. A PCFB is therefore tics of these have been studied for significantly different from the many years both theoretically and by brickwall filter bank. A detailed ex- extensive measurements [2]. Figure 19 ample presented in [25] for typical shows typical power spectra of the ADSL downstream service shows that NEXT and FEXT noise sources in a the difference in performance can be 50-pair cable. In addition to the NEXT significant for a small number of and FEXT, DSL services also suffer subchannels M. For example, assume from AM radio interference and ama- M = 8 and probability of error 10–9 in teur radio (HAM) interference as each subchannel. With typical num- shown schematically in the figure. The bers chosen for various noise compo- main point of this discussion is that the nents for the ADSL downstream sce- total noise spectrum is quite compli- nario, the required power can be calcu- cated and is far from a constant or a lated. For an overall bit rate of 3.2 Mb/s monotone decreasing function. And it is verified in [25] that the required since the channel gain has dips due to power is 4.68 mW for traditional DFT bridge taps, the effective noise spec- type of multitone, and only 0.94 mW for

v (n) y (n) x 0 (n) 0 x(n) 0 N C(z) + D(z) N v (n) z –1 x1 (n) 1 z y1 (n) N channel equalizer N e(n) z –1 z R(z) E (z)

y (n) x M – 1 (n) M –1

–1 vN – 1(n) z z N N

interleaving deinterleaving

Figure 20. A DMT system with redundancy.

20 Filter Banks in Digital Communications

PCFB using ideal filters.The interme- outputs of R(z). The expanders ↑ N diate value of 2.76 mW is achieved and the set of delay elements follow- when traditional DFT is replaced by ing them simply interleave the outputs the KLT. Finally, a traditional ideal of R(z) to produce the composite chan- brickwall filter bank uses 1.28 mW, nel signal x(n) (see Fig. 6(a)). Using slightly worse than the ideal PCFB. standard multirate identities [17] we can draw the system of Fig. 20 in terms Filter Banks with Redundancy of transmitting and receiving filters as We now mention some generaliza- shown in Fig. 21. If the samples of tions of the DMT structure that have xk(n) are separated by 1 second, for received significant attention recently, example, then the samples of x(n) are especially in the signal processing separated by 1/N seconds, instead of community. First consider Fig. 20 and 1/M seconds as before. This introduces

u 0 (n) y 0 (n) x 0 (n) y(n) N F0 (z) H 0 (z) N

y 1 (n) x (n) u 1 (n) 1 N>M N F 1 (z) H 1 (z) N

x (n) y M–1 (n) M–1 u M–1(n) x(n) F (z) H (z) N M–1 C(z) + D(z) M–1 N Transmitter Receiving filters filters e(n) Figure 21. Redrawing the DMT system with redundancy. The only difference from Fig. 7(a) is that the expander ratio is N which is larger than M. compare with the DFT based DMT redundancy because the actual symbol system of Fig. 14. The DFT and IDFT rate of M per second has been increased matrices have been replaced with the to N per second. The factor N/M is matrices R(z) and E(z) which can de- called the bandwidth expansion factor. pend on z. This means that the filters There are good practical reasons

Fk(z) and Hk(z) can have arbitrarily for the incorporation of such redun- large orders. This freedom can be ex- dancy. For example channel equaliza- ploited to design DMT systems with tion is easier [3]; there is no need to better performance (e.g., smaller total directly approximate the channel in- power for a given set of requirements). verse 1/C(z), which is undesirable Second, and more importantly, when C(z) is a filter of high order. In a there is a new integer N > M in the fig- DMT system with redundancy, it is ure which represents the number of customary to use a simple FIR or IIR

21 equalizer D(z) such that the product of incorporating redundancy into a sig- D(z)C(z) is a good approximation of an nal before putting it on the channel. FIR filter of small length, say L. This The system shown in the figure is is called the channel shortening step. called a filter bank precoder. By us- Now, if the integer N is chosen as N = ing standard multirate identities, the M + L – 1, we have L – 1 extra rows in system of Fig. 22 can be redrawn as the matrix R(z). It is possible to choose shown in Fig. 23(a) or equivalently as these appropriately in such a way that in Fig. 23(b). a simple set of M multipliers at the out- There are many applications which put of E(z) can equalize the channel can be described with the help of the practically completely. One special filter bank precoder configuration. If case of this idea is where the first M s(n) comes from a finite field and all rows of R(z) are chosen from the DFT the arithmetic operations in the filter- matrix and the last L – 1 rows are rep- ing are finite field operations, then we etitions of the first L – 1 rows. This can derive traditional channel coders results in a scheme called the cyclic [14] such as block coders and convo- prefix explained in detail in [22]. Fur- lutional coders as special cases of this ther interesting extensions and deeper system. These introduce redundancy in results can be found in [26]. order to make the best use of the noisy channel. A more recent application is Filter Bank Precoders the use of such redundancy in chan-

In a DMT system the symbols xk(n) nel equalization. The channel C(z) are obtained by parsing binary data can usually be approximated well by s(n) as shown earlier in Fig. 8(a). In- an FIR or stable IIR filter. The zero stead of this, imagine that the symbols forcing equalizer 1/C(z) is in general

xk(n) are obtained by blocking a sca- IIR and could even be unstable (poles lar signal s(n) which itself belongs to outside the unit circle). When we in- a PAM or QAM constellation. Then troduce redundancy as above, the use the structure of Fig. 20 can be redrawn of an IIR filter to approximate 1/C(z) as in Fig. 22. In this system, a sequence is unnecessary. In a beautiful paper [7] of symbols s(n) is converted to another Xia has shown that for almost any sequence x(n) before being fed into the channel (FIR or IIR) there exist FIR

channel. filters Ak(z) and Bk(z) such that the If we assume that two successive channel is completely equalized (i.e.,

samples of xk(n) are spaced apart by the received signal r(n) is equal to s(n) one second, then the samples of s(n) in absence of channel noise e(n)). In are separated only by 1/M seconds and fact the well known class of fraction- the channel input samples x(n) are ally spaced equalizers (FSE) [27] is a separated by the even smaller duration special case of the filter bank precoder of 1/N seconds (see the blue, green and with M =1 and uses N-fold redun- red signals in the figure). This is a way dancy. The fascinating fact about the

22 x (n) v (n) y (n) s(n) 0 0 x(n) y(n) 0 r(n) M N C(z) + N M v (n) z –1 y (n) z –1 z x 1(n) 1 z 1 M N channel N M e(n) –1 z z z z –1 R(z) E(z)

x (n) y (n) –1 z M –1 M–1 z M M –1 v N – 1 (n) z z N N deinterleaving interleaving

interleaving deinterleaving Figure 22. A modification of the redundant filter bank of Fig. 20. This is called the filter bank precoder. The only conceptual difference is in

the interpretation of the vector {x0(n), x1(n), … xM–1(n)}. The precoder allows us to equalize FIR channels with FIR filters. It also opens up the possibility of blind equalization. See text.

s(n) x(n) y(n) r(n) A 0 (z) M N C(z) + detector N M B0 (z) z –1 z A (z) M N e(n) 1 N M B1 (z) z –1 z (a) N>M z –1 z A (z) N–1 M N N M B N–1(z) Transmitter Receiving filters filters

s(n) x(n) y(n) r(n) M N F0 (z) C(z) + H0 (z) N M detector z z –1 N e(n) H (z) M F1 (z) 1 N M z z –1

(b) z z –1 M N FM– 1 (z) HM –1 (z) N M Transmitter Receiving filters filters

Figure 23. (a) The filter bank precoder redrawn to show the transmitting and receiving filters. Notice that there are N transmitting filters. (b) An equivalent drawing where the number of filters is M instead of N. Note that the delay and advance operators have been relocated as well. Both configurations have been used in the literature [7, 8].

23 Filter Banks in Digital Communications

filter bank precoder is that even if M = ample we could impose a power con- N – 1 it is almost always possible to straint and optimize these filters for have such FIR equalizers; by making maximization of signal to noise ratio N arbitrarily large we can therefore at the detector input. Some of these reduce the bandwidth expansion fac- ideas have been pursued in the context tor N / (N – 1) to almost unity and still of filter bank precoders in [8]. A fur- have FIR equalization! ther generalization of the DMT system All the above discussions assume can be obtained by using nonuniform that the channel transfer function C(z) filter banks (i.e., systems where the is known. There are many situations decimator and expander are not the where this is not true. In these cases the same in all bands). removal of intersymbol interference from the received signal falls under the References category of blind equalization. It has [1] K. Shenoi, Digital Signal Processing in been shown by Giannakis [12] that the Telecommunications. Prentice Hall, Inc., redundancy introduced by filter bank 1995. [2] T. Starr, J. M. Cioffi, and P. J. Silverman, precoders can actually be exploited to Understanding Digital Subscriber Line perform blind equalization. Further Technology. Prentice Hall, Inc., 1999. detailed results on blind as well as non- [3] G. Cherubini, E. Eleftheriou, S. Olcer, and blind equalization with filter banks can J. M. Cioffi, “Filter Bank Modulation Tech- be found in [8] and [13]. niques for Very High Speed Digital Sub- scriber Lines”, IEEE Communications Concluding Remarks Magazine, pp. 98–104, May 2000. [4] R. E. Crochiere and L. R. Rabiner, Filter banks have solved a number Multirate Digital Signal Processing. of problems in communications, but Prentice Hall, Inc., 1983. many new questions and ideas have [5] M. Vetterli, “Perfect Transmultiplexers”, Proceedings ICASSP, pp. 2567–2570, been opened up as well. In the section 1986. Biorthogonality and Perfect DMT Sys- [6] A. N. Akansu, P. Duhamel, X. Lin, and M. tems, we imposed the perfect symbol de Courville, “Orthogonal Trans- recovery condition (or biorthogonality multiplexers in Communications: A Re- IEEE Transactions on Signal Pro- condition) on the DMT filter bank and view”, cessing, April 1998. furthermore assumed a zero forcing [7] X.-G. Xia, “New Precoding for Inter- equalizer. The filters were optimized symbol Interference Cancellation Using under these two conditions. However, Nonmaximally Decimated Multirate Fil- neither of these conditions is actually ter Banks with Ideal FIR Equalizers”, IEEE Transactions on Signal Processing necessary. In fact, a priori imposition , vol. 45, no. 10, pp. 2431–2441, October of these conditions is a loss of gener- 1997. ality. It is more appropriate to optimize [8] A. Scaglione, G. B. Giannakis, and S. the transmitter and receiver filters Barbarossa, “Redundant Filter Bank jointly (the equalizer being regarded as Precoders and Equalizers Part I: IEEE part of the receiver filters). For ex- Unification and Optimal Designs”,

24 Transactions on Signal Processing, vol. [20] M. Vetterli and J. Kovaˇcevíc, Wavelets and 47, no. 7, pp. 1988–2006, July 1999. Subband Coding. Prentice Hall, Inc., 1995. [9] B. Vrcelj and P. P. Vaidyanathan, “Theory [21] P. P. Vaidyanathan and S. Akkarakaran, of MIMO Biorthogonal Partners and “A Review of the Theory and Applica- Their Application in Channel Equaliza- tions of Optimal Subband and Transform tion”, Proceedings ICC, Helsinki, Fin- Coders”, Journal of Applied and Compu- land, June 2001. tational Harmonic Analysis, to appear. [10] M. A. Tzannes, M. C. Tzannes, J. G. [22]A. Peled and A. Ruiz, “Frequency Do- Proakis, and P. N. Heller, “DMT Systems, main Data Transmission Using Reduced DWMT Systems, and Digital Filter Computational Complexity Algorithms”, Banks”, Proceedings ICC, pp. 311–315, Proceedings IEEE ICASSP, Denver, 1994. Colorado, pp. 964–967, April 1980. [11]P. P. Vaidyanathan, Y.-P. Lin, S. [23] S. Akkarakaran and P. P. Vaidyanathan, Akkarakaran, and S.-M. Phoong, “Filter Bank Optimization with Convex “Optimalilty of Principal Component Fil- Objectives, and the Optimality of Princi- ter Banks for Discrete Multitone Commu- pal Component Forms”, IEEE Transac- nication Systems”, Proceedings IEEE tions on Signal Processing, pp. 100–114, ISCAS, Geneva, May 2000. January 2001. [12] G. B. Giannakis, “Filter Banks for Blind [24] S. Akkarakaran and P. P. Vaidyanathan, Channel Identification and Equalization”, “Discrete Multitone Communication with IEEE Signal Processing Letters, vol. 4, Principal Component Filter Banks”, Pro- no. 6, pp. 184–187, June 1997. ceedings of the ICC, Helsinki, Finland, [13] A. Scaglione, G. B. Giannakis, and S. June 2001. Barbarossa, “Redundant Filter Bank [25] P. P. Vaidyanathan, Y.-P. Lin, S. Precoders and Equalizers Part II: Syn- Akkarakaran, and S.-M. Phoong, “Dis- chronization and Direct Equalization”, crete Multitone Modulation with Princi- IEEE Transactions on Signal Processing, pal Component Filter Banks”, Technical vol. 47, no. 7, pp. 2007–2022, July 1999. Report, California Institute of Technol- [14] J. G. Proakis, Digital Communications. ogy, December 2000. McGraw Hill 1995. [26] Y.-P. Lin and S.-M. Phoong, “Minimum [15] J. A. C. Bingham, “Multicarrier Modula- Redundancy ISI Free FIR Filter Bank tion for Data Transmission: An Idea Transceivers”, Proceedings of the SPIE, Whose Time Has Come”, IEEE Commu- San Diego, California, vol. 4119, pp. 745– nications Magazine, pp. 5–14, May 1990. 755, July 2000. [16] I. Kalet, “The Multitone Channel”, IEEE [27] J. R. Treichler, I. Fijalkow, and C. R. Transactions on Communications, pp. Johnson, Jr., “Fractionally Spaced Equal- 119–124, February 1989. izers: How Long Should They Be?”, [17] P. P. Vaidyanathan, Multirate Systems and IEEE Signal Processing Magazine, pp. Filter Banks. Prentice Hall, Inc., 1993. 65–81, May 1996. [18] M. J. Narasimha and A. M. Peterson, “De- [28] Y.-P. Lin and S.-M. Phoong, “Perfect Dis- sign of a 24-Channel Transmultiplexer”, crete Multitone Modulation with Optimal IEEE Transactions on Acoustics, Speech, Transceivers”, IEEE Transactions on Sig- and Signal Processing, vol. 27, Decem- nal Processing, vol. 48, pp. 1702–1712, ber 1979. June 2000. [19] J. S. Chow, J. C. Tu, and J. M. Cioffi, “A [29] M. K. Tsatsanis and G. B. Giannakis, Discrete Multitone Transreceiver System “Principal Component Filter Banks for for HDSL Applications”, IEEE Journal Optimal Multiresolution Analysis”, IEEE on Selected Areas of Communications, Transactions on Signal Processing, vol. vol. 9, no. 6, pp. 895–908, August 1991. 43, pp. 1766–1777, August 1995.

25 AA MicrophoneMicrophone ArrayArray forfor HearingHearing AidsAids

by Bernard Widrow

1531-636X/06/$10.00©2001IEEE31-636X/06 0.00© 26 This method enables the de- sign of highly-directive-hear-

Abstract—A directional acoustic re- ing instruments which are ceiving system is constructed in the comfortable, inconspicuous, form of a necklace including an array of two or more microphones mounted and convenient to use. The ar- on a housing supported on the chest of ray provides the user with a a user by a conducting loop encircling the user’s neck. Signal processing elec- dramatic improvement in tronics contained in the same housing speech perception over exist- receive and combine the microphone signals in such a manner as to provide ing hearing aid designs, par- an amplified output signal which em- ticularly in the presence of phasizes sounds of interest arriving in a direction forward of the user. The background noise, reverbera- amplified output signal drives the sup- tion, and feedback. porting conducting loop to produce a representative magnetic field. An elec- troacoustic transducer including a magnetic field pick up coil for receiv- ing the magnetic field is mounted in or Microphone Array on the user’s ear and generates an for Hearing Aids acoustic signal representative of the There is a big difference between sounds of interest. The microphone hearing speech and understanding output signals are weighted (scaled) speech. Most hearing-impaired people and combined to achieve desired spa- will be able to hear speech when given tial directivity responses. The weight- sufficient amplification from their ing coefficients are determined by an hearing aids. In many cases, however, optimization process. By bandpass fil- they will hear but will not understand. tering the weighted microphone sig- The benefits of amplification alone nals, with a set of filters covering the are limited. In a noisy place, hearing audio frequency range, and summing aids will amplify the noise as well as the filtered signals, a receiving micro- the desired speech signal. In a rever- phone array with a small aperture size berant place, hearing aids will amplify is caused to have a directivity pattern late multipath arrivals as well as the that is essentially uniform over fre- direct first-arrival signal. Furthermore, quency in two or three dimensions. feedback associated with high output This method enables the design of hearing aids distorts the frequency re- highly-directive-hearing instruments sponse of the hearing aid, which was which are comfortable, inconspicuous, carefully tuned to compensate for the and convenient to use. The array pro- individual’s hearing loss; and some- vides the user with a dramatic im- times causes oscillation. provement in speech perception over We describe a microphone array existing hearing aid designs, particu- for hearing aids that overcomes some larly in the presence of background of these limitations and has the capa- noise, reverberation, and feedback. bility of enhancing speech understand-

27 Hearing Aid with Tele-coil

hears the usual sound received by the hearing aid’s own microphone. Neck Loop The original purpose of the telecoil was to enable the hearing aid wearer to converse over the telephone. A hear- ing-aid compatible telephone receiver radiates a time varying magnetic field corresponding to the telephone signal. This is generally leakage flux from the receiver. Using the telecoil, many pa- Directional Array tients can hear over the telephone and Signal Processing Microphone much more effectively. We are able to Electronics take advantage of the telecoil, which is commonly available in the most powerful behind-the-ear hearing aid Figure 1. Chest-mounted “necklace”, directional microphone array types, to provide a wireless link be- with neck loop, and hearing aid with telecoil. tween the chest-mounted array and the hearing aid. Telecoils can be fitted to ing for hearing-impaired patients. The almost all hearing aids. microphone array is worn on the chest Use of the array enhances the as part of a necklace, in accord with the patient’s hearing in the following three diagram of Fig. 1. A processed signal ways. from the array drives current through a conducting neck loop thus creating Signal-to-noise ratio a time-variable magnetic field that is The array enhances signal-to-noise representative of the received sound. ratio. The patient aims his or her body The magnetic field provides a wireless toward the person who is speaking. means for carrying the sound signal to The array beam is 60° wide in both azi- conventional hearing aid devices lo- muth and elevation. The sound in the cated in the ears of the wearer. In or- beam is enhanced relative to the sur- der to receive the signal, the hearing rounding sound. The speech of inter- aid must be equipped with a “telecoil”, est is enhanced relative to omnidirec- a small induction coil contained within tional background noise by about 10 the hearing aid whose output can be dB, from about 200 Hz to 6 kHz. The switch selected by the wearer to serve gains of the array sidelobes vary be- in place of the hearing aid’s micro- tween 20-35 db below the gain at the phone signal. When switching the center of the main beam. hearing aid to telecoil position, the wearer hears the sound received by the Effects of reverberation array. When switching the hearing aid The array reduces the effects of re- to the microphone position, the wearer verberation. Because the array is gen-

28 A Microphone Array for Hearing Aids erally steered toward the sound of in- ments and the requirements for low terest, the direct primary path is thus battery drain. aligned with the beam. The secondary In this device, the audio spectrum paths for the most part arrive at angles from 200 Hz to 6 kHz is divided into outside the beam and are thus attenu- twelve bands, each with its own digi- ated by the array. Reducing reverbera- tal gain control. The six microphone tion enhances sound clarity since the signals are amplified and weighted and ear and the brain are somewhat con- then fed to each of the twelve band- fused by multiple arrivals. This is spe- pass filters. Different microphone-sig- cially the case with hearing-impaired nal weightings were designed for each individuals. frequency band so that the beam width was able to be held at approximately Feedback 60° over the entire frequency range of Use of the array reduces feedback interest. The microphone weights were by about 15 dB, because the chest is designed off-line by using adaptive at a much greater distance from the beamforming techniques to achieve hearing-aid loudspeaker than is the mi- the desired beam shape and to achieve crophone on the hearing aid itself. Re- a specified robustness to inherent duction of feedback makes available variations in microphone characteris- louder sound for the patient, without tics. A least square error criterion was oscillation, and allows the hearing aid used for the design. Anechoic cham- to function with a frequency response ber testing was used to verify the de- closer to the desired compensation curve. The current array design and ge- ometry are shown in Fig. 2. The device is comprised of an array of six micro- phones, four pushbuttons for control, and a plastic case designed to fit both the adult male and female torso. The plastic case was designed by computer, completely specified in software. It contains batteries and all of the signal processing electronics. Two custom ASIC chips were designed for this de- vice, one for signal processing and the other to serve as an interface between a PC computer and the signal process- ing chip when this chip is being pro- grammed. Custom chips were needed Figure 2: The current array design geometry. because of the tight space require-

29 sign. Theoretical and measured beam spectrally weighted bandpass noise. patterns turned out to be remarkably Four additional loudspeakers in the close. four corners at the ceiling were also U.S. Patent number 5,793,857 has used to carry the same noise. The room been granted to Michael A. Lehr and was not anechoic but had some sound Bernard Widrow for this technology. damping. The noise carried by the Canadian, European, and Australian eight corner loudspeakers produced a patents have been granted, and patents noise field that was approximately iso- are pending in other countries. tropic. Patient testing was performed to The test voice and the test noise evaluate the effectiveness of the micro- were stored in a PC computer. The phone array and to compare listening voice and noise data were obtained with the hearing aid alone with listen- from Dr. Sig Soli of the House Ear In- ing to the array and hearing aid in stitute in Los Angeles. We performed telecoil mode. Fig. 3 shows the floor a modified version of his HINT test arrangement of the test room. The pa- (hearing and noise test). tient was seated before a loudspeaker With the patient seated at a pre- that carried the sound of a male test scribed location marked on the floor, voice. Four loudspeakers on the floor the volume control of the hearing aid in the four corners of the room carried and the volume control of the array

Door

Noise Speakers 1 and 2

Noise Voice Speakers Speaker 5’1” 7 and 8

2’ Figure 3. Arrangement of test room.

Patient

6’6”

5’1”

Noise Noise Table Speakers Amplifier Computer Speakers 3 and 4 5 and 6

13’ 30 100 18 90 16 80 14 70 12 60 HA Only 10 50 8 40 HA+Array 6 30 4 20 Words Correct (Percent) Correct Words

SNR Improvement (dB) SNR Improvement 2 10 0 0 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 9 Subjects Subjects (a) (b)

Figure 4. Results obtained on HINT and Intelligibility tests. (a) Magnitude of improvement in sentence speech recognition threshold in noise (HINT) with the microphone array, in comparison with the hearing aid alone. (b) Percent speech intelligibility in the presence of noise. For each test subject, the purple bar is the result obtained using the hearing aid and the blue bar is the result using the array. were set so that the measured volume when using the array is plotted in delivered to the patient’s ear would be Fig. 4(a) for nine test patients. This im- the same when listening to the test provement averages more than 10 dB, voice through the hearing aid and which is consistent with anechoic through the array. The volume level of chamber measurements and theoretical the test voice was set to be comfortable calculations. for the patient, in the absence of noise. Other testing was done with the Word phrases were spoken to the noise volume fixed and the volume patient by the test voice, with some level of the test voice fixed. Individual noise applied. The patient was asked words randomly selected were pre- to repeat the words. If any word in the sented by the test voice. The responses phrase was repeated incorrectly, the of the patients were observed when response was considered to be incor- using the hearing aid, and when using rect. The noise level was reduced the array. The results are shown for the by 2 dB, and another randomly chosen same nine patients, in Fig. 4(b). Patient phrase was read. If the response was #1 had a 25% correct response with the incorrect again, the noise was lowered hearing aid, and a 95% correct re- by another 2 dB and so forth. When a sponse with the array. Patient #2 had correct response was obtained, the a 15% correct response with the hear- noise level for the next phase was ing aid, and an 80% correct response raised by 2 dB. If another correct re- with the array. And so forth. These sponse was obtained, the noise level improvements are rather dramatic. was raised by another 2 dB and so One young woman in Palo Alto, forth. The noise level went up and California, has been wearing one of down, and the average noise level was these devices on a daily basis over the observed over ten or twenty phrases. past five years. As the design evolved, The average noise level when us- she always had the latest for testing. ing the hearing aid was compared to She is totally deaf in one ear and is 95– that when using the array. The im- 105 db below normal in her “good” provement in signal-to-noise ratio ear. Using her hearing aid and with

31 A Microphone Array for Hearing Aids

One young woman in Palo Alto, California, has been wearing one of these devices on a daily basis over the past five years. As the design evolved, she always had the latest for testing. She is totally deaf in one ear and is 95-105 db below normal in her “good” ear. Using her hearing aid and with good lip reading, she can correctly recognize zero to two words in a typical long sentence. With her hearing aid and an array, she gets essentially every word. She can do very well even with her eyes closed. Her hearing loss is in the profound range. Hearing losses Bernard Widrow Stanford University are generally characterized as mild, moderate, se- Dept. of Electrical Eng. Packard Bldg., Rm. 273 vere and profound. The array will find its best ap- Stanford, CA 94305 phone: (650) 723–4949 plication with the difficult cases, the severe and pro- fax: (650) 723–8473 [email protected] found ones.

good lip reading, she can correctly rec- profound ones. ognize zero to two words in a typical The microphone array devices are long sentence. With her hearing aid now being manufactured and marketed and an array, she gets essentially ev- by Starkey Laboratories, 6700 Wash- ery word. She can do very well even ington Ave, Eden Prairie, MN, 55344, with her eyes closed. Her hearing loss U.S.A. The trade name for the device is in the profound range. Hearing is Radiant Beam Array (RBA). It is the losses are generally characterized as most powerful hearing device on the mild, moderate, severe and profound. market. It remains to be seen how well The array will find its best application it will be accepted by the hearing-im- with the difficult cases, the severe and paired community.

Bernard Widrow received the S.B., M.S. and Sc.D. degrees in electrical engineering from the Massachusetts Institute of Technology in 1951,1953, and 1956, respectively. He joined the MIT faculty and taught there from 1956–1959. In 1959, he joined the faculty of Stanford University, where he is currently Professor of Electrical Engineering. Dr. Widrow is a Life Fellow of the IEEE and a fellow of AAAS. He received the IEEE Centennial Medal in 1984, the IEEE Alexander Graham Bell medal in 1986, the IEEE Neural Networks Pioneer Medal in 1991, the IEEE Signal Processing Society Medal in 1998, the IEEE Millennium Medal in 2000, and the Benjamin Franklin Medal of the Franklin Institute in 2001. He was inducted into the National Academy of Engineering in 1995, and into the Silicon Valley Engineering Council Hall of fame in 1999. Dr. Widrow is a past president and currently a member of the Governing Board of the international Neural Network Society. He is asso- ciate editor of several journals, and is the author of about 100 technical papers and 15 patents. He is co-author of Adaptive Signal Processing and Adaptive Inverse Control, both Prentice Hall books. A new book, Quantization Noise, is in preparation. 32 TransactionsTransactions

speech recognition into many environments. Two demonstra- Prize Paper Awards tions are described: one is PC-based and is near-real-time, the s noted in the ISCAS’01 report in this issue (see page 39), second is digital signal processing based and is real-time. Asix best paper awards were announced at the International Sample results are described to quantify the performance of Symposium on Circuits and Systems in Sydney, Australia. Listed the overall systems. below are the awards and the papers which received them. IEEE Transactions on Circuits and Systems, Part II: Analog and Digital Signal Processing, July 1999, pp. GUILLEMIN-CAUER AWARD 915–922. “Filtering through Combination OUTSTANDING YOUNG AUTHOR AWARD of Positive Filters Criteria” Luca Benvenuti, Lorenzo Farina, “Exploiting Symmetry When Verifying and Brian D. O. Anderson Transistor-Level Circuits by Symbolic Trajectory Evaluation” Abstract—The linear filters characterized by a state-variable realization given by Manish Pandey (with Randal E. Bryant) matrices with nonnegative entries (called Abstract—We describe the use of symmetry for verifi- positive filters) are heavily restricted in cation of transistor-level circuits by symbolic trajectory evalu- their achievable performance. Neverthe- ation (STE). We present a new formulation of STE which al- less, such filters are the only choice when lows a succinct description of symmetry properties in circuits. dealing with the charged coupled device Symmetries in circuits are classified as structural symmetries, MOS technology of charge routing net- arising from similarities in circuit structure, data symmetries, works (CRN’S), since nonnegativity is a arising from similarities in the handling of data values, and consequence of the underlying physical mixed structural-data symmetries.We use graph isomorphism mechanism. In order to exploit the advantages offered by this testing and symbolic simulation to verify the symmetries in technology, the authors try to overcome the above-mentioned the original circuit. Using conservative approximations, we limitation by realizing an arbitrary transfer function as a dif- partition a circuit to expose the symmetries in its components, ference of two positive filters. and construct reduced system models which can be verified IEEE Transactions on Circuits and Systems, Part I: Fun- efficiently. Introducing X-drivers into switch-level circuits damental Theory and Applications, December 1999, pp. simplifies the task of creating conservative approximations of 1431–1440. switch-level circuits. Our empirical results show that exploit- ing symmetry with conservative approximations can allow one DARLINGTON AWARD to verify systems several orders of magnitude larger than oth- erwise possible. We present results of verifying Static Random “Voice Extraction by Access Memory circuits with up to 1.5 Million transistors. On-Line Signal IEEE Transactions on Computer-Aided Design of In- Separation and Recovery” tegrated Circuits and Systems, July 1999, pp. 918–935. Gail Erten and Fathi M. Salam CSVT TRANSACTIONS Abstract—The paper presents a for- mulation and an implementation of a sys- BEST PAPER AWARD tem for voice output extraction (VOX) in real-time and near-real-time realistic real- “On End-to-End Architecture for Transporting world applications. A key component includes voice-signal MPEG-4 Video over the Internet” separation and recovery from a mixture in practical environ- Dapeng Wu, Yiwei Thomas Hou, Wenwu Zhu, ments. The signal separation and extraction component in- cludes several algorithmic modules with a variety of sophisti- Hung-Ju Lee, Tihao Chiang, Ya-Qin Zhang, cation levels, which include dynamic processing neural net- and H. Jonathan Chao works in tandem with (dynamic) adaptive methods. These Abstract—With the success of the Internet and flexibility adaptive methods make use of optimization theory subject to of MPEG-4, transporting MPEG-4 video over the Internet is the dynamic network constraints to enable practical algorithms. expected to be an important component of many multimedia The underlying technology platforms used in the compiled applications in the near future. Video applications typically VOX software can significantly facilitate the embedding of have delay and loss requirements, which cannot be adequately

33 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS

A PUBLICATION OF THE IEEE CIRCUITS AND SYSTEMS SOCIETY WWW.IEEE-CAS.ORG

IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING

A PUBLICATION OF THE IEEE CIRCUITS AND SYSTEMS SOCIETY WWW.IEEE-CAS.ORG

IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY A PUBLICATION OF THE IEEE CIRCUITS AND SYSTEMS SOCIETY WWW.IEEE-CAS.ORG

IEEE TRANSACTIONS ON COMPUTER-AIDED DESIGN OF INTEGRATED CIRCUITS AND SYSTEMS

A PUBLICATION OF THE IEEE CIRCUITS AND SYSTEMS SOCIETY WWW.IEEE-CAS.ORG

IEEE TRANSACTIONS ON VERY LARGE SCALE INTEGRATION (VLSI) SYSTEMS A JOINT PUBLICATION OF THE IEEE CIRCUITS AND SYSTEMS SOCIETY THE IEEE COMPUTER SOCIETY o1 THE IEEE SOLID-STATE CIRCUITS SOCIETY SSCS

34 TransactionsTransactions supported by the current Internet. Thus, it is a challenging prob- lem to design an efficient MPEG-4 video delivery system that VLSI TRANSACTIONS can maximize the perceptual quality while achieving high re- BEST PAPER AWARD source utilization. This paper addresses this problem by pre- senting an end-to-end architecture for transporting MPEG-4 “Toward Achieving Energy Efficiency in video over the Internet. We present a framework for transport- Presence of Deep Submicron Noise” ing MPEG-4 video, which includes source rate adaptation, Rajamohana Hegde and Naresh R. Shanbhag packetization, feedback control, and error control. The main contributions of this paper are: 1) a feedback control algorithm Abstract—Presented in this paper are 1) information- based on Real Time Protocol (RTP) and Real Time Control Pro- theoretic lower bounds on energy consumption of noisy digi- tocol (RTCP); 2) an adaptive source-encoding algorithm for tal gates and 2) the concept of noise tolerance via coding for MPEG-4 video which is able to adjust the output rate of MPEG- achieving energy efficiency in the presence of noise. In par- 4 video to the desired rate; and 3) an efficient and robust ticular, lower bounds on a) circuit speed fc and supply voltage packetization algorithm for MPEG video bit-streams at the sync Vdd; b) transition activity t in presence of noise; c) dynamic layer for Internet transport. Simulation results show that our end- energy dissipation; and d) total (dynamic and static) energy to-end transport architecture achieves good perceptual picture dissipation are derived. A surprising result is that in a scenario quality for MPEG-4 video under low bit-rate and varying net- where dynamic component of power dissipation dominates, work conditions and efficiently utilizes network resources. the supply voltage for minimum energy operation (Vdd, opt) is greater than the minimum supply voltage (Vdd, min) for reliable IEEE Transactions on Circuits and Systems for Video operation. We then propose noise tolerance via coding to ap- Technology, September 2000, pp. 923–941. proach the lower bounds on energy dissipation. We show that the lower bounds on energy for an off-chip I/O signaling ex- CAD TRANSACTIONS ample are a factor of 24x below present day systems. A very BEST PAPER AWARD simple Hamming code can reduce the energy consumption by a factor of 3x, while Reed-Muller (RM) codes give a 4x re- “SPFD: A New Method to duction in energy dissipation. Express Functional Flexibility” IEEE Transactions on Very Large Scale Integration Shigeru Yamashita, Hiroshi Sawada, (VLSI) Systems, August 2000, pp. 379–391. and Akira Nagoya Abstract—In this paper, we propose a unique way to ex- CALL FOR PAPERS press functional flexibility by using sets of pairs of functions called “Sets of Pairs of Functions to be Distinguished” (SPFDs) IEEE TVLSI Special SLIP Issue rather than traditional incompletely specified functions. This method was very naturally derived from a unique concept for he IEEE Transactions on VLSI Systems will publish a distinguishing two logic functions, which we explain in detail Tspecial issue on System Level Interconnect Prediction in in this paper. the summer of 2002. This special issue focuses on system-level The flexibility represented by an SPFD assumes that the in- predictions and a priori and on-line (i.e., pre-layout and dur- ternal logic of a node in a circuit can be freely changed. SPFDs ing layout) estimations of interconnect related design param- make good use of this assumption, and they can express larger eters and characteristics (delay, power, yield, wiring layers, ...). flexibility than incompletely specified functions in some cases. It intends to build on the research work presented in the Although the main subject of this paper is to explain the first special issue on System-Level Interconnect Prediction concept of SPFDs, we also present an efficient method for cal- (T-VLSI December 2000) and also on work that will be pre- culating the functional flexibilities by SPFDs because the con- sented in the second special issue (targeted for the fall of 2001). cept becomes useful only if there is an efficient calculation Authors are encouraged to submit high-quality research con- method for it. Moreover, we present a method to use SPFDs tributions that will not require major revisions. Extensions of for circuit transformation along with a proof of the correct- papers presented at the IEEE/ACM International Workshop on ness of the method. System-Level Interconnect Prediction (SLIP 2001 http:// We further make a comparison between SPFDs and com- www.SLIPonline.org) are welcomed, but work not presented patible sets of permissible functions (CSPFs), which express at the workshop is also encouraged. functional flexibility by incompletely specified functions. More information can be downloaded from As an application of SPFDs, we show a method to opti- http://www.ece.udel.edu/~christie/Research/Slip mize LUT (look-up table) networks and experimental results. The submission deadline is August 1st, 2001. IEEE Transactions on Computer-Aided Design of Inte- —Phillip Christie, University of Delaware grated Circuits and Systems, August 2000, pp. 840–849. Guest editor TVLSI SLIP issue 35 TechnicalTechnical CommitteesCommittees

with the supply voltage being near 0.6 V. The transistor cutoff New Adventure by frequency will be beyond 500 GHz. While the silicon semiconductor technology progresses IEEE CAS Society toward the nano-scale regime, new nano-structures or nano- devices in the scale of 0.1–10 nm are under aggressive devel- Introduction to Nanoelectronics opment with the possibility to supplement the capability of and Gigascale Systems silicon technologies. Many innovative nano-structures/nano- Technical Committee devices such as resonant tunneling diodes, single-electron tun- neling (SET) transistors, quantum dots, quantum wires, tera- by Chung-Yu Wu, Founding Chair herz devices, and DNA bio-devices, have been proposed and [email protected] investigated. These nano-devices and nano-structures make delicate use of quantum phenomena or interaction to perform ith a vast amount of research and development from creative physical functions. Whigh-tech industry and strong demands from personal Although technology and physics of nano-devices and computers, cellular phones, and Internet, the silicon semicon- nano-structures in silicon and other materials are increasingly ductor technology continues its marvelous progress into the understood, it is still a great challenge to integrate millions and nanometer regime in the 21st century. The great advances of billions of these devices into gigascale micro-systems. For the semiconductor mass-production technology closely follow 18-nm technology node with BJT (bipolar junction transistor)– Moor’s Law in scaling down the effective transistor channel like gate current characteristics and 0.1-V threshold voltage, length. It is expected that in the 130–100 nm CMOS genera- over a 1 Tera-Bit DRAM chip or logic processor chip could tion for mass production in around years 2003–2006, the gate- be built. At such an unprecedented integration level, the num- oxide thickness will be 2 nm with the gate leakage current of ber of transistors on a single chip will be ten times the number 0.05 nA/µm2, which is about 1/20 of the off-state leakage cur- of neurons in a human brain. There are many research challenges rent through the transistor channel. The corresponding thresh- and opportunities from the circuit and system perspective to har- old voltage of the transistor is about 0.3 V with the power sup- vest the gigascale integration. Bio-inspired or neuromorphic ply voltage close to 1 V. The cutoff frequency is beyond 100 networks might provide some stimulating information-process- GHz [1]. A significant barrier will be encountered at around ing architectures to alleviate severe problems of internal data the 18-nm technology node in the year 2020 as shown in Fig.1. transfer and the I/O function at the microchip peripheral. In that technology generation, the gate-oxide thickness is be- Similarly, there exist great challenges in using nano-struc- low 1 nm with the gate leakage current of 0.05µA/µm2, which tures/nano-devices to form nanoelectronic integration systems. is about 50 times off the off-state leakage current through the In nanoelectronics, quantum phenomena could be used for new transistor channel. The threshold voltage will be below 0.1 V signal processing operations. Electrostatic coupling or Cou- 10 15

10000 Number of transistors on a single chip DRAM 14 Figure 1. Semiconductor roadmap prediction toward 2020 10 1.4 time/year 10 13 1000 800 1TB 12 500 10 350 64GB (2023) 250 11 180 (2015) 10 130 100 Number of Human 10 10 100 70 50 Neurons 35 10 9 25 18 10 8

Device dimension (nm) 10 Very Difficult in Technology 10 7 10 6 5 1 10 1989 1990 1993 1996 1999 2002 2005 2008 2011 2014 2017 2020 2023 2026 Year 36 TechnicalTechnical CommitteesCommittees

Chung-Yu Wu Bing Sheu Leon Chua Josef Nossek

Tamas Roska Rui de Figueiredo Eby Friedman Mona Zaghloul

George Moschytz Fathi Salam David Chung-Laung Liu Mohammed Ismail

Andreas Andreou Robert Newcomb Keshab Parhi Cliff Lau 37 TechnicalTechnical CommitteesCommittees

formed in the IEEE Circuits and Systems Society in early 2001. The Technical Committee on Nanoelectronics and Gigascale Sys- tems has been actively participating in the publications, confer- ences, and meetings of this fast-growing profession. Active and interested researchers and engineers are all welcome. Many enthusiastic researchers and engineers are actively involved. They include the following: Prof. Chung-Yu Wu (Technical Committee Founding Chair), National Chiao Tung University; Prof. Leon Chua (1977 President of CASS), Uni- versity of California at Berkeley; Prof. Rui de Figueiredo (1998 President of CASS), University of California at Irvine; Dr. Bing Sheu (Past President of CASS), Nassda Corporation; Arpad Csurgay Wolfgang Porod Prof. Mona Zaghloul (VP-Technical Activities of CASS), George Washington University; Prof. George Moschytz (1999 President of CASS), Swiss Federal Institute of Technology; Prof. Josef Nossek (President-Elect of CASS), Munich Uni- versity of Technology; Prof./President David Chung-Laung Liu, National Tsing Hua University, (Contributing Researcher) Prof./President Chun-Yen Chang, National Chiao Tung Uni- versity; Prof. Robert Newcomb, University of Maryland; Dr. Cliff Lau, Office of Naval Research; Prof. Tamas Roska, Hun- garian Academy of Sciences; Prof. Arpad Csurgay, Technical University of Budapest; Prof. Eby Friedman, University of Rochester; Prof. Fathi Salam, Michigan State University; Prof. Mohammed Ismail, Ohio State University; Prof. Andreas Andreou, John Hopkins University; Prof. Keshab Parhi, Uni- versity of Minnesota; Prof. Wolfgang Porod, University of Liang-Gee Chen Sethuraman Panchanathan Notre Dame; Prof. Liang-Gee Chen, National Univer- sity; Prof. Ramalingam Sridhar, State University of New York at Buffalo; Prof. Sethuraman Panchanathan, Arizona State University; Dr. Jeong-Taek Kong, Samsung Electronics Co.; Dr. Pedro Julian, University of California at Berkeley. Many more will join the Technical Committee later. References [1] Y. Taur, “The Incredible Shrinking Transistor”, IEEE Spectrum, vol. 36, pp. 25–29, July 1999. [2] W. Porod, “Towards Nanoelectronics: Possible CNN Implemen- tations Using Nanoelectronics Devices”, Proceedings of IEEE International Workshop on Cellular Neural Networks and Their Applications, pp. 20–25, London, England, April 1998. Jeong-Taek Kong Pedro Julian

lomb interactions could be applied for interconnection or I/O access. The promising architecture for nanoelectronic proces- sors is quantum device or molecule arrays, such as quantum- dot cellular nonlinear networks, quantum-dot cellular automata [2], and DNA array chips. The development of nanoelectronics could lead to the applications of quantum computers, molecular computers, or DNA computers in industrial/business work and in daily life. This will bring us into a new Era in human history. With the exciting opportunities described above, it is clear that researches on nanoelectronics and gigascale systems are tightly coupled together. To promote and excel in the areas of nanoelectronics and gigascale systems for the circuits and sys- tems research/design community, a new technical committee was Chun-Yen Chang Ramalingam Sridhar

38 MeetingsMeetings

ISCAS 2001 in Sydney – president, Ibrahim Hajj as VP-Administration, M.N.S. Swamy as VP-Publications, Mona Zaghloul as VP-Technical Activi- A Great Success ties, Geert de Veirman as VP-Conferences, plus regional VPs: Wasfy Mikhael, Anthony Davies, Juan Cousseau, and Nobuo he 2001 IEEE International Symposium on Circuits and Fujii. Great service contributions by Society administrator Bar- TSystems (ISCAS) was very successfully held at the Con- bara Wehner and her staff Tom Wehner were publicly recognized vention Center near Darling Harbor of Sydney, Australia dur- with a plaque presented by Graham Hellestrand. ing May 5–9, 2001. More than one thousand and two hundred The Symposium banquet on Tuesday evening was held researchers, engineers, and students from the scientific/aca- in the new ballroom at the Convention Center. Truly consid- demic and industrial community gathered together and pre- erate arrangements, warm friendship, and best services from sented their newest and most exciting technical results/find- the Symposium team deeply touched the hearts of all the par- ings and exchanged ideas with their peers. As the historic first ticipants. The awards ceremony was held before the banquet. IEEE ISCAS conference to be held in the Southern Hemi- Many awards were presented, including Society Awards, prize- sphere, ISCAS-2001 shines like a bright star in the history of paper awards, and IEEE Fellow awards. For the Society IEEE Circuits and Systems Society. Awards, only Prof. Robert Newcomb attended the ceremony All the participants greatly enjoyed the keynote presen- to receive the Society Education Award. The others will have tations that highlighted industrial contributions to the modern to be presented in other IEEE meetings. Figure 1 shows the circuits and systems field: “The Impact of SOC Methodolo- photo of Prof. Robert Newcomb, with Prof. Fathi Salam who gies on Circuits and Systems” prepared by Fred Shlapak of is recipient of the Darlington Award for a T—CAS II paper, Motorola Inc., and “Cochlear Implants: Technology and Ap- Dr. Bing Sheu who serves as the 2001 Society past president, plications” by John Parker of Cochlear Ltd. Barrie Gilbert of and Prof. Chung-Yu Wu who serves as Society technical edi- Analog Devices Inc. also prepared the material on “Analog tor to the Circuits & Devices Magazine. Design at the Millennial Crossroads”, but could not attend due The IEEE Circuits and Systems Society continues to to illness. One popular panel session was held on May 9 in maintain a strong technical and financial position. The ISCAS the afternoon to address “Successful Entrepreneurs: Generat- attendees all plan on coming to Phoenix, Arizona to enjoy the ing Wealth from Circuits and Systems R&D (Not Necessarily warm weather and outdoor activities in May 2002. in Silicon Valley)”. The panelists from successful start-up com- panies included Neil Weste of Radiata Inc., Antoni Cantoni —Chung-Yu Wu, CAS Technical Editor of Atmoshpere Inc., and Chris Roberts of ResMed Corp. Circuits & Devices Magazine A unique feature in this year’s Technical Program was the —Michael Sain, Editor-in-Chief Interactive Sessions that included brief oral presentations by IEEE CAS Magazine the authors for the poster papers. It helped attract large audiences to the poster papers room. Many authors indicated strong in- terest in participating in similar arrange- ments in future ISCAS conferences. The tutorial sessions on May 6 were heavily sponsored by a special grant from the IEEE Circuits and Systems Society to al- low more audience to attend the various talks. Conference general co-chairs, David Skellern of Radiata Inc. and Macquarie University, with Graham Hellestrand of VaST Systems Technology Corporation and University of New South Wales, the technical program chair Chris Toumazou of Imperial College, and their organizing committee deserve recognition for the wonderful success of the confer- ence. This year’s Society executive com- mittee strongly supports the ISCAS con- Figure 1. Photo taken after awards ceremony. (From left: Dr. Bing Sheu, ference, and the members include Hari Society past president; Prof. Robert Newcomb, recipient of Society Education Reddy as Society president, Josef Nossek Award; Prof. Fathi Salam, recipient of Darlington Award; and Chung-Yu Wu, as president-elect, Bing Sheu as past Society technical editor to the Circuits & Devices Magazine.)

39 SocietySociety

www.ieee-cas.org/html/student-promo.html, fill it out com- Swamy Assumes pletely, provide necessary signatures and fax or mail it to: Publication Post IEEE Circuits & Systems Society P.O. Box 265 M.N.S. Swamy is serving as Vice 15 West Marne Ave. Beverly Shores, IN 46301 President—Publications for the USA Circuits and Systems Society dur- Fax: +1 219 871 0211 ing the years 2001–2002. Applicants should make sure that ALL information is pro- vided completely and legibly. INCOMPLETE applications will M.N.S. Swamy is presently a research professor and the be rejected. Students with additional questions may send an director of the Center for Signal Processing and Communica- email message to [email protected]. tions in the Department of Electrical and Computer Engineer- Terms & Limitations ing at Concordia University, Montreal, Canada, where he • This offer is valid until 1 August 2001. Any applica- served as chair of the Department of Electrical Engineering tions received after this date will be rejected. from 1970 to 1977, and dean of Engineering and Computer • This free student membership is being offered on a first- Science from 1977 to 1993. He has published extensively in come, first-served basis until the allocated number of the areas of number theory, circuits, systems and signal pro- memberships has been exhausted. cessing, and holds four patents. He is the co-author of two book • The IEEE CAS Society reserves the right to withdraw chapters and three books, one of which has been translated into this offer at any time. Chinese and Russian. He is the editor-in-chief of the IEEE Transactions on Circuits and Systems—I and associate editor of the Fibonacci Quarterly. Dr. Swamy is a fellow of a number of professional soci- Trajkovic Begins Three-Year eties including the IEEE, IEE (UK), and the Engineering In- Term on Society Board stitute of Canada. He has served the IEEE in various capaci- ties such as vice president of the IEEE CAS Society in 1976, Ljiljana Trajkovic is now serving the first program chair for the 1973 IEEE CAS Symposium, general chair for the 1984 IEEE CAS Symposium, vice-chair for the year of a three-year term as an elected mem- 1999 IEEE CAS Symposium, and a member of the Board of ber of the CAS Society Board of Governors. Governors. He is a co-recipient of the IEEE CAS 1986 Guillemin-Cauer Best Paper Award. Recently, he received the She presently serves as a member of the Con- IEEE-CAS Society Golden Jubilee Medal, as well as the year ference and Regional Activities Divisions. 2000 IEEE-CAS Society Education Award in recognition of Ljiljana Trajkovic received the Dipl. Ing. degree in her his contributions. native Yugoslavia in 1974, M.Sc. degrees in electrical engi- neering and computer engineering from Syracuse University in 1979 and 1981, respectively, and a Ph.D. degree in electri- FREE Student Membership cal engineering from UCLA in 1986. Dr. Trajkovic is currently Promotion for 2001 an Associate Professor in the School of Engineering Science at Simon Fraser University, Burnaby, British Columbia, ntil August 1, 2001, the CAS Society is offering FREE Canada. From 1995 to 1997, she was an NSF Visiting Pro- Umembership in IEEE and the CAS Society for students. fessor in the Electrical Engineer- The only stipulations are that the student must be an active ing and Computer Sciences De- full time graduate or senior undergraduate student in an elec- partment at UC Berkeley. She trical/electronic engineering or computer engineering/sciences was a Research Scientist at Bell discipline, must be a first time applicant to the IEEE and CAS Communications Research, Society, and must have his or her application signed by a fac- Morristown, NJ from 1990 to ulty member who is also an IEEE member. 1997, and a Member of the This offer is NOT valid for renewals, previous student Technical Staff at AT&T Bell members whose membership has lapsed, recent graduates, or Laboratories, Murray Hill, NJ full IEEE/CAS members (i.e. employed in an engineering pro- from 1988 to 1990. fession). Her research interests in- Students who meet the above requirements should down- clude high-performance com- load the special Student Membership Application at http:// munication networks, com-

40 SocietySociety puter-aided circuit analysis and design, and theory of nonlin- CHAPTER-OF-THE-YEAR AWARD ear circuits and dynamical systems. Russian CAS Chapter in Moscow Dr. Trajkovic was recently elected to the IEEE Circuits Chair: Vagan V. Shakhguildian and Systems Board of Governors. She is also associate editor Vice Chair: Alexander S. Dmitriev of the IEEE Transactions on Circuits and Systems (Part II). Justification for selection: For outstanding contributions She was chair of the IEEE Technical Committee on Nonlin- to promoting circuits and systems related technical activities ear Circuits and Systems (1998), and associate editor of the in the Moscow area as well as in other regions of the country. IEEE Transactions on Circuits and Systems (Part I) from 1993–1995. She also serves on various technical committees of IEEE conferences. She is the vice chair of the upcoming THE ADVENTURES OF … ISCAS 2004, to be held in Vancouver, British Columbia. …THE 'UMBLE OHM Editor’s Note: Beginning with this issue, Circuits and Systems …Shlomo Karni Magazine is profiling recently elected Society officers and members of the Board of Governors. Other electees will appear in future is- sues. Circuits and Systems Society Names 2001 Awardees lso noted in the ISCAS’01 report in this issue (see page 39), Asix Society awards were announced at the International Symposium on Circuits and Systems in Sydney, Australia. Listed below are the awards and the persons who received them. M.E. VAN VALKENBURG AWARD Alfred Fettweis ohminy: hulled corn Citation: For significant, creative and innovative research contributions to circuits and systems theory and technology, and particularly to the discipline of digital signal processing. Society Announces Divisions TECHNICAL ACHIEVEMENT AWARD resident Hari C. Reddy has announced the organization Ruey-Wen Liu Pof the Circuits and Systems Society into its four constitu- Citation: For his leadership and fundamental contributions ent divisions: Conference, Publications, Technical Activities, in nonlinear circuit theory and blind signal processing; and for and Regional Activities. Membership and contact information his unceasing leadership in enhancing and sustaining the tech- follows. nical vibrancy of the Circuits and Systems Society. SOCIETY EDUCATION AWARD Conference Division Robert W. Newcomb Chair — Geert De Veirman Citation: For his outstanding contributions to engineer- Tel: (714) 573–6829; Fax: (714) 573–6916 ing education, the publication of pioneering textbooks in cir- Ian Galton cuit theory, VLSI, and control, and his classic work on Linear Tel: 858–822–1332; Fax: 858–822–3425 Multiport Synthesis, and for being a creative and inspiring teacher and researcher in the areas of circuits and systems, mi- Ljiljana Trajkovic croelectronics, and neural systems for close to four decades. Tel: (604) 291–3998; Fax: (604) 291–4951 Martin Hasler MERITORIOUS SERVICE AWARD Tel: +41 21 693 2622; Fax: +41 21 693 6700 Philip V. Lopresti Peter Pirsch Citation: For providing a quarter century of outstanding Tel: +49–511–762–19640; Fax: +49–511–762–19601 and faithful leadership and service to the CAS Society. Georges Gielen INDUSTRIAL PIONEER AWARD Tel: ++32–16–321047; Fax: ++32–16–321975 Aart De Geus Ex-Officio — Hari C. Reddy Citation: For pioneering logic synthesis. Tel: (562) 985–5106; Fax: (562) 985–8887

41 SocietySociety RecognitionsRecognitions Publications Division Chair — M.N.S. Swamy IEEE CAS Fellow Profiles 2001 Tel: (514) 848–3091; Fax: (514) 848–2802 Massoud Pedram Tel:213 740 4458; Fax:213 740 9803 Hon-Sum Philip Wong Robert J. Marks For contributions to solid-state image sensors and Tel: (206) 543–6990; Fax: (206) 543–3842 nanoscale CMOS devices. Eby Friedman H.-S. Philip Wong received the B.Sc. (Hons.) degree from Tel: (716) 275–1022; Fax: (716) 275–2073 the University of Hong Kong in 1982, the M.S. degree from Rajesh K. Gupta the State University of New Tel: 949 824 8052; Fax: 949 824 8019 York, Stony Brook, in 1983, and the Ph.D. degree in electri- Hiroto Yasuura cal engineering from Lehigh Tel:+81–92–583–7620; Fax:+81–92–583–1338 University, Bethlehem, Penn- Chris Toumazou sylvania, in 1988. He joined Tel: +44–171–594–6255; Fax: +44–171–581–4419 the IBM Thomas J. Watson Ex-Officio — Hari C. Reddy Research Center, Yorktown Tel: (562) 985–5106; Fax: (562) 985–8887 Heights, in 1988, where he is Technical Activities Division now senior manager of explor- atory devices and integration Chair — Mona Zaghloul technology. Tel: (202) 994–3772; Fax: (202) 994–0227 Since 1993, he has been Michael Green working on the device physics, Tel: 949 585 6854; Fax: 949 450 1079 fabrication, and applications of nanoscale CMOS devices. His Jason Cong recent work has been on the physics and fabrication technol- Tel: 310–206–2775; Fax: 310–825–2773 ogy of double-gate and back-gate MOSFETs for CMOS tech- Ellen Yoffa nologies toward the 25 nm channel length regime. Tel: (914) 945–3270; Fax: (914) 945–4349 In the applications arena, his work has been on solid-state Chris Toumazou imaging. His recent work has been imaging devices using Tel: +44–171–594–6255; Fax: +44–171–581–4419 CMOS technologies and impact of device scaling on CMOS Andreas G. Andreou imaging systems. His interest in solid-state imaging began when he joined IBM in 1988. From 1988 to 1992, he was a Ex-Officio — Hari C. Reddy member of a team that worked on the design, fabrication, and Tel: (562) 985–5106; Fax: (562) 985–8887 characterization of a high resolution, high color-fidelity CCD Regional Activities Division image scanner for art work archiving. These scanners are now Chair—Josef Nossek in use at several premier museums around the world. Tel: +49–89–2892–8501; Fax: +49–89–2892–8504 VP R1–7 — Wasfy Mikhael Tel: (407) 823–3210; Fax: (407) 823–5835 Masao Hotta VP R8 — Anthony C. Davies For contributions to the development of low-power video-fre- Tel: +44 20 7848 2441; Fax: +44 20 7836 4781 quency analog-to-digital converters for mixed-signal system VP R9 — Juan Cousseau large scale integrated circuits. Tel: +549 155 8679; Fax: +549 155 3122 Masao Hotta received the VP R10 — Nobuo Fujii M.S. and Ph.D. degrees in elec- Tel: +81 3 5734 2561; Fax: +81 3 5734 2909 tronics from Hokkaido Univer- sity, Sapporo, Japan, in 1973 Ellen Yoffa and 1976, respectively. Tel: (914) 945–3270; Fax: (914) 945–4349 Since 1976 he has been Ljiljana Trajkovic with the Central Research Tel: (604) 291–3998; Fax: (604) 291–4951 Laboratory, Hitachi Ltd., To- Georges Gielen kyo, Japan. He initially en- Tel: ++32–16–321047; Fax: ++32–16–321975 gaged in research and develop- Ex-Officio — Hari C. Reddy ment of high precision mono- Tel: (562) 985–5106; Fax: (562) 985–8887 lithic D/A converters with reso-

42 RecognitionsRecognitions lution of 14 bits. From 1981 he started research and develop- Middle East Technical Univer- ment of high-speed A/D converters for video use, and ultra- sity (METU) in Ankara and af- high-speed D/A converters for high definition display appli- ter 1974 at ITU. He was also cations. Since 1986, he was senior researcher conducting re- responsible for conducting re- search on analog circuits and high precision DACs and ADCs. search at TUBITAK Marmara From 1995 to 1999, Dr. Hotta was manager of the Advanced Scientific and Industrial Re- Device Development Department, Semiconductor Technology search Institute (MRE) where Development Center, Semiconductor & Integrated Circuits Di- he was head of the electronics vision, Hitachi Ltd., conducting development on microproces- research department and then sors, memories, RF devices and DA/CAD systems. He is pres- became the director of the in- ently senior chief engineer and senior manager of the Advanced stitute. After serving seventeen Analog Technology Center, Semiconductor & Integrated Cir- years at MRE he asked for cuits, Hitachi Ltd. He is conducting development on RF power early retirement from ITU. He amplifier modules, RF transceiver LSIs and mixed-signal LSIs. was at Bilkent University, An- Dr. Hotta served on the technical program committees of kara (1989–1992) , and Eastern Mediterranean University at CICC, BCTM and ASIC/SOC Conferences. North Cyprus (1992-1997). He is now a professor at ISIK Uni- versity in Istanbul. In 1982 he received the Science Award from Chai Wah Wu the Turkish Scientific and Technical Research Council For contributions to synchronization of chaotic systems and (TUBITAK) for his contributions in network theory. He has its applications. co-authored the book Analysis of Discrete Physical Systems (McGraw-Hill Book Co., 1967) and also ten books on circuits Chai Wah Wu received the and systems, seven of them in Turkish. B.S. in computer engineering and the B.A. in cognitive sci- Masakazu Sengoku ence from Lehigh University, Bethlehem, Pennsylvania in For contributions to graph theoretic research on circuits 1990, and the M.S. in electri- and communication network systems. cal engineering, the M.A. in Masakazu Sengoku received the B.E. degree from Niigata mathematics, and the Ph.D. in University, Japan, in 1967, and the M.E. and Dr. Eng. degrees electrical engineering from the from Hokkaido University in University of California, Ber- 1969 and 1972, respectively. In keley in 1991, 1994 and 1995 1972, he joined the staff at respectively. He was a Hokkaido University as re- postdoctoral researcher at the search associate. In 1978, he University of California, Berkeley before he joined IBM T. J. was associate professor in the Watson Research Center where he is currently research staff Department of Information En- member in the mathematical sciences department. He has co- gineering, Niigata University, taught graduate level seminars on nonlinear dynamics at Co- where he is presently professor. lumbia University and City University of New York. His re- His research interests include search interests include synchronization in chaotic systems, digi- network theory, graph theory, tal halftoning, image watermarking and multimedia security. transmission of information and Dr. Wu was associate editor of IEEE Transactions on Cir- mobile communications. He re- cuits and Systems, Part I. He is the author and co-author of more ceived the 1992, 1996, 1997 and than 50 refereed publications and the holder of 15 US patents. 1998 Best Paper Awards from IEICE, and IEEE ICNNSP Best Paper Award in 1996. He was chair of the IEICE Technical Group on Circuits and Systems in 1995. He is a member of the edito- Yilmaz Tokad rial boards of ACM, URSI, Wireless Networks, and Baltzer Sci- For lifetime contributions to research and education in ence Pub. He is vice-president of the Communication Society, system science and circuit theory. IEICE, for 2000–2001. He is author and coauthor of more than Dr. Yilmaz Tokad received the B.S. and M.S. degrees in 170 scientific publications and four books and holds several pat- electrical engineering from the Technical University of Istanbul ents. He was technical program member the 51st IEEE VTC, (ITU) in 1951 and 1952. He worked for Turkish PTT as an 2000; technical program co-chair, IEEE Region 10 Conference engineer for about a year and then taught at ITU. He received (TENCON’99); technical program chairman, IEEE/IEICE/IEEK the Ph.D. in electrical engineering from Michigan State Uni- ITC-CSCC’99; and general vice chair, IEEE International Sym- versity (MSU) in 1959 and served at MSU until 1971 as a pro- posium on Multi-Dimensional Mobile Communications fessor. He continued his research and teaching activities at (MDMC’98), November 1998, among others.

43 CallsCalls forfor PapersPapers andand ParticipationParticipation ICME2001 IEEE International Conference on Multimedia and Expo August 22–25, 2001 , Tokyo, Japan First IEEE Conference on Nanotechnology

October 28–30, 2001 The International Conference on Multimedia & Expo in Tokyo is the first ICME in the 21st century and the 2nd ICME co-organized by the four IEEE societies—the Circuits Outrigger Wailea Resort, Maui, Hawaii, USA and Systems Society, the Computer Society, the Signal Processing Society, and the Com- munications Society—in joint collaboration, so as to unify all IEEE as well as non-IEEE The First IEEE Conference on Nanotechnology will be held in Maui, Hawaii, USA multimedia related activities under the common umbrella of ICME. from October 28 (Sunday) to October 30 (Tuesday), 2001. The state-of-the-art technical achievements on all aspects of nanotechnology will be reported. General Chair Nobuyoshi Terashima IEEE/RSJ International Conference on Intelligent Robots and Systems: IROS 2001 ICME2001 office will be held at the same place from October 29 to November 3, 2001 (http://www.icra- Waseda University iros.com/iros2001). Global Information and [email protected] General Co-Chairs: Telecommunication Institute, Japan Toshio Fukuda, Nagoya University, [email protected] [email protected] Conference Webpage Robert D. Shull, NIST, [email protected] Technical Program Chair Jun Ohya (Waseda U./ Japan) http://www.giti.waseda.ac.jp/ICME2001/ Program Chair: Clifford Lau, ONR, [email protected] For more information, please check our World Wide Web site. The 2002 IEEE Asia-Pacific Conference on Circuits and Systems http://www.mein.nagoya-u.ac.jp/IEEE-NANO 2001 IEEE International Workshop on Behavioral Modeling and Simulation BMAS 2001 October 11–12, 2001 October 28–31, 2002 Denpasar, Bali, Indonesia Fountaingrove Inn Santa Rosa, California, USA http://apccas2002.itb.ac.id http://www.vlsi.ss.titech.ac.jp/apccas2002/ http://www.esat.kuleuven.ac.be/~bmas The 2002 IEEE Asia-Pacific Conference on Circuits and Systems (APCCAS’02) is For general information, contact: the sixth in the series of biennial Asia-Pacific Conferences sponsored by the IEEE Cir- Prof. Alan Mantooth, BMAS 2001 General Chair cuits and Systems Society. Department of Electrical Engineering University of Arkansas Technical Papers Author’s Schedule Fayetteville, AR 72701 Deadline for Submission of Papers: March 15, 2002 E-mail: [email protected] Notification of Acceptance: June 15, 2002 Deadline for Submission of Final Papers: August 15, 2002 CALL FOR PARTICIPATION Please use our website above for detailed list of areas, specification of the camera- ready format and electronic submission of papers. For further information contact Prof. Soegijardjo Soegijoko the APCCAS2002 Gen- A Joint Conference in 2002 eral Chair at the Conference Secretariat, Electronics Laboratory, Department of Electri- cal Engineering, Institut Teknologi Bandung, Jalan Ganesha 10 Bandung 40132 INDO- The Seventh Asia and South Pacific The Fifteenth International NESIA; Tel: (62–22) 253 4117; Fax: (62–22) 250 1895; E-mail: [email protected] Design Automation Conference Conference on VLSI Design January 7–11 2002, Bangalore, India ISLPED’01 GENERAL CHAIR International Symposium on Low Power Electronics and Design Sunil D. Sherlekar Sasken Communication Technologies Ltd, India Huntington Beach, California August 6–7, 2001 [email protected] General Chair: Program Co-Chairs: For more information, please visit the website: Enrico Macii, Politecnico di Torino Vivek De, Intel Corporation http://aspdac.com, or http://vlsi.ccrl.nj.nec.com, or [email protected] [email protected] http://vlsidesignindia.com Mary Jane Irwin, Penn State University [email protected] ACM/IEEE ISLPED Homepage: http://www.cse.psu.edu/~islped ASICON 2001 International Conference on Computer Aided Design A Conference for the EE CAD Professional The 4th INTERNATIONAL CONFERENCE ON ASIC November 4–8, 2001 Doubletree Hotel San José, California October 22–25, 2001, Hotel Equatorial, Shanghai, China General Co-Chairs: Qianling Zhang, Richard.M.M.Chen, Walter H.Ku General Chair: Rolf Ernst Conference Manager: Kathy MacLennan Dr. Albert Wang (Overseas) (P. R. China) Ms. Huihua Yu Tech. Univ. of Braunschweig MP Associates, Inc. 3255-6 Scott Blvd 220 Han Dan Road Inst. Furdatenverarbeitungsanlagn 5305 Spine Rd., Ste. A Santa Clara, CA 95054 Department of Electronic Engineering Hans-sommer-str. 66 Boulder, CO 80301 USA Fudan University, Shanghai, 200433 P. R. China Braunschweig D38106, Germany 303-530-4562 Tel: 408-327-7306 Tel: +86-21-65643761 49-531-391-3730; fax: 49-531-391-4587 fax: 303-530-4334 Fax: 408-986-8919 Fax: +86-21-65648267 [email protected] [email protected] [email protected] [email protected] Conference Web site: http://www.iccad.com For further information please visit our web site: http://www.cie-china.org 44 2001 IEEE CIRCUITS AND SYSTEMS SOCIETY ROSTER** (CAS-04, Division I) Division I Director President President – Elect Past President R. W. Wyndrum, Jr. H. C. Reddy J. A. Nossek B. J. Sheu +1 732 949 7409 + 1 562 985 5106 +49 89 2892 8501 + 1 408 562 9168 Ext 223 Administrative Vice President Vice President – Conferences Vice President – Publications Vice President – Technical Activities I. N. Hajj G. A. De Veirman M. N. S. Swamy M. E. Zaghloul + 1 961 347 952 +1 714 890 4118 +1 514 848 3091 +1 202 994 3772 Vice President – Regions 1-7 Vice President – Region 8 Vice President – Region 9 Vice President – Region 10 Advisor to President W. B. Mikhael A. C. Davies J. E. Cousseau N. Fujii G. De Micheli +1 407 823 3210 +44 20 7848 2441 +54 291 459 5153 +81 3 5734 2561 +1 650 725 3632 Elected Members of the Board of Governors 2001 2002 2003 Ex Officio J. Cong G. G. E. Gielen A. G. Andreou Chairman of Sponsored or E. G. Friedman R. Gupta R. J. Marks II Co-Sponsored Conferences I. Galton M. Hasler M. Pedram Standing Committee Chairs M. Green P. Pirsch L. Trajkovic Technical Committee Chairs C. Toumazou H. Yasuura E. Yoffa Chapter Chairpersons Editors of Society Sponsored Transactions & Magazines Society Administrator Division I Director B. Wehner D. Senese +1 219 871 0210 M. Ward-Callan

STANDING COMMITTEE CHAIRS Sensors and Micromachining M. E. Zaghloul Awards B. J. Sheu Visual Signal Processing and Communication J. Ostermann Constitution and Bylaws K. Thulasiraman VLSI Systems and Applications R. Sridhar Fellows W.-K. Chen REPRESENTATIVES Nominations G. S. Moschytz Neural Networks Council M. Ahmadi, A.R. Stubberud Parliamentarian R. J. Marks II Sensors Council M. E. Zaghloul, G. Barrows CONFERENCE ACTIVITIES G. A. De Veirman Solid-State Circuits Society E. Sanchez-Sinencio 2002 APCCAS S. Soegijoko Society on Social Implications of Technology R. W. Newcomb 2001 DAC J. Rabaey Society Education Chairs Committee J. Choma 2001 ICCAD R. Ernst IEEE Press M. A. Bayoumi 2001 ICCD S. Kundu IEEE TAB Magazines Committee W. H. Wolf 2001 ICECS J. Micallef IEEE TAB Nanotechnology Committee B. J. Sheu, C.-Y. Wu 2002 ISCAS S. Panchanathan/D. Allstot IEEE TAB New Technology 2001 MWSCAS R. L. Ewing Directions Committee C.-Y. Wu IEEE TAB Newsletters Committee M. Haenggi DISTINGUISHED LECTURER PROGRAM E. Yoffa IEEE TAB Nominations & MEMBERSHIP DEVELOPMENT AFFAIRS J. A. Nossek Appointments Committee G. S. Moschytz PUBLICATION ACTIVITIES M. N. S. Swamy IEEE TAB Periodicals Committee M. A. Bayoumi CAS Magazine Editor M. K. Sain IEEE TAB Transactions Committee M. Hasler Circuits and Devices Magazine Editor R. W. Waynant IEEE-USA R. de Figueiredo - Editor for CAS Society C.-Y. Wu IEEE-USA Professional Activities Committee Transactions on Circuits and Systems Part I: for Engineers (PACE) P. K. Rajan Fundamental Theory and Applications M. N. S. Swamy CORRESPONDING MEMBERS Transactions on Circuits and Systems Part II: TAB Awards and Recognition Committee L.-G. Chen Analog and Digital Signal Processing C. Toumazou Conference Publications Committee L. Goras Transactions on Computer Aided Design of P2SB/TAB Electronic Products and Integrated Circuits & Systems G. De Micheli Services Committee A. Ioinovici Transactions on Circuits and Systems for TAB Finance Committee K. Thulasiraman Video Technology W. Li TAB Periodicals Committee C.-S. Li Transactions on Multimedia M.-T. Sun TAB Periodicals Packages Committee J. Silva-Martinez Transactions on VLSI Systems E. G. Friedman TAB Periodicals Review Committee J. Vandewalle REGIONAL ACTIVITIES J. A. Nossek TAB Products Committee M. Ogorzalek TAB Society Review Committee J. Cousseau TECHNICAL ACTIVITIES M. E. Zaghloul TAB Strategic Planning and Review Committee R. Gupta Analog Signal Processing G. Cauwenberghs Cellular Neural Networks and Array Computing T. Roska 2001 Advisory and Consulting Group of Past CAS Society Presidents Circuits and Systems for Communications M. A. Bayoumi Professor Belle Shenoi Professor Ruey-Wen Liu Computer-Aided Network Design R. Gupta Professor George Moschytz Professor Wai-Kai Chen Digital Signal Processing P. S. R. Diniz Professor Leon Chua Professor Rolf Schaumann Multimedia Systems and Applications C.-S. Li Professor Rui de Figueiredo Professor Ming Liou Nanoelectronics and Gigascale Systems C.-Y. Wu Professor John Choma Professor Sanjit K. Mitra Neural Systems and Applications M. Ahmadi Professor Michael Lightner Nonlinear Circuits and Systems G. R. Chen ** To obtain full contact information for any one of the volunteers, please Power Systems and Power Electronics Circuits A. Ioinivici contact the Society Adminstrator, Barbara Wehner at [email protected].

45 Call for Papers 2001 Midwest Symposium on 2002 IEEE International Symposium Circuits and Systems On Circuits and Systems August 14–17, 2001 Fairborn, Ohio The Midwest Symposium on Circuits and Systems (MWSCAS) is one of the premier IEEE confer- ences presenting research in all as- pects of theory, design, and applica- tions of circuits and systems. In ad- dition to presenting new research results in core circuits and systems areas, this 2001 year conference will also be focusing on innovative information technology, intelligent exploitation http://www.iscas2002.com in microsystems and wireless information exchange. Six short courses, two plenary speakers, several invited sessions and a student paper contest are Sunday, May 26, 2002 – Wednesday, May 29, 2002 planned. Please check the symposium website for current information. Scottsdale Princess Resort Scottsdale, Arizona ISCAS 2002 Theme: General Chair General Co-Chair Circuits and Systems for Ubiquitous Computing Robert L. Ewing Harold W. Carter Air Force Research Laboratory University of Cincinnati The 2002 IEEE International Symposium on Circuits and Systems will be held in Scottsdale (Metropolitan Phoenix), Arizona on May 26–29, 2002. The Air Force Institute of Tech. [email protected] Symposium is sponsored by the IEEE Circuits and Systems Society (CASS). [email protected] Authors are expected to present their paper at the Symposium. At least one Technical Program Chair Technical Program Co-Chair author of each paper MUST register for the Symposium by January 2002 for Gary B. Lamont Tuna Tarim papers to be included in the program. Air Force Institute of Tech. Texas Instruments Author Deadlines gary.lamont@afit.mil [email protected] Submission of Tutorial proposal Friday, September 21, 2001 Submission of papers Monday, October 15, 2001 This symposium is for researchers, practitioners, and students inter- Notification of paper acceptance Friday, December 14, 2001 ested in circuits and systems, and their impact on computing, informa- Author registration Friday, January 25, 2002 tion, processing, design, and applications. This year's conference will Submission of FINAL paper Monday, January 28, 2002 include expanded focus on developing or using advanced aspects of theory, design, implementation, testing, and application of analog, digi- General Co-Chairs tal, and microwave circuits and systems. Besides the general presenta- Sethuraman (Panch) Panchanathan David J. Allstot tion session, six short courses, two plenary speakers, several invited ses- Arizona State University University of Washington sions and a student paper contest are planned. [email protected] [email protected] Dr. M. Ismail, Ohio State University, 205 Dreese Laboratory, Conference Management Services 2015 Neil Avenue, Columbus, OH 43210–1272. Billene Mercer Conference Management Services For further updated information, please see the symposium website. [email protected] http://www.ececs.uc.edu/~hcarter/mwscas

THE INSTITUTE OF ELECTRICAL & ELECTRONICS ENGINEERS, INC. 445 HOES LANE PISCATAWAY, NJ 08855