LINEARIZATION OF NONLINEAR EQUATIONS By Dominick Andrisani. dg x. ( ) ( ) dx

Similar documents
Partial Differential Equations

too many conditions to check!!

MAT2400 Assignment 2 - Solutions

Notes on iteration and Newton s method. Iteration

THE CONSERVATIVE DIFFERENCE SCHEME FOR THE GENERALIZED ROSENAU-KDV EQUATION

Numerical Methods for Finding Multiple Solutions of a Superlinear Problem

Castiel, Supernatural, Season 6, Episode 18

Appendix: The Laplace Transform

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

The z Transform. The Discrete LTI System Response to a Complex Exponential

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

In this document, if A:

CS321. Numerical Analysis and Computing

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

Diagonalization of Quadratic Forms. Recall in days past when you were given an equation which looked like

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

CS537. Numerical Analysis and Computing

Apply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.

Mathematical Series (You Should Know)

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

Coordinate Systems. Things to think about:

Find a formula for the exponential function whose graph is given , 1 2,16 1, 6

Inverse Matrix. A meaning that matrix B is an inverse of matrix A.

4. Linear Classification. Kai Yu

Fast Power Flow Methods 1.0 Introduction

Math 257: Finite difference methods

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Topic 9: Sampling Distributions of Estimators

x x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,

Chapter 2: Numerical Methods

September 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1

Math 312 Lecture Notes One Dimensional Maps

Chapter 4. Fourier Series

Lecture Notes for Analysis Class

Definition 2.1 (The Derivative) (page 54) is a function. The derivative of a function f with respect to x, represented by. f ', is defined by

Chapter 4 : Laplace Transform

Three-Step Iterative Methods with Sixth-Order Convergence for Solving Nonlinear Equations

INEQUALITIES BJORN POONEN

Special Modeling Techniques

Approximation of the Likelihood Ratio Statistics in Competing Risks Model Under Informative Random Censorship From Both Sides

( a) ( ) 1 ( ) 2 ( ) ( ) 3 3 ( ) =!

Sequences and Series of Functions

Chimica Inorganica 3

and the sum of its first n terms be denoted by. Convergence: An infinite series is said to be convergent if, a definite unique number., finite.

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

x x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula

Taylor Series (BC Only)

Lesson 10: Limits and Continuity

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

f t dt. Write the third-degree Taylor polynomial for G

HWA CHONG INSTITUTION JC1 PROMOTIONAL EXAMINATION Wednesday 1 October hours. List of Formula (MF15)

MATHEMATICS I COMMON TO ALL BRANCHES

wavelet collocation method for solving integro-differential equation.

Stability of Solution for Nonlinear Singular Systems with Delay

μ are complex parameters. Other

Lecture 3. Digital Signal Processing. Chapter 3. z-transforms. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev. 2016

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

Elementary Linear Algebra

Exponential Moving Average Pieter P

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

ECONOMICS 200A MATHEMATICAL HANDOUT

8. Applications To Linear Differential Equations

1 Approximating Integrals using Taylor Polynomials

Asymptotic Behaviour of Solutions of Certain Third Order Nonlinear Differential Equations via Phase Portrait Analysis

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

Precalculus MATH Sections 3.1, 3.2, 3.3. Exponential, Logistic and Logarithmic Functions

Infinite Sequences and Series

Chapter 10: Power Series

SNAP Centre Workshop. Basic Algebraic Manipulation

EXISTENCE OF CONCENTRATED WAVES FOR THE DIFFERENT TYPE OF NONLINEARITY. V. Eremenko and N. Manaenkova

Linear Regression Demystified

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.

Continuous and Discrete Modeling for IBIS-AMI

Curve Sketching Handout #5 Topic Interpretation Rational Functions

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

Probability and Statistics

CHAPTER 5. Theory and Solution Using Matrix Techniques

