Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Vankka J. - Digital Synthesizers and Transmitters for Software Radio (2000)(en)

.pdf
Скачиваний:
68
Добавлен:
15.08.2013
Размер:
14.63 Mб
Скачать

234

Chapter 11

resulting in a limited range of filter characteristics.

Filter sharpening can be used to improve the response of a CIC filter. This technique applies the same filter several times to an input to improve both passband and stopband characteristics. The interested reader is referred to [Kwe97] for more details.

11.11 Pipelining/Interleaving

The pipelining/interleaving (P/I) technique can be used to reduce the amount of hardware when the same filtering operation is performed for several independent data streams [Par89], [Jia97]. The idea is to use the same hardware for all K data streams, as shown in Figure 11-11. This may be achieved with some additional logic. The data streams are interleaved in a serial form in time. The filtering operations are made with a K times higher clock frequency; after filtering, the data stream is de-interleaved back into K parallel data streams. This serial-in-time computation requires K register elements instead of one register in the filter structure. In addition to that, extra logic is needed in the interleave/de-interleave interfaces. The expense of this technique is the increased power consumption. This increase accrues because the number of delay elements remains same, but the clock frequency is doubled. The increased power consumption due to the doubled clock frequency in the other elements such as in the taps and adders is compensated by the halved number of elements.

REFERENCES

[Che82] P. R. Chevillat, and G. Ungerboeck, "Optimum FIR Transmitter and Receiver Filters for Data Transmission over Band-Limited Channels," IEEE Trans. Commun., Vol. 30, No. 8, pp. 1909-1915, Aug. 1982.

[Che99] C. Chen, and A. Willson, "A Trellis Search Algorithm for the Design of FIR Filters with Signed-Powers-of-Two Coefficients," IEEE Transactions on Circuits and Systems-II, Vol. 46, No. 1, pp. 29-39, Jan. 1999.

[Con99] G. A. Constantinides, P. Y. K. Cheung, and W. Luk, "Truncation noise in Fixed-Point SFGs [digital filters]," Electron. Lett., Vol. 35, No. 23,

Inter-

leaver K : 1

H(ZK)

K × Fs

Deinter- leaver 1 : K

Figure 11-11 Digital filtering of K independent signal sequences using a single filter.

Filter Design Algorithms

235

pp. 2012 -2014, Nov. 1999.

[Dua00] J. N. Duan, L. Ko, and B. Daneshrad, "A Highly Versatile Beamforming ASIC for Application in Broad-Band Fixed Wireless Access Systems," IEEE J. of Solid State Circuits, Vol. 35, No. 3, pp. 391-400, Mar. 2000.

[Fli94] N. J. Fliege, "Multirate Digital Signal Processing," John Wiley & Sons, 1994.

[Gra01] E. Grayver, and B. Daneshrad, "VLSI implementation of a 100- W Multirate FSK Receiver," IEEE J. of Solid State Circuits, Vol. 36, No. 11, pp. 1821-1828, Nov. 2001.

[Har96] R. Hartley, "Subexpression Sharing in Filters Using Canonic Signed Digit Multipliers," IEEE Transactions on Circuits and Systems-II, Vol. 43, No. 10, pp. 677-688, Oct. 1996.

[Hat01] M. Hatamian, "Efficient FIR Filter for High-Speed Communication," U. S. Patent 6,272,173, Aug. 7, 2001.

[Hau93] J. C. Hausman, R. R. Harnden, E. G. Cohen, and H. G. Mills, "Programmable Canonic Signed Digit Filter Chip," U. S. Patent 5,262,974, Nov. 16, 1993.

[Haw96] R. Hawley, B. Wong, T. Lin, J. Laskowski, and H. Samueli, "Design Techniques for Silicon Compiler Implementations of High-Speed FIR Digital Filters," IEEE Journal of Solid-State Circuits, Vol. 31, No. 5, pp.

656-667, May 1996.

