Advanced methods for ODEs and DAEs
|
|
- Samuel Holland
- 5 years ago
- Views:
Transcription
1 Lecture : Implicit Runge Kutta method Bojana Roić, 9. April 7
2 What you need to know before thi lecture numerical integration: Lecture from ODE iterative olver: Lecture 5-8 from ODE 9. April 7 Bojana Roić Seite
3 Why not explicit Runge Kutta method? Explicit Runge Kutta method can be uccefully ued for integration of non-tiff problem x = x fail to capture the tiff phenomena x = 6 x even when the tep ize i decreaed dratically. (Plauibilize (Sketch) with RHS a vector field and Euler-FW and BW-Method!) 9. April 7 Bojana Roić Seite
4 General implicit Runge Kutta method Runge Kutta method i preented by Butcher table: c a a c a a c a a b b a a..., a b or c A bt Here, A may be uch that aij =, j > i explicit Runge Kutta method (ERK) uch that aij =, j > i and at leat one aii 6= diagonal implicit Runge Kutta method (DIRK) uch that aij =, j > i and all aii 6= and identical ingly diagonal implicit Runge Kutta method (SDIRK) full matrix fully implicit Runge Kutta method (IRK) 9. April 7 Bojana Roić Seite 4
5 Derivation Given firt order differential equation x = f (t, x) the olution hall be numerically evaluated dx = f (t, x) x(tn ) = x(tn )+ dt Z tn f (τ, x)dτ = x(tn )+h tn uch that the local truncation error atifie maxk loc h k 6 Chp 9. April 7 Bojana Roić Seite 5 i= bi k i
6 Implicit Euler method The implet poible cae i when = and b =, k = f (tn + h, x(tn + h)) i.e., the area under the graph f (τ, x) from tn up to tn i approximated by a right rectangle Z tn f (τ, x)dτ = hf (tn + h, x(tn + h)) = hf (tn, x(tn )) tn uch that xn = xn + hf (tn, xn ). Thu, the local error i loc = xa (tn ) x(tn ) = xa (tn ) xa (tn ) in which xa (tn ) = xa (tn ) + dxa dt tn (tn 9. April 7 Bojana Roić Seite 6 dxa (tn tn ) = O(h ) dt tn tn ) + O(h )
7 Implicit Euler method 9. April 7 Bojana Roić Seite 7
8 Generaliation: collocation baed RKM To numerically evaluate the integral in Z tn x(tn ) = x(tn ) + f (τ, x)dτ, tn let u define the collocation polynomial P(t) of degree N whoe derivative dp dt coincide with the function f at the ditinct collocation point (tn + ci h) for i =,..., and c,..., c being ditinct real number between and. That i, dp (tn + ci h) = f (tn + ci h, P(tn + ci h)) =: ki, dt and P(tn ) = x(tn ), P(tn + h) = x(tn ). 9. April 7 Bojana Roić Seite 8 i =,...,
9 Generaliation: collocation baed RKM Note that dp dt i the polynomial of mot degree characteried by the mot of coefficient. Thee coefficient can be uniquely determined from given point. In Lecture, ODE, we have hown that the polynomial can be written dp dp (tn +τh) = (tn +ci h)`i (τ) = ki `i (τ) dt dt i= i= in term of Lagrange polynomial `j (τ) = Y i=,i6=j τ ci, cj ci τ [, ] with property `i (cj ) = δij in which δij i the Kronecker delta ymbol. 9. April 7 Bojana Roić Seite 9
10 Generaliation: collocation baed RKM One may write that lope at t + hci i that of P, hence dp (t + cj h) = ki `i (cj ) dt kj = and i= Z tn Z x(tn ) = x(tn ) + f (τ, x)dτ = x(tn ) + h tn dp (tn + τh)dτ dt Integrating Z dp (t + τh)dτ = dt Z ki `i (τ)dτ = i= i= Z ki `i (τ)dτ = i= Z bi = `i (τ)dτ, ki = f (tn + ci h, P(tn + ci h)) 9. April 7 Bojana Roić Seite ki bi
11 Generaliation: collocation baed RKM Furthermore, P(tn + ci h) i unknown in ki = f (tn + ci h, P(tn + ci h)) and can be numerically etimated Z ci P(tn + ci h) = P(tn ) + h = x(tn ) + h Z ci dp (tn + τh)dτ dt kj `j (τ)dτ j= = x(tn ) + h j= Z ci kj `j (τ)dτ = x(tn ) + h j= Z ci aij := `j (τ)dτ aij are the RK coeff! 9. April 7 Bojana Roić Seite kj aij
12 Generaliation: collocation baed RKM Altogether: ki = f (tn + ci h, xn + h i =,..., kj aij ), j= x(tn ) = x(tn ) + h bi ki i= Z ci in which aij = `j (τ)dτ Z bi = `i (τ)dτ Thu, for collocation baed RKM one only need to know ci (quadrature rule). 9. April 7 Bojana Roić Seite
13 Propertie of collocation baed RKM Interpolation of time polynomial up to degree hall be error free τq = cjq `j (τ), q =,..., Hence, q =, Z ci τ dτ = ci = Rci j= τ dτ = ci, and alo Z ci j= P j= aij : `j (τ)dτ = Z ci cj `j (τ)dτ = j= invariance to autonomiation! 9. April 7 Bojana Roić Seite Z ci j= `j (τ)dτ = j= aij,
14 Propertie of collocation baed RKM Similarly, one may how aij cjm = j= a well a cim, m i =,...,, m =,... bj cjm i= = Z i= i= Z = ( `j (τ)dτ)cjm i= Z cjm `j (τ)dτ bi cim aij = bj ( cjm ) m 9. April 7 Bojana Roić Seite 4, τm dτ = = m j =,...,, m =,...,
15 Order condition Thee three equalitie can be gathered to B(p) : bi cim =, m aij cjm = cim, m i= C (q) : j= D(r ) : i= bi cim aij = bj k =,..., p b T c m = m i =,...,, m =,...q ( cjm ) m, j =,...,, m =,..., r with r=q = and p =. D commonly derived by comparion with Taylor. 9. April 7 Bojana Roić Seite 5
16 Order condition The condition B(p) i equivalent to a quadrature rule with the node ci and the weight bi which integrate polynomial of degree p exact. The condition C (q) have the following meaning: the intermediate value ki are integrated exactly by a quadrature rule with the weight aij and the node ci. Thi rule integrate polynomial of degree q exactly. The condition D(r ) mean that the olution x i integrated exactly up to the polynomial degree r. 9. April 7 Bojana Roić Seite 6
17 Order condition Theorem An RK-method with internal tage ha the convergence order p, if the implifying condition B(p), C (l), and D(m) with p 6 min{l + m +, l + } are atified. 9. April 7 Bojana Roić Seite 7
18 Lobatto-IIIA method 9. April 7 Bojana Roić Seite
19 Lobatto-IIIB method 9. April 7 Bojana Roić Seite
20 Lobatto-IIIC method 9. April 7 Bojana Roić Seite
21 General RKM A natural generalization of collocation method i obtained by allowing the cofficient ci,bi, and aij, i =,..., to take on arbitrary value, not necearily related to quadrature rule. In fact, one no longer aume ci to be ditinct. The reult i the cla of general Runge-Kutta method. However, thi tory i not currently the cope of thi ubject. 9. April 7 Bojana Roić Seite
22 Implementation Let u implement RKM on the imple general nonlinear ODE x = f (x, t), x() = x Hence, one obtain kj = f (tn +cj h, xn +h ki aij ), i= x(tn ) = x(tn ) + h bi k i i= 9. April 7 Bojana Roić Seite j =,...,
23 Implementation Implicit Runge Kutta method have nice propertie. However, the main drawback i that one obtain the nonlinear ytem of equation k = f (tn + c h, xn + ha k + + ha k )... k = f (tn + c h, xn + ha k + + ha k ) The dimenion of thi ytem i. 9. April 7 Bojana Roić Seite
24 Implementation In pecial cae, for example in the cae of diagonally implicit Runge Kutta method (A i lower triangular including diagonal), thi ytem plit into non-linear decoupled equation: k = f (tn + c h, xn + ha k )... k = f (tn + c h, xn + ha k + + ha k ) 9. April 7 Bojana Roić Seite 4
25 Implementation To olve thee ytem one require nonlinear iterative olver: fixedpoint iteration (Lecture 5 and 6 from lat emeter) Newton or quai-newton method (Lecture 8) 9. April 7 Bojana Roić Seite 5
26 Fixed point iteration To compute unknown k,..., k one may ue fixed point iteration (ν+) (ν) (ν) k = f tn + c h, xn + ha k + + ha k =: f (k,..., k )... (ν+) k (ν) (ν) = f tn + c h, xn + ha k + + ha k =: f (k,..., k ) () () tarting from k,..., k. 9. April 7 Bojana Roić Seite 6
27 Fixed point iteration Gathering all ki into a vector k = (k,..., k )T, one may rewrite previou equation in a matrix form a f (k,..., k ) f (k,..., k ) k(ν+) = F (k(ν) ) =... f (k,..., k ) 9. April 7 Bojana Roić Seite 7
28 Banach fixed point theorem Theorem Let V be a Banach pace, let I V be a cloed ubet, and let F : I I be a contraction ( 6 q < ) uch that kf (k) F (w)k 6 qkk wk for all k, w I. hold. Then the following concluion hold:.) F ha a unique fixed point k I, i. e. F (k ) = k..) For any k I the equence xn+ = F (xn ) converge to x. Furthermore for differentiable function we aid that (Lecture 6, ODE) q := up kf (ξ)k, ξ = w + η(k w ), η [,] 9. April 7 Bojana Roić Seite 8 w 6ξ6k
29 Convergence Hence, in our cae f / k..... F =. f / k... f / k a fx.... = h.. f / k a fx a fx... a fx in which f (tn + ci h, xn + hai k + + haij kj + + hai k ) fi = = haij fx kj kj Thu, for convergence mut hold q := up kf k < 9. April 7 Bojana Roić Seite 9
30 Exercie Let u integrate ODE x = f (t, x) = te x, x = x(.) =, h = by uing Radau IA method: /4 /4 / /4 5/ /4 /4 i.e. k = f (tn + h, xn + h k h k ) = tn e xn +h 4 k h 4 k 4 4 k = f (tn + h, xn + h k + h k ) = (tn + h)e xn +h 4 k +h 4 k April 7 Bojana Roić Seite
31 Exercie and xn+ = xn + h bi ki = xn + h( k + k ) 4 4 In the firt tep n = one ha k = t e x +h 4 k h 4 k =.e + 4 k 4 k k = (t + h)e x +h 4 k +h 4 k = (.7667)e + 4 k + 4 k Thi i nonlinear ytem of equation which can be olved by fixed point iteration (ν) (ν) (ν+) k =.e + 4 k 4 k (ν+) k (ν) = (.7667)e + 4 k () in which the initial etimate k () = k 9. April 7 Bojana Roić Seite (ν) + 4 k = f (t, x ) =.e.
32 Exercie Thu, the firt fixedpoint iteration i () k () k () =.e + 4 k () 4 k () = (.7667)e + 4 k =.78 () + 4 k =.75 The econd i () k () k () 4 k () =.e + 4 k () = (.7667)e + 4 k =.47 () + 4 k = and o on until m th iteration in which (m+) fix = ki (m) ki (m) ki 9. April 7 Bojana Roić Seite 6 tol, i =,
33 Exercie (m) Once the converged k the firt tep value (m) and k are evaluated, one may etimate (m) (m) bi ki = x + h( k + k ) 4 4 Now the whole proce ha to be repeated for the econd time tep: x = x + h bi ki = x + h( k + k 4 4 in which again k and k are unknown. Hence, again one applie fixed point iteration (ν) (ν) (ν+) k =.e + 4 k 4 k x = x + h (ν+) k (ν) = (.7667)e + 4 k () in which the initial etimate k () = k 9. April 7 Bojana Roić Seite (ν) + 4 k = f (t, x ).
34 Newton method To olve nonlinear equation: k = f (tn + c h, xn + ha k + + ha k )... k = f (tn + c h, xn + ha k + + ha k ) by Newton method, let u rewrite it a f (k ) := k f (tn + c h, xn + ha k + + ha k ) =... f (k ) := k f (tn + c h, xn + ha k + + ha k ) = 9. April 7 Bojana Roić Seite 4
35 Newton method Hence, f (k,..., k ) f (k,..., k ) = G (k) =... f (k,..., k ) The olution i obtained a k(i+) = k(i) (Gk ) (k(i) )G (k(i) ) in which the Jacobian i obtained a f / k.. Gk (k) =. f / k f / k ha fx.... =.. f / k ha fx 9. April 7 Bojana Roić Seite ha fx... ha fx
36 Convergence To invetigate the convergence one may oberve the mapping k(i+) = Φ(k(i) ) = k(i) (Gk (k(i) )) G (k(i) ) and apply fixed point theorem q := up kφ k < See lide from Lecture 8, ODE. 9. April 7 Bojana Roić Seite 6
37 Exercie Let u integrate ODE x = f (t, x) = te x, x = x(.) =, h = by uing Radau IA method: /4 /4 / /4 5/ /4 /4 i.e. k = f (tn + h, xn + h k h k ) = tn e xn +h 4 k h 4 k 4 4 k = f (tn + h, xn + h k + h k ) = (tn + h)e xn +h 4 k +h 4 k April 7 Bojana Roić Seite 7
38 Exercie and xn+ = xn + h bi ki = xn + h( k + k ) 4 4 In the firt tep n = one ha k = t e x +h 4 k h 4 k =.e + 4 k 4 k k = (t + h)e x +h 4 k +h 4 k = (.7667)e + 4 k + 4 k Thi i nonlinear ytem of equation which can be rewritten a: f (k, k ) = k.e + 4 k 4 k = f (k, k ) = k (.7667)e + 4 k + 4 k = and olved by Newton method. 9. April 7 Bojana Roić Seite 8
39 Exercie i.e. (i+) k (i+) k! (i) = k (i) k! (G (k(i) ) with G= = (i) f / k f / k () Starting value are again k (i) () 9. April 7 Bojana Roić Seite 9 (i) = k (i) + 4 k 4 k ) (.e k + 4 k + 4 k ) ((.7667)e k = f (t, x ).! 4 k k (.7667)e + 4 k f / k f / k (.e + 4 k 4 k ) k k +k ) ((.7667)e k (i) k.e + 4 k (i) + 4 k
40 Exercie (m) (m) Once the Newton method ha converged to value k and k, one may etimate the firt tep value (m) (m) x = x + h bi ki = x + h( k + k ) 4 4 Now the whole proce ha to be repeated for the econd time tep: x = x + h bi ki = x + h( k + k 4 4 in which again k and k are unknown. Hence, again one applie Newton method to olve k = t e x +h 4 k h 4 k k = (t + h)e x +h 4 k +h 4 k 9. April 7 Bojana Roić Seite 4
41 Simplified Newton method In Newton method one ha to etimate Jacobian, which i more than expenive. In order to avoid thi, one may ue implified Newton method which etimate Jacobian only at the old time tep The olution i obtained a k(i+) = k(i) (J) G (k(i) ) in which the Jacobian i obtained a ha fx (tn, xn ).. J=. ha fx (tn, xn ) ha fx (tn, xn ).... ha fx (tn, xn ) compute Jacobian firt time and keep it contant for each Newton iteration 9. April 7 Bojana Roić Seite 4
Advanced methods for ODEs and DAEs
Advanced methods for ODEs and DAEs Lecture 13: Overview Bojana Rosic, 13. Juli 2016 Dynamical system is a system that evolves over time possibly under external excitations: river flow, car, etc. The dynamics
More informationAdvanced methods for ODEs and DAEs
Advanced methods for ODEs and DAEs Lecture 8: Dirk and Rosenbrock Wanner methods Bojana Rosic, 14 Juni 2016 General implicit Runge Kutta method Runge Kutta method is presented by Butcher table: c1 a11
More informationReformulation of Block Implicit Linear Multistep Method into Runge Kutta Type Method for Initial Value Problem
International Journal of Science and Technology Volume 4 No. 4, April, 05 Reformulation of Block Implicit Linear Multitep Method into Runge Kutta Type Method for Initial Value Problem Muhammad R., Y. A
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationProblem Set 8 Solutions
Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem
More informationMulti-dimensional Fuzzy Euler Approximation
Mathematica Aeterna, Vol 7, 2017, no 2, 163-176 Multi-dimenional Fuzzy Euler Approximation Yangyang Hao College of Mathematic and Information Science Hebei Univerity, Baoding 071002, China hdhyywa@163com
More informationMATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:
MATEMATIK Datum: 20-08-25 Tid: eftermiddag GU, Chalmer Hjälpmedel: inga A.Heintz Telefonvakt: Ander Martinon Tel.: 073-07926. Löningar till tenta i ODE och matematik modellering, MMG5, MVE6. Define what
More informationSolving Differential Equations by the Laplace Transform and by Numerical Methods
36CH_PHCalter_TechMath_95099 3//007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the
More information18.03SC Final Exam = x 2 y ( ) + x This problem concerns the differential equation. dy 2
803SC Final Exam Thi problem concern the differential equation dy = x y ( ) dx Let y = f (x) be the olution with f ( ) = 0 (a) Sketch the iocline for lope, 0, and, and ketch the direction field along them
More informationPractice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions
For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationChapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem
Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationLecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)
Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationA SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho
A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the
More informationMultistage Methods I: Runge-Kutta Methods
Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.
More informationFactor Analysis with Poisson Output
Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu
More informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationLaplace Transformation
Univerity of Technology Electromechanical Department Energy Branch Advance Mathematic Laplace Tranformation nd Cla Lecture 6 Page of 7 Laplace Tranformation Definition Suppoe that f(t) i a piecewie continuou
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationHybrid Control and Switched Systems. Lecture #6 Reachability
Hbrid Control and Switched Stem Lecture #6 Reachabilit João P. Hepanha Univerit of California at Santa Barbara Summar Review of previou lecture Reachabilit tranition tem reachabilit algorithm backward
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationAccuracy of Symmetric Partitioned Runge-Kutta Methods for Differential Equations on Lie-Groups
Accuracy of Symmetric Partitioned Runge-Kutta Method for Differential Equation on Lie-Group WUB / 11-23 BUW-IMACM 11-19, Michael Günther, Franceco Knechtli and Michèle Wandelt Bergiche Univerität Wuppertal,
More informationMarch 18, 2014 Academic Year 2013/14
POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of
More informationThe machines in the exercise work as follows:
Tik-79.148 Spring 2001 Introduction to Theoretical Computer Science Tutorial 9 Solution to Demontration Exercie 4. Contructing a complex Turing machine can be very laboriou. With the help of machine chema
More informationCHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More informationLecture 8: Period Finding: Simon s Problem over Z N
Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationUNIT 15 RELIABILITY EVALUATION OF k-out-of-n AND STANDBY SYSTEMS
UNIT 1 RELIABILITY EVALUATION OF k-out-of-n AND STANDBY SYSTEMS Structure 1.1 Introduction Objective 1.2 Redundancy 1.3 Reliability of k-out-of-n Sytem 1.4 Reliability of Standby Sytem 1. Summary 1.6 Solution/Anwer
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationNONISOTHERMAL OPERATION OF IDEAL REACTORS Plug Flow Reactor
NONISOTHERMAL OPERATION OF IDEAL REACTORS Plug Flow Reactor T o T T o T F o, Q o F T m,q m T m T m T mo Aumption: 1. Homogeneou Sytem 2. Single Reaction 3. Steady State Two type of problem: 1. Given deired
More informationName: Solutions Exam 3
Intruction. Anwer each of the quetion on your own paper. Put your name on each page of your paper. Be ure to how your work o that partial credit can be adequately aeed. Credit will not be given for anwer
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationSome Sets of GCF ϵ Expansions Whose Parameter ϵ Fetch the Marginal Value
Journal of Mathematical Reearch with Application May, 205, Vol 35, No 3, pp 256 262 DOI:03770/jin:2095-26520503002 Http://jmredluteducn Some Set of GCF ϵ Expanion Whoe Parameter ϵ Fetch the Marginal Value
More informationChapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog
Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou
More informationLTV System Modelling
Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction
More informationDigital Control System
Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)
More informationNotes on Phase Space Fall 2007, Physics 233B, Hitoshi Murayama
Note on Phae Space Fall 007, Phyic 33B, Hitohi Murayama Two-Body Phae Space The two-body phae i the bai of computing higher body phae pace. We compute it in the ret frame of the two-body ytem, P p + p
More informationCHE302 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHE3 ECTURE V APACE TRANSFORM AND TRANSFER FUNCTION Profeor Dae Ryook Yang Fall Dept. of Chemical and Biological Engineering Korea Univerity CHE3 Proce Dynamic and Control Korea Univerity 5- SOUTION OF
More informationNew bounds for Morse clusters
New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität
More informationmultistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):
MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem
More informationCorrection for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002
Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in
More informationControl Systems Analysis and Design by the Root-Locus Method
6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If
More informationSolutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam
BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning
More informationMA 266 FINAL EXAM INSTRUCTIONS May 2, 2005
MA 66 FINAL EXAM INSTRUCTIONS May, 5 NAME INSTRUCTOR. You mut ue a # pencil on the mark ene heet anwer heet.. If the cover of your quetion booklet i GREEN, write in the TEST/QUIZ NUMBER boxe and blacken
More informationTHE DIVERGENCE-FREE JACOBIAN CONJECTURE IN DIMENSION TWO
THE DIVERGENCE-FREE JACOBIAN CONJECTURE IN DIMENSION TWO J. W. NEUBERGER Abtract. A pecial cae, called the divergence-free cae, of the Jacobian Conjecture in dimenion two i proved. Thi note outline an
More informationL 2 -transforms for boundary value problems
Computational Method for Differential Equation http://cmde.tabrizu.ac.ir Vol. 6, No., 8, pp. 76-85 L -tranform for boundary value problem Arman Aghili Department of applied mathematic, faculty of mathematical
More informationManprit Kaur and Arun Kumar
CUBIC X-SPLINE INTERPOLATORY FUNCTIONS Manprit Kaur and Arun Kumar manpreet2410@gmail.com, arun04@rediffmail.com Department of Mathematic and Computer Science, R. D. Univerity, Jabalpur, INDIA. Abtract:
More informationOptimal Coordination of Samples in Business Surveys
Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationMoment of Inertia of an Equilateral Triangle with Pivot at one Vertex
oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.
More informationAdvanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis
Advanced Digital ignal Proceing Prof. Nizamettin AYDIN naydin@yildiz.edu.tr Time-Frequency Analyi http://www.yildiz.edu.tr/~naydin 2 tationary/nontationary ignal Time-Frequency Analyi Fourier Tranform
More informationSampling and the Discrete Fourier Transform
Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at
More informationSTOCHASTIC DIFFERENTIAL GAMES:THE LINEAR QUADRATIC ZERO SUM CASE
Sankhyā : he Indian Journal of Statitic 1995, Volume 57, Serie A, Pt. 1, pp.161 165 SOCHASIC DIFFERENIAL GAMES:HE LINEAR QUADRAIC ZERO SUM CASE By R. ARDANUY Univeridad de Salamanca SUMMARY. hi paper conider
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar
DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:
More informationHyperbolic Partial Differential Equations
Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a
More informationProblem 1. Construct a filtered probability space on which a Brownian motion W and an adapted process X are defined and such that
Stochatic Calculu Example heet 4 - Lent 5 Michael Tehranchi Problem. Contruct a filtered probability pace on which a Brownian motion W and an adapted proce X are defined and uch that dx t = X t t dt +
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationList coloring hypergraphs
Lit coloring hypergraph Penny Haxell Jacque Vertraete Department of Combinatoric and Optimization Univerity of Waterloo Waterloo, Ontario, Canada pehaxell@uwaterloo.ca Department of Mathematic Univerity
More informationAdvanced methods for ODEs and DAEs. Lecture 9: Multistep methods/ Differential algebraic equations
Advanced methods for ODEs and DAEs Lecture 9: Multistep methods/ Differential algebraic equations Bojana Rosic, 22. Juni 2016 PART I: MULTISTEP METHODS 22. Juni 2016 Bojana Rosić Advanced methods for ODEs
More informationOnline Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat
Online Appendix for Managerial Attention and Worker Performance by Marina Halac and Andrea Prat Thi Online Appendix contain the proof of our reult for the undicounted limit dicued in Section 2 of the paper,
More informationChapter Landscape of an Optimization Problem. Local Search. Coping With NP-Hardness. Gradient Descent: Vertex Cover
Coping With NP-Hardne Chapter 12 Local Search Q Suppoe I need to olve an NP-hard problem What hould I do? A Theory ay you're unlikely to find poly-time algorithm Mut acrifice one of three deired feature
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More informationUnbounded solutions of second order discrete BVPs on infinite intervals
Available online at www.tjna.com J. Nonlinear Sci. Appl. 9 206), 357 369 Reearch Article Unbounded olution of econd order dicrete BVP on infinite interval Hairong Lian a,, Jingwu Li a, Ravi P Agarwal b
More informationSimulation Study on the Shock Properties of the Double-Degree-of-Freedom Cushioning Packaging System
Proceeding of the 7th IAPRI World Conference on Packaging Simulation Study on the Shock Propertie of the Double-Degree-of-Freedom Cuhioning Packaging Sytem Xia Zhu, Qiaoqiao Yan, Xiaoling Yao, Junbin Chen,
More informationLaplace Adomian Decomposition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernels
Studie in Nonlinear Science (4): 9-4, ISSN -9 IDOSI Publication, Laplace Adomian Decompoition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernel F.A. Hendi Department of Mathematic
More informationTRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL
GLASNIK MATEMATIČKI Vol. 38583, 73 84 TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL p-laplacian Haihen Lü, Donal O Regan and Ravi P. Agarwal Academy of Mathematic and Sytem Science, Beijing, China, National
More informationLecture 9: Shor s Algorithm
Quantum Computation (CMU 8-859BB, Fall 05) Lecture 9: Shor Algorithm October 7, 05 Lecturer: Ryan O Donnell Scribe: Sidhanth Mohanty Overview Let u recall the period finding problem that wa et up a a function
More information1 Basic Equations of the PLLs
1 Baic Equation of the PLL 1.1 INTRODUCTION Phae lock loop (PLL) belong to a larger et of regulation ytem. A an independent reearch and deign field it tarted in the 1950 [1] and gained major practical
More informationNAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE
POLITONG SHANGHAI BASIC AUTOMATIC CONTROL June Academic Year / Exam grade NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE Ue only thee page (including the bac) for anwer. Do not ue additional
More informationDiagonally Implicit Symplectic Runge-Kutta Methods with Special Properties
Appl. Math. Inf. Sci. 9, No. 1L, 11-17 (2015) 11 Applied Mathematic & Information Science An International Journal http://dx.doi.org/10.12785/ami/091l02 Diagonally Implicit Symplectic Runge-Kutta Method
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More informationCHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang
CHBE3 ECTURE V APACE TRANSFORM AND TRANSFER FUNCTION Profeor Dae Ryook Yang Spring 8 Dept. of Chemical and Biological Engineering 5- Road Map of the ecture V aplace Tranform and Tranfer function Definition
More informationMath 201 Lecture 17: Discontinuous and Periodic Functions
Math 2 Lecture 7: Dicontinuou and Periodic Function Feb. 5, 22 Many example here are taken from the textbook. he firt number in () refer to the problem number in the UA Cutom edition, the econd number
More informationCS 170: Midterm Exam II University of California at Berkeley Department of Electrical Engineering and Computer Sciences Computer Science Division
1 1 April 000 Demmel / Shewchuk CS 170: Midterm Exam II Univerity of California at Berkeley Department of Electrical Engineering and Computer Science Computer Science Diviion hi i a cloed book, cloed calculator,
More informationExercise set 1 Solution
Exercie et 1 Solution 1. What are the main poible advantage of DNA-baed computing? How about diadvantage? Advantage: - More parallel computation - Able to tore tremendou amount of information in very mall
More informationEnergy-diminishing integration of gradient systems
Verion of 16 May 213 Energy-diminihing integration of gradient ytem ERNST HAIRER Section de mathématique, Univerité de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland AND CHRISTIAN LUBICH Mathematiche
More informationEfficient Methods of Doppler Processing for Coexisting Land and Weather Clutter
Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr
More informationSTABILITY OF A LINEAR INTEGRO-DIFFERENTIAL EQUATION OF FIRST ORDER WITH VARIABLE DELAYS
Bulletin of Mathematical Analyi and Application ISSN: 1821-1291, URL: http://bmathaa.org Volume 1 Iue 2(218), Page 19-3. STABILITY OF A LINEAR INTEGRO-DIFFERENTIAL EQUATION OF FIRST ORDER WITH VARIABLE
More informationRunge-Kutta and Collocation Methods Florian Landis
Runge-Kutta and Collocation Methods Florian Landis Geometrical Numetric Integration p.1 Overview Define Runge-Kutta methods. Introduce collocation methods. Identify collocation methods as Runge-Kutta methods.
More informationWeber Schafheitlin-type integrals with exponent 1
Integral Tranform and Special Function Vol., No., February 9, 47 53 Weber Schafheitlin-type integral with exponent Johanne Kellendonk* and Serge Richard Univerité de Lyon, Univerité Lyon, Intitut Camille
More informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More informationLecture 13. Thermodynamic Potentials (Ch. 5)
Lecture 13. hermodynamic Potential (Ch. 5) So far we have been uing the total internal energy U and ometime the enthalpy H to characterize variou macrocopic ytem. hee function are called the thermodynamic
More informationSOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU. I will collect my solutions to some of the exercises in this book in this document.
SOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU CİHAN BAHRAN I will collect my olution to ome of the exercie in thi book in thi document. Section 2.1 1. Let A = k[[t ]] be the ring of
More informationSOME RESULTS ON INFINITE POWER TOWERS
NNTDM 16 2010) 3, 18-24 SOME RESULTS ON INFINITE POWER TOWERS Mladen Vailev - Miana 5, V. Hugo Str., Sofia 1124, Bulgaria E-mail:miana@abv.bg Abtract To my friend Kratyu Gumnerov In the paper the infinite
More informationPHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS B - HW #6 Spring 4, Solution by David Pace Any referenced equation are from Griffith Problem tatement are paraphraed. Problem. from Griffith Show that the following, A µo ɛ o A V + A ρ ɛ o Eq..4 A
More informationEE Control Systems LECTURE 6
Copyright FL Lewi 999 All right reerved EE - Control Sytem LECTURE 6 Updated: Sunday, February, 999 BLOCK DIAGRAM AND MASON'S FORMULA A linear time-invariant (LTI) ytem can be repreented in many way, including:
More informationLINEAR ALGEBRA METHOD IN COMBINATORICS. Theorem 1.1 (Oddtown theorem). In a town of n citizens, no more than n clubs can be formed under the rules
LINEAR ALGEBRA METHOD IN COMBINATORICS 1 Warming-up example Theorem 11 (Oddtown theorem) In a town of n citizen, no more tha club can be formed under the rule each club have an odd number of member each
More informationMODERN CONTROL SYSTEMS
MODERN CONTROL SYSTEMS Lecture 1 Root Locu Emam Fathy Department of Electrical and Control Engineering email: emfmz@aat.edu http://www.aat.edu/cv.php?dip_unit=346&er=68525 1 Introduction What i root locu?
More information