(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:

which are generalizations of Ceva s theorem on the triangle

CDS 101: Lecture 5.1 Controllability and State Space Feedback

Calculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.

Application of Digital Filters

A.1 Algebra Review: Polynomials/Rationals. Definitions:

PRACTICE FINAL/STUDY GUIDE SOLUTIONS

Introduction. Question: Why do we need new forms of parametric curves? Answer: Those parametric curves discussed are not very geometric.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Microscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions

LESSON 2: SIMPLIFYING RADICALS

Definition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.

ECONOMICS 100ABC MATHEMATICAL HANDOUT

Digital signal processing: Lecture 5. z-transformation - I. Produced by Qiangfu Zhao (Since 1995), All rights reserved

L 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!

Transcription:

LINEAIZATION OF NONLINEA EQUATIONS By Domiick Adrisai A. Liearizatio of Noliear Fctios A. Scalar fctios of oe variable. We are ive the oliear fctio (). We assme that () ca be represeted si a Taylor series epasio abot some poit as follows d d ( ) ( ) = ( ) = + = + = ( ) ( ) ( ) d! d + hiher order terms A liear approimatio for () ivolves taki oly the first two terms d( ) ( ) ( ) = + ( ) d This approimatio is most accrate if ( ) order terms are eliible. = is small so that the elected hiher

Eample: ( )= Epadi () abot = ives ( ) d ( ) + d ( ) = = + ( ) = = 4+ 4( ) = 4+ 4 Notice that this is a liear fctio of. The simplificatio reslted becase we evalated all oliear terms at the mber riht had side of the eqatio above at = =. Becase we evalated the terms o the = =, the oly term that depeds o is the - term, ad this is a liear term. This approimatio to () is valid ear =. We ca eerate aother approimatio to () by epadi () abot = ives ( ) = + ( ) = + ( ) = = + This secod approimatio is valid ear =. Clearly the liear approimatio depeds o the choice of referece poit. A. Scalar fctio of variables Give the oliear fctio (, ). This fctio ca be represeted by a Taylor series epasio abot, as follows

(, ) = (, ) +, ( ) +, ( = = = =! ( ) + ( ) =, = =, = + ( )( ) + h. o. t. =, = A liear approimatio of ca be obtaied by retaii the first three terms above (derlied). The two variables i this problem ca be associated toether i a vector as follows ( ) where = Eample: (, ) = cos ca be approimated abot =, = as follows,, = = = = (, ) = ( cos ) + ( cos ) ( ) ( si ) ( ) =, = = 4+ 4( ) + = 4+ 4 The same fctio ca be approimated abot = =π/ 4 π (, ) = cos( ) + ( cos ) ( ) π 4 =, = 4 π ( si ) ( ) π =, = 4 4 π = + 4( ) = π 4 4 Notice aai how importat the liearizatio poit or referece poit is to the liearized reslt.

A.3 Vector fctio of a vector of variables. Let ( ) be a vector of oliear fctios. Let be a vector of variables = = M M, A liear approimatio abot is + = ( ) ( ) where ( ) (, ) (, ) = M L L = vector = M L L L = Jacobia Matri = matri A.4 Accracy of liearized soltio. Whe we approimate ( ) by retaii oly the liear terms, we mst aratee that the deleted terms, i.e., the h.o.t. are eliible. This is tre oly whe is small, i.e. whe the pertrbatios from the referece poit are small. B. Liearizatio o Noliear Differetial Eqatios i First Order Form B. First order form Noliear differetial eqatios i first order form ca be writte as

where = M = (, ), ( ) =, M, =, M = M Note that represets specified forci fctios ad ( ) is a specified iitial coditio vector. m Eample B.a =, = + = + B. The referece or trim soltio Whe we were liearizi oliear fctios, we saw how importat the choice of referece poit was. I liearizi oliear differetial eqatios, we are also cocered with the referece abot which we liearize. However, we are ow iterested i obtaii a liearized soltio valid for all time. This reqires that we liearize arod a referece soltio, which is valid for all time. Let () t be a kow soltio to the oliear differetial eqatio with specified forci fctio t. i.e., () ad specified iitial coditio ( ) t ( ) = ( (), t ()) t ( ) () t is said to be the referece soltio to the oliear differetial eqatio. Eample B.b For the differetial eqatios ive i Eample B.a () t =, () t, = () t =

Eample B.c Eample B.d is a costat soltio to the oliear differetial eqatio. Verify this fact for yorself by sbstitti this soltio ito the differetial eqatio ive i Eample B.a. Please keep straiht i yor mid the differece betwee a differetial eqatio (e.. = ) ad a soltio to a differetial eqatio (e.. = for = ). For the differetial eqatios ive i Eample B.a () t = () t =, ẋ = [ ] is aother costat soltio to the oliear differetial eqatios. For the differetial eqatios ive i Eample B.a =± = =± = cost = cost ẋ = is a costat soltio to the oliear differetial eqatios for ay costat. B.3 Liearizatio abot a referece soltio Let (), t () t be a referece soltio. We ow wat to fid a liearized soltio to the oliear differetial eqatio abot this referece soltio. We aai epad ( ) i a Taylor series epasio abot ad i.e., = (, ) +, ( ) + =, ( ) = = = +h.o.t. The liear approimatio is obtaied by assri that eoh that the h.o.t. ca be elected. B.4 Defiitio of small distribtio variables ad are small Defie = = = For the liearized soltio to be valid, these pertrbatios mst be small.

B.5 Separatio of the liearized differetial eqatios ito two parts Assmi that the pertrbatios are small, we ca write the approimatio to the differetial eqatios as = (, ) + ( ) + ( ) we ca ow sbstitte the small pertrbatio variables + = (, ) + + I the eqatio above we have simplified the otatio with to deote, = =. Notice that the derlied terms are merically eqal from the defiitio of referece soltio. Sice they are eqal, they ca be cacelled ot leavi = + This is a set of liear small pertrbatio differetial eqatios. I smmary, the oriial oliear problem = (, ), ( ) with soltio t () for specified ipt (t) has bee decomposed ito two separate problems. The referece problem r = (, ) with iitial coditio ( ) with soltio () t to ipt () t The small pertrbatio problem with iitial coditio with soltio () t to ipt (). t = + ( ) = ( ) ( ) Fially the total approimate soltio is ive by the etire soltio procedre is show i Fire. t () = () t +t ()

B.6 O picki a referece soltio Ay soltio to (, ) = makes a ood referece soltio bt these soltios ca be hard to fid. A easier set of soltios are costat soltios i.e., soltios so that () t = ad ()= t costat for t ()= costat. For costat referece soltios, fidi the referece soltio to a oliear differetial eqatio becomes a problem of fidi the soltio to a oliear alebraic eqatio (, )= B.7 Liearizatio Eample t ( ) = = + a) Choice of eferece Soltio To simplify or choice, assme that the referece soltio is costat, i.e., = =. This reqires that = ad + =. These eqatios ca be satisfied wheever = ad = Vales of ad which satisfy these eqatios are =± where is ay costat =±

Eact Noliear Soltio ((t)) (e.., from merical iteratio) Noliear Differetial Eqatios = (, ), ( ) () t = specified Noliear eferece Problem r = (, ) ( ), specified where typically ad are costats Liear Small Pertrbatio Problem = + ( ) = ( ) ( ) () t = t () Noliear eferece Soltio (e.., costat soltio), r Liear Soltio (i.e., from Laplace Trasforms) () t Total Approimate Soltio t () = + t () ( ) = + ( ) t () = + t () Fire Soltio Procedres for Noliear Differetial Eqatios

We will cosider two differet referece soltios e.# (), () ( ) e.# (), () ( ) f t t f t t = + + =+ = = = = b) Small Pertrbatio Eqatios of Motio = + where = = A B = = =, = = =, Usi ef. # = =, = + + Usi ef. # = =, = + + + c. The Liear Soltio for eferece # () ẋ =

() ẋ = take Laplace trasforms () s () s () = () s s () () s () s () = () s mltiply () by s s () s s () = s () s mltiply () by set these eqal s () s = 4 () s + 4() s () s () s s() = 4() s + 4() s () ( s) [ s + 4]= s( ) ( ) + 4( s) s ( ) ( ) 4 () s = + ( ) s ( ) s + 4 s + 4 The first term o the riht ives iitial coditio respose. The secod term o the riht () s 4 cotais the trasfer fctio s () = s + 4. To fid () t take the iverse Laplace trasform. From () s () s () = () s () s = s() s () [ ] To fid the soltios () t ad () t yo mst be ive the ipt t () ad the iitial coditios ( ( ), ( ) ). The the soltios ca be fod si iverse Laplace trasforms. d) Total Soltio for eferece # () t () t = () t () t + () t () t = + + t () = () t + t () = + t () + ( ) ( ) = + ( ) () t () t