[Hog81] E. B. Hogenauer, "An Economial Class of Digital Filters for Decimation and Interpolation," IEEE Trans. Acoust., Speech, Signal Process., Vol. ASSP-29, No. 2, pp. 155-162, Apr. 1981.

[Hwa75] S. Hwang, "On Monotonicity of Lp and lp Norms," IEEE Transactions on Acoustics, Speech, and Signal Processing, Dec. 1975, pp. 593-594. [Jia97] Z. Jiang, and A. Willson, "Efficient Digital Filtering Architectures Using Pipelining/Interleaving," IEEE Transactions on Circuits and SystemsII, Vol. 44, No. 2, pp. 110-119, Feb. 1997.

[Kho01] K. Y. Khoo, Z. Yu, and A. Willson, "Design of Optimal Hybrid Form FIR Filter," IEEE International Symposium on Circuits and Systems 2001, Vol. 2, pp. 621-624.

[Kho96] K. Y. Khoo, A. Kwentus, and A. N. Willson Jr., "A Programmable FIR Digital Filter using CSD Coefficients," IEEE J. of Solid State Circuits, Vol. 31, No. 6, pp. 869-874, June 1996.

[Kwe97] A. Y. Kwentus, Z. Jiang, and Alan N. Wilson, Jr. , "Application of Filter Sharpening to Cascaded Integrator-Comb Decimation Filters," IEEE Transactions on Signal Processing, Vol. 45, No. 2, pp. 457-467, 1997. [Lee96] H. R. Lee, C. W. Jen, and C. M. Liu, "A New Hardware-Efficient Architecture for Programmable FIR Filters," IEEE Transactions on Circuits and Systems-II, Vol. 43, No. 9, pp. 637-644, Sept. 1996.

236

Chapter 11

[Li93] D. Li, J. Song, and Y. Lim, "Polynomial-Time Algorithm for Designing Digital Filters with Power-of-Two Coefficients," In Proc. IEEE International symposium on Circuits and Systems 1993, Vol. 1, May 1993, pp. 84-

87.

[Man89] E. D. Man, and T. Noll, "Arrangement for Bit-Parallel Addition of Binary Numbers with Carry-Save Overflow Correction," U. S. Patent 4,888,723, Dec. 19, 1989.

[Mor95] F. Moreau de Saint-Martin, and P. Siohan, "Design of Optimal Lin- ear-Phase Transmitter and Receiver Filters for Digital Systems,'' IEEE International Symposium on Circuits and Systems, 1995, Vol. 2, pp. 885 -888.

[Nol91] T. Noll, "Carry-Save Architectures for High-Speed Digital Signal Processing," Journal of VLSI Signal Processing, Vol. 3, pp. 121-140, June 1991.

[Oh95] W. J. Oh, and Y. H. Lee, "Implementation of Programmable Multiplierless FIR Filters with Powers-of-Two Coefficients," IEEE Trans. Circuits and Systems II, Vol. 42, No. 8, pp. 553 - 556, Aug. 1995.

[Par89] K. Parhi, and D. Messerschmitt, "Pipeline Interleaving and Parallelism in Recursive Digital Filters-Part I: Pipelining Using Scattered LookAhead and Decomposition," IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. 37, No. 7, pp. 1099-1117, July 1989.

[Sam88] H. Samueli, "On the Design of Optimal Equiripple FIR Digital Filters for Data Transmission Applications," IEEE Transactions on Circuits and Systems, Vol. 35, No. 12, pp. 1542-1546, Dec. 1988.

[Sam89] H. Samueli, "An Improved Search Algorithm for the Design of Multiplierless FIR Filters with Powers-of-Two Coefficients," IEEE Transactions on Circuits and Systems, Vol. 36, No. 7, pp. 1044-1047, Jul. 1989.

[Sam91] H. Samueli, "On the Design of FIR Digital Data Transmission Filters with Arbitrary Magnitude Specifications," IEEE Transactions on Circuits and Systems, Vol. 38, No. 12, pp. 1563-1567, Dec. 1991.

