通信工程畢業(yè)設(shè)計-基于matlab的iir數(shù)字濾波器設(shè)計(更新版)

　　

【正文】 step in the development of a digital filter is the determination of a realizable transfer function G(z) approximating the given frequency response specifications. If an IIR filter is desired,it is also necessary to ensure that G(z) is stable. The process of deriving the transfer function G(z) is called digital filter design. After G(z) has been obtained, the next step is to realize it in the form of a suitable filter structure. In chapter 8,we outlined a variety of basic structures for the realization of FIR and IIR transfer functions. In this chapter,we consider the IIR digital filter design problem. The design of FIR digital filters is treated in chapter 10. First we review some of the issues associated with the filter design problem. A widely used approach to IIR filter design based on the conversion of a prototype analog transfer function to a digital transfer function is discussed next. Typical design examples are included to illustrate this approach. We then consider the transformation of one type of IIR filter transfer function into another type, which is achieved by replacing the plex variable z by a function of z. Four monly used transformations are summarized. Finally we consider the puteraided design of IIR digital filter. To this end, we restrict our discussion to the use of matlab in determining the transfer functions. preliminary considerations There are two major issues that need to be answered before one can develop the digital transfer function G(z). The first and foremost issue is the development of a reasonable filter frequency response specification from the requirements of the overall system in which the digital filter is to be employed. The second issue is to determine whether an FIR or IIR digital filter is to be designed. In the section ,we examine these two issues first . Next we review the basic analytical approach to the design of IIR digital filters and then consider the determination of the filter order that meets the prescribed specifications. We a34 lso discuss appropriate scaling of the transfer function. Digital Filter Specifications As in the case of the analog filter,either the magnitude and/or the phase(delay) response is specified for the design of a digital filter for most applications. In some situations, the unit sample response or step response may be specified. In most practical applications, the problem of interest is the development of a realizable approximation to a given magnitude response specification. As indicated in section , the phase response of the designed filter can be corrected by cascading it with an allpass section. The design of allpass phase equalizers has received a fair amount of attention in the last few years. We restrict our attention in this chapter to the magnitude approximation problem only. We pointed out in section that there are four basic types of filters,whose magnitude responses are shown in Figure . Since the impulse response corresponding to each of these is noncausal and of infinite length, these ideal filters are not realizable. One way of developing a realizable approximation to these filter would be to truncate the impulse response as indicated in Eq.() for a lowpass filter. The magnitude response of the FIR lowpass filter obtained by truncating the impulse response of the ideal lowpass filter does not have a sharp transition from passband to stopband but, rather, exhibits a gradual rolloff. Thus, as in the case of the analog filter design problem outlined in section , the magnitude response specifications of a digital filter in the passband and in the stopband are given with some acceptable tolerances. In addition, a transition band is specified between the passband and the stopband to permit the magnitude to drop off smoothly. For example, the magnitude )( ?jeG of a lowpass filter may be given as shown in Figure . As indicated in the figure, in the passband defined by 0 p???? , we require that the magnitude approximates unity with an error of p?? ,., ppjp fo reG ???? ? ????? ,1)(1 . In the stopband, defined by ??? ??s ,we require that the magnitude approximates 35 zero with an error of is,? .e., ,)( sjeG ?? ? for ??? ??s . The frequencies p? and s? are , respectively, called the passband edge frequency and the stopband edge frequency. The limits of the tolerances in the passband and stopband, p? and s? , are usually called the peak ripple values. Note that the frequency response )( ?jeG of a digital filter is a periodic function of ? ,and the magnitude response of a realcoefficient digital filter is an even function of ? . As a result, the digital filter specifications are given only for the range ????0 . Digital filter specifications are often given in terms of the loss function,)(lo g20)( 10 ??? jeG?? , in dB. Here the peak passband ripple p? and the minimum stopband attenuation s? are given in dB,., the loss specifications of a digital filter are given by dBpp )1(lo g20 10 ?? ??? , dBss )(lo g20 10 ?? ?? . Preliminary Considerations As in the case of an analog lowpass filter, the specifications for a digital lowpass filter may alternatively be given in terms of its magnitude response, as in Figure . Here the maximum value of the magnitude in the passband is assumed to be unity, and the maximum passband deviation, denoted as 1/ 21 ?? ,is given by the minimum value of the magnitude in the passband. The maximum stopband magnitude is denoted by 1/A. For the normalized specification, the maximum value of the gain function or the minimum value of the loss function is therefore 0 dB. The quantity max? given by dB)1(lo g20 210m a x ?? ?? Is called the maximum passband attenuation. For p? ?? 1, as is typically the case, it36 can be shown that pp ??? 2)21(lo g20 10m a x ???? The passband and stopband edge frequencies,