Numerical Analysis II. Problem Sheet 8
|
|
- Magdalen Hall
- 5 years ago
- Views:
Transcription
1 P. Grohs S. Hosseini Ž. Kereta Spring Term 15 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 8 Problem 8.1 The dierential equation Autonomisation and Strang-Splitting ẏ = A(t)y, y(t ) = y, A(t) = ta, A = can be autonomised by introducing a variable z := ż = F (z) = }{{} =:(z) [ ] 1 1 [ ] y, which yields the ODE t (8.1.1) [ ] A(t)y. (8.1.) 1 (8.1a) Formulate the discrete evolution o a Strang splitting method or (8.1.) corresponding to the decomposition [ ] [ ] A(t)y F (z) = + 1 }{{} =:g(z) based on the exact evolutions or and g. HINT: Use the representation o the exact solution o an IVP or a homogeneous linear dierential equation. Solution: The exact evolutions to the linear homogeneous dierential equations ż = (z) and ż = g(z) are given by [ ] [ ] [ ] [ ] y y y e Φ h = and Φ h g = ha(t) y. t t + h t t The discrete evolutions o the corresponding Strang splitting method is now given by Ψ h Fz = Φ h/ Φ h g Φ h/ z. We thus obtain [ ] [ ] [ ] y Ψ h F = Φ h/ t Φ h g Φ h/ y = Φ h/ y t Φ h g t + 1h [ ] [ ] = Φ h/ e ha(t+ 1 h) y e ha(t+ 1 t + 1h = h) y. t + h Problem Sheet 8 Page 1 Problem 8.1
2 (8.1b) Implement the resulting method by completing the MATLAB template [t, y] = strangaut(y, t, T, h) HINT: The matrix exponential unction can be computed using expm. Solution: See strangaut.m, Listing 8.1. Listing 8.1: Implementation o strangaut.m. 1 u n c t i o n [t, y] = strangaut(y, t, T, h) % Strang splitting method or the solution o an autonomized IVP 3 % 4 % Input: 5 % y: initial value 6 % t: initial time 7 % T: end time 8 % h: step size 9 % 1 % output: 11 % t: time grid 1 % y: solution at the points o the grid % initialization 15 A = [, 1; -1, ]; % adapt step size, so that end time is attained precisely 18 N = c e i l ((T-t) / h); 19 h = (T-t) / N; 1 % allocate memory or results t = z e r o s (N+1, 1); 3 y = z e r o s (, N+1); 4 5 % set starting values 6 t(1) = t; 7 y(:, 1) = y; 8 9 % Timestepping 3 o r k = 1:N 31 3 % apply Strang splitting method y = expm(h * (t +.5*h) * A) * y; 35 t = t + h; 36 % save results 37 t(k+1) = t; 38 y(:, k+1) = y; Problem Sheet 8 Page Problem 8.1
3 39 end (8.1c) Plot the numerical solution obtained by applying the splitting method to the IVP (8.1.1) where y = [1, ], t =, T = 1, h =.1. To do so, complete the template strangautplot.m. Use time t as your z-coordinate. HINT: You might ind the unction plot3 to be o use. Solution: See strangautplot.m, Listing 8.. The resulting plot is shown in Figure 8.1. Listing 8.: Implementation o strangautplot.m. 1 u n c t i o n strangautplot % plot results rom strangaut 3 4 % initialization 5 t = ; 6 T = 1; 7 h =.1; 8 y = [1; ]; 9 1 % compute numerical solution 11 [t, y] = strangaut(y, t, T, h); 1 13 % plot numerical solution 14 i g u r e 15 p l o t 3 (y(1, :), y(, :), t, b- ) 16 x l a b e l ( y_1 ) 17 y l a b e l ( y_ ) 18 z l a b e l ( t, Rotation, ) 19 t i t l e ( numerical solution o the autonomized IVP via the strang splitting method ) Problem 8. Strang-Splitting For continuous dierentiable unction g : D R d R d consider the autonomous dierential equation ( ) ( ) ( ) u u g(v) ẏ = (y) with y =, (y) = :=. (8..1) v v g(u) (8.a) Give the method equations o the Strang-splitting 1-step method [NUMODE, Eq. (.5.3)] or (8..1) which is based on the decomposition ( ) ( ) g(v) (y) = (y) + 1 (y) := +. (8..) g(u) Solution: Exact evolution o the problems: ( ) ( u u + tg(v) Φ t 1 = v v ), Φ t ( ) ( u = v u v + tg(u) Problem Sheet 8 Page 3 Problem 8. ).
4 numerical solution o the autonomized IVP via the strang splitting method 1 8 t y y Figure 8.1: Numerical solution o the autonomized linear IVP using the Strang splitting method. Plugging this into the ormula [NUMODE, Eq. (.5.3)] rom the lecture gives the discrete evolution o the Strang-splitting 1-step method ( ) ( ) ( u Ψ h = Φ h/ Φ h 1 Φ h/ u u + hg(v + 1 = hg(u)) ) v v v + 1hg(u) + 1hg(u + hg(v + 1hg(u))). (8..3) (8.b) Is the 1-step method rom subproblem (8.a) an explicit or an implicit method? Solution: The method is explicit, because Ψ t can be evaluted without solving an equation. (8.c) What is the minimal order o this method in the general case? Solution: Order by [NUMODE, Thm..5.5] rom the lecture. (8.d) Implement a MATLAB-unction u = sped(g,u,t,n), that realizes N equidistant iterations o the Strang-splitting 1-step method or (8..1) with initial value ( u ) u on the time interval [, T ] using a minimal number o g( ) evaluations. The unction argument g is a The return value is a d (N + 1) matrix o the values o the u-component o the approximation y k, k =,..., N. Solution: Let ṽ := v + h g(u) then (8..3), turns into ( ) ( ) u1 u + hg(ṽ = ), (8..4) ṽ + hg(u 1 ) ṽ 1 Problem Sheet 8 Page 4 Problem 8.
5 which requires N evaluations o the unction g. See Listing 8.3. Listing 8.3: subproblem (8.d) 1 u n c t i o n u = sped(g,u,t,n) % unction sped (Strang Splitting) 3 % g : handle to right hand side unction 4 % u: initial value 5 % T : inal time 6 % N : number o timesteps 7 % return value: values or u-components h = T/N; 11 u = z e r o s ( l e n g t h (u),n+1); 1 u(:,1) = u; 13 v = u +.5*h*g(u(:,1)); o r k=1:n-1 16 u(:,k+1) = u(:,k) + h*g(v); 17 v = v + h*g(u(:,k+1)); 18 end 19 u(:,n+1) = u(:,n) + h*g(v); (8.e) Write a MATLAB-unction spedcirc.m that numerically solves the equation o motion or circular motion ẋ = y, ẏ = x or x() = 1, y() = using the method sped rom subproblem (8.d) with T = 6, N =. It should also plot the solution. Solution: In order to solve the circular motion using the method sped rom subproblem (8.d) we set u = (x, y) and v = (x, y), see Listing 8.4. For a plot o the solution see Figure 8.. Listing 8.4: subproblem (8.e) 1 % MATLAB script: circular movement with sped 3 % right hand side 4 g [-y();y(1)]; 5 6 % initial value 7 y = [1;]; 8 9 % Final time 1 T = 6; 11 1 % Number o timesteps Problem Sheet 8 Page 5 Problem 8.
6 sped: T = 6., N = y Figure 8.: Approximated trajectory o the solution in subproblem (8.e). x 13 N = ; y = sped(g,y,t,n); % Plot trajectory i g u r e ( name, spedcirc ); p l o t (y(1,:),y(,:), r. ); hold on; 1 p l o t (cos(:.1:*pi), s i n (:.1:*pi), g- ); p l o t ([ ],[ ], k- ); 3 p l o t ([ ],[ ], k- ); 4 x l a b e l ( {\b x}, ontsize,14); 5 y l a b e l ( {\b y}, ontsize,14); 6 t i t l e ( s p r i n t ( sped: T = %, N = %d,t,n)); 7 a x i s ([ ]); 8 9 p r i n t -depsc spedcirc.eps ; 3 31 i g u r e ( name, length deviation ); 3 p l o t ((:T/N:T), (y(1,:).ˆ+y(,:).ˆ)-1, r+ ); hold on; 34 x l a b e l ( {\b t}, ontsize,14); 35 y l a b e l ( {\b u_k ˆ-1}, ontsize,14); 36 t i t l e ( s p r i n t ( sped: T = %, N = %d,t,n)); p r i n t -depsc spedcirclen.eps ; Problem Sheet 8 Page 6 Problem 8.
7 (8.) Investigate numerically i the approximations u k rom sped in subproblem (8.e) are on a circle. Solution: The length is not preserved, as can be seen in Figure 8.3. One can also show this by perorming one step o the method with sped and subtracts the norm o the resulting vector by the norm o the initial vector ( 1 )..5 sped: T = 6., N =..15 u k Figure 8.3: Length variation in the numerical experiment rom subproblem (8.) t Problem 8.3 Stability o a Decomposed Method For the autonomous ODE ẏ = F(y) we will consider the ollowing decomposition o the right hand side, which we assume to be suiciently smooth: F(y) = DF(y )(y y ) + F(y) DF(y )(y y ). (8.3.1) }{{}}{{} =:(y) =:g(y) Taking y = y we deine the discrete evolution y 1 = Ψ h y o a single-step method as ollows: Compute ỹ by perorming one step o the explicit trapezoid rule with step-size h the initial value problem ẏ = g(y), y() = y. applied to Compute ŷ as the exact solution o the initial value problem ẏ = (y), y() = ỹ at the point t = h. Compute y 1 by applying the explicit trapezoid rule with step-size h problem ẏ = g(y), y() = ŷ. to the initial vlaue Problem Sheet 8 Page 7 Problem 8.3
8 (8.3a) Determine the stability unction o the method described above. HINT: To check your results, use the MATLAB unction S=stabn(z) (it is given as a pcode). It evaluates the stability unction at a given point z. Solution: To ind the stability unction, we apply the method to the model problem ẏ = λy =: F (y), y() = y. (8.3.) Because DF (y) = λ, the corresponding decomposition is given by g(y) = λy and (y) = λ(y y ). Thus, we receive ỹ = y + h 4 [g(y ) + g(y + h g(y ))] = y + h 4 [λy + λy ] = (1 + 1 λh)y ẏ = λ(y y ) d dt (y y ) = λ(y y ) y(h) y = (y() y )e λh ŷ = y + (ỹ y )e λh = (1 + 1 λheλh )y y 1 = ŷ + h 4 [g(ŷ) + g(ŷ + h g(ŷ))] = ŷ + h 4 [λy + λy ] = (1 + 1 λh + 1 λheλh )y The method s stability unction is now given as S(z) = z + 1 zez. (8.3b) Complete the MATLAB template plotstabdom which plots the boundary o the stability domain. Save the resulting plot as stabdom.eps. HINT: Use the MATLAB unction contour to draw the boundary. Take 5 and 5 as bounds or both axes. Solution: See Listing u n c t i o n plotstabdom % PLOTSTABDOM Listing 8.5: plotstabdom.m 3 % plots the boundary o the stability domain 4 % o the given method 5 % 6 % Input: 7 % none 8 % 9 % Output: 1 % none 11 1 % deine stability unction 13 S 1 + (1 + exp(z)).* z/; % create grid 16 [X, Y] = meshgrid( l i n s p a c e (-5, 5, 1)); % evaluate stability unction on grid 19 Z = abs(s(x + 1i*Y)); Problem Sheet 8 Page 8 Problem 8.3
9 Figure 8.4: Method s Stability Function 1 % plot stability domain contour(x, Y, Z, [1 1]); 3 4 % save plot 5 p r i n t -depsc stabdom.eps The stability domain is depicted in Figure 8.4 (8.3c) Show that the method described above has a convergence rate o at least. HINT: Notice that the previously described method can be written as a well known (and well studied) method. A theorem rom the lectures will then yield the statement. Solution: Let Φ and Φ g respectively be the exact evolutions to the autonomous dierential equations with right hand side and g respectively. Further, let Ψ g be the discrete evolution or ẏ = g(y) that is calculated by the explicit trapezoid rule. With this notation, the (discrete) evolution Ψ h o the given method can be written as Ψ h = Ψ h g Φ h Ψ h g. In other words, it is a Strang splitting method, in which the exact evolution Φ h g is replaced by Ψ h g. As we know, Strang splitting methods as well as the explicit trapezoid rule have consistency order, meaning, Φ h y = (Φ h g Φ h Φ h g )y + O(h 3 ), (8.3.3) Φ h gy = Ψ h gy + O(h 3 ) (8.3.4) holds, where Φ h denotes the exact evolution or the autonomous dierential equation ẏ = F(y). Problem Sheet 8 Page 9 Problem 8.3
10 As it ollows that Φ h i (y + y) = Φ h i y + O( y ), i {, g}, (8.3.5) Φ h y (8.3.3) = Φ h g ( Φ h (Φ h g y) ) + O(h 3 ) (8.3.4) = Φ h g ( Φ h (Ψ h g y + O(h 3 )) ) + O(h 3 ) (8.3.5) = Φ h g ( Φ h (Ψ h g y) + O(h 3 ) ) + O(h 3 ) (8.3.5) = Φ h g ( Φ h (Ψ h g y) ) + O(h 3 ) (8.3.4) = Ψ h g ( Φ h (Ψ h g y) ) + O(h 3 ) = Ψ h y + O(h 3 ). The consistency and thereby the convergence order o the method is thus. Published on 14 April 15. To be submitted by 1 April 15. Reerences [NUMODE] Lecture Slides or the course Numerical Methods or Ordinary Dierential Equations, SVN revision # Last modiied on April 1, 15 Problem Sheet 8 Page 1 Reerences
Numerical Analysis II. Problem Sheet 12
P. Grohs S. Hosseini Ž. Kereta Spring Term 2015 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 12 Problem 12.1 The Exponential Runge-Kutta Method Consider the exponential Runge Kutta single-step
More informationNumerical Analysis II. Problem Sheet 9
H. Ammari W. Wu S. Yu Spring Term 08 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 9 Problem 9. Extrapolation of the Implicit Mid-Point Rule. Taking the implicit mid-point rule as a basis method,
More informationTelescoping Decomposition Method for Solving First Order Nonlinear Differential Equations
Telescoping Decomposition Method or Solving First Order Nonlinear Dierential Equations 1 Mohammed Al-Reai 2 Maysem Abu-Dalu 3 Ahmed Al-Rawashdeh Abstract The Telescoping Decomposition Method TDM is a new
More informationNumerical Methods for CSE. Examination - Solutions
R. Hiptmair A. Moiola Fall term 2010 Numerical Methods for CSE ETH Zürich D-MATH Examination - Solutions February 10th, 2011 Total points: 65 = 8+14+16+16+11 (These are not master solutions but just a
More informationNumerical Methods for CSE. Examination Solutions
R. Hiptmair A. Moiola Fall term 010 Numerical Methods for CSE ETH Zürich D-MATH Examination Solutions August 11th, 011 Total points: 90 = 16+5+16+13+0 Problem 1. Structured matrix-vector product [16=+6+++4
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationDynamical Systems. August 13, 2013
Dynamical Systems Joshua Wilde, revised by Isabel Tecu, Takeshi Suzuki and María José Boccardi August 13, 2013 Dynamical Systems are systems, described by one or more equations, that evolve over time.
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationNumerical Methods for CSE. Examination
R. Hiptmair A. Moiola Fall term 2010 Numerical Methods for CSE ETH Zürich D-MATH Examination August 11th, 2011 Instructions. Duration of examination: 180 minutes. Total points: 90. Concise answers are
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationScattering of Solitons of Modified KdV Equation with Self-consistent Sources
Commun. Theor. Phys. Beijing, China 49 8 pp. 89 84 c Chinese Physical Society Vol. 49, No. 4, April 5, 8 Scattering o Solitons o Modiied KdV Equation with Sel-consistent Sources ZHANG Da-Jun and WU Hua
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationOn linear and non-linear equations.(sect. 2.4).
On linear and non-linear equations.sect. 2.4). Review: Linear differential equations. Non-linear differential equations. Properties of solutions to non-linear ODE. The Bernoulli equation. Review: Linear
More informationOn linear and non-linear equations. (Sect. 1.6).
On linear and non-linear equations. (Sect. 1.6). Review: Linear differential equations. Non-linear differential equations. The Picard-Lindelöf Theorem. Properties of solutions to non-linear ODE. The Proof
More informationMEAN VALUE THEOREM. Section 3.2 Calculus AP/Dual, Revised /30/2018 1:16 AM 3.2: Mean Value Theorem 1
MEAN VALUE THEOREM Section 3. Calculus AP/Dual, Revised 017 viet.dang@humbleisd.net 7/30/018 1:16 AM 3.: Mean Value Theorem 1 ACTIVITY A. Draw a curve (x) on a separate sheet o paper within a deined closed
More informationLesson 15: Oregonator, Chemical Reactions and ode15s
Lesson 15: Oregonator, Chemical Reactions and ode15s 15.1 Applied Problem. The chemical reaction with bromous acid, bromide ion and cerium ion exhibits a remarkable chemical attributes, which oscillate
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear and
More informationSolutions to Math 53 First Exam April 20, 2010
Solutions to Math 53 First Exam April 0, 00. (5 points) Match the direction fields below with their differential equations. Also indicate which two equations do not have matches. No justification is necessary.
More informationReceived: 30 July 2017; Accepted: 29 September 2017; Published: 8 October 2017
mathematics Article Least-Squares Solution o Linear Dierential Equations Daniele Mortari ID Aerospace Engineering, Texas A&M University, College Station, TX 77843, USA; mortari@tamu.edu; Tel.: +1-979-845-734
More informationThu June 16 Lecture Notes: Lattice Exercises I
Thu June 6 ecture Notes: attice Exercises I T. Satogata: June USPAS Accelerator Physics Most o these notes ollow the treatment in the class text, Conte and MacKay, Chapter 6 on attice Exercises. The portions
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationVector Fields and Solutions to Ordinary Differential Equations using MATLAB/Octave
Vector Fields and Solutions to Ordinary Differential Equations using MATLAB/Octave Andreas Stahel 5th December 27 Contents Vector field for the logistic equation 2 Solutions of ordinary differential equations
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationSection 1.4 Circles. Objective #1: Writing the Equation of a Circle in Standard Form.
1 Section 1. Circles Objective #1: Writing the Equation of a Circle in Standard Form. We begin by giving a definition of a circle: Definition: A Circle is the set of all points that are equidistant from
More information1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal
. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More informationResearch Article Diagonally Implicit Block Backward Differentiation Formulas for Solving Ordinary Differential Equations
International Mathematics and Mathematical Sciences Volume 212, Article ID 767328, 8 pages doi:1.1155/212/767328 Research Article Diagonally Implicit Block Backward Differentiation Formulas for Solving
More informationLast name Department Computer Name Legi No. Date Total. Always keep your Legi visible on the table.
First name Last name Department Computer Name Legi No. Date 02.02.2016 1 2 3 4 5 Total Grade Fill in this cover sheet first. Always keep your Legi visible on the table. Keep your phones, tablets and computers
More informationLecture : Feedback Linearization
ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical
More informationMath Review and Lessons in Calculus
Math Review and Lessons in Calculus Agenda Rules o Eponents Functions Inverses Limits Calculus Rules o Eponents 0 Zero Eponent Rule a * b ab Product Rule * 3 5 a / b a-b Quotient Rule 5 / 3 -a / a Negative
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationMath 2412 Activity 1(Due by EOC Sep. 17)
Math 4 Activity (Due by EOC Sep. 7) Determine whether each relation is a unction.(indicate why or why not.) Find the domain and range o each relation.. 4,5, 6,7, 8,8. 5,6, 5,7, 6,6, 6,7 Determine whether
More informationMathematical Notation Math Calculus & Analytic Geometry III
Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and should be emailed to the instructor at james@richland.edu.
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationExample: When describing where a function is increasing, decreasing or constant we use the x- axis values.
Business Calculus Lecture Notes (also Calculus With Applications or Business Math II) Chapter 3 Applications o Derivatives 31 Increasing and Decreasing Functions Inormal Deinition: A unction is increasing
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationOrdinary Differential Equations
Ordinary Differential Equations Introduction: first order ODE We are given a function f(t,y) which describes a direction field in the (t,y) plane an initial point (t 0,y 0 ) We want to find a function
More informationPractice Midterm Solutions
Practice Midterm Solutions Math 4B: Ordinary Differential Equations Winter 20 University of California, Santa Barbara TA: Victoria Kala DO NOT LOOK AT THESE SOLUTIONS UNTIL YOU HAVE ATTEMPTED EVERY PROBLEM
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationChapter 6: Functions with severable variables and Partial Derivatives:
Chapter 6: Functions with severable variables and Partial Derivatives: Functions o several variables: A unction involving more than one variable is called unction with severable variables. Eamples: y (,
More informationMatrix Solutions to Linear Systems of ODEs
Matrix Solutions to Linear Systems of ODEs James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 3, 216 Outline 1 Symmetric Systems of
More informationNON-AUTONOMOUS INHOMOGENEOUS BOUNDARY CAUCHY PROBLEMS AND RETARDED EQUATIONS. M. Filali and M. Moussi
Electronic Journal: Southwest Journal o Pure and Applied Mathematics Internet: http://rattler.cameron.edu/swjpam.html ISSN 83-464 Issue 2, December, 23, pp. 26 35. Submitted: December 24, 22. Published:
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More information0,0 B 5,0 C 0, 4 3,5. y x. Recitation Worksheet 1A. 1. Plot these points in the xy plane: A
Math 13 Recitation Worksheet 1A 1 Plot these points in the y plane: A 0,0 B 5,0 C 0, 4 D 3,5 Without using a calculator, sketch a graph o each o these in the y plane: A y B 3 Consider the unction a Evaluate
More informationHomogeneous Equations with Constant Coefficients
Homogeneous Equations with Constant Coefficients MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 General Second Order ODE Second order ODEs have the form
More informationDON T PANIC! If you get stuck, take a deep breath and go on to the next question. Come back to the question you left if you have time at the end.
Math 307A, Midterm 1 Spring 2013 Name: Instructions. DON T PANIC! If you get stuck, take a deep breath and go on to the next question. Come back to the question you left if you have time at the end. There
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationLecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material.
Lecture 1 1 1 Department of Mathematics Dartmouth College 3/28/07 Outline Course Overview http://www.math.dartmouth.edu/~m23s07 Matlab Ordinary differential equations Definition An ordinary differential
More informationMath 310 Introduction to Ordinary Differential Equations Final Examination August 9, Instructor: John Stockie
Make sure this exam has 15 pages. Math 310 Introduction to Ordinary Differential Equations inal Examination August 9, 2006 Instructor: John Stockie Name: (Please Print) Student Number: Special Instructions
More information8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).
8. THEOREM I the partial derivatives and eist near (a b) and are continuous at (a b) then is dierentiable at (a b). For a dierentiable unction o two variables z= ( ) we deine the dierentials d and d to
More informationConstrained Optimal Control I
Optimal Control, Guidance and Estimation Lecture 34 Constrained Optimal Control I Pro. Radhakant Padhi Dept. o Aerospace Engineering Indian Institute o Science - Bangalore opics Motivation Brie Summary
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI097 TU) Thursday January 014, 18:0-1:0
More informationIntersection of Ellipsoids
Intersection of Ellipsoids David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationMathematical Notation Math Calculus & Analytic Geometry III
Name : Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and can e printed and given to the instructor
More informationThe Direct Transcription Method For Optimal Control. Part 2: Optimal Control
The Direct Transcription Method For Optimal Control Part 2: Optimal Control John T Betts Partner, Applied Mathematical Analysis, LLC 1 Fundamental Principle of Transcription Methods Transcription Method
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
More informationAPPM 2360: Midterm exam 1 February 15, 2017
APPM 36: Midterm exam 1 February 15, 17 On the front of your Bluebook write: (1) your name, () your instructor s name, (3) your recitation section number and () a grading table. Text books, class notes,
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationProblem set 6 for Quantum Field Theory course
Problem set 6 or Quantum Field Theory course 2018.03.13. Toics covered Scattering cross-section and decay rate Yukawa theory and Yukawa otential Scattering in external electromagnetic ield, Rutherord ormula
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
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 informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationLecture 9. Systems of Two First Order Linear ODEs
Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 1 Dianne P. O Leary c 2008 1 The
More informationDynamical Systems & Scientic Computing: Homework Assignments
Fakultäten für Informatik & Mathematik Technische Universität München Dr. Tobias Neckel Summer Term Dr. Florian Rupp Exercise Sheet 3 Dynamical Systems & Scientic Computing: Homework Assignments 3. [ ]
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationFirst In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 20 September 2018
First In-Class Exam Solutions Math 246, Professor David Levermore Thursday, 20 September 208 () [6] In the absence of predators the population of mosquitoes in a certain area would increase at a rate proportional
More informationImage Enhancement (Spatial Filtering 2)
Image Enhancement (Spatial Filtering ) Dr. Samir H. Abdul-Jauwad Electrical Engineering Department College o Engineering Sciences King Fahd University o Petroleum & Minerals Dhahran Saudi Arabia samara@kupm.edu.sa
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationFeedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control
More informationWhy are Discrete Maps Sufficient?
Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well
More informationdy dt = ty, y(0) = 3. (1)
2. (10pts) Solve the given intial value problem (IVP): dy dt = ty, y(0) = 3. (1) 3. (10pts) A plot of f(y) =y(1 y)(2 y) of the right hand side of the differential equation dy/dt = f(y) is shown below.
More informationLinear algebra and differential equations (Math 54): Lecture 20
Linear algebra and differential equations (Math 54): Lecture 20 Vivek Shende April 7, 2016 Hello and welcome to class! Last time We started discussing differential equations. We found a complete set of
More informationME8230 Nonlinear Dynamics
ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass
More informationMOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1
MOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1 OPTIONAL question 1. Below, several phase portraits are shown.
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationSolutions to Math 53 Math 53 Practice Final
Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points
More informationChristoffel symbols and Gauss Theorema Egregium
Durham University Pavel Tumarkin Epiphany 207 Dierential Geometry III, Solutions 5 (Week 5 Christoel symbols and Gauss Theorema Egregium 5.. Show that the Gauss curvature K o the surace o revolution locally
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
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 informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationLecture 31. Basic Theory of First Order Linear Systems
Math 245 - Mathematics of Physics and Engineering I Lecture 31. Basic Theory of First Order Linear Systems April 4, 2012 Konstantin Zuev (USC) Math 245, Lecture 31 April 4, 2012 1 / 10 Agenda Existence
More informationMecE 390 Final examination, Winter 2014
MecE 390 Final examination, Winter 2014 Directions: (i) a double-sided 8.5 11 formula sheet is permitted, (ii) no calculators are permitted, (iii) the exam is 80 minutes in duration; please turn your paper
More informationBasic mathematics of economic models. 3. Maximization
John Riley 1 January 16 Basic mathematics o economic models 3 Maimization 31 Single variable maimization 1 3 Multi variable maimization 6 33 Concave unctions 9 34 Maimization with non-negativity constraints
More information1 The relation between a second order linear ode and a system of two rst order linear odes
Math 1280 Spring, 2010 1 The relation between a second order linear ode and a system of two rst order linear odes In Chapter 3 of the text you learn to solve some second order linear ode's, such as x 00
More informationMesa College Math SAMPLES
Mesa College Math 6 - SAMPLES Directions: NO CALCULATOR. Write neatly, show your work and steps. Label your work so it s easy to ollow. Answers without appropriate work will receive NO credit. For inal
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
More informationOutline. Approximate sampling theorem (AST) recall Lecture 1. P. L. Butzer, G. Schmeisser, R. L. Stens
Outline Basic relations valid or the Bernstein space B and their extensions to unctions rom larger spaces in terms o their distances rom B Part 3: Distance unctional approach o Part applied to undamental
More information