MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

Similar documents
Control System Engineering (EE301T) Assignment: 2

DYNAMICS and CONTROL

CSE 245: Computer Aided Circuit Simulation and Verification

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

Elementary Differential Equations and Boundary Value Problems

Lecture 4: Laplace Transforms

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

CONTROL SYSTEM DESIGN OF GUIDED MISSILE

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

C From Faraday's Law, the induced voltage is, C The effect of electromagnetic induction in the coil itself is called selfinduction.

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

PFC Predictive Functional Control

Microscopic Flow Characteristics Time Headway - Distribution

CPSC 211 Data Structures & Implementations (c) Texas A&M University [ 259] B-Trees

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

Midterm exam 2, April 7, 2009 (solutions)

Transfer function and the Laplace transformation

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

A Backstepping Simple Adaptive Control Application to Flexible Space Structures

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

Double Slits in Space and Time

EE 350 Signals and Systems Spring 2005 Sample Exam #2 - Solutions

Chapter 7 Stead St y- ate Errors

Decline Curves. Exponential decline (constant fractional decline) Harmonic decline, and Hyperbolic decline.

Continous system: differential equations

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

Ratio-Product Type Exponential Estimator For Estimating Finite Population Mean Using Information On Auxiliary Attribute

Lecture 2: Current in RC circuit D.K.Pandey

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

First Lecture of Machine Learning. Hung-yi Lee

5. An object moving along an x-coordinate axis with its scale measured in meters has a velocity of 6t

whereby we can express the phase by any one of the formulas cos ( 3 whereby we can express the phase by any one of the formulas

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

Problem 2. Describe the following signals in terms of elementary functions (δ, u,r, ) and compute. x(t+2) x(2-t) RT_1[x] -3-2 = 1 2 = 1

Chemistry 988 Part 1

UNSTEADY FLOW OF A FLUID PARTICLE SUSPENSION BETWEEN TWO PARALLEL PLATES SUDDENLY SET IN MOTION WITH SAME SPEED

Math 3301 Homework Set 6 Solutions 10 Points. = +. The guess for the particular P ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos 2 t : 4D= 2

Supplementary Materials

3.4 Repeated Roots; Reduction of Order

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

The transition:transversion rate ratio vs. the T-ratio.

LaPlace Transform in Circuit Analysis

( ) ( ) + = ( ) + ( )

I) Title: Rational Expectations and Adaptive Learning. II) Contents: Introduction to Adaptive Learning

On the Speed of Heat Wave. Mihály Makai

Chapter 9 Transient Response

Linear Systems Analysis in the Time Domain

Transient Performance Analysis of Serial Production Lines

Ministry of Education and Science of Ukraine National Technical University Ukraine "Igor Sikorsky Kiev Polytechnic Institute"

Midterm Examination (100 pts)

Wave Equation (2 Week)

General Article Application of differential equation in L-R and C-R circuit analysis by classical method. Abstract

Revisiting what you have learned in Advanced Mathematical Analysis

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

A Method of Performance Assessment of PID Controller with Actuator Saturation

The Variance-Covariance Matrix

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

Chapter 6. PID Control

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Linear Systems

REPETITION before the exam PART 2, Transform Methods. Laplace transforms: τ dτ. L1. Derive the formulas : L2. Find the Laplace transform F(s) if.

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

where: u: input y: output x: state vector A, B, C, D are const matrices

Coherence and interactions in diffusive systems. Cours 4. Diffusion + e-e interations

PS#4 due today (in class or before 3pm, Rm ) cannot ignore the spatial dimensions

Lie Groups HW7. Wang Shuai. November 2015

EXERCISE - 01 CHECK YOUR GRASP

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional

Failure Load of Plane Steel Frames Using the Yield Surface Method

CHE302 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES. Professor Dae Ryook Yang

Problem Set #3: AK models

Economics 302 (Sec. 001) Intermediate Macroeconomic Theory and Policy (Spring 2011) 3/28/2012. UW Madison

symmetric/hermitian matrices, and similarity transformations

The Science of Monetary Policy

Copyright 2012 Pearson Education, Inc. Publishing as Prentice Hall.

XV Exponential and Logarithmic Functions

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

Phys463.nb Conductivity. Another equivalent definition of the Fermi velocity is

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

u r du = ur+1 r + 1 du = ln u + C u sin u du = cos u + C cos u du = sin u + C sec u tan u du = sec u + C e u du = e u + C

Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

( ) 2! l p. Nonlinear Dynamics for Gear Fault Level. ( ) f ( x) ( ),! = sgn % " p. Open Access. Su Xunwen *,1, Liu Jinhao 1 and Wang Shaoping 2. !

Chapter 17 Handout: Autocorrelation (Serial Correlation)

Consider a system of 2 simultaneous first order linear equations

Chapter 8 The Complete Response of RL and RC Circuits

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

Chapter 12 Introduction To The Laplace Transform

ERROR ANALYSIS A.J. Pintar and D. Caspary Department of Chemical Engineering Michigan Technological University Houghton, MI September, 2012

SOLUTIONS. 1. Consider two continuous random variables X and Y with joint p.d.f. f ( x, y ) = = = 15. Stepanov Dalpiaz

A Condition for Stability in an SIR Age Structured Disease Model with Decreasing Survival Rate

Sundials and Linear Algebra

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 3: Ideal Nozzle Fluid Mechanics

CHAPTER 9 Compressible Flow

H is equal to the surface current J S

Transcription:

