Stability Analysis of Continuous-Time Linear Time-Invariant Systems Kam Modjtahedzadeh Department of Physics, UC Santa Barbara

Total Page:16

File Type:pdf, Size:1020Kb

Stability Analysis of Continuous-Time Linear Time-Invariant Systems Kam Modjtahedzadeh Department of Physics, UC Santa Barbara Stability Analysis of Continuous-Time Linear Time-Invariant Systems Kam Modjtahedzadeh Department of Physics, UC Santa Barbara Abstract This paper focuses on the mathematical approaches to the analysis of stability that is a crucial step in the design of dynamical systems. Three methods are presented, namely, absolutely integrable impulse response, Fourier integral, and Laplace transform. The superiority of Laplace transform over the other methods becomes clear for several reasons that include the following: 1) It allows for the analysis of the stable, as well as, the unstable systems. 2) It not only determines absolute stability (a yes/no answer), but also shines light on the relative stability (how stable/unstable the system is), allowing for a design with a good degree of stability. 3) Its algebraic and convolution properties significantly simplify the mathematical manipulations involved in the analysis, especially when tackling a complex system composed of several simpler ones. A brief relevant introduction to the subject of systems is presented for the unfamiliar reader. Additionally, appropriate physical concepts and examples are presented for better clarity. Introduction later work of Kalman. Current research includes 퐻∞ (H infinity) and Hardy space [6]. However, Stability is a vast subject in the field of signals coverage of such advanced methods is beyond the and systems, in general, and, controls, in scope of this paper. particular. All controlled dynamical systems, being biological, chemical, physical, nuclear, or System Classifications of any other domain, must have a good degree of stability. The Three Mile Island and Chernobyl Strictly speaking, there are no linear-time catastrophes that occurred due to the process’s invariant (LTI) systems in reality. Elements of the temperatures going out of control are still fresh in systems are functions of the dependent, our memories. independent (e.g., time), or both variables. However, such effects become minimal when the Seldom, dynamical systems are intentionally systems operate near their equilibrium points, and made unstable. That adds an extra layer of temperature variation and other environmental difficulty in the design of the systems for effects are kept within a range. Thus the LTI acceptable overall stability [5]. Examples are models of the systems results in ordinary Segway, essentially an inverted pendulum for fun differential equations (ODE) with constant and space efficiency, and jet fighter planes, for coefficients. high maneuverability. This paper only concerns the continuous-time LTI Considering the degree of complexity of today’s systems. Although digital controllers have been systems, stability issue is a huge liability and a replacing the analog ones, the concepts headache for the control engineer; hence, the developed for continuous domain are extendable need for relevant and powerful mathematical to the discrete domain in the framework of z- tools. transform [3], [4]. A lot has been developed by control theoreticians Another classification of the systems is causality [4] such as earlier works of Bode, Nichol, and (causal versus non-causal). Almost all real Nyquist, building upon Laplace transform, and systems are causal as the output depends on the 1 1 1 past and present inputs but not on the future ℎ(푡) = 푒−푎푡 푢(푡), 푎 = inputs. However, time advance (positive time 푅 푅퐶 1 1 shift), such as preview, and anticipation is 퐻(푠) = 픏{ℎ(푡)} = ( ) , 퐼퐶 = 0 actually employed in practice and is non-causal. 푅 푠+푎 푄(푠) 1 1 Also, 퐻(푠) = = ( ) , 퐼퐶 = 0 Impulse Response and Transfer 푉(푠) 푅 푠+푎 Function An alternative name for the transfer function is As the name implies, the impulse response, ℎ(푡), system function. Note that the impulse response of an LTI system is defined as the output when and the transfer function characterize the system the input is a unit-impulse function, 푢(푡). and are neither input nor IC dependent; hence the Practical impulse resembling inputs include, name system function. force and voltage gendered by a hammer stroke A transfer function is a rational function when and fast on-off switching respectively. both the numerator and the denominator are The transfer function, 퐻(푠), of an LTI system is polynomials. The roots of the denominator and defined as the Laplace transform of the impulse the numerator are called poles and zeros response, with all the initial conditions (IC) set to respectively. The denominator polynomial set to zero. Alternatively, it can be defined as the zero forms the characteristic equation of the Laplace transform of the output over the Laplace system (the roots, or poles, determine the types of transform of the input, when 퐼퐶 = 0. functions in the partial fraction expansion). Laplace transform is covered in a dedicated Summarizing the system of figure 1, it is a causal section later in this paper. However, some of its continuous-time LTI system with a rational basic operations are employed in this section, in transfer function. Throughout the rest of this order to accomplish rudimentary system concepts paper, the assumptions of continuous-time linear all at once in the beginning. time-invariant will remain. Consider the RC circuit of figure 1 when a voltage source, 푣(푡), is placed in the loop, as the input, and the capacitor charge, 푞(푡), as the output. The parameters are constants for the range of operation and with respect to time. The voltage Fig. 2, Block Diagram of a Basic System source provides a continuous signal. The system is causal as it is a real practical circuit. Systems are often depicted by their corresponding block diagrams, showing their transfer functions, impulse responses, or even Fourier integrals in the blocks. For instance, figure 2 represents a simple system such as that of figure 1. The solution for the output , 푦(푡), of figure 2 via Fig. 1, RC Circuit impulse response is through the convolution The corresponding ODE and Laplace transform integral of the input, 푥(푡), and the impulse equations, when 퐼퐶 = 0, are: response, ℎ(푡): 푑 1 ( ) ∞ ( ) ( ) 푅 푞(푡) + 푞(푡) = 푣(푡) 푦 푡 = ∫−∞ 푥 푡 − 휏 ℎ 휏 푑휏 푑푡 퐶 ∞ 1 = 푥(휏)ℎ(푡 − 휏)푑휏 푅푆푄(푠) + 푄(푠) = 푉(푠) ∫−∞ 퐶 The integration becomes cumbersome for The impulse response and the transfer function complicated functions. A much more efficient are obtained according to their definitions: and systematic solution is by the Laplace 2 ∞ transform method, applying its convolution 푦(푡) = ∫ ℎ(휏)푥(푡 − 휏)푑휏 property: −∞ ∞ 푦(0) = ∫ ℎ(휏)푠푔푛(ℎ(휏))푑휏 푌(푠) = 퐻(푠)푋(푠) −∞ ∞ | | Inverse Laplace transform can be performed = ∫−∞ ℎ(휏) 푑휏 = ∞ subsequently with ease, with the aid of partial fraction expansion and an inverse table. As an example, the exponentially decaying However, the concept behind convolution is impulse response of the RC circuit is absolutely insightful, and also is employed in derivations integrable and stable. However, had the exponent of the impulse response been positive, then the and proofs. system would be neither absolutely integrable nor Absolutely Integrable Impulse stable. A negative resistor, i.e., physically a Response voltage source in phase with the input voltage source, would put energy into the system rather Bounded input bounded output (BIBO) stability than damping out energy, so causing the definition has been vastly accepted as the exponentially growing output. criterion for the stability of systems [3], [4]. It states that: a system is stable, if for every bounded Note that the absolutely integrable method only input the corresponding output is bounded. In reveals the absolute stability, but not the relative other words, the system is stable if the output is stability. finite for all possible finite inputs. Fourier integral For the particular case of continuous-time LTI Spectrum, 푋(푗휔), of a continuous-time systems, it can be proven that a system is (BIBO) aperiodic signal, 푥(푡), is obtained by the Fourier stable, if and only if, the impulse response ℎ(푡) is integral, also known as the Fourier transform absolutely integrable. of 푥(푡), if the integral converges: ∞ |ℎ(푡)|푑푡 < ∞ ∞ ∫−∞ ( ) { ( )} ( ) −푗휔푡 푋 푗휔 = ℱ 푥 푡 = ∫−∞ 푥 푡 푒 푑푡 The sufficiency part of the proof is obtained by The transformation is valid, provided that 푥(푡) manipulating the absolute value of the can be reconstructed accurately by an integral convolution integral: called the synthesis equation [1]: ∞ |푦(푡)| = |∫ ℎ(휏)푥(푡 − 휏)푑휏| 1 ∞ −∞ 푥(푡) = ∫ 푋(푗휔)푒푗휔푡푑휔 2휋 −∞ ∞ ≤ ∫ |ℎ(휏)푥(푡 − 휏)푑휏| −∞ The synthesis and the spectrum equations form ∞ the Fourier transform pair. A sufficient (but not ≤ ∫ |ℎ(휏)||푥(푡 − 휏)|푑휏 −∞ necessary) condition for the validity of the ∞ transformation is, if 푥(푡) is square integrable ≤ 퐵 ∫ |ℎ(휏)|푑휏 ≤ 퐵 푥 −∞ 푦 [1]: Where 퐵 and 퐵 are the input and output bounds ∞ 푥 푦 ∫ |푥(푡)|2푑푡 < ∞ respectively and are positive numbers. −∞ For the necessity part of the proof, it can be Alternatively, 푥(푡) should be absolutely shown that if the integral is not bounded, then integrable [1] (as one of the 3 Dirichlet there exists at least a bounded input that will drive conditions, where the other 2 are normally the output out of bound. satisfied for practical signals and systems). ∞ 1 ℎ(−푡) > 0 ∫−∞ |푥(푡)|푑푡 < ∞ 푥(푡) = 푠푔푛(ℎ(−푡)) = { 0 ℎ(−푡) = 0 −1 ℎ(−푡) < 0 3 The alternative condition also implies that any working in frequency domain (i.e., Bode, Nichol, stable system possesses Fourier transform of its and Nyquist methods not covered here) [4]. impulse response. Bilateral Laplace Transform The RC circuit of Figure 1 has an absolutely integrable, as well as a square integrable, impulse Bilateral Laplace transform of a continuous-time response. Thus, its Fourier transform exists and signal, 푥(푡), is defined as: is: { ( )} ∞ −푠푡 푋(푠) = 픏 푥 푡 = ∫−∞ 푥(푡)푒 푑푡 1 1 1 퐻(푗휔) = ( ) , 푎 = > 0 푅 푗휔+푎 푅퐶 It converts 푥(푡) into the complex function 푋(푠) of complex variable 푠 = 휎 + 푗휔, provided that The Fourier transform would not exist for if the the integral exists.
Recommended publications
  • Lecture 3: Transfer Function and Dynamic Response Analysis
    Transfer function approach Dynamic response Summary Basics of Automation and Control I Lecture 3: Transfer function and dynamic response analysis Paweł Malczyk Division of Theory of Machines and Robots Institute of Aeronautics and Applied Mechanics Faculty of Power and Aeronautical Engineering Warsaw University of Technology October 17, 2019 © Paweł Malczyk. Basics of Automation and Control I Lecture 3: Transfer function and dynamic response analysis 1 / 31 Transfer function approach Dynamic response Summary Outline 1 Transfer function approach 2 Dynamic response 3 Summary © Paweł Malczyk. Basics of Automation and Control I Lecture 3: Transfer function and dynamic response analysis 2 / 31 Transfer function approach Dynamic response Summary Transfer function approach 1 Transfer function approach SISO system Definition Poles and zeros Transfer function for multivariable system Properties 2 Dynamic response 3 Summary © Paweł Malczyk. Basics of Automation and Control I Lecture 3: Transfer function and dynamic response analysis 3 / 31 Transfer function approach Dynamic response Summary SISO system Fig. 1: Block diagram of a single input single output (SISO) system Consider the continuous, linear time-invariant (LTI) system defined by linear constant coefficient ordinary differential equation (LCCODE): dny dn−1y + − + ··· + _ + = an n an 1 n−1 a1y a0y dt dt (1) dmu dm−1u = b + b − + ··· + b u_ + b u m dtm m 1 dtm−1 1 0 initial conditions y(0), y_(0),..., y(n−1)(0), and u(0),..., u(m−1)(0) given, u(t) – input signal, y(t) – output signal, ai – real constants for i = 1, ··· , n, and bj – real constants for j = 1, ··· , m. How do I find the LCCODE (1)? .
    [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]
  • Step Response of Series RLC Circuit ‐ Output Taken Across Capacitor
    ESE 271 / Spring 2013 / Lecture 23 Step response of series RLC circuit ‐ output taken across capacitor. What happens during transient period from initial steady state to final steady state? 1 ESE 271 / Spring 2013 / Lecture 23 Transfer function of series RLC ‐ output taken across capacitor. Poles: Case 1: ‐‐two differen t real poles Case 2: ‐ two identical real poles ‐ complex conjugate poles Case 3: 2 ESE 271 / Spring 2013 / Lecture 23 Case 1: two different real poles. Step response of series RLC ‐ output taken across capacitor. Overdamped case –the circuit demonstrates relatively slow transient response. 3 ESE 271 / Spring 2013 / Lecture 23 Case 1: two different real poles. Freqqyuency response of series RLC ‐ output taken across capacitor. Uncorrected Bode Gain Plot Overdamped case –the circuit demonstrates relatively limited bandwidth 4 ESE 271 / Spring 2013 / Lecture 23 Case 2: two identical real poles. Step response of series RLC ‐ output taken across capacitor. Critically damped case –the circuit demonstrates the shortest possible rise time without overshoot. 5 ESE 271 / Spring 2013 / Lecture 23 Case 2: two identical real poles. Freqqyuency response of series RLC ‐ output taken across capacitor. Critically damped case –the circuit demonstrates the widest bandwidth without apparent resonance. Uncorrected Bode Gain Plot 6 ESE 271 / Spring 2013 / Lecture 23 Case 3: two complex poles. Step response of series RLC ‐ output taken across capacitor. Underdamped case – the circuit oscillates. 7 ESE 271 / Spring 2013 / Lecture 23 Case 3: two complex poles. Freqqyuency response of series RLC ‐ output taken across capacitor. Corrected Bode GiGain Plot Underdamped case –the circuit can demonstrate apparent resonant behavior.
    [Show full text]
  • Frequency Response and Bode Plots
    1 Frequency Response and Bode Plots 1.1 Preliminaries The steady-state sinusoidal frequency-response of a circuit is described by the phasor transfer function Hj( ) . A Bode plot is a graph of the magnitude (in dB) or phase of the transfer function versus frequency. Of course we can easily program the transfer function into a computer to make such plots, and for very complicated transfer functions this may be our only recourse. But in many cases the key features of the plot can be quickly sketched by hand using some simple rules that identify the impact of the poles and zeroes in shaping the frequency response. The advantage of this approach is the insight it provides on how the circuit elements influence the frequency response. This is especially important in the design of frequency-selective circuits. We will first consider how to generate Bode plots for simple poles, and then discuss how to handle the general second-order response. Before doing this, however, it may be helpful to review some properties of transfer functions, the decibel scale, and properties of the log function. Poles, Zeroes, and Stability The s-domain transfer function is always a rational polynomial function of the form Ns() smm as12 a s m asa Hs() K K mm12 10 (1.1) nn12 n Ds() s bsnn12 b s bsb 10 As we have seen already, the polynomials in the numerator and denominator are factored to find the poles and zeroes; these are the values of s that make the numerator or denominator zero. If we write the zeroes as zz123,, zetc., and similarly write the poles as pp123,, p , then Hs( ) can be written in factored form as ()()()s zsz sz Hs() K 12 m (1.2) ()()()s psp12 sp n 1 © Bob York 2009 2 Frequency Response and Bode Plots The pole and zero locations can be real or complex.
    [Show full text]
  • Linear Time Invariant Systems
    UNIT III LINEAR TIME INVARIANT CONTINUOUS TIME SYSTEMS CT systems – Linear Time invariant Systems – Basic properties of continuous time systems – Linearity, Causality, Time invariance, Stability – Frequency response of LTI systems – Analysis and characterization of LTI systems using Laplace transform – Computation of impulse response and transfer function using Laplace transform – Differential equation – Impulse response – Convolution integral and Frequency response. System A system may be defined as a set of elements or functional blocks which are connected together and produces an output in response to an input signal. The response or output of the system depends upon transfer function of the system. Mathematically, the functional relationship between input and output may be written as y(t)=f[x(t)] Types of system Like signals, systems may also be of two types as under: 1. Continuous-time system 2. Discrete time system Continuous time System Continuous time system may be defined as those systems in which the associated signals are also continuous. This means that input and output of continuous – time system are both continuous time signals. For example: Audio, video amplifiers, power supplies etc., are continuous time systems. Discrete time systems Discrete time system may be defined as a system in which the associated signals are also discrete time signals. This means that in a discrete time system, the input and output are both discrete time signals. For example, microprocessors, semiconductor memories, shift registers etc., are discrete time signals. LTI system:- Systems are broadly classified as continuous time systems and discrete time systems. Continuous time systems deal with continuous time signals and discrete time systems deal with discrete time system.
    [Show full text]
  • Simplified, Physically-Informed Models of Distortion and Overdrive Guitar Effects Pedals
    Proc. of the 10th Int. Conference on Digital Audio Effects (DAFx-07), Bordeaux, France, September 10-15, 2007 SIMPLIFIED, PHYSICALLY-INFORMED MODELS OF DISTORTION AND OVERDRIVE GUITAR EFFECTS PEDALS David T. Yeh, Jonathan S. Abel and Julius O. Smith Center for Computer Research in Music and Acoustics (CCRMA) Stanford University, Stanford, CA [dtyeh|abel|jos]@ccrma.stanford.edu ABSTRACT retained, however, because intermodulation due to mixing of sub- sonic components with audio frequency components is noticeable This paper explores a computationally efficient, physically in- in the audio band. formed approach to design algorithms for emulating guitar distor- tion circuits. Two iconic effects pedals are studied: the “Distor- Stages are partitioned at points in the circuit where an active tion” pedal and the “Tube Screamer” or “Overdrive” pedal. The element with low source impedance drives a high impedance load. primary distortion mechanism in both pedals is a diode clipper This approximation is also made with less accuracy where passive with an embedded low-pass filter, and is shown to follow a non- components feed into loads with higher impedance. Neglecting linear ordinary differential equation whose solution is computa- the interaction between the stages introduces magnitude error by a tionally expensive for real-time use. In the proposed method, a scalar factor and neglects higher order terms in the transfer func- simplified model, comprising the cascade of a conditioning filter, tion that are usually small in the audio band. memoryless nonlinearity and equalization filter, is chosen for its The nonlinearity may be evaluated as a nonlinear ordinary dif- computationally efficient, numerically robust properties.
    [Show full text]
  • Control System Design Methods
    Christiansen-Sec.19.qxd 06:08:2004 6:43 PM Page 19.1 The Electronics Engineers' Handbook, 5th Edition McGraw-Hill, Section 19, pp. 19.1-19.30, 2005. SECTION 19 CONTROL SYSTEMS Control is used to modify the behavior of a system so it behaves in a specific desirable way over time. For example, we may want the speed of a car on the highway to remain as close as possible to 60 miles per hour in spite of possible hills or adverse wind; or we may want an aircraft to follow a desired altitude, heading, and velocity profile independent of wind gusts; or we may want the temperature and pressure in a reactor vessel in a chemical process plant to be maintained at desired levels. All these are being accomplished today by control methods and the above are examples of what automatic control systems are designed to do, without human intervention. Control is used whenever quantities such as speed, altitude, temperature, or voltage must be made to behave in some desirable way over time. This section provides an introduction to control system design methods. P.A., Z.G. In This Section: CHAPTER 19.1 CONTROL SYSTEM DESIGN 19.3 INTRODUCTION 19.3 Proportional-Integral-Derivative Control 19.3 The Role of Control Theory 19.4 MATHEMATICAL DESCRIPTIONS 19.4 Linear Differential Equations 19.4 State Variable Descriptions 19.5 Transfer Functions 19.7 Frequency Response 19.9 ANALYSIS OF DYNAMICAL BEHAVIOR 19.10 System Response, Modes and Stability 19.10 Response of First and Second Order Systems 19.11 Transient Response Performance Specifications for a Second Order
    [Show full text]
  • Mathematical Modeling of Control Systems
    OGATA-CH02-013-062hr 7/14/09 1:51 PM Page 13 2 Mathematical Modeling of Control Systems 2–1 INTRODUCTION In studying control systems the reader must be able to model dynamic systems in math- ematical terms and analyze their dynamic characteristics.A mathematical model of a dy- namic system is defined as a set of equations that represents the dynamics of the system accurately, or at least fairly well. Note that a mathematical model is not unique to a given system.A system may be represented in many different ways and, therefore, may have many mathematical models, depending on one’s perspective. The dynamics of many systems, whether they are mechanical, electrical, thermal, economic, biological, and so on, may be described in terms of differential equations. Such differential equations may be obtained by using physical laws governing a partic- ular system—for example, Newton’s laws for mechanical systems and Kirchhoff’s laws for electrical systems. We must always keep in mind that deriving reasonable mathe- matical models is the most important part of the entire analysis of control systems. Throughout this book we assume that the principle of causality applies to the systems considered.This means that the current output of the system (the output at time t=0) depends on the past input (the input for t<0) but does not depend on the future input (the input for t>0). Mathematical Models. Mathematical models may assume many different forms. Depending on the particular system and the particular circumstances, one mathemati- cal model may be better suited than other models.
    [Show full text]
  • 12: Resonance
    12: Resonance • QuadraticFactors + • Damping Factor and Q • Parallel RLC • Behaviour at Resonance • Away from resonance • Bandwidth and Q • Power and Energy at Resonance + • Low Pass Filter • Resonance Peak for LP filter • Summary 12: Resonance E1.1 Analysis of Circuits (2017-10213) Resonance: 12 – 1 / 11 Quadratic Factors + 12: Resonance 2 • QuadraticFactors + A quadratic factor in a transfer function is: F (jω) = a (jω) + b (jω) + c. • Damping Factor and Q • Parallel RLC • Behaviour at Resonance • Away from resonance • Bandwidth and Q • Power and Energy at Resonance + • Low Pass Filter • Resonance Peak for LP filter • Summary E1.1 Analysis of Circuits (2017-10213) Resonance: 12 – 2 / 11 Quadratic Factors + 12: Resonance 2 • QuadraticFactors + A quadratic factor in a transfer function is: F (jω) = a (jω) + b (jω) + c. • Damping Factor and Q • Parallel RLC 2 • Behaviour at Resonance Case 1: If b 4ac then we can factorize it: • Away from resonance ≥ • Bandwidth and Q F (jω) = a(jω p1)(jω p2) • Power and Energy at − − Resonance + • Low Pass Filter • Resonance Peak for LP filter • Summary E1.1 Analysis of Circuits (2017-10213) Resonance: 12 – 2 / 11 Quadratic Factors + 12: Resonance 2 • QuadraticFactors + A quadratic factor in a transfer function is: F (jω) = a (jω) + b (jω) + c. • Damping Factor and Q • Parallel RLC 2 • Behaviour at Resonance Case 1: If b 4ac then we can factorize it: • Away from resonance ≥ • Bandwidth and Q F (jω) = a(jω p1)(jω p2) • Power and Energy at − − Resonance + b √b2 4ac • Low Pass Filter where p = − ± − . • Resonance Peak for LP i 2a filter • Summary E1.1 Analysis of Circuits (2017-10213) Resonance: 12 – 2 / 11 Quadratic Factors ++ 12: Resonance 2 • QuadraticFactors + A quadratic factor in a transfer function is: F (jω) = a (jω) + b (jω) + c.
    [Show full text]
  • Discussion 3B
    EECS 16B Designing Information Devices and Systems II Spring 2018 J. Roychowdhury and M. Maharbiz Discussion 3B 1 Transfer Functions When we analyzed circuits in the phasor domain, we always told you what the input voltage, or the input sinusoid to the circuit was. However, sometimes we have many input sinusoids, and we want to look at how a circuit (or system) generically responds to a sinusoid input of frequency w. We want to see how an input sinusoid “transfers” into an output sinusoid. How do we do this? Let’s start with a simple RC circuit. R i(t) + + vin C vout − − Figure 1: First Order RC Low Pass Filter 1 In the phasor domain, the impedance of the capacitor is ZC = jwC and the impedance of the resistor is ZR = R. Because we treat impedances the same as resistances, this circuit looks like a voltage divider in the phasor domain. Remember we must also represent vin as a phasor V˜in; transfer functions are in the phasor domain only, not the time domain. 1 Z 1 V˜ = C V˜ = jwC V˜ = V˜ out Z + Z in 1 in jwRC + 1 in R C R + jwC We define the frequency response as V˜ 1 H(w) = out = V˜in jwRC + 1 Now, given an arbitrary input sinusoid, if we multiply it by the frequency response, we can get the output sinusoid. What this allows us to do is to model any arbitrary circuit as a 4-port (two input, two output) black box. The transfer function completely defines how our circuit works.
    [Show full text]
  • Resonance of Fractional Transfer Functions of the Second Kind Rachid Malti, Xavier Moreau, Firas Khemane
    Resonance of fractional transfer functions of the second kind Rachid Malti, Xavier Moreau, Firas Khemane To cite this version: Rachid Malti, Xavier Moreau, Firas Khemane. Resonance of fractional transfer functions of the second kind. The 3th IFAC Workshop on Fractional Differentiation and its Applications, FDA08, Nov 2008, Ankara, Turkey. pp.1-6. hal-00326418 HAL Id: hal-00326418 https://hal.archives-ouvertes.fr/hal-00326418 Submitted on 26 Jan 2009 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Resonance of fractional transfer functions of the second kind Rachid MALTI, Xavier MOREAU, and Firas KHEMANE ∗ Bordeaux University – IMS, 351, cours de la Lib´eration, 33405 Talence Cedex, France. firstname.lastname @ims-bordeaux.fr { } Abstract: Canonical fractional transfer function of the second kind is studied in this paper. Stability and resonance conditions are determined in terms of pseudo-damping factor and commensurable order. Keywords: Resonance, fractional system, second order transfer function, canonical form. 1. INTRODUCTION F (s) + + 1 1 Σ K1 ν Σ K2 Commensurable fractional systems can be represented in s sν a transfer function form as: mB − − b sνj ν j T (s ) j=0 H (s)= = , (1) ν PmA R(s ) νi 1+ ais Fig.
    [Show full text]
  • The Z-Transform
    CHAPTER 33 The z-Transform Just as analog filters are designed using the Laplace transform, recursive digital filters are developed with a parallel technique called the z-transform. The overall strategy of these two transforms is the same: probe the impulse response with sinusoids and exponentials to find the system's poles and zeros. The Laplace transform deals with differential equations, the s-domain, and the s-plane. Correspondingly, the z-transform deals with difference equations, the z-domain, and the z-plane. However, the two techniques are not a mirror image of each other; the s-plane is arranged in a rectangular coordinate system, while the z-plane uses a polar format. Recursive digital filters are often designed by starting with one of the classic analog filters, such as the Butterworth, Chebyshev, or elliptic. A series of mathematical conversions are then used to obtain the desired digital filter. The z-transform provides the framework for this mathematics. The Chebyshev filter design program presented in Chapter 20 uses this approach, and is discussed in detail in this chapter. The Nature of the z-Domain To reinforce that the Laplace and z-transforms are parallel techniques, we will start with the Laplace transform and show how it can be changed into the z- transform. From the last chapter, the Laplace transform is defined by the relationship between the time domain and s-domain signals: 4 X(s) ' x(t) e &st dt m t '&4 where x(t) and X(s) are the time domain and s-domain representation of the signal, respectively.
    [Show full text]