Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.
|
|
- Regina Black
- 5 years ago
- Views:
Transcription
1 Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need to establish that the largest invariant set in E is the origin, that is: = { }. his is done by showing that no solution an stay identially in E other than the trivial solution x() t. heorem 4.1 (due to Barbashin and Krasovskii) Let x = be an equilibrium point for (3.1). Let V : D R be a ontinuously n differentiable positive definite funtion on a domain D R ontaining the origin, suh that V x() t S = x D: V x = and suppose that no other solution in D. Let { } an stay in S, other than the trivial solution x( t). hen the origin is loally asymptotially stable. If, in addition, V ( x ) is radially unbounded then the origin is globally asymptotially stable. Note that if V ( x) is negative definite then { } S = and the above theorem oinides with the Lyapunov nd theorem. Also note that the LaSalle s invariant set theorems are appliable to autonomous system only. Example 4. Consider the 1 st order system together with its adaptive ontrol law he dynamis of the adaptive gain ˆk ( t ) is x = ax+ u ˆ u = k t x ˆk = γ x where γ > is alled the rate of adaptation. hen the losed-loop system beomes: ( ˆ x = k() t a) x ˆ k = γ x he line x = represents the system equilibrium set. We want to show that the trajetories approah this equilibrium set, as t, whih means that the adaptive ontroller regulates x() t to zero in the presene of onstant unertainty in a. Consider the Lyapunov funtion andidate: 13
2 (, ˆ) = ( ˆ ) V x k x k b γ where b> a. ime derivative of V along the trajetories of the system is given by 1 V ( ˆ) ( ˆ x, k = xx + k b) kˆ= x kˆ a + kˆ b x = x ( b a) γ V x, k ˆ is positive definite and radially unbounded funtion, whose derivative Sine V (, ˆ x k) is semi-negative, the set (, ˆ) : (, ˆ) positively invariant set. hus taking satisfied. he set E is given by ( ˆ) { } Ω = x k R V x k is ompat, Ω =Ω, all the onditions of LaSalle s heorem are {, : } E = x k Ω x=. Beause any point on the line x = is an equilibrium point, E is an invariant set. herefore, in this example = E. From LaSalle s heorem we onlude that every trajetory starting in Ω approahes E, as t, that is () xt as t. oreover, sine (, ˆ) (, ˆ ) Ω an be hosen large enough that ( ˆ) V x k is radially unbounded, the onlusion is global, that is it holds for all initial onditions onstant in the definition of x k beause the x, k Ω. 5. Boundedness and Ultimate Boundedness Consider the nonautonomous system where :, [ ) [, ) D, and f D R (, ) x = f t x (5.1) n is pieewise ontinuous in t, loally Lipshitz in x on D R n is a domain that ontains the origin x =. Note that if the origin is an equilibrium point for (5.1) then by definition: f ( t,) =, t. On the other hand, even if there is no equilibrium at the origin, Lyapunov analysis an still be used to show boundedness of the system trajetories. We begin with a motivating example. Example 5.1 Consider the IVP with nonautonomous salar dynamis x = x+ sin t (5.) xt ( ) = a> > he system has no equilibrium points. he IVP expliit solution an be easily found and shown to be bounded for all t t, uniformly in t, that is with a bound b independent of t. In this ase, the solution is said to be uniformly ultimately bounded (UUB), and b is alled the ultimate bound, (prove it). 14
3 he UUB property of (5.) an be established via Lyapunov analysis and without using x the expliit solution of the state equation. In fat, starting with V ( x ) =, we an alulate the time derivative of V, along the system trajetories. V ( x) = xx = x( x+ sin t) = x + xsin t x + x = x ( x ) It immediately follows that V ( x) <, x > In other words, the time derivative of V is negative outside of the ompat set B = { x }, or equivalently, all the solutions that start outside of B will enter this set within a finite time, and will remain inside the set for all future times. Formally, it an be stated as follows. Choose >. hen all the solutions that start in the set B = { V( x) } B x will remain therein, for all future time. he latter takes plae sine V is negative on the boundary. Hene, the solutions are uniformly bounded. oreover, an ultimate bound of these solutions an also be found. Choose ε suh that < ε < hen V is negative in the annulus set { V ( x) } ( ()) { V( x) ε} ε, whih implies that in this set V x t will derease monotonially in time, until the solution enters the set. From that time on, the trajetory will remain in the set, beause V is x negative on the boundary: V ( x) = ε. Sine V ( x ) =, we an onlude that this solution is UUB with the ultimate bound x ε. Definition 5.1 he solutions of (5.1) are uniformly bounded if there exists a positive onstant, independent of t, and for every a (, ), there is β β ( a) = >, independent of t, suh that x t a x t β, t t (5.3) globally uniformly bounded if (5.3) holds for arbitrarily large a 15
4 uniformly ultimately bounded with ultimate bound b if there exist positive onstants b and, independent of a,, there is (, ) = a b, independent of t, suh that t, and for every x t a x t b, t t + (5.4) globally uniformly ultimately bounded if (5.4) holds for arbitrarily large a. Figure 5.1: UUB Conept In the definition above, the term uniform indiates that the bound b does not depend on t. he term ultimate indiates that boundedness holds after the lapse of a ertain time. he onstant defines a neighborhood of the origin, independent of t, suh that all trajetories starting in the neighborhood will remain bounded in time. If an be hosen arbitrarily large then the UUB notion beomes global. Basially, the notion of UUB an be onsidered as a milder form of stability in the sense of Lyapunov (SISL). A omparison between SISL and UUB onepts is given below. SISL is defined with respet to an equilibrium, while UUB is not. Asymptoti SISL is a strong property that is very diffiult to ahieve in pratial dynamial systems. SISL requires the ability to keep the state arbitrarily lose to the system equilibrium by starting suffiiently lose to it. his is still too strong a requirement for pratial systems operating in the presene of unknown disturbanes. 16
5 he main differene between UUB and SISL is that the UUB bound b annot be made arbitrarily small by starting loser to the equilibrium or the origin. In pratial systems, the bound b depends on disturbanes and system unertainties. o demonstrate how Lyapunov analysis an be used to study UUB, onsider a ontinuously differentiable positive definite funtion V ( x ). Choose < ε <. Suppose that the sets Ω ε = { V ( x) ε} and Ω = { V ( x) } are ompat. Let Λ= ε V ( x) =Ω Ω ε { } and suppose that it is known that the time derivative of ( ) V x t along the trajetories of the nonautonomous dynamial system (5.1) is negative definite inside Λ, that is V ( x() t ) W( x( t) ) <, x Λ, t t where W ( x() t ) is a ontinuous positive definite funtion. Sine V is negative in Λ, a trajetory starting in Λ must move in the diretion of dereasing V ( x() t ). In fat, it an be shown that in the set Λ the trajetory behaves as if the origin was uniformly asymptotially stable, (whih it does not have to be in this ase). Consequently, the funtion V ( x() t ) will ontinue dereasing until the trajetory enters the set Ω ε in finite time and stays there for all future time. Hene, the solutions of (5.1) are UUB with the ultimate bound b= max x. A sketh of the sets Λ, Ω, Ω ε is shown in Figure 5.. x Ω ε Figure 5.: UUB by Lyapunov Analysis In many problems, the relation V ( t, x) W( x) is derived and shown to be valid on a domain whih is speified in terms of x. In suh ases, the UUB analysis involves finding the orresponding domains of attration and an ultimate bound. 6. Barbalat s Lemma and Invariane like heorems 17
6 For autonomous systems, LaSalle s invariane set theorems allow asymptoti stability onlusions to be drawn even when V is only negative semi-definite in a domain Ω. In that ase, the system trajetory approahes the largest invariant set E, whih is a subset of all points x Ω where V ( x ) =. However the invariant set theorems are not appliable to nonautonomous systems. In the ase of the latter, it may not even be lear how to define a set E, sine V may expliitly depend on both t and x. Even when V = V ( x) does not expliitly depend on t the nonautonomous nature of the system dynamis preludes the use of the LaSalle s invariant set theorems. Example 6.1 he losed-loop error dynamis of an adaptive ontrol system for 1 st order plant with one unknown parameter is e = e+ θ w( t) θ = ew t () where e represents the traking error and wt is a bounded funtion of time t. Due to the presene of wt (), the system dynamis is nonautonomous. Consider the Lyapunov funtion andidate V ( e, θ ) = e + θ Its time derivative along the system trajetories is V e, θ = ee + θθ = e e+ θ w t + θ ew t = e ( ) ( ) his inequality implies that V is a dereasing funtion of time, and therefore, both et and θ () t are bounded signals of time. But due to the nonautonomous nature of the system dynamis, the LaSalle s invariane set theorems annot be used to onlude the onvergene of et () to the origin., then we may expet that the trajetory of the system In general, if V ( t x) W( x) approahes the set { W( x ) = }, as t. Before we formulate main results, we state a lemma that is interesting in its own sake. he lemma is an important result about asymptoti properties of funtions and their derivatives and it is known as the Barbalat s lemma. Definition 6.1 (uniform ontinuity) A funtion f (): t R R is said to be uniformly ontinuous if ε > = ( ε) > t t1 f ( t) f ( t1) ε (6.1) Note that t 1 and t play a symmetri role in the definition of the uniform ontinuity. Lemma 6.1 (due to Barbalat) 18
7 Let f (): t R R be differentiable and has a finite limit as t. If f () t is uniformly ontinuous then f () t as t. Lemma 6. If f () t is bounded then f () t is uniformly ontinuous. An immediate and pratial orollary of Barbalat s lemma an now be stated. Corollary 6.1 If f (): t R R is twie differentiable, has a finite limit, and its nd derivative is bounded then f () t as t. In general, the fat that derivative tends to zero does not imply that the funtion has a limit. Also, the onverse is not true. f t C f t Example 6. As t, f () t sin ( ln t) os ln t = does not have a limit, while f () t =. t t t f t e sin e f t = e t sin e t + e t sin e t. As t, () =, while Example 6.3 Consider an LI system x = Ax+ Bu with a Hurwitz matrix A and a uniformly bounded in time input u( t ). hese two fats imply that the state x() t is bounded. hus, the state derivative x ( t) is bounded. Let y = Cx represent the system output. hen y = Cx is bounded and, onsequently the system output y () t is a uniformly ontinuous funtion of time. oreover, if the input u() t = u is onstant then the output y ( t ) tends to a limit, as t. he latter ombined with the fat that y is uniformly ontinuous implies, (by Barbalat s Lemma), that the output time derivative y asymptotially approahes zero. o apply Barbalat s lemma to the analysis of nonautonomous dynami systems we state the following immediate orollary. Corollary 6. (Lyapunov-like Lemma) V = V t, x is suh that If a salar funtion V ( t, x ) is lower bounded V ( t, x) x = f ( t, x ) 19
8 V ( t, x) then V ( t x) is uniformly ontinuous in time,, as t. Notie that the first two assumptions imply that V ( t, x ) tends to a limit. he latter oupled with the 3 rd assumption proves (using Barbalat s lemma) the orollary. Example 6.4 Consider again the losed-loop error dynamis of an adaptive ontrol system from Example 6.1. Choosing V ( e, θ ) = e + θ, it was shown that along the system trajetories: V ( e, θ ) = e. he nd time derivative of V is V ( e, θ) = 4ee = 4e( e+ θ w( t) ) Sine wt () is bounded by hypothesis, and etand θ ( t) were shown to be bounded, it is lear that V is bounded.. Hene, V is uniformly ontinuous and by the Barbalat s lemma (or the Lyapunov-like lemma), V whih in turn indiates that the traking error et tends to zero, as t. Often, the Barbalat Lemma is formulated with respet to an integral of a funtion, rather than to the funtion itself. Lemma 6.3 (due to Barbalat) If f : R + R is uniformly ontinuous for all t, and if the limit of the integral exists and finite, then t ( τ ) lim f dτ (6.) t lim f t = (6.3) t Definition 6. he L norm of a funtion f : R + R is defined as: p for all p = [ 1, ), while denotes the L norm of f. 1 p p τ (6.4) f () t = f ( τ ) d p max f ( t) f t = (6.5) t
9 We say that f Lp if () p f t exists and finite. In partiular, f L dτ <. if f ( τ) he next statement diretly follows from Lemma 6.3, if one hooses f () t g() t p Corollary 6. If gg, L and g() t Lp, for some p = [ 1, ), then g( t) lim =. t =. 1
10 Part II: System ID ethods 1. Introdution Read [Ioannou], Chapters 1. he goal of any parameter estimation method is to infer the values of unertain system parameters from the measurements of input and output signals of the system. Parameter estimation an be done either on-line or off-line. he latter may be preferable if the parameters are onstant and there is a suffiient time for estimation to be performed. However, for parameters that vary during operation, on-line parameter estimation is neessary to keep trak of the parameter values. Sine problems in robust and adaptive ontrol ontext usually involve time-varying parameters, on-line estimation methods are thus more relevant. In what follows, several basi methods of on-line estimation will be presented. Note that, although the main purpose of the on-line estimators may be to provide parameter estimates for adaptive / self-tuning ontrol, they an also be used for other purposes, suh as system health monitoring and failure detetion.. Stati Linear in Parameters odel An estimation model that relates the available data to the unknown parameters forms the basis of any parameter estimation method. A quite general model for parameter estimation appliations is in the linear parameterization form: y t =Θ Φ x t + ε t (.1) where () m n y t R, x() t R are the system output and input signals, the ( N m) N m Θ R ontains the unknown parameters to be estimated, x( t) R N matrix Φ (alled the regressor ) represents the N dimensional vetor of hosen basis funtions, and ε t R denotes the non-parametri unertainties in the system (suh as noise, m () modeling errors, et). Note that in (.1) both y ( t ) and x( t) Φ are known quantities. his implies that (.1) is simply is a system of linear equations in terms of the unknown Θ. For every time instant t there is suh a system. So if ontinuous measurements of y () t and Φ ( x() t ) are available on-line then, in general, there will be more equations than the unknowns. However, this will failitate parameter estimation in the presene of inevitable noise and modeling error. Φ for a period of time, and solves the equations one and for all. In on-line estimation, equations are In off-line estimation, one ollets the data points of y ( t ) and x( t)
11 solved reursively, implying that the estimated value of ˆ ˆ ( t) set of data y () t and Φ x() t beomes available. Θ=Θ is updated one a new Stati linear in parameters model (.1), although simple, is atually quite general, as seen from the following example. Example.1: Linear parameterization of 1 st order dynamis Consider the 1 st order system x = x+ ax+ bu+ξ ( t) (.) x( ) = where x, u are the salar state and input respetively, a and b are the unknown onstants to be identified on-line using the available measurements of x and u. Also in (.), ξ () t represent bounded unknown disturbane, (suh as proess noise). he very first step in the design of an on-line parameter identifier is to lump the unknown parameters a and b in a vetor, that is one needs to rewrite (.) in the form of (.1). oward that end, one an easily integrate (.) and get: () ( ) t t t t τ t τ t τ x t = x + a e x τ dτ + b e u τ dτ + e ξ τ dτ () () ε () xf t uf t t (.3) or, equivalently: x() t = axf ( t) + buf ( t) + ε ( t) (.4) Finally, the system an be written in the desired form of (.1) : Θ ( t) () t xf x = ( a b) + ε () t u (.5) f Example.: Linear parameterization of airraft pith dynamis Consider the following pith dynamis of a rigid-body, fixed-wing airraft. IY q = ( α) + q q+ + pr( Iz Ix) + ( r p ) Ixz I is the pith inertia, ( q, q) Φ (, ) aero IC p r (.6) where Y are the pith aeleration and angular rate. Also in (.6), aero denotes the airraft total aerodynami pithing moment whih is omprised of the baseline omponent ( α ), pith damping q q, and pith inrement due to elevator defletion. Finally, (, ) IC pr represents the sum of ross-oupling terms 3
12 that are due to the angular roll and yaw rates, ( p, r ), saled by the airraft inertia omponents ( x, z, xz) I I I. Suppose that the aerodynami pithing moment omponents α, q, q, are unknown and need to be estimated on-line using available input-output measurements ( α, p, q, r, ) and q, respetively. It is further assumed that the airraft inertia terms are onstant and known. Consequently, the inertial oupling term IC ( pr, ) in the right hand side of (.6) is known, as well. Rewriting (.6) using aerodynami oeffiients yields: I Y IC q q C C C q QS = α + + IY V (.7) In (.7), Q denotes the dynamis pressure, S is the airraft wing referene area, and is the mean geometri hord. For simpliity, the unknown pith ontrol effetiveness C and the pith rate damping oeffiient C are assumed to be onstant. On the other q hand, the unertain aerodynami pithing moment omponent C ( α ) depends on the airraft angle of attak α. Suppose that the pithing moment an be written as: C ( α ) = θ Φ α + ε α (.8) N where θ R is the N dimensional vetor of unknown onstant parameters, and N ( α ) R basis funtions, suh as splines, polynomials, RBF-s, et. Also in (.8), ε Φ is the N dimensional regressor vetor that onsists of appropriately hosen α is the non-parametri approximation error that an be made small by hoosing a suffiiently large number of basis funtions N. Substituting (.8) into (.7) gives the desired linear in parameters model. I Y IC q q C C q QS = θ Φ α ε α IY V (.9) or equivalently: Φ ( α ) I Y IC q q ( θ C C = ) q + ε (.1) QS IY V Θ yt () Ψ ( α, qv,, ) Note that the error term ε may ontain other non-parametri unertainties in the system, suh as measurement noise, numerial inauraies et. As a result, the airraft pith dynamis relation (.6) is represented using the linear-in-parameters model: y t =Θ Ψ t + ε t (.11) 4
13 Note that the stati model (.11) requires knowledge of the pith aeleration signal q. ost often in real-life appliations, vehile rotational aelerations are not available as on board measurements. In this ase, an alternative modeling approah an be utilized. We rewrite (.6) in the form, aero ( α, q, ) + IC ( p, r) q = af q+ af q+ (.1) IY where a f >, and integrate both sides, (assuming zero initial onditions). t t t aero ( α( τ), q( τ), ( τ) ) IC ( p( τ), r t t t ( τ τ τ τ )) q() t = a e q( τ ) dτ + e dτ + e dτ (.13) hen, f I Y I Y qf () t aero f () t IC f () t () IC f f f aero f yt () q t t a q t = t (.14) where y () t is the known output, while aero f ( t ) represents the vehile filtered pith aerodynamis. his term an be modeled as a linear in unknown parameters funtion. Φ ( α ) q QS aero f = ( θ C C ) q + ε f (.15) V IY Θ Ψ ( α, qv,, ) hus, again we arrive at a stati model in the form of (.11). Remark 1: In some pratial appliations, using a filtered version of the derivative, suh f s as q = q, in plae of the derivative q works well beause the signal q that is τ q s + 1 being filtered is ontinuous. herefore, the signal energy is at low frequenies, where the speified filter is a good approximation to the derivative, ( 1 defines the bandwidth of τ q the derivative filter). It is evident that the derivative approximation error an be folded into the bounded non-parametri unertainty vetor ε. 5
Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties
Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationCHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS. Professor Dae Ryook Yang
CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 208 Dept. of Chemial and Biologial Engineering 0- Road Map of the Leture X Stability of losed-loop ontrol system
More informationA Characterization of Wavelet Convergence in Sobolev Spaces
A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat
More informationEvaluation of effect of blade internal modes on sensitivity of Advanced LIGO
Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationA NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS
Vietnam Journal of Mehanis, VAST, Vol. 4, No. (), pp. A NONLILEAR CONTROLLER FOR SHIP AUTOPILOTS Le Thanh Tung Hanoi University of Siene and Tehnology, Vietnam Abstrat. Conventional ship autopilots are
More informationMost results in this section are stated without proof.
Leture 8 Level 4 v2 he Expliit formula. Most results in this setion are stated without proof. Reall that we have shown that ζ (s has only one pole, a simple one at s =. It has trivial zeros at the negative
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationmax min z i i=1 x j k s.t. j=1 x j j:i T j
AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be
More informationSensitivity Analysis in Markov Networks
Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores
More information). In accordance with the Lorentz transformations for the space-time coordinates of the same event, the space coordinates become
Relativity and quantum mehanis: Jorgensen 1 revisited 1. Introdution Bernhard Rothenstein, Politehnia University of Timisoara, Physis Department, Timisoara, Romania. brothenstein@gmail.om Abstrat. We first
More informationSQUARE ROOTS AND AND DIRECTIONS
SQUARE ROOS AND AND DIRECIONS We onstrut a lattie-like point set in the Eulidean plane that eluidates the relationship between the loal statistis of the frational parts of n and diretions in a shifted
More informationBäcklund Transformations: Some Old and New Perspectives
Bäklund Transformations: Some Old and New Perspetives C. J. Papahristou *, A. N. Magoulas ** * Department of Physial Sienes, Helleni Naval Aademy, Piraeus 18539, Greee E-mail: papahristou@snd.edu.gr **
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationLECTURE NOTES FOR , FALL 2004
LECTURE NOTES FOR 18.155, FALL 2004 83 12. Cone support and wavefront set In disussing the singular support of a tempered distibution above, notie that singsupp(u) = only implies that u C (R n ), not as
More informationAircraft CAS Design with Input Saturation Using Dynamic Model Inversion
International Journal of Control, Automation, and Systems Vol., No. 3, September 003 35 Airraft CAS Design with Input Saturation sing Dynami Model Inversion Sangsoo Lim and Byoung Soo Kim Abstrat: This
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationImproved Extended Kalman Filter for Parameters Identification
Improved Extended Kalman Filter for Parameters Identifiation Peipei Zhao. Zhenyu Wang. ao Lu. Yulin Lu Institute of Disaster Prevention, Langfang,Hebei Provine, China SUMMARY: he extended Kalman filter
More informationWave Propagation through Random Media
Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationEECS 120 Signals & Systems University of California, Berkeley: Fall 2005 Gastpar November 16, Solutions to Exam 2
EECS 0 Signals & Systems University of California, Berkeley: Fall 005 Gastpar November 6, 005 Solutions to Exam Last name First name SID You have hour and 45 minutes to omplete this exam. he exam is losed-book
More informationChapter Review of of Random Processes
Chapter.. Review of of Random Proesses Random Variables and Error Funtions Conepts of Random Proesses 3 Wide-sense Stationary Proesses and Transmission over LTI 4 White Gaussian Noise Proesses @G.Gong
More informationEstimating the probability law of the codelength as a function of the approximation error in image compression
Estimating the probability law of the odelength as a funtion of the approximation error in image ompression François Malgouyres Marh 7, 2007 Abstrat After some reolletions on ompression of images using
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More informationRelativity in Classical Physics
Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of
More informationON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS
ON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS MARIO LEFEBVRE and JEAN-LUC GUILBAULT A ontinuous-time and ontinuous-state stohasti proess, denoted by {Xt), t }, is defined from a proess known as
More informationRelativistic Dynamics
Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationModeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers
A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation
More informationECE-320 Linear Control Systems. Winter 2013, Exam 1. No calculators or computers allowed, you may leave your answers as fractions.
ECE-320 Linear Control Systems Winter 2013, Exam 1 No alulators or omputers allowed, you may leave your answers as frations. All problems are worth 3 points unless noted otherwise. Total /100 1 Problems
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationChapter 2 Lecture 8 Longitudinal stick fixed static stability and control 5 Topics
Flight dynamis II Stability and ontrol hapter 2 Leture 8 Longitudinal stik fied stati stability and ontrol 5 Topis 2.6 ontributions of power plant to mg and mα 2.6.1 Diret ontributions of powerplant to
More informationSimple FIR Digital Filters. Simple FIR Digital Filters. Simple Digital Filters. Simple FIR Digital Filters. Simple FIR Digital Filters
Simple Digital Filters Later in the ourse we shall review various methods of designing frequeny-seletive filters satisfying presribed speifiations We now desribe several low-order FIR and IIR digital filters
More informationModes are solutions, of Maxwell s equation applied to a specific device.
Mirowave Integrated Ciruits Prof. Jayanta Mukherjee Department of Eletrial Engineering Indian Institute of Tehnology, Bombay Mod 01, Le 06 Mirowave omponents Welome to another module of this NPTEL mok
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationSupporting Information
Supporting Information Olsman and Goentoro 10.1073/pnas.1601791113 SI Materials Analysis of the Sensitivity and Error Funtions. We now define the sensitivity funtion Sð, «0 Þ, whih summarizes the steepness
More informationQuasi-Monte Carlo Algorithms for unbounded, weighted integration problems
Quasi-Monte Carlo Algorithms for unbounded, weighted integration problems Jürgen Hartinger Reinhold F. Kainhofer Robert F. Tihy Department of Mathematis, Graz University of Tehnology, Steyrergasse 30,
More informationSensor management for PRF selection in the track-before-detect context
Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationLongitudinal Static Stability
ongitudinal Stati Stability Some definitions C m M V S pithing moment without dimensions (so without influene of ρ, V and S) it is a shape parameter whih varies with the angle of attak. Note the hord in
More informationMeasuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach
Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La
More informationThe Laws of Acceleration
The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the
More informationConformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices
AAS 03-190 Conformal Mapping among Orthogonal, Symmetri, and Skew-Symmetri Matries Daniele Mortari Department of Aerospae Engineering, Texas A&M University, College Station, TX 77843-3141 Abstrat This
More informationMethods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationChapter 8 Hypothesis Testing
Leture 5 for BST 63: Statistial Theory II Kui Zhang, Spring Chapter 8 Hypothesis Testing Setion 8 Introdution Definition 8 A hypothesis is a statement about a population parameter Definition 8 The two
More informationSome facts you should know that would be convenient when evaluating a limit:
Some fats you should know that would be onvenient when evaluating a it: When evaluating a it of fration of two funtions, f(x) x a g(x) If f and g are both ontinuous inside an open interval that ontains
More informationUniversity of Groningen
University of Groningen Port Hamiltonian Formulation of Infinite Dimensional Systems II. Boundary Control by Interonnetion Mahelli, Alessandro; van der Shaft, Abraham; Melhiorri, Claudio Published in:
More informationCALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS
International Journal of Modern Physis A Vol. 24, No. 5 (2009) 974 986 World Sientifi Publishing Company CALCULATION OF NONLINEAR TUNE SHIFT USING BEAM POSITION MEASUREMENT RESULTS PAVEL SNOPOK, MARTIN
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E')
22.54 Neutron Interations and Appliations (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') Referenes -- J. R. Lamarsh, Introdution to Nulear Reator Theory (Addison-Wesley, Reading, 1966),
More informationGreen s function for the wave equation
Green s funtion for the wave equation Non-relativisti ase January 2019 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 43 and 44): 1 2 A 2 2 2 A = µ 0
More informationSpeed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion
Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3
More informationarxiv: v2 [math.pr] 9 Dec 2016
Omnithermal Perfet Simulation for Multi-server Queues Stephen B. Connor 3th Deember 206 arxiv:60.0602v2 [math.pr] 9 De 206 Abstrat A number of perfet simulation algorithms for multi-server First Come First
More informationSimplified Buckling Analysis of Skeletal Structures
Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation
More informationModal Horn Logics Have Interpolation
Modal Horn Logis Have Interpolation Marus Kraht Department of Linguistis, UCLA PO Box 951543 405 Hilgard Avenue Los Angeles, CA 90095-1543 USA kraht@humnet.ula.de Abstrat We shall show that the polymodal
More informationThe Lorenz Transform
The Lorenz Transform Flameno Chuk Keyser Part I The Einstein/Bergmann deriation of the Lorentz Transform I follow the deriation of the Lorentz Transform, following Peter S Bergmann in Introdution to the
More informationCritical Reflections on the Hafele and Keating Experiment
Critial Refletions on the Hafele and Keating Experiment W.Nawrot In 1971 Hafele and Keating performed their famous experiment whih onfirmed the time dilation predited by SRT by use of marosopi loks. As
More informationthe following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined
How should a snake turn on ie: A ase study of the asymptoti isoholonomi problem Jianghai Hu, Slobodan N. Simić, and Shankar Sastry Department of Eletrial Engineering and Computer Sienes University of California
More informationThe transition between quasi-static and fully dynamic for interfaces
Physia D 198 (24) 136 147 The transition between quasi-stati and fully dynami for interfaes G. Caginalp, H. Merdan Department of Mathematis, University of Pittsburgh, Pittsburgh, PA 1526, USA Reeived 6
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationHILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES
HILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES L ERBE, A PETERSON AND S H SAKER Abstrat In this paper, we onsider the pair of seond-order dynami equations rt)x ) ) + pt)x
More informationAn Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems
An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,
More information3 Tidal systems modelling: ASMITA model
3 Tidal systems modelling: ASMITA model 3.1 Introdution For many pratial appliations, simulation and predition of oastal behaviour (morphologial development of shorefae, beahes and dunes) at a ertain level
More informationOn Component Order Edge Reliability and the Existence of Uniformly Most Reliable Unicycles
Daniel Gross, Lakshmi Iswara, L. William Kazmierzak, Kristi Luttrell, John T. Saoman, Charles Suffel On Component Order Edge Reliability and the Existene of Uniformly Most Reliable Uniyles DANIEL GROSS
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationDIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS
CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More informationWord of Mass: The Relationship between Mass Media and Word-of-Mouth
Word of Mass: The Relationship between Mass Media and Word-of-Mouth Roman Chuhay Preliminary version Marh 6, 015 Abstrat This paper studies the optimal priing and advertising strategies of a firm in the
More informationANTI-SATURATION DYNAMIC SURFACE CONTROL FOR SPACECRAFT TERMINAL SAFE APPROACH BASED ON COMMAND FILTER. Guanqun Wu, Shenmin Song and Jingguang Sun
International Journal of Innovative Computing, Information and Control ICIC International 2018 ISSN 1349-4198 Volume 14, Number 1, February 2018 pp. 33 52 ANTI-SATURATION DYNAMIC SURFACE CONTROL FOR SPACECRAFT
More informationSPLINE ESTIMATION OF SINGLE-INDEX MODELS
SPLINE ESIMAION OF SINGLE-INDEX MODELS Li Wang and Lijian Yang University of Georgia and Mihigan State University Supplementary Material his note ontains proofs for the main results he following two propositions
More informationSensitivity analysis for linear optimization problem with fuzzy data in the objective function
Sensitivity analysis for linear optimization problem with fuzzy data in the objetive funtion Stephan Dempe, Tatiana Starostina May 5, 2004 Abstrat Linear programming problems with fuzzy oeffiients in the
More informationPlanning with Uncertainty in Position: an Optimal Planner
Planning with Unertainty in Position: an Optimal Planner Juan Pablo Gonzalez Anthony (Tony) Stentz CMU-RI -TR-04-63 The Robotis Institute Carnegie Mellon University Pittsburgh, Pennsylvania 15213 Otober
More informationThe gravitational phenomena without the curved spacetime
The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,
More informationDirectional Coupler. 4-port Network
Diretional Coupler 4-port Network 3 4 A diretional oupler is a 4-port network exhibiting: All ports mathed on the referene load (i.e. S =S =S 33 =S 44 =0) Two pair of ports unoupled (i.e. the orresponding
More informationA Functional Representation of Fuzzy Preferences
Theoretial Eonomis Letters, 017, 7, 13- http://wwwsirporg/journal/tel ISSN Online: 16-086 ISSN Print: 16-078 A Funtional Representation of Fuzzy Preferenes Susheng Wang Department of Eonomis, Hong Kong
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More informationFinite-time stabilization of chaotic gyros based on a homogeneous supertwisting-like algorithm
OP Conferene Series: Materials Siene Engineering PAPER OPEN ACCESS Finite-time stabilization of haoti gyros based on a homogeneous supertwisting-like algorithm To ite this artile: Pitha Khamsuwan et al
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationBrazilian Journal of Physics, vol. 29, no. 3, September, Classical and Quantum Mechanics of a Charged Particle
Brazilian Journal of Physis, vol. 9, no. 3, September, 1999 51 Classial and Quantum Mehanis of a Charged Partile in Osillating Eletri and Magneti Fields V.L.B. de Jesus, A.P. Guimar~aes, and I.S. Oliveira
More informationCase I: 2 users In case of 2 users, the probability of error for user 1 was earlier derived to be 2 A1
MUTLIUSER DETECTION (Letures 9 and 0) 6:33:546 Wireless Communiations Tehnologies Instrutor: Dr. Narayan Mandayam Summary By Shweta Shrivastava (shwetash@winlab.rutgers.edu) bstrat This artile ontinues
More informationCMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017
CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,
More informationMOLECULAR ORBITAL THEORY- PART I
5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.
More informationBerry s phase for coherent states of Landau levels
Berry s phase for oherent states of Landau levels Wen-Long Yang 1 and Jing-Ling Chen 1, 1 Theoretial Physis Division, Chern Institute of Mathematis, Nankai University, Tianjin 300071, P.R.China Adiabati
More informationDifferential Equations 8/24/2010
Differential Equations A Differential i Equation (DE) is an equation ontaining one or more derivatives of an unknown dependant d variable with respet to (wrt) one or more independent variables. Solution
More information9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry
9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived
More informationDiscrete Bessel functions and partial difference equations
Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat
More informationSURFACE WAVES OF NON-RAYLEIGH TYPE
SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh
More informationTime Domain Method of Moments
Time Domain Method of Moments Massahusetts Institute of Tehnology 6.635 leture notes 1 Introdution The Method of Moments (MoM) introdued in the previous leture is widely used for solving integral equations
More informationON-LINE ESTIMATION OF THE VENTILATION RATE OF GREENHOUSES
ON-LINE ESTIMATION OF THE VENTILATION RATE OF GREENHOUSES J. Bontsema *, E.J. van Henten *, J.G. Kornet *, J. Budding ** and T. Rieswijk ** * Agrotehnology and Food Innovations B.V., Greenhouse Tehnology
More informationSimplification of Network Dynamics in Large Systems
Simplifiation of Network Dynamis in Large Systems Xiaojun Lin and Ness B. Shroff Shool of Eletrial and Computer Engineering Purdue University, West Lafayette, IN 47906, U.S.A. Email: {linx, shroff}@en.purdue.edu
More informationThe Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge
The Conept of Mass as Interfering Photons, and the Originating Mehanism of Gravitation D.T. Froedge V04 Formerly Auburn University Phys-dtfroedge@glasgow-ky.om Abstrat For most purposes in physis the onept
More informationReview of Force, Stress, and Strain Tensors
Review of Fore, Stress, and Strain Tensors.1 The Fore Vetor Fores an be grouped into two broad ategories: surfae fores and body fores. Surfae fores are those that at over a surfae (as the name implies),
More informationADAPTIVE POISSON DISORDER PROBLEM
ADAPTIVE POISSON DISORDER PROBLEM ERHAN BAYRAKTAR, SAVAS DAYANIK, AND IOANNIS KARATZAS Abstrat. We study the quikest detetion problem of a sudden hange in the arrival rate of a Poisson proess from a known
More informationNormative and descriptive approaches to multiattribute decision making
De. 009, Volume 8, No. (Serial No.78) China-USA Business Review, ISSN 57-54, USA Normative and desriptive approahes to multiattribute deision making Milan Terek (Department of Statistis, University of
More informationTHE TWIN PARADOX A RELATIVISTIC DOMAIN RESOLUTION
THE TWIN PARADOX A RELATIVISTIC DOMAIN RESOLUTION Peter G.Bass P.G.Bass www.relativitydomains.om January 0 ABSTRACT This short paper shows that the so alled "Twin Paradox" of Speial Relativity, is in fat
More informationController Design Based on Transient Response Criteria. Chapter 12 1
Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors
More information13.Prandtl-Meyer Expansion Flow
3.Prandtl-eyer Expansion Flow This hapter will treat flow over a expansive orner, i.e., one that turns the flow outward. But before we onsider expansion flow, we will return to onsider the details of the
More information