Moving Horizon Estimation for Staged QP Problems
|
|
- Mildred Berry
- 6 years ago
- Views:
Transcription
1 51s IEEE Conference on Decision and Conrol, Maui, HI, December 212 Moving Horizon Esimaion for Saged QP Problems Eric Chu, Arezou Keshavarz, Dimiry Gorinevsky, and Sephen Boyd Absrac his paper considers moving horizon esimaion (MHE) approach o soluion of saged quadraic programming (QP) problems. Using an insigh ino he consrained soluion srucure for he growing horizon, we develop a very accurae ieraive updae of he arrival cos in he MHE soluion. he updae uses a quadraic approximaion of he arrival cos and informaion abou he previously acive or inacive consrains. In he absence of consrains, he updae is he familiar Kalman filer in informaion form. In he presence of he consrains, he updae requires solving a sequence of linear sysems wih varying size. he proposed MHE updae provides very good performance in numerical examples. his includes problems wih l 1 regularizaion where opimal esimaion allows us o perform online segmenaion of sreaming daa. A. Problem saemen I. INRODUCION Given a series of quadraic forms g (u), g (u,v), where u,v R n and is an ineger index, consider a series of ime saged quadraic programming (QP) opimizaion problems for increasing g (z )+ =1 g (z 1,z ) F eq z = G eq z 1 +h eq F ineq z G ineq z 1 +h ineq = 1,...,, where z R n are decision variables; he marices F eq, F ineq have sizes compaible wih z and he vecors h eq, h ineq ; and denoes componen-wise for wo vecors. he opimizaion inerval [1, ] is furher called full horizon. Our goal is o provide a sequence of opimal soluions Z () = (z,...,z ) as increases. he sage coss g and he iniial cos funcion, g, are convex quadraic funcions ha can be wrien in he form [ ] [ ][ ] [ ] [ ] u R Q g (u,v) = u s u 2, (2) v M v r v Q (1) g (v) = v P v 2q v, (3) where marices R, Q, M, P R n,n provide he convexiy. he moivaion for considering he QP formulaion and nonlinear esimaion problem examples ha lead o such formulaion are presened below. In he absence of consrains, he ime saged QP is a sligh generalizaion of he opimal linear Gaussian esimaion problem. In linear Gaussian esimaion g (z 1,z ) = 1 2 y C z 2 Φ z A z 1 2 Ψ, (4) he auhors are wih he Informaion Sysems Laboraory, Elecrical Engineering Deparmen, Sanford Universiy, Sanford, CA , USA. {echu58,arezou,gorin,boyd}@sanford.edu. where y is he observed daa; z is he esimaed sae; A and C are he sae updae and he observaion marices, respecively; we use he noaion x 2 Q = x Qx, and Φ, Ψ are he inverse covariance marices of he observaion and process noise, respecively. As grows, he soluion for problem (4) can be recursively compued using he wellknown Kalman filer approach. his paper pursues a moving horizon esimaion (MHE) approach o solve he ime-saged QP (1). For each he MHE solves he problem V N (z N )+ = N+1 g (z 1,z ) F eq z = G eq z 1 +h eq F ineq z G ineq z 1 +h ineq = N +1,...,, (5) where N is fixed, and he ime inerval [ N +1, ] is called he receding horizon of lengh N. he arrival cos V N (z N ) in (5) is compued hrough a recursive updae. he arrival cos updae is he key conribuion of his paper and is described below. Unlike he recursive Kalman filer soluion, which is exac, he MHE soluion is approximae, because he arrival cos updae is approximae. his paper proposes a very accurae approximaion of he arrival cos. B. Moivaion ime-saged QP problems of he form (1) arise in many nonlinear esimaion applicaions, such as roboics, image processing, process monioring, navigaion, and oher filering and rend esimaion applicaions. his QP formulaion includes sae and observaion consrains, which makes i more versaile han he sandard Kalman filer. An imporan source of he ime-saged QP formulaions are problems relaed o compressed sensing or robus esimaion. Such problem impose penaly funcions o shape he error residuals. Common examples include he Huber and l 1 penaly funcion. One well-known example is oal-variaion denoising in which an l 1 penaly is used o promoe sparsiy in labeling daa segmens. We consider he oal-variaion denoising example in Secion III-B of his paper. ime-saged QPs also arise in order-resriced inference problems and bioequivalence problems in drug esing. For increasing, hese problems evenually require runcaing he size of he QP problem solved. he MHE offers a sysemaic approach o such sreaming daa processing. A brief discussion of some earlier work on MHE is presened in he nex subsecion. I is worh noing ha many (perhaps even mos) examples in he published MHE work consider QP formulaions ha can be described by (1).
2 C. Previous work Some of he imporan earlier work on MHE includes [1], [2], [3], [4], where he MHE problems and soluions are esablished as he (esimaion) dual o receding horizon conrol, or model predicive conrol (MPC). In [5] and laer work, he sabiliy of MPC is esablished by inroducing a quadraic approximaion for he cos-o-go a he erminal sae of he receding horizon. As a dual problem, he MHE abiliy o rack he full horizon soluion (sabiliy) requires adding a quadraic arrival cos penaly a he beginning of he receding horizon. For MPC, his line of reasoning is well discussed in [6], where furher references can be found. For MHE, a quadraic approximaion of he arrival cos leads o an exended Kalman filer (EKF) updae for he arrival cos. Such an updae is used in much of he exising work in MHE. In he presence of consrains, however, he quadraic approximaion of he arrival cos used in EKF migh be inaccurae. his is especially rue if he inequaliy consrains swich beween being acive and inacive a he opimal soluion. As a resul, uning he EKF updae in MHE is more of an ar han science. I is well recognized ha improperly uned arrival cos updae can render he MHE updae unsable [7]. his paper proposes a new approach o approximaing he arrival cos ha is aware of he consrains changing from acive o inacive as he opimizaion horizon moves. he earlier paper [8] discussed a special case of he MHE approach proposed in his paper in applicaion o isoonic regression esimaion. he MHE arrival cos updae in [8] uses knowledge abou he consrains being acive or inacive based on he previous MHE updae. his paper develops an exension of his idea o a general class of QP-represenable MHE problems. D. Conribuion of he paper his paper esablishes a class of saged QP problems (1) (3). his class includes many MHE problems in paricular he problems wih l 1 regularizaion where he opimal MHE esimaion allows us o perform on-line segmenaion of sreaming daa. wo examples of such problems are considered in Secion III. he conribuion of he paper is wofold: firs, we inroduce he approximaion hypohesis ha he opimal soluion of a ime-saged QP problem deep inside he horizon (i.e., far back enough in he filer hisory) is insensiive o he las sage of he problem. In paricular, he acive se of he inequaliy consrains becomes fixed afer a paricular ime as he problem horizon grows. his hypohesis holds for mos problems in pracice. Second, we inroduce an ieraive updae of he arrival cos in he MHE soluion. he ieraive updae uses a quadraic approximaion of he arrival cos and acive/inacive consrain informaion from he previous sep. In he absence of he consrains, he updae becomes he familiar Kalman filer in informaion form. In he presence of he consrains, he updae requires solving linear sysems of varying size. his paper demonsraes ha our approximaion hypohesis holds well in pracice and he proposed arrival cos updae leads o good esimaes in our numerical examples. A heoreical proof of our approximaion hypohesis is no included in his paper; noe ha sabiliy proofs for MPC appear well afer MPC s adopion. II. MOVING HORIZON ESIMAION A. Exac arrival cos funcion Since minimizing he full problem is equivalen o minimizing over a parially d problem, he full horizon problem (1) (3) can be presened in he form (5). his requires using an arrival cosv N (z N ) = V N (z N) obained by parial minimizaion wih respec o he decision variables z,z 1,...,z N 1. V N (z N) = ( min g (z )+ ) = N z =1 g (z 1,z ),z 1,...,z N 1 F eq z = G eq z 1 +h eq F ineq z G ineq z 1 +h ineq = 1,..., N We will furher call V N (z N) in (6) he exac arrival cos a N o disinguish i from he approximae arrival cos inroduced below. A recursive updae for he exac arrival cos from 1 o can be obained by exending he parial minimizaion o he nex decision variable z. his yields he Bellman equaion V (u) = min v V 1(v)+g (u,v) F eq v = G eq v +h eq F ineq v G ineq v +h ineq, where u corresponds o z and v o z 1. B. Arrival cos and acive consrains Because of he inequaliy consrains in (6), he exac arrival cos V (u) is a complicaed piecewise quadraic funcion of is argumen. ransiions beween he quadraic segmens corresponds o he inequaliy consrains in (1) swiching beween being acive and inacive. Calculaing he arrival cos funcion in (6) involves solving an inequaliy consrained QP in he decision variablesz,z 1,...,z N 1. As he number N of decision variables grows, in mos cases we can no longer effecively calculae V N (u) in (6) for all u R n. he proposed approach is o use quadraic approximaion of he arrival cos in he viciniy of he opimal soluion Z () = (z,...,z ). If Z () were known ahead of ime, hen we could replace he affine inequaliy consrains in (6) wih equaliy consrains whenever we know ha he consrain is igh for Z (). When he inequaliy consrain is no igh for he soluion Z (), we could drop he consrain. For he consrains in (6), consider he index se A of he inequaliy consrains ha are acive a sep for he (6) (7) 2
3 soluion Z () of he full-horizon problem wih he horizon. he acive consrain se A can be described as A = {i F ineq i z = Gineq i z 1 +hineq i }, (8) for = 1,...,. Noe ha he he se of acive consrains A depends on he horizon. As grows, he se A of he consrains acive a ime sep may change. By replacing he consrains in (6) wih he affine equaliy consrains, he arrival cos becomes a soluion of equaliy consrained QP - a quadraic funcion. Unforunaely, finding he se of acive consrains A exacly requires solving he full-horizon inequaliy consrained QP in he firs place. We will nex discuss how a quadraic approximaion of he arrival cos funcion and he approximae se of he acive consrains A can be compued recursively wihou he need o solve he full horizon problem. C. Approximaion hypohesis he MHE updae relies on represening he full horizon problem (1) in he form (5), wheren is fixed and he horizon is growing. In wha follows, we assume ha N is fixed and > N. For each, he arrival cos is compued based on he MHE soluion and he arrival cos ha have been obained a he previous sep, for 1. Afer ha, he new MHE soluion o he N-sep QP problem is compued. he proposed MHE updae is based on he following approximaion hypohesis: for N large enough and N, A = A +d for any d >. his approximaion implies ha, as grows, he acive consrain se for N does no change and ha A = A, where A is independen of for large enough. Assuming ha A is known, he combined equaliy and inequaliy consrains in (6) for = 1,..., N can be changed o he equaliy consrains F z = G z 1 +h, (9) where he marices F, G, and h are defined as [ ] [ ] [ ] F eq G eq h eq F =, G =, h =, (1) F ineq A G ineq A h ineq A and F A denoes he sub-marix of F formed from he rows F i for all i A. he approximaion hypohesis is a sylized fac ha is no exacly rue, bu appears o hold in pracice for reasonably large N. As is shown below, i allows us o obain very accurae MHE approximaions of he exac full-horizon soluion. D. Updae of he approximae arrival cos Assuming ha he inroduced approximaion hypohesis holds and A is known, he approximae arrival cos V can be found via an approximae Bellman ieraion. By replacing he consrains in (7) wih he equaliy consrains (9), we ge he approximae Bellman updae as V (v) = min u V 1 (u)+g (u,v) F 1 v = G 1 u+h 1 (11) A each horizon, he proposed MHE approach solves he problem (5), where he (approximae) arrival cos V N is updaed according o (11). he arrival cos updae (11) depends on he acive consrain se A N hrough (1). A each horizon, we ake A N = A N = A N 1, he acive se for he soluion ha was compued a horizon 1. Such updae of he acive consrain se follows he approximaion hypohesis described in Secion II-C. Unlike he exac Bellman updae (7) ha has inequaliy consrains, he approximae updae (11) has affine equaliy consrains only. Saring wih V (v) = g (v), he approximae arrival cos funcion obained in such updae is always quadraic and can be presened in he form V (v) = v P v 2q v +α, (12) where marix P and vecor q fully define he approximae arrival cos; he scalar consan α has no impac on he opimizaion problem. he approximae arrival cos updae (11) can be formulaed as an updae for P, q in (12). Consider he parial minimizaion problem for he approximae Bellman ieraion (11), wherev 1 is given by (12) andg (u,v) is given by (2). he Karush-Kuhn-ucker (KK) condiions for opimaliy yield P 1 +R Q F Q M G u v = q 1 +s r, F G µ h Eliminaing he variables u and µ from hese condiions yields a quadraic form in v ha describes V (v) (12). he expressions forp andq can be exraced from his quadraic form. he obained recursive updae for P and q has he form ] Q [ P = M [ P 1 +R G F G ] Q [ q = r [ P 1 +R G F G ] 1 [ ] Q F ] 1 [ ] q 1 +s. h his, along wih (1) and he acive consrain se A N updae, fully specifies he arrival cos updae (11). In conras o he familiar Kalman filer, he marices updaed and invered in our equaions have varying size. he size of he marices F, G, and h changes depending on he number of he inequaliy consrains ha are acive a ime. E. MHE updae summary o recap, he proposed MHE approach approximaely solves he ime-saged QP problem (1) (3). A each horizon, i provides a soluion Z mhe () ha approximaes he opimal soluion Z (). For N, we solve he full problem (1). For > N, he N-sep QP problem (5) is solved wih he approximae arrival cos funcion V N (z N ) of he form (12), where marices P N and q N are defined recursively as described in Secion II-D. he soluion o he QP problem (5) can be obained using any suiable QP solver. For > N he N-sep soluion o (5) defines he las N decision vecors z mhe in Z mhe () = 3
4 (z mhe,zmhe 1,...,zmhe ). he decision vecors for < N are aken from Z mhe ( 1) such ha z mhe = zmhe 1. A. Kalman filer III. EXAMPLES Consider a linear Gaussian esimaion problem for a dynamical sysem x +1 = A x +w y = C x +v, where x is he ime-dependen sae vecor, y is he observaion vecor, and A and C are marices of appropriae sizes. he noise vecors w and v are assumed o be Gaussian wih disribuion N(,Ψ 1 ) and N(,Φ 1 ), respecively. he iniial sae x is assumed o have disribuion N(P 1 q,p 1 ). he problem of maximum likelihood esimaion of he sae z, given he observaions y has he form (1) wih no consrains, where g (z 1,z ) is described by (4), and can be presened in he form (2) wih M = Ψ +C Φ C, Q = Ψ A, R = A Ψ A, r = C Φ y, s =. In his example, here are no consrains and all approximaions are exac. he updae (7) for he (exac) arrival cos funcion V can be performed by he recursive ieraion P = Ψ +C Φ C Ψ A (P 1 +A Ψ A ) 1 A Ψ q = C Φ y +Ψ A (P 1 +A Ψ A ) 1 q 1. In his case, our approximaion hypohesis holds exacly (since here are no consrains), and any choice of N in he MHE formulaion (5) wih he exac arrival cos will yield he same soluion. hese recursive esimaion equaions provide informaion form of he Kalman filer, where P is he inverse covariance marix (informaion marix). he sandard Kalman filer corresponds o he updae wih N =, and provides he las poin of he opimal full-horizon soluion as z = P 1 q. B. oal variaion denoising Consider he esimaion problem x Px 2q x + =1 λ x x = y Cx 2 Φ, (13) where x R n, =,..., are he decision variables, y R m, =,..., are he raw daa, C R m n is he observaion marix, and Φ R m m is a posiive definie weigh marix. he l 1 -norm is used o induce sparsiy of changes in he soluion x and force segmenaion of he soluion ino several inervals wih consan x. Problem (13) is known as oal-variaion denoising and is imporan for many applicaions. o implemen he MHE updae of Secion II-E, we formulae an equivalen problem x Px 2q x + =1 λ1 a + = (y Cx ) Φ(y Cx ) a x x 1 a x x 1 = 1,...,. Using decision variable vecorsz = (x,a ), we can express his as he ime-saged QP problem z P z 2q z + =1 g (z 1,z ) F ineq z G ineq z 1, = 1,...,, where [ ] [ ] P +C P = ΦC q +C, q = Φy, [ ] [ ] I I I F ineq =, G I I ineq =. I he cos funcion g has he form (2) wih [ ] C M = ΦC, Q = R =, [ ] C r = Φy, s (λ/2)1 =. Using his represenaion, we solve an insance of problem (13) wih scalar sae x and scalar daa y, wih C = 1, Φ = 1,P =,q =, and = 2. We use he moving horizon of N = 2 and λ = 2. he ime series y, =,..., was obained by adding random noise o a piecewise consan signal. We implemen he MHE updae of Secion II-E o solve his oal variaion denoising problem. Figure 1 shows he raw daa y as dos and he sae x in he full-horizon soluion Z () for = 2 as a piece-wise consan dashed line. he uneven solid line in he plo corresponds o he x-componen of he final sae z mhe esimaed a each sep of he MHE updae. In his plo, he MHE soluion overlaps wih he full horizon soluion, we have z mhe = z wihin he accuracy of he compuaions. he verical lines in Figure 1 delineae change poins where some of he consrains change beween acive and inacive. C. l 1 rend filering We now consider he l 1 rend filering problem described in [9]. he problem can be formulaed as x Px 2q x + =2 λ x 2x 1 +x 2 + = (y x ) 2, (14) where x is a scalar decision variable and ime series y, =,..., represens he daa o be filered. his problem places an l 1 penaly on he second difference x 2x 1 +x 2 o make i sparse. he soluions of (14) follow a piecewise linear rend wih sparse kink poins. he segmenaion of he soluion ino affine inervals (sraigh 4
5 Fig. 1. MHE filering resuls for oal variaion denoising line segmens) faciliaes daa inerpreaion, compression, and forecasing. Similar o he previous subsecion, he problem can be ransformed ino a ime-saged QP. An equivalen form of problem (14) is x Px 2q x + =2 λa + = (y x ) 2 x 2x 1 +x 2 a = 2,...,. Wih he decision variable z = (x,x 1,a ), we can express he problem as z P z 2q z + =1 g (z 1,z ) F ineq z G ineq z 1 F eq z = G eq z 1, = 2,...,, where P = P +1, q = q +y [ ] [ ] F ineq =, G ineq 1 1 =, F eq = [ 1 ], G eq = [ 1 ]. he represenaion (2) for he cos funcion g is given by M = 1, Q = R =, r = y, s =. (λ/2) Figure 2 shows he resul of applying he MHE algorihm of Secion II-E o problem (14) wih λ = 5, P =, q =. he daa y, =,..., is shown as dos. he solid line shows he opimal filering soluion z. In his case, he MHE soluion z mhe = z (wihin he accuracy of he compuaions) for all =,...,2. he verical lines in Figure 2 delineae he change poins Fig. 2. MHE filering resuls for l 1 filer. = 5 = x x x x 2 4 Fig. 3. he rue V (x) (dashed) and he approximae V(x) (solid) arrival coss. IV. ACCURACY OF APPROXIMAION A. Accuracy of arrival cos approximaion Consider he oal-variaion denoising problem (13) discussed as Example 2 in Secion III-B. his is a onedimensional example wih he scalar sae x, scalar observaion y, C = 1, Φ = 1, and λ = 2. In his example, an MHE updae wih window of N = 2 is applied o daa ime series of he lengh = 2. We calculae he exac arrival cos funcion V ( ), =,...,2 by gridding, and compare i o he approximae arrival cos funcions, obained using our mehod. Figure 3 compares he exac arrival cos funcion V (x) and he proposed quadraic approximaion of he arrival cos funcion V (x) used in he MHE. he plos show hese wo value funcions for = 5 (lef) and = 15 (righ). On he x-axis, we plo he deviaion of he sae variable from he exac full horizon soluion. he y-axis is he arrival cos funcion value offse by a consan (which is inconsequenial for he opimizaion problem). We observe ha he exac arrival cos funcion V and he approximae arrival cos funcion V mach very well in he viciniy of he operaing poin x =. B. Accuracy of segmenaion Consider he l 1 rend filering problem (14), discussed in Secion III-C. In his formulaion, he proposed MHE updae algorihm provides a piecewise linear segmenaion of he underlying rend from sreaming daa. he segmenaion is 5
6 x mhe x Fig. 5. Comparison of he error beween he full horizon esimaes (x ) and he moving horizon esimaes (x mhe ) Fig. 4. Segmenaion for l 1 rend filering: a (blue) do indicaes a change poin where our MHE soluion agrees wih he full horizon soluion, a (red) x represens false negaives where our MHE soluion failed o deec a change poin, and a (green) o represens false posiives where our MHE soluion erroneously deeced a change poin. described by he change poins (kinks) in he soluion. Using he proposed MHE updae, we are able o obain online segmenaion of sreaming daa as he problem size grows ou of bounds for a bach soluion. Figure 4 shows he change poins deeced using he proposed MHE algorihm for he same simulaion resuls as shown in Figure 2. he change poins were deeced as poins where x 2x 1 + x 2 > ǫ, where ǫ is defined by he accuracy of he QP opimizaion soluion. he dos show he deeced change poins. he poins where he full-horizon soluion Z () and our MHE algorihm do no agree are indicaed wih an x (false negaive) or an o (false posiive). he horizonal axis represens he horizon lengh, while he verical axis shows ime inside he full horizon [1, ]. he doed diagonal represens he causaliy boundary: he soluion is obained for. he MHE filering soluion in Figure 4 does mach he full horizon soluion exremely accuraely. he differing change poins have he soluion ha is very close o he consrain, close o he QP solver accuracy limi. C. Accuracy of opimal soluion approximaion Checking he accuracy of he arrival cos approximaion in Secion IV-A needs o be complemened by verificaion of opimal soluion approximaion. Furhermore, if decision vecor z has dimension larger han 2, visualizing he value funcion is no pracical. We compare he performance of he MHE wih proposed arrival cos approximaion agains he performance of an esimaor ha solves he full-horizon problem a each sep. We again consider he oal variaion denoising example (13), bu we use a larger problem dimension, wih x and y having size of n = m = 5, C = I, Φ = I, and λ = 2. We consider an MHE updae wih window of N = 5 applied o he daa ime series y of he lengh = 2. Figure 5 shows he resuls. For each, we plo all n = 5 enries of he esimae obained. he las sae esimaed wih our MHE updae x mhe is ploed along he horizonal axis. he esimaion error x mhe x is shown on he verical axis. he small error observed is on he order of accuracy of he QP solver. Furhermore, he numerical QP soluion is more accurae for he MHE soluion ha requires o solve he problem of small size (5 seps) compared o he fullhorizon soluion (2 seps). his resul means ha he MHE performance wih he proposed approximaion of he arrival cos is pracically idenical o he full-horizon esimaor, a leas in his example. V. CONCLUSION his paper has presened a moving horizon esimaor (MHE) for saged QP problems suiable for sreaming daa processing. he MHE updae is based on he approximaion assumpion ha he acive consrain se remains fixed far enough in he pas from he las ime sage. his assumpion allows us o derive an explici Bellman updae for he arrival cos. he updae is a generalizaion of Kalman filer updae; wih a difference ha a each sep he dimensions of he updaed marices change depending on he number of acive inequaliy consrains. REFERENCES [1] H. Michalska and D. Mayne, Moving horizon observers and observerbased conrol, IEEE ransacions on Auomaic Conrol, no. 6, pp , Jun [2] P. Moraal and J. Grizzle, Observer design for nonlinear sysems wih discree-ime measuremens, IEEE ransacions on Auomaic Conrol, no. 3, pp , Mar [3] C. Rao, J. Rawlings, and D. Mayne, Consrained sae esimaion for nonlinear discree-ime sysems: Sabiliy and moving horizon approximaions, IEEE ransacions on Auomaic Conrol, no. 2, pp , 23. [4] G. C. Goodwin, M. M. Seron, and J. A. D. Dona, Consrained Conrol and Esimaion - An Opimisaion Approach. Springer Verlag, 25. [5] H. Michalska and D. Mayne, Receding horizon conrol of nonlinear sysems, IEEE ransacions on Auomaic Conrol, no. 5, pp , May 199. [6] J. Primbs and V. Nevisic, A new approach o sabiliy analysis of finie receding horizon conrol wihou end consrain, IEEE ransacions on Auomaic Conrol, no. 8, pp , Aug. 2. [7] J. H. Lee and N. L. Ricker, Exended kalman filer based nonlinear model predicive conrol, Indusrial and Engineering Chemisry Research, vol. 33, no. 6, pp , [8] D. Gorinevsky, Efficien filering using monoonic walk model, in American Conrol Conference, Jun. 28, pp [9] S. J. Kim, K. Koh, S. Boyd, and D. Gorinevsky, l 1 rend filering, SIAM Review, no. 2, pp , 29. 6
Notes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationAugmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004
Augmened Realiy II Kalman Filers Gudrun Klinker May 25, 2004 Ouline Moivaion Discree Kalman Filer Modeled Process Compuing Model Parameers Algorihm Exended Kalman Filer Kalman Filer for Sensor Fusion Lieraure
More informationOnline Convex Optimization Example And Follow-The-Leader
CSE599s, Spring 2014, Online Learning Lecure 2-04/03/2014 Online Convex Opimizaion Example And Follow-The-Leader Lecurer: Brendan McMahan Scribe: Sephen Joe Jonany 1 Review of Online Convex Opimizaion
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationTwo Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017
Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =
More informationTechnical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.
Technical Repor Doc ID: TR--203 06-March-203 (Las revision: 23-Februar-206) On formulaing quadraic funcions in opimizaion models. Auhor: Erling D. Andersen Convex quadraic consrains quie frequenl appear
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationLecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples
EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,
More information1 Review of Zero-Sum Games
COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any
More informationSequential Importance Resampling (SIR) Particle Filter
Paricle Filers++ Pieer Abbeel UC Berkeley EECS Many slides adaped from Thrun, Burgard and Fox, Probabilisic Roboics 1. Algorihm paricle_filer( S -1, u, z ): 2. Sequenial Imporance Resampling (SIR) Paricle
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationComparing Means: t-tests for One Sample & Two Related Samples
Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationTimed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1.
Timed Circuis Asynchronous Circui Design Chris J. Myers Lecure 7: Timed Circuis Chaper 7 Previous mehods only use limied knowledge of delays. Very robus sysems, bu exremely conservaive. Large funcional
More informationLinear Gaussian State Space Models
Linear Gaussian Sae Space Models Srucural Time Series Models Level and Trend Models Basic Srucural Model (BSM Dynamic Linear Models Sae Space Model Represenaion Level, Trend, and Seasonal Models Time Varying
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationOn-line Adaptive Optimal Timing Control of Switched Systems
On-line Adapive Opimal Timing Conrol of Swiched Sysems X.C. Ding, Y. Wardi and M. Egersed Absrac In his paper we consider he problem of opimizing over he swiching imes for a muli-modal dynamic sysem when
More informationOutline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests
Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange
More informationArticle from. Predictive Analytics and Futurism. July 2016 Issue 13
Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning
More informationApplication of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing
Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3
More informationReferences are appeared in the last slide. Last update: (1393/08/19)
SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationMost Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation
Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,
More informationOptimal Path Planning for Flexible Redundant Robot Manipulators
25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering
More informationA Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs
PROC. IEEE CONFERENCE ON DECISION AND CONTROL, 06 A Primal-Dual Type Algorihm wih he O(/) Convergence Rae for Large Scale Consrained Convex Programs Hao Yu and Michael J. Neely Absrac This paper considers
More informationDecentralized Stochastic Control with Partial History Sharing: A Common Information Approach
1 Decenralized Sochasic Conrol wih Parial Hisory Sharing: A Common Informaion Approach Ashuosh Nayyar, Adiya Mahajan and Demoshenis Tenekezis arxiv:1209.1695v1 [cs.sy] 8 Sep 2012 Absrac A general model
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 0.038/NCLIMATE893 Temporal resoluion and DICE * Supplemenal Informaion Alex L. Maren and Sephen C. Newbold Naional Cener for Environmenal Economics, US Environmenal Proecion
More informationSliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game
Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp
More informationProbabilistic Robotics
Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae
More informationA Shooting Method for A Node Generation Algorithm
A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan
More informationACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin
ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model
More informationDistribution of Estimates
Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion
More informationA DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS
A DELAY-DEPENDENT STABILITY CRITERIA FOR T-S FUZZY SYSTEM WITH TIME-DELAYS Xinping Guan ;1 Fenglei Li Cailian Chen Insiue of Elecrical Engineering, Yanshan Universiy, Qinhuangdao, 066004, China. Deparmen
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More informationAnnouncements. Recap: Filtering. Recap: Reasoning Over Time. Example: State Representations for Robot Localization. Particle Filtering
Inroducion o Arificial Inelligence V22.0472-001 Fall 2009 Lecure 18: aricle & Kalman Filering Announcemens Final exam will be a 7pm on Wednesday December 14 h Dae of las class 1.5 hrs long I won ask anyhing
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More information2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006
2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)
More informationProbabilistic Robotics SLAM
Probabilisic Roboics SLAM The SLAM Problem SLAM is he process by which a robo builds a map of he environmen and, a he same ime, uses his map o compue is locaion Localizaion: inferring locaion given a map
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationRandom Walk with Anti-Correlated Steps
Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationIntroduction to Mobile Robotics
Inroducion o Mobile Roboics Bayes Filer Kalman Filer Wolfram Burgard Cyrill Sachniss Giorgio Grisei Maren Bennewiz Chrisian Plagemann Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel
More informationProbabilistic Robotics SLAM
Probabilisic Roboics SLAM The SLAM Problem SLAM is he process by which a robo builds a map of he environmen and, a he same ime, uses his map o compue is locaion Localizaion: inferring locaion given a map
More informationInstitute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler
MULTIVARIATE TIME SERIES ANALYSIS AND FORECASTING Manfred Deisler E O S Economerics and Sysems Theory Insiue for Mahemaical Mehods in Economics Universiy of Technology Vienna Singapore, May 2004 Inroducion
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationSEIF, EnKF, EKF SLAM. Pieter Abbeel UC Berkeley EECS
SEIF, EnKF, EKF SLAM Pieer Abbeel UC Berkeley EECS Informaion Filer From an analyical poin of view == Kalman filer Difference: keep rack of he inverse covariance raher han he covariance marix [maer of
More informationEKF SLAM vs. FastSLAM A Comparison
vs. A Comparison Michael Calonder, Compuer Vision Lab Swiss Federal Insiue of Technology, Lausanne EPFL) michael.calonder@epfl.ch The wo algorihms are described wih a planar robo applicaion in mind. Generalizaion
More informationAn recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes
WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationSupplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence
Supplemen for Sochasic Convex Opimizaion: Faser Local Growh Implies Faser Global Convergence Yi Xu Qihang Lin ianbao Yang Proof of heorem heorem Suppose Assumpion holds and F (w) obeys he LGC (6) Given
More informationThe Arcsine Distribution
The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationBook Corrections for Optimal Estimation of Dynamic Systems, 2 nd Edition
Boo Correcions for Opimal Esimaion of Dynamic Sysems, nd Ediion John L. Crassidis and John L. Junins November 17, 017 Chaper 1 This documen provides correcions for he boo: Crassidis, J.L., and Junins,
More informationChristos Papadimitriou & Luca Trevisan November 22, 2016
U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream
More informationScheduling of Crude Oil Movements at Refinery Front-end
Scheduling of Crude Oil Movemens a Refinery Fron-end Ramkumar Karuppiah and Ignacio Grossmann Carnegie Mellon Universiy ExxonMobil Case Sudy: Dr. Kevin Furman Enerprise-wide Opimizaion Projec March 15,
More informationInventory Control of Perishable Items in a Two-Echelon Supply Chain
Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationDesigning Information Devices and Systems I Spring 2019 Lecture Notes Note 17
EES 16A Designing Informaion Devices and Sysems I Spring 019 Lecure Noes Noe 17 17.1 apaciive ouchscreen In he las noe, we saw ha a capacior consiss of wo pieces on conducive maerial separaed by a nonconducive
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationA Specification Test for Linear Dynamic Stochastic General Equilibrium Models
Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models
More informationEnergy Storage Benchmark Problems
Energy Sorage Benchmark Problems Daniel F. Salas 1,3, Warren B. Powell 2,3 1 Deparmen of Chemical & Biological Engineering 2 Deparmen of Operaions Research & Financial Engineering 3 Princeon Laboraory
More informationOptimality Conditions for Unconstrained Problems
62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x
More informationLecture 33: November 29
36-705: Inermediae Saisics Fall 2017 Lecurer: Siva Balakrishnan Lecure 33: November 29 Today we will coninue discussing he boosrap, and hen ry o undersand why i works in a simple case. In he las lecure
More informationKalman filtering for maximum likelihood estimation given corrupted observations.
alman filering maimum likelihood esimaion given corruped observaions... Holmes Naional Marine isheries Service Inroducion he alman filer is used o eend likelihood esimaion o cases wih hidden saes such
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationSome Ramsey results for the n-cube
Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain
More informationTime series model fitting via Kalman smoothing and EM estimation in TimeModels.jl
Time series model fiing via Kalman smoohing and EM esimaion in TimeModels.jl Gord Sephen Las updaed: January 206 Conens Inroducion 2. Moivaion and Acknowledgemens....................... 2.2 Noaion......................................
More informationObject tracking: Using HMMs to estimate the geographical location of fish
Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationThe General Linear Test in the Ridge Regression
ommunicaions for Saisical Applicaions Mehods 2014, Vol. 21, No. 4, 297 307 DOI: hp://dx.doi.org/10.5351/sam.2014.21.4.297 Prin ISSN 2287-7843 / Online ISSN 2383-4757 The General Linear Tes in he Ridge
More informationPade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol
Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne
More information0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED
0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. XI Control of Stochastic Systems - P.R. Kumar
CONROL OF SOCHASIC SYSEMS P.R. Kumar Deparmen of Elecrical and Compuer Engineering, and Coordinaed Science Laboraory, Universiy of Illinois, Urbana-Champaign, USA. Keywords: Markov chains, ransiion probabiliies,
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More information12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =
1: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME Moving Averages Recall ha a whie noise process is a series { } = having variance σ. The whie noise process has specral densiy f (λ) = of
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More information