FIR的英文资料.docx

FIR Digital Filter Design In chapter 9 we considered the design of IIR digital filters. For such filters, it is also necessary to ensure that the derived transfer function G(z) is stable. On the other hand, in the case of FIR digital filter design,the stability is not a design issue as the transfer function is a polynomial in z-1 and is thus always guaranteed stable. In this chapter, we consider the FIR digital filter design problem. Unlike the IIR digital filter design problem, it is always possible to design FIR digital filters with exact linear-phase. First ,we describe a popular approach to the design of FIR digital filters with linear-phase. We then consider the computer-aided design of linear-phase FIR digital filters. To this end, we restrict our discussion to the use of matlab in determining the transfer functions. Since the order of the FIR transfer function is usually much higher than that of an IIR transfer function meeting the same frequency response specifications, we outline two methods for the design of computationally efficient FIR digital filters requiring fewer multipliers than a direct form realization. Finally, we present a method of designing a minimum-phase FIR digital filter that leads to a transfer function with smaller group delay than that of a linear-phase equivalent. The minimum-phase FIR digital filter is thus attractive in applications where the linear-phase requirement is not an issue. 10.1 preliminary considerations In this section,we first review some basic approaches to the design of FIR digital filters and the determination of the filter order to meet the prescribed specifications. 10.1.1 Basic Approaches to FIR Digital Filter Design Unlike IIR digital filter design, FIR filter design does not have any connection with the design of analog filters. The design of FIR filters is therefore based on a direct approximation of the specified magnitude response,with the often added requirement that the phase response be linear. Recall a causal FIR transfer function H(z) of length N+1 is a polynomial in z（-1） of degree N: H(z) = hnzn n=0N (10.1)The corresponding frequency response is given byH (ej) = hne(-jn) n=0 N (10.2)It has been shown in section 5.3.1 that any finite duration sequence xn of length N+1 is completely characterized by N+1 samples of its discrete-time Fourier transform X(ej). As a result, the design of an FIR filter of length N+1 can be accomplished by finding either the impulse response sequence hn or N+1 samples of its frequency response H(e j). Also ,to ensure a linear-phase design, the conditionhn = h N-n ,must be satisfied. Two direct approaches to the design of FIR filters are the windowed Fourier series approach and the frequency sampling approach. We describe the former approach in Section 10.2. The second approach is treated in Problems 10.31 and 10.32. In section 10.3, we outline computer-based digital filter design methods. 10.1.2 Estimation of the Filter Order After the type of the digital filter has selected, the next step in the filter design process is to estimate the filter order should be the smallest integer greater than or equal to the estimated value. FIR Digital Filter Order Estimation For the design of lowpass FIR digital filters, several authors have advanced formulas for estimating the minimum value of the filter order N directly from the digital filter specifications: normalized passband edge angular frequency p , normalizef stopband edge angulars,bandpass peak ripplep，stopband peak ripples. Kaisers Formula. A rather simple formula developed by Kaiser Kai74 is given by N-20log10(ps)-13/14.6(s-p)/2We illustrate the application of the above formula in Example 10.1.Bellangers Formula. Another simple formula advanced by Bellanger is given by Bel81 N-2log10(10ps) /3(s-p)/ 2-1 .Its application is considered in Example 10.2. Hermanns Formula. The formula due to Hermann et al.Her73 gives a slightly more accurate value for the order and is given by ND(p,s)-F(p,s)*(s-p)/ 2 2/(s-p)/ 2Where D（p,s）= a1(log10p) 2+a2(log10p)+a3 log10s- a4(log10p) 2+a5(log10p)+a6And F(p,s) = b1 + b2log10p-log10s , With a1=0.005309, a2=0.07114, a3=-0.4761,a4=0.00266, a5=0.5941, a6=0.4278,b1=11.01217, b2=0.51244. The formula given in Eq.(10.5) is valid for ps . If p s , then the filter order formula to be used is obtained by interchanging p and s in Eq.(10.6a) and (10.6b).For small values of p and s , all of the above formulas provide reasonably close and p and s are large, Eq.(10.5) yields accurate results. On the other hand, when the values of p and s are large,Eq.(10.5) a more accurate value for the order. A Comparison of FIR Filter Order Formulas Note that the filter order computed in Examples 10.1, 10.2 and 10.3, using Eqs.(10.3),(10.3),and (10.5). Respectively ,are all different. Each of these three formulas provide only an estimate of the required filter order. The frequency response of the FIR filter designed using this estimated order may or may not meet the given specifications. If the specifications are not met, it is recommended that the filter order be gradually increased until the specifications are met. Estimation of the FIR filter order using MATLAB is discussed in Section 10.5.1. An important property of each of the above three formulas is that the estimated filter order N of the FIR filter is inversely proportional to the transition band width (s-p) and does not depend on the actual location of the transition band. This implies that a sharp cutoff FIR filter with a narrow transition band would be of very high order, whereas an FIR filter with a wide transition band will have a very low order. Another interesting property of Kaisers and Bellangers formulas is that the order depends on the product ps . This implies that if the values of p and s are interchanged, the order remains the same. To compare the accuracy of the the above formulas, we estimate using each formula the order of three linear-phase lowpass FIR filters of known order, bandedges, and ripples. The specifications of the three filters are as follows: Filter No.1: p = 0.10625,s=0.14375,p=0.0224,s = 0.0001