CHAPTER 2 MODELING OF THREE-TANK SYSTEM

Size: px
Start display at page:

Download "CHAPTER 2 MODELING OF THREE-TANK SYSTEM"

Transcription

1 7 CHAPTER MODELING OF THREE-TANK SYSTEM. INTRODUCTION Te interacting tree-tank system is a typical example of a nonlinear MIMO system. Heiming and Lunze (999) ave regarded treetank system as a bencmark problem for reconfigurable control and observer based fault diagnosis. Te linear discrete state space models are essential for model based control tecniques. In tis capter, developments of linear discrete state space models for fault-free and faulty systems against actuator failure are presented. Te tree-tank system proposed by Hou et al (005) is used for analysis and syntesis of single objective and multi-objective reconfigurable control systems explained in Capter and Capter 4.. THREE-TANK SYSTEM DESCRIPTION Te scematic diagram of te coupled tree-tank system is sown in Figure.. It is composed of tree identical tanks wit a circular cross section of area S. Te tanks are interconnected by two cylindrical pipes wit a circular cross-section of area S C and outflow coefficients of tank and tank are a z and a z respectively. Te nominal inflows (q and q ) are located at tank and tank respectively. Te inflow rate can be continuously manipulated from 0 to a maximum flow rate of q max to maintain te tank level of max. Te measured variables are te level of tank ( ), tank ( ) and tank ( ). Te nominal outflow pipe as a cross section S C wit an outflow

2 8 coefficient a z and located at tank. Te control objective is to control level of tank and tank by manipulating te inflow rates q and q. Pump Pump q S q S c a z a z a z Tank Tank Tank Figure. Scematic diagram of tree-tank system in Equation (.) Te tree-tank system represented using te mass balance is given d = dt d = dt q - S a sgn ( - ) g( - ) z S S a sgn ( - ) g( - ) - S a sgn ( - ) g( - ) z z d q + S a sgn ( - ) = g( - ) - S a g dt S z z S (.) te Table.. Te pysical parameters of te tree tank system are presented in

3 9 Table. Pysical parameters of te tree-tank system Parameters Values Tank cross-section area S = 0.07 m Pipe cross- section area S C = S = S = S = m Pipe outflow coefficients a z = 0.5, a z = 0.579, a z = 0.7 Maximum level max = 0.68 m Maximum in-flow rate q max =. 0-4 m /s Te tree-tank system equations involve square-root nonlinearities and te flow-rates become proportional to te square root of te tank level. In control engineering, a normal operation of te system may be around an equilibrium point and te signals may be considered as small signals around te equilibrium. However, if te system operates around an equilibrium point and if te signals involved are small signals, ten it is possible to approximate te nonlinear system by a linear system. Suc a linear system is equivalent to te nonlinear system considered witin a limited operating range (Ogata 004). Te linearization procedure is presented in te following section.. DEVELOPMENT OF FAULT-FREE AND FAULTY MODELS Te discrete state space model development tecnique involves te following steps to obtain fault-free and faulty models:. Linearization of nonlinear equations around te operating point using Taylor s series expansion metod.. Discretization of continuous model for fault-free system.

4 0. Modeling of faulty system by introducing loss in control effectiveness in actuator. 4. Linearization and discretization of faulty system... Linear State Space Model for Fault-free System Equation (.) Te linearized state space model in continuous form is given in x(t) = Acx(t) + Bcu(t) y(t) = Ccx(t) (.) were, b b f f f 0 b b b b f f f, A C b b b f f f 0 d d d f, f, f dt dt dt f f 0 q q S f f B C 0 0 q q f f 0 q q S

5 x o x= x = o x o u q -q o u = = u q -q o were [ 0, 0, 0 ] and [q 0, q 0 ] are steady state operating points of level and flow rate respectively. a z S g a z S g a z S g b =, b =, b = S S S Te linearization tecnique is valid in te vicinity of te operating point. Te above nonlinear system is linearized around te following steady state operating points T [ ] = [ ] m o o o and [ q T T -4 o q o ] = [ ] 0 m / s as given below Te continuous state space model for te parameters in Table. is AC , BC 0 0 and CC Most practical systems are continuous-time systems. However, tey use discrete-time controller to obtain optimum performance.

