Design of Digital Filters

Total Page:16

File Type:pdf, Size:1020Kb

Design of Digital Filters Chapter 8 Design of Digital Filters Contents Overview . 8.2 General Considerations . 8.2 Design of FIR Filters . 8.4 Design of linear-phase FIR filters using windows . 8.7 Time-delay in desired response . 8.9 Sidelobes . 8.10 ◦ Hilbert transform (90 phase shift) . 8.12 Design of FIR filters by frequency sampling . 8.12 “Optimum” equiripple linear-phase FIR filters . 8.14 Comparison of FIR methods . 8.17 Design of IIR Filters from Analog Filters . 8.18 IIR filter design by bilinear transformation . 8.19 Frequency Transformations . 8.24 Design of Digital Filters Based on Least-Squares Method . 8.24 Summary . 8.24 8.1 8.2 c J. Fessler, May 27, 2004, 13:18 (student version) So far our treatment of DSP has focused primarily on the analysis of discrete-time systems. Now we finally have the analytical tools to begin to design discrete-time systems. All LTI systems can be thought of as filters, so, at least for LTI systems, to “design” a system means to design a digital filter. (The design of nonlinear or time-varying systems is generally more complicated, and often more case specific.) Goal: given desired • magnitude response j Hd(!) j • phase response \ Hd(!) • tolerance specifications (how far from ideal?), we want to choose the filter parameters N, M, fakg, fbkg such that the system function M b z−1 (z − z ) H(z) = k=0 k = g i i ; PN −1 Q (z − p ) k=0 akz j j P Q where a0 = 1, yields a frequency response H(!) ≈ Hd(!). Rational H(z), so LTI system described by a constant-coefficient difference equation, so can be implemented with finite # of adds, multiplies, and delays. In other words, filter design means choosing the number and locations of the zeros and poles, or equivalently the number and values of the filter coefficients, and thus H(z), h[n], H(!). Overview • N = 0, FIR or all-zero. linear passband phase. • N > 0, IIR. lower sidelobes for same number of coefficients. 8.1 General Considerations • Ideally would like N and M or N + M to be as small as possible for minimal computation / storage. • Causal (for now) • Poles inside unit circle for stability 8.1.1 Causality We will focus on designing causal digital filters, since those can be implemented in real time. Noncausal filter design (e.g., for off-line applications) is much easier and many of the same principles apply anyway. An LTI system is causal iff • input/output relationship: y[n] depends only on current and past input signal values. • impulse response: h[n] = 0 for n < 0 • system function: number of finite zeros ≤ number of finite poles • frequency response: What can we say about H(!)? Fact: if h[n] is causal, then • Paley-Wiener Theorem: H(!) cannot be exactly zero over any band of frequencies. (Except in the trivial case where h[n] = 0.) • Furthermore, jH(!)j cannot be flat (constant) over any finite band. • HR(!) and HI (!) are Hilbert transform pairs. Therefore they are not independent. Hence magnitude and phase response are interdependent. Thus those ideal filters with finite bands of zero response cannot be implemented with a causal filter. Instead, we must design filters that approximate the desired frequency response Hd(!). c J. Fessler, May 27, 2004, 13:18 (student version) 8.3 8.1.2 Characteristics of practical frequency-selective filters No perfectly flat regions Fact: since causal filters cannot have a band of frequencies with zero response, nor can they have any band of frequencies over which the frequency response is a constant. Proof by contradiction. Suppose H(!) = c for !1 ≤ ! ≤ !2, with corresponding impulse response h[n]. Now define a new filter g[n] = h[n] −c δ[n]. Then certainly g[n] is also causal. But G(!) = H(!) −c = 0 for !1 ≤ ! ≤ !2, which is impossible if g[n] is causal. Causal filters cannot have a band of frequencies with zero response. Nor can they have an infinitely sharp transition between the passband and the stopband. Nor can they have perfectly flat passband. So a typical realistic magnitude response looks like the following. jH(!)j 1 + δ1 1 Passband Ripple 1 − δ1 Transition Band Passband PSfrag replacements Stopband δ2 Stopband Ripple !p !s π ! Practical filter design means choosing δ1, δ2, !c and !s, and then designing N, M, fakg, fbkg to satisfy those requirements. Often one must iterate. Often we plot jH(!)j using dB, i.e., 20 log10 jH(!)j, and express the ripple in dB as well. Note that stopband ripple is not defined peak-to-peak, since the highest magnitude response in the stopband is more important than how wiggly the response is in the stopband. Example application: CD digital crossover Separate woofer and tweeter signals digitally inside CD player, rather in analog at speaker. Why is passband ripple tolerable in this context (distortion)? Room acoustics act as another filter. Speaker response is not perfectly flat. Keep filter ripple below these other effects. Now imagine that we have specified δ1, δ2, !c and !s, and we wish to design N, M, fakg, fbkg to satisfy those requirements. As mentioned above, we have two broad choices: FIR and IIR. We focus first on FIR. 8.4 c J. Fessler, May 27, 2004, 13:18 (student version) 8.2 Design of FIR Filters An FIR filter of length M is an LTI system with the following difference equation1: M−1 y[n] = bk x[n − k] : kX=0 Note that the book changes the role of M here. Earlier, when discussing rational system functions, M was the number of zeros. Now M is the number of “nonzero” elements of h[n], which corresponds to at most M − 1 zeros. (More precisely, we assume bM−1 6= 0 and b0 6= 0, but some of the coefficients in between could be zero.) M−1 The problem: given δ1, δ2, !c, and !s, we wish to choose M and fbkgk=0 to achieve that goal. We focus on linear-phase FIR filters, because if linear phase is not needed, then IIR is probably preferable anyway. We focus on lowpass filters, since transformations can be made to form highpass, bandpass from lowpass, as discussed previously. Impulse response Clearly for an FIR filter b ; n = 0; : : : ; M − 1 h[n] = n 0; otherwise: So rather than writing everything in terms of bk's, we write it directly in terms of the impulse response h[n]. In fact, for FIR filter design we usually design h[n] directly, rather than starting from a pole-zero plot. (An exception would be notch filters.) 8.2.1 Symmetric and antisymmetric FIR filters I focus on the symmetric case. System function: M−1 H(z) = h[n] z−n: nX=0 How do we make a filter have linear phase? We previously answered this in the pole-zero domain. Now we examine it in the time domain. An FIR filter has linear phase if h[n] = h[M − 1 − n]; n = 0; 1; : : : ; M − 1. Example. For M = 5: h[n] = fb0; b1; b2; b1; b0g : Such an FIR filter is called symmetric. Caution: this is not “even symmetry” though in the sense we discussed previously. This is related, but not exactly the same as circular symmetry. Example. For M = 3 and h[n] = 1=2; 1; 1=2 . Does this filter have linear phase? Is it lowpass or highpass? n o 1 1 1 1 H(!) = + e−|! + e−|2! = e−|! e|! + 1 + e−|! = e−|! (1 + cos !); 2 2 2 2 so since 1 + cos ! ≥ 0, \ H(!) = −!; which is linear phase. So it works for this particular example, but why does the symmetry condition ensure linear phase in general? 1 M M−1 Caution. At this point the book switches from k=0 to k=0 apparently. This inconsistent with MATLAB, so there are “M − 1” factors that appear frequently in the MATLAB calls. I think that MATLABPis consistentP and the book makes an undesirable switch of convention here. c J. Fessler, May 27, 2004, 13:18 (student version) 8.5 Symmetric real FIR filters, h[n] = h[M − 1 − n]; n = 0; : : : ; M − 1, are linear phase. Proof. Suppose M is even: M−1 H(z) = h[n] z−n Xn=0 M=2−1 M−1 = h[n] z−n + h[n] z−n (split sum) Xn=0 n=XM=2 M=2−1 M=2−1 0 = h[n] z−n + h[M − 1 − n0] z−(M−1−n ) (n0 = M − 1 − n) Xn=0 nX0=0 M=2−1 M=2−1 = h[n] z−n + h[n] z−(M−1−n) (symmetry of h[n]) Xn=0 Xn=0 M=2−1 = h[n] z−n + z−(M−1−n) (combine) Xn=0 h i M=2−1 = z−(M−1)=2 h[n] z(M−1)=2−n + z−((M−1)=2−n) (split phase). Xn=0 h i Thus the frequency response is H(!) = H(z) z=e|! M=2−1 = e−|!(M−1)=2 h[n] e|!((M−1)=2−n) + e−|!((M−1)=2−n) Xn=0 h i M=2−1 M − 1 = 2 e−|!(M−1)=2 h[n] cos ! − n 2 Xn=0 −|!(M−1)=2 = e Hr(!) M=2−1 M − 1 H (!) = 2 h[n] cos ! − n : (Real since h[n] is real.) r 2 Xn=0 Phase response: ! M−1; ! > \ − 2 Hr( ) 0 H(!) = M−1 −! 2 + π; Hr(!) < 0: The case for M odd is similar, and leads to the same phase response but with a slightly different Hr(!).
Recommended publications
  • J-DSP Lab 2: the Z-Transform and Frequency Responses
    J-DSP Lab 2: The Z-Transform and Frequency Responses Introduction This lab exercise will cover the Z transform and the frequency response of digital filters. The goal of this exercise is to familiarize you with the utility of the Z transform in digital signal processing. The Z transform has a similar role in DSP as the Laplace transform has in circuit analysis: a) It provides intuition in certain cases, e.g., pole location and filter stability, b) It facilitates compact signal representations, e.g., certain deterministic infinite-length sequences can be represented by compact rational z-domain functions, c) It allows us to compute signal outputs in source-system configurations in closed form, e.g., using partial functions to compute transient and steady state responses. d) It associates intuitively with frequency domain representations and the Fourier transform In this lab we use the Filter block of J-DSP to invert the Z transform of various signals. As we have seen in the previous lab, the Filter block in J-DSP can implement a filter transfer function of the following form 10 −i ∑bi z i=0 H (z) = 10 − j 1+ ∑ a j z j=1 This is essentially realized as an I/P-O/P difference equation of the form L M y(n) = ∑∑bi x(n − i) − ai y(n − i) i==01i The transfer function is associated with the impulse response and hence the output can also be written as y(n) = x(n) * h(n) Here, * denotes convolution; x(n) and y(n) are the input signal and output signal respectively.
    [Show full text]
  • Moving Average Filters
    CHAPTER 15 Moving Average Filters The moving average is the most common filter in DSP, mainly because it is the easiest digital filter to understand and use. In spite of its simplicity, the moving average filter is optimal for a common task: reducing random noise while retaining a sharp step response. This makes it the premier filter for time domain encoded signals. However, the moving average is the worst filter for frequency domain encoded signals, with little ability to separate one band of frequencies from another. Relatives of the moving average filter include the Gaussian, Blackman, and multiple- pass moving average. These have slightly better performance in the frequency domain, at the expense of increased computation time. Implementation by Convolution As the name implies, the moving average filter operates by averaging a number of points from the input signal to produce each point in the output signal. In equation form, this is written: EQUATION 15-1 Equation of the moving average filter. In M &1 this equation, x[ ] is the input signal, y[ ] is ' 1 % y[i] j x [i j ] the output signal, and M is the number of M j'0 points used in the moving average. This equation only uses points on one side of the output sample being calculated. Where x[ ] is the input signal, y[ ] is the output signal, and M is the number of points in the average. For example, in a 5 point moving average filter, point 80 in the output signal is given by: x [80] % x [81] % x [82] % x [83] % x [84] y [80] ' 5 277 278 The Scientist and Engineer's Guide to Digital Signal Processing As an alternative, the group of points from the input signal can be chosen symmetrically around the output point: x[78] % x[79] % x[80] % x[81] % x[82] y[80] ' 5 This corresponds to changing the summation in Eq.
    [Show full text]
  • Designing Filters Using the Digital Filter Design Toolkit Rahman Jamal, Mike Cerna, John Hanks
    NATIONAL Application Note 097 INSTRUMENTS® The Software is the Instrument ® Designing Filters Using the Digital Filter Design Toolkit Rahman Jamal, Mike Cerna, John Hanks Introduction The importance of digital filters is well established. Digital filters, and more generally digital signal processing algorithms, are classified as discrete-time systems. They are commonly implemented on a general purpose computer or on a dedicated digital signal processing (DSP) chip. Due to their well-known advantages, digital filters are often replacing classical analog filters. In this application note, we introduce a new digital filter design and analysis tool implemented in LabVIEW with which developers can graphically design classical IIR and FIR filters, interactively review filter responses, and save filter coefficients. In addition, real-world filter testing can be performed within the digital filter design application using a plug-in data acquisition board. Digital Filter Design Process Digital filters are used in a wide variety of signal processing applications, such as spectrum analysis, digital image processing, and pattern recognition. Digital filters eliminate a number of problems associated with their classical analog counterparts and thus are preferably used in place of analog filters. Digital filters belong to the class of discrete-time LTI (linear time invariant) systems, which are characterized by the properties of causality, recursibility, and stability. They can be characterized in the time domain by their unit-impulse response, and in the transform domain by their transfer function. Obviously, the unit-impulse response sequence of a causal LTI system could be of either finite or infinite duration and this property determines their classification into either finite impulse response (FIR) or infinite impulse response (IIR) system.
    [Show full text]
  • Control Theory
    Control theory S. Simrock DESY, Hamburg, Germany Abstract In engineering and mathematics, control theory deals with the behaviour of dynamical systems. The desired output of a system is called the reference. When one or more output variables of a system need to follow a certain ref- erence over time, a controller manipulates the inputs to a system to obtain the desired effect on the output of the system. Rapid advances in digital system technology have radically altered the control design options. It has become routinely practicable to design very complicated digital controllers and to carry out the extensive calculations required for their design. These advances in im- plementation and design capability can be obtained at low cost because of the widespread availability of inexpensive and powerful digital processing plat- forms and high-speed analog IO devices. 1 Introduction The emphasis of this tutorial on control theory is on the design of digital controls to achieve good dy- namic response and small errors while using signals that are sampled in time and quantized in amplitude. Both transform (classical control) and state-space (modern control) methods are described and applied to illustrative examples. The transform methods emphasized are the root-locus method of Evans and fre- quency response. The state-space methods developed are the technique of pole assignment augmented by an estimator (observer) and optimal quadratic-loss control. The optimal control problems use the steady-state constant gain solution. Other topics covered are system identification and non-linear control. System identification is a general term to describe mathematical tools and algorithms that build dynamical models from measured data.
    [Show full text]
  • Finite Impulse Response (FIR) Digital Filters (II) Ideal Impulse Response Design Examples Yogananda Isukapalli
    Finite Impulse Response (FIR) Digital Filters (II) Ideal Impulse Response Design Examples Yogananda Isukapalli 1 • FIR Filter Design Problem Given H(z) or H(ejw), find filter coefficients {b0, b1, b2, ….. bN-1} which are equal to {h0, h1, h2, ….hN-1} in the case of FIR filters. 1 z-1 z-1 z-1 z-1 x[n] h0 h1 h2 h3 hN-2 hN-1 1 1 1 1 1 y[n] Consider a general (infinite impulse response) definition: ¥ H (z) = å h[n] z-n n=-¥ 2 From complex variable theory, the inverse transform is: 1 n -1 h[n] = ò H (z)z dz 2pj C Where C is a counterclockwise closed contour in the region of convergence of H(z) and encircling the origin of the z-plane • Evaluating H(z) on the unit circle ( z = ejw ) : ¥ H (e jw ) = åh[n]e- jnw n=-¥ 1 p h[n] = ò H (e jw )e jnwdw where dz = jejw dw 2p -p 3 • Design of an ideal low pass FIR digital filter H(ejw) K -2p -p -wc 0 wc p 2p w Find ideal low pass impulse response {h[n]} 1 p h [n] = H (e jw )e jnwdw LP ò 2p -p 1 wc = Ke jnwdw 2p ò -wc Hence K h [n] = sin(nw ) n = 0, ±1, ±2, …. ±¥ LP np c 4 Let K = 1, wc = p/4, n = 0, ±1, …, ±10 The impulse response coefficients are n = 0, h[n] = 0.25 n = ±4, h[n] = 0 = ±1, = 0.225 = ±5, = -0.043 = ±2, = 0.159 = ±6, = -0.053 = ±3, = 0.075 = ±7, = -0.032 n = ±8, h[n] = 0 = ±9, = 0.025 = ±10, = 0.032 5 Non Causal FIR Impulse Response We can make it causal if we shift hLP[n] by 10 units to the right: K h [n] = sin((n -10)w ) LP (n -10)p c n = 0, 1, 2, ….
    [Show full text]
  • The Scientist and Engineer's Guide to Digital Signal Processing Properties of Convolution
    CHAPTER 7 Properties of Convolution A linear system's characteristics are completely specified by the system's impulse response, as governed by the mathematics of convolution. This is the basis of many signal processing techniques. For example: Digital filters are created by designing an appropriate impulse response. Enemy aircraft are detected with radar by analyzing a measured impulse response. Echo suppression in long distance telephone calls is accomplished by creating an impulse response that counteracts the impulse response of the reverberation. The list goes on and on. This chapter expands on the properties and usage of convolution in several areas. First, several common impulse responses are discussed. Second, methods are presented for dealing with cascade and parallel combinations of linear systems. Third, the technique of correlation is introduced. Fourth, a nasty problem with convolution is examined, the computation time can be unacceptably long using conventional algorithms and computers. Common Impulse Responses Delta Function The simplest impulse response is nothing more that a delta function, as shown in Fig. 7-1a. That is, an impulse on the input produces an identical impulse on the output. This means that all signals are passed through the system without change. Convolving any signal with a delta function results in exactly the same signal. Mathematically, this is written: EQUATION 7-1 The delta function is the identity for ( ' convolution. Any signal convolved with x[n] *[n] x[n] a delta function is left unchanged. This property makes the delta function the identity for convolution. This is analogous to zero being the identity for addition (a%0 ' a), and one being the identity for multiplication (a×1 ' a).
    [Show full text]
  • Exploring Anti-Aliasing Filters in Signal Conditioners for Mixed-Signal, Multimodal Sensor Conditioning by Arun T
    Texas Instruments Incorporated Amplifiers: Op Amps Exploring anti-aliasing filters in signal conditioners for mixed-signal, multimodal sensor conditioning By Arun T. Vemuri Systems Architect, Enhanced Industrial Introduction Figure 1. Multimodal, mixed-signal sensor-signal conditioner Some sensor-signal conditioners are used to process the output Analog Digital of multiple sense elements. This Domain Domain processing is often provided by multimodal, mixed-signal condi- tioners that can handle the out- puts from several sense elements Sense Digital Amplifier 1 ADC 1 at the same time. This article Element 1 Filter 1 analyzes the operation of anti- Processed aliasing filters in such sensor- Intelligent Output Compensation signal conditioners. Sense Digital Basics of sensor-signal Amplifier 2 ADC 2 Element 2 Filter 2 conditioners Sense elements, or transducers, convert a physical quantity of interest into electrical signals. Examples include piezo resistive bridges used to measure pres- sure, piezoelectric transducers used to detect ultrasonic than one sense element is processed by the same signal waves, and electrochemical cells used to measure gas conditioner is called multimodal signal conditioning. concentrations. The electrical signals produced by sense Mixed-signal signal conditioning elements are small and exhibit nonidealities, such as tem- Another aspect of sensor-signal conditioning is the electri- perature drifts and nonlinear transfer functions. cal domain in which the signal conditioning occurs. TI’s Sensor analog front ends such as the Texas Instruments PGA309 is an example of a device where the signal condi- (TI) LMP91000 and sensor-signal conditioners such as TI’s tioning of resistive-bridge sense elements occurs in the PGA400/450 are used to amplify the small signals produced analog domain.
    [Show full text]
  • 2913 Public Disclosure Authorized
    WPS A 13 POLICY RESEARCH WORKING PAPER 2913 Public Disclosure Authorized Financial Development and Dynamic Investment Behavior Public Disclosure Authorized Evidence from Panel Vector Autoregression Inessa Love Lea Zicchino Public Disclosure Authorized The World Bank Public Disclosure Authorized Development Research Group Finance October 2002 POLIcy RESEARCH WORKING PAPER 2913 Abstract Love and Zicchino apply vector autoregression to firm- availability of internal finance) that influence the level of level panel data from 36 countries to study the dynamic investment. The authors find that the impact of the relationship between firms' financial conditions and financial factors on investment, which they interpret as investment. They argue that by using orthogonalized evidence of financing constraints, is significantly larger in impulse-response functions they are able to separate the countries with less developed financial systems. The "fundamental factors" (such as marginal profitability of finding emphasizes the role of financial development in investment) from the "financial factors" (such as improving capital allocation and growth. This paper-a product of Finance, Development Research Group-is part of a larger effort in the group to study access to finance. Copies of the paper are available free from the World Bank, 1818 H Street NW, Washington, DC 20433. Please contact Kari Labrie, room MC3-456, telephone 202-473-1001, fax 202-522-1155, email address [email protected]. Policy Research Working Papers are also posted on the Web at http://econ.worldbank.org. The authors may be contacted at [email protected] or [email protected]. October 2002. (32 pages) The Policy Research Working Paper Series disseminates the findmygs of work mn progress to encouirage the excbange of ideas about development issues.
    [Show full text]
  • Transformations for FIR and IIR Filters' Design
    S S symmetry Article Transformations for FIR and IIR Filters’ Design V. N. Stavrou 1,*, I. G. Tsoulos 2 and Nikos E. Mastorakis 1,3 1 Hellenic Naval Academy, Department of Computer Science, Military Institutions of University Education, 18539 Piraeus, Greece 2 Department of Informatics and Telecommunications, University of Ioannina, 47150 Kostaki Artas, Greece; [email protected] 3 Department of Industrial Engineering, Technical University of Sofia, Bulevard Sveti Kliment Ohridski 8, 1000 Sofia, Bulgaria; mastor@tu-sofia.bg * Correspondence: [email protected] Abstract: In this paper, the transfer functions related to one-dimensional (1-D) and two-dimensional (2-D) filters have been theoretically and numerically investigated. The finite impulse response (FIR), as well as the infinite impulse response (IIR) are the main 2-D filters which have been investigated. More specifically, methods like the Windows method, the bilinear transformation method, the design of 2-D filters from appropriate 1-D functions and the design of 2-D filters using optimization techniques have been presented. Keywords: FIR filters; IIR filters; recursive filters; non-recursive filters; digital filters; constrained optimization; transfer functions Citation: Stavrou, V.N.; Tsoulos, I.G.; 1. Introduction Mastorakis, N.E. Transformations for There are two types of digital filters: the Finite Impulse Response (FIR) filters or Non- FIR and IIR Filters’ Design. Symmetry Recursive filters and the Infinite Impulse Response (IIR) filters or Recursive filters [1–5]. 2021, 13, 533. https://doi.org/ In the non-recursive filter structures the output depends only on the input, and in the 10.3390/sym13040533 recursive filter structures the output depends both on the input and on the previous outputs.
    [Show full text]
  • ELEG 5173L Digital Signal Processing Ch. 5 Digital Filters
    Department of Electrical Engineering University of Arkansas ELEG 5173L Digital Signal Processing Ch. 5 Digital Filters Dr. Jingxian Wu [email protected] 2 OUTLINE • FIR and IIR Filters • Filter Structures • Analog Filters • FIR Filter Design • IIR Filter Design 3 FIR V.S. IIR • LTI discrete-time system – Difference equation in time domain N M y(n) ak y(n k) bk x(n k) k 1 k 0 – Transfer function in z-domain N M k k Y (z) akY (z)z bk X (z)z k 1 k 0 M k bk z Y (z) k 0 H (z) N X (z) k 1 ak z k 1 4 FIR V.S. IIR • Finite impulse response (FIR) – difference equation in the time domain M y(n) bk x(n k) k 0 – Transfer function in the Z-domain M Y (z) k H (z) bk z X (z) k 0 – Impulse response h(n) [b ,b ,,b ] 0 1 M • The impulse response is of finite length finite impulse response 5 FIR V.S. IIR • Infinite impulse response (IIR) – Difference equation in the time domain N M y(n) ak y(n k) bk x(n k) k1 k0 – Transfer function in the z-domain M k bk z Y(z) k0 H (z) N X (z) k 1 ak z k1 – Impulse response can be obtained through inverse-z transform, and it has infinite length 6 FIR V.S. IIR • Example – Find the impulse response of the following system. Is it a FIR or IIR filter? Is it stable? 1 y(n) y(n 2) x(n) 4 7 FIR V.S.
    [Show full text]
  • Time Alignment)
    Measurement for Live Sound Welcome! Instructor: Jamie Anderson 2008 – Present: Rational Acoustics LLC Founding Partner & Systweak 1999 – 2008: SIA Software / EAW / LOUD Technologies Product Manager 1997 – 1999: Independent Sound Engineer A-1 Audio, Meyer Sound, Solstice, UltraSound / Promedia 1992 – 1997: Meyer Sound Laboratories SIM & Technical Support Manager 1991 – 1992: USC – Theatre Dept Education MFA: Yale School of Drama BS EE/Physics: Worcester Polytechnic University Instructor: Jamie Anderson Jamie Anderson [email protected] Rational Acoustics LLC Who is Rational Acoustics LLC ? 241 H Church St Jamie @RationalAcoustics.com Putnam, CT 06260 Adam @RationalAcoustics.com (860)928-7828 Calvert @RationalAcoustics.com www.RationalAcoustics.com Karen @RationalAcoustics.com and Barb @RationalAcoustics.com SmaartPIC @RationalAcoustics.com Support @RationalAcoustics.com Training @Rationalacoustics.com Info @RationalAcoustics.com What Are Our Goals for This Session? Understanding how our analyzers work – and how we can use them as a tool • Provide system engineering context (“Key Concepts”) • Basic measurement theory – Platform Agnostic Single Channel vs. Dual Channel Measurements Time Domain vs. Frequency Domain Using an analyzer is about asking questions . your questions Who Are You?" What Are Your Goals Today?" Smaart Basic ground rules • Class is informal - Get comfortable • Ask questions (Win valuable prizes!) • Stay awake • Be Courteous - Don’t distract! TURN THE CELL PHONES OFF NO SURFING / TEXTING / TWEETING PLEASE! Continuing
    [Show full text]
  • Digital Signal Processing Filter Design
    2065-27 Advanced Training Course on FPGA Design and VHDL for Hardware Simulation and Synthesis 26 October - 20 November, 2009 Digital Signal Processing Filter Design Massimiliano Nolich DEEI Facolta' di Ingegneria Universita' degli Studi di Trieste via Valerio, 10, 34127 Trieste Italy Filter design 1 Design considerations: a framework |H(f)| 1 C ıp 1 ıp ıs f 0 fp fs Passband Transition Stopband band The design of a digital filter involves five steps: Specification: The characteristics of the filter often have to be specified in the frequency domain. For example, for frequency selective filters (lowpass, highpass, bandpass, etc.) the specification usually involves tolerance limits as shown above. Coefficient calculation: Approximation methods have to be used to calculate the values hŒk for a FIR implementation, or ak, bk for an IIR implementation. Equivalently, this involves finding a filter which has H.z/ satisfying the requirements. Realisation: This involves converting H.z/ into a suitable filter structure. Block or flow diagrams are often used to depict filter structures, and show the computational procedure for implementing the digital filter. 1 Analysis of finite wordlength effects: In practice one should check that the quantisation used in the implementation does not degrade the performance of the filter to a point where it is unusable. Implementation: The filter is implemented in software or hardware. The criteria for selecting the implementation method involve issues such as real-time performance, complexity, processing requirements, and availability of equipment. 2 Finite impulse response (FIR) filter design A FIR filter is characterised by the equations N 1 yŒn D hŒkxŒn k kXD0 N 1 H.z/ D hŒkzk: kXD0 The following are useful properties of FIR filters: They are always stable — the system function contains no poles.
    [Show full text]