Podlubny I., Skovranek T, Verbickij V., Vinagre B., Chen YQ., Petras I.
|
|
- Rudolf Barnett
- 5 years ago
- Views:
Transcription
1 DISCRETE FRACTIONAL CALCULUS: NON-EQUIDISTANT GRIDS AND VARIABLE STEP LENGTH Podlubny I., Skovranek T, Verbickij V., Vinagre B., Chen YQ., Petras I. Technical University of Kosice, Slovakia Odessa National University, Odessa, Ukraine University of Extremadura, Badajoz, Spain Utah State University, Logan, USA
2 History: OFDEs: 2000, PFDEs: 2008, SOFDEs: 2009 MATRIX APPROACH TO DISCRETE FRACTIONAL CALCULUS Igor Podlubny 2
3 Implemented in MATLAB (2008, updated 2009)
4 Triangular strip matrices (TSMs)
5 Generating functions for TSMs
6 Operations with TSMs
7 Left-sided fractional derivatives
8 Discrete analog of left-sided fractional derivatives
9 Solving ordinary FDEs Bagley-Torvik equation: ay (t)+b 0 D 3/2 t y(t)+cy(t) =f(t) add zero initial conditions: ab (2) + bb n n (3/2) + cb n (0) Yn = F n, 9
10 Solving partial FDEs β U x β β U x β 10
11 Left-sided fractional integrals: inverse of fractional derivatives
12 Non-equidistant grid: discretization is cumbersome even in simplest cases... w (x) : 1 ĥ i ( wi+1 w i h i w i w i 1 h i 1 ) h i = x i+1 x i, ĥ i =ˆx i ˆx i 1, e difference proble
13 Recall that I α N = BN α 1 Change the viewpoint: Left-sided fractional derivatives: inverse of fractional integrals; then B α N =(I α N ) 1 Any approximation of fractional integration after inversion gives an approximation for fractional differentiation on the same grid!
14 The simplest approach: approximation of a function under integration by a piecewise constant function B α N =(I α N ) 1 Coefficients of I α N as I k,j = (t k t j 1 ) α (t k t j ) α Γ(α + 1) j =1,...,k; k =1,...,N., For non-equidistant grids, the matrix is not a TSM.
15 Example 1: fractional integrals of sin(x) y = 0.1 = 0.3 = 0.5 = 0.7 = [ ] with non equidistant step t
16 Example 2: fractional derivatives of sin(x) y = 0.1 = 0.3 = 0.5 = 0.7 = [0.1, 0.3, 0.5, 0.7], with non equidistant step t
17 Example 3: two-term ordinary FDE y (α) (t)+y(t) =1, Exact solution: y(0) = 0, y (0) = 0. y(t) =t α E α,α+1 ( t α ) y analytical solution 0.2 numerical solution (non equidistant step) t α = 1.8, number of discretization nodes N = 500
18 Example 4(a): Bagley-Torvik equation Ay (t)+by 3/2 (t)+cy(t) =F (t), F (t) = 8, (0 t 1) 0, (t >1), y(0) = y (0) = equidistant step non equidistant step 3 y t A = 1, B = 1, C = 1.
19 Example 4(b): Bagley-Torvik equation Ay (t)+by 3/2 (t)+cy(t) =F (t), F (t) = 8, (0 t 1) 0, (t >1), y(0) = y (0) = equidistant step non equidistant step 4 y t A = 1, B = 0.5, C = 0.5.
20 Can we have variable step length? As seen in MATLAB: ode23.m and ode45.m solvers optimized for a variable step using a variable step ensures that a large step size is used for low frequencies and a small step size is used for high frequencies make a step, estimate the error at this step, check if the value is greater than or less than the tolerance, and adapt the step size accordingly There was nothing like this available for FDEs so far...
21 Method of large steps 0Dt α y(t) =f(y(t),t), (t >0), y(0) = 0, Suppose we obtained its solution in the interval (0,a) (and the final value y a at t =a), then we can use this for transforming the above problem to adt α y(t) =f(y(t),t) 0 Ra α y(t), (t >a), y(a) =y a, where 0R α a y(t) = 1 Γ(1 α) 0R α a y(t) = 0 D α t a 0 (t τ) α 1 y(τ)dτ, (t >a). (1 H(t a))y(t)
22 First large step in [0, a] 0D α t y(t) =f(y(t),t), (t >0), y(0) = 0, ad α t y(t) =f(y(t),t) 0 R α a y(t), (t >a), Auxiliary function: y(a) =y a, Second large step in [a, b] y(t) = u(t) +y a, the function u(t), w ad α t u(t) =f(u(t)+y a,t) 0 R α a y(t) y a, u(a) =0. (t >a),
23 Method of large steps : example (1) 0D 1/2 t y(t)+y(t) = t1.5 Γ(1.5) y(0) = 0. + t, (t >0), Exact solution: y(t) =t. First large step : interval [0, 1]: clear all h = 0.01; Using the matrix approach t = 0:h:1; N = 1/h + 1; M = zeros(n,n); M = ban(0.5, N, h) + eye(n,n); F = (t.^(0.5)/gamma(1.5) + t) ; M = eliminator(n,[1])*m*eliminator(n, [1]) ; F = eliminator(n,[1])*f; Y = M\F; Y0 = [0; Y]; plot (t,y0, b ) set(gca, xlim, [0 2], ylim, [0 2] ) grid on, hold on
24 Method of large steps : example (2) Second large step : interval [1, 2] 0D 1/2 t y(t) = 1 D 1/2 t y(t)+ 1 Γ(0.5) 1 0 y (τ)dτ, (t >1) (t τ) 1/2 1D 1/2 t y(t)+y(t) = t1.5 Γ(1.5) + t 1 Γ(0.5) 1 0 dτ (t τ) 1/2 (t>1). 1D α t y(t)+y(t) = t1.5 Γ(1.5) + t 2t0.5 Γ(0.5) + 2(t 1)0.5 Γ(0.5) ; (t>1) y(1) = 1. y(t) =u(t)+1,
25 Method of large steps : example (3) The problem to solve in [1, 2]: 1D α t u(t)+u(t) = t1.5 Γ(1.5) + t 2t0.5 Γ(0.5) + 2(t 1)0.5 Γ(0.5) 1; (t >1) u(1) = 0. Using the matrix approach clear all h = 0.01; t = 1:h:2; N = 1/h + 1; M = zeros(n,n); M = ban(0.5, N, h) + eye(n,n); F = (t.^(0.5)/gamma(1.5) + t - 2*t.^(0.5)/gamma(0.5) *(t-1).^(0.5)/gamma(0.5) - 1) ; M = eliminator(n,[1])*m*eliminator(n, [1]) ; F = eliminator(n,[1])*f; U = M \F; U0 = [0; U]; Y0 = U0 + 1; plot(t, Y0, g )
26 Conclusions Easy, transparent, uniform approach to fractional integration and differentiation on non-equidistant grids No need for cumbersome formulas for finite differences on nonequidistant grids Any numerical approximation of a fractional-order integral on a non-equidistant mesh gives a matrix for fractional-order differentiation on the same grid Works for left-sided, right-sided, and symmetric (Riesz) fractional integrals and derivatives Can be used for solving PDEs on non-uniform spatial-temporal meshes Moving towards variable-length and adaptive-length step methods in numerical solution of ordinary and partial FDEs.
27 Thank you!
MATRIX APPROACH TO DISCRETE FRACTIONAL CALCULUS. Igor Podlubny. Abstract
MATRIX APPROACH TO DISCRETE FRACTIOAL CALCULUS Igor Podlubny Dedicated to Prof Rudolf Gorenflo, on the occasion of his 70-th birthday Abstract A matrix form representation of discrete analogues of various
More informationA ROBUST STABILITY TEST PROCEDURE FOR A CLASS OF UNCERTAIN LTI FRACTIONAL ORDER SYSTEMS
International Carpathian Control Conference ICCC 22 MALENOVICE, CZECH REPUBLIC May 27-, 22 A ROBUST STABILITY TEST PROCEDURE FOR A CLASS OF UNCERTAIN LTI FRACTIONAL ORDER SYSTEMS Ivo PETRÁŠ, YangQuan CHEN
More informationA computationally effective predictor-corrector method for simulating fractional order dynamical control system
ANZIAM J. 47 (EMA25) pp.168 184, 26 168 A computationally effective predictor-corrector method for simulating fractional order dynamical control system. Yang F. Liu (Received 14 October 25; revised 24
More informationNUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 6(2) (217), 515 523 Palestine Polytechnic University-PPU 217 NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX Raghvendra
More informationx(t) = t[u(t 1) u(t 2)] + 1[u(t 2) u(t 3)]
ECE30 Summer II, 2006 Exam, Blue Version July 2, 2006 Name: Solution Score: 00/00 You must show all of your work for full credit. Calculators may NOT be used.. (5 points) x(t) = tu(t ) + ( t)u(t 2) u(t
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
More informationUsing MATLAB with the Convolution Method
ECE 350 Linear Systems I MATLAB Tutorial #5 Using MATLAB with the Convolution Method A linear system with input, x(t), and output, y(t), can be described in terms of its impulse response, h(t). x(t) h(t)
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Pseudo-Time Integration 1 / 10 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 10 Outline 1
More informationComputational methods for large-scale linear matrix equations and application to FDEs. V. Simoncini
Computational methods for large-scale linear matrix equations and application to FDEs V. Simoncini Dipartimento di Matematica, Università di Bologna valeria.simoncini@unibo.it Joint work with: Tobias Breiten,
More informationChapter 3 Convolution Representation
Chapter 3 Convolution Representation DT Unit-Impulse Response Consider the DT SISO system: xn [ ] System yn [ ] xn [ ] = δ[ n] If the input signal is and the system has no energy at n = 0, the output yn
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 informationLecture 5 Least-squares
EE263 Autumn 2008-09 Stephen Boyd Lecture 5 Least-squares least-squares (approximate) solution of overdetermined equations projection and orthogonality principle least-squares estimation BLUE property
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts Signals A signal is a pattern of variation of a physical quantity as a function of time, space, distance, position, temperature, pressure, etc. These quantities are usually
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
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 informationTheory, Solution Techniques and Applications of Singular Boundary Value Problems
Theory, Solution Techniques and Applications of Singular Boundary Value Problems W. Auzinger O. Koch E. Weinmüller Vienna University of Technology, Austria Problem Class z (t) = M(t) z(t) + f(t, z(t)),
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationThe Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers
The Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers Chengbin Ma and Yoichi Hori Electrical Control System Engineering, Information & System Division, Institute of Industrial
More informationNumerical Differential Equations: IVP
Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential
More informationHomework Solutions:
Homework Solutions: 1.1-1.3 Section 1.1: 1. Problems 1, 3, 5 In these problems, we want to compare and contrast the direction fields for the given (autonomous) differential equations of the form y = ay
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (211) 219 223 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Laplace transform and fractional differential
More informationSWARM INTELLIGENCE OPTIMIZED NEURAL NETWORKS FOR SOLVING FRACTIONAL DIFFERENTIAL EQUATIONS. Received March 2010; revised August 2010
International Journal of Innovative Computing, Information and Control ICIC International c 2011 ISSN 1349-4198 Volume 7, Number 11, November 2011 pp. 6301 6318 SWARM INTELLIGENCE OPTIMIZED NEURAL NETWORKS
More informationTHE UNIVERSITY OF WESTERN ONTARIO. Applied Mathematics 375a Instructor: Matt Davison. Final Examination December 14, :00 12:00 a.m.
THE UNIVERSITY OF WESTERN ONTARIO London Ontario Applied Mathematics 375a Instructor: Matt Davison Final Examination December 4, 22 9: 2: a.m. 3 HOURS Name: Stu. #: Notes: ) There are 8 question worth
More informationNumerical solutions of ordinary fractional differential equations with singularities
arxiv:86.3269v [mat.na] 8 Jun 28 Numerical solutions of ordinary fractional differential equations wit singularities Yuri Dimitrov, Ivan Dimov, Venelin Todorov Abstract Te solutions of fractional differential
More informationDETC DISCRETE FRACTIONAL CALCULUS: NON-EQUIDISTANT GRIDS AND VARIABLE STEP LENGTH
Proceedings of he ASME Inernaional Design Engineering Technical onferences & ompuers and Informaion in Engineering onference IDET/IE Augus 9-3,, Washingon, D, USA DET-4763 DISRETE FRATIONAL ALULUS: NON-EQUIDISTANT
More informationHomework 3 Solutions
18-290 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 2018 Homework 3 Solutions Part One 1. (25 points) The following systems have x(t) or x[n] as input and y(t) or y[n] as output. For each
More informationA BRIEF INTRODUCTION INTO SOLVING DIFFERENTIAL EQUATIONS
MATTHIAS GERDTS A BRIEF INTRODUCTION INTO SOLVING DIFFERENTIAL EQUATIONS Universität der Bundeswehr München Addresse des Autors: Matthias Gerdts Institut für Mathematik und Rechneranwendung Universität
More informationSOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 3(2) (214), 281 288 Palestine Polytechnic University-PPU 214 SOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX Naresh Berwal, Dinesh Panchal and
More information1.17 : Consider a continuous-time system with input x(t) and output y(t) related by y(t) = x( sin(t)).
(Note: here are the solution, only showing you the approach to solve the problems. If you find some typos or calculation error, please post it on Piazza and let us know ).7 : Consider a continuous-time
More informationDifferential Equations Class Notes
Differential Equations Class Notes Dan Wysocki Spring 213 Contents 1 Introduction 2 2 Classification of Differential Equations 6 2.1 Linear vs. Non-Linear.................................. 7 2.2 Seperable
More informationEN Applied Optimal Control Lecture 8: Dynamic Programming October 10, 2018
EN530.603 Applied Optimal Control Lecture 8: Dynamic Programming October 0, 08 Lecturer: Marin Kobilarov Dynamic Programming (DP) is conerned with the computation of an optimal policy, i.e. an optimal
More informationUniversity of Kentucky Department of Electrical and Computer Engineering. EE421G: Signals and Systems I Fall 2007
University of Kentucky Department of Electrical an Computer Engineering EE4G: Signals an Systems I Fall 7 Issue: October 4, 7 Problem Set 6 Due: October, 7 (In class) Reaing Assignments: Rea Chapter 3.7,
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 4, 6, 2M Open access books available International authors and editors Downloads Our authors are
More informationFigure 1 A linear, time-invariant circuit. It s important to us that the circuit is both linear and time-invariant. To see why, let s us the notation
Convolution In this section we consider the problem of determining the response of a linear, time-invariant circuit to an arbitrary input, x(t). This situation is illustrated in Figure 1 where x(t) is
More informationMULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS
MULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS Hossein Jafari & M. A. Firoozjaee Young Researchers club, Islamic Azad University, Jouybar Branch, Jouybar, Iran
More informationLinear Multistep Methods I: Adams and BDF Methods
Linear Multistep Methods I: Adams and BDF Methods Varun Shankar January 1, 016 1 Introduction In our review of 5610 material, we have discussed polynomial interpolation and its application to generating
More informationDifferential Equations
Differential Equations Math 341 Fall 21 MWF 2:3-3:25pm Fowler 37 c 21 Ron Buckmire http://faculty.oxy.edu/ron/math/341/1/ Worksheet 29: Wednesday December 1 TITLE Laplace Transforms and Introduction to
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More information/ / MET Day 000 NC1^ INRTL MNVR I E E PRE SLEEP K PRE SLEEP R E
05//0 5:26:04 09/6/0 (259) 6 7 8 9 20 2 22 2 09/7 0 02 0 000/00 0 02 0 04 05 06 07 08 09 0 2 ay 000 ^ 0 X Y / / / / ( %/ ) 2 /0 2 ( ) ^ 4 / Y/ 2 4 5 6 7 8 9 2 X ^ X % 2 // 09/7/0 (260) ay 000 02 05//0
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationMA Ordinary Differential Equations and Control Semester 1 Lecture Notes on ODEs and Laplace Transform (Parts I and II)
MA222 - Ordinary Differential Equations and Control Semester Lecture Notes on ODEs and Laplace Transform (Parts I and II Dr J D Evans October 22 Special thanks (in alphabetical order are due to Kenneth
More informationLecture 29. Convolution Integrals and Their Applications
Math 245 - Mathematics of Physics and Engineering I Lecture 29. Convolution Integrals and Their Applications March 3, 212 Konstantin Zuev (USC) Math 245, Lecture 29 March 3, 212 1 / 13 Agenda Convolution
More informationNumerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)
Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation
More informationMathematical Methods and its Applications (Solution of assignment-12) Solution 1 From the definition of Fourier transforms, we have.
For 2 weeks course only Mathematical Methods and its Applications (Solution of assignment-2 Solution From the definition of Fourier transforms, we have F e at2 e at2 e it dt e at2 +(it/a dt ( setting (
More informationAn Analytical Scheme for Multi-order Fractional Differential Equations
Tamsui Oxford Journal of Mathematical Sciences 26(3) (2010) 305-320 Aletheia University An Analytical Scheme for Multi-order Fractional Differential Equations H. M. Jaradat Al Al Bayt University, Jordan
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 4, pp. 373-38, 23. Copyright 23,. ISSN 68-963. ETNA A NOTE ON THE RELATION BETWEEN THE NEWTON HOMOTOPY METHOD AND THE DAMPED NEWTON METHOD XUPING ZHANG
More informatione st f (t) dt = e st tf(t) dt = L {t f(t)} s
Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationDiscussion Section #2, 31 Jan 2014
Discussion Section #2, 31 Jan 2014 Lillian Ratliff 1 Unit Impulse The unit impulse (Dirac delta) has the following properties: { 0, t 0 δ(t) =, t = 0 ε ε δ(t) = 1 Remark 1. Important!: An ordinary function
More informationAsymptotic Behaviour of Solutions of n-order Differential Equations with Regularly Varying Nonlinearities
International Workshop QUALITDE 206, December 24 26, 206, Tbilisi, Georgi4 Asymptotic Behaviour of Solutions of n-order Differential Equations with Regularly Varying Nonlinearities K. S. Korepanova Odessa
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 informationKevin James. MTHSC 3110 Section 2.1 Matrix Operations
MTHSC 3110 Section 2.1 Matrix Operations Notation Let A be an m n matrix, that is, m rows and n columns. We ll refer to the entries of A by their row and column indices. The entry in the i th row and j
More informationMAT292 - Calculus III - Fall Solution for Term Test 2 - November 6, 2014 DO NOT WRITE ON THE QR CODE AT THE TOP OF THE PAGES.
MAT9 - Calculus III - Fall 4 Solution for Term Test - November 6, 4 Time allotted: 9 minutes. Aids permitted: None. Full Name: Last First Student ID: Email: @mail.utoronto.ca Instructions DO NOT WRITE
More informationSOME APPROXIMATIONS OF FRACTIONAL ORDER OPERATORS USED IN CONTROL THEORY AND APPLICATIONS. Abstract
SOME APPROXIMATIONS OF FRACTIONAL ORDER OPERATORS USED IN CONTROL THEORY AND APPLICATIONS B. M. Vinagre, I. Podlubny, A. Hernández, V. Feliu Abstract In this paper, starting from the formulation of some
More informationLTI Systems (Continuous & Discrete) - Basics
LTI Systems (Continuous & Discrete) - Basics 1. A system with an input x(t) and output y(t) is described by the relation: y(t) = t. x(t). This system is (a) linear and time-invariant (b) linear and time-varying
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationMath 341 Fall 2008 Friday December 12
FINAL EXAM: Differential Equations Math 341 Fall 2008 Friday December 12 c 2008 Ron Buckmire 1:00pm-4:00pm Name: Directions: Read all problems first before answering any of them. There are 17 pages in
More informationDISCRETE-TIME SYSTEMS AND DIGITAL CONTROLLERS
DISCRETE-TIME SYSTEMS AND DIGITAL CONTROLLERS In practice controllers are nowadays almost exclusively implemented digitally. This means that the controller operates in discrete time, although the controlled
More informationAnalytic solution of fractional integro-differential equations
Annals of the University of Craiova, Mathematics and Computer Science Series Volume 38(1), 211, Pages 1 1 ISSN: 1223-6934 Analytic solution of fractional integro-differential equations Fadi Awawdeh, E.A.
More informationNumerical Methods - Boundary Value Problems for ODEs
Numerical Methods - Boundary Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Boundary Value Problems for ODEs 2013 1 / 14 Outline 1 Boundary Value
More informationSection 3: Summary. Section 4. Evaluating Convolution Integrals. Convolution integral (memorise this):
Section 3: Summary Convolution integral (memorise this): f(t) input g(t) impulse response y(t) output y(t) g(t τ) f(τ) dτ Way to find the output of a linear system, described by a differential equation,
More informationNumerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)
Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationSolving linear equations with Gaussian Elimination (I)
Term Projects Solving linear equations with Gaussian Elimination The QR Algorithm for Symmetric Eigenvalue Problem The QR Algorithm for The SVD Quasi-Newton Methods Solving linear equations with Gaussian
More informationLecture 7 - Separable Equations
Lecture 7 - Separable Equations Separable equations is a very special type of differential equations where you can separate the terms involving only y on one side of the equation and terms involving only
More informationarxiv:math/ v1 [math.oc] 11 Apr 2000
The fractional - order controllers: Methods for their synthesis and application arxiv:math/0004064v1 [math.oc] 11 Apr 2000 Abstract This paper deals with fractional-order controllers. We outline mathematical
More informationWinfried Auzinger, Ernst Karner, Othmar Koch, Ewa Weinmüller
Opuscula Mathematica Vol. 26 No. 2 26 Winfried Auzinger, Ernst Karner, Othmar Koch, Ewa Weinmüller COLLOCATION METHODS FOR THE SOLUTION OF EIGENVALUE PROBLEMS FOR SINGULAR ORDINARY DIFFERENTIAL EQUATIONS
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
More informationEE 210. Signals and Systems Solutions of homework 2
EE 2. Signals and Systems Solutions of homework 2 Spring 2 Exercise Due Date Week of 22 nd Feb. Problems Q Compute and sketch the output y[n] of each discrete-time LTI system below with impulse response
More information1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is
ECE 55, Fall 2007 Problem Set #4 Solution The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ Ax + Bu is x(t) e A(t ) x( ) + e A(t τ) Bu(τ)dτ () This formula is extremely important
More informationSection 2.4 Linear Equations
Section 2.4 Linear Equations Key Terms: Linear equation Homogeneous linear equation Nonhomogeneous (inhomogeneous) linear equation Integrating factor General solution Variation of parameters A first-order
More informationEssential Ordinary Differential Equations
MODULE 1: MATHEMATICAL PRELIMINARIES 10 Lecture 2 Essential Ordinary Differential Equations In this lecture, we recall some methods of solving first-order IVP in ODE (separable and linear) and homogeneous
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 informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationECE 301 Division 1 Exam 1 Solutions, 10/6/2011, 8-9:45pm in ME 1061.
ECE 301 Division 1 Exam 1 Solutions, 10/6/011, 8-9:45pm in ME 1061. Your ID will be checked during the exam. Please bring a No. pencil to fill out the answer sheet. This is a closed-book exam. No calculators
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More informationComputing inverse Laplace Transforms.
Review Exam 3. Sections 4.-4.5 in Lecture Notes. 60 minutes. 7 problems. 70 grade attempts. (0 attempts per problem. No partial grading. (Exceptions allowed, ask you TA. Integration table included. Complete
More informationAnalysis of Fractional-order Telegraph Model for Neutron Transport in Nuclear Reactor with Slab Geometry
23 European Control Conference (ECC July 7-9, 23, Zürich, Switzerland. Analysis of Fractional-order Telegraph Model for Neutron Transport in Nuclear Reactor with Slab Geometry Vishwesh A. Vyawahare,2 and
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 informatione jωt y (t) = ω 2 Ke jωt K =
BME 171, Sec 2: Homework 2 Solutions due Tue, Sep 16 by 5pm 1. Consider a system governed by the second-order differential equation a d2 y(t) + b dy(t) where a, b and c are nonnegative real numbers. (a)
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 informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationNew computational method for solving fractional Riccati equation
Available online at www.isr-publications.com/jmcs J. Math. Computer Sci., 17 2017), 106 114 Research Article Journal Homepage: www.tjmcs.com - www.isr-publications.com/jmcs New computational method for
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationComplete Identification of a Dynamic Fractional Order System Under Non-ideal Conditions Using Fractional Differintegral Definitions
Complete Identification of a Dynamic Fractional Order System Under Non-ideal Conditions Using Fractional Differintegral Definitions Deepyaman Maiti # Ayan Acharya #2 R. Janarthanan *3 Amit Konar #4 # Department
More informationUsing Fractional Calculus for Lateral and Longitudinal Control of Autonomous Vehicles
Using Fractional Calculus for Lateral and Longitudinal Control of Autonomous Vehicles J. I. Suárez 1,B.M.Vinagre 1,A.J.Calderón 1,C.A.Monje 1,andY.Q.Chen 2 1 EII, Universidad de Extremadura, Badajoz -
More informationAnalytical solutions for the fractional Klein-Gordon equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 2, No. 2, 214, pp. 99-114 Analytical solutions for the fractional Klein-Gordon equation Hosseni Kheiri Faculty of Mathematical
More informationHomework 4. May An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt
Homework 4 May 2017 1. An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt Determine the impulse response of the system. Rewriting as y(t) = t e (t
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 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationto have roots with negative real parts, the necessary and sufficient conditions are that:
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 543 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 7 Sebastian A. Nugroho November 6, 7 Due date of the homework is: Sunday, November 6th @ :59pm.. The following
More informationLecture Notes on Numerical Differential Equations: IVP
Lecture Notes on Numerical Differential Equations: IVP Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu URL:
More informationFinal Exam Solutions
EE55: Linear Systems Final Exam SIST, ShanghaiTech Final Exam Solutions Course: Linear Systems Teacher: Prof. Boris Houska Duration: 85min YOUR NAME: (type in English letters) I Introduction This exam
More informationNumerical Solution of Fuzzy Differential Equations
Applied Mathematical Sciences, Vol. 1, 2007, no. 45, 2231-2246 Numerical Solution of Fuzzy Differential Equations Javad Shokri Department of Mathematics Urmia University P.O. Box 165, Urmia, Iran j.shokri@mail.urmia.ac.ir
More informationRemark on the Sensitivity of Simulated Solutions of the Nonlinear Dynamical System to the Used Numerical Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 55, 2749-2754 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.59236 Remark on the Sensitivity of Simulated Solutions of
More informationNUMERICAL SOLUTION OF FRACTIONAL DIFFUSION-WAVE EQUATION WITH TWO SPACE VARIABLES BY MATRIX METHOD. Mridula Garg, Pratibha Manohar.
NUMERICAL SOLUTION OF FRACTIONAL DIFFUSION-WAVE EQUATION WITH TWO SPACE VARIABLES BY MATRIX METHOD Mridula Garg, Pratibha Manohar Abstract In the present paper we solve space-time fractional diffusion-wave
More information