6 Te discrete state space model is obtained by discretizing te continuous state space model wit sampling period T s = s as given in Equation (.) x(k +) = Ax(k) + Bu (k) y(k) = Cx(k) (.) were A , B and C x(k), u(k) and y(k) are te state, input and output vectors of discrete model respectively... Linear State Space Model for Faulty System Te actuator fault is modeled as a bias fault troug te control effectiveness factor wic represents te gain of actuator. If te gain is zero, ten te actuator is 00% effective. During normal operation, te actuator delivers te control signal witout any loss. Terefore, =0. Wen a fault occurs in te actuator, due to eiter partial blockage or aging, te actuator cannot deliver te control signal witout loss. Suc a failure condition can be represented by a reduced control effectiveness factor. Te magnitude of reflects te severity of actuator fault. Te actuator fault wit 80% loss of control effectiveness in actuator (pump) is considered for simulation. Te continuous and discrete faulty models ave te form given in Equations (.4) and (.5)

7 x (t) = A x (t) + Bcf u (t) f cf f f y (t) = C x (t) f cf f (.4) were x f (t) and y f (t) are te state and output vectors of faulty system in continuous form respectively and faulty system state matrix A cf = A c, faulty system input matrix B cf = (-) B c and faulty system output matrix C cf = C c. x (k+) = A x (k)+b u (k) f f f f f y (k) = C x (k) f f f (.5) were x f (k) and y f (k) are te state and output vectors of faulty system in discrete form respectively. A f, B f and C f are te state, input and output matrices of faulty system in discrete form respectively. given below: Model parameters for linearized system in continuous form are as A AC , B 0 0 and C C cf cf cf C given below: Model parameters for linearized system in discrete form are as A A , B and C C f f f

8 4 From te continuous model parameters, it is clear tat 80% loss of control effectiveness factor (actuator fault) on pump as caused canges in te input matrix, by noticing te value of te first element in te first column canges from (B c ) to.7 (B cf ). Te corresponding canges are reflected in all te elements of input matrix B f..4 SUMMARY Te detailed description of te interacting tree-tank system and first principle model are presented. Te discrete state space models for faultfree and faulty systems are derived using Taylor s series linearization tecnique. In te following capter, te conventional control tecnique using te derived discrete model will be discussed.

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

Runge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))

Runge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t)) Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001

Investigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001 Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x

More information

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Chapter 5 FINITE DIFFERENCE METHOD (FDM) MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential

More information

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1 Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.

More information

158 Calculus and Structures

158 Calculus and Structures 58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you

More information

Logarithmic functions

Logarithmic functions Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic

More information

MVT and Rolle s Theorem

MVT and Rolle s Theorem AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state

More information

LECTURE 14 NUMERICAL INTEGRATION. Find

LECTURE 14 NUMERICAL INTEGRATION. Find LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use

More information

De-Coupler Design for an Interacting Tanks System

De-Coupler Design for an Interacting Tanks System IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System

More information

Chapters 19 & 20 Heat and the First Law of Thermodynamics

Chapters 19 & 20 Heat and the First Law of Thermodynamics Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,

More information

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow 1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

More information

Continuity. Example 1

Continuity. Example 1 Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)

More information

1. State whether the function is an exponential growth or exponential decay, and describe its end behaviour using limits.

1. State whether the function is an exponential growth or exponential decay, and describe its end behaviour using limits. Questions 1. State weter te function is an exponential growt or exponential decay, and describe its end beaviour using its. (a) f(x) = 3 2x (b) f(x) = 0.5 x (c) f(x) = e (d) f(x) = ( ) x 1 4 2. Matc te

More information

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means

More information

Chapter 4 Derivatives [ ] = ( ) ( )= + ( ) + + = ()= + ()+ Exercise 4.1. Review of Prerequisite Skills. 1. f. 6. d. 4. b. lim. x x. = lim = c.

Chapter 4 Derivatives [ ] = ( ) ( )= + ( ) + + = ()= + ()+ Exercise 4.1. Review of Prerequisite Skills. 1. f. 6. d. 4. b. lim. x x. = lim = c. Capter Derivatives Review of Prerequisite Skills. f. p p p 7 9 p p p Eercise.. i. ( a ) a ( b) a [ ] b a b ab b a. d. f. 9. c. + + ( ) ( + ) + ( + ) ( + ) ( + ) + + + + ( ) ( + ) + + ( ) ( ) ( + ) + 7

More information

Section 15.6 Directional Derivatives and the Gradient Vector

Section 15.6 Directional Derivatives and the Gradient Vector Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,

More information

Definition of the Derivative

Definition of the Derivative Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

Polynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions

Polynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations

More information

HOMEWORK HELP 2 FOR MATH 151

HOMEWORK HELP 2 FOR MATH 151 HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,

More information

5.1 We will begin this section with the definition of a rational expression. We

5.1 We will begin this section with the definition of a rational expression. We Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat

More information

Physically Based Modeling: Principles and Practice Implicit Methods for Differential Equations