[Tat90] L. R. Tate, "Truncated Product Partial Canonical Signed Digit Multiplier," U. S. Patent 4,967,388, Oct. 30, 1990.

[Vai87] P. Vaidyanathan, and T. Nguyen, "A "Trick'' for the Design of FIR Half-Band Filters," IEEE Transactions on Circuits and Systems, Vol. 34, No. 3, pp. 297-300, Mar. 1987.

[Vai93] P. P. Vaidyanathan, "Multirate Systems and Filter Banks," PrenticeHall, 1993.

[Whi00] B. A. White, and M. I Elmasry, "Low-Power Design of Decimation Filters for a Digital IF Receiver," IEEE Transactions on Very Large Scale Integration (VLSI), Vol. 8, No. 3, pp. 339-345, June 2000.

[Yag96] M. Yagu, A. Nishihara, and N. Fujii, "Design of FIR Digital Filters Using Estimates of Error Function over CSD Coefficient Space," IEICE Transactions on Fundamentals, Vol. E79-A, No. 3, pp. 283-290, Mar. 1996.

Filter Design Algorithms

237

[Zha01] T. Zhangwen, Z. Z. Zhanpeng, Z. Jie, and M. Hao "A High-Speed, Programmable, CSD Coefficient FIR Filter," in Proc. 13th Annual IEEE International ASIC/SOC Conference, Sept. 2001, pp. 397-400.

[Zha88] Q. Zhao, and Y. Tadokoro, "A Simple Design of FIR Filters with Powers-of-Two Coefficients," IEEE Transactions on Circuits and Systems, Vol. 35, No. 5, pp. 566-570, May 1988.

Chapter 12

12. RE-SAMPLING

The multi-standard modulator has to be able to accept data with different symbol rates. This fact leads to the need for a variable interpolator that performs a conversion between variable sampling frequencies. Furthermore, because the symbol rates in the different systems are not multiples of each other, the interpolation ratio has to be a rational number. An interpolator that is used for the sampling rate conversion is often called a re-sampler.

For example, the symbol rate of the GSM is 270.833 ksym/s and the symbol rate of the WCDMA is 3.84Msym/s. These rates are not multiples of each other. If these rates need to be converted to the sample rate 76.8MHz, the interpolation ratios for the GSM and WCDMA modes have to be 283.56 and 20, respectively. The ratio of the sampling frequency to the signal bandwidth determines, in part, the complexity of an interpolator with a rational number [Ves99]. Therefore, in these particular cases, an integer rate interpolator should be used to increase the pre-oversampling ratio [Ves99], [Sam00], [Cho01], [Ram84]. The sampling rate can first be raised using power-of-two ratios, for example. An efficient integer ratio interpolator is fairly easy to implement using a polyphase decomposition, as discussed in Section 11.8.

There are several methods of realizing a re-sampler with an arbitrary sampling rate conversion. In this chapter, the design of the polynomial-based interpolation filter using the Lagrange method is presented. Some other poly- nomial-based methods are also discussed.

Even though the concentration in this chapter is focused on the polyno- mial-based approaches, these are not the only ways to implement an interpolator with a rational number ratio. One well-known approach presented in most of the digital signal processing (DSP) books, in [Mit98], for example, is to cascade an interpolator of an integer ratio L and a decimator of an integer ratio M to form an interpolator with a rational number ratio L/M. Usually

240

Chapter 12

the multistage techniques are most efficient for these kind of solutions [Mit98]. In [Gra01], a similar approach in which the decimation ratio M can be a mixed integer/fractional number is used. In this system, the nearest interpolated value is selected as an output value and the exact value at the time instant is not calculated, which introduces some timing jitter. The methods described in [Che93] and [Che98] allow the decimation and interpolation factors to be prime numbers while preserving the multistage efficiency. Programmable cascaded integrator-comb filters are used for interpolation with a ratio of variable rational numbers in [Hen99] and [Ana01].

12.1 Interpolation for Timing Adjustment

