260 likes | 509 Views
Filtros Digitais. Diagramas de bloco e grafos de fluxo de sinal Estruturas de filtros IIR Projeto de filtro FIR. 6.1.1 Filter Transfer Function. Linear constant-coefficient difference equations (LCCDEs). Taking the two sided Z-transform we have
E N D
Filtros Digitais • Diagramas de bloco e grafos de fluxo de sinal • Estruturas de filtros IIR • Projeto de filtro FIR
6.1.1 Filter Transfer Function • Linear constant-coefficient difference equations (LCCDEs). • Taking the two sided Z-transform we have • H(z) is a filter transfer function, and can be used to describe any digital filter. • Expressed as a diagram or signal flow graph • Implemented as a digital circuit
6.1.2 Block Diagram and Signal Flow Graph • A signal flow graph representation of LCCDE is same as a block diagram, and it is a network of directed branches that connect at nodes which are variables. Block diagram representation of a first- order digital filter. Structure of the signal flow graph. Structure of the signal flow graph with the delay branch indicated by z-1.
6.1.3 Block Diagram: Direct Form I Realization This structure (non-canonic, direct form I) can be too sensitive to finite word-length errors (quantization errors) – errors are summed, fed back and re-amplified over and over.
6.1.4 Block Diagram: Direct Form II Realization This structure (canonic direct from, direct form II) requires less delay elements. The minimum number of delays is max(N, M).
6.1.5 Example of LTI Implementation Consider the LTI system with system (transfer) function we have two implementation as follows
6.1.6 Signal Flow Graph: Direct Forms • Given the LCCDE Signal Flow Chart of Direct From I Signal Flow Chart of Direct From II
6.2.1 Structure of IIR: Cascade Form • By factoring the numerator and denominator we can write which can be drawn as a cascade of smaller sections: • Advantages: Smaller sections – less feedback error. • Disadvantages: Errors fed from section-to-section.
6.2.2 Parallel Realization • By performing a partial fraction expansion we can write which can be drawn as a parallel sum of smaller sections • Advantages: smaller sections- less feedback error, and error confined to each section.
6.2.3 Structure of IIR: Example (Cascade) Given a two-order system • Cascade Structure (Not unique)
6.2.4 Structure of IIR: Example (Parallel) Given a two-order system • Parallel Structure (Not unique) Parallel-form structure using second- order system (form I) Parallel-form structure using first-order system
6.3.1 Digital Filter Design • Given a set specifications or stated constraints on • magnitude spectrum • phase spectrum • Find where, • The constraints may include • zero, small, or linear phase • specific bandlimit within a passband • amount of ripple within a passband • amount of ripple within a stopband • sharpness of transitions between passband/stopband • filter order K, M.
6.3.2 Finite Impulse Response (FIR) Filter Design • Here it is assumed that • Hence • And so the unit pulse response of the filter is clearly: • Problem: Given specifications on and , • find • FIR filters are often called non-recursive for obvious reasons.
6.3.3 FIR Filter: Advantages and Disadvantages • Advantages: • Always stable (assume non-recursive implementation). • Quantization noise is not much of a problem. • Can be designed to have exact linear phase even when causal, while meeting a prescribed phase to arbitrary accuracy. • Design complexity generally linear. • Transients have a finite duration. • Disadvantages: • A high-order filter is generally needed to satisfy the stated specification – so more coefficients are needed with more storage and computation.
6.3.4 FIR Filter Design: Linear Phase Condition • Definition: The digital filter is linear phase if for some real number C. If C=0, then the filter has zero- phase, which is only possible when the filter is non-causal. • Achieving linear phase is quite important in applications where is desirable not to distort the signal phase much –i.e., where the frequency locations are critical, such as speech signals. • Many applications benefit be the linear phase thought as • shaping frequencies according to the magnitude spectrum. • Time-shifting the response by an amount -C
6.3.5 Linear Phase Condition • Theorem: a causal FIR filter with unit pulse response is linear phase if h(n) is even symmetric: Proof: Suppose M is odd. Then which finishes the proof, why? • Consider case of Meven to be an exercise.
6.4.1 FIR Filter Design: Windowing • Goal: Design an FIR digital filter with M+1 coefficients that approximates a desired frequency response with • Usually d(n) cannot be realized for some reasons. • d(n) has infinite duration if contains discontinuities; • If d(n) is non-causal and we want it causal; • If d(n) is longer than can be computed efficiently; • It’s generally desirable to have few coefficients; • Windowing is the simplest approach to FIR filter design. One can proceed naively, and thus obtain poor results. But with little care (basic windowing strategies), windowing can be very effective.
6.4.2 General Windowing Approach • Define where • Then designed filter then has frequency response • Observations: We desire conflicting goals • be time-limited to • be spectrally localized – impulse-like, if
6.4.3 Truncation Windowing • Rectangular window: • The designed filter has frequency response where • Is the frequency response a good approximation to the desired frequency response ? • Actually, for M give is optimal in the mean square sense (MSE). • However, while rectangular windowing does the best global MSE job, it suffers dramatically at frequencies.
6.4.4 Triangle (Bartlett) Windowing • Suppose that which • Note that (ignoring the shift) is a positive function, hence must rise monotonically at a jump discontinuity (why?). • In the prior example, using the triangular window gives an approximation with smooth, but wider transition.
6.4.5 Windowing: Trade-off Ripples vs. Transition Width • Rectangular window has a sharp transition but severe ripple. • Triangular window has no ripple but a very wide transition.
6.4.6 Other Windows • Other windows attempt to optimize this trade-off. Widely used windows that give intermediate results are: • Hamming Window: • Hanning Window: • Blackman Window:
6.4.7 Windowing Comparisons Rectangular: transition width is optimized. Blackman: Ripple is minimized.
6.4.8 Kaiser Window Design • Here Where and represents the zeroth-order modified Bessel function of the first kind, and there are two important parameters: M, . • For M held constant, increasing reduces sidelobe but increase mainlobe width. • For held constant, increasing M reduces mainlobe width but does not affect sidelobes much. • Kaiser developed an empirical but careful design procedure for windowing a filter having sharp discontinuity (e.g. an ideal LPF).