Learning-based Adaptive Control for Nonlinear Systems

Size: px
Start display at page:

Download "Learning-based Adaptive Control for Nonlinear Systems"

Transcription

1 MITSUBISHI ELECTRIC RESEARCH LABORATORIES Lerning-bsed Adptive Control for Nonliner Systems Benosmn, M. TR14-64 June 14 Abstrct We present in this pper preliminry result on lerning-bsed dptive trjectory trcking control for nonliner systems. We propose, for the clss of nonliner systems with prmetric uncertinties which cn be rendered integrl Input-to-Stte stble w.r.t. the prmeter estimtion errors input, tht it is possible to merge together the integrl Input-to-Stte stbilizing feedbck controller nd model-free extremum seeking ES lgorithm to relize lerning-bsed dptive controller. We show the efficiency of this pproch on mechtronic exmple. Europen Control Conference ECC 14 This work my not be copied or reproduced in whole or in prt for ny commercil purpose. Permission to copy in whole or in prt without pyment of fee is grnted for nonprofit eductionl nd reserch purposes provided tht ll such whole or prtil copies include the following: notice tht such copying is by permission of Mitsubishi Electric Reserch Lbortories, Inc.; n cknowledgment of the uthors nd individul contributions to the work; nd ll pplicble portions of the copyright notice. Copying, reproduction, or republishing for ny other purpose shll require license with pyment of fee to Mitsubishi Electric Reserch Lbortories, Inc. All rights reserved. Copyright c Mitsubishi Electric Reserch Lbortories, Inc., 14 1 Brodwy, Cmbridge, Msschusetts 139

2

3 14 Europen Control Conference ECC June 4-7, 14. Strsbourg, Frnce Lerning-bsed Adptive Control for Nonliner Systems Mouhcine Benosmn Abstrct We present in this pper preliminry result on lerning-bsed dptive trjectory trcking control for nonliner systems. We propose, for the clss of nonliner systems with prmetric uncertinties which cn be rendered integrl Input-to-Stte stble w.r.t. the prmeter estimtion errors input, tht it is possible to merge together the integrl Input-to-Stte stbilizing feedbck controller nd model-free extremum seeking ES lgorithm to relize lerning-bsed dptive controller. We show the efficiency of this pproch on mechtronic exmple. I. INTRODUCTION Extremum seeking ES is well known pproch by which one cn serch for the extremum of cost function ssocited with given process performnce under some conditions without the need of detiled modelling of the process, e.g. [1], [], [3]. Severl ES lgorithms with their stbility nlysis hve been proposed, e.g. [4], [5], [], [6], [3], [6], [1], [7], [8], nd mny pplictions of ES lgorithms hve been reported, e.g [9], [1], [11], [1], [13]. On the other hnd, clssicl dptive control dels with controlling prtilly unknown processes bsed on their uncertin model, i.e., controlling plnts with prmeters uncertinties. Clssicl dptive methods cn be clssified s direct, where the controller is updted to dpt to the process, or indirect, where the model is updted to better reflect the ctul process. Mny dptive methods hve been proposed over the yers for liner nd nonliner systems, we could not possibly cite here ll the design nd nlysis results tht hve been reported, insted we refer the reder to e.g. [14], [15] nd the references therein for more detils. Wht we wnt to underline here is tht these results in clssicl dptive control re minly bsed on the structure of the model of the system, e.g. liner vs. nonliner, with liner uncertinties prmetriztion vs. nonliner prmeteriztions, etc. Another dptive control prdigm is the one which uses lerning schemes to estimte the uncertin prt of the process. Indeed, in this prdigm the lerning-bsed controller, bsed either on mchine lerning theory, neurl networks, fuzzy systems, etc. is trying either to estimte the prmeters of n uncertin model, or the structure of deterministic or stochstic function representing prt or Mouhcine Benosmn m benosmn@ieee.org is with Mitsubishi Electric Reserch Lbortories, Cmbridge, MA 139, USA. totlity of the model. Severl results hve been proposed in this re s well, nd we refer the reder to e.g. [16] nd the references therein for more detils. We wnt to concentrte in this pper on the use of ES theory in the lerning-bsed dptive control prdigm. Indeed, severl results were recently developed in this direction, e.g. [17], [18], [19], [], [9], [1], [1], [13]. For instnce in [17], [18] the uthors used model-free ES, i.e., only bsed on desired cost function, to estimte prmeters of liner stte feedbck to compenste for unknown prmeters for liner systems. In [19], [] n extremum seeking-bsed controller for nonliner ffine systems with liner prmeters uncertinties ws proposed. The controller drives the sttes of the system to unknown optiml sttes tht optimize desired objective function. The ES controller used in [19], [] is not model-free in the sense tht it is bsed on the known prt of the model, i.e., it is designed bsed on the objective function nd the nonliner model structure. A similr pproch is used in [9], [1] when deling with more specific exmples. In [1], the uthors used for the cse of electromgnetic ctutors, model-free ES, i.e., only bsed on the cost function without the use of the systems model, to lern the best feedbck gins of pssive robust stte feedbck. Similrly, in [1], bckstepping controller ws merged with model-free ES to estimte the uncertin prmeters of nonliner model for electromgnetic ctutors. Although, no stbility nlysis ws presented for the full controller i.e., bckstepping plus ES estimtor, very promising numericl results where reported. In this work we propose to generlize the ide of [1], for the clss of nonliner system with prmetric uncertinties which cn be rendered iiss w.r.t. the prmeter estimtion errors. The ide is bsed on modulr design, where we first design feedbck controller which mkes the closed-loop trcking error dynmic iiss w.r.t. the estimtion errors nd then complement this iiss-controller with mode-free ES lgorithm tht cn minimize desired cost function, by tuning, i.e., estimting, the unknown prmeters of the model. The modulr design simplifies the nlysis of the totl controller, i.e., iiss-controller plus ES estimtion lgorithm. We first propose this formultion in the generl cse of nonliner systems nd then show cse-study on mechtronic exmple. This pper is orgnized s follows: Section II is used to EUCA 9