Physically Based Modeling: Principles and Practice Implicit Methods for Differential Equations Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

Taylor Series and the Mean Value Theorem of Derivatives

Taylor Series and the Mean Value Theorem of Derivatives 1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential

More information

Chapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3.

Chapter 1 Functions and Graphs. Section 1.5 = = = 4. Check Point Exercises The slope of the line y = 3x+ 1 is 3. Capter Functions and Graps Section. Ceck Point Exercises. Te slope of te line y x+ is. y y m( x x y ( x ( y ( x+ point-slope y x+ 6 y x+ slope-intercept. a. Write te equation in slope-intercept form: x+

More information

Theoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump

Theoretical Analysis of Flow Characteristics and Bearing Load for Mass-produced External Gear Pump TECHNICAL PAPE Teoretical Analysis of Flow Caracteristics and Bearing Load for Mass-produced External Gear Pump N. YOSHIDA Tis paper presents teoretical equations for calculating pump flow rate and bearing

More information

Chapter 1D - Rational Expressions

Chapter 1D - Rational Expressions - Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere

More information

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x

More information

Seepage Analysis through Earth Dam Based on Finite Difference Method

Seepage Analysis through Earth Dam Based on Finite Difference Method J. Basic. Appl. Sci. Res., (11)111-1, 1 1, TetRoad Publication ISSN -44 Journal of Basic and Applied Scientific Researc www.tetroad.com Seepage Analysis troug Eart Dam Based on Finite Difference Metod

More information

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY (Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative

More information

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as

More information

Recall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if

Recall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions

More information

1.5 Functions and Their Rates of Change

1.5 Functions and Their Rates of Change 66_cpp-75.qd /6/8 4:8 PM Page 56 56 CHAPTER Introduction to Functions and Graps.5 Functions and Teir Rates of Cange Identif were a function is increasing or decreasing Use interval notation Use and interpret

More information

Key Concepts. Important Techniques. 1. Average rate of change slope of a secant line. You will need two points ( a, the formula: to find value

Key Concepts. Important Techniques. 1. Average rate of change slope of a secant line. You will need two points ( a, the formula: to find value AB Calculus Unit Review Key Concepts Average and Instantaneous Speed Definition of Limit Properties of Limits One-sided and Two-sided Limits Sandwic Teorem Limits as x ± End Beaviour Models Continuity

More information

Lesson 6: The Derivative

Lesson 6: The Derivative Lesson 6: Te Derivative Def. A difference quotient for a function as te form f(x + ) f(x) (x + ) x f(x + x) f(x) (x + x) x f(a + ) f(a) (a + ) a Notice tat a difference quotient always as te form of cange

More information

MA119-A Applied Calculus for Business Fall Homework 4 Solutions Due 9/29/ :30AM

MA119-A Applied Calculus for Business Fall Homework 4 Solutions Due 9/29/ :30AM MA9-A Applied Calculus for Business 006 Fall Homework Solutions Due 9/9/006 0:0AM. #0 Find te it 5 0 + +.. #8 Find te it. #6 Find te it 5 0 + + = (0) 5 0 (0) + (0) + =.!! r + +. r s r + + = () + 0 () +

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x.

1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x. Problem. Let f x x. Using te definition of te derivative prove tat f x x Solution. Te function f x is only defined wen x 0, so we will assume tat x 0 for te remainder of te solution. By te definition of

More information

3.4 Algebraic Limits. Ex 1) lim. Ex 2)

3.4 Algebraic Limits. Ex 1) lim. Ex 2) Calculus Maimus.4 Algebraic Limits At tis point, you sould be very comfortable finding its bot grapically and numerically wit te elp of your graping calculator. Now it s time to practice finding its witout

More information

DYNAMIC MODELING OF ORGANIC RANKINE CYCLE (ORC) SYSTEM FOR FAULT DIAGNOSIS AND CONTROL SYSTEM DESIGN

DYNAMIC MODELING OF ORGANIC RANKINE CYCLE (ORC) SYSTEM FOR FAULT DIAGNOSIS AND CONTROL SYSTEM DESIGN DYNAMIC MODELING OF ORGANIC RANKINE CYCLE (ORC SYSTEM FOR FAULT DIAGNOSIS AND CONTROL SYSTEM DESIGN Sungjin Coi and Susan Krumdieck University of Canterbury, Private Bag 48, Cristcurc 84 New Zealand sungjin.coi@pg.canterbury.ac.nz

More information

Practice Problem Solutions: Exam 1

Practice Problem Solutions: Exam 1 Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical

More information

Chapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods

Chapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary

More information

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014 Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.

More information

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y

More information

The derivative function

The derivative function Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

The Verlet Algorithm for Molecular Dynamics Simulations

The Verlet Algorithm for Molecular Dynamics Simulations Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical

More information

Differential Calculus (The basics) Prepared by Mr. C. Hull

Differential Calculus (The basics) Prepared by Mr. C. Hull Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit

More information

The Derivative as a Function

The Derivative as a Function Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips

More information

MATH CALCULUS I 2.1: Derivatives and Rates of Change

MATH CALCULUS I 2.1: Derivatives and Rates of Change MATH 12002 - CALCULUS I 2.1: Derivatives and Rates of Cange Professor Donald L. Wite Department of Matematical Sciences Kent State University D.L. Wite (Kent State University) 1 / 1 Introduction Our main

More information

Use of fin analysis for determination of thermal conductivity of material

Use of fin analysis for determination of thermal conductivity of material RESEARCH ARTICLE OPEN ACCESS Use of fin analysis for determination of termal conductivity of material Nea Sanjay Babar 1, Saloni Suas Desmuk 2,Sarayu Dattatray Gogare 3, Snea Barat Bansude 4,Pradyumna

More information

Chapter 1. Density Estimation

Chapter 1. Density Estimation Capter 1 Density Estimation Let X 1, X,..., X n be observations from a density f X x. Te aim is to use only tis data to obtain an estimate ˆf X x of f X x. Properties of f f X x x, Parametric metods f

More information

MATH1151 Calculus Test S1 v2a

MATH1151 Calculus Test S1 v2a MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums

More information

Lecture 21. Numerical differentiation. f ( x+h) f ( x) h h

Lecture 21. Numerical differentiation. f ( x+h) f ( x) h h Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However

More information

Chemical Engineering & Process Techniques

Chemical Engineering & Process Techniques emical Engineering & Process Tecniques eview Article eedback ontrol for Liquid Level in a Gravity-Drained Multi-Tank System Larry K Jang* Department of emical Engineering, alifornia State University, USA

More information

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =

Test 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 = Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS INTODUCTION ND MTHEMTICL CONCEPTS PEVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips of sine,

More information

Finding and Using Derivative The shortcuts

Finding and Using Derivative The shortcuts Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex

More information

Chemical Engineering & Process Techniques

Chemical Engineering & Process Techniques emical Engineering & Process Tecniques eview Article eedback ontrol for Liquid Level in a Gravity-Drained Multi-Tank System Larry K. Jang* Department of emical Engineering, alifornia State University,

More information

ETNA Kent State University

ETNA Kent State University Electronic Transactions on Numerical Analysis. Volume 34, pp. 14-19, 2008. Copyrigt 2008,. ISSN 1068-9613. ETNA A NOTE ON NUMERICALLY CONSISTENT INITIAL VALUES FOR HIGH INDEX DIFFERENTIAL-ALGEBRAIC EQUATIONS

More information

Chapter 9 - Solved Problems

Chapter 9 - Solved Problems Capter 9 - Solved Problems Solved Problem 9.. Consider an internally stable feedback loop wit S o (s) = s(s + ) s + 4s + Determine weter Lemma 9. or Lemma 9. of te book applies to tis system. Solutions

More information

AVL trees. AVL trees

AVL trees. AVL trees Dnamic set DT dnamic set DT is a structure tat stores a set of elements. Eac element as a (unique) ke and satellite data. Te structure supports te following operations. Searc(S, k) Return te element wose

More information

Section 3.1: Derivatives of Polynomials and Exponential Functions

Section 3.1: Derivatives of Polynomials and Exponential Functions Section 3.1: Derivatives of Polynomials and Exponential Functions In previous sections we developed te concept of te derivative and derivative function. Te only issue wit our definition owever is tat it

More information

The distance between City C and City A is just the magnitude of the vector, namely,

The distance between City C and City A is just the magnitude of the vector, namely, Pysics 11 Homework III Solutions C. 3 - Problems 2, 15, 18, 23, 24, 30, 39, 58. Problem 2 So, we fly 200km due west from City A to City B, ten 300km 30 nort of west from City B to City C. (a) We want te

More information

(4.2) -Richardson Extrapolation

(4.2) -Richardson Extrapolation (.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as

More information

Exercise 19 - OLD EXAM, FDTD

Exercise 19 - OLD EXAM, FDTD Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms

More information

Section 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is

Section 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit

More information

An optimum design of robotic food handling by using Burger model

An optimum design of robotic food handling by using Burger model DOI 10.1007/s11370-008-0032-5 ORIGINAL RESEARCH PAPER An optimum design of robotic food andling by using Burger model Naoki Sakamoto Mitsuru Higasimori Tosio Tsuji Makoto Kaneko Received: 28 February 2008

More information

1 Introduction to Optimization

1 Introduction to Optimization Unconstrained Convex Optimization 2 1 Introduction to Optimization Given a general optimization problem of te form min x f(x) (1.1) were f : R n R. Sometimes te problem as constraints (we are only interested

More information

Notes on wavefunctions II: momentum wavefunctions

Notes on wavefunctions II: momentum wavefunctions Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles

More information

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c) Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''

More information

Manipulator Dynamics (1) Read Chapter 6

Manipulator Dynamics (1) Read Chapter 6 Manipulator Dynamics (1) Read Capter 6 Wat is dynamics? Study te force (torque) required to cause te motion of robots just like engine power required to drive a automobile Most familiar formula: f = ma

More information

lim 1 lim 4 Precalculus Notes: Unit 10 Concepts of Calculus

lim 1 lim 4 Precalculus Notes: Unit 10 Concepts of Calculus Syllabus Objectives: 1.1 Te student will understand and apply te concept of te limit of a function at given values of te domain. 1. Te student will find te limit of a function at given values of te domain.

More information

Derivation Of The Schwarzschild Radius Without General Relativity

Derivation Of The Schwarzschild Radius Without General Relativity Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:

More information

NCCI: Simple methods for second order effects in portal frames

NCCI: Simple methods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames Tis NCC presents information

More information

Pre-lab Quiz/PHYS 224 Earth s Magnetic Field. Your name Lab section

Pre-lab Quiz/PHYS 224 Earth s Magnetic Field. Your name Lab section Pre-lab Quiz/PHYS 4 Eart s Magnetic Field Your name Lab section 1. Wat do you investigate in tis lab?. For a pair of Helmoltz coils described in tis manual and sown in Figure, r=.15 m, N=13, I =.4 A, wat

More information

An Order Quantity Decision System for the Case of Approximately Level Demand

An Order Quantity Decision System for the Case of Approximately Level Demand Capter 5 An Order Quantity Decision System for te Case of Approximately Level Demand Demand Properties Inventory problems exist only because tere are demands; oterwise, we ave no inventory problems. Inventory

More information

Math 34A Practice Final Solutions Fall 2007

Math 34A Practice Final Solutions Fall 2007 Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen

More information

Digital Filter Structures

Digital Filter Structures Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical

More information

Exam 1 Solutions. x(x 2) (x + 1)(x 2) = x

Exam 1 Solutions. x(x 2) (x + 1)(x 2) = x Eam Solutions Question (0%) Consider f() = 2 2 2 2. (a) By calculating relevant its, determine te equations of all vertical asymptotes of te grap of f(). If tere are none, say so. f() = ( 2) ( + )( 2)

More information

Section 3: The Derivative Definition of the Derivative

Section 3: The Derivative Definition of the Derivative Capter 2 Te Derivative Business Calculus 85 Section 3: Te Derivative Definition of te Derivative Returning to te tangent slope problem from te first section, let's look at te problem of finding te slope

More information

3 Parabolic Differential Equations

3 Parabolic Differential Equations 3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded

More information

Consider the element shown in Figure 2.1. The statement of energy conservation applied to this element in a time period t is that:

Consider the element shown in Figure 2.1. The statement of energy conservation applied to this element in a time period t is that: . Conduction. e General Conduction Equation Conduction occurs in a stationary medium wic is most liely to be a solid, but conduction can also occur in s. Heat is transferred by conduction due to motion

More information

Pre-Calculus Review Preemptive Strike

Pre-Calculus Review Preemptive Strike Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly

More information

Chapter 2 Limits and Continuity

Chapter 2 Limits and Continuity 4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(

More information

Problem Solving. Problem Solving Process

Problem Solving. Problem Solving Process Problem Solving One of te primary tasks for engineers is often solving problems. It is wat tey are, or sould be, good at. Solving engineering problems requires more tan just learning new terms, ideas and

More information

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0 3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

The Priestley-Chao Estimator

The Priestley-Chao Estimator Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are

More information

The Laplace equation, cylindrically or spherically symmetric case

The Laplace equation, cylindrically or spherically symmetric case Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,

More information

Derivatives. By: OpenStaxCollege

Derivatives. By: OpenStaxCollege By: OpenStaxCollege Te average teen in te United States opens a refrigerator door an estimated 25 times per day. Supposedly, tis average is up from 10 years ago wen te average teenager opened a refrigerator

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information