MEM 355 Prformanc Enhancmn of Dynamical Sysms A Firs Conrol Problm - Cruis Conrol Harry G. Kwany Darmn of Mchanical Enginring & Mchanics Drxl Univrsiy

Cruis Conrol ( ) mv = F mg sinθ cv v +.2v= u 9.8θ v[ m/ s] sd ( m/s=36 km/h=22 mils/hr) u normalizd hrol u θ [ rad] roadway slo sd command sd rror - conrol driving forc ngin disurbing forc - vhicl sd

Criria Wha characrisics do w xc in a good cruis conrol sysm? Ulima racking rror in rsons o sd command Ulima racking rror in rsons o disurbanc Transin rsons im Transin ovrshoo, undrshoo Robusnss olranc o unmodld dynamics or aramr variaion

Cruis Conrol: driv rror rsons from diffrnial quaions Assumions: 'roorional' lus 'ingral' conrol, ignor ngin dynamics ( ) i ( ( ) ( )) ( ) ( ) ( ) u = k v v + k v τ v τ dτ ( ) = ( ) ( ) Dfin: v v v +.2v= u 9.8θ +.2 = u + 9.8 θ + v +.2v ( ) ( ( τ) ) ( ) ( ) u = k v v + k v v dτ u = k + k i i + (.2 + k) + k i = 9.8 θ + v +.2v ( s + ).2 9.8s E( s) = V ( s) + Θ s s +.2 + k s+ k s +.2 + k s+ k ( ) ( ) 2 2 i i Error rsons o command ( ) sd rror 2 Inus: v θ sd command road disurbanc Error rsons o disurbanc

Cruis Conrol: driv rror rsons using block diagrams Engin Comnsaor θ 9.8 ( ) Vhicl v Rmmbr h formula Or do h algbra ks+ k s i s +.2 ( s +.2) v : Gv = G = ks+ k s + ( + k ) s+ k + s s+.2 ks + ki θ : E = ( 9.8Θ) E s+.2 s ks + ki 9.8 9.8s + E = Θ Gθ = 2 s s+.2 s +.2 + k s+ ki v 2 i.2 i ( ) v Closd loo ransfr funcions wih wo dsign aramrs

Cruis Conrol 4.8.6.4 Rsons o uni s disurbanc 9.8s E( s) = Θ( s), Θ ( s) = s (.2 ) 2 s + + k s+ ki.2 2 4 6 8 2 4 ( s + ).2 E( s) = V ( s), V ( s) = s (.2 ) 2 s + + k s+ ki.75.5.25 Rsons o uni s command k =, k =,.5,, 2, 6 i -.25 -.5 2 4 6 8 2 4

Obsrvaions Boh ransfr funcions hav sam dnominaor (sam ols), bu diffrn numraors (diffrn zros) Whn = (roorional conrol) h ulima rror is no zro, k i in fac h ulima rror in disurbanc is larg. rsons o command is vry small, bu o For sabiliy w can look a ihr ransfr funcion, bu for rformanc w nd o considr boh. To valua k, k i is hlful o mak h associaion i s + (.2 + k ) s+ k s + 2ρω s+ ω 2 2 2 i

Rfin h Conrol Noic ha w can scify, ρω, and choos k i = ω, k = 2ρω.2 2 L us fix ρ =.77 and look a ω =, 2, 2,3.4.3.2..8.6.4.2 Disurbanc rsons 2 4 6 8 2 4 Command rsons -.2 2 4 6 8 2 4

Effc of Engin Dynamics G ng Error rsons o command s using h fass conrollr. ρ =.77, ω = 3 Engin im consan. sc, and.2 sc. 2 =, τ =.,.2 τ s +.75.5.25 -.25 -.5 -.75 6 4 2-2 -4-6 τ =.sc 2 4 6 8 2 4 τ =.2 sc 2 4 6 8 2 4 Dramaic rducion in daming raion Unsabl

Effc of Engin Dynamics ~ 2 Suos w us h slows conrollr ρ =.77, ω = Hr w s how h rsons dgrads whn slow ngin is includd, a las i is sill sabl Pushing for high rformanc ofn lads o non-robus dsign..6 Wih slow ngin.5.4.3.2. -. Wih idal ngin 2 4 6 8 2 4 Error rsons o disurbanc

Add Fdforward v Gc ks+ k s i G θ ( ) 9.8 Us il snsor (acclromrs) o sima roadway slo. Gv s +.2 v.6.4 E ( G ) Gv = + GGG v c D.2 Slow conrollr, slow ngin 2 4 6 8 2 4 -.2 Wih fdforward

Summary Any on of h closd loo ransfr funcions can b usd for sabiliy analysis (all hav sam ols) Prformanc analysis usually rquirs considring wo or mor closd loo ransfr funcions. Ulima rror dnds on conrollr y,.g. PI conrollr rsuld in zro rror vnually, bu P conrollr lf som rsidual rror nonrivial in h cas of disurbanc. W can choos conrol aramrs o sha ransin rsons (loca closd loo ols) in his scial cas w usd our knowldg of 2 nd ordr sysm bhavior. Th sysm may b snsiiv o modl accuracy, including nglcd dynamics vn o h oin of insabiliy.

Nx Ss Conrolling h ulima rror Shaing h ransin rsons: closd loo ol locaion via roo locus Evaluaing closd loo sysm sabiliy robusnss h abiliy o rmain sabl whn h modl is uncrain Conrollr dsign o achiv robus rformanc Ohr (mor dirc) mhods for shaing ransin rsons