Efficient quadratic penalization through the partial minimization technique
|
|
- Isabella Malone
- 5 years ago
- Views:
Transcription
1 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control Efficient qadratic penalization throgh the partial imization techniqe Aleksandr Y Aravkin Dmitriy Drsvyatskiy Tristan van Leewen Abstract Common comptational problems, sch as parameter estimation in dynamic models and PDE constrained optimization, reqire data fitting over a set of axiliary parameters sbject to physical constraints over an nderlying state Naive qadratically penalized formlations, commonly sed in practice, sffer from inherent ill-conditioning We show that srprisingly the partial imization techniqe reglarizes the problem, making it well-conditioned This viewpoint sheds new light on variable projection techniqes, as well as the penalty method for PDE constrained optimization, and motivates robst extensions In addition, we otline an inexact analysis, showing that the partial imization sbproblem can be solved very loosely in each iteration We illstrate the theory and algorithms on bondary control, optimal transport, and parameter estimation for robst dynamic inference I INTRODUCTION In this work, we consider a strctred class of optimization problems having the form y, fy + g sbject to Ay q Here, f : R n R is convex and smooth, q R n is a fixed vector, and A is a smoothly varying invertible matrix For now, we make no assmptions on g : R d R, thogh in practice, it is typically either a smooth or a simple nonsmooth fnction Optimization problems of this form often appear in PDE constrained optimization [6], [2], [23], Kalman filtering [2], [4], [5], bondary control problems [], [9], and optimal transport [], [3] Typically, encodes axiliary variables while y encodes the state of the system; the constraint Ay q corresponds to a discretized PDE describing the physics Since the discretization Ay q is already inexact, it is appealing to relax it in the formlation A seegly naive relaxation approach is based on the qadratic penalty: y, F y, : fy + g + Ay q 2 2 Here > is a relaxation parameter for the eqality constraints in, corresponding to relaxed physics The classical qadratic penalty method in nonlinear programg proceeds by applying an iterative optimization algorithm to the nconstrained problem 2 ntil some teration criterion is satisfied, then increasing, and repeating the procedre with the previos iterate sed as a warm start For a detailed discssion, see eg [7, Section Department of Applied Mathematics, University of Washington, Seattle, WA saravkin@wed Department of Mathematics, University of Washington, Seattle, WA ddrsv@wed Department of Mathematics, Utrecht University, Utrecht, Nethelands TvanLeewen@nl logf-f * logf-f * iteration iteration Fig : Top panel shows logf k F of L-BFGS with partial imization applied to 2 for the bondary control problem in Section III Each partial imization solves a least sqares problem in y Bottom panel shows logf k F of L-BFGS withot partial imization applied to 2, for the same vales of Both methods are initialized at random, and y that imizes F, Performance of L-BFGS withot partial imization degrades as increases, while performance of L-BFGS with partial imization is insensitive to 7] The athors of [22] observe that this strategy helps to avoid extraneos local ima, in contrast to the original formlation From this consideration alone, the formlation 2 appears to be sefl Conventional wisdom teaches s that the qadratic penalty techniqe is rarely appropriate The difficlty is that one mst allow to tend to infinity in order to force near-feasibility in the original problem ; the residal error Ay q at an optimal pair, y for the problem 2 is at best on the c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
2 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control 2 order of O/ Conseqently, the maximal eigenvale of the Hessian of the penalty term scales linearly with and the problems 2 become comptationally difficlt Indeed, the maximal eigenvale of the Hessian deteres the behavior of nmerical methods gradient descent, qasi-newton far away from the soltion a regime in which most hge scale problems are solved Figre is a simple nmerical illstration of this inherent difficlty on a bondary control problem; see Section III for more details on the problem formlation The bottom panel in the figre tracks progress of the objective fnction in 2 when an L-BFGS method is applied jointly in the variables, y After iterations of L-BFGS, the objective vale significantly increases with increasing, while the actal imal vale of the objective fnction converges to that of, and so hardly changes In other words, performance of the method scales poorly with, illstrating the ill-conditioning In this paper, we show that by sing a simple partial imization step, this complexity blow-p can be avoided entirely The reslting algorithm is perfectly sited for many large-scale problems, where satisfying the constraint Ay q to high accracy is not reqired or even possible The strategy is straightforward: we rewrite 2 as ϕ + g, 3 where the fnction ϕ is defined implicitly by { ϕ fy + Ay q 2} 4 y We will call ϕ the redced fnction Thogh this approach of imizing ot the variable y, sometimes called variable projection, is widely sed eg [7], [], [22], little theoretical jstification for its speriority is known In this work, we show that not only does partial imization perform well nmerically for the problem class 2, bt is also theoretically gronded We prove that srprisingly the Lipschitz constant of ϕ is bonded by a constant independent of Therefore, iterative methods can be applied directly to the formlation 3 The performance of the new method is illstrated sing a toy example top panel of Figre We se L-BFGS to attack the oter problem; solving it within 35 iterations The inner solver for the toy example simply solves the least sqares problem in y The inner problem 4 can be solved efficiently since its condition nmber is nearly independent of When f is a convex qadratic and A is sparse, one can apply sparse direct solvers or iterative methods sch as LSQR [8] More generally, when f is an arbitrary smooth convex fnction, one can apply first-order methods, which converge globally linearly with the rate governed by the condition nmber of the strongly convex objective in 4 Qasi-newton methods or variants of Newton s method are also available; even if g is nonsmooth, BFGS methods can still be applied The otline of the paper is as follows In Section II, we present complexity garantees of the partial imization techniqe In Section III, we nmerically illstrate the overall approach on bondary control and optimal transport problems, and on tning an oscillator from very noisy measrements II THEORY In this section, we show that the proposed framework is insensitive to the parameter Throghot we assme that f and A are C 2 -smooth, f is convex, and A is invertible for every R n To shorten the formlas, in this section, we will se the symbol A instead of A throghot Setting the stage, define the fnction A qick comptation shows ϕ, y fy + 2 A y q 2 y ϕ, y fy + A T A y q, ϕ, y G, y T A y q, where G, y is the Jacobian with respect to of the map A y Clearly the Lipschitz constant Lip ϕ scales with This can be detrimental to nmerical methods For example, basic gradient descent will find a point k, y k satisfying ϕ k, y k 2 < ɛ after at most O Lip ϕϕ,x ϕ ɛ iterations [6, Section 23] As discssed in the introdction, imizing ϕ amonts to the imization problem ϕ for the redced fnction ϕ defined in 4 Note since f is convex and A is invertible, the fnction ϕ, admits a niqe imizer, which we denote by y Appealing to the classical implicit fnction theorem eg [2, Theorem 58] we dedce that ϕ is differentiable with ϕ ϕ, y G, y T A y q We aim to pper bond the Lipschitz constant of ϕ by a qantity independent of We start by estimating the residal A y q Throghot the paper, we se the following simple identity Given an invertible map F : R n R n and invertible matrix C, for any points x R n and nonzero R, we have F x F x, and 2 C F C T C T F C We often apply this observation to the invertible map F x fx + Bx, where B is a positive definite matrix Lemma Residal bond For any point, the ineqality holds: f A q A y q Proof Note that the imizers y of ϕ, are characterized by first order optimality conditions 5 fy + A T A y q 6 Applying the implicit fnction theorem, we dedce that y depends C 2 -smoothly on with y given by 2 fy + A T A A T A y q On the other hand, from the eqality 6 we have y f + A T A A T f q + AT A A T q c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
3 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control 3 Therefore, we dedce A y q A f + AT A A T q q A T f A + I I q Define now the operator F : A T and the point z : F q Note that F x x f A + I f A x Letting L be a Lipschitz constant of F, we obtain A y q F q F z L q z L q F q L A T fa q Now the inverse fnction theorem yields for any point y the ineqality F y F F y 2 f A F y + I, where the last ineqality follows from the fact that by convexity of f A all eigenvales of 2 f A are nonnegative Ths we may set L, completing the proof For ease of reference, we record the following direct corollary Corollary For any point, we have y A q A A y q A f A q Next we will compte the Hessian of ϕ, y, and se it to show that the norm of the Hessian of ϕ is bonded by a constant independent of Defining ] ] R, y, v [G, y T v and K, v [A T v, we can partition the Hessian as follows: [ ] 2 ϕ ϕ ϕ y where ϕ, y ϕ y ϕ yy G, y T G, y + R, y, A y q, ϕ yy, y 2 fy + A T A, ϕ y, y K, A y q + A T G, y See [22, Section 4] for more details Moreover, it is known that the Hessian of the redced fnction ϕ admits the expression [22, Eqation 22] 2 ϕ ϕ, y ϕ y, y ϕ yy, y ϕ y, y, 7 which is simply the Schr complement of ϕ yy, y in 2 ϕ, y We define the operator norms G, y T : sp [G, y T v], v A T : sp [A T v] v Using this notation, we can prove the following key bonds Corollary 2 For any points and y, the ineqalities hold: 2 ϕ yy, y A, f A q R, y, A y q G, y, f A q A T K, A y q Proof The first bond follows by the ineqality ϕ yy, y A A T 2 fya + I A T A 2, and the remaining bonds are immediate from Lemma Next, we need the following elementary linear algebraic fact Lemma 2 For any positive semidefinite matrix B and a real >, we have I I + B 2 B Proof Define the matrix F I I + B and consider an arbitrary point z Observing the ineqality I + B and defining the point p : I + Bz, we obtain F z z I + B z I + B p I + B z I + B p z Bz B z Since this holds for all z, the reslt follows Ptting all the pieces together, we can now prove the main theorem of this section Theorem Norm of the redced Hessian The operator norm of 2 ϕ is bonded by a qantity CA, f, q, independent of Proof To simplify the proof, define G : G, y, R : R, y, A y q, K : K, A y q, and ϕ yy, y When, the operator norm of 2 φ has a trivial bond c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
4 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control 4 directly from eqation 7 For large, after rearranging 7, we can write 2 ϕ as follows: R 2 K T K + K T A T G + G T A K + G T G 2 G T A A T G 8 Corollary 2 implies that the operator norm of the first row of 8 is bonded above by the qantity L G + L2 A T 2 A 2 + 2L A T A 2 A T G, where we set L : f A q Notice that the expression in 9 is independent of We rewrite the second row of 8 sing the explicit expression for : G T G 2 G T A A T G 9 G T G G T A T 2 fy A + I G G T I Applying Lemma 2 with B A T G T G 2 G T A A T G Setting A T 2 fy A + I G 2 fy A, we have G 2 A T 2 fy A CA, f, q, : L G + 2L A T A 2 A T G + G 2 A 2 2 fy + L2 A T 2 A 2 For >, the last term is always trivially bonded by L 2 A T 2 A 2 and the reslt follows A Inexact analysis of the projection sbproblem In practice, one can rarely evalate ϕ exactly It is therefore important to nderstand how inexact soltions of the inner problems 4 impact iteration complexity of the oter problem The reslts presented in the previos section form the fondation for sch an analysis For simplicity, we assme that g is smooth, thogh the reslts can be generalized, as we comment on shortly In this section, we compte the overall complexity of the partial imization techniqe when the oter nonconvex imization problem 3 is solved by an inexact gradient descent algorithm When g is nonsmooth, a completely analogos analysis applies to the prox-gradient method We only focs here on gradient descent, as opposed to more sophisticated methods, since the analysis is straightforward We expect qasi- Newton methods and limited memory variants to exhibit exactly the same behavior eg Figre We do not perform a similar analysis here for inexact qasi-newton methods, as the global efficiency estimates even for exact qasi-newton methods for nonconvex problems are poorly nderstood Define the fnction H : g + ϕ Let β > be the Lipschitz constant of the gradient H g + ϕ Fix a constant c >, and sppose that in each iteration k, we compte a vector v k with v k H k c k Consider then the inexact gradient descent method k+ k β v k Then we dedce H k+ H k H k, β v k + β 2 β v k 2 v 2β k H k 2 H k 2 Hence we obtain the convergence garantee: i,,k H i 2 k k H i 2 i 2β H H k 2β H H k + k k v k H k 2 i + c2 π 2 6k β2 2 + c 2 π 2 /6 k where is sed to go from line to line 2 Now, if we compte g exactly, the qestion is how many inner iterations are needed to garantee v k H k c k For fixed k, the inner objective is ϕ k, y fy + 2 A ky q 2 The condition nmber ratio of Lipschitz constant of the gradient over the strong convexity constant of ϕ k, y in y is κ k : Lipf A k 2 + A k 2 A k 2 Notice that κ k converges to the sqared condition nmber of A k as Gradient descent on the fnction ϕ k, garantees y i y 2 ɛ after κ k log y y 2 ɛ iterations Then we have ϕ k, y i ϕ k, y G, y y i y Since we want the left hand side to be bonded by c k, we simply need to ensre y i y 2 c 2 k 2 2 G k, y 2 Therefore the total nmber of inner iterations is no larger than y κ k log y 2 G k, y 2 c 2 k 2 2, which grows very slowly with k and with In particlar, the nmber of iterations to solve the inner problem scales as logk to achieve a global k rate in H 2 If instead we se a fast-gradient method [6, Section 22] for imizing ϕ k,, we can replace κ k with the mch better qantity κ k throghot c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
5 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control 5 III NUMERICAL ILLUSTRATIONS In this section, we present two representative examples of PDE constrained optimization bondary control and optimal transport and a problem of robst dynamic inference In each case, we show that practical experience spports theoretical reslts from the previos section In particlar, in each nmerical experiment, we stdy the convergence behavior of the proposed method as increases A Bondary control In bondary control, the goal is to steer a system towards a desired state by controlling its bondary conditions Perhaps the simplest example of sch a problem is the following Given a sorce qx, defined on a domain Ω, we seek bondary conditions sch that the soltion to the Poisson problem y q y Ω for x Ω is close to a desired state y d Discretizing the PDE yields the system Ay + B q where A is a discretization of the Laplace operator on the interior of the domain and B coples the interior gridpoints to the bondary The corresponding PDE-constrained optimization problem is given by,y 2 y y d 2 2 sbject to Ay + B q, whereas the penalty formlation reads,y 2 y y d Ay + B q 2 2 Since both terms are qadratic in y, we can qickly solve for y explicitly Nmerical experiments: In this example, we consider an L-shaped domain with a sorce q shaped like a Gassian bell, as depicted in figre 2 Or goal is to get a constant distribtion y d in the entire domain The soltion for is shown in figre 3 a To solve the optimization problem we se a steepest-descent method with a fixed stepsize, detered from the Lipschitz constant of the gradient The reslt of the constrained formlation is shown in figre 3 b We see that by adapting the bondary conditions we get a more even distribtion The convergence behavior for varios vales of is shown in figre 4 a We see that as, the behavior tends towards that of the constrained formlation, as expected The Lipschitz constant of the gradient evalated at the initial point, as a fnction of is shown in figre 4 b; the crve levels off as the theory predicts B Optimal transport The second class of PDE-constrained problems we consider comes from optimal transport, where the goal is to detere a mapping, or flow, that optimally transforms one mass density fnction into another Say we have two density fnctions y x and y T x, with x Ω, we can formlate the problem as Fig 2: L-shaped domain with the sorce fnction q a b Fig 3: Bondary vales,, and soltion in the interior for the initial and optimized bondary vales are depicted in a and b respectively gradient - -2 inf iteration a Lipschitz constant Fig 4: The convergence plots for varios vales of are depicted in a, while b shows the dependence of the nmerically compted Lipschitz constant on t, x finding a flowfield, t, x, sch that y 2 t, x T x yt, x and y x y, x, where yt, x solves y t + y Discretizing sing an implicit Lax-Friedrichs scheme [2], the PDE reads Ay q, b c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
6 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control 6 y data Fig 8: Left: Densities, Gassian black dash, Hber red solid, and Stdent s t ble dot Right: Negative Log Likelihoods gradient Fig 5: The initial left and desired right mass density are shown yt Fig 6: The mass density obtained after optimization left and the corresponding time-averaged flow field right are shown where q contains the initial condition and we have I + tb M I + tb 2 A M M I + tb N with M a for-point averaging matrix and B containing the discretization of the derivative terms Adding reglarization to promote smoothness of and y in time [2], we obtain the problem ,y iteration a 2 P y y T α2 2 yt Ldiag inf sbject to Ay q Lipschitz constant Fig 7: Convergence plots for varios vales of are shown in a, while b shows the nmerically compted Lipschitz constant as a fnction of b, Here, P restricts the soltion y to t T, α is a reglarization parameter and L is a block matrix with I on the main and pper diagonal The penalized formlation is,y 2 P y y T α2 2 yt Ldiag + 2 Ay q 2 2 Again the partial imization in y amont to imizing a qadratic fnction Nmerical experiments: For the nmerical example we consider the domain Ω [, ] 2, discretized with x /6 and T /32 with a stepsize of t /8 The initial and desired state are depicted in figre 5 The reslting state obtained at time T and the corresponding time-averaged flowfield are depicted in figre 6 The initial flow was generated by iid samples from a standard Gassian random variable To imize 2, we sed a steepest-descent method with constant step size, sing the largest eigenvale of the Gass-Newton Hessian at the initial as an estimate of the Lipschitz constant The convergence behavior for varios vales of as well as the corresponding estimates of the Lipschitz constant at the final soltion are shown in figre 7 C Robst dynamic inference with the penalty method In many settings, data is natrally very noisy, and a lot of effort mst be spent in pre-processing and cleaning before applying standard inversion techniqes To narrow the scope, consider dynamic inference, where we wish to infer both hidden states and nknown parameters driven by an nderlying ODE Recent efforts have focsed on developing inference formlations that are robst to otliers in the data [3], [4], [9], sing convex penalties sch as l, Hber [4] and non-convex penalties sch as the Stdent s t log likelihood in place of the least sqares penalty The goal is to develop formlations and estimators that achieve adeqate performance when faced with otliers; these may arise either as gross measrement errors, or real-world events that are not modeled by the dynamics Figre 8 shows the probability density fnctions and penalties corresponding to Gassian, Hber, and Stdent s t densities Qadratic tail growth corresponds to extreme decay of the Gassian density for large inpts, and linear growth of the inflence of any measrement on the fit In contrast, Hber and Stdent s t have linear and sblinear tail growth, respectively, which ensres every observation has bonded inflence We focs on the Hber fnction [4], since it is both C - smooth and convex In particlar, the fnction fy in and 4 is chosen to be a composition of the Hber with an observation model Note that Hber is not C 2, so this case is c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
7 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control 7 not immediately captred by the theory we propose However, Hber can be closely approximated by a C 2 fnction [8], and then the theory flly applies For or nmerical examples, we apply the algorithm developed in this paper directly to the Hber formlation We illstrate robst modeling sing a simple representative example Consider a 2-dimensional oscillator, governed by the following eqations: [ ] [ y ] [ ] [ ] y sinωt + 3 y 2 y 2 where we can interpret as the freqency, and 2 is the damping Discretizing in time, we have [ ] k+ [ y 2 t 2 t 2 ] [ ] k [ ] y sinωtk + t y 2 t y 2 We now consider direct observations of the second component, z k [ ] [ ] k y + w y k, w k N, R k 2 We can formlate the joint inference problem on states and measrements as follows: φ, y :,y 2 ρ R /2 Hy z st Gy v, 4 with v k sinωt k, v the initial condition for y, and ρ is either the least sqares or Hber penalty, and we se the following definitions: R diag{r k } H diag{h k } y vec{y k } z vec{z,, z N } I G G 2 I G N I Note in particlar that there are only two nknown parameters, ie R 2, while the state y lies in R 2N, with N the nmber of modeled time points The redced optimization problem for is given by f : 2 ρr /2 HG v z To compte the derivative of the ODE-constrained problem, we can se the adjoint state method Defining the Lagrangian Ly, x, R 2 ρ /2 Hy z + x, Gy v, we write down the optimality conditions L and obtain y G v x G T H T R /2 ρr /2 HG v z f x, Gy The inexact penalized problem is given by ϕ, y,y 2 ρ R /2 Hy z + 2 Gy v 2 5 TABLE I: Reslts for Kalman experiment Penalty method for both least sqares and hber achieves the same reslts for moderate vales of as does the projected formlation While Hber reslts converge to nearly the tre parameters, least sqares reslts converge to an incorrect parameter estimate ρ Iter Opt 3 l , 98 5 l , 43 7 l , 9 l , 8 l , 8 3 h h , 7 h , 9 h , h , We then immediately find y arg y 2 ρ R /2 Hy z + Gy v 2 2 φ Gy Gy v When ρ is the least sqares penalty, y is available in closed form However, when ρ is the Hber, y reqires an iterative algorithm Rather than solving for y sing a first-order method, we se IPsolve, an interior point method well sited for Hber [5] Even thogh each iteration reqires inversions of systems of size ON, these systems are very sparse, and the complexity of each iteration to compte y is ON for any piecewise linear qadratic fnction [5] Once again, we see that the comptational cost does not scale with Nmerical Experiments: We simlate a data contaation scenario by solving the ODE 3 for the particlar parameter vale 2, The second component of the reslting state y is observed, and the observations are contaated In particlar, in addition to Gassian noise with standard deviation σ, for % of the measrements niformly distribted errors in [, 2] are added The state y R 24 is finely sampled over 4 periods For the least sqares and the Hber penalty with κ, we solved both the ODE constrained problem 4 and the penalized version 5 for { 3, 5, 7, 9 } The reslts are presented in Table I The Hber formlation behaves analogosly to the formlation sing least sqares; in particlar the oter projected fnction in is no more difficlt to imize And, as expected, the robst Hber penalty finds the correct vales for the parameters A state estimate generated from -estimates corresponding to large is shown in Figre 9 The hberized approach is able to ignore the otliers, and recover both better estimates of the nderlying dynamics parameters, and the tre observed and hidden components of the state y IV CONCLUSIONS In this paper, we showed that, contrary to conventional wisdom, the qadratic penalty techniqe can be sed effectively c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
8 This article has been accepted for pblication in a ftre isse of this jornal, bt has not been flly edited Content may change prior to final pblication Citation information: DOI 9/TAC , IEEE Transactions on Atomatic Control Tre L2 recovery Hber recovery measrements penalty composed with a linear model as the fnction fy Nmerical reslts illstrated the overall approach, inclding convergence behavior of the penalized projected scheme, as well as modeling advantages of robst penalized formlations Acknowledgment Research of A Aravkin was partially spported by the Washington Research Fondation Data Science Professorship Research of D Drsvyatskiy was partially spported by the AFOSR YIP award FA The research of T van Leewen was in part financially spported by the Netherlands Organisation of Scientific Research NWO as part of research programme Tre L2 recovery Hber recovery Fig 9: Top panel: first samples of state y 2 solid black, least sqares recovery red dash-dot and hber recovery ble dash Noisy measrements z appear as ble diamonds on the plot, with otliers shown at the top of the panel Bottom panel: first samples of state y solid black, least sqares recovery red dash-dot and Hber recovery ble dash for control and PDE constrained optimization problems, if done correctly In particlar, when combined with the partial imization techniqe, we showed that the penalized projected scheme g + {fy + 2 } Ay b 2 y has the following advantages: The Lipschitz constant of the gradient of the oter fnction in is bonded as, and hence we can effectively analyze the global convergence of first-order methods 2 Convergence behavior of the data-reglarized convex inner problem is controlled by parametric matrix A, a fndamental qantity that does not depend on 3 The inner problem can be solved inexactly, and in this case, the nmber of inner iterations of a first-order algorithm needs to grow only logarithmically with and the oter iterations conter, to preserve the natral rate of gradient descent As an immediate application, we extended the penalty method in [22] to convex robst formlations, sing the Hber REFERENCES [] L Ambrosio and N Gigli A User s Gide to Optimal Transport In Modelling and Optimisation of Flows on Networks, pages 55 Springer, 23 [2] B D Anderson and J B Moore Optimal filtering 979, 979 [3] A Aravkin, B Bell, J Brke, and G Pillonetto An l -Laplace robst Kalman smoother IEEE Transactions on Atomatic Control, 2 [4] A Aravkin, J Brke, and G Pillonetto Robst and trend-following Stdent s t Kalman smoothers SIAM Jornal on Control and Optimization, 525: , 24 [5] A Aravkin, J V Brke, and G Pillonetto Sparse/robst estimation and Kalman smoothing with nonsmooth log-concave densities: Modeling, comptation, and theory Jornal of Machine Learning Research, 4: , 23 [6] A Aravkin, M Friedlander, F Herrmann, and T van Leewen Robst inversion, dimensionality redction, and randomized sampling Mathematical Programg, 34: 25, 22 [7] A Y Aravkin and T van Leewen Estimating nisance parameters in inverse problems Inverse Problems, 28:56, 22 [8] K P Bbe and R T Langan Hybrid l/l2 imization with applications to tomography Geophysics, 624:83 95, 997 [9] S Farahmand, G B Giannakis, and D Angelosante Dobly Robst Smoothing of Dynamical Processes via Otlier Sparsity Constraints IEEE Transactions on Signal Processing, 59: , 2 [] G Golb and V Pereyra Separable nonlinear least sqares: the variable projection method and its applications Inverse Problems, 92:R R26, Apr 23 [] M Ggat Control problems with the wave eqation SIAM Jornal on Control and Optimization, 485:326 35, 29 [2] E Haber and L Hanson Model Problems in PDE-Constrained Optimization Technical report, 27 [3] S Haker, L Zh, A Tannenbam, and S Angenent Optimal Mass Transport for Registration and Warping International Jornal of Compter Vision, 63:225 24, dec 24 [4] P J Hber Robst Statistics John Wiley and Sons, 24 [5] R E Kalman A New Approach to Linear Filtering and Prediction Problems Transactions of the AMSE - Jornal of Basic Engineering, 82D:35 45, 96 [6] Y Nesterov Introdctory lectres on convex optimization, volme 87 of Applied Optimization Klwer Academic Pblishers, Boston, MA, 24 A basic corse [7] J Nocedal and S Wright Nmerical optimization Springer Series in Operations Research and Financial Engineering Springer, New York, second edition, 26 [8] C C Paige and M A Sanders LSQR: An algorithm for sparse linear eqations and sparse least sqares ACM Transactions on Mathematical Software TOMS, 8:43 7, 982 [9] P Philip Optimal Control and Partial Differential Eqations Technical report, LMU Berlin, 29 [2] R Rockafellar and R Wets Variational Analysis, volme 37 Springer, 998 [2] A Tarantola Inverse Problem Theory SIAM, 25 [22] T van Leewen and F J Herrmann A penalty method for pdeconstrained optimization in inverse problems Inverse Problems, 32:57, 25 [23] J Viriex and S Operto An overview of fll-waveform inversion in exploration geophysics Geophysics, 746:WCC WCC26, c 27 IEEE Personal se is permitted, bt repblication/redistribtion reqires IEEE permission See for more information
Modelling by Differential Equations from Properties of Phenomenon to its Investigation
Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationB-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Quadratic Optimization Problems in Continuous and Binary Variables
B-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Qadratic Optimization Problems in Continos and Binary Variables Naohiko Arima, Snyong Kim and Masakaz Kojima October 2012,
More informationDirect linearization method for nonlinear PDE s and the related kernel RBFs
Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationA Model-Free Adaptive Control of Pulsed GTAW
A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,
More informationStudy of the diffusion operator by the SPH method
IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationNonlinear parametric optimization using cylindrical algebraic decomposition
Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More informationA Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model
entropy Article A Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model Hyenkyn Woo School of Liberal Arts, Korea University of Technology and Edcation, Cheonan
More informationOPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1
OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body
More informationStability of Model Predictive Control using Markov Chain Monte Carlo Optimisation
Stability of Model Predictive Control sing Markov Chain Monte Carlo Optimisation Elilini Siva, Pal Golart, Jan Maciejowski and Nikolas Kantas Abstract We apply stochastic Lyapnov theory to perform stability
More informationDownloaded 07/06/18 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 4, No., pp. A4 A7 c 8 Society for Indstrial and Applied Mathematics Downloaded 7/6/8 to 8.83.63.. Redistribtion sbject to SIAM license or copyright; see http://www.siam.org/jornals/ojsa.php
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
More informationWhen are Two Numerical Polynomials Relatively Prime?
J Symbolic Comptation (1998) 26, 677 689 Article No sy980234 When are Two Nmerical Polynomials Relatively Prime? BERNHARD BECKERMANN AND GEORGE LABAHN Laboratoire d Analyse Nmériqe et d Optimisation, Université
More informationA Survey of the Implementation of Numerical Schemes for Linear Advection Equation
Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear
More informationAn Investigation into Estimating Type B Degrees of Freedom
An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More informationSources of Non Stationarity in the Semivariogram
Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary
More informationFRÉCHET KERNELS AND THE ADJOINT METHOD
PART II FRÉCHET KERNES AND THE ADJOINT METHOD 1. Setp of the tomographic problem: Why gradients? 2. The adjoint method 3. Practical 4. Special topics (sorce imaging and time reversal) Setp of the tomographic
More informationSubstructure Finite Element Model Updating of a Space Frame Structure by Minimization of Modal Dynamic Residual
1 Sbstrctre Finite Element Model Updating of a Space Frame Strctre by Minimization of Modal Dynamic esidal Dapeng Zh, Xinn Dong, Yang Wang, Member, IEEE Abstract This research investigates a sbstrctre
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationJoint Transfer of Energy and Information in a Two-hop Relay Channel
Joint Transfer of Energy and Information in a Two-hop Relay Channel Ali H. Abdollahi Bafghi, Mahtab Mirmohseni, and Mohammad Reza Aref Information Systems and Secrity Lab (ISSL Department of Electrical
More informationOn Multiobjective Duality For Variational Problems
The Open Operational Research Jornal, 202, 6, -8 On Mltiobjective Dality For Variational Problems. Hsain *,, Bilal Ahmad 2 and Z. Jabeen 3 Open Access Department of Mathematics, Jaypee University of Engineering
More informationDecision Oriented Bayesian Design of Experiments
Decision Oriented Bayesian Design of Experiments Farminder S. Anand*, Jay H. Lee**, Matthew J. Realff*** *School of Chemical & Biomoleclar Engineering Georgia Institte of echnology, Atlanta, GA 3332 USA
More informationReducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation
JOURNAL OF AIRCRAFT Vol. 42, No. 4, Jly Agst 2005 Redcing Conservatism in Fltterometer Predictions Using Volterra Modeling with Modal Parameter Estimation Rick Lind and Joao Pedro Mortaga University of
More informationLecture 8: September 26
10-704: Information Processing and Learning Fall 2016 Lectrer: Aarti Singh Lectre 8: September 26 Note: These notes are based on scribed notes from Spring15 offering of this corse. LaTeX template cortesy
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationNonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator
Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational
More informationAssignment Fall 2014
Assignment 5.086 Fall 04 De: Wednesday, 0 December at 5 PM. Upload yor soltion to corse website as a zip file YOURNAME_ASSIGNMENT_5 which incldes the script for each qestion as well as all Matlab fnctions
More informationDiscussion Papers Department of Economics University of Copenhagen
Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,
More informationDiscussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli
1 Introdction Discssion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen Department of Economics, University of Copenhagen and CREATES,
More informationAN IMPROVED VARIATIONAL MODE DECOMPOSITION METHOD FOR INTERNAL WAVES SEPARATION
3rd Eropean Signal Processing Conference (EUSIPCO AN IMPROVED VARIATIONAL MODE DECOMPOSITION METHOD FOR INTERNAL WAVES SEPARATION Jeremy Schmitt, Ernesto Horne, Nelly Pstelni, Sylvain Jobad, Philippe Odier
More informationThe spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam
J.L. Pan W.Y. Zh Nonlinear Sci. Lett. Vol.8 No. pp.- September The spreading reside harmonic balance method for nonlinear vibration of an electrostatically actated microbeam J. L. Pan W. Y. Zh * College
More informationL 1 -smoothing for the Ornstein-Uhlenbeck semigroup
L -smoothing for the Ornstein-Uhlenbeck semigrop K. Ball, F. Barthe, W. Bednorz, K. Oleszkiewicz and P. Wolff September, 00 Abstract Given a probability density, we estimate the rate of decay of the measre
More informationA Macroscopic Traffic Data Assimilation Framework Based on Fourier-Galerkin Method and Minimax Estimation
A Macroscopic Traffic Data Assimilation Framework Based on Forier-Galerkin Method and Minima Estimation Tigran T. Tchrakian and Sergiy Zhk Abstract In this paper, we propose a new framework for macroscopic
More informationDEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1. C. Silvestre A. Pascoal
DEPTH CONTROL OF THE INFANTE AUV USING GAIN-SCHEDULED REDUCED-ORDER OUTPUT FEEDBACK 1 C. Silvestre A. Pascoal Institto Sperior Técnico Institte for Systems and Robotics Torre Norte - Piso 8, Av. Rovisco
More informationA fundamental inverse problem in geosciences
A fndamental inverse problem in geosciences Predict the vales of a spatial random field (SRF) sing a set of observed vales of the same and/or other SRFs. y i L i ( ) + v, i,..., n i ( P)? L i () : linear
More informationLecture 13: Duality Uses and Correspondences
10-725/36-725: Conve Optimization Fall 2016 Lectre 13: Dality Uses and Correspondences Lectrer: Ryan Tibshirani Scribes: Yichong X, Yany Liang, Yanning Li Note: LaTeX template cortesy of UC Berkeley EECS
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationFREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS
7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,
More informationThe Dual of the Maximum Likelihood Method
Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All
More informationFrequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1
Freqency Estimation, Mltiple Stationary Nonsinsoidal Resonances With Trend 1 G. Larry Bretthorst Department of Chemistry, Washington University, St. Lois MO 6313 Abstract. In this paper, we address the
More informationarxiv: v1 [physics.flu-dyn] 4 Sep 2013
THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv:1309.1728v1 physics.fl-dyn]
More informationActive Flux Schemes for Advection Diffusion
AIAA Aviation - Jne, Dallas, TX nd AIAA Comptational Flid Dynamics Conference AIAA - Active Fl Schemes for Advection Diffsion Hiroaki Nishikawa National Institte of Aerospace, Hampton, VA 3, USA Downloaded
More informationFUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS
MODELOWANIE INŻYNIERSKIE ISNN 896-77X 3, s. 433-438, Gliwice 6 FUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS JERZY SKRZYPCZYK HALINA WITEK Zakład
More informationFRTN10 Exercise 12. Synthesis by Convex Optimization
FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be
More informationMove Blocking Strategies in Receding Horizon Control
Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it
More informationLinear and Nonlinear Model Predictive Control of Quadruple Tank Process
Linear and Nonlinear Model Predictive Control of Qadrple Tank Process P.Srinivasarao Research scholar Dr.M.G.R.University Chennai, India P.Sbbaiah, PhD. Prof of Dhanalaxmi college of Engineering Thambaram
More informationOptimal Control, Statistics and Path Planning
PERGAMON Mathematical and Compter Modelling 33 (21) 237 253 www.elsevier.nl/locate/mcm Optimal Control, Statistics and Path Planning C. F. Martin and Shan Sn Department of Mathematics and Statistics Texas
More informationNetwork Coding for Multiple Unicasts: An Approach based on Linear Optimization
Network Coding for Mltiple Unicasts: An Approach based on Linear Optimization Danail Traskov, Niranjan Ratnakar, Desmond S. Ln, Ralf Koetter, and Mriel Médard Abstract In this paper we consider the application
More informationThe Real Stabilizability Radius of the Multi-Link Inverted Pendulum
Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract
More informationarxiv: v1 [cs.sy] 22 Nov 2018
AAS 18-253 ROBUST MYOPIC CONTROL FOR SYSTEMS WITH IMPERFECT OBSERVATIONS Dantong Ge, Melkior Ornik, and Ufk Topc arxiv:1811.09000v1 [cs.sy] 22 Nov 2018 INTRODUCTION Control of systems operating in nexplored
More informationChapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS
Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system
More informationRemarks on strongly convex stochastic processes
Aeqat. Math. 86 (01), 91 98 c The Athor(s) 01. This article is pblished with open access at Springerlink.com 0001-9054/1/010091-8 pblished online November 7, 01 DOI 10.1007/s00010-01-016-9 Aeqationes Mathematicae
More informationA New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration
Proceedings of the 5th WSEAS Int. Conf. on Power Systems and Electromagnetic Compatibility, Corf, Greece, Agst 3-5, 005 (pp43-48) A New Method for Calclating of Electric Fields Arond or Inside Any Arbitrary
More informationKrauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.
Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore
More informationA Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units
A Reglator for Continos Sedimentation in Ideal Clarifier-Thickener Units STEFAN DIEHL Centre for Mathematical Sciences, Lnd University, P.O. Box, SE- Lnd, Sweden e-mail: diehl@maths.lth.se) Abstract. The
More informationBREGMAN ITERATIVE ALGORITHMS FOR l 1 -MINIMIZATION WITH APPLICATIONS TO COMPRESSED SENSING
BREGMAN ITERATIVE ALGORITHMS FOR l 1 -MINIMIZATION WITH APPLICATIONS TO COMPRESSED SENSING WOTAO YIN, STANLEY OSHER, DONALD GOLDFARB, AND JEROME DARBON Abstract. We propose simple and extremely efficient
More informationAn extremum seeking approach to parameterised loop-shaping control design
Preprints of the 19th World Congress The International Federation of Atomatic Control An extremm seeking approach to parameterised loop-shaping control design Chih Feng Lee Sei Zhen Khong Erik Frisk Mattias
More informationTrustSVD: Collaborative Filtering with Both the Explicit and Implicit Influence of User Trust and of Item Ratings
Proceedings of the Twenty-Ninth AAAI Conference on Artificial Intelligence TrstSVD: Collaborative Filtering with Both the Explicit and Implicit Inflence of User Trst and of Item Ratings Gibing Go Jie Zhang
More informationNon-collinear upconversion of infrared light
Non-collinear pconversion of infrared light Christian Pedersen, Qi H, Lasse Høgstedt, Peter Tidemand-Lichtenberg, and Jeppe Seidelin Dam * DTU Fotonik, Frederiksborgvej 399, 4000 Roskilde, Denmark * jdam@fotonik.dt.dk
More informationTheoretical and Experimental Implementation of DC Motor Nonlinear Controllers
Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,
More informationChapter 2 Difficulties associated with corners
Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces
More informationAn optimizing reduced order FDS for the tropical Pacific Ocean reduced gravity model
An optimizing redced order FDS for the tropical Pacific Ocean redced gravity model Zhendong Lo a,, Jing Chen b,, Jiang Zh c,, Riwen Wang c,, and I. M. Navon d,, a School of Science, Beijing Jiaotong University,
More informationOil production optimization of several wells subject to choke degradation
Proceedings of the 3rd IFAC Workshop on Atomatic Control in Offshore Oil and Gas Prodction, Esbjerg, Denmark, May 30 - Jne 1, 2018 We_A_Reglar_Talk.1 Oil prodction optimization of several wells sbject
More informationAffine Invariant Total Variation Models
Affine Invariant Total Variation Models Helen Balinsky, Alexander Balinsky Media Technologies aboratory HP aboratories Bristol HP-7-94 Jne 6, 7* Total Variation, affine restoration, Sobolev ineqality,
More informationSensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters
Sensitivity Analysis in Bayesian Networks: From Single to Mltiple Parameters Hei Chan and Adnan Darwiche Compter Science Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwiche}@cs.cla.ed
More information1 Undiscounted Problem (Deterministic)
Lectre 9: Linear Qadratic Control Problems 1 Undisconted Problem (Deterministic) Choose ( t ) 0 to Minimize (x trx t + tq t ) t=0 sbject to x t+1 = Ax t + B t, x 0 given. x t is an n-vector state, t a
More informationLinear System Theory (Fall 2011): Homework 1. Solutions
Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts
More informationComputational Geosciences 2 (1998) 1, 23-36
A STUDY OF THE MODELLING ERROR IN TWO OPERATOR SPLITTING ALGORITHMS FOR POROUS MEDIA FLOW K. BRUSDAL, H. K. DAHLE, K. HVISTENDAHL KARLSEN, T. MANNSETH Comptational Geosciences 2 (998), 23-36 Abstract.
More informationTIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS
International Jornal of Wavelets, Mltiresoltion and Information Processing Vol. 4, No. (26) 65 79 c World Scientific Pblishing Company TIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL
More informationarxiv: v3 [gr-qc] 29 Jun 2015
QUANTITATIVE DECAY RATES FOR DISPERSIVE SOLUTIONS TO THE EINSTEIN-SCALAR FIELD SYSTEM IN SPHERICAL SYMMETRY JONATHAN LUK AND SUNG-JIN OH arxiv:402.2984v3 [gr-qc] 29 Jn 205 Abstract. In this paper, we stdy
More informationSafe Manual Control of the Furuta Pendulum
Safe Manal Control of the Frta Pendlm Johan Åkesson, Karl Johan Åström Department of Atomatic Control, Lnd Institte of Technology (LTH) Box 8, Lnd, Sweden PSfrag {jakesson,kja}@control.lth.se replacements
More informationCurves - Foundation of Free-form Surfaces
Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections
More informationSolving large-scale PDE-constrained Bayesian inverse problems with Riemann manifold Hamiltonian Monte Carlo
Inverse Problems PAPE Solving large-scale PDE-constrained Bayesian inverse problems with iemann manifold Hamiltonian Monte Carlo To cite this article: 04 Inverse Problems 30 404 elated content - Likelihood-informed
More informationDevelopment of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model
Development of Second Order Pls Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Lemma D. Tfa*, M. Ramasamy*, Sachin C. Patwardhan **, M. Shhaimi* *Chemical Engineering Department, Universiti
More informationEVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE
13 th World Conference on Earthqake Engineering Vancover, B.C., Canada Agst 1-6, 2004 Paper No. 3099 EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE Ellen M. RATHJE 1, Wen-Jong CHANG 2, Kenneth
More informationOutline. Model Predictive Control: Current Status and Future Challenges. Separation of the control problem. Separation of the control problem
Otline Model Predictive Control: Crrent Stats and Ftre Challenges James B. Rawlings Department of Chemical and Biological Engineering University of Wisconsin Madison UCLA Control Symposim May, 6 Overview
More informationHybrid modelling and model reduction for control & optimisation
Hybrid modelling and model redction for control & optimisation based on research done by RWTH-Aachen and TU Delft presented by Johan Grievink Models for control and optimiation market and environmental
More informationElements of Coordinate System Transformations
B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and
More informationA State Space Based Implicit Integration Algorithm. for Differential Algebraic Equations of Multibody. Dynamics
A State Space Based Implicit Integration Algorithm for Differential Algebraic Eqations of Mltibody Dynamics E. J. Hag, D. Negrt, M. Ianc Janary 28, 1997 To Appear Mechanics of Strctres and Machines Abstract.
More informationA generalized Alon-Boppana bound and weak Ramanujan graphs
A generalized Alon-Boppana bond and weak Ramanjan graphs Fan Chng Department of Mathematics University of California, San Diego La Jolla, CA, U.S.A. fan@csd.ed Sbmitted: Feb 0, 206; Accepted: Jne 22, 206;
More informationBayes and Naïve Bayes Classifiers CS434
Bayes and Naïve Bayes Classifiers CS434 In this lectre 1. Review some basic probability concepts 2. Introdce a sefl probabilistic rle - Bayes rle 3. Introdce the learning algorithm based on Bayes rle (ths
More informationLab Manual for Engrd 202, Virtual Torsion Experiment. Aluminum module
Lab Manal for Engrd 202, Virtal Torsion Experiment Alminm modle Introdction In this modle, o will perform data redction and analsis for circlar cross section alminm samples. B plotting the torqe vs. twist
More informationStudy on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom
EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical
More informationWe automate the bivariate change-of-variables technique for bivariate continuous random variables with
INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,
More informationConvergence analysis of ant colony learning
Delft University of Technology Delft Center for Systems and Control Technical report 11-012 Convergence analysis of ant colony learning J van Ast R Babška and B De Schtter If yo want to cite this report
More informationRobust Tracking and Regulation Control of Uncertain Piecewise Linear Hybrid Systems
ISIS Tech. Rept. - 2003-005 Robst Tracking and Reglation Control of Uncertain Piecewise Linear Hybrid Systems Hai Lin Panos J. Antsaklis Department of Electrical Engineering, University of Notre Dame,
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationAdaptive Regularization with the Structure Tensor
1 Adaptive Reglarization with the Strctre Tensor V. Estellers, S. Soato, X. Bresson Abstract Natral images exhibit geometric strctres that are informative of the properties of the nderlying scene. Modern
More informationEffects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations
Geotechnical Safety and Risk V T. Schweckendiek et al. (Eds.) 2015 The athors and IOS Press. This article is pblished online with Open Access by IOS Press and distribted nder the terms of the Creative
More informationControl Performance Monitoring of State-Dependent Nonlinear Processes
Control Performance Monitoring of State-Dependent Nonlinear Processes Lis F. Recalde*, Hong Ye Wind Energy and Control Centre, Department of Electronic and Electrical Engineering, University of Strathclyde,
More informationExecution time certification for gradient-based optimization in model predictive control
Exection time certification for gradient-based optimization in model predictive control Giselsson, Ponts Pblished in: [Host pblication title missing] 2012 Link to pblication Citation for pblished version
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationJ.A. BURNS AND B.B. KING redced order controllers sensors/actators. The kernels of these integral representations are called fnctional gains. In [4],
Jornal of Mathematical Systems, Estimation, Control Vol. 8, No. 2, 1998, pp. 1{12 c 1998 Birkhaser-Boston A Note on the Mathematical Modelling of Damped Second Order Systems John A. Brns y Belinda B. King
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More information