When a sampling rate is to be converted, new samples using the known input samples need to be calculated. An interpolator calculates one interpolant y(kTout) at a time using N input samples. Let's first consider a continuoustime signal that is sampled with a sampling frequency Fs = 1/Ts. The output interpolants of an ideal interpolator are the same as if the original continu- ous-time signal had been sampled at the sampling frequency Fout = 1/Tout [Ves96a].

The time of the interpolant is governed by the fractional interval k in Figure 12-1 and defined by

x(nTs)

= y(nTout)

x(t)

k+ Ts

k Ts

(n - 2) Ts

(n - 1) Ts

n Ts

(n + 1) Ts

Figure 12-1 Interpolation in time-domain.

Re-Sampling

 

241

µk

(kTout nTs

)

(12.1)

Ts

 

 

 

 

where the fractional interval is the output sample interval and n output of the interpolator is

k [0,1), Ts the input sample interval, Tout the largest integer for which nTs kTout. The

 

 

 

N / 2

 

 

 

 

 

 

 

 

 

 

 

y(kT

out

)

¦

 

[(

 

 

i)T ] h

 

 

[(

i

µ

 

k

)T ]

(12.2)

 

 

 

 

 

 

 

 

 

 

 

 

i

N / 2+1

 

 

 

 

 

 

 

 

 

 

 

where x[(n i)Ts] and ha[( i+ k)Ts] for i = -N/2+1,...,N/2 are the N input samples, i.e. basepoints, and the N samples of the continuous impulse response of the interpolation filter, respectively.

12.2 Interpolation Filter with Polynomial-Based Impulse Response

As can be seen from Figure 12-1, in case of the arbitrary re-sampling rate, the fractional interval k has to be controlled on-line during the computation. Equation (12.2) states that the continuous-time impulse response ha(t) is a function of k. This means that the filter coefficients of the interpolation filter must change for every new value of k. This can be achieved by calculating the coefficients for all possible values of k and storing these into a lookup table. A new set of coefficients is then read into the filter structure every time k changes. Of course, this leads to a need for a huge coefficient memory, particularly when using long filters and a high resolution of the fractional interval. To avoid the very inefficient implementation using the coefficient memory, the filter coefficients can be calculated on-line using a polynomial-based interpolation.

The idea of the polynomial-based interpolator is that the continuous-time impulse response ha(t) of the interpolation filter is expressed in each interval of length Ts by means of polynomial [Ves99]

M

 

ha [(i + µk )Ts ] ¦cm (i) (µk ) m

(12.3)

m 0

 

for i = -N/2+1, -N/2+2,..., N/2 and for k [0,1). cm(i) denotes the polynomial coefficients for ha[(i + k)Ts] and M is the degree of polynomial. The coefficients cm(i) are constants and independent of k. This leads to the most attractive feature of the polynomial-based interpolation filters; they can be implemented very efficiently using the so-called Farrow structure [Far98] described in Section 12.3.

242

Chapter 12

 

12.2.1 Lagrange Interpolation

One of the best known polynomial-based interpolation methods is the Lagrange interpolation. An advantage of the Lagrange interpolation is that the polynomial coefficients can be determined in a closed form, and therefore no optimization is needed [Ves99]. On the other hand, this leads to a fixed frequency response that cannot be controlled by any other design parameter than the degree of the polynomial M. All the same, the Lagrange interpolators provide a fairly good approximation of the continuous-time signal at low frequencies, i.e. far away from half the sampling rate [Ves99]. It is also a maximally flat design of the interpolation filters or the fractional delay finite impulse response (FIR) filters [Väl95].

The approximating polynomial of degree M, denoted by ya(t), can be ob-

tained by using the following M + 1 time-domain conditions

 

ya

(iTs

) x(i), for n N + 1 ≤ i n + N

(12.4)

 

 

2

2

 

