Application of a Generalized Budan-Fourier Algorithm to Stability Analysis of Time-Delay Systems
|
|
- Christine Owen
- 5 years ago
- Views:
Transcription
1 Application of a Generalized Budan-Fourier Algorithm to Stability Analysis of Time-Delay Systems Bernard Hanzon 1 Joint Work With Oliver Mason 2 1 Edgeworth Centre for Financial Mathematics, Department of Mathematics, University College Cork, Ireland 2 Hamilton Institute, NUI Maynooth, Ireland 10 July / 32
2 Introduction Introduction Differential-delay systems have many applications in various fields including financial mathematics, control etc. Here we consider classical LTI differential-delay systems and investigate how to determine whether the system is stable, and if not, where the unstable poles are located (in the complex plane). 1 / 32
3 Introduction Introduction This problem has been around for a long time. We will follow an exposition of [Luc Habets 1992] of how one could approach this problem. However we combine the method suggested there with a relatively recent Generalized Budan Fourier (GBF) algorithm, which allows one to determine the zeros of any real Exponential-Polynomial-Trigonometric (EPT) function on any finite interval of the real axis with arbitrary precision under the assumption that we can determine the (-,0 or +) of an EPT function in any given real point. The GBF algorithm will also allow us to isolate all the (finitely many) unstable poles with arbitrary precision, in the complex plane, if there are any such poles. 2 / 32
4 Stability Analysis for Linear Differential-Delay Systems Stability Analysis for Linear Differential-Delay Systems Consider a linear differential-delay system with state equation: ẋ(t) = Fx(t) + k j=1 G jx(t τ j ) + Lu(t), x(t) R n, τ j > 0, j = 1, 2,..., k This system is stable if the poles of the transfer function all lie in the open left half plane (so all have negative real part). The poles are the zeros of the denominator of the transfer function, which is given by (assuming minimality of the state space realization) f (z) = det(zi F k j=1 G je zτ j ) This function is also the characteristic function of the relevant differential delay operator. Consider the homogeneous polynomial q(u, v, w 1,..., w k ) = det(ui vf k j=1 G jw j ) Then f (z) = q(z, 1, e zτ 1,..., e zτ k ) 3 / 32
5 Stability Analysis for Linear Differential-Delay Systems Stability Analysis for Linear Differential-Delay Systems We can also write this in the form f (z) = z n + n 1 h=0 p h(e τ1z,..., e τkz ).z h, z C, where the p h are polynomials in k variables. Habets[1992] shows that f (z) has only a finite number of zeros in the closed right half plane. Stability holds if this number is actually zero. As f (z) is an entire function (analytic on the whole complex plane) we can count the number of zeros inside any Jordan curve in the complex plane by counting the number of windings around the origin that f (z) makes when z traverses (once) along the Jordan curve in counterclockwise direction. 4 / 32
6 Stability Analysis for Linear Differential-Delay Systems Stability Analysis for Linear Differential-Delay Systems By taking a Jordan curve consisting of an interval on the imaginary axis and a large semi-circle in the right half plane, one can reduce the problem to counting the total increase in the argument of f (z) along a sufficiently large interval on the imaginary axis, as the contribution to the increase in the argument of moving along the semi-circle is known to be nπ if the semi-circle is sufficiently big. This leads to Theorem 2.7 of [Habets 1992] which says that the number of right half plane zeros is n 2 1 π.totarg(f, [0, ik]) where totarg(f, [0, ik]) denotes the net increase in argument of f (z) when z traverses along the imaginary axis from the origin to the number ik, for some sufficiently large K. 5 / 32
7 Stability Analysis for Linear Differential-Delay Systems Stability Analysis for Linear Differential-Delay Systems Finally the net increase in argument can actually be found by counting the number of crossings by f (z) of the real and of the imaginary axis, when z runs from 0 to ik. This in turn boils down to locating the zeros of Re[f (z)] and Im[f (z)], and determining the sign of the other at any location where one of these is zero. See the picture from [Habets,1992] 6 / 32
8 EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computing Science Memorandum CaSaR A Reliable Stability Test for Exponential Polynomials L.C.G.J.M. Habets Eindhoven, NoveJIlber 1992 The Netherlands
9 1m - Qll.is ""'!I""' Re ~ Qxis Figure 1: Jordan-curve JR Definition 2.3 Let R E R+. Then the half circle GR is defined as GR := {z EC Ilzl = RARe (z) ~ O}, the part IR of the imaginary axis as IR := {z Eel Re(z) = 0 A Izl < R}, and the Jordan-curve JR as JR:= GRUIR. This Jordan-curve is traversed in counter clockwise direction as depicted in Figure 1. Now suppose that the exponential polynomial f has no zeros on the imaginary axis. Let NJ denote the number of zeros of f in the open RHP. When R becomes large enough, JR will enclose all the NJ zeros of f. Hence NJ = lim _1_i f'(z) dz. R-+oo 21r~.7R f(z) (5) The computation of the integral (5) can be split into two parts: the integral over the half circle GR, and over the imaginary axis IR. The first term can be determined analytically, but before we can state this result we need some preliminary lemmas. Lemma 2.4 Let f(z) be an exponential polynomial of the form n-l f(z) = zn + ~ Pi(e-1"I Z,...,e-1""z)zi. i=o (6) Define A() d ( ( -1"IZ -1""Z)) Z := dz pn-l e,...,e. (7) Then for large values of Izl, such that Re (z) ~ 0 we have f'(z) = ~+ A(z) +O(~). f(z) z Z z2 (8) 3
10 1m -axis "Va. 'II. - '/a. 'r '/l... Ill. 'Iz..1/'1..'/ Figure 2: Counting the number of encirclements The only problem left is to track the curve r of I(f.W) for w E [0, K max ] in an accurate way that ascertains that all the intersections with the real axis (when n is odd) or the imaginary axis (when n is even) are detected. However, in an algorithm only a finite number of points on the curve r can be calculated. Therefore the tracking problem can be reformulated as the question to find a method for the selection of a finite number of points on the curve r in such a way that all the intersections of r with the real and imaginary axis can be detected from this finite set of points. For this purpose, linear search in the parameter space with constant step length I is until now most commonly used. Unfortunately, this method is not always reliable as will be illustrated by an example in the next section. Therefore we propose another, more reliable method to overcome this problem. The main idea is to make the step length variable depending on the curvature of the curve r. For a curve r = {(u(w),v(w)) IwE [O,K]} in a two-dimensional plane, parametrized by the variable w, the curvature in a point x(wo) = (u(wo), v(wo)) at w = wo, is given by u(wo)v(wo) - ij(wo)u(wo) (20) (u2(wo) + ij2(wo))3/2 (See for example [8, pp ]; the formula can be found literally in [2, p. 590]). The curvature in a point x(wo) on r is a measure for the rate of change of the tangent in x(wo), when proceeding along the curve. For example, a circle with radius R has in every point a curvature of *. It is easily seen that in a small neighborhood of the point x on the curve r, with curvature kx, the curve r itself behaves like a circle with radius Ik~ I. This implies that to track the curve accurately we have to take small steps along the curve when the curvature is large in absolute value, and we can take somewhat larger steps when the absolute value of the curvature is small. In this way we obtain the following rule: IcurvatureI X step length along the curve = constant (21) Note that the curve {(u(w), v(w)) I w E [0, K]} is parametrized by wand not by its length along the curve. The length of the curve between Wo and Wo + ~w is given by l wo +oo.w vu2(w) +ij2(w)dw. wo (22) For small values of ~w, this integral is approximately vu2(wo) + ij2(wo)~w. Substitution of this formula and (20) in (21) yields lu(wo)v(wo) - ij(wo)u(wo) I. ~ = C u2(wo) + ij2(wo) W. 8
11 Realization of Exponential-Polynomial-Trigonometric functions Realization of Exponential-Polynomial-Trigonometric functions So what s new? Firstly note that f (z) is an exponential-polynomial-trigonometric (EPT) function. Therefore it can be realized as the impulse response function of a linear system! So there exists a triple (A, b, c), A : N N, b : N 1, c : 1 N such that f (z) = ce Az b. Actually it is not too hard to show that an upper bound of the McMillan degree (i.e N in a minimal realization) is given by ( ) m + k + 1. k / 32
12 Realization of Exponential-Polynomial-Trigonometric functions Realization of Exponential-Polynomial-Trigonometric functions The real as well as the imaginary part of f (z) on the imaginary axis are again EPT functions and we can show that their McMillan degree is bounded by ( ) m + k (m + 1). k + 1 A realization for f (z) can be found in various ways. For instance one can use well-known realization formulas for the sum and the product of realizations of EPT functions (see e.g. [H.-Ober 2002] and references given there) and the well-known realizations of polynomials and scalar exponential functions. 8 / 32
13 It follows that the problem can be solved once we have a method to find the zeros of any given real EPT function on a finite real interval. This can be done using the GBF algorithm [H.-Holland 2009;2010] which was first used in the context of non-negativity of impulse-response functions of Finite Dimensional LTI systems and applications in financial matematics (non-negativity of Nelson-Siegel forward rate curves, and one-sided and two-sided EPT probability densities). As this method has not been used before in the context of locating roots in the complex plane we will give an overview of it here. 9 / 32
14 Consider the class of EPT functions k y(t) = Re p j (t)e µ j t, j=0 where p j (t) C[t] polynomials with complex coefficients, j = 0, 1, 2,..., k, µ j C, j = 0, 1, 2,..., k. This class contains the polynomials, the exponential functions and the scaled trigonometric polynomials (such as sin(νt), cos(νt), ν R but not tan(νt),..), and products and sums of such functions (the set of functions is a ring!). 10 / 32
15 A (real) EPT function can be written as y(t) = ce At b, (A, b, c) minimal and real. Distinguish three cases: (1) y(t) : polynomial iff σ(a) = {0}; then A nilpotent, A n = 0, e At = I + At A2 t (n 1)! An 1 t n 1 (2) If σ(a) R, real, then by bringing A in Jordan canonical form, one can show y(t) = ce At b = k h=0 p h(t)e λht, p h (t) R[t], real polynomial, λ 0 = 0 and λ h R, real eigenvalue of A, h=1,2,...,k. These will be called real exponential polynomials ( EP class). If deg(p h ) = 0, h = 0, 1, 2,..., K, we will speak of the exponential sums class ( E ) 11 / 32
16 (3) If σ(a) C the function will be called a function of the exponential-polynomial-trigonometric EPT class. In this case: { k } ce At b = Re h=0 p h(t)e λ ht, p h (t) C[t], λ h C, h = 0, 1, 2,..., k. This class/these classes go under many different names: quasi-exponential matrix-exponential exponential-polynomial; in systems theory they were sometimes called Bohl functions: this is now fading out. 12 / 32
17 Now consider the problem of determining the sign-changing zeros of such an EPT function on a finite interval (a, u). sign-changing zeros : zeros where the sign of the function changes from + to or vice versa. (EPT functions are real analytic, hence the zeros are isolated). CAN WE CONSTRUCT AN ALGORITHM THAT PRODUCES ALL SIGN-CHANGING ZEROS OF AN EPT FUNCTION ON A GIVEN FINITE INTERVAL? 13 / 32
18 For polynomials: Yes! Consider polynomial of degree n 1. Use Budan-Fourier sequence of higher order derivatives: p (0) (t) = p(t) = c.e At.b, σ(a) = {0} p (1) (t) = c.ae At.b, p (2) (t) = c.a 2.e At.b,. p (n) (t) = c.a n.e At.b 0, (A n = 0) We know k = n is the smallest integer value for which p (k) 0 so p (n 1) is a non-zero constant and p (n 2) is either strictly increasing (if p (n 1) > 0) or strictly decreasing (if p (n 1) < 0). 14 / 32
19 Let a < u. If p (n 2) (a) and p (n 2) (u) have the same (non-zero) sign then p (n 2) (t) 0, t [a, u]. If p (n 2) (a) and p (k 2) (u) have opposite (non-zero) signs then there is exactly one sign-changing zero of p (n 2) (t) on (a, u)]. This zero can be calculated using bisection, with arbitrary precision! One way of viewing this is as follows: Definition An open interval (a, u) R will be called simple for the function f : R R, f continuous, if f has at most one sign-changing zero on (a, u). 15 / 32
20 Remark. Suppose f has simple interval (a, u). Then if lim Sign [f (a + ɛ)] = lim Sign [f (u + ɛ)] 0 (1) ɛ 0 ɛ 0 then there is no sign-changing zero of f in (a, u). If lim ɛ 0 Sign [f (a + ɛ)] lim ɛ 0 Sign [f (u + ɛ)] < 0, then bisection gives the unique sign-changing zero c (a, u) such that f (c) = 0. (Of course the case f 0 on (a, u) can be handled in a straightforward manner. ) 16 / 32
21 Definition A grid {a 0, a 1,..., a N }, a 0 = a, a N = u, a 0 < a 1 <... < a N, is called simple for f if each interval (a k 1, a k ), k = 1, 2,..., N is simple for f. Remark. Given a simple grid the sign-changing zeros of f on the interval (a, u) can be found all and with arbitrary precision, using a bisection algorithm. 17 / 32
22 Now note that h {1, 2,..., n} the zeros of p (h) (t), together with the boundary points a, u, form a simple grid for p (h 1) (t), where p (0) (t) p(t) is the original polynomial. Therefore the Budan-Fourier sequence gives a guaranteed method to find all the real zeros of the original polynomial p with arbitrary precision. 18 / 32
23 Simple grid properties and a Generalized Budan-Fourier sequence for EP functions Simple grids have some nice properties: (i) If {a = a 0, a 1,..., a N = u} is a simple grid for f, then also for f.g, where g(t) 0, t (a, u)\{a 0, a 1,..., a N } and g continuous on (a, u). (ii) If {a 0, a 1,..., a N } consists of the boundary points together with the sign-changing zeros of a function h, h continuous on (a, u); then the same grid is obtained if d is replaced by d.k for any continuous function d with d(t) 0, t (a, u). 19 / 32
24 We can now construct a GBF sequence for the EP class! Let g(t) = ce At b, σ(a) R, (A, b, c) minimal and real. Let λ 1, λ 2,..., λ n denote the (real) eigenvalues of A : n n. Let g (0) (t) := g(t) = c.e At.b, g (h) := c(a λ h.i )... (A λ 1.I ).e At.b, h = 1, 2,..., n 1, g (n) (t) := c(a λ n.i )... (A λ 1.I )e At b 0 (Cayley-Hamilton!) 20 / 32
25 Theorem This is a Generalized Budan-Fourier sequence, for any pair of boundary points {a, u}, a < u. 21 / 32
26 Proof Note that for every h {0, 1, 2,..., n 1}, g (h+1) is obtained from g (h) by applying the formula: g (h+1) = e λ h+1.t d [ ] e λh+1.t.g (h) (t) dt hence d [ ] e λh+1.t.g (h) (t) = e λh+1.t.g (h+1) (t) dt Using (i),(ii) and the fact that e λ h+1.t > 0, t R, it follows that the sign-changing zeros of g h+1) (t) (which are also the sign-changing zeros of e λ h+1.t g (h+1) (t)), together with the boundary points {a, u}, form a simple grid for e λ h+1.t g (h) (t), hence for g (h) (t). 22 / 32
27 As g (n) 0 because (A λ n I )(A λ n 1 I )... (A λ 1 I ) = 0, it follows that g (n 1) has constant sign (no zeros). The theorem now follows by induction. QED 23 / 32
28 Now consider a general EPT function g(t) = c.e At.b, where σ(a) can contain complex elements, but g(t) R, t R. Assuming (A, b, c) minimal, it follows that if λ σ(a) with multiplicity m(λ), then λ σ(a), with the same multiplicity m(ˆλ) = m(λ). Let us order the eigenvalues λ 1, λ 2,..., λ N such that the complex conjugates come in pairs λ h = θ h + iν h, λ h+1 = θ h iν h, with ν h > / 32
29 Theorem A GBF sequence with boundary points {0, u}, for g(t) = ce At b, (A, b, c) minimal, real, is given by: g (0) (t) = g(t), with boundary points {0, u}; if λ h R : g (h) = c(a λ h I )(A λ h 1 I )... (A λ 1 I ).e At.b with boundary points {0, u}; if λ h = θ h + [ iν h, λ h+1 = θ h iν h, ν h > 0 then] g (h) := Im e λ h t c(a λ h I )... (A λ 1 I )e At b with extended set of boundary points {0, π ν h, 2 π ν h,..., u π/ν h π ν h, u}; and g (h+1) := c(a λ h+1 I )(A λ h I )... (A λ 1 I )e At b, with (same) extended set of boundary points {0, π ν h, 2 π ν h,..., u π/ν h π ν h, u} 25 / 32
30 Note that (A λ h+1 I )(A λ h I ) = (A λ h )(A λ h ) is a real matrix. Notation: x = max{n N; n x}. 26 / 32
31 Proof. Consider g (h 1) (t) = c(a λ h 1 I )... (A λ 1 I )e At b, and consider an open interval J l = (l π ν h, (l + 1) π ν h ), l {0,..., u φ/ν h }. On J l the function sin(ν h t) has no zero, hence does not change sign. On I l (0, u) a simple grid for g (h 1) can be obtained by calculating the ( sign-changing ) zeros of sin 2 (ν h t) d g (h 1) (t)e θ h t dt = sin(ν h t) 27 / 32
32 sin(ν h t).( d dt g (h 1) (t))e θht + sin(ν h t).g (h 1) (t)( θ h )e θht + g (h 1) (t).e θht ν h cos(ν h t) = Im [ e θht (cos(ν h t) + i sin(ν h t)).( d dt g (h 1) (t) + g (h 1) (t)( θ h iν h )) ] = Im [ e θ ht+iν h t ( d dt λ h)(g (h 1) (t)) ] [ = ] Im e λ h t c.(a λ h I )(A λ h 1 I )... e At.b = g (h) (t). 28 / 32
33 Now, perhaps surprisingly, e θ h t d sin(ν h t) dt [ g (h) (t) ] = d2 dt 2 g (h 1) (t) 2θ h d dt g (h 1) (t) + (θ 2 + ν 2 h )g (h 1) (t) = ( d dt λ h I )( d dt λ h)g (h 1) (t) = c(a λ h I )(A λ h I )(A λ h 1 I )... (A λ 1 I )e At.b =: g (h+1) (t), where we use λ h+1 = λ h. Therefore the sign-changing zeros of g (h+1) (t), together with the (extended set of) boundary points {0, π ν h, 2 π ν h,..., u π/ν h π ν h, u} form a simple grid for g (h) (t) and the sign-changing zeros and the sign-changing zeros of g (h) (t), together with the (extended set of) boundary points {0, π ν h, 2 π ν h,..., u π/ν h π ν h, u} form a simple grid for g (h 1) (t). The theorem now follows by induction (in combination with the proof of the previous theorem). QED 29 / 32
34 Final remarks and conclusions Final remarks and conclusions Applying the GBF bisection method to the real and imaginary parts of the characteristic function of a differential-delay system one can determine stability. This provides an alternative to existing algebraic methods. A further application is to isolate the right-half plane zeros of the characteristic function. This can be done by using the fact that the increase of the argument over any Jordan curve that takes the form of a rectangle with horizontal and vertical sides, can be obtained again by determining the zero crossings of the real and imaginary parts of f (z) over the sides of such a rectangle. That will reveal the number of zeros inside the rectangle. By systematically cutting the rectangle into smaller ones one can isolate the complex zeros with arbitrary precision (as there only finitely many). We are not aware of other methods that can do this. 30 / 32
35 Final remarks and conclusions Final remarks and conclusions We are developing software that can handle arbitrary matrices F, G j, j = 1, 2,..., k as inputs. Presently we are considering possible generalizations to more general differential-delay systems (such as neutral systems) based on results of [H-Holland IWOTA 2010] concerning positivity of the tail of an EPT function. 31 / 32
36 Thank You Thank you! 1 For references and links to software contact us at b.hanzon@ucc.ie 1 Research support for B.Hanzon by Science Foundation Ireland, grant nrs RFP2007-MATF802 and 07/MI/ / 32
The Generalized Laplace Transform: Applications to Adaptive Control*
The Transform: Applications to Adaptive * J.M. Davis 1, I.A. Gravagne 2, B.J. Jackson 1, R.J. Marks II 2, A.A. Ramos 1 1 Department of Mathematics 2 Department of Electrical Engineering Baylor University
More informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationTopic 4 Notes Jeremy Orloff
Topic 4 Notes Jeremy Orloff 4 auchy s integral formula 4. Introduction auchy s theorem is a big theorem which we will use almost daily from here on out. Right away it will reveal a number of interesting
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationLinear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form
Linear Planar Systems Math 246, Spring 2009, Professor David Levermore We now consider linear systems of the form d x x 1 = A, where A = dt y y a11 a 12 a 21 a 22 Here the entries of the coefficient matrix
More informationMid Term-1 : Solutions to practice problems
Mid Term- : Solutions to practice problems 0 October, 06. Is the function fz = e z x iy holomorphic at z = 0? Give proper justification. Here we are using the notation z = x + iy. Solution: Method-. Use
More informationChapter 6: The Laplace Transform. Chih-Wei Liu
Chapter 6: The Laplace Transform Chih-Wei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace
More informationSolutions Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du
Solutions hapter 9 dode 9 asic Solution Techniques 9 hoose one or more of the following differential equations, and then: (a) Solve the equation directly (b) Write down its phase plane equivalent, and
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationComplex Variables. Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit.
Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit. 1. The TI-89 calculator says, reasonably enough, that x 1) 1/3 1 ) 3 = 8. lim
More informationav 1 x 2 + 4y 2 + xy + 4z 2 = 16.
74 85 Eigenanalysis The subject of eigenanalysis seeks to find a coordinate system, in which the solution to an applied problem has a simple expression Therefore, eigenanalysis might be called the method
More informationThe Residue Theorem. Integration Methods over Closed Curves for Functions with Singularities
The Residue Theorem Integration Methods over losed urves for Functions with Singularities We have shown that if f(z) is analytic inside and on a closed curve, then f(z)dz = 0. We have also seen examples
More informationHandout 1 - Contour Integration
Handout 1 - Contour Integration Will Matern September 19, 214 Abstract The purpose of this handout is to summarize what you need to know to solve the contour integration problems you will see in SBE 3.
More informationKatznelson Problems. Prakash Balachandran Duke University. June 19, 2009
Katznelson Problems Prakash Balachandran Duke University June 9, 9 Chapter. Compute the Fourier coefficients of the following functions (defined by their values on [ π, π)): f(t) { t < t π (t) g(t) { t
More informationComplex Variables Notes for Math 703. Updated Fall Anton R. Schep
Complex Variables Notes for Math 703. Updated Fall 20 Anton R. Schep CHAPTER Holomorphic (or Analytic) Functions. Definitions and elementary properties In complex analysis we study functions f : S C,
More informationPHYS 3900 Homework Set #03
PHYS 3900 Homework Set #03 Part = HWP 3.0 3.04. Due: Mon. Feb. 2, 208, 4:00pm Part 2 = HWP 3.05, 3.06. Due: Mon. Feb. 9, 208, 4:00pm All textbook problems assigned, unless otherwise stated, are from the
More informationQualifying Exam Complex Analysis (Math 530) January 2019
Qualifying Exam Complex Analysis (Math 53) January 219 1. Let D be a domain. A function f : D C is antiholomorphic if for every z D the limit f(z + h) f(z) lim h h exists. Write f(z) = f(x + iy) = u(x,
More informationChapter 14: Vector Calculus
Chapter 14: Vector Calculus Introduction to Vector Functions Section 14.1 Limits, Continuity, Vector Derivatives a. Limit of a Vector Function b. Limit Rules c. Component By Component Limits d. Continuity
More informationChapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =
Chapter 3 The Laplace Transform 3. The Laplace Transform The Laplace transform of the function f(t) is defined L[f(t)] = e st f(t) dt, for all values of s for which the integral exists. The Laplace transform
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationComplex Analysis Homework 1: Solutions
Complex Analysis Fall 007 Homework 1: Solutions 1.1.. a) + i)4 + i) 8 ) + 1 + )i 5 + 14i b) 8 + 6i) 64 6) + 48 + 48)i 8 + 96i c) 1 + ) 1 + i 1 + 1 i) 1 + i)1 i) 1 + i ) 5 ) i 5 4 9 ) + 4 4 15 i ) 15 4
More informationSolutions to practice problems for the final
Solutions to practice problems for the final Holomorphicity, Cauchy-Riemann equations, and Cauchy-Goursat theorem 1. (a) Show that there is a holomorphic function on Ω = {z z > 2} whose derivative is z
More informationPutzer s Algorithm. Norman Lebovitz. September 8, 2016
Putzer s Algorithm Norman Lebovitz September 8, 2016 1 Putzer s algorithm The differential equation dx = Ax, (1) dt where A is an n n matrix of constants, possesses the fundamental matrix solution exp(at),
More information3 + 4i 2 + 3i. 3 4i Fig 1b
The introduction of complex numbers in the 16th century was a natural step in a sequence of extensions of the positive integers, starting with the introduction of negative numbers (to solve equations of
More informationLecture 16 and 17 Application to Evaluation of Real Integrals. R a (f)η(γ; a)
Lecture 16 and 17 Application to Evaluation of Real Integrals Theorem 1 Residue theorem: Let Ω be a simply connected domain and A be an isolated subset of Ω. Suppose f : Ω\A C is a holomorphic function.
More informationHere are brief notes about topics covered in class on complex numbers, focusing on what is not covered in the textbook.
Phys374, Spring 2008, Prof. Ted Jacobson Department of Physics, University of Maryland Complex numbers version 5/21/08 Here are brief notes about topics covered in class on complex numbers, focusing on
More informationINTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES PHILIP FOTH 1. Cauchy s Formula and Cauchy s Theorem 1. Suppose that γ is a piecewise smooth positively ( counterclockwise ) oriented simple closed
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
More informationMAE 105 Introduction to Mathematical Physics HOMEWORK 1. Due on Thursday October 1st 2015
MAE 5 Introduction to Mathematical Physics HOMEWORK Due on Thursday October st 25 PROBEM : Evaluate the following integrals (where n =, 2, 3,... is an integer) and show all your steps: (a) x nπx We use
More informationMTHE 227 Problem Set 2 Solutions
MTHE 7 Problem Set Solutions 1 (Great Circles). The intersection of a sphere with a plane passing through its center is called a great circle. Let Γ be the great circle that is the intersection of the
More informationNeed for transformation?
Z-TRANSFORM In today s class Z-transform Unilateral Z-transform Bilateral Z-transform Region of Convergence Inverse Z-transform Power Series method Partial Fraction method Solution of difference equations
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationThe Laplace Transform
The Laplace Transform Introduction There are two common approaches to the developing and understanding the Laplace transform It can be viewed as a generalization of the CTFT to include some signals with
More informationSTABILITY ANALYSIS. Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated using cones: Stable Neutral Unstable
ECE4510/5510: Feedback Control Systems. 5 1 STABILITY ANALYSIS 5.1: Bounded-input bounded-output (BIBO) stability Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated
More informationFIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS II: Homogeneous Linear Systems with Constant Coefficients
FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS II: Homogeneous Linear Systems with Constant Coefficients David Levermore Department of Mathematics University of Maryland 28 January 2012 Because
More informationEigenpairs and Diagonalizability Math 401, Spring 2010, Professor David Levermore
Eigenpairs and Diagonalizability Math 40, Spring 200, Professor David Levermore Eigenpairs Let A be an n n matrix A number λ possibly complex even when A is real is an eigenvalue of A if there exists a
More informationOnthe Role of the Critical Value Polynomial in Algebraic Optimization
Onthe Role of the Critical Value Polynomial in Algebraic Optimization Bernard Hanzon 1 Joint Work With Andrei Mustata 2 1 Edgeworth Centre for Financial Mathematics, Department of Mathematics, University
More informationSOLUTIONS OF VARIATIONS, PRACTICE TEST 4
SOLUTIONS OF VARIATIONS, PRATIE TEST 4 5-. onsider the following system of linear equations over the real numbers, where x, y and z are variables and b is a real constant. x + y + z = 0 x + 4y + 3z = 0
More informationGetting Some Big Air
Getting Some Big Air Control #10499 February 14, 2011 Abstract In this paper we address the problem of optimizing a ramp for a snowboarder to travel. Our approach is two sided. We first address the forward
More informationFormal power series rings, inverse limits, and I-adic completions of rings
Formal power series rings, inverse limits, and I-adic completions of rings Formal semigroup rings and formal power series rings We next want to explore the notion of a (formal) power series ring in finitely
More information18.03 LECTURE NOTES, SPRING 2014
18.03 LECTURE NOTES, SPRING 2014 BJORN POONEN 7. Complex numbers Complex numbers are expressions of the form x + yi, where x and y are real numbers, and i is a new symbol. Multiplication of complex numbers
More informationME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms
I/O Stability Preliminaries: Vector and function norms 1. Sup norms are used for vectors for simplicity: x = max i x i. Other norms are also okay 2. Induced matrix norms: let A R n n, (i stands for induced)
More informationSPRING 2008: POLYNOMIAL IMAGES OF CIRCLES
18.821 SPRING 28: POLYNOMIAL IMAGES OF CIRCLES JUSTIN CURRY, MICHAEL FORBES, MATTHEW GORDON Abstract. This paper considers the effect of complex polynomial maps on circles of various radii. Several phenomena
More informationMATH 135: COMPLEX NUMBERS
MATH 135: COMPLEX NUMBERS (WINTER, 010) The complex numbers C are important in just about every branch of mathematics. These notes 1 present some basic facts about them. 1. The Complex Plane A complex
More informationControl Systems I. Lecture 9: The Nyquist condition
Control Systems I Lecture 9: The Nyquist condition Readings: Åstrom and Murray, Chapter 9.1 4 www.cds.caltech.edu/~murray/amwiki/index.php/first_edition Jacopo Tani Institute for Dynamic Systems and Control
More informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using
More informationLaplace Transforms and use in Automatic Control
Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
More informationLecture 1 From Continuous-Time to Discrete-Time
Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................
More informationSubdiagonal pivot structures and associated canonical forms under state isometries
Preprints of the 15th IFAC Symposium on System Identification Saint-Malo, France, July 6-8, 29 Subdiagonal pivot structures and associated canonical forms under state isometries Bernard Hanzon Martine
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More informationNotes on Complex Analysis
Michael Papadimitrakis Notes on Complex Analysis Department of Mathematics University of Crete Contents The complex plane.. The complex plane...................................2 Argument and polar representation.........................
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationIntroduction. Performance and Robustness (Chapter 1) Advanced Control Systems Spring / 31
Introduction Classical Control Robust Control u(t) y(t) G u(t) G + y(t) G : nominal model G = G + : plant uncertainty Uncertainty sources : Structured : parametric uncertainty, multimodel uncertainty Unstructured
More informationMATH 311: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE
MATH 3: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE Recall the Residue Theorem: Let be a simple closed loop, traversed counterclockwise. Let f be a function that is analytic on and meromorphic inside. Then
More informationMAT389 Fall 2016, Problem Set 11
MAT389 Fall 216, Problem Set 11 Improper integrals 11.1 In each of the following cases, establish the convergence of the given integral and calculate its value. i) x 2 x 2 + 1) 2 ii) x x 2 + 1)x 2 + 2x
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationComplex varibles:contour integration examples
omple varibles:ontour integration eamples 1 Problem 1 onsider the problem d 2 + 1 If we take the substitution = tan θ then d = sec 2 θdθ, which leads to dθ = π sec 2 θ tan 2 θ + 1 dθ Net we consider the
More informationNovember 18, 2013 ANALYTIC FUNCTIONAL CALCULUS
November 8, 203 ANALYTIC FUNCTIONAL CALCULUS RODICA D. COSTIN Contents. The spectral projection theorem. Functional calculus 2.. The spectral projection theorem for self-adjoint matrices 2.2. The spectral
More informationProblem 1A. Suppose that f is a continuous real function on [0, 1]. Prove that
Problem 1A. Suppose that f is a continuous real function on [, 1]. Prove that lim α α + x α 1 f(x)dx = f(). Solution: This is obvious for f a constant, so by subtracting f() from both sides we can assume
More informationControl Systems I. Lecture 7: Feedback and the Root Locus method. Readings: Jacopo Tani. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 7: Feedback and the Root Locus method Readings: Jacopo Tani Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 2, 2018 J. Tani, E. Frazzoli (ETH) Lecture 7:
More informationMathematical Methods for Engineers and Scientists 1
K.T. Tang Mathematical Methods for Engineers and Scientists 1 Complex Analysis, Determinants and Matrices With 49 Figures and 2 Tables fyj Springer Part I Complex Analysis 1 Complex Numbers 3 1.1 Our Number
More informationLinear and Bilinear Algebra (2WF04) Jan Draisma
Linear and Bilinear Algebra (2WF04) Jan Draisma CHAPTER 3 The minimal polynomial and nilpotent maps 3.1. Minimal polynomial Throughout this chapter, V is a finite-dimensional vector space of dimension
More informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
More informationLeast Squares Based Self-Tuning Control Systems: Supplementary Notes
Least Squares Based Self-Tuning Control Systems: Supplementary Notes S. Garatti Dip. di Elettronica ed Informazione Politecnico di Milano, piazza L. da Vinci 32, 2133, Milan, Italy. Email: simone.garatti@polimi.it
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More informationVectors Coordinate frames 2D implicit curves 2D parametric curves. Graphics 2008/2009, period 1. Lecture 2: vectors, curves, and surfaces
Graphics 2008/2009, period 1 Lecture 2 Vectors, curves, and surfaces Computer graphics example: Pixar (source: http://www.pixar.com) Computer graphics example: Pixar (source: http://www.pixar.com) Computer
More informationControl Systems. Frequency Method Nyquist Analysis.
Frequency Method Nyquist Analysis chibum@seoultech.ac.kr Outline Polar plots Nyquist plots Factors of polar plots PolarNyquist Plots Polar plot: he locus of the magnitude of ω vs. the phase of ω on polar
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationπ 1 = tr(a), π n = ( 1) n det(a). In particular, when n = 2 one has
Eigen Methods Math 246, Spring 2009, Professor David Levermore Eigenpairs Let A be a real n n matrix A number λ possibly complex is an eigenvalue of A if there exists a nonzero vector v possibly complex
More informationIntegration in the Complex Plane (Zill & Wright Chapter 18)
Integration in the omplex Plane Zill & Wright hapter 18) 116-4-: omplex Variables Fall 11 ontents 1 ontour Integrals 1.1 Definition and Properties............................. 1. Evaluation.....................................
More informationMTH 3102 Complex Variables Final Exam May 1, :30pm-5:30pm, Skurla Hall, Room 106
Name (Last name, First name): MTH 02 omplex Variables Final Exam May, 207 :0pm-5:0pm, Skurla Hall, Room 06 Exam Instructions: You have hour & 50 minutes to complete the exam. There are a total of problems.
More informationMath 273 (51) - Final
Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More informationPart IB. Further Analysis. Year
Year 2004 2003 2002 2001 10 2004 2/I/4E Let τ be the topology on N consisting of the empty set and all sets X N such that N \ X is finite. Let σ be the usual topology on R, and let ρ be the topology on
More informationCanonical lossless state-space systems: staircase forms and the Schur algorithm
Canonical lossless state-space systems: staircase forms and the Schur algorithm Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics School of Mathematical Sciences Projet APICS Universiteit
More informationCMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation
CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2014
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2014 1. Please write your 1- or 2-digit exam number on
More informationHomogeneous Constant Matrix Systems, Part II
4 Homogeneous Constant Matrix Systems, Part II Let us now expand our discussions begun in the previous chapter, and consider homogeneous constant matrix systems whose matrices either have complex eigenvalues
More informationChapter 7: The z-transform
Chapter 7: The -Transform ECE352 1 The -Transform - definition Continuous-time systems: e st H(s) y(t) = e st H(s) e st is an eigenfunction of the LTI system h(t), and H(s) is the corresponding eigenvalue.
More information1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =
Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values
More informationELLIPTIC FUNCTIONS AND THETA FUNCTIONS
ELLIPTIC FUNCTIONS AND THETA FUNCTIONS LECTURE NOTES FOR NOV.24, 26 Historically, elliptic functions were first discovered by Niels Henrik Abel as inverse functions of elliptic integrals, and their theory
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Spring Semester 2019
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Spring Semester 2019 1. Please write your 1- or 2-digit exam number
More informationDr. Allen Back. Sep. 8, 2014
in R 3 Dr. Allen Back Sep. 8, 2014 in R 3 in R 3 Def: For f (x, y), the partial derivative with respect to x at p 0 = (x 0, y 0 ) is f x = lim f (x 0 + h, y 0 ) f (x 0, y 0 ) h 0 h or f x = lim f (p 0
More informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationProblem Set 7 Solution Set
Problem Set 7 Solution Set Anthony Varilly Math 3: Complex Analysis, Fall 22 Let P (z be a polynomial Prove there exists a real positive number ɛ with the following property: for all non-zero complex numbers
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationWe are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero
Chapter Limits of Sequences Calculus Student: lim s n = 0 means the s n are getting closer and closer to zero but never gets there. Instructor: ARGHHHHH! Exercise. Think of a better response for the instructor.
More informationCalculus C (ordinary differential equations)
Calculus C (ordinary differential equations) Lesson 9: Matrix exponential of a symmetric matrix Coefficient matrices with a full set of eigenvectors Solving linear ODE s by power series Solutions to linear
More informationMoment of a force (scalar, vector ) Cross product Principle of Moments Couples Force and Couple Systems Simple Distributed Loading
Chapter 4 Moment of a force (scalar, vector ) Cross product Principle of Moments Couples Force and Couple Systems Simple Distributed Loading The moment of a force about a point provides a measure of the
More informationSecond In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 2011
Second In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 31 March 211 (1) [6] Give the interval of definition for the solution of the initial-value problem d 4 y dt 4 + 7 1 t 2 dy dt
More informationSolution of Linear State-space Systems
Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationODEs Cathal Ormond 1
ODEs Cathal Ormond 2 1. Separable ODEs Contents 2. First Order ODEs 3. Linear ODEs 4. 5. 6. Chapter 1 Separable ODEs 1.1 Definition: An ODE An Ordinary Differential Equation (an ODE) is an equation whose
More informationMATH 2433 Homework 1
MATH 433 Homework 1 1. The sequence (a i ) is defined recursively by a 1 = 4 a i+1 = 3a i find a closed formula for a i in terms of i.. In class we showed that the Fibonacci sequence (a i ) defined by
More informationCONSTRAINED PERCOLATION ON Z 2
CONSTRAINED PERCOLATION ON Z 2 ZHONGYANG LI Abstract. We study a constrained percolation process on Z 2, and prove the almost sure nonexistence of infinite clusters and contours for a large class of probability
More informationDigital Control: Part 2. ENGI 7825: Control Systems II Andrew Vardy
Digital Control: Part 2 ENGI 7825: Control Systems II Andrew Vardy Mapping the s-plane onto the z-plane We re almost ready to design a controller for a DT system, however we will have to consider where
More information