Interior point algorithm for the optimization of a space shuttle re-entry trajectory
|
|
- Frank Dickerson
- 5 years ago
- Views:
Transcription
1 Interior point algorithm for the optimization of a space shuttle re-entry trajectory. Julien Laurent-Varin, CNES-INRIA-ONERA Common work with Julien Laurent-Varin J. F. Bonnans INRIA, N. Bérend ONERA, M. Haddou U. Orléans, C. Talbot CNES julien.laurent-varin@inria.fr
2 Context Ultimate aims of this work Survey Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 2
3 Context Context Ultimate aims of this work Survey A typical mission Orbital station Rendez-vous with station End of mission Booster separation Orbiter et Booster on ground Retour booster A challenge for optimal control software Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 3
4 Ultimate aims of this work Context Ultimate aims of this work Survey Efficient method for solving optimal control problems with path constraints Compute the ascent trajectory of classical launcher Compute reentry trajectory of an orbitor Compute the whole trajectory of a future space launcher (with possible rendez-vous) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 4
5 Survey Context Ultimate aims of this work Survey [1] J.T. Betts. Survey of numerical methods for trajectory optimization. AIAA J. of Guidance, Control and Dynamics, 21: , [2] J.T. Betts. Practical methods for optimal control using nonlinear programming. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 5
6 Problem descrition (P) Discretized Problem (DP) Optimal conditions of (DP) Optimal conditions of (P) Commutative diagram Summary Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 6
7 Problem descrition (P) Consider the optimal control problem (P) : Problem descrition (P) Discretized Problem (DP) Optimal conditions of (DP) Optimal conditions of (P) Commutative diagram Summary Min Φ(y(T )); ẏ(t) = f(u(t), y(t)), t [0, T ]; y(0) = y 0. (P ) Discretize dynamic with Runge-Kutta method Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 7
8 Discretized Problem (DP) Problem descrition (P) Discretized Problem (DP) Optimal conditions of (DP) Optimal conditions of (P) Commutative diagram Summary Min Φ(y N ); y k+1 = y k + h s k i=1 b if(u ki, y ki ), y ki = y k + h s k j=1 a ijf(u kj, y kj ), y 0 = y 0. (DP ) [1] W. Hager. Runge-Kutta methods in optimal control and the transformed adjoint system. Numer. Math., 87(2): , Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 8
9 Optimal conditions of (DP) Problem descrition (P) Discretized Problem (DP) Optimal conditions of (DP) Optimal conditions of (P) Commutative diagram Summary y k+1 = y k + h k s i=1 b if(u ki, y ki ), y ki = y k + h k s j=1 a ijf(u kj, y kj ), p k+1 = p k + h k s i=1 ˆb i f y (y ki, u ki ) T p ki, p ki = p k + h k s j=1 âijf y (y kj, u kj ) T p kj, 0 = f u (y k, u k ) T p k, 0 = f u (y ki, u ki ) T p ki, y 0 = y 0, p N = Φ (y N ). Where : ˆb = b and â ij = b ib j b j a ij b i for all i, j. Partitioned RK methods (DOC) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 9
10 Optimal conditions of (P) Problem descrition (P) Discretized Problem (DP) Optimal conditions of (DP) Optimal conditions of (P) Commutative diagram Summary ẏ(t) = f(u(t), y(t)), ṗ(t) = f y (u(t), y(t)) T p(t), p(t ) = Φ (y(t )), y(0) = y 0, 0 = f u (u(t), y(t)) T p(t). (OC) Previous scheme = (partitioned) symplectic schemes for (OC) deeply described in [1] [1] E. Hairer, C. Lubich, and G. Wanner. Geometric numerical integration. Springer-Verlag, Berlin, Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 10
11 Commutative diagram Problem descrition (P) Discretized Problem (DP) Optimal conditions of (DP) Optimal conditions of (P) Commutative diagram Summary optimality conditions (P ) (OC) discretization (DP ) optimality conditions discretization (DOC) (D) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 11
12 Summary Problem descrition (P) Discretized Problem (DP) Optimal conditions of (DP) Optimal conditions of (P) Commutative diagram Summary Complete error analysis for : unconstrained problems with strongly convex Hamiltonians Existing theory for : control constraint order one state constraints (Dontchev, Hager) Open problems : Singular arcs high order state constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 12
13 - I - II - III Symplectic schemes - order 4 Symplectic schemes - order 4 Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 13
14 - I - I - II - III Symplectic schemes - order 4 Symplectic schemes - order 4 For partitionned Runge-Kutta Schemes, we have order conditions based on bi-colored rooted tree. For example : Bi-coloured tree t φ(t) = 1/γ(t) s k i,j,k,l=1 j i l b i â ij a jk a jl = 1 12 Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 14
15 - II - I - II - III Symplectic schemes - order 4 Symplectic schemes - order 4 But our partitionned schemes is particular and have the following propertise : (1) ˆbi = b i â ij = b j b j b i a ji i = 1,..., s i = 1,..., s j = 1,..., s This propertise leads to computation that can be interpreted in term of graph computation. Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 15
16 - III - I - II - III Symplectic schemes - order 4 Symplectic schemes - order 4 Bi-coloured tree t φ(t) = 1/γ(t) φ(t) = φ(t 1 ) φ(t 2 ) t 1 t 2 k l j i s b i â ij a jk a jl = 1 12 s s b i b j a jk a jl b j a ji a jk a jl i,j,k,l=1 i,j,k,l=1 i,j,k,l=1 k l k l j j i i Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 16
17 Symplectic schemes - order 4 - I - II - III Symplectic schemes - order 4 Symplectic schemes - order 4 Graph Condition Graph Condition 1 b k a lk d k d l = 1 8 bi b k a ik c i d k = 5 24 c 2 j d j = 1 12 ajk d j c k = 1 24 bi a ij c i c j = 1 8 bi c 3 i = b k c k d 2 k = b 2 l d 3 l = 1 4 Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 17
18 Symplectic schemes - order 4 In the previous table, we use the usual notations d j = i b i a ij - I - II - III Symplectic schemes - order 4 Symplectic schemes - order 4 and c i = j a ij. Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 18
19 History Basic idea For our problem For our problem Combination with mesh refinement Sparsity structure Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 19
20 History History Basic idea For our problem For our problem Combination with mesh refinement Sparsity structure [1] A.V. Fiacco and G.P. McCormick. Nonlinear programming. John Wiley and Sons, Inc., New York-London-Sydney, [2] L. G. Hačijan. A polynomial algorithm in linear programming. Dokl. Akad. Nauk SSSR, 244(5): , [3] N. Karmarkar. A new polynomial-time algorithm for linear programming. Combinatorica, 4(4): , Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 20
21 Basic idea Transformation of constrained problem (forget equality constraints in order to simplify presentation) History Basic idea For our problem For our problem Combination with mesh refinement Sparsity structure Logarithmic penalty Min f(x); g(x) 0 (P ) Min f(x) ε ln g(x) (P ε ) Unconstrained problem solved by Newton steps on optimality conditions, with basic safeguards Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 21
22 For our problem Min T 0 l(u(t), y(t), u p, t)dt + l f (y(t ), u p, T ); History Basic idea For our problem ẏ(t) = f(u(t), y(t), u p, t), t [0, T ]; For our problem Combination with mesh refinement Sparsity structure y i (0) = y 0 i i I 0, y i (T ) = y T i i I T a g(u(t), y(t), u p, t) b, t [0, T ]. (P ) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 22
23 For our problem History Basic idea For our problem For our problem Combination with mesh refinement Sparsity structure With Min T 0 l ε (u(t), y(t), u p, t)dt + l f (y(t ), u p, T ); ẏ(t) = f(u(t), y(t), u p, t), t [0, T ]; y i (0) = yi 0 i I 0, y i (T ) = yi T i I T. (P ε ) l ε (u, y, u p, t) := l(u, y, u p, t) ε X i [log(g i (u, y, u p, t) a i ) + log(b i g i (u, y, u p, t)) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 23
24 Combination with mesh refinement We must have ε 0 (IP parameter) and E 0 (Integration error) History Basic idea For our problem For our problem Combination with mesh refinement Sparsity structure Good News : Both can be reduced simultaneously : interpolate values at additional points. Our (somewhat arbitrary) choice : For a fixed ε, we refine mesh until E cε; c = 1/10 in our tests. Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 24
25 Sparsity structure History Basic idea For our problem For our problem Combination with mesh refinement Sparsity structure A C B D Band matrix: distributed variables, band QR available Full pieces: static parameters Solve by elimination of distributed variables Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 25
26 Modelisation of error Optimal refinement problem Maximal gain index Resolution of (ORP) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 26
27 Modelisation of error Modelisation of error Optimal refinement problem Maximal gain index Resolution of (ORP) k : index of time step k : error model (valid for h k small) estimated by variable order (symplectic) schemes. Adding (q k 1) points in k th step estimate of error becomes : e k q p k e k = C k h p+1 Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 27
28 Optimal refinement problem Modelisation of error Optimal refinement problem Maximal gain index Resolution of (ORP) Reach a specific error estimation Minimise the number of added points Integer programming problem with a single nonlinear constraints Min q N N N q k ; k=1 N k=1 e k q p k E (ORP ) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 28
29 Maximal gain index Modelisation of error Optimal refinement problem Maximal gain index Resolution of (ORP) Definition 1. Maximal marginal gain g and maximal gain index k g for which the maximum error reduction is obtained by adding only one point : (2) (3) g(q) := max k k g := argmax k e k ( 1/q p k 1/(q k + 1) p) e k ( 1/q p k 1/(q k + 1) p) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 29
30 Resolution of (ORP) Modelisation of error Optimal refinement problem Maximal gain index Resolution of (ORP) Algorithm 1. AORP For k = 1,..., N do q k := 1. End for While N k=1 e k/q p k > E do Compute k g, the maximal gain index. q kg := q kg + 1. End While This algorithm solves (ORP) in O(((E 0 /E) p + 1)N log N) operations Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 30
31 Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 31
32 Control variables µ Control variables α State variables Dynamic I Dynamic II γ Cost Constraints Result without heating constraints Result with heating constraints Result (state) Bank angle µ V Flightpath angle (γ) Angle of attack (α) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 32
33 State variables Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints z altitude V velocity λ longitude γ flightpath angle φ latitude ψ azimuth Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 33
34 Dynamic I ż = V sin γ Control variables State variables Dynamic I λ = V cos γ sin ψ (z + R e ) cos φ Dynamic II Cost Constraints Result without heating constraints φ = V cos γ cos ψ z + R e Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 34
35 Dynamic II V = D m g sin γ Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints γ = L cos µ mv + ψ = L sin µ mv cos γ + ( V 2 ) cos γ g z + R e V V z + R e cos γ sin ψ tan φ Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 35
36 Cost Max crossrange (final latitude) J := φ(t ) Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 36
37 Constraints Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Control constraints : { 0 α µ 0 Mixed control-state contraint (Heating flux) : with q(α) ρv 3.07 Q max q(α) = c 0 + c 1 α + c 2 α 2 + c 3 α 3 Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 37
38 Result without heating constraints Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints ε Integration error Time steps Newton steps Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 38
39 Result (control) Angle of attack (deg) Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Bank angle (deg) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 39
40 Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 40
41 Result with heating constraints We add the constraint of heating and we compare the results Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 41
42 Result (control) Angle of attack (deg) Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Bank angle (deg) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 42
43 Result (state) Altitude (ft) Longitude (deg) Latitude (deg) Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Velocity (ft/s) Flight path angle (deg) Azimute (deg) Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 43
44 Result (constraint:heating) 180 Heating 160 Control variables State variables Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 44
45 3D view: no heating constraints Control variables State variables 26 Dynamic I Dynamic II Cost Constraints Result without heating constraints Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints h/10^ theta phi Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 45
46 3D view: heating constraints Control variables State variables Dynamic I 26 Dynamic II Cost Constraints Result without heating constraints h/10^ Result with heating constraints Result (state) Result (constraint:heating) 3D view: no heating constraints 3D view: heating constraints theta phi 40 Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 46
47 Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 47
48 Scenario graph Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references V set of vertices E set of edges Edges: Distributed in time cost and constraints Vertices: Junction conditions, cost and constraints Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 48
49 Edges: variables Edge E e = (i, j) y(t) IR n e : state u(t) IR m e : control π IR r e : static optimization parameters Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 49
50 A view after discretization Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 50
51 Edges: cost and constraints l e (t, y(t), u(t), π): Distributed cost ẏ(t) = f e (t, y(t), u(t), π): Dynamics g e (t, y(t), u(t), π) 0: Distributed constraints Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references The value function over arc = integral of cost function. Feasibility = dynamics are satisfied and distributed constraints hold. Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 51
52 Vertices: variables Vertice j V Variables z j IR n j Includes copy of boundary (initial of final) state y B e and parameters π e of connected edges Possibility of other variables Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 52
53 Vertices: cost and constraints l i (z): Cost y B e = z e : Copy constraints g i (z)) = 0: Equality constraints g i (z)) 0: Inequality constraints Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 53
54 Cost of factorization and solve Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references Arcs: O(n 2 N T ) - n number of state variables - N T average number of time steps per arc Vertices (junctions): O(n) size linear systems - Recursive elimination, starting from leaves - O(n 3 ) operations for each (elimination of) vertex Total number of operations: n 2 O( E N T + n V ) Seems to be the least possible number! Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 54
55 Perspectives Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references Implementation of multi-arc methodology Test on various examples Analysis of convergence (difficult) - Vanishing parameter ε - Mesh refinement Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 55
56 Some references Scenario graph Edges: variables A view after discretization Edges: cost and constraints Vertices: variables Vertices: cost and constraints Cost of factorization and solve Perspectives Some references [1] J. Laurent-Varin, N. Bérend, F. Bonnans, M. Haddou, and C. Talbot. On the refinement of discretization for optimal control problems. 16 th IFAC SYMPOSIUM Automatic Control in Aerospace, june, St. Petersburg, Russia. [2] J Laurent-Varin and F Bonnans. Computation of order conditions for symplectic partitioned Runge-Kutta schemes with application to optimal control. Technical Report RR 5398, INRIA, Julien Laurent-Varin 7 th French-Latin American Congress on Applied Mathematics - 56
57 56-1
Numerical Resolution of optimal control problems by a Piecewise Linear continuation method
Numerical Resolution of optimal control problems by a Piecewise Linear continuation method Pierre Martinon november 25 Shooting method and homotopy Introduction We would like to solve numerically the
More informationIntroduction to Optimization Problems and Methods
Introduction to Optimization Problems and Methods wjch@umich.edu December 10, 2009 Outline 1 Linear Optimization Problem Simplex Method 2 3 Cutting Plane Method 4 Discrete Dynamic Programming Problem Simplex
More informationAlgorithms for convex optimization
Algorithms for convex optimization Michal Kočvara Institute of Information Theory and Automation Academy of Sciences of the Czech Republic and Czech Technical University kocvara@utia.cas.cz http://www.utia.cas.cz/kocvara
More informationMathematics in Orbit
Mathematics in Orbit Dan Kalman American University Slides and refs at www.dankalman.net Outline Basics: 3D geospacial models Keyhole Problem: Related Rates! GPS: space-time triangulation Sensor Diagnosis:
More informationThe Design and Implementation of a Modeling Package
The Design and Implementation of a Modeling Package p.1/20 The Design and Implementation of a Modeling Package Hossein ZivariPiran hzp@cs.toronto.edu Department of Computer Science University of Toronto
More informationResearch Interests Optimization:
Mitchell: Research interests 1 Research Interests Optimization: looking for the best solution from among a number of candidates. Prototypical optimization problem: min f(x) subject to g(x) 0 x X IR n Here,
More informationMATH3016: OPTIMIZATION
MATH3016: OPTIMIZATION Lecturer: Dr Huifu Xu School of Mathematics University of Southampton Highfield SO17 1BJ Southampton Email: h.xu@soton.ac.uk 1 Introduction What is optimization? Optimization is
More informationRapid Design Space Exploration for Conceptual Design of Hypersonic Missions
Rapid Design Space Exploration for Conceptual Design of Hypersonic Missions Michael J. Grant, Ian G. Clark, and Robert D. Braun Georgia Institute of Technology, Atlanta, GA, 3332 During conceptual design,
More informationPanagiotis Tsiotras. Dynamics and Control Systems Laboratory Daniel Guggenheim School of Aerospace Engineering Georgia Institute of Technology
Panagiotis Tsiotras Dynamics and Control Systems Laboratory Daniel Guggenheim School of Aerospace Engineering Georgia Institute of Technology ICRAT 12 Tutotial on Methods for Optimal Trajectory Design
More informationConvex Optimization M2
Convex Optimization M2 Lecture 1 A. d Aspremont. Convex Optimization M2. 1/49 Today Convex optimization: introduction Course organization and other gory details... Convex sets, basic definitions. A. d
More informationA fast algorithm for sparse reconstruction based on shrinkage, subspace optimization and continuation [Wen,Yin,Goldfarb,Zhang 2009]
A fast algorithm for sparse reconstruction based on shrinkage, subspace optimization and continuation [Wen,Yin,Goldfarb,Zhang 2009] Yongjia Song University of Wisconsin-Madison April 22, 2010 Yongjia Song
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-4: Constrained optimization Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428 June
More informationValue function and optimal trajectories for a control problem with supremum cost function and state constraints
Value function and optimal trajectories for a control problem with supremum cost function and state constraints Hasnaa Zidani ENSTA ParisTech, University of Paris-Saclay joint work with: A. Assellaou,
More informationPRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING. 1. Introduction
PRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING KELLER VANDEBOGERT AND CHARLES LANNING 1. Introduction Interior point methods are, put simply, a technique of optimization where, given a problem
More informationLecture Note 3: Rotational Motion
ECE5463: Introduction to Robotics Lecture Note 3: Rotational Motion Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio, USA Spring 2018 Lecture 3 (ECE5463
More informationarxiv: v2 [math.oc] 12 Feb 2019
A CENTRAL PATH FOLLOWING METHOD USING THE NORMALIZED GRADIENTS Long Chen 1, Wenyi Chen 2, and Kai-Uwe Bletzinger 1 arxiv:1902.04040v2 [math.oc] 12 Feb 2019 1 Chair of Structural Analysis, Technical University
More informationOn the graphical display of Powell-Sabin splines: a comparison of three piecewise linear approximations
On the graphical display of Powell-Sabin splines: a comparison of three piecewise linear approximations Hendrik Speleers Paul Dierckx Stefan Vandewalle Report TW515, January 008 Ò Katholieke Universiteit
More informationA Simplex-Cosine Method for Solving Hard Linear Problems
A Simplex-Cosine Method for Solving Hard Linear Problems FEDERICO TRIGOS 1, JUAN FRAUSTO-SOLIS 2 and RAFAEL RIVERA-LOPEZ 3 1 Division of Engineering and Sciences ITESM, Campus Toluca Av. Eduardo Monroy
More informationLecture 9. Introduction to Numerical Techniques
Lecture 9. Introduction to Numerical Techniques Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering May 27, 2015 1 / 25 Logistics hw8 (last one) due today. do an easy problem
More informationScalable Algorithms in Optimization: Computational Experiments
Scalable Algorithms in Optimization: Computational Experiments Steven J. Benson, Lois McInnes, Jorge J. Moré, and Jason Sarich Mathematics and Computer Science Division, Argonne National Laboratory, Argonne,
More informationComputational Methods. Constrained Optimization
Computational Methods Constrained Optimization Manfred Huber 2010 1 Constrained Optimization Unconstrained Optimization finds a minimum of a function under the assumption that the parameters can take on
More informationMeasuring Lengths The First Fundamental Form
Differential Geometry Lia Vas Measuring Lengths The First Fundamental Form Patching up the Coordinate Patches. Recall that a proper coordinate patch of a surface is given by parametric equations x = (x(u,
More informationARTICLE IN PRESS. Topology and its Applications ( )
S0166-8641(05)0018-7/FLA AID:2822 Vol. ( ) [DTD5] P.1 (1-13) TOPOL:m1a v 1.42 Prn:11/08/2005; 10:03 top2822 by:rita p. 1 Topology and its Applications ( ) www.elsevier.com/locate/topol The topological
More informationEXERCISES SHORTEST PATHS: APPLICATIONS, OPTIMIZATION, VARIATIONS, AND SOLVING THE CONSTRAINED SHORTEST PATH PROBLEM. 1 Applications and Modelling
SHORTEST PATHS: APPLICATIONS, OPTIMIZATION, VARIATIONS, AND SOLVING THE CONSTRAINED SHORTEST PATH PROBLEM EXERCISES Prepared by Natashia Boland 1 and Irina Dumitrescu 2 1 Applications and Modelling 1.1
More information3D Kinematics. Consists of two parts
D Kinematics Consists of two parts D rotation D translation The same as D D rotation is more complicated than D rotation (restricted to z-ais) Net, we will discuss the treatment for spatial (D) rotation
More informationPerformance Evaluation of an Interior Point Filter Line Search Method for Constrained Optimization
6th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING, Venice, Italy, November 21-23, 2007 18 Performance Evaluation of an Interior Point Filter Line Search Method for Constrained
More informationComputer Algebra Investigation of Known Primitive Triangle-Free Strongly Regular Graphs
Computer Algebra Investigation of Known Primitive Triangle-Free Strongly Regular Graphs Matan Ziv-Av (Jointly with Mikhail Klin) Ben-Gurion University of the Negev SCSS 2013 RISC, JKU July 5, 2013 Ziv-Av
More informationDavid G. Luenberger Yinyu Ye. Linear and Nonlinear. Programming. Fourth Edition. ö Springer
David G. Luenberger Yinyu Ye Linear and Nonlinear Programming Fourth Edition ö Springer Contents 1 Introduction 1 1.1 Optimization 1 1.2 Types of Problems 2 1.3 Size of Problems 5 1.4 Iterative Algorithms
More informationLeaderless Formation Control for Multiple Autonomous Vehicles. Wei Ren
AIAA Guidance, Navigation, and Control Conference and Exhibit - 4 August 6, Keystone, Colorado AIAA 6-669 Leaderless Formation Control for Multiple Autonomous Vehicles Wei Ren Department of Electrical
More informationOn the deviation of a parametric cubic spline interpolant from its data polygon
Computer Aided Geometric Design 25 (2008) 148 156 wwwelseviercom/locate/cagd On the deviation of a parametric cubic spline interpolant from its data polygon Michael S Floater Department of Computer Science,
More informationSurrogate Gradient Algorithm for Lagrangian Relaxation 1,2
Surrogate Gradient Algorithm for Lagrangian Relaxation 1,2 X. Zhao 3, P. B. Luh 4, and J. Wang 5 Communicated by W.B. Gong and D. D. Yao 1 This paper is dedicated to Professor Yu-Chi Ho for his 65th birthday.
More informationCourse Introduction. Scheduling: Terminology and Classification
Outline DM87 SCHEDULING, TIMETABLING AND ROUTING Lecture 1 Course Introduction. Scheduling: Terminology and Classification 1. Course Introduction 2. Scheduling Problem Classification Marco Chiarandini
More informationBilinear Programming
Bilinear Programming Artyom G. Nahapetyan Center for Applied Optimization Industrial and Systems Engineering Department University of Florida Gainesville, Florida 32611-6595 Email address: artyom@ufl.edu
More informationThe Cross-Entropy Method for Mathematical Programming
The Cross-Entropy Method for Mathematical Programming Dirk P. Kroese Reuven Y. Rubinstein Department of Mathematics, The University of Queensland, Australia Faculty of Industrial Engineering and Management,
More informationCONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION USED AS A TEST CASE FOR THE INVISCID PART OF RANS SOLVERS
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION
More informationNumerical Method in Optimization as a Multi-stage Decision Control System
Numerical Method in Optimization as a Multi-stage Decision Control System B.S. GOH Institute of Mathematical Sciences University of Malaya 50603 Kuala Lumpur MLYSI gohoptimum@gmail.com bstract: - Numerical
More informationCompactness Theorems for Saddle Surfaces in Metric Spaces of Bounded Curvature. Dimitrios E. Kalikakis
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 51, 2005 (45 52) Compactness Theorems for Saddle Surfaces in Metric Spaces of Bounded Curvature Dimitrios E. Kalikakis Abstract The notion of a non-regular
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationROC-HJ: Reachability analysis and Optimal Control problems - Hamilton-Jacobi equations
ROC-HJ: Reachability analysis and Optimal Control problems - Hamilton-Jacobi equations Olivier Bokanowski, Anna Désilles, Hasnaa Zidani February 2013 1 Introduction The software ROC-HJ is a C++ library
More informationIntroduction to Optimization
Introduction to Optimization Second Order Optimization Methods Marc Toussaint U Stuttgart Planned Outline Gradient-based optimization (1st order methods) plain grad., steepest descent, conjugate grad.,
More informationStructure in Optimization: Factorable Programming and Functions
Structure in Optimization: Factorable Programming and Functions Laurent Hascoët, Shahadat Hossain and Trond Steihaug Abstract It is frequently observed that effective exploitation of problem structure
More information3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis
Adv. Theor. Appl. Mech., Vol. 1, 2008, no. 3, 131-139 3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis Yunhua Luo Department of Mechanical & Manufacturing Engineering, University of
More informationTMA946/MAN280 APPLIED OPTIMIZATION. Exam instructions
Chalmers/GU Mathematics EXAM TMA946/MAN280 APPLIED OPTIMIZATION Date: 03 05 28 Time: House V, morning Aids: Text memory-less calculator Number of questions: 7; passed on one question requires 2 points
More informationNATCOR Convex Optimization Linear Programming 1
NATCOR Convex Optimization Linear Programming 1 Julian Hall School of Mathematics University of Edinburgh jajhall@ed.ac.uk 5 June 2018 What is linear programming (LP)? The most important model used in
More informationRandomized rounding of semidefinite programs and primal-dual method for integer linear programming. Reza Moosavi Dr. Saeedeh Parsaeefard Dec.
Randomized rounding of semidefinite programs and primal-dual method for integer linear programming Dr. Saeedeh Parsaeefard 1 2 3 4 Semidefinite Programming () 1 Integer Programming integer programming
More informationarxiv: v1 [math.co] 24 Aug 2009
SMOOTH FANO POLYTOPES ARISING FROM FINITE PARTIALLY ORDERED SETS arxiv:0908.3404v1 [math.co] 24 Aug 2009 TAKAYUKI HIBI AND AKIHIRO HIGASHITANI Abstract. Gorenstein Fano polytopes arising from finite partially
More informationAn augmented Lagrangian method for equality constrained optimization with fast infeasibility detection
An augmented Lagrangian method for equality constrained optimization with fast infeasibility detection Paul Armand 1 Ngoc Nguyen Tran 2 Institut de Recherche XLIM Université de Limoges Journées annuelles
More informationAdaptive and Smooth Surface Construction by Triangular A-Patches
Adaptive and Smooth Surface Construction by Triangular A-Patches Guoliang Xu Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences, Beijing, China Abstract
More informationPOLYHEDRAL GEOMETRY. Convex functions and sets. Mathematical Programming Niels Lauritzen Recall that a subset C R n is convex if
POLYHEDRAL GEOMETRY Mathematical Programming Niels Lauritzen 7.9.2007 Convex functions and sets Recall that a subset C R n is convex if {λx + (1 λ)y 0 λ 1} C for every x, y C and 0 λ 1. A function f :
More informationInstitutionen för matematik, KTH.
Institutionen för matematik, KTH. Chapter 10 projective toric varieties and polytopes: definitions 10.1 Introduction Tori varieties are algebraic varieties related to the study of sparse polynomials.
More informationInterpolation by Spline Functions
Interpolation by Spline Functions Com S 477/577 Sep 0 007 High-degree polynomials tend to have large oscillations which are not the characteristics of the original data. To yield smooth interpolating curves
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationLecture 9. Semidefinite programming is linear programming where variables are entries in a positive semidefinite matrix.
CSE525: Randomized Algorithms and Probabilistic Analysis Lecture 9 Lecturer: Anna Karlin Scribe: Sonya Alexandrova and Keith Jia 1 Introduction to semidefinite programming Semidefinite programming is linear
More informationCooperative Task Planning of Multi-Robot Systems with Temporal Constraints 1
Paper Number: ICRA-2003 Cooperative Task Planning of Multi-Robot Systems with Temporal Constraints 1 Feng-Li Lian a and Richard Murray b (a) Electrical Engineering, National Taiwan University (b) Control
More informationINTRODUCTION TO CLUSTER ALGEBRAS
INTRODUCTION TO CLUSTER ALGEBRAS NAN LI (MIT) Cluster algebras are a class of commutative ring, introduced in 000 by Fomin and Zelevinsky, originally to study Lusztig s dual canonical basis and total positivity.
More informationSPIRAL TRANSITION CURVES AND THEIR APPLICATIONS. Zulfiqar Habib and Manabu Sakai. Received August 21, 2003
Scientiae Mathematicae Japonicae Online, e-2004, 25 262 25 SPIRAL TRANSITION CURVES AND THEIR APPLICATIONS Zulfiqar Habib and Manabu Sakai Received August 2, 200 Abstract. A method for family of G 2 planar
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More informationDecomposition of log-linear models
Graphical Models, Lecture 5, Michaelmas Term 2009 October 27, 2009 Generating class Dependence graph of log-linear model Conformal graphical models Factor graphs A density f factorizes w.r.t. A if there
More informationComparison of Interior Point Filter Line Search Strategies for Constrained Optimization by Performance Profiles
INTERNATIONAL JOURNAL OF MATHEMATICS MODELS AND METHODS IN APPLIED SCIENCES Comparison of Interior Point Filter Line Search Strategies for Constrained Optimization by Performance Profiles M. Fernanda P.
More informationWhat is linear programming (LP)? NATCOR Convex Optimization Linear Programming 1. Solving LP problems: The standard simplex method
NATCOR Convex Optimization Linear Programming 1 Julian Hall School of Mathematics University of Edinburgh jajhall@ed.ac.uk 14 June 2016 What is linear programming (LP)? The most important model used in
More informationOn Hierarchical Combination of Discrete Wavelet Transform with Canny Edge Detector
On Hierarchical Combination of Discrete Wavelet Transform with Canny Edge Detector Tiehang Duan Department of Computer Science and Engineering The State University of New York at Buffalo Abstract Discrete
More informationQuaternions and Rotations
CSCI 520 Computer Animation and Simulation Quaternions and Rotations Jernej Barbic University of Southern California 1 Rotations Very important in computer animation and robotics Joint angles, rigid body
More informationA ray-based IPDG method for high-frequency time-domain acoustic wave propagation in inhomogeneous media
A ray-based IPDG method for high-frequency time-domain acoustic wave propagation in inhomogeneous media Eric T. Chung Chi Yeung Lam Jianliang Qian arxiv:1704.06916v1 [math.na] 23 Apr 2017 October 15, 2018
More informationRobust time-varying shortest path with arbitrary waiting time at vertices
Croatian Operational Research Review 525 CRORR 8(2017), 525 56 Robust time-varying shortest path with arbitrary waiting time at vertices Gholamhassan Shirdel 1, and Hassan Rezapour 1 1 Department of Mathematics,
More information2008 by authors and 2008 American Institute of Aeronautics and Astronautics (AIAA)
Karelahti, J., Virtanen, K., and Öström, J., Automated Generation of Realistic Near Optimal Aircraft Trajectories, Journal of Guidance, Control, and Dynamics, accepted for publication. 28 by authors and
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn School of Mathematical Sciences Tel Aviv University Michael S. Floater Department of Informatics University of
More informationNumerical schemes for Hamilton-Jacobi equations, control problems and games
Numerical schemes for Hamilton-Jacobi equations, control problems and games M. Falcone H. Zidani SADCO Spring School Applied and Numerical Optimal Control April 23-27, 2012, Paris Lecture 2/3 M. Falcone
More informationConvexity Theory and Gradient Methods
Convexity Theory and Gradient Methods Angelia Nedić angelia@illinois.edu ISE Department and Coordinated Science Laboratory University of Illinois at Urbana-Champaign Outline Convex Functions Optimality
More informationarxiv: v1 [cs.ne] 22 Mar 2016
Adaptive Parameter Selection in Evolutionary Algorithms by Reinforcement Learning with Dynamic Discretization of Parameter Range arxiv:1603.06788v1 [cs.ne] 22 Mar 2016 ABSTRACT Arkady Rost ITMO University
More informationAlgorithm That Mimics Human Perceptual Grouping of Dot Patterns
Algorithm That Mimics Human Perceptual Grouping of Dot Patterns G. Papari and N. Petkov Institute of Mathematics and Computing Science, University of Groningen, P.O.Box 800, 9700 AV Groningen, The Netherlands
More informationComputation of Voronoi Diagrams and Delaunay Triangulation via Parametric Linear Programming
Computation of Voronoi Diagrams and Delaunay Triangulation via Parametric Linear Programming Saša V. Raković, Pascal Grieder and Colin Jones September 14, 2004 Abstract This note illustrates how Voronoi
More informationCS612 - Algorithms in Bioinformatics
Fall 2017 Structural Manipulation November 22, 2017 Rapid Structural Analysis Methods Emergence of large structural databases which do not allow manual (visual) analysis and require efficient 3-D search
More informationTripod Configurations
Tripod Configurations Eric Chen, Nick Lourie, Nakul Luthra Summer@ICERM 2013 August 8, 2013 Eric Chen, Nick Lourie, Nakul Luthra (S@I) Tripod Configurations August 8, 2013 1 / 33 Overview 1 Introduction
More informationApplied Lagrange Duality for Constrained Optimization
Applied Lagrange Duality for Constrained Optimization Robert M. Freund February 10, 2004 c 2004 Massachusetts Institute of Technology. 1 1 Overview The Practical Importance of Duality Review of Convexity
More informationAnimation Lecture 10 Slide Fall 2003
Animation Lecture 10 Slide 1 6.837 Fall 2003 Conventional Animation Draw each frame of the animation great control tedious Reduce burden with cel animation layer keyframe inbetween cel panoramas (Disney
More informationOptimal boundary control of a tracking problem for a parabolic distributed system using hierarchical fuzzy control and evolutionary algorithms
Optimal boundary control of a tracking problem for a parabolic distributed system using hierarchical fuzzy control and evolutionary algorithms R.J. Stonier, M.J. Drumm and J. Bell Faculty of Informatics
More informationLEAST COST ROUTING ALGORITHM WITH THE STATE SPACE RELAXATION IN A CENTRALIZED NETWORK
VOL., NO., JUNE 08 ISSN 896608 00608 Asian Research Publishing Network (ARPN). All rights reserved. LEAST COST ROUTING ALGORITHM WITH THE STATE SPACE RELAXATION IN A CENTRALIZED NETWORK Y. J. Lee Department
More informationQuaternion to Euler Angle Conversion for Arbitrary Rotation Sequence Using Geometric Methods
uaternion to Euler Angle Conversion for Arbitrary Rotation Sequence Using Geometric Methods ê = normalized Euler ation axis i Noel H. Hughes Nomenclature = indices of first, second and third Euler
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationGeneralized barycentric coordinates
Generalized barycentric coordinates Michael S. Floater August 20, 2012 In this lecture, we review the definitions and properties of barycentric coordinates on triangles, and study generalizations to convex,
More informationState Plane Coordinates and Computations using them GISC Spring 2013
State Plane Coordinates and Computations using them GISC-3325 - Spring 2013 Map Projections From UNAVCO site hosting.soonet.ca/eliris/gpsgis/lec2geodesy.html Taken from Ghilani, SPC State Plane Coordinate
More informationDIRECT SEQUENTIAL DYNAMIC OPTIMIZATION WITH AUTOMATIC SWITCHING STRUCTURE DETECTION. Martin Schlegel, Wolfgang Marquardt 1
DIRECT SEQUENTIAL DYNAMIC OPTIMIZATION WITH AUTOMATIC SWITCHING STRUCTURE DETECTION Martin Schlegel, Wolfgang Marquardt 1 Lehrstuhl für Prozesstechnik, RWTH Aachen University D 52056 Aachen, Germany Abstract:
More informationInertial Measurement Units II!
! Inertial Measurement Units II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 10! stanford.edu/class/ee267/!! wikipedia! Polynesian Migration! Lecture Overview! short review of
More informationMA4254: Discrete Optimization. Defeng Sun. Department of Mathematics National University of Singapore Office: S Telephone:
MA4254: Discrete Optimization Defeng Sun Department of Mathematics National University of Singapore Office: S14-04-25 Telephone: 6516 3343 Aims/Objectives: Discrete optimization deals with problems of
More informationNonlinear Programming
Nonlinear Programming SECOND EDITION Dimitri P. Bertsekas Massachusetts Institute of Technology WWW site for book Information and Orders http://world.std.com/~athenasc/index.html Athena Scientific, Belmont,
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationOPTIMUM FLIGHT TRAJECTORIES AND SENSITIVITY ANALYSIS FOR TERRAIN COLLISION AVOIDANCE SYSTEMS
2 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES OPTIMUM FLIGHT TRAJECTORIES AND SENSITIITY ANALYSIS FOR TERRAIN COLLISION AOIDANCE SYSTEMS T Sharma*, C Bil*, A Eberhard** * Sir Lawrence Wackett
More information43rd International Mathematical Olympiad
43rd International Mathematical Olympiad Glasgow, United Kingdom, July 2002 1 Let n be a positive integer let T be the set of points (x, y) in the plane where x and y are non-negative integers and x +
More informationA Feasible Region Contraction Algorithm (Frca) for Solving Linear Programming Problems
A Feasible Region Contraction Algorithm (Frca) for Solving Linear Programming Problems E. O. Effanga Department of Mathematics/Statistics and Comp. Science University of Calabar P.M.B. 1115, Calabar, Cross
More informationSensor-Target and Weapon-Target Pairings Based on Auction Algorithm
Proceedings of the 11th WSEAS International Conference on APPLIED MATHEMATICS, Dallas, Texas, USA, March 22-24, 2007 92 Sensor-Target and Weapon-Target Pairings Based on Auction Algorithm Z. R. BOGDANOWICZ,
More informationCOLUMN GENERATION IN LINEAR PROGRAMMING
COLUMN GENERATION IN LINEAR PROGRAMMING EXAMPLE: THE CUTTING STOCK PROBLEM A certain material (e.g. lumber) is stocked in lengths of 9, 4, and 6 feet, with respective costs of $5, $9, and $. An order for
More informationREINER HORST. License or copyright restrictions may apply to redistribution; see
MATHEMATICS OF COMPUTATION Volume 66, Number 218, April 1997, Pages 691 698 S 0025-5718(97)00809-0 ON GENERALIZED BISECTION OF n SIMPLICES REINER HORST Abstract. A generalized procedure of bisection of
More informationNelson Functions 11 Errata
Nelson Functions 11 Errata 1: Introduction to Functions Location Question Correct Answer Getting Started 6a Graph is correct but vertex and axis of symmetry are not labelled. Add blue point labelled (in
More informationVisual Recognition: Image Formation
Visual Recognition: Image Formation Raquel Urtasun TTI Chicago Jan 5, 2012 Raquel Urtasun (TTI-C) Visual Recognition Jan 5, 2012 1 / 61 Today s lecture... Fundamentals of image formation You should know
More informationAn Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm
An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm P.A. Sherar, C.P. Thompson, B. Xu, B. Zhong Abstract A new method is presented to deal with shape optimization problems.
More informationDetermination of Angle of Attack (AOA) for Rotating Blades
Downloaded from orbit.dtu.dk on: Sep 1, 218 Determination of Angle of Attack (AOA) for Rotating Blades Shen, Wen Zhong; Hansen, Martin Otto Laver; Sørensen, Jens Nørkær Published in: Wind Energy Publication
More informationEffect of Uncertainties on UCAV Trajectory Optimisation Using Evolutionary Programming
2007 Information, Decision and Control Effect of Uncertainties on UCAV Trajectory Optimisation Using Evolutionary Programming Istas F Nusyirwan 1, Cees Bil 2 The Sir Lawrence Wackett Centre for Aerospace
More informationConvex Optimization CMU-10725
Convex Optimization CMU-10725 Ellipsoid Methods Barnabás Póczos & Ryan Tibshirani Outline Linear programs Simplex algorithm Running time: Polynomial or Exponential? Cutting planes & Ellipsoid methods for
More informationA Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete Data
Applied Mathematical Sciences, Vol. 1, 16, no. 7, 331-343 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/1.1988/ams.16.5177 A Cumulative Averaging Method for Piecewise Polynomial Approximation to Discrete
More informationLATIN SQUARES AND THEIR APPLICATION TO THE FEASIBLE SET FOR ASSIGNMENT PROBLEMS
LATIN SQUARES AND THEIR APPLICATION TO THE FEASIBLE SET FOR ASSIGNMENT PROBLEMS TIMOTHY L. VIS Abstract. A significant problem in finite optimization is the assignment problem. In essence, the assignment
More information