which means that the reconstructed signal value is exactly the value of the input signal at the sampling instants. To minimize the approximation error, only the middle interval nTs ≤ t (n + 1)Ts of the approximating polynomial is normally used. In order to uniquely determine the middle interval, the length of the filter N M + 1 has to be even. Even though not recommended in [Eru93], for instance, it is also possible to use the Lagrange interpolator with an odd N. In this case, the locations of the filter taps have to be changed when k passes the value 0.5 to maintain k within half a sample from the center point of the filter [Väl95]. This is practically difficult to implement and it also causes a discontinuity to the output signal [Väl95].

The impulse response ha(t) of the interpolation filter based on the Lagrange interpolation can be derived by using the conditions of Equation (12.4). It is a piecewise polynomial and can be expressed as (12.3),

M

 

ha [(i + µk )Ts ] ¦cm (i)(µk )m

(12.5)

m 0

 

for i = -N/2+1, -N/2+2,..., N/2. The so-called Lagrangian coefficient functions can be determined from the following equation

M

 

 

N / 2

j

x

 

¦

cm

(i)xm

(12.6)

 

i + j

m 0

 

 

j N / 2+1

 

 

 

 

 

 

ji

 

 

 

One of the possible derivations of (12.6) is shown in Appendix [Väl95]. The coefficients cm(i) of a cubic, i.e. third order M = 3, Lagrange interpolator are calculated in Table 12-1. The impulse and magnitude responses of the linear (M = 1) and cubic (M = 3) Lagrange interpolators are depicted in

Re-Sampling

243

Table 12-1 Filter coefficients for a cubic Lagrange interpolator.

 

m = 0

m = 1

m = 2

m = 3

cm(-1)

0

-1/3

1/2

-1/6

cm(0)

1

-1/2

-1

1/2

cm(1)

0

1

1/2

-1/2

cm(2)

0

-1/6

0

1/6

Figure 12-2 and Figure 12-3

respectively. One important property of the

interpolator, when using it as a sampling rate converter, is the phase linearity. Even though the underlying continuous time impulse response of the Lagrange interpolator is symmetric and thus linear phase, the phase response of the interpolator is not linear for all k. Fortunately, it is phase linear on average. The average phase linearity can also be seen from the phase responses of the cubic Lagrange interpolation filter for k =0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1 depicted in Figure 12-4.

12.3 Farrow Structure

The Farrow structure is a very efficient implementation structure for a poly- nomial-based interpolator. It suits well in real-time systems because it is fast and it has fixed coefficient values. A Farrow structure can be found for any polynomial-based filter [Ves99].

Substituting (12.5) into (12.2), the following is obtained

h(t)

IMPULSE RESPONSES

1

Linear

Cubic

0.8

0.6

0.4

0.2

0

-0.2

-2

-1

0

1

2

 

 

t/Ts

 

 

Figure 12-2 Impulse responses of linear and cubic Lagrange interpolators.

244

Chapter 12

 

FREQUENCY RESPONSES

 

 

 

 

 

 

 

0

 

 

 

Linear

 

 

 

 

Cubic

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−10

H(f) dB

−20

−30

−40

−50

−60

−70

−80

−90

0

1

2

f Ts

Figure 12-3 Magnitude responses of linear and cubic Lagrange interpolators.

 

 

 

 

N / 2

M

 

 

 

 

 

 

 

 

 

M

 

 

y(kT

out

)

¦

 

 

[(

 

 

 

 

i)T ] ¦

 

 

 

 

 

 

(i) (

µ

k

)m = ¦

 

(

 

) (

µ

 

k

)m

(12.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

N / 2+1

m 0

 

 

 

 

 

 

 

 

 

m 0

 

where

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N / 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

)

¦

 

 

[(

 

 

 

 

 

i)T ]

 

 

 

(i)

(12.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i N / 2+1

(12.7) and (12.8) indicate that the polynomial-based interpolation filter

Phase (degrees)

0

-50

-100

-150

-200

-250

-300

-350

-4000

0.1

0.2

0.3

0.4

0.5

Normalized frequency

Figure 12-4 Phase response of cubic Lagrange interpolator for different k values.

Соседние файлы в предмете Химия