AN OPTIMAL CONTROL PROBLEM FOR SPLINES ASSOCIATED TO LINEAR DIFFERENTIAL OPERATORS
|
|
- Gloria McCormick
- 6 years ago
- Views:
Transcription
1 CONTROLO 6 7h Poruguese Conference on Auomaic Conrol Insiuo Superior Técnico, Lisboa, Porugal Sepember -3, 6 AN OPTIMAL CONTROL PROBLEM FOR SPLINES ASSOCIATED TO LINEAR DIFFERENTIAL OPERATORS Rui C. Rodrigues, F. Silva Leie, Insiuo Superior de Engenharia de Coimbra Deparameno de Física e Maemáica Coimbra, Porugal and Insiuo de Sisemas e Robóica Universidade de Coimbra - Pólo II 33-9 Coimbra, Porugal ruicr@isec.p Deparameno de Maemáica Universidade de Coimbra Coimbra, Porugal and Insiuo de Sisemas e Robóica Universidade de Coimbra - Pólo II 33-9 Coimbra, Porugal fleie@ma.uc.p Absrac: We formulae an opimal conrol problem, where he dynamics are given by an inpu-oupu equaion, in order o emphasize he opimal properies of some inerpolaing curves, known in he lieraure as L-splines. We presen he soluion of he opimal conrol problem and illusrae he heory wih some examples. Keywords: Opimal conrol, Linear sysems, Splines.. INTRODUCTION In many engineering applicaions one has o find rajecories ha are consrained o pass hrough some specific poins wih prescribed derivaives and some required degree of smoohness. Examples of hese are he cubic splines ha also minimize acceleraion and oher polynomial funcions, such as quinic polynomials ha minimize jerk. Parially suppored by PRODEP/3 (Medida 5, Acção 5.3), Projec NCT4 (FCT) and ISR-Coimbra. Parially suppored by Projec NCT4 (FCT) and ISR- Coimbra. One common feaure o hese examples is ha he number of inerpolaion condiions, as well as he number of smoohness condiions, is he same a each poin (excep possibly he firs and he las). Less known are he curves named L-splines. These are obained by joining smoohly pieces of soluions of cerain linear differenial equaions. Conrary o he previous examples, hey allow he prescripion of uneven condiions a each poin. In his aricle we demonsrae he opimal behavior of L-splines, by showing ha hese curves
2 can be realized as he opimal oupu of a cerain opimal conrol problem wih piecewise conrols. The formulaion of his problem appears in secion. The dynamics are given by an appropriae inpu-oupu equaion, in opposiion o he usual classical formulaion, where he dynamics appear in sae-space form. The connecion wih L-splines is presened in secion 3, and he properies of hese curves are used, in secion 4, o presen he soluion of he proposed opimal conrol problem. An equivalen sae-space formulaion appears in secion 5. Finally, in he las secion, we include hree examples o illusrae he heory. This paper highlighs several aspecs of he opimal behavior of splines associaed o linear differenial operaors reaed in (Rodrigues and Silva Leie, ) and correcs some inaccuracies and imprecise saemens conained in his aricle.. FORMULATION OF THE OPTIMAL CONTROL PROBLEM Le : a = < < < m = b be any pariion of he ime inerval [a, b] R, p N, and Z = (z, z, z,...,z m, z m ), wih z = z m = p and z i p, any incidence vecor for. L and L c are wo linear differenial operaors wih real consan coefficiens given by L D p + a p D p + + a D + a, L c c p D p + + c D + c. We do no impose c p and define σ as σ = maxj : c j }. We assume ha he characerisic polynomials associaed o L and L c, have no common facors. The admissible conrols are piecewise coninuous funcions, bounded over [a, b], which are also C p+σ -smooh in each subinerval [ i, i+ ]. As usual, he se of all admissible conrols is denoed by U. We consider he following opimal conrol problem, denoed by (P c ). b a L y() = L c u(), (L c u()) d u U he inerpolaion condiions y (k) ( i ) = α k i, i =,...,m, k =,...,z i, he boundary condiions y (k) ( ) = η k, y(k) ( m ) = η k m, k =,...,p, min, and he smoohness condiions y (k) ( i ) = y (k) ( i + ), i =,...,m, k =,...,p z i, where α k i, ηk and η k m are prescribed real numbers. This is no a classical formulaion of an opimal conrol problem since he dynamics are described by an inpu-oupu equaion and he opimal oupu is required o saisfy a cerain number of uneven condiions a he kno poins i. Noe ha here is a direc connecion beween he incidence vecor Z and he number of condiions prescribed a each kno poin of he pariion. Indeed, here are z = p iniial condiions a he ime insan, z m = p final condiions a he ime insan m, and z i inerpolaion condiions a i, i =,...,m. Taking ino accoun he number of prescribed smoohness condiions, one ges a oal of p condiions a each i, i =,...,m, and p condiions a and m. For deails concerning linear conrol sysems see, for insance, (Kailah, 98), (Ribeiro, a), (Ribeiro, b) and (Rugh, 996). 3. CONNECTIONS WITH L-SPLINES In his secion we show ha he opimal oupu of he problem (P c ) is a curve named L-spline, known in he lieraure for many years (Schulz and Varga, 967). Esablishing he connecion beween our opimal conrol problem and hese curves, will highligh he opimal behavior of L- splines already demonsraed in (Rodrigues and Silva Leie, ). This complemens wha has already been done for oher kinds of splines, namely for he so called generalized splines (see (Ahlberg e al., 964)), in (Zhang e al., 997) and (Rodrigues e al., 999). The opimal conrol problem (P c ) has a naural variaional formulaion as: J(y( )) = b a (Ly()) d y Ω subjec o he inerpolaion condiions min, y (k) ( i ) = α k i, i =,...,m, k =,...,z i, he boundary condiions y (k) ( ) = η k, y(k) ( m ) = η k m, k =,...,p, and he smoohness condiions y (k) ( i ) = y (k) ( i + ), i =,...,m, k =,...,p z i, where Ω is he se of funcions which are C p - smooh in each subinerval [ i, i+ ].
3 We denoe by Ω A he se of all admissible curves for his variaional problem, and by L he adjoin of he operaor L, ha is, L ( ) p D p + ( ) p a p D p + +a. Nex, we sae a necessary and sufficien condiion for an admissible curve o minimize he given funcional, which can be derived from resuls in (Schulz and Varga, 967). Theorem y Ω A is a soluion of he variaional problem above if and only if y is a soluion of he linear differenial equaion L L x = in every subinerval [ i, i+ ]. Remark () The differenial equaion in he previous heorem is nohing bu he Euler-Lagrange equaion associaed o he variaional problem. () I urns ou ha here is only one curve which is a soluion of his differenial equaion (of order p) on each subinerval and saisfies he prescribed inerpolaion, boundary and smoohness condiions. Thus, he variaional problem has one and only one soluion. This soluion is an L-spline (of ype I) ha can be easily compued by solving a consisen se of pm linear algebraic equaions in pm unknowns. 4. SOLVING THE OPTIMAL CONTROL PROBLEM We now discuss how o solve he opimal conrol problem (P c ) using he facs presened in he previous secion. Theorem When σ, problem (P c ) has several soluions, which are piecewise coninuous funcions composed by m segmens, each one depending on σ free parameers. When σ =, problem (P c ) has a unique soluion. Proof. Since he variaional problem has a soluion, i is clear ha problem (P c ) has, a leas, a soluion, which generaes an L-spline curve as opimal oupu. Le ỹ denoe his L-spline curve, which is he unique soluion of he variaional problem, and le u denoe any soluion of problem (P c ). Since ỹ is composed by m segmens (which are C p -smooh) and y and u are such ha L ỹ() = L c u() in inerval [a, b], hen u is a piecewise funcion composed by m segmens (which are C p+σ -smooh). Define u i and ỹ i o be he segmens of u and ỹ in he inerval [ i, i+ ]. I is clear ha u i is a soluion of he complee linear differenial equaion L c x = φ i where φ i = L ỹ i. The segmen u i is an elemen of he affine space v + kerl c (where v is any paricular soluion of he differenial equaion L c x = φ i ) and hus u i depends on σ free parameers. Therefore, when σ problem (P c ) has several soluions wih he saed properies. Clearly, here is only one soluion oherwise. Remark () As seen, when σ problem (P c ) has many soluions. To compue a specific soluion one has o fix a oal of mσ real consans. When σ =, he opimal conrol is given by ũ() = c L ỹ(), [a, b]. () In general he opimal soluions are no coninuous in [a, b]. However, when σ = and he incidence vecor Z = (p, z,...,z m, p) is such ha z i < p, he unique opimal soluion is always coninuous. The nex resul shows an imporan propery of he opimal soluions. Theorem 3 If u is an opimal conrol for problem (P c ), hen each segmen u i belongs o he kernel of he linear differenial operaor L L c. Proof. This follows immediaely from he fac ha L L ỹ i = and L ỹ i = L c u i. In paricular, when σ =, ũ i is a soluion of L x =. 5. STATE-SPACE APPROACH We now presen an alernaive formulaion of our opimal conrol problem, which is classical in he sense ha he dynamics are in sae-space form. In his formulaion one expecs o find he opimal conrol wihou previous knowledge of he corresponding oupu. Wih his approach, previous resuls abou he opimal behavior of generalized splines follow easily. The sae-space represenaion (inernal represenaion) corresponding o he inpu-oupu equaion (exernal represenaion) L y() = L c u() is ẋ() = Ax() + Bu(), y() = Cx() + du(), where [a, b] and x() X = R p. This single inpu single oupu linear ime invarian sysem is a minimal realizaion since he operaors L and L c have no common facors. Therefore, wihou loss of generaliy, we may assume ha he marices A, B, C and scalar d are he following:
4 . A =..... a a a p p p B =. C = [ c a c p c p a p c p ] p d = c p. p This paricular sae-space represenaion is clearly conrollable and observable. Remark 3 If x is he firs componen of he sae vecor, he following can be easily derived from he saespace represenaion above: L x () = u(), y() = L c x (). The inpu-oupu equaion L y() = L c u() follows from here immediaely. sense ha hey share he same se of soluions and also produce he same opimal oupu. () We poin ou ha if he operaors L and L c are no relaively prime, hen all he resuls in he previous secion are sill valid. However, in his case, he problems (P c ) and (P cs ) are no longer equivalen. Indeed, he se of opimal conrols of he problem (P cs ) is a proper subse of he se of opimal conrols of he problem (P c ). In spie of his discrepancy, boh problems produce he same oupu, even when he operaors L and L c have common facors. For he paricular case when he incidence vecor is Z = (p,,...,, p) and σ =, he resuls coincide wih hose obained in (Zhang e al., 997) and (Rodrigues e al., 999). Now, le (P cs ) denoe he following opimal conrol problem: b a (L c u()) d u U ẋ() = Ax() + Bu(), y() = Cx() + du(), he inerpolaion condiions y (k) ( i ) = α k i, i =,...,m, k =,...,z i, he boundary condiions y (k) ( ) = η k, y (k) ( m ) = η k m, k =,...,p, and he smoohness condiions min, y (k) ( i ) = y (k) ( i + ), i =,...,m, k =,...,p z i, where α k i, ηk and ηk m are prescribed real numbers. Since we prescribe inerpolaion condiions on he oupu, i would be naural o assume ha he sysem is oupu conrollable, ha is, he marix Γ = [ CB CAB CA p B d ] (p+) has full rank. Neverheless, in our case his follows from he conrollabiliy and observabiliy of he sysem. This opimal conrol problem may be solved direcly by applying he Ponryagin Maximum Principle. Remark 4 () As formulaed, boh opimal conrol problems (P c ) and (P cs ) are equivalen in he 6. EXAMPLES We consider hree examples for p =, he pariion : a = < / < < < 9/4 < 3 = b, and he associaed incidence vecor Z = (,,,,, ). We prescribe he inerpolaion condiions y(/) = 3/, ẏ(/) =, y() =, y() = /, ẏ() = /, y(9/4) = /6, and he boundary condiions y() = 3, ẏ() =, y(3) =, ẏ(3) =. The smoohness condiions are immediaely clear from he incidence vecor above. The fis example is adaped from he opimal conrol problem The cheapes sop of a rain conained in (Agrachev and Sachkov, 4). Example 3 ÿ() = u() (u()) d min, and all he inerpolaion, boundary and smoohness condiions saed above.
5 Noice ha he conrol sysem is given in saespace form by ẋ = x, y = x. ẋ = u, This opimal conrol problem has a unique soluion (because σ = ). Each segmen of he opimal u is a sraigh line and u is clearly coninuous a = and a = 9/4. The opimal oupu is a cubic polynomial in each subinerval [ i, i+ ], alhough no a classical cubic spline on [, 3]. Example 3 ( u() + u()) d min, ÿ() + 36y() = u() + u() and all he inerpolaion, boundary and smoohness condiions saed above. The conrol sysem is given in sandard form by ẋ = x, y = x + x. ẋ = 36x + u, The soluion of he problem is no unique, since σ =. The opimal soluion presened bellow corresponds o a paricular choice of a free parameer in each subinerval. The opimal oupu funcion is in each subinerval a linear combinaion of ψ () = cos(6), ψ () = cos(6), ψ 3 () = sin(6) and ψ 4 () = sin(6). Example 3 3 (ü() u() u()) d min, ACKNOWLEDGEMENT We would like o hank E. Marques de Sá for helpful hins and criicisms ha helped us o clarify some of he ideas conained here. REFERENCES Agrachev, Andrei A. and Yuri L. Sachkov (4). Conrol heory from he geomeric viewpoin. Vol. 87 of Encyclopaedia of Mahemaical Sciences. Springer-Verlag. Berlin. Ahlberg, J. H., E. N. Nilson and J. L. Walsh (964). Fundamenal properies of generalized splines. Proc. Na. Acad. Sci. U.S.A. 5, Kailah, T. (98). Linear sysems. Prenice-Hall. Prenice-Hall Informaion and Sysem Sciences Series. Ribeiro, Maria Isabel (a). Análise de sisemas lineares. Vol. of Colecção ensino da ciência e da ecnologia. IST Press. Ribeiro, Maria Isabel (b). Análise de sisemas lineares. Vol. of Colecção ensino da ciência e da ecnologia. IST Press. Rodrigues, Rui C. and F. Silva Leie (). L- splines a manifesaion of opimal conrol. IMA J. Mah. Conrol Inform. 9(3), Rodrigues, Rui C., F. Silva Leie and C. Simões (999). Generalized splines and opimal conrol. In: Proceedings of he European Conrol Conference, ECC 99. Karlsruhe, Germany. CD-ROM paper F59.PDF. Rugh, Wilson J. (996). Linear sysem heory. nd ed.. Prenice-Hall. Schulz, M. H. and R. S. Varga (967). L-splines. Numer. Mah., Zhang, Zhimin, John Tomlinson and Clyde Marin (997). Splines and linear conrol heory. Aca Appl. Mah. 49(), 34. ÿ() + 9ẏ() y() = ü() u() u() and all he inerpolaion, boundary and smoohness condiions saed above. The conrol sysem is given in sandard form by ẋ = x, y = 8x x +u. ẋ = x 9x + u, The soluion of he problem is no unique, since σ =. The opimal soluion presened bellow corresponds o a paricular choice of free parameer in each subinerval. The opimal oupu funcion is in each subinerval a linear combinaion of exponenial funcions. The opimal curves for hese problems are presened a he end of he paper. 3 y() Fig.. Firs example - opimal oupu
6 u() Fig.. Firs example - opimal inpu u() Fig. 6. Third example - an opimal inpu 3 y() Fig. 3. Second example - opimal oupu 6 u() Fig. 4. Second example - an opimal inpu 3 y() Fig. 5. Third example - opimal oupu
Lecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationKEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationMath 315: Linear Algebra Solutions to Assignment 6
Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationA New Perturbative Approach in Nonlinear Singularity Analysis
Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The
More informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationKEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More information4. Advanced Stability Theory
Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium
More informationExistence of positive solution for a third-order three-point BVP with sign-changing Green s function
Elecronic Journal of Qualiaive Theory of Differenial Equaions 13, No. 3, 1-11; hp://www.mah.u-szeged.hu/ejqde/ Exisence of posiive soluion for a hird-order hree-poin BVP wih sign-changing Green s funcion
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
More informationdy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page
Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,
More informationSingle and Double Pendulum Models
Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double
More informationarxiv: v1 [math.ca] 15 Nov 2016
arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationContinuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.
Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More informationINDEX. Transient analysis 1 Initial Conditions 1
INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationSolutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore
Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s
More informationAnalytical Solutions of an Economic Model by the Homotopy Analysis Method
Applied Mahemaical Sciences, Vol., 26, no. 5, 2483-249 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/.2988/ams.26.6688 Analyical Soluions of an Economic Model by he Homoopy Analysis Mehod Jorge Duare ISEL-Engineering
More informationCONTRIBUTION TO IMPULSIVE EQUATIONS
European Scienific Journal Sepember 214 /SPECIAL/ ediion Vol.3 ISSN: 1857 7881 (Prin) e - ISSN 1857-7431 CONTRIBUTION TO IMPULSIVE EQUATIONS Berrabah Faima Zohra, MA Universiy of sidi bel abbes/ Algeria
More informationCERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS
SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien
More informationUndetermined coefficients for local fractional differential equations
Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More informationTheory of! Partial Differential Equations!
hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationSignal and System (Chapter 3. Continuous-Time Systems)
Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b
More informationNotes for Lecture 17-18
U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up
More informationProperties Of Solutions To A Generalized Liénard Equation With Forcing Term
Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More information11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu
ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning
More informationExistence Theory of Second Order Random Differential Equations
Global Journal of Mahemaical Sciences: Theory and Pracical. ISSN 974-32 Volume 4, Number 3 (22), pp. 33-3 Inernaional Research Publicaion House hp://www.irphouse.com Exisence Theory of Second Order Random
More informationMATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence
MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationChapter Three Systems of Linear Differential Equations
Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n
More informationEssential Maps and Coincidence Principles for General Classes of Maps
Filoma 31:11 (2017), 3553 3558 hps://doi.org/10.2298/fil1711553o Published by Faculy of Sciences Mahemaics, Universiy of Niš, Serbia Available a: hp://www.pmf.ni.ac.rs/filoma Essenial Maps Coincidence
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationTrajectory planning in Cartesian space
Roboics 1 Trajecory planning in Caresian space Prof. Alessandro De Luca Roboics 1 1 Trajecories in Caresian space in general, he rajecory planning mehods proposed in he join space can be applied also in
More informationLecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples
EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,
More informationSystem of Linear Differential Equations
Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y
More informationMathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation
Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion
More informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationOmega-limit sets and bounded solutions
arxiv:3.369v [mah.gm] 3 May 6 Omega-limi ses and bounded soluions Dang Vu Giang Hanoi Insiue of Mahemaics Vienam Academy of Science and Technology 8 Hoang Quoc Vie, 37 Hanoi, Vienam e-mail: dangvugiang@yahoo.com
More informationMean-square Stability Control for Networked Systems with Stochastic Time Delay
JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationOscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,
More informationTheory of! Partial Differential Equations-I!
hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationNEW APPROACH TO DIFFERENTIAL EQUATIONS WITH COUNTABLE IMPULSES
1 9 NEW APPROACH TO DIFFERENTIAL EQUATIONS WITH COUNTABLE IMPULSES Hong-Kun ZHANG Jin-Guo LIAN Jiong SUN Received: 1 January 2007 c 2006 Springer Science + Business Media, Inc. Absrac This paper provides
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationMath 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.
1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be
More informationKinematics and kinematic functions
Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion
More informationIn this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should
Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion
More informationMODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE
Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion
More informationA DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS
A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen
More informationSOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM
SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM FRANCISCO JAVIER GARCÍA-PACHECO, DANIELE PUGLISI, AND GUSTI VAN ZYL Absrac We give a new proof of he fac ha equivalen norms on subspaces can be exended
More informationMath Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.
Mah 50-004 Week 5: Secion 7.4, mass-spring sysems. These are noes for Monday. There will also be course review noes for Tuesday, posed laer. Mon Apr 3 7.4 mass-spring sysems. Announcemens: Warm up exercise:
More informationSolution of Integro-Differential Equations by Using ELzaki Transform
Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationConcourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations
Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationMath 333 Problem Set #2 Solution 14 February 2003
Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial
More informationOnline Appendix to Solution Methods for Models with Rare Disasters
Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,
More informationPOSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION
Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial
More informationThe motions of the celt on a horizontal plane with viscous friction
The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationWEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x
WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile
More informationIMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013
IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher
More informationChapter 7 Response of First-order RL and RC Circuits
Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial
More informationA Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero
JOURNAL OF MAEMAICAL ANALYSIS AND APPLICAIONS 24, 7887 1997 ARICLE NO. AY965143 A Necessary and Sufficien Condiion for he Soluions of a Funcional Differenial Equaion o Be Oscillaory or end o Zero Piambar
More informationLecture Notes 2. The Hilbert Space Approach to Time Series
Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship
More informationdi Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.
di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x
More informationLie Derivatives operator vector field flow push back Lie derivative of
Lie Derivaives The Lie derivaive is a mehod of compuing he direcional derivaive of a vecor field wih respec o anoher vecor field We already know how o make sense of a direcional derivaive of real valued
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More information