4 recll some nottions nd definitions. In Section III we present the min result of this pper, nmely the ES-bsed lerning dptive controller. Section IV is dedicted to n ppliction exmple, nd the pper ends with Conclusion in Section V. II. PRELIMINARIES Throughout the pper we will use. to denote the Eucliden norm; i.e., for x R n we hve x = x T x. We will use. for the short nottion of time derivtive. We denote by C k functions tht re k times differentible. A function is sid nlytic in given set, if it dmits convergent Tylor series pproximtion in some neighborhood of every point of the set. A continuous function α : [, [, is sid to belong to clss K if it is strictly incresing nd α =. A continuous function β : [, [, [, is sid to belong to clss KL if, for ech fixed s, the mpping βr, s belongs to clss K with respect to r nd, for ech fixed r, the mpping βr, s is decresing with respect to s nd βr, s s s. Let us now introduce the definition of Locl Integrl Inputto-Stte Stbility iiss. Definition 1 Locl Integrl Input-to-Stte Stbility [1]: Consider the system ẋ = ft, x, u 1 where x D R n such tht D, nd f : [, D D u R n is piecewise continuous in t nd loclly Lipschitz in x nd u, uniformly in t. The inputs re ssumed to be mesurble nd loclly bounded functions u : R D u R m. Given ny control u D u nd ny ξ D D, there is unique mximl solution of the initil vlue problem ẋ = ft, x, u, xt = ξ. Without loss of generlity, ssume t =. The unique solution is defined on some mximl open intervl, nd it is denoted by x, ξ, u. System 1 is loclly integrl input-to-stte stble LiISS if there exist functions α, γ K nd β KL such tht, for ll ξ D nd ll u D u, the solution xt, ξ, u is defined for ll t nd α xt, ξ, u β ξ, t + t γ us ds for ll t. Equivlently, system 1 is LiISS if nd only if there exist functions β KL nd γ 1, γ K such tht t xt, ξ, u β ξ, t + γ 1 γ us ds 3 for ll t, ll ξ D nd ll u D u. III. LEARNING-BASED ADAPTIVE CONTROLLER Consider the system 1, with n dditionl rgument representing prmetric uncertinties R p ẋ = ft, x,, u 4 We ssocite with 4, the output vector y = hx 5 where h : R n R h. The control objective here is for y to symptoticlly trck desired smooth vector time-dependent trjectory y ref : [, R h. Let us now define the output trcking error vector s e y t = yt y ref t. We then ssume the following Assumption 1: There exists robust control feedbck u iss t, x, ˆ : R R n R p R m, with ˆ being the dynmic estimte of the uncertin vector, such tht, the closed-loop error dynmics ė y = f ey t, e y, e 6 is iiss from the input vector e = ˆ to the stte vector e y. Remrk 1: Assumption 1 might seem too generl, however, severl control pproches cn be used to design controller u iss rendering n uncertin system iiss, for instnce bckstepping control pproch hs been shown to chieve such property for prmetric strict-feedbck systems, e.g. [15]. This is preliminry report, nd we do not pretend here to present detiled solution for ll the cses. A more detiled study of how to chieve Assumption 1 for specific clsses of systems nd how to use it in the context of ES lerning-bsed dptive control, will be presented in our future reports. Let us define now the following cost function Q ˆ = F e y ˆ 7 where F : R h R, F =, F e y > for e y. We need the following ssumptions on Q. Assumption : The cost function Q hs locl minimum t ˆ =. Assumption 3: The initil error e t is sufficiently smll, i.e., The originl prmeters estimtes vector ˆ is close enough to the ctul prmeters vector. Assumption 4: The cost function is nlytic nd its vrition with respect to the uncertin vribles is bounded in the neighborhood of, i.e., Q ξ, ξ >, V, where V denotes compct neighborhood of. Remrk : Assumption simply mens tht we cn consider tht Q hs t lest locl minimum t the true vlues of the uncertin prmeters. Remrk 3: Assumption 3 indictes tht our result will be of locl nture, mening tht our nlysis holds in smll neighborhood of the ctul vlues of the prmeters. We cn now present the following Lemm. 91

5 Lemm 1: Consider the system 4, 5, with the cost function 7, then under Assumptions 1,, 3, nd 4, the controller u iss, where ˆ is estimted with the multi-prmeter extremum seeking lgorithm ẋ i = i sinω i t + π Q ˆ ˆ i = x i + i sinω i t π, i {1,..., p} 8 with ω i ω j, ω i + ω j ω k, i, j, k {1,..., p}, nd ω i > ω, i {1,..., p}, with ω lrge enough, ensures tht the norm of the error vector e y dmits the following bound e yt β e y, t+α t γ β e, t+ e mxds, i=p where e mx = ξ1 ω + i=1 i, ξ 1, ξ >, e D e, ω = mx i {1,...,p} ω i, α K, β KL, β KL nd γ K. Proof: Consider the system 4, 5, then under Assumption 1, the controller u iss ensures tht the trcking error dynmic 6 is iiss between the input e nd the stte vector e y, which by Definition 1, implies tht there exist functions α K, β KL nd γ K, such tht, for ll e D e nd e D e, the norm of the error vector e dmits the following bound e y t β e y, t + α t γ e ds 9 for ll t. Now, we need to evlute the bound on the estimtion vector, to do so we use the results presented in []. First, bsed on Assumption 4, the cost function is loclly Lipschitz, i.e., η 1 >, s.t. Q 1 Q η 1 1, 1, V. Furthermore, since Q is nlytic it cn be pproximted loclly in V with qudrtic function, e.g. Tylor series up to second order. Bsed on this nd on Assumptions nd 3, we cn write the following bound [], pges : e t dt e t dt β e, t + ξ1 e t β e, t + ξ1 ω e t β e, t + ξ1 ω + + dt i=p i=1 i, with β KL, ξ 1 >, t, ω = mx i {1,...p} ω i, dt = [ 1 sinω 1 t + π,..., psinω p t + π ]T. which together with the bound 9 completes the proof. Remrk 4: The estimted prmeters upper bounds used in Lemm 1 re correlted to the choice of the first order multi-vrible extremum seeking MES 7 nd 8, however, these bounds cn be esily chnged by using other MES lgorithms, e.g. [3], [4], which is due to the modulr design of the controller, tht uses the iiss robust prt to ensure boundedness of the error dynmics nd the lerning prt to improve the trcking performnce. Remrk 5: We point out here tht ISS cn be substituted for iiss if we re deling with time-invrint systems nd ω, solving regultion problem insted of time-vrying trjectory trcking. A. Controller design IV. CASE STUDY We study here the exmple of electromgnetic ctutors modelled by the nonliner equtions [4] m d x dt = kx x η dx u = Ri + di b+x dt dt i b+x + f d i dx b+x dt, x x f 1 where, x represents the rmture position physiclly constrined between the initil position of the rmture, nd the mximl position of the rmture x f, dx dt represents the rmture velocity, m is the rmture mss, k the spring constnt, x the initil spring length, η the dmping coefficient, b+x i represents the electromgnetic force EMF generted by the coil,, b being constnt prmeters of the coil, f d constnt term modelling unknown disturbnce force, e.g. sttic friction, R the resistnce of the coil, L = b+x the coil inductnce ssumed to be dependent on the i dx position of the rmture, b+x dt represents the bck EMF. Finlly, i denotes the coil current, di dt its time derivtive nd u represents the control voltge pplied to the coil. We consider there the control problem of the electromgnetic system ssuming uncertinties on the spring constnt k, the dmping coefficient η, nd the dditive disturbnce f d. Let us define the stte vector z := [z 1 z z 3 ] T = [x ẋ i] T. The objective of the control is to mke the vribles z 1, z robustly trck sufficiently smooth i.e. C time-vrying position nd velocity trjectories z ref 1 t, z ref t = dzref 1 t dt tht stisfy the following constrints: z ref 1 t = z 1int, z ref 1 t f = z 1f, ż ref 1 t = ż ref 1 t f =, z ref 1 t = z ref 1 t f =, where t is the strting time of the trjectory, t f is the finl time, z 1int is the initil position nd z 1f is the finl position. To strt, we first write the system 1 in the following form ż 1 = z ż = k m x z 1 η m z ż 3 = R b+z 1 z 3 + z3 b+z 1 z + mb+z 1 z 3 11 u b+z 1 The uthors in [5] hve shown tht under the following 9

6 feedbck controller u= Rb+z 1 b+z 1 z 3 z z 3 b+z z 3 mb+z 1 z zref c z 3 ũ + mz ˆk z 3 m x z1 ˆη m z+ ˆf dm +c 3z 1 z ref 1 +c 1z z ref ż ref +κ 1z z ref ψ ż ref + mb+z 1 z 3 ˆk c1+κ 1 ψ m ˆη ˆη m żref z ˆk x z 1 z m + κ z 3 ũ mb+z 1 z 3 m x z1 ˆη m z+ f ˆ dm mb+z 1 z 3 + mb+z 1 z 3 κ 1z z ref z 3 mb+z 1 m m x z 1 ˆη m z + ˆf dm [ c1+κ 1 ψ ˆη m + κ1z z ref ] z m κ 3z 3 ũ mb+z 1 z ψ 3 +mb+z 1 z ˆk 3 m z zref +c 3z z ref with ũ= mb+z 1 +c 1z z ref ż ref ˆk m x z1 ˆη m z+ ˆf dm +c 3z 1 z ref 1 + mb+z 1 κ 1z z ref ψ, ψ, 1 13 where ˆk, ˆη, ˆfd re the system prmeter estimtes, nd ψ x z1 m, z m, 1 m T ; the system 1, 1 nd 13, is loclly iiss- LiISS. Next, we define the cost function tf tf Q ˆ = q 1z 1s z 1s ref ds+ q z s z ref s ds, 14 where q 1, q >, nd ˆ = ˆ k, ˆ η, ˆ fd T represents the vector of the lerned prmeters, defined such tht ˆkt = k nominl + ˆ k t ˆηt = η nominl + ˆ η t ˆf d t = f d nominl + ˆ fd t 15 where k nominl, η nominl, f d nominl re the nominl vlues of the prmeters, nd the ˆ i s re computed using discrete version of 8, given by x k k + 1 = x k k + k t f sinω k k t f + π Q ˆ k k + 1 = x k k k sinω k k t f π, x η k + 1 = x η k + η t f sinω η k t f + π Q ˆ η k + 1 = x η k η sinω η k t f π, x fd k + 1 = x fd k + fd t f sinω fd k t f + π Q ˆ fd k + 1 = x fd k fd sinω fd k t f π, 16 eventully, we conclude bsed on Lemm 1, tht the controller 1, 13, 15 nd 16, ensures tht the trcking error norm is bounded with decresing function of the estimtion error. B. Numericl results In this section, we illustrte our pproch for the nonliner electromgnetic ctutor modelled by 1, using the system prmeters given in Tble I [6]. The reference trjectory is designed to be 5 th order polynomil, x ref t = 5 i= i t t f i where the coefficients i re selected such tht the following conditions re stisfied: 93 Prmeter m R η x k b Vlue.7 [kg] 6 [Ω] 7.53 [ kg s ] 8 [mm] 158 [ N mm ] [ Nm [m] TABLE I SYSTEM PARAMETER VALUES A ] x ref =. mm, x ref.5 =.7 mm, ẋ ref =, ẋ ref.5 =, ẍ ref =, ẋ ref.5 =. We consider the uncertinties given by k = 4.5, η =.7 nd f d = 7.5. To mke the simultion cse more chllenging we lso introduced n initil error x =.1 mm on the rmture position. We implemented the controller 1 nd 13 with the coefficients c 1 = 1, c = 1, c 3 = 5, κ 1 = κ = κ 3 =.5, together with the lerning lgorithm 14, 15 nd 16 with the coefficients k =.5, ω k = 7.5, η =., ω η = 7.4, fd = 1, ω fd = 7.3, q 1 = q = 1. For more detils bout the tuning of the MES coefficients we refer the reder to [], [], [1], however, we underline here tht the frequencies ω i, ı = 1,, 3 hve been selected high enough to ensure efficient explortion on the serch spce nd ensure convergence nd tht the mplitudes i, i = 1,, 3 of the dither signls, hve been chosen such tht the serch is fst enough for this ppliction. Here due to the cyclic nture of the problem, i.e., cyclic motion of the rmture between nd x f, the uncertin prmeters estimte vector ˆk, ˆη, ˆfd T is updted for ech cycle, i.e., t the end of ech cycle t t = t f, the cost function Q is updted, nd the new estimte of the prmeters is computed for the next cycle. The purpose of using MES scheme long with iiss-bckstepping controller is to improve the performnce of the iiss-bckstepping controller by better estimting the system prmeters over mny cycles, hence decresing the error in the prmeters over time to provide better trjectory following for the ctutor. As cn be seen in Figures 1 nd 1b, the robustifiction of the bckstepping control vi extremum seeking gretly improves the trcking performnce. Figure b shows tht the cost function decreses below 1 within itertions. It cn be seen in Figure tht the cost strts t n initil vlue round 9, nd decreses rpidly fterwrds. Moreover, the estimted prmetric uncertinties ˆ k, ˆ η nd ˆ fd converge to regions round the ctul prmeter vlues, s shown on Figure 3. The number of itertions for the estimte to rech the ctul vlue of the prmeters my pper to be high. The reson behind tht is tht the llowed uncertinties in the prmeters re lrge, hence the extremum seeking

7 position [mm] x w/o ESC x ref x w ESC Cost function Q time [s] Obtined Armture Position vs. Reference Trjectory Number of itertions Cost function- zoom 1.5 v w/o ESC v ref v w ESC 5 velocity [mm/s] 1.5 Cost function Q time [s] b Obtined Armture Velocity vs. Reference Trjectory Number of itertions b Cost function Fig. 1. Obtined trjectories vs. Reference Trjectory- Cse with uncertin k, η, f d Fig.. Cost function- Cse with uncertin k, η, f d scheme requires lot of itertions to improve performnce. Furthermore, we purposely tested the chllenging cse of three simultneous uncertinties, which mkes the spce serch for the lerning lgorithm lrge note tht this cse of multiple uncertinties could not be solved with other clssicl model-bsed dptive controller [7], due to some intrinsic limittions of model-bsed dptive controller. However, In rel-life pplictions uncertinties ccumulte grdully over long period of time, while the lerning lgorithm keeps trcking these chnges continuously. Thus, the extremum seeking lgorithm will be ble to improve the controller performnce in much fewer itertions. Finlly, the control voltge is depicted on Figure 4, which shows n initil high vlue due to the reltively lrge simulted initil condition error on the rmture position. V. CONCLUSION In this pper we hve studied the problem of lerningbsed dptive control for nonliner systems with prmetric uncertinties. We rgued tht for the clss of nonliner systems which cn be rendered iiss w.r.t. the prmeter estimtion errors, by robust feedbck controller, it is possible to combine the iiss feedbck controller with model-free ES lgorithm to obtin lerning-bsed dptive controller. We showed detiled ppliction of this pproch on mechtronic exmple nd reported encourging numericl results. In this preliminry pper, we introduced the ide in generl setting, however, further investigtion re needed to nlyze specific nonliner systems clsses which cn be stbilized in the iiss sense w.r.t. to the estimtion error of the uncertin prmeters, nd show for these specific clsses constructive control design pproch, in the context of the lerning-bsed dptive control presented here. Further work will lso del with using different ES lgorithms with less restrictive conditions on the dither signls mplitude nd frequencies, e.g.[8], [4], together with investigting other model-free lerning lgorithms such us reinforcement lerning lgorithms, nd compring the obtined controllers to the vilble clssicl dptive controllers. REFERENCES [1] K. Ariyur nd M. Krstić, Rel-time optimiztion by extremum-seeking control. Wiley-Blckwell, 3. [] K. B. Ariyur nd M. Krstic, Multivrible extremum seeking feedbck: Anlysis nd design, in Proc. of the Mthemticl Theory of Networks nd Systems, South Bend, IN, August. [3] D. Nesic, Extremum seeking control: Convergence nlysis, Europen Journl of Control, vol. 15, no. 34, pp , 9. 94

8 ˆ η ˆ k Number of itertions Prmeter k estimte Number of itertions ˆ f d b Prmeter η estimte Number of itertions Fig. 3. U [v] c Prmeter f d estimte Prmeters estimtes- Cse with uncertin k, η, f d time [s] Fig. 4. Control voltge- Cse with uncertin k, η, f d [4] A. Scheinker, Simultneous stbiliztion of nd optimiztion of unkown time-vrying systems, in Americn Control Conference, June 13, pp [5] M. Krstic, Performnce improvement nd limittions in extremum seeking, Systems & Control Letters, vol. 39, pp ,. [6] Y. Tn, D. Nesic, nd I. Mreels, On non-locl stbility properties of extremum seeking control, Automtic, no. 4, pp , 6. [7] M. Rote, Anlysis of multivrible extremum seeking lgorithms, in Proceedings of the Americn Control Conference, vol. 1, no. 6. IEEE,, pp [8] M. Guy, S. Dhliwl, nd D. Dochin, A time-vrying extremumseeking control pproch, in Americn Control Conference, 13, pp [9] T. Zhng, M. Guy, nd D. Dochin, Adptive extremum seeking control of continuous stirred-tnk biorectors, AIChE J., no. 49, p , 3. [1] N. Hudon, M. Guy, M. Perrier, nd D. Dochin, Adptive extremumseeking control of convection-rection distributed rector with limited ctution, Computers & Chemicl Engineering, vol. 3, no. 1, pp , 8. [11] C. Zhng nd R. Ordez, Extremum-Seeking Control nd Applictions. Springer-Verlg, 1. [1] M. Benosmn nd G. Atinc, Multi-prmetric extremum seekingbsed lerning control for electromgnetic ctutors, in Americn Control Conference, 13, pp [13], Nonliner lerning-bsed dptive control for electromgnetic ctutors, in Europen Control Conference, 13, pp [14] I. D. Lndu, R. Lozno, M. M Sd, nd A. Krimi, Adptive Control Algorithms, Anlysis nd Applictions, ser. Communictions nd Control Engineering. Springer-Verlg, 11. [15] M. Krstic, I. Knellkopoulos, nd P. Kokotovic, Nonliner nd dptive control design. John Wiley & Sons New York, [16] C. Wng nd D. Hill, Deterministic Lerning Theory for Identifiction, Recognition, nd Control, ser. Automtion nd control engineering series. Tylor & Frncis Group, 6. [17] P. Hghi nd K. Ariyur, On the extremum seeking of model reference dptive control in higher-dimensionl systems, in Americn Control Conference, 11, pp [18] K. B. Ariyur, S. Gnguli, nd D. F. Enns, Extremum seeking for model reference dptive control, in Proc. of the AIAA Guidnce, Nvigtion, nd Control Conference, 9, doi: 1.514/ [19] M. Guy nd T. Zhng, Adptive extemum seeking control of nonliner dynmic systems with prmetric uncertinties, Automtic, no. 39, p , 3. [] V. Adetol nd M. Guy, Prmeter convergence in dptive extremum-seeking control, Automtic, no. 43, p. 1511, 7. [1] H. Ito nd Z. Jing, Necessry nd sufficient smll gin conditions for integrl input-to-stte stble systems: A Lypunov perspective, IEEE Trnsctions on Automtic Control, vol. 54, no. 1, pp , 9. [] M. A. Rote, Anlysis of multivrible extremum seeking lgorithms, in Americn Control Conference, June, pp [3] W. Nose, Y. Tn, D. Nesic, nd C. Mnzie, Non-locl stbility of multi-vrible extremum-seeking scheme, in IEEE, Austrlin Control Conference, November 11, pp [4] K. Peterson nd A. Stefnopoulou, Extremum seeking control for soft lnding of electromechnicl vlve ctutor, Automtic, vol. 4, pp , 4. [5] G. Atinc nd M. Benosmn, Nonliner lerning-bsed dptive control for electromgnetic ctutors with proof of stbility, in IEEE, Conference on Decision nd Control, 13, pp [6] N. Khveci nd I. Kolmnovsky, Control design for electromgnetic ctutors bsed on bckstepping nd lnding reference governor, in 5th IFAC Symposium on Mechtronic Systems, Cmbridge, September 1, pp [7] M. Benosmn nd G. Atinc, Nonliner dptive control of electromgnetic ctutors, in SIAM Conference on Control nd Applictions, 13, pp

Multi-parametric Extremum Seeking-based Learning Control for Electromagnetic Actuators

Multi-parametric Extremum Seeking-based Learning Control for Electromagnetic Actuators MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Multi-prmetric Extremum Seeking-bsed Lerning Control for Electromgnetic Actutors Benosmn, M.A.; Atinc, G.M. TR213-137 June 213 Abstrct We study

More information

Extremum Seeking-based Indirect Adaptive Control for Nonlinear Systems

Extremum Seeking-based Indirect Adaptive Control for Nonlinear Systems MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Extremum Seeking-based Indirect Adaptive Control for Nonlinear Systems Benosman, M. TR14-85 August 14 Abstract We present in this paper a preliminary

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

On Stability of Nonlinear Slowly Time-Varying and Switched Systems

On Stability of Nonlinear Slowly Time-Varying and Switched Systems 2018 IEEE Conference on Decision nd Control CDC) Mimi Bech, FL, USA, Dec. 17-19, 2018 On Stbility of Nonliner Slowly Time-Vrying nd Switched Systems Xiobin Go, Dniel Liberzon, nd Tmer Bşr Abstrct In this

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

NUMERICAL INTEGRATION

NUMERICAL INTEGRATION NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

Taylor Polynomial Inequalities

Taylor Polynomial Inequalities Tylor Polynomil Inequlities Ben Glin September 17, 24 Abstrct There re instnces where we my wish to pproximte the vlue of complicted function round given point by constructing simpler function such s polynomil

More information

Generation of Lyapunov Functions by Neural Networks

Generation of Lyapunov Functions by Neural Networks WCE 28, July 2-4, 28, London, U.K. Genertion of Lypunov Functions by Neurl Networks Nvid Noroozi, Pknoosh Krimghee, Ftemeh Sfei, nd Hmed Jvdi Abstrct Lypunov function is generlly obtined bsed on tril nd

More information

1 Online Learning and Regret Minimization

1 Online Learning and Regret Minimization 2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in

More information

1.9 C 2 inner variations

1.9 C 2 inner variations 46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for

More information

Reinforcement learning II

Reinforcement learning II CS 1675 Introduction to Mchine Lerning Lecture 26 Reinforcement lerning II Milos Huskrecht milos@cs.pitt.edu 5329 Sennott Squre Reinforcement lerning Bsics: Input x Lerner Output Reinforcement r Critic

More information

Linear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System

Linear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System Pure nd Applied Mthemtics Journl 017; 6(1): 5-13 http://www.sciencepublishinggroup.com/j/pmj doi: 10.11648/j.pmj.0170601.1 ISSN: 36-9790 (Print); ISSN: 36-981 (Online) Liner nd Non-liner Feedbck Control

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Czechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction

Czechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCK-KURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When rel-vlued

More information

Department of Mechanical Engineering MECE 551 Final examination Winter 2008 April 16, 9:00 11:30. Question Value Mark

Department of Mechanical Engineering MECE 551 Final examination Winter 2008 April 16, 9:00 11:30. Question Value Mark Deprtment of Mechnicl Engineering MECE 55 Finl exmintion Winter 8 April 6, 9: :3 Notes: You my hve your text book nd one pge formul sheet Electronic devices re not llowed except n pproved clcultor NAME:

More information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information

A Bernstein polynomial approach for solution of nonlinear integral equations

A Bernstein polynomial approach for solution of nonlinear integral equations Avilble online t wwwisr-publictionscom/jns J Nonliner Sci Appl, 10 (2017), 4638 4647 Reserch Article Journl Homepge: wwwtjnscom - wwwisr-publictionscom/jns A Bernstein polynomil pproch for solution of

More information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

Improved Results on Stability of Time-delay Systems using Wirtinger-based Inequality

Improved Results on Stability of Time-delay Systems using Wirtinger-based Inequality Preprints of the 9th World Congress he Interntionl Federtion of Automtic Control Improved Results on Stbility of ime-dely Systems using Wirtinger-bsed Inequlity e H. Lee Ju H. Prk H.Y. Jung O.M. Kwon S.M.

More information

The final exam will take place on Friday May 11th from 8am 11am in Evans room 60.

The final exam will take place on Friday May 11th from 8am 11am in Evans room 60. Mth 104: finl informtion The finl exm will tke plce on Fridy My 11th from 8m 11m in Evns room 60. The exm will cover ll prts of the course with equl weighting. It will cover Chpters 1 5, 7 15, 17 21, 23

More information

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Partial Derivatives. Limits. For a single variable function f (x), the limit lim Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles

More information

Non-linear Adaptive Control for Electromagnetic Actuators

Non-linear Adaptive Control for Electromagnetic Actuators MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.erl.co Non-liner Adptive Control for Electrognetic Actutors Benosn, M.; Atinc, G.M. TR15-4 Jnury 15 Abstrct The uthors study here the proble of dptive

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

Introduction to the Calculus of Variations

Introduction to the Calculus of Variations Introduction to the Clculus of Vritions Jim Fischer Mrch 20, 1999 Abstrct This is self-contined pper which introduces fundmentl problem in the clculus of vritions, the problem of finding extreme vlues

More information

KRASNOSEL SKII TYPE FIXED POINT THEOREM FOR NONLINEAR EXPANSION

KRASNOSEL SKII TYPE FIXED POINT THEOREM FOR NONLINEAR EXPANSION Fixed Point Theory, 13(2012), No. 1, 285-291 http://www.mth.ubbcluj.ro/ nodecj/sfptcj.html KRASNOSEL SKII TYPE FIXED POINT THEOREM FOR NONLINEAR EXPANSION FULI WANG AND FENG WANG School of Mthemtics nd

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

THIELE CENTRE. Linear stochastic differential equations with anticipating initial conditions

THIELE CENTRE. Linear stochastic differential equations with anticipating initial conditions THIELE CENTRE for pplied mthemtics in nturl science Liner stochstic differentil equtions with nticipting initil conditions Nrjess Khlif, Hui-Hsiung Kuo, Hbib Ouerdine nd Benedykt Szozd Reserch Report No.

More information

ADAPTIVE CONTROL AND SYNCHRONIZATION OF A HIGHLY CHAOTIC ATTRACTOR

ADAPTIVE CONTROL AND SYNCHRONIZATION OF A HIGHLY CHAOTIC ATTRACTOR ADAPTIVE CONTROL AND SYNCHRONIZATION OF A HIGHLY CHAOTIC ATTRACTOR Sundrpndin Vidynthn 1 1 Reserch nd Development Centre, Vel Tech Dr. RR & Dr. SR Technicl University Avdi, Chenni-600 062, Tmil Ndu, INDIA

More information

Adaptive Stabilization of Nonholonomic Mobile Robots with Unknown Kinematic Parameters

Adaptive Stabilization of Nonholonomic Mobile Robots with Unknown Kinematic Parameters WSEAS TRANSACTIONS on SYSTEMS nd CONTROL Deheng Hou, Fngzheng Go Adptive Stbiliztion of Nonholonomic Mobile Robots with Unknown Kinemtic Prmeters DEHENG HOU FANGZHENG GAO School of Softwre Engineering

More information

Best Approximation. Chapter The General Case

Best Approximation. Chapter The General Case Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Frobenius numbers of generalized Fibonacci semigroups

Frobenius numbers of generalized Fibonacci semigroups Frobenius numbers of generlized Fiboncci semigroups Gretchen L. Mtthews 1 Deprtment of Mthemticl Sciences, Clemson University, Clemson, SC 29634-0975, USA gmtthe@clemson.edu Received:, Accepted:, Published:

More information

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space. Clculus 3 Li Vs Spce Curves Recll the prmetric equtions of curve in xy-plne nd compre them with prmetric equtions of curve in spce. Prmetric curve in plne x = x(t) y = y(t) Prmetric curve in spce x = x(t)

More information

ON PREY-PREDATOR MODEL WITH HOLLING-TYPE II AND LESLIE-GOWER SCHEMES AHMED BUSERI ASHINE

ON PREY-PREDATOR MODEL WITH HOLLING-TYPE II AND LESLIE-GOWER SCHEMES AHMED BUSERI ASHINE Journl of Mthemtics nd Computer Applictions Reserch (JMCAR) ISSN(P): 5-48; ISSN(E): Applied Vol 3, Issue 1, Dec 16, 7-34 TJPRC Pvt Ltd ON PREY-PREDATOR MODEL WITH HOLLING-TYPE II AND LESLIE-GOWER SCHEMES

More information

Numerical Analysis: Trapezoidal and Simpson s Rule

Numerical Analysis: Trapezoidal and Simpson s Rule nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =

More information

Deteriorating Inventory Model for Waiting. Time Partial Backlogging

Deteriorating Inventory Model for Waiting. Time Partial Backlogging Applied Mthemticl Sciences, Vol. 3, 2009, no. 9, 42-428 Deteriorting Inventory Model for Witing Time Prtil Bcklogging Nit H. Shh nd 2 Kunl T. Shukl Deprtment of Mthemtics, Gujrt university, Ahmedbd. 2

More information

Modeling and Controller Design for the Air-to-Air Missile Uncertain System

Modeling and Controller Design for the Air-to-Air Missile Uncertain System Interntionl ournl of Computer (IC) ISSN 07-45 (Print & Online) Globl Society of Scientific Reserch nd Reserchers http://ijcjournl.org/ Modeling nd Controller Design for the Air-to-Air Missile Uncertin

More information

Reinforcement Learning

Reinforcement Learning Reinforcement Lerning Tom Mitchell, Mchine Lerning, chpter 13 Outline Introduction Comprison with inductive lerning Mrkov Decision Processes: the model Optiml policy: The tsk Q Lerning: Q function Algorithm

More information

Some estimates on the Hermite-Hadamard inequality through quasi-convex functions

Some estimates on the Hermite-Hadamard inequality through quasi-convex functions Annls of University of Criov, Mth. Comp. Sci. Ser. Volume 3, 7, Pges 8 87 ISSN: 13-693 Some estimtes on the Hermite-Hdmrd inequlity through qusi-convex functions Dniel Alexndru Ion Abstrct. In this pper

More information

dt. However, we might also be curious about dy

dt. However, we might also be curious about dy Section 0. The Clculus of Prmetric Curves Even though curve defined prmetricly my not be function, we cn still consider concepts such s rtes of chnge. However, the concepts will need specil tretment. For

More information

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer 5 Proceedings of the Interntionl Conference on Informtion nd Automtion, December 58, 5, Colombo, Sri Ln. Appliction of Exct Discretiztion for Logistic Differentil Equtions to the Design of Discrete-ime

More information

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient * Interntionl Mthemticl Forum, 4, 9, no., 7-3 Appliction of Exp-Function Method to Huxley Eqution with Vrible Coefficient * Li Yo, Lin Wng nd Xin-Wei Zhou. Deprtment of Mthemtics, Kunming College Kunming,Yunnn,

More information

Ordinary differential equations

Ordinary differential equations Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Passive Control of Mechanical Systems

Passive Control of Mechanical Systems 1 Perry Y. Li Deprtment of Mechnicl Engineering University of Minnesot Minnepolis MN 55455 pli@me.umn.edu http://www.me.umn.edu/~pli/ 2 Outline Introduction / Motivtion Pssive Velocity Field Control (PVFC)

More information

An approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95

An approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95 An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech

More information

Positive Solutions of Operator Equations on Half-Line

Positive Solutions of Operator Equations on Half-Line Int. Journl of Mth. Anlysis, Vol. 3, 29, no. 5, 211-22 Positive Solutions of Opertor Equtions on Hlf-Line Bohe Wng 1 School of Mthemtics Shndong Administrtion Institute Jinn, 2514, P.R. Chin sdusuh@163.com

More information

Recitation 3: More Applications of the Derivative

Recitation 3: More Applications of the Derivative Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech

More information

AMATH 731: Applied Functional Analysis Fall Some basics of integral equations

AMATH 731: Applied Functional Analysis Fall Some basics of integral equations AMATH 731: Applied Functionl Anlysis Fll 2009 1 Introduction Some bsics of integrl equtions An integrl eqution is n eqution in which the unknown function u(t) ppers under n integrl sign, e.g., K(t, s)u(s)

More information

5.3 Nonlinear stability of Rayleigh-Bénard convection

5.3 Nonlinear stability of Rayleigh-Bénard convection 118 5.3 Nonliner stbility of Ryleigh-Bénrd convection In Chpter 1, we sw tht liner stbility only tells us whether system is stble or unstble to infinitesimlly-smll perturbtions, nd tht there re cses in

More information

Chapter 3 Solving Nonlinear Equations

Chapter 3 Solving Nonlinear Equations Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

1 The Lagrange interpolation formula

1 The Lagrange interpolation formula Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x

More information

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b. Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn

More information

Realistic Method for Solving Fully Intuitionistic Fuzzy. Transportation Problems

Realistic Method for Solving Fully Intuitionistic Fuzzy. Transportation Problems Applied Mthemticl Sciences, Vol 8, 201, no 11, 6-69 HKAR Ltd, wwwm-hikricom http://dxdoiorg/10988/ms20176 Relistic Method for Solving Fully ntuitionistic Fuzzy Trnsporttion Problems P Pndin Deprtment of

More information

CBE 291b - Computation And Optimization For Engineers

CBE 291b - Computation And Optimization For Engineers The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn

More information

Research Article Moment Inequalities and Complete Moment Convergence

Research Article Moment Inequalities and Complete Moment Convergence Hindwi Publishing Corportion Journl of Inequlities nd Applictions Volume 2009, Article ID 271265, 14 pges doi:10.1155/2009/271265 Reserch Article Moment Inequlities nd Complete Moment Convergence Soo Hk

More information

Time Optimal Control of the Brockett Integrator

Time Optimal Control of the Brockett Integrator Milno (Itly) August 8 - September, 011 Time Optiml Control of the Brockett Integrtor S. Sinh Deprtment of Mthemtics, IIT Bomby, Mumbi, Indi (emil : sunnysphs4891@gmil.com) Abstrct: The Brockett integrtor

More information

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004 Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

Lecture 1. Functional series. Pointwise and uniform convergence.

Lecture 1. Functional series. Pointwise and uniform convergence. 1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is

More information

Solution to Fredholm Fuzzy Integral Equations with Degenerate Kernel

Solution to Fredholm Fuzzy Integral Equations with Degenerate Kernel Int. J. Contemp. Mth. Sciences, Vol. 6, 2011, no. 11, 535-543 Solution to Fredholm Fuzzy Integrl Equtions with Degenerte Kernel M. M. Shmivnd, A. Shhsvrn nd S. M. Tri Fculty of Science, Islmic Azd University

More information

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), ) Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s

More information

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE HANS RINGSTRÖM. Questions nd exmples In the study of Fourier series, severl questions rise nturlly, such s: () (2) re there conditions on c n, n Z, which ensure

More information

QUADRATURE is an old-fashioned word that refers to

QUADRATURE is an old-fashioned word that refers to World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd

More information

Expected Value of Function of Uncertain Variables

Expected Value of Function of Uncertain Variables Journl of Uncertin Systems Vol.4, No.3, pp.8-86, 2 Online t: www.jus.org.uk Expected Vlue of Function of Uncertin Vribles Yuhn Liu, Minghu H College of Mthemtics nd Computer Sciences, Hebei University,

More information

Chapter 1. Basic Concepts

Chapter 1. Basic Concepts Socrtes Dilecticl Process: The Þrst step is the seprtion of subject into its elements. After this, by deþning nd discovering more bout its prts, one better comprehends the entire subject Socrtes (469-399)

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

S. S. Dragomir. 2, we have the inequality. b a

S. S. Dragomir. 2, we have the inequality. b a Bull Koren Mth Soc 005 No pp 3 30 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Abstrct Compnions of Ostrowski s integrl ineulity for bsolutely

More information

The Riemann-Lebesgue Lemma

The Riemann-Lebesgue Lemma Physics 215 Winter 218 The Riemnn-Lebesgue Lemm The Riemnn Lebesgue Lemm is one of the most importnt results of Fourier nlysis nd symptotic nlysis. It hs mny physics pplictions, especilly in studies of

More information

A New Grey-rough Set Model Based on Interval-Valued Grey Sets

A New Grey-rough Set Model Based on Interval-Valued Grey Sets Proceedings of the 009 IEEE Interntionl Conference on Systems Mn nd Cybernetics Sn ntonio TX US - October 009 New Grey-rough Set Model sed on Intervl-Vlued Grey Sets Wu Shunxing Deprtment of utomtion Ximen

More information

REGULARITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2

REGULARITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2 EGULAITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2 OVIDIU SAVIN AND ENICO VALDINOCI Abstrct. We show tht the only nonlocl s-miniml cones in 2 re the trivil ones for ll s 0, 1). As consequence we obtin tht

More information

On the Continuous Dependence of Solutions of Boundary Value Problems for Delay Differential Equations

On the Continuous Dependence of Solutions of Boundary Value Problems for Delay Differential Equations Journl of Computtions & Modelling, vol.3, no.4, 2013, 1-10 ISSN: 1792-7625 (print), 1792-8850 (online) Scienpress Ltd, 2013 On the Continuous Dependence of Solutions of Boundry Vlue Problems for Dely Differentil

More information

Chaos in drive systems

Chaos in drive systems Applied nd Computtionl Mechnics 1 (2007) 121-126 Chos in drive systems Ctird Krtochvíl, Mrtin Houfek, Josef Koláčný b, Romn Kříž b, Lubomír Houfek,*, Jiří Krejs Institute of Thermomechnics, brnch Brno,

More information

FUZZY HOMOTOPY CONTINUATION METHOD FOR SOLVING FUZZY NONLINEAR EQUATIONS

FUZZY HOMOTOPY CONTINUATION METHOD FOR SOLVING FUZZY NONLINEAR EQUATIONS VOL NO 6 AUGUST 6 ISSN 89-668 6-6 Asin Reserch Publishing Networ (ARPN) All rights reserved wwwrpnjournlscom FUZZY HOMOTOPY CONTINUATION METHOD FOR SOLVING FUZZY NONLINEAR EQUATIONS Muhmmd Zini Ahmd Nor

More information

Math 270A: Numerical Linear Algebra

Math 270A: Numerical Linear Algebra Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner

More information

On Error Sum Functions Formed by Convergents of Real Numbers

On Error Sum Functions Formed by Convergents of Real Numbers 3 47 6 3 Journl of Integer Sequences, Vol. 4 (), Article.8.6 On Error Sum Functions Formed by Convergents of Rel Numbers Crsten Elsner nd Mrtin Stein Fchhochschule für die Wirtschft Hnnover Freundllee

More information

Remark on boundary value problems arising in Ginzburg-Landau theory

Remark on boundary value problems arising in Ginzburg-Landau theory Remrk on boundry vlue problems rising in Ginzburg-Lndu theory ANITA KIRICHUKA Dugvpils University Vienibs Street 13, LV-541 Dugvpils LATVIA nit.kiricuk@du.lv FELIX SADYRBAEV University of Ltvi Institute

More information

An optimal 3-point quadrature formula of closed type and error bounds

An optimal 3-point quadrature formula of closed type and error bounds Revist Colombin de Mtemátics Volumen 8), págins 9- An optiml 3-point qudrture formul of closed type nd error bounds Un fórmul de cudrtur óptim de 3 puntos de tipo cerrdo y error de fronter Nend Ujević,

More information

A Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System

A Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System SPIE Aerosense 001 Conference on Signl Processing, Sensor Fusion, nd Trget Recognition X, April 16-0, Orlndo FL. (Minor errors in published version corrected.) A Signl-Level Fusion Model for Imge-Bsed

More information

Predict Global Earth Temperature using Linier Regression

Predict Global Earth Temperature using Linier Regression Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id

More information

New Integral Inequalities for n-time Differentiable Functions with Applications for pdfs

New Integral Inequalities for n-time Differentiable Functions with Applications for pdfs Applied Mthemticl Sciences, Vol. 2, 2008, no. 8, 353-362 New Integrl Inequlities for n-time Differentible Functions with Applictions for pdfs Aristides I. Kechriniotis Technologicl Eductionl Institute

More information

Variational Problems

Variational Problems Vritionl Problems Com S 477/577 Notes Yn-Bin Ji Dec 7, 017 ThevritionlderivtiveoffunctionlJ[y]cnbedefinedsδJ/δy = F y x,y,y d dx F y x,y,y [1, pp. 7 9]. Euler s eqution essentilly sttes tht the vritionl

More information

Chapters 4 & 5 Integrals & Applications

Chapters 4 & 5 Integrals & Applications Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions

More information

Definite integral. Mathematics FRDIS MENDELU

Definite integral. Mathematics FRDIS MENDELU Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the

More information

Overview of Calculus I

Overview of Calculus I Overview of Clculus I Prof. Jim Swift Northern Arizon University There re three key concepts in clculus: The limit, the derivtive, nd the integrl. You need to understnd the definitions of these three things,

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

Line Integrals. Partitioning the Curve. Estimating the Mass

Line Integrals. Partitioning the Curve. Estimating the Mass Line Integrls Suppose we hve curve in the xy plne nd ssocite density δ(p ) = δ(x, y) t ech point on the curve. urves, of course, do not hve density or mss, but it my sometimes be convenient or useful to

More information

An instructive toy model: two paradoxes

An instructive toy model: two paradoxes Tel Aviv University, 2006 Gussin rndom vectors 27 3 Level crossings... the fmous ice formul, undoubtedly one of the most importnt results in the ppliction of smooth stochstic processes..j. Adler nd J.E.

More information