e) Commet: For this procedre to be valid the pertrbatios mst be small, i.e., all mst be small. Sppose we have the oliear problem with ( ) =. ( ) =. 99 t ( ) =. +. siωt the we ca se ef.# ad the we ca have ( ) =. ( ) =. t ( ) =. siωt O the other had if for the oliear problem we have ( ) =. ( ) =. 99 t ( ) =. siωt We wold se ef. # with ( ) =. ( ) =. t ( ) =. siωt

C. Otpt Eqatios Ofte oliear differetial eqatios are associated with oliear otpt eqatios. This may come abot i modeli the sesors aboard a aircraft. The sesors are ofte oliear fctios of the state vector ad cotrol vector. Otpt eqatios ca be epressed as follows. y = h(, ) where is y a p vector This ca also be liearized abot referece soltio ad as follows. h h y = y + y h(, ) + ( ) + ( ) The derlied terms are eqal by defiitio of y o the referece ad ca be cacelled ot o both sides of the eqatio. That leaves the liear small pertrbatio otpt eqatios i terms of small pertrbatio variables. h h y = + = C + D The total otpt eqatio i liear form is the ive by the followi. yt () y +yt () D. Stability Oe of may possible defiitios of dyamic stability for oliear systems is ive i terms of the eievales of the Jacobia matri, A = = If all the eievales of A have eative real parts we say that the referece soltio, (, ) is stable.., If at least oe of the eievale of A has a positive real part we say that the referece soltio, (, ), is stable. If at least oe eievales of A has a zero real part, ad if all the other eievales have eative real parts, we ca draw o coclsio abot the stability of the referece soltio, (, ). E. Cocldi Commets We have see how the soltio to oliear differetial eqatios ca be fod by decomposi the problem ito two simpler parts. The referece part is simpler becase it is ofte a oliear alebraic problem. The secod small pertrbatio part is simpler becase it ofte ivolves solvi liear differetial eqatios with costat coefficiets. The total approimate soltio to the oriial oliear differetial eqatio was show to be the sm of the two simpler parts.