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

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

Elementary Differential Equations and Boundary Value Problems

Transfer function and the Laplace transformation

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

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

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

Chapter 12 Introduction To The Laplace Transform

Midterm exam 2, April 7, 2009 (solutions)

Charging of capacitor through inductor and resistor

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!

On the Speed of Heat Wave. Mihály Makai

CSE 245: Computer Aided Circuit Simulation and Verification

fiziks Institute for NET/JRF, GATE, IIT JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES MATEMATICAL PHYSICS SOLUTIONS are

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

Lecture 4: Laplace Transforms

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

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

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

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

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

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

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

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

Wave Equation (2 Week)

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

LaPlace Transform in Circuit Analysis

H is equal to the surface current J S

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

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

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

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

EXERCISE - 01 CHECK YOUR GRASP

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

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

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

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

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

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

Chap.3 Laplace Transform

Review Lecture 5. The source-free R-C/R-L circuit Step response of an RC/RL circuit. The time constant = RC The final capacitor voltage v( )

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

CHAPTER. Linear Systems of Differential Equations. 6.1 Theory of Linear DE Systems. ! Nullcline Sketching. Equilibrium (unstable) at (0, 0)

3(8 ) (8 x x ) 3x x (8 )

Control System Engineering (EE301T) Assignment: 2

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

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

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

Institute of Actuaries of India

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

Microscopic Flow Characteristics Time Headway - Distribution

EE 434 Lecture 22. Bipolar Device Models

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

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:

DE Dr. M. Sakalli

Logistic equation of Human population growth (generalization to the case of reactive environment).

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!

Poisson process Markov process

Final Exam : Solutions

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

A Propagating Wave Packet Group Velocity Dispersion

Impulsive Differential Equations. by using the Euler Method

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

On General Solutions of First-Order Nonlinear Matrix and Scalar Ordinary Differential Equations

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

Why Laplace transforms?

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

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

Inverse Fourier Transform. Properties of Continuous time Fourier Transform. Review. Linearity. Reading Assignment Oppenheim Sec pp.289.

Lagrangian for RLC circuits using analogy with the classical mechanics concepts

Revisiting what you have learned in Advanced Mathematical Analysis

Chapter 6 Differential Equations and Mathematical Modeling

3.4 Repeated Roots; Reduction of Order

Continous system: differential equations

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

A MATHEMATICAL MODEL FOR NATURAL COOLING OF A CUP OF TEA

Circuit Transients time

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

Discussion 06 Solutions

46. Let y = ln r. Then dy = dr, and so. = [ sin (ln r) cos (ln r)

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Nikesh Bajaj. Fourier Analysis and Synthesis Tool. Guess.? Question??? History. Fourier Series. Fourier. Nikesh Bajaj

Physics 160 Lecture 3. R. Johnson April 6, 2015

Circuits and Systems I

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

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

( ) C R. υ in RC 1. cos. ,sin. ω ω υ + +

Double Slits in Space and Time

Lecture 26: Leapers and Creepers

Mundell-Fleming I: Setup

AN INTRODUCTION TO FOURIER ANALYSIS PROF. VEDAT TAVSANOĞLU

14.02 Principles of Macroeconomics Fall 2005 Quiz 3 Solutions

10. If p and q are the lengths of the perpendiculars from the origin on the tangent and the normal to the curve

Frequency Response. Lecture #12 Chapter 10. BME 310 Biomedical Computing - J.Schesser

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

Basic Polyhedral theory

University of Kansas, Department of Economics Economics 911: Applied Macroeconomics. Problem Set 2: Multivariate Time Series Analysis

The Science of Monetary Policy

Asymptotic Solutions of Fifth Order Critically Damped Nonlinear Systems with Pair Wise Equal Eigenvalues and another is Distinct

Chapter 1 Fundamental Concepts

Transcription:

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw. conx and dircion Imagin a sysm ha varis in im; w migh plo is oupu vs. im. A plo migh imply an quaion, and h quaion is usually an ODE (ordinary diffrnial quaion. Thrfor, w will rviw h mah of h firs-ordr ODE whil mphasizing how i can rprsn a namic sysm. W xamin how h sysm is affcd by is iniial condiion and by disurbancs, whr h disurbancs may b non-smooh, mulipl, or dlayd.. firs-ordr, linar, variabl-cofficin ODE Th dpndn variabl y( dpnds on is firs drivaiv and forcing funcion x(. Whn h indpndn variabl is, y is y. a ( + y( = Kx( y( = y (.- d In wriing (.- w hav arrangd a cofficin of + for y. Thrfor a( mus hav dimnsions of indpndn variabl, and K has dimnsions of y/x. W solv (.- by dfining h ingraing facor p( d p ( = xp (.- a( Noic ha p( is dimnsionlss, as is h quoin undr h ingral. Th soluion p( y( K p(x( y ( = + d (.-3 p( p( a( compriss conribuions from h iniial condiion y( and h forcing funcion Kx(. Ths ar known as h homognous (as if h righ-hand sid wr zro and paricular (dpnds on h righ-hand sid soluions. In h languag of namic sysms, w can hink of y( as h rspons of h sysm o inpu disurbancs Kx( and y(.. firs-ordr ODE, spcial cas for procss conrol applicaions Th indpndn variabl will rprsn im. For many procss conrol applicaions, a( in (.- will b a posiiv consan; w call i h im consan. + y( = Kx( y( = y (.- d Th ingraing facor (.- is rvisd 5 Jan

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw d p ( = xp = (.- and h soluion (.-3 bcoms ( K y( = y + x(d (.-3 Th iniial condiion affcs h sysm rspons from h bginning, bu is ffc dcays o zro according o h magniud of h im consan - largr im consans rprsn slowr dcay. If no furhr disurbd by som x(, h firs ordr sysm rachs quilibrium a zro. Howvr, mos pracical sysms ar disurbd. K is a propry of h sysm, calld h gain. By is magniud and sign, h gain influncs how srongly y rsponds o x. Th form of h rspons dpnds on h naur of h disurbanc. Exampl: suppos x is a uni sp funcion a im. Bfor w procd formally, l us hink inuiivly. From (.-3 w xpc h rspons y o dcay oward zro from IC y. A im, h sysm will rspond o bing hi wih a sp disurbanc. Afr a long im, hr will b no mmory of h iniial condiion, and h sysm will rspond only o h disurbanc inpu. Bcaus his is consan afr h sp, w guss ha h rspons will also bcom consan. Now h mah: from (.-3 y( = y = y ( K + ( ( + KU( U( d (.-4 Figur.- shows h soluion. Noic ha h paricular soluion maks no conribuion bfor im. Th iniial condiion dcays, and wih no disurbanc would coninu o zro. A, howvr, h sysm rsponds o h sp disurbanc, approaching consan valu K as im bcoms larg. This immdia rspons, followd by asympoic approach o h nw sa sa, is characrisic of firs-ordr sysms. Bcaus h rspons dos no rack h sp inpu faihfully, h rspons is said o lag bhind h inpu; h firs-ordr sysm is somims calld a firs-ordr lag. rvisd 5 Jan

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw disurbanc.5.5 3 4 5 6 y rspons.5 K.5 3 4 5 6 im Figur.- firs-ordr rspons o iniial condiion and sp disurbanc.3 picwis ingraion of non-smooh disurbancs Th soluion (.-3 is applid ovr succding im inrvals, ach fauring an iniial condiion (from h prcding inrval and disurbanc inpu. ( K y( + ( K ( = y( + c. < < < < y x(d x(d (.3- Exampl: suppos rvisd 5 Jan 3

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw d x = + y = x ( y( = < < < < < (.3- In his problm, variabls, x, and y should b prsumd o hav appropria, if unsad, unis; in hs unis, boh gain and im consan ar of magniud. From (.3-, y( = ( ( + ( < < < < < (.3-3 Wih a zro iniial condiion and no disurbanc, h sysm rmains a quilibrium unil h ramp disurbanc bgins a =. Thn h oupu immdialy riss in rspons, lagging bhind h linar ramp. A =, h disurbanc cass, and h oupu dcays back oward quilibrium. rvisd 5 Jan 4

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw disurbanc.5.5.5 rspons.8.7.6.5.4.3.. 3 4 5 6 3 4 5 6 im.4 mulipl disurbancs and suprimposiion Sysms can hav mor han on inpu. Considr a firs-ordr sysm wih wo disurbanc funcions. + y( = Kx( + K x ( y( = y (.4- d Applying (.-3 and disribuing h ingral across h disurbancs, w find ha h ffcs of h disurbancs on y ar addiiv. ( y K K ( = y + x(d + x (d (.4- This addiiv bhavior is a happy characrisic of linar sysms. Thus anohr way o viw problm (.4- is o dcompos i ino componn problms. Tha is, dfin y = y + y + y (.4-3 H rvisd 5 Jan 5

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw and wri (.4- in hr quaions. W pu h iniial condiion wih no disurbancs, and ach disurbanc wih a zro iniial condiion. H + yh ( = d + y( = Kx( d + y( = K x ( d y ( = y y ( = H y ( = (.4-4 Equaions and iniial condiions (.4-4 can b summd o rcovr h original problm spcificaion (.4-. Th soluions ar y H ( = y K y( = K y( = ( x x (d (d (.4-5 and of cours hs soluions can b addd o rcovr original soluion (.4-. Thus w can viw h problm of mulipl disurbancs as a sysm rsponding o h sum of h disurbancs, or as h sum of rsponss from svral idnical sysms, ach rsponding o a singl disurbanc. Exampl: considr 3 + y = U( U( 3 y( = (.4-6 4 d 4 4 W firs plac h quaion in sandard form, in which h cofficin of y is +. + y = 3U( 4U( 3 y( = (.4-7 d Equaion (.4-7 shows us ha h im consan is, and ha h sysm rsponds o h firs disurbanc wih a gain of 3, and o h scond wih a gain of -4. Th soluion is y ( ( 3 ( 4U( 3 ( = + 3U( (.4-8 rvisd 5 Jan 6

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw In Figur.4-, h individual soluion componns ar plod as solid racs; hir sum, which is h sysm rspons, is a dashd rac. Noic how h firs-ordr lag rsponds o ach nw disurbanc as i occurs. disurbancs.5 rspons 4 3 - - -3-4 -5 4 6 4 6 im 8 8 Figur.4- firs-ordr rspons o mulipl disurbancs Wriing h sp funcions xplicily in soluion (.4-8 mphasizs ha paricular disurbancs do no influnc h soluion unil h im of hir occurrnc. For xampl, if hy wr omid, som dcpivly corrc bu inappropria rarrangmn would lad o rrors. y = + 3 = + 3 3 = + 3 ( ( ( ( ( 4 4 + 4 + 4 ( 3 ( 3 ( 3 (do no do his! (.4-9 This noaion a las implis ha wo of h xponnial funcions hav dlayd onss. Howvr, furhr corrc-bu-inappropria rarrangmn maks hings vn wors. rvisd 5 Jan 7

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw y = + = + = + 3 3 ( + 4 + 4 3 ( 3 + 4 ( 3 3 (do no do his! (.4- Th incorrc soluions ar plod wih (.4-8 in Figur.4-. Equaion (.4-9 has bcom disconinuous - h rspons aks non-physical laps a h ons of ach nw disurbanc. Equaion (.4- has los all dpndnc on h disurbancs and dcays from a non-physical iniial condiion. Evn wih h misaks, boh incorrc soluions lad o h corrc long-rm condiion. 5 rspons 5-5 4 6 im 8 soluion q.4.9 q.4. Figur.4- comparison of corrc and incorrc soluions.5 dlayd rspons o disurbancs Considr a sysm ha racs o a disurbanc, bu only afr som inrvning im inrval θ has passd. Tha is + y( = Kx( θ y( = y (.5- d Equaion (.5- shows h dpndnc of y, a any im, on h valu of x a arlir im - θ. Th soluion is wrin dircly from (.-3. ( K y( = y + x( θd (.5- W mus ingra h disurbanc considring h im dlay. Tak as an xampl a disurbanc x( occurring a im. Th plo shows h rvisd 5 Jan 8

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw disurbanc, as wll as h disurbanc as h sysm xprincs i, which bgins a im + θ. W could xprss his disurbanc-as-xprincd as som nw funcion x (, occurring a im + θ. disurbanc as i occurs x( disurbanc as xprincd by sysm x( - θ = x ( + θ x( - θ - θ Alrnaivly, w could dfin a nw im variabl = θ (.5-3 and wri h inpu as x(. Th ingral in (.5- bcoms, hn, x( θd = x(d = θ +θ x( d (.5-4 Thrfor, soluion (.5- bcoms ( K θ y( = y + x( d (.5-5 θ Exampl: considr a sp disurbanc a im = ha affcs h sysm 3 im unis lar. Using (.5-5 + y = x( 3 d x( = U( y( = (.5-6 rvisd 5 Jan 9

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw y = = = 3 3 3 3 3 U( = U( 3 = U( 5 = U( 5 U( d U( d + 3 3 [ ] 3 + 3 + 3 [ ] ( 5 [ ] U( d (.5-7 Figur.5- shows ha a ypical firs-ordr lag sp rspons occurs 3 im unis afr bing disurbd a =. disurbancs.5 4 6 8.8 rspons.6.4. 4 6 8 im Figur.5- sp rspons of firs ordr sysm wih dad im Th im dlay in rsponding o a disurbanc is ofn calld dad im. Dad im is diffrn from lag. Lag occurs bcaus of h combinaion of y and is drivaiv on h lf-hand sid of h quaion. Dad im rvisd 5 Jan

Spring 6 Procss Dynamics, Opraions, and Conrol.45 Lsson : Mahmaics Rviw occurs bcaus of a im dlay in procssing a disurbanc on h righhand sid..6 conclusion Plas bcom comforabl wih handling ODEs. Viw hm as sysms; idnify hir inpus and oupus, hir gains and im paramrs. rvisd 5 Jan