Second Order Transfer Function Discrete Equations
|
|
- Brice Preston
- 6 years ago
- Views:
Transcription
1 Second Order Transfer Function Discrete Equations J. Riggs 23 Aug 2017
2 Transfer Function Equations pg 1 1 Introduction The objective of this paper is to develop the equations for a discrete implementation of a generic second order transfer function for use in digital simulation. The general form for the second order transfer function is given by where 2 Cs is the function output Rs is the function input. is the oscillation frequency when is zero is the damping ratio. As increases from zero, oscillations dampen exponentially. When 1 oscillation does not occur. The equation represents a system with the following differential equation: This second order system has two states. Let 2 Where is the function input. State might represent position, and therefore would be the first derivative of with respect to time, which is velocity. The original differential equation can now be written in terms of two linear equations; 2 The time history of the transfer function is obtained by solving the initial value problem represented by these two simultaneous differential equations.
3 Transfer Function Equations pg 2 2 Discrete Implementation The transfer function can be represented in state space, canonical form which can be solved (approximately) using numerical techniques. Take the general form for the linear function: The companion matrix, P, represents the system dynamics and input matrix, B, represents how the inputs couple into the system. Together, they define the change in system states. Propagating these changes over time step gives the value for the system states at time. The details of P and B are obtained from the system differential equations. Using the two linear differential equations, and the definitions of the system states, above, we can determine P and B for our generic second order system; companion matrix Input matrix Euler Integration The derivative can be approximated numerically (using the forward difference method) by Using this representation for the derivative in the state space equation gives This equation can be solved for the future value of, which is : (or) The above Equation gives the formula for the value of the state vector, x, for the next sample in the time series (updated at sample time)
4 Transfer Function Equations pg 3 Applying these substitutions and performing the algebraic expansion gives the following equations for the future value of the states; Δ Δ Runge Kutta 2 nd Order Integration For increased accuracy (compared to Euler integration), the transfer function equations can also be solved using a higher order method to approximate the function derivative, such as the second order Runge Kutta method. Since this is a multi pass method, the solution approach is a little different than the previous derivation. The implementation of 2 nd order Runge Kutta integration is performed as follows; Given; initial values for the two state variables x and x function input r(t), Equations for the state derivatives, x, x (i.e. x fx,x,t, x fx,x,t Apply 2 nd order Runge-Kutta integration to propagate the value of x and x forward in time. The algorithm is as follows; 1. Compute 1 x f x,x,r x f x,x,r 2. Compute x x Δt x x x Δt x tt Δt 3. Compute x f x,x,r x f x,x,r 4. Compute x x x x x x Estimated state values at time t Δt 1 Note that the time dependency of the function is contained within the function input, r. For example, r = k*t (ramp), r = sin(omega*t) (sine wave), etc.
5 Transfer Function Equations pg 4 Using the second order transfer function equations, define f x f ω x r 2ζω x The function equations can be expanded algebraically, in accordance with the 4 step solution algorithm, above, to compute the future value for the state variables. Step 1: x x x ω x r 2ζω x Step 2: x x Δt x x x Δt ω x r 2ζω x tt Δt Step 3: x x x ω x r 2ζω x x x Δt ω x r 2ζω x x ω x Δt x r 2ζω x Δt ω x r 2ζω x Step 4: x x Δt x x x Δt x Applying substitution and some algebraic manipulation, the updated states can be expressed in terms of the previous state values. These two equations give the future value at time of the two system states based on the current system state values at time t and the function input, r, using a second order Runge Kutta approximation to propagate the state derivative. x x Δt x ω x r 2 ζω x x x Δt 2 ω 2 Δt ζ ω 1x r x r 2 2Δt ζ ω 2Δt ζω Δt 2 ω x Note that in the above equation; r is the commanded value of the function at time t, which is the time of validity for x and x. r is the commanded value of the function at time Δ.
6 Transfer Function Equations pg 5 3 Performance Characterization This section illustrates the performance of the two numerical solutions to the second order transfer function using a side by side comparison of the two sets of equations. Figure 2 illustrates the discrete implementation of the transfer function equation for the forward difference (Euler) implementation, and Figure 2 shows the implementation using the second order Runge Kutta integration technique. The key difference between these two is that the second order method requires a memory state from the previous time step. For initialization, the two function states, x 1 and x 2 are assigned initial values. The function input, r(t), is the desired (i.e. commanded) value of the x 1 state, and the function outputs are the current state (x 1 ) and its first derivative (x 2 ). For the performance evaluation and comparison, we choose to use a nominal 2 Hz period (i.e. 4) and will propagate the solution of the transfer function for a unit step response, therefore x 1 initial value = 0 (initial position), and r(t) = 1 (unit step, from 0). We also set the initial value of x 2 (initial velocity) to zero. The performance of the two numerical solutions are compared by varying the damping ratio,, and the calculation time step, Δt. Figure 1 Transfer Function Discrete Implementation Euler Integration Figure 2 Transfer Function Discrete Implementation RK2 Integration
7 Transfer Function Equations pg Critically Damped Function The first comparison is for a function with critical damping ( 1.0) and the two function implementations will be compared at multiple calculation time steps. = 2 Hz (4 ) = 1.0 t = 0.05 (20 Hz) Figure 3 Unit Step response Comparison; Omega = 2 Hz, Zeta = 1, dt =.05 = 2 Hz (4 ) = 1.0 t = (40 Hz) Figure 4 Unit Step response Comparison; Omega = 2 Hz, Zeta = 1, dt =.025
8 Transfer Function Equations pg 7 = 2 Hz (4 ) = 1.0 t = 0.01 (100 Hz) Figure 5 Unit Step response Comparison; Omega = 2 Hz, Zeta = 1, dt =.01 = 2 Hz (4 ) = 1.0 t = (200 Hz) Figure 6 Unit Step response Comparison; Omega = 2 Hz, Zeta = 1, dt =.005
9 Transfer Function Equations pg 8 The preceding sequence of four charts illustrates the effect of increasing the calculation frequency (i.e. reducing the calculation time step). The first figure, above, shows the effect of a 20 Hz calculation frequency. Since the fundamental frequency of the function represented is 2 Hz, this represents a calculation rate of only 10 samples per period. There is a significant difference in the peak velocity between the two methods. The true peak is at a value of (at T= ), indicated by the short dashed line. The Euler method over estimates the peak velocity by more than 70%, while the Runge Kutta method under estimates the peak velocity by 27.7%. Comparing the first and second figures, above, shows the effect of doubling the calculation frequency from 20 Hz o 40 Hz. Notice that the first order method shows much greater change in the peak velocity as the time step is reduced, indicating that it is less accurate at the larger time step. The Euler method is now 20% high, and the Runge Kutta method is 5% low. Reducing the time step further improves the accuracy of both calculations. Table 1 Effect of Calculation Rate on Peak Velocity Error ( = 2Hz, = 1) Calculation Peak Velocity Error Freq (Hz) t Euler RK % 27.7% % 5.0% % 0.61% % 0.14% % 0.09% % 0.035% % 0.023% 1, % % 2, % % 4, % % 5, % % 6, % % 10, % % 20, % % Table 1 shows the improvement in the calculation as the calculation rate is increased. Clearly the second order Runge Kutta function improves at a much faster rate compared to the Euler integration implementation. For example, from increasing the calculation rate from 100 Hz to 10,000 Hz (two orders of magnitude) The Euler function improves by two orders of magnitude, from 6.8% to.063%), where the RK2 function improves by four orders of magnitude, from 0.6% to %. Another metric of comparison is to observe the calculation frequency required for a particular level of accuracy, say 0.1% or better. The RK2 method provides this level of accuracy at a calculation rate of 250 Hz, where the Euler method requires a rate of 6250 Hz, (25 times higher).
10 Transfer Function Equations pg Lightly damped Function Next we look at the two functions when the damping ratio is small. This is an area where the two functions will show some of the greatest differences. It must be borne in mind that the Euler implementation is a first order approximation to the function, therefore, it will eventually break down when trying to follow a second order function. In contrast, the 2 nd order Runge Kutta implementation will have a significant advantage in reproducing the second order transfer function. This is observed in the case of very low damping. To begin, we use a damping ratio of zero, resulting in a function that is in undamped oscillation at 2 Hz. Figure 7 shows the comparison of the two function implementations. Note that even at a calculation rate of 100 Hz, the first order Euler implementation is unstable, and shows a negatively damped characteristic. But the second order Runge Kutta implementation exhibits pure oscillation. By increasing the damping slightly to 0.063, the first order function becomes stable as seen in Figure 8. Note that in this figure, the first order function exhibits no damping under these conditions, while the second order function shows positive damping, more faithfully reproducing the expected result. Thus, a small increase in the damping ratio provides an improvement in the numerical stability of the Euler integration method = 2 Hz (4 ) = 0.0 t = 0.01 (100 Hz) Vel (EU) Vel (RK2) Pos (EU) Pos (RK2) Figure 7 Unit Step response Comparison; Omega = 2 Hz, Zeta = 0, dt =.001
11 Transfer Function Equations pg = 2 Hz (4 ) = t = 0.01 (100 Hz) Vel (EU) Vel (RK2) Pos (EU) Pos (RK2) Figure 8 Unit Step response Comparison; Omega = 2 Hz, Zeta = 0.063, dt =.001 Increasing the time step for the zero damping case will eventually make the second order implementation go unstable, as is illustrated in the following two figures. = 2 Hz (4 ) = 0.0 t = 0.02 (50 Hz) Figure 9 Unit Step response RK2 Implementation; Omega = 2 Hz, Zeta = 0.0, dt =.02
12 Transfer Function Equations pg 11 = 2 Hz (4 ) = 0.0 t = 0.04 (25Hz) Figure 10 Unit Step response RK2 Implementation; Omega = 2 Hz, Zeta = 0.0, dt =.04 Figure 9 shows that the second order Runge Kutta method is slightly unstable at a calculation rate of 50 Hz, as is evident by the slightly increasing amplitude over time in the velocity curve. Comparing Figure 10 to Figure 9 shows that reducing the calculation rate from 50 Hz to 25 Hz further increases the instability. Conclusions For the zero damping case, a calculation rate of about 100 Hz is necessary when using the second order Runge Kutta integrator, but the Euler integration method requires a significantly faster calculation rate (several orders of magnitude faster) in order to maintain numerical stability. Increasing the damping in the transfer function helps improve the stability of the Euler method tremendously.
13 Transfer Function Equations pg 12 4 Applying a Rate Limit It is herein assumed that the transfer function state x 1 represents a position and the user desires to apply a rate limit to the function output. The limit process is applied as follows; let be the rate limit. Compute an incremental limit for x 1 The position state, x, is expressed in terms of the previous value plus an increment; A limit of ± is placed on x x x x A limit of ± is placed on the value of x 4.1 Euler Method with Rate Limit Δ Δ 2 x Δ x,, nd Order Runge Kutta Method with Rate Limit x x Δt x ω x r 2 ζω x x x Δt 2 ω 2 Δt ζ ω 1x r x r 2Δt ζ ω 2Δt ζω Δt 2 ω x x Δt x ω x r 2 ζω x x,,
14 Transfer Function Equations pg 13 5 Observations and Conclusions Two numerical implementations of the generic second order transfer function have been developed and compared, one using a first order forward difference (Euler) method to propagate the state derivatives, and the second using a second order Runge Kutta algorithm. The RK2 based algorithm provides significant improvement in accuracy compared to the forward difference Euler method. Both methods break down in their ability to accurately represent the function under certain conditions, which depend on the calculation rate, the function natural frequency and damping ratio. Break down occurs at very small, and very large values of the damping ratio. Both numerical methods become unstable if the calculation time step is too large. The allowable time step magnitude is a function of the frequency and damping. For a damping ratio of 1, the largest allowable time step is approximately 1/(7 x ) ( in Hz), i.e. Δ 1 7 One interpretation of this is that we require at least 7 samples per cycle to maintain numeric stability in the calculation (with a damping ratio of 1). Reducing the damping ratio to less than 1 will require a smaller time step, with the Euler method requiring an even greater reduction in the time step compared to the Runge Kutta method. The Euler method is completely unsuitable in order to represent a function with zero damping. Because both methods have limitations, users should explore the function behavior under the conditions of their intended application in order to verify that the method will provide accurate results, especially if very high frequencies or low damping is desired.
6.1 Sketch the z-domain root locus and find the critical gain for the following systems K., the closed-loop characteristic equation is K + z 0.
6. Sketch the z-domain root locus and find the critical gain for the following systems K (i) Gz () z 4. (ii) Gz K () ( z+ 9. )( z 9. ) (iii) Gz () Kz ( z. )( z ) (iv) Gz () Kz ( + 9. ) ( z. )( z 8. ) (i)
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationSTABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential
More informationImplementation Issues for the Virtual Spring
Implementation Issues for the Virtual Spring J. S. Freudenberg EECS 461 Embedded Control Systems 1 Introduction One of the tasks in Lab 4 is to attach the haptic wheel to a virtual reference position with
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationMechatronics Assignment # 1
Problem # 1 Consider a closed-loop, rotary, speed-control system with a proportional controller K p, as shown below. The inertia of the rotor is J. The damping coefficient B in mechanical systems is usually
More informationProblem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013
EE 56: Digital Control Systems Problem Set 3: Solution Due on Mon 7 th Oct in class Fall 23 Problem For the causal LTI system described by the difference equation y k + 2 y k = x k, () (a) By first finding
More informationCourse Summary. The course cannot be summarized in one lecture.
Course Summary Unit 1: Introduction Unit 2: Modeling in the Frequency Domain Unit 3: Time Response Unit 4: Block Diagram Reduction Unit 5: Stability Unit 6: Steady-State Error Unit 7: Root Locus Techniques
More informationRunga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t)
Runga-Kutta Schemes Exact evolution over a small time step: x(t + t) =x(t)+ t 0 d(δt) F (x(t + δt),t+ δt) Expand both sides in a small time increment: x(t + t) =x(t)+x (t) t + 1 2 x (t)( t) 2 + 1 6 x (t)+
More informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationA Brief Introduction to Numerical Methods for Differential Equations
A Brief Introduction to Numerical Methods for Differential Equations January 10, 2011 This tutorial introduces some basic numerical computation techniques that are useful for the simulation and analysis
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationThursday, August 4, 2011
Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationEE402 - Discrete Time Systems Spring Lecture 10
EE402 - Discrete Time Systems Spring 208 Lecturer: Asst. Prof. M. Mert Ankarali Lecture 0.. Root Locus For continuous time systems the root locus diagram illustrates the location of roots/poles of a closed
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More information13 Numerical Solution of ODE s
13 NUMERICAL SOLUTION OF ODE S 28 13 Numerical Solution of ODE s In simulating dynamical systems, we frequently solve ordinary differential equations. These are of the form dx = f(t, x), dt where the function
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 8: Response Characteristics Overview In this Lecture, you will learn: Characteristics of the Response Stability Real
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. Åström and Murray, Exercise 1.3 2. Åström and Murray, Exercise 1.4 3. Åström and Murray, Exercise 2.6, parts (a) and (b) CDS 110a 1. Åström and Murray, Exercise 1.4 2. Åström and Murray, Exercise
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
More informationComputational project: Modelling a simple quadrupole mass spectrometer
Computational project: Modelling a simple quadrupole mass spectrometer Martin Duy Tat a, Anders Hagen Jarmund a a Norges Teknisk-Naturvitenskapelige Universitet, Trondheim, Norway. Abstract In this project
More informationThe Nonlinear Pendulum
The Nonlinear Pendulum - Pádraig Ó Conbhuí - 08531749 TP Monday 1. Abstract This experiment was performed to examine the effects that linearizing equations has on the accuracy of results and to find ways
More informationMEG6007: Advanced Dynamics -Principles and Computational Methods (Fall, 2017) Lecture DOF Modeling of Shock Absorbers. This lecture covers:
MEG6007: Advanced Dynamics -Principles and Computational Methods (Fall, 207) Lecture 4. 2-DOF Modeling of Shock Absorbers This lecture covers: Revisit 2-DOF Spring-Mass System Understand the den Hartog
More informationChapter 5 HW Solution
Chapter 5 HW Solution Review Questions. 1, 6. As usual, I think these are just a matter of text lookup. 1. Name the four components of a block diagram for a linear, time-invariant system. Let s see, I
More informationTime and Spatial Series and Transforms
Time and Spatial Series and Transforms Z- and Fourier transforms Gibbs' phenomenon Transforms and linear algebra Wavelet transforms Reading: Sheriff and Geldart, Chapter 15 Z-Transform Consider a digitized
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationDamped Oscillations *
OpenStax-CNX module: m58365 1 Damped Oscillations * OpenStax This wor is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 4.0 By the end of this section, you will be
More informationAN INTRODUCTION TO THE CONTROL THEORY
Open-Loop controller An Open-Loop (OL) controller is characterized by no direct connection between the output of the system and its input; therefore external disturbance, non-linear dynamics and parameter
More informationDynamic System Response. Dynamic System Response K. Craig 1
Dynamic System Response Dynamic System Response K. Craig 1 Dynamic System Response LTI Behavior vs. Non-LTI Behavior Solution of Linear, Constant-Coefficient, Ordinary Differential Equations Classical
More informationIn this Lecture. Frequency domain analysis
In this Lecture Frequency domain analysis Introduction In most cases we want to know the frequency content of our signal Why? Most popular analysis in frequency domain is based on work of Joseph Fourier
More informationTable of Laplacetransform
Appendix Table of Laplacetransform pairs 1(t) f(s) oct), unit impulse at t = 0 a, a constant or step of magnitude a at t = 0 a s t, a ramp function e- at, an exponential function s + a sin wt, a sine fun
More informationDiscrete Systems. Step response and pole locations. Mark Cannon. Hilary Term Lecture
Discrete Systems Mark Cannon Hilary Term 22 - Lecture 4 Step response and pole locations 4 - Review Definition of -transform: U() = Z{u k } = u k k k= Discrete transfer function: Y () U() = G() = Z{g k},
More informationApplication Note #3413
Application Note #3413 Manual Tuning Methods Tuning the controller seems to be a difficult task to some users; however, after getting familiar with the theories and tricks behind it, one might find the
More informationSecond Order Systems
Second Order Systems independent energy storage elements => Resonance: inertance & capacitance trade energy, kinetic to potential Example: Automobile Suspension x z vertical motions suspension spring shock
More informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationController Design using Root Locus
Chapter 4 Controller Design using Root Locus 4. PD Control Root locus is a useful tool to design different types of controllers. Below, we will illustrate the design of proportional derivative controllers
More informationProblem Weight Score Total 100
EE 350 EXAM IV 15 December 2010 Last Name (Print): First Name (Print): ID number (Last 4 digits): Section: DO NOT TURN THIS PAGE UNTIL YOU ARE TOLD TO DO SO Problem Weight Score 1 25 2 25 3 25 4 25 Total
More informationA Novel Approach for Complete Identification of Dynamic Fractional Order Systems Using Stochastic Optimization Algorithms and Fractional Calculus
5th International Conference on Electrical and Computer Engineering ICECE 2008, 20-22 December 2008, Dhaka, Bangladesh A Novel Approach for Complete Identification of Dynamic Fractional Order Systems Using
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 21 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University 1 All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationHomework Assignment 3
ECE382/ME482 Fall 2008 Homework 3 Solution October 20, 2008 1 Homework Assignment 3 Assigned September 30, 2008. Due in lecture October 7, 2008. Note that you must include all of your work to obtain full
More informationSAMPLE EXAMINATION PAPER (with numerical answers)
CID No: IMPERIAL COLLEGE LONDON Design Engineering MEng EXAMINATIONS For Internal Students of the Imperial College of Science, Technology and Medicine This paper is also taken for the relevant examination
More informationAPPLICATIONS FOR ROBOTICS
Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table
More informationCompensator Design to Improve Transient Performance Using Root Locus
1 Compensator Design to Improve Transient Performance Using Root Locus Prof. Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning
More information2 Background: Fourier Series Analysis and Synthesis
Signal Processing First Lab 15: Fourier Series Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in the Pre-Lab section before
More informationAdders, subtractors comparators, multipliers and other ALU elements
CSE4: Components and Design Techniques for Digital Systems Adders, subtractors comparators, multipliers and other ALU elements Adders 2 Circuit Delay Transistors have instrinsic resistance and capacitance
More informationThermodynamic Energy Equation
Thermodynamic Energy Equation The temperature tendency is = u T x v T y w T z + dt dt (1) where dt/dt is the individual derivative of temperature. This temperature change experienced by the air parcel
More informationChapter 7. Digital Control Systems
Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,
More informationFEEDBACK CONTROL SYSTEMS
FEEDBAC CONTROL SYSTEMS. Control System Design. Open and Closed-Loop Control Systems 3. Why Closed-Loop Control? 4. Case Study --- Speed Control of a DC Motor 5. Steady-State Errors in Unity Feedback Control
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationOrdinary differential equations. Phys 420/580 Lecture 8
Ordinary differential equations Phys 420/580 Lecture 8 Most physical laws are expressed as differential equations These come in three flavours: initial-value problems boundary-value problems eigenvalue
More informationMaple in Differential Equations
Maple in Differential Equations and Boundary Value Problems by H. Pleym Maple Worksheets Supplementing Edwards and Penney Differential Equations and Boundary Value Problems - Computing and Modeling Preface
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationPoles, Zeros and System Response
Time Response After the engineer obtains a mathematical representation of a subsystem, the subsystem is analyzed for its transient and steady state responses to see if these characteristics yield the desired
More informationChapter 7: Time Domain Analysis
Chapter 7: Time Domain Analysis Samantha Ramirez Preview Questions How do the system parameters affect the response? How are the parameters linked to the system poles or eigenvalues? How can Laplace transforms
More informationP441 Analytical Mechanics - I. RLC Circuits. c Alex R. Dzierba. In this note we discuss electrical oscillating circuits: undamped, damped and driven.
Lecture 10 Monday - September 19, 005 Written or last updated: September 19, 005 P441 Analytical Mechanics - I RLC Circuits c Alex R. Dzierba Introduction In this note we discuss electrical oscillating
More informationTransient response via gain adjustment. Consider a unity feedback system, where G(s) = 2. The closed loop transfer function is. s 2 + 2ζωs + ω 2 n
Design via frequency response Transient response via gain adjustment Consider a unity feedback system, where G(s) = ωn 2. The closed loop transfer function is s(s+2ζω n ) T(s) = ω 2 n s 2 + 2ζωs + ω 2
More information2.004 Dynamics and Control II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.004 Dynamics and Control II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts Institute
More informationOSE801 Engineering System Identification. Lecture 05: Fourier Analysis
OSE81 Engineering System Identification Lecture 5: Fourier Analysis What we will study in this lecture: A short introduction of Fourier analysis Sampling the data Applications Example 1 Fourier Analysis
More informationChapter 3 HW Solution
ME 48/58 Chapter 3 HW February 6, Chapter 3 HW Solution Problem. Here you re given a lead network often used in control systems to improve the transient response which adds around 6 of phase angle at about
More informationThe Nonlinear Pendulum
The Nonlinear Pendulum Evan Sheridan 11367741 Feburary 18th 013 Abstract Both the non-linear linear pendulum are investigated compared using the pendulum.c program that utilizes the trapezoid method for
More informationIndex. Index. More information. in this web service Cambridge University Press
A-type elements, 4 7, 18, 31, 168, 198, 202, 219, 220, 222, 225 A-type variables. See Across variable ac current, 172, 251 ac induction motor, 251 Acceleration rotational, 30 translational, 16 Accumulator,
More informationQuanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control. DC Motor Control Trainer (DCMCT) Student Manual
Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control DC Motor Control Trainer (DCMCT) Student Manual Table of Contents 1 Laboratory Objectives1 2 References1 3 DCMCT Plant
More informationNotes for ECE-320. Winter by R. Throne
Notes for ECE-3 Winter 4-5 by R. Throne Contents Table of Laplace Transforms 5 Laplace Transform Review 6. Poles and Zeros.................................... 6. Proper and Strictly Proper Transfer Functions...................
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Numerics Michael Bader SCCS Technical University of Munich Summer 018 Recall: Molecular Dynamics System of ODEs resulting force acting on a molecule: F i = j
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationAnalysis and Design of Control Systems in the Time Domain
Chapter 6 Analysis and Design of Control Systems in the Time Domain 6. Concepts of feedback control Given a system, we can classify it as an open loop or a closed loop depends on the usage of the feedback.
More informationCHAPTER 80 NUMERICAL METHODS FOR FIRST-ORDER DIFFERENTIAL EQUATIONS
CHAPTER 8 NUMERICAL METHODS FOR FIRST-ORDER DIFFERENTIAL EQUATIONS EXERCISE 33 Page 834. Use Euler s method to obtain a numerical solution of the differential equation d d 3, with the initial conditions
More informationThe (Fast) Fourier Transform
The (Fast) Fourier Transform The Fourier transform (FT) is the analog, for non-periodic functions, of the Fourier series for periodic functions can be considered as a Fourier series in the limit that the
More informationTransient Response of a Second-Order System
Transient Response of a Second-Order System ECEN 830 Spring 01 1. Introduction In connection with this experiment, you are selecting the gains in your feedback loop to obtain a well-behaved closed-loop
More informationLecture V : Oscillatory motion and spectral analysis
Lecture V : Oscillatory motion and spectral analysis I. IDEAL PENDULUM AND STABILITY ANALYSIS Let us remind ourselves of the equation of motion for the pendulum. Remembering that the external torque applied
More informationDate: _15 April (1) The only reference material you may use is one 8½x11 crib sheet and a calculator.
PH1140: Oscillations and Waves Name: SOLUTIONS Conference: Date: _15 April 2005 EXAM #2: D2005 INSTRUCTIONS: (1) The only reference material you may use is one 8½x11 crib sheet and a calculator. (2) Show
More informationDamped Oscillation Solution
Lecture 19 (Chapter 7): Energy Damping, s 1 OverDamped Oscillation Solution Damped Oscillation Solution The last case has β 2 ω 2 0 > 0. In this case we define another real frequency ω 2 = β 2 ω 2 0. In
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationINTRODUCTION TO DIGITAL CONTROL
ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall K(s +1)(s +2) G(s) =.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering. Dynamics and Control II Fall 7 Problem Set #7 Solution Posted: Friday, Nov., 7. Nise problem 5 from chapter 8, page 76. Answer:
More informationReinforcement Learning In Continuous Time and Space
Reinforcement Learning In Continuous Time and Space presentation of paper by Kenji Doya Leszek Rybicki lrybicki@mat.umk.pl 18.07.2008 Leszek Rybicki lrybicki@mat.umk.pl Reinforcement Learning In Continuous
More information(b) A unity feedback system is characterized by the transfer function. Design a suitable compensator to meet the following specifications:
1. (a) The open loop transfer function of a unity feedback control system is given by G(S) = K/S(1+0.1S)(1+S) (i) Determine the value of K so that the resonance peak M r of the system is equal to 1.4.
More informationChapter 2 PARAMETRIC OSCILLATOR
CHAPTER- Chapter PARAMETRIC OSCILLATOR.1 Introduction A simple pendulum consists of a mass m suspended from a string of length L which is fixed at a pivot P. When simple pendulum is displaced to an initial
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationu (t t ) + e ζωn (t tw )
LINEAR CIRCUITS LABORATORY OSCILLATIONS AND DAMPING EFFECT PART I TRANSIENT RESPONSE TO A SQUARE PULSE Transfer Function F(S) = ω n 2 S 2 + 2ζω n S + ω n 2 F(S) = S 2 + 3 RC ( RC) 2 S + 1 RC ( ) 2 where
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationVALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur 603 203. DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING SUBJECT QUESTION BANK : EC6405 CONTROL SYSTEM ENGINEERING SEM / YEAR: IV / II year
More informationControl Systems I Lecture 10: System Specifications
Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture
More informationNumerical Integration of Ordinary Differential Equations for Initial Value Problems
Numerical Integration of Ordinary Differential Equations for Initial Value Problems Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu These slides are a
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
More informationSpontaneous Speed Reversals in Stepper Motors
Spontaneous Speed Reversals in Stepper Motors Marc Bodson University of Utah Electrical & Computer Engineering 50 S Central Campus Dr Rm 3280 Salt Lake City, UT 84112, U.S.A. Jeffrey S. Sato & Stephen
More informationChecking the Radioactive Decay Euler Algorithm
Lecture 2: Checking Numerical Results Review of the first example: radioactive decay The radioactive decay equation dn/dt = N τ has a well known solution in terms of the initial number of nuclei present
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More information7.2 Controller tuning from specified characteristic polynomial
192 Finn Haugen: PID Control 7.2 Controller tuning from specified characteristic polynomial 7.2.1 Introduction The subsequent sections explain controller tuning based on specifications of the characteristic
More informationx(t+ δt) - x(t) = slope δt t+δt
Techniques of Physics Worksheet 2 Classical Vibrations and Waves Introduction You will have encountered many different examples of wave phenomena in your courses and should be familiar with most of the
More informationElectric Circuit Theory
Electric Circuit Theory Nam Ki Min nkmin@korea.ac.kr 010-9419-2320 Chapter 8 Natural and Step Responses of RLC Circuits Nam Ki Min nkmin@korea.ac.kr 010-9419-2320 8.1 Introduction to the Natural Response
More informationLaboratory handout 5 Mode shapes and resonance
laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful
More informationLecture 7:Time Response Pole-Zero Maps Influence of Poles and Zeros Higher Order Systems and Pole Dominance Criterion
Cleveland State University MCE441: Intr. Linear Control Lecture 7:Time Influence of Poles and Zeros Higher Order and Pole Criterion Prof. Richter 1 / 26 First-Order Specs: Step : Pole Real inputs contain
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationComputational Methods in Plasma Physics
Computational Methods in Plasma Physics Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Purpose of Talk Describe use of numerical methods to solve simple problem in plasma
More information