For later use we dene a constraint-set Z k utilizing an expression for the state predictor Z k (x k )=fz k =( k k ): k 2 ^ k 2 ^ k = F k ( k x k )g 8
|
|
- Colleen Shields
- 6 years ago
- Views:
Transcription
1 IDENTIFICATION AND CONVEITY IN OPTIMIZING CONTROL Bjarne A. Foss Department of Engineering Cybernetics, Norwegian University of Science and Technology, N-734 Trondheim, Norway. Tor A. Johansen SINTEF Electronics and Cybernetics, Automatic Control Department, N-734 Trondheim, Norway. Abstract. The interplay between identication of a predictor and the convexity of an associated optimization-based control problem is studied. The reason for the appearance of nonconvexity is investigated and a method to favour convexity of the optimization-based control problem in the identication problem is suggested. Further, a parametrization which makes the identication problem tractable is proposed. Simulation results suggest a merit to method. Keywords. Optimization, nonconvex problems, identication, nonlinear systems, model predictive control 1. INTRODUCTION The combined use of dynamic models and optimization oers a concept in which process knowledge can be linked to operational goals formulated by some optimization criterion. This concept has seen widespread use. One such area is process control and the use of model predictive control (MPC). MPC refers to a class of algorithms where an optimization problem is solved repetitively, at every new time-instant, see for example Lee (1996). The interaction between control and optimization is discussed in an illuminating way in Mayne (199). In this paper the divide between convex and nonconvex problems is emphasized. In particular problems arise by using a nonlinear predictor as opposed to a linear predictor since the optimization problem in the former cases in general becomes nonconvex. This latter observation forms the departure point for this study. The paper is structured as follows: First, a problem denition is provided. Thereafter, we investigate convexity in nonlinear optimization-based control. In particular, we discuss design tradeos between identifying a predictor, and convexity of an associated optimization problem. Section 4 highlights some of the ndings by a simulation example, while a discussion and some conclusions ends the paper. 2. PROBLEM DEFINITION The optimality criterion is dened in discrete-time k (z k )= k ( k k )= i2i k l i (x i+1 u i ) (1) where k = fu T k ::: ut k+n;1 gt =U ::: U k = fx T k+1 ::: x T N g T = ::: I k = fk : : : k + N ; 1g u i 2 U R mu x i 2 R mx u i denotes the control input of a system. Since we are assuming a time-discrete model the control input u i is constant during the time span [i i + 1 >. x i denotes the system state at time i. Above, the sets and (or U and ) invoke constraints on k and k. The states and control inputs are, however, not independent, since the state sequence k in a dynamic system is uniquely dened by the initial value, ie. x k, and the control input sequence k.we dene a state-space model, valid on at least I k, and a predictor by h i (x i+1 x i u i )=x i+1 ; f i (x i u i )= (2) k = F k ( k x k ) (3) = ff k1 ( k x k ) T ::: F kn ( k x k ) T g T which means that F k performs the mapping F k :!. The predictor can be realized by multiple use of a state-space model. The minimization problem can be stated as k o = arg min k 2 k ( k k ) (4) subject to k = F k ( k x k ) 2 x k ; given 1
2 For later use we dene a constraint-set Z k utilizing an expression for the state predictor Z k (x k )=fz k =( k k ): k 2 ^ k 2 ^ k = F k ( k x k )g 8 x k 2 () 3. CONVEITY ISSUES IN OPTIMIZATION BASED CONTROL We will in the following categorize optimization problems arising in control and consider the tradeos which should be made. First, we separate the convex from the nonconvex cases. 3.1 Convex problems We rst note that problem (4) is convex if k is a convex functional and Z k is a (nonempty) convex set. The issue of k being a convex functional does not impose severe restrictions in optimizing control. The optimization criterion will typically be de- ned by some distance measure based on a norm, and any norm is a convex functional. The problem of limiting Z k (x k ) 8 x k 2 to a convex set imposes severe restrictions from a control point of view. To discuss this we rst derive sucient conditions for Z k (x k ) to be a convex set. Sucient conditions for convexity The set Z k (x k ) in () is convex if the following assumptions are satised. A1: is a convex set. A2: is a convex set. A3: The mapping F k ( x k ):! R mx :::R mx is ane for all x k 2. A4: Z k (x k ) is nonempty for all x k 2. The proof is omitted here due to page limitations. The conditions are sucient, since there exists nonane mappings F k that renders a convex set Fk (x k ). This is readily seen in the 1-dimensional case where any continuous mapping provides a convex set Fk (x k ). We will in the following discuss assumptions A1 to A4 above with particular emphasis on how they restrict the control problem formulation. Since is convex the set U must be convex, cf. (4). The control input at a given time instant is chosen from U, ie. u i 2 U R mu. If the restrictions on the m u dierent control inputs are independent U will be convex if each of the control inputs are continuous in the sense that u ij 2 [u j u j ] R 8j 2 f1 ::: m u g. Dening U by including dependency between the control inputs will in most cases not introduce problems. Typical limitations like ju ij j + ju il j1, ju ij + u il j1, or ku ij + u il k 2 1, where u ij and u il denote two dierent control inputs at a time i, do not by themselves impose nonconvexity. As a conclusion Assumption A1 does not impose severe restriction on the formulation of the control problem as long as the control inputs are continuous. The discussion of Assumption A2 is very similar to the discussion of A1 above. Hence, if the states are continuous A2 does not impose severe restriction on the problem formulation. A3 is the problematic assumption since it severely limits the class of possible predictors. Viewing linear predictors, they always satisfy this assumption. Nonlinear predictors also satisfy this condition as long as they remain linear in the control input k. Assumption A4 is related to feasibility, ie. is it possible to nd a control sequence that satises the state constraint?we will assume that this assumption always is satis- ed. From a control design point of view this should be interpreted as a controllability constraint. 3.2 Nonconvex problems We will now discuss nonconvex problems, and assume that Assumptions A1, A2 and A4 are satis- ed. This means that nonconvexity is caused by the predictor alone. By the former discussion we claim that this is an important cause for nonconvexity, hence this type of nonconvex problems is highly relevant to optimization-based control Identication method Let us rst formalize an identication problem tailored for optimization-based control. The issue is to construct a predictor which produces small prediction errors on the optimization horizon, eg. I k in (1). First we dene a set of predictors parameterized by some parameter, cf. (3). F = ff ( ) : 2 g R m (6) For the sake of convenience and notation we have limited the set of predictors above in the sense that they are time-invariant as compared to the time-varying predictors in (3). We assume to have a set of input and output data for the system in question and structure these in an information vector. The output is assumed to be a function of the current state. Further, we also dene a truncated information vector.
3 I M = fu y 1 ::: u M;1 y M g (7) y i 2 Y R my y i = g(x i ) I i = fu y 1 ::: u i;1 y i g (8) This allows us to dene a causal predictor for the system output. ^y i+jji = g(^x i+jji )=g(f j ( i ^x iji )) (9) 8i 2 I j 2 I N ^y i+jji ^x i+jji ; predicted output or state given I i ^x iji ; state estimate of x i byuseofi i I = f ::: M ; Ng j 2 I N = f ::: Ng F j refers to the j-step ahead predictor, ie. the j'th componentoff j dened in (3). Note that the predictor is time-invariant. We observe that the prediction is based on a current state estimate ^x iji and the use of the predictor F j. In the following we will make the following two simplications. We will assume that the state is known in the sense that ^x iji = x i. Hence, we do not address the state-estimation problem. We further assume that the system output function g is known which means that the system output is known, ie. ^y iji = y i. The identication problem is posed as o = arg min 2 V () (1) V ( I M F g)= l pj (y i+j ^y i+jji )(11) i2i j2i N l pj : Y Y! R Let us consider the inner part of the optimality criterion. At a given time step i 2 I, the inner sum computes the eect of the prediction error on the horizon fi+1 ::: i+ng. l pj is some metric. It may vary on the prediction horizon. In particular, it can be of interest to fade away the eect of prediction errors for large values of j. Ultimately, we would like to identify predictors which render a convex optimization problem (4). Two means of attacking this problem are (i) limiting the class of predictors F, and (ii) modifying modify the identication criterion (11). Limiting the class of predictors to comply with convexity leaves us with ane predictors, or even more limited, linear predictors. Viewing the identication problem the minimum criterion value V o = V ( o I M F g) will increase if the size of the set of predictor F decreases. Hence, there is a trade-o between the quality of the predictions and the advantage of having a convex optimization problem. Referring to (ii), limiting the class of predictor can be viewed as inferring hard constraints on the optimization problem. An alternative is to include soft constraints by modifying the identication criterion in (1). In this work we do this by using ideas from (Johansen 1996) in which the identication criterion is modied to account for prior knowledge or desired properties on a problem. The modied identication criterion is given by: V ( I M F g)= i2i i2i N l pj (y i+j ^y i+jji ) + ks F F k + kf ; F lin k (12) 2 [ 1) S F is an operator that projects the nonlinearities of F onto a normed function space. S F will typically contain 2nd order derivative terms which vanish if F is linear. kf ; F lin k is a distance measure between an identied predictor F and the best linear predictor, F lin. The norm is typically a weighted 2-norm. Let us view two possibilities. If > = the predictor nonlinearities are limited by a trade-o between the prediction error and the degree of nonlinearity of F. Second, if = > it is necessary to obtain an apriori linear predictor. The distance from this linear model is then penalized. A problem with the aboveapproach is to compute the norms of S F F or F ; F lin, in particular their gradients with respect to. It is highly desirable to have a linear parameterization since this greatly reduces the computational eort. The proposed method may imply a signicant complexity increase compared to the standard identication formulation. The identication problem is, however, an o-line problem. Further, the complexity increase is motivated by the fact that we are seeking predictors that render a convex or close-to-convex optimization-based control problem. These are on-line problems in which timeecient algorithms can be critical for stability and performance. The method adds soft constraints in the sense that an identied predictor may well produce a nonconvex problem. We do, however, know that a linear predictor renders a convex problem given Assumptions A1, A2 and A4. Hence, if the deviation from a linear predictor is penalized it seems reasonable to believe that this convexies the optimization problem. Alternative methods can be suggested by including convexity conditions as hard constraints in the identication problem.
4 3.3 Convexity measure To the authors' knowledge there does not exist any measures of the degree of convexity in the literature. To evaluate our results in quantitative terms, however, we dene a convexity measure based on the epigraph of over, cf. (Luenberger 1969), p.192. We dene the set (v x k )=f 2 :( k ) v (13) ^ k = F k ( x k ) 2g v 2 R Further we dene the convex hull of following Floudas (Floudas 199), p.21. Denition of a convexity measure Let (v x k ) be a set (convex or nonconvex) in R m. The convex hull, denoted ~(v x k ), of (v x k ) is the intersection of all convex sets in R m which contains (v x k ) as a subset. Hence, the convex hull is the smallest extension ~(v x k )of(v x k )such that ~(v x k )isconvex. If (v x k ) is a convex set, it is equal to its convex hull. The convexity measure is based on the volume (v x k ) compared to its convex hull. = 1 v ; v Z v v R ~(v x k );(v x k ) d R ~ (v xk ) d dv (14) By choosing v equal to the global minimum of ( k ) the denominator will always be larger than. v is given by subject to v = max k 2 k ( k k ) (1) k 2 h i (x i+1 x i u i )= 8 i 2 I k x k ; given Because of the normalization the convexity measure is bounded by and 1, ie. 2 [ 1). The convexity measure will be computed by discretizing uniformly, and discretizing v. The computation is extensive, in particular the calculation of the convex hull. 3.4 Parametrization The identication problems posed using V can easily become intractable. It is particularly important tochoose a parametrization which limits the complexity. There are multiple choices. In this work we choose a very simple parametrization. We assume that a linear predictor and a nonlinear predictor has been identied by means of (1). A new predictor is dened as a convex combination of these two predictors. F = p F l +(1; p )F n p 2 [ 1] (16) F l and F n are the linear and nonlinear predictors identied by (1). Hence, one parameter, p, the weighting between the linear and nonlinear predictor, needs to be estimated. p will take onvalues close to if the inuence of the terms ks F F k and kf ;F l k is small compared to the prediction error term, provided F n is a better predictor than F l. 4. EAMPLE In this example we investigate properties related to a model of a continuously stirred tank reactor with the following three reactions: A! k1 B! k2 C 2A k3x3! D E g(u)! 2E, and state-space model formulated in continuous time. _x 1 = ;k 1 (u)x 1 ; k 3 (u)x 3 x v V (x 1 ; x 1 ) _x 2 = k 1 (u)x 1 ; k 2 (u)x 2 ; v V x 2 _x 3 = g(u)x 3 + v V (x 3 ; x 3 ) (17) y = x 2 k i (u)=k i e ; E i u i 2f1 2 3g g(u)=be ;(u;us)2 Parameter values are given in Table 1. k 1 k 2 k 3 1:287e +12, 1:287e + 12, 9:43e +9 E 1 E 2 E 3 978:3, 978:3, 86: us 387:3K v=v 14:18 b 1: x 1 x 3 :1 1: Table 1. Example 2: Parameter values. x 1 x 2 x 3 u are the concentrations of A, B and E, and temperature, respectively, while y constitutes the process output. E is some biochemical substance which grows according to temperature conditions given by g(u). It uses some inert substance to support this growth. E catalysis the reaction 2A! D. We assume that B is the valuable product and hence try to maximize the productionofthisover some time-horizon. If the reaction E! 2E is removed the model resembles a Van de Vusse reactor. We dene the following criterion using a sampling time of :1. ( )=; i=1 jy i j (18) = fu u u u 1 u 1 u 1 u 1 u 1 u 1 u 1 g = fx 1 ::: x 1 g U =[u s ; : u s +:] =[1: 2:] [: 2:] [:8 2:] x ; given
5 The graph of, ie. a surface in R, using the state-space model (17) for prediction, is shown in the upper left part of Fig.3 for dierent control inputs and initial conditions x =(2:143 1:9 1:). Nonconvexity is readily observed in this case. The nonconvexity is mainly caused by g(u). The global minimum for is obtained for (u u 1 )=(2:1 :). The surface changes with dierent values x 2. The nonconvexity is, however, maintained. We will assume that the nonlinear process model is known. This forms the basis for the nonlinear predictor. This means that the optimal control input, in this case, is searched for on a graph like the upper left part of Fig.4. A linear predictor is identied on the basis of the two solid line input data sequences shown in Fig.1 and the two corresponding output sequences shown in Fig.2. This gives 2 data points. The initial condition is x =(2:143 1:9 1:) T. A constant level is removed from each of these data sequences. No noise is added to the data. The identication criterion is identical to (11) with l pj as a quadratic function, ^y iji = y i and I N = f1 ::: 1g. The linear predictor is based on multiple use of the 1-step linear predictor. We use the Akaike Final Prediction Error to choose the model structure. y i+1 =1:3y i ; :48y i;1 (19) +:9u i ; :68u i;1 Since the model is based on data in which the bias is removed the input and output in (19) must be adjusted when used for prediction. In practice the input/output data to (19) can be corrected by subtracting a lowpass-ltered value of the input and output. This is done in the prediction evaluation below. To illustrate the prediction properties of the linear model, the 1-step ahead prediction, ie. :1 time units, is shown in Fig.3 for the initial condition x =(2:1 1:84 1:2) and a control input sequence dierent from the identication data set, see the validation data in Fig.1. The prediction is quite good. It should be remarked that if the control input gets close to u s the dynamics change signicantly. No identication data for control inputs in this region is available. Hence, the prediction using the linear model will deteriorate significantly for validation data in this region. The graph of using the linear model (19) is a linear plane, thus it is convex. We use the parametrization (16) and apply (12) with =. In this case V can be written as V ( p I M F g)= l pj (y i+j ^y i+jji )(2) i2i i2i N + kf l ; F n k(1 ; p ) = i2i i2i N l pj (y i+j ^y i+jji )+ (1 ; p ) The optimal choice of p increases with an increase in (or equivalently an increase in ). Further, the predictor becomes linear, ie. p = 1, for large values of. To study convexity we compute the convexity measure (14) for 4 dierent values of. The measure is computed by discretizing the -set uniformly into 2121 blocks, and discretizing the v-axis into 1 equidistant intervals. The results using this approximation are shown in Table 2. In this case is a monotonically decreasing function of p.further, the computational load and the number of failures to reach the global minimum decreases with an increase in p. V 1 provides the values for the prediction error term of V using the identication data in Figs. 1 and 2. It increases quadratically with p since there is an exact match tothe identication data for p =. p V Average no. of functions calls %failto global minimum Table 2. Example 2: The eect of different p -values on the convexity measure, the computational load, and the search for the minimizing control inputs,, isgiven. The decision surfaces for dierent p values are shown in Fig. 4. As can be expected the minimizing control sequence changes with varying predictors. For large p values the minimizing control sequence is f;: :g instead of f2:1 :g as for the correct predictor.. DISCUSSION The simulation example shows that it is feasible to include convexity into the identication problem. The identication problem is aected in two ways. First, the prediction quality of the identied model decreases. This is clearly observed by viewing row V 1 in Table 2. Second, the identication problem becomes more complex. The complexity increase is, however, very modest because of the suggested parametrization (16). The example gives a merit to the proposed method since all the tested properties of the optimizationbased control problem, ie. the convexity measure,
6 u u_s Process output y time Fig. 1. Example 2: Identication input data (solid lines), validation data (broken line). Sampling time :1. the computational load, and the ability to reach the global minimum, improves with an increase in (or equivalently p ). This result can, however, not be generalized since there does exist examples where the proposed parametrization does not improve the convexity. The results on the computational load and the failure rate to reach a global minimum will depend on the optimization method used. In this work we use a Sequential Quadratic Programming method which is the basis for the Matlab routine constr (Grace 1994). The motivation for this choice is the fact that this is an extensively used mainstream method time Fig. 2. Example 2: Identication process output data. Solid and broken line corresponds to input data given by the upper and lower solid lines in the previous gure. Sampling time : CONCLUSIONS This study shows how it is possible to take into account the convexity of an on-line optimization problem when identifying a predictor. 7. REFERENCES Floudas, C. (199). Nonlinear and Mixed-Integer Optimization. Oxford University Press. Grace, A. (1994). Optimization Toolbox - User's guide. The MathWorks Inc. Johansen, T. A. (1996). Identication of nonlinear systems using empirical data and prior knowledge - an optimization approach. Automatica 32, 337{36. Lee, J. H. (1996). Recent advances in model predictive control and other related areas. In: Preprints CPC-V, Lake Tahoe, USA. Luenberger, D. G. (1969). Optimization by Vector Space Methods. JohnWiley. Mayne, D. Q. (199). Optimization in model based control. In: Preprints IFAC Symposium DYCORD, Helsingor, Denmark time Fig. 3. Example 2: Process output data. Solid and broken line correspond to process output and 1-step ahead prediction using the linear model. Control input data is given by the broken line data in Fig.3. Sampling time : \theta_p=. \theta_p= \theta_p=.2 \theta_p=.7 Fig. 4. Example 2: Decision surfaces, ie. (u u 1 ) in (18), for p = f: :2: :7g are shown. x =(2:143 1:9 1:) T.
1 Introduction When the model structure does not match the system, is poorly identiable, or the available set of empirical data is not suciently infor
On Tikhonov Regularization, Bias and Variance in Nonlinear System Identication Tor A. Johansen SINTEF Electronics and Cybernetics, Automatic Control Department, N-7034 Trondheim, Norway. Email: Tor.Arne.Johansen@ecy.sintef.no.
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More information1 Introduction Determining a model from a nite sample of observations without any prior knowledge about the system is an ill-posed problem, in the sen
Identication of Non-linear Systems using Empirical Data and Prior Knowledge { An Optimization Approach Tor A. Johansen 1 Department of Engineering Cybernetics, Norwegian Institute of Technology, 73 Trondheim,
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationThe Uniformity Principle: A New Tool for. Probabilistic Robustness Analysis. B. R. Barmish and C. M. Lagoa. further discussion.
The Uniformity Principle A New Tool for Probabilistic Robustness Analysis B. R. Barmish and C. M. Lagoa Department of Electrical and Computer Engineering University of Wisconsin-Madison, Madison, WI 53706
More information1 I A Q E B A I E Q 1 A ; E Q A I A (2) A : (3) A : (4)
Latin Squares Denition and examples Denition. (Latin Square) An n n Latin square, or a latin square of order n, is a square array with n symbols arranged so that each symbol appears just once in each row
More informationOutline Introduction: Problem Description Diculties Algebraic Structure: Algebraic Varieties Rank Decient Toeplitz Matrices Constructing Lower Rank St
Structured Lower Rank Approximation by Moody T. Chu (NCSU) joint with Robert E. Funderlic (NCSU) and Robert J. Plemmons (Wake Forest) March 5, 1998 Outline Introduction: Problem Description Diculties Algebraic
More informationNovel determination of dierential-equation solutions: universal approximation method
Journal of Computational and Applied Mathematics 146 (2002) 443 457 www.elsevier.com/locate/cam Novel determination of dierential-equation solutions: universal approximation method Thananchai Leephakpreeda
More informationThe model reduction algorithm proposed is based on an iterative two-step LMI scheme. The convergence of the algorithm is not analyzed but examples sho
Model Reduction from an H 1 /LMI perspective A. Helmersson Department of Electrical Engineering Linkoping University S-581 8 Linkoping, Sweden tel: +6 1 816 fax: +6 1 86 email: andersh@isy.liu.se September
More informationGlucose. Oxygen Temperature
NONLINEAR PREDICTIVE CONTROL USING LOCAL MODELS { APPLIED TO A BATCH FERMENTATION PROCESS Bjarne A. Foss, Tor A. Johansen, and Aage V. Srensen Department of Engineering Cybernetics, Norwegian Institute
More informationRelative Irradiance. Wavelength (nm)
Characterization of Scanner Sensitivity Gaurav Sharma H. J. Trussell Electrical & Computer Engineering Dept. North Carolina State University, Raleigh, NC 7695-79 Abstract Color scanners are becoming quite
More informationContents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces
Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v 250) Contents 2 Vector Spaces 1 21 Vectors in R n 1 22 The Formal Denition of a Vector Space 4 23 Subspaces 6 24 Linear Combinations and
More informationMerging Optimization and Control
Merging Optimization and Control Bjarne Foss and Tor Aksel N. Heirung February 18, 2014 Department of Engineering Cybernetics Norwegian University of Science and Technology Contents Contents ii 1 Introduction
More informationConstrained State Estimation Using the Unscented Kalman Filter
16th Mediterranean Conference on Control and Automation Congress Centre, Ajaccio, France June 25-27, 28 Constrained State Estimation Using the Unscented Kalman Filter Rambabu Kandepu, Lars Imsland and
More informationLearning with Ensembles: How. over-tting can be useful. Anders Krogh Copenhagen, Denmark. Abstract
Published in: Advances in Neural Information Processing Systems 8, D S Touretzky, M C Mozer, and M E Hasselmo (eds.), MIT Press, Cambridge, MA, pages 190-196, 1996. Learning with Ensembles: How over-tting
More informationIE 5531: Engineering Optimization I
IE 5531: Engineering Optimization I Lecture 15: Nonlinear optimization Prof. John Gunnar Carlsson November 1, 2010 Prof. John Gunnar Carlsson IE 5531: Engineering Optimization I November 1, 2010 1 / 24
More informationLinear Systems with Saturating Controls: An LMI Approach. subject to control saturation. No assumption is made concerning open-loop stability and no
Output Feedback Robust Stabilization of Uncertain Linear Systems with Saturating Controls: An LMI Approach Didier Henrion 1 Sophie Tarbouriech 1; Germain Garcia 1; Abstract : The problem of robust controller
More informationAbstract Details and results of a combined moving-horizon fault detection/isolation and parameter estimation technique will be presented. Parameters a
Paper No. 231b Moving Horizon Parameter Estimation using Qualitative Knowledge and Multiple Linear Models Edward P.Gatzke, Francis J. Doyle III Department of Chemical Engineering University of Delaware
More informationEXPERIMENTAL VALIDATION OF A MARINE PROPELLER THRUST ESTIMATION SCHEME. Luca Pivano yvind N. Smogeli Thor Inge Fossen Tor Arne Johansen
EXPERIMENTAL VALIDATION OF A MARINE PROPELLER THRUST ESTIMATION SCHEME Luca Pivano yvind N. Smogeli Thor Inge Fossen Tor Arne Johansen Department of Engineering Cybernetics, Norwegian University of Science
More information58 Appendix 1 fundamental inconsistent equation (1) can be obtained as a linear combination of the two equations in (2). This clearly implies that the
Appendix PRELIMINARIES 1. THEOREMS OF ALTERNATIVES FOR SYSTEMS OF LINEAR CONSTRAINTS Here we consider systems of linear constraints, consisting of equations or inequalities or both. A feasible solution
More informationProblem Description The problem we consider is stabilization of a single-input multiple-state system with simultaneous magnitude and rate saturations,
SEMI-GLOBAL RESULTS ON STABILIZATION OF LINEAR SYSTEMS WITH INPUT RATE AND MAGNITUDE SATURATIONS Trygve Lauvdal and Thor I. Fossen y Norwegian University of Science and Technology, N-7 Trondheim, NORWAY.
More informationIdentication and Control of Nonlinear Systems Using. Neural Network Models: Design and Stability Analysis. Marios M. Polycarpou and Petros A.
Identication and Control of Nonlinear Systems Using Neural Network Models: Design and Stability Analysis by Marios M. Polycarpou and Petros A. Ioannou Report 91-09-01 September 1991 Identication and Control
More informationRECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK
RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK TRNKA PAVEL AND HAVLENA VLADIMÍR Dept of Control Engineering, Czech Technical University, Technická 2, 166 27 Praha, Czech Republic mail:
More information1 Introduction A large proportion of all optimization problems arising in an industrial or scientic context are characterized by the presence of nonco
A Global Optimization Method, BB, for General Twice-Dierentiable Constrained NLPs I. Theoretical Advances C.S. Adjiman, S. Dallwig 1, C.A. Floudas 2 and A. Neumaier 1 Department of Chemical Engineering,
More informationCONTROL. J.S. Senent M. Martnez J. Sanchis. Departamento de Ingeniera de Sistemas y Automatica. Universidad Politecnica de Valencia.
QUADRATIC INDEX ANALYSIS IN PREDICTIVE CONTROL J.S. Senent M. Martnez J. Sanchis Departamento de Ingeniera de Sistemas y Automatica. Universidad Politecnica de Valencia. Camino de Vera, 14 E-46022 Valencia
More informationWidths. Center Fluctuations. Centers. Centers. Widths
Radial Basis Functions: a Bayesian treatment David Barber Bernhard Schottky Neural Computing Research Group Department of Applied Mathematics and Computer Science Aston University, Birmingham B4 7ET, U.K.
More informationBlowup for Hyperbolic Equations. Helge Kristian Jenssen and Carlo Sinestrari
Blowup for Hyperbolic Equations Helge Kristian Jenssen and Carlo Sinestrari Abstract. We consider dierent situations of blowup in sup-norm for hyperbolic equations. For scalar conservation laws with a
More informationNational Institute of Standards and Technology USA. Jon W. Tolle. Departments of Mathematics and Operations Research. University of North Carolina USA
Acta Numerica (1996), pp. 1{000 Sequential Quadratic Programming Paul T. Boggs Applied and Computational Mathematics Division National Institute of Standards and Technology Gaithersburg, Maryland 20899
More informationMerging Optimization and Control
Merging Optimization and Control Bjarne Foss and Tor Aksel N. Heirung September 4, 2013 Contents 1 Introduction 2 2 Optimization 3 2.1 Classes of optimization problems.................. 3 2.2 Solution
More information368 XUMING HE AND GANG WANG of convergence for the MVE estimator is n ;1=3. We establish strong consistency and functional continuity of the MVE estim
Statistica Sinica 6(1996), 367-374 CROSS-CHECKING USING THE MINIMUM VOLUME ELLIPSOID ESTIMATOR Xuming He and Gang Wang University of Illinois and Depaul University Abstract: We show that for a wide class
More informationChapter 1. Comparison-Sorting and Selecting in. Totally Monotone Matrices. totally monotone matrices can be found in [4], [5], [9],
Chapter 1 Comparison-Sorting and Selecting in Totally Monotone Matrices Noga Alon Yossi Azar y Abstract An mn matrix A is called totally monotone if for all i 1 < i 2 and j 1 < j 2, A[i 1; j 1] > A[i 1;
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More information= w 2. w 1. B j. A j. C + j1j2
Local Minima and Plateaus in Multilayer Neural Networks Kenji Fukumizu and Shun-ichi Amari Brain Science Institute, RIKEN Hirosawa 2-, Wako, Saitama 35-098, Japan E-mail: ffuku, amarig@brain.riken.go.jp
More informationON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS. Abstract. We introduce a wide class of asymmetric loss functions and show how to obtain
ON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS FUNCTIONS Michael Baron Received: Abstract We introduce a wide class of asymmetric loss functions and show how to obtain asymmetric-type optimal decision
More informationDomain and range of exponential and logarithmic function *
OpenStax-CNX module: m15461 1 Domain and range of exponential and logarithmic function * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License
More informationComputer Science Dept.
A NOTE ON COMPUTATIONAL INDISTINGUISHABILITY 1 Oded Goldreich Computer Science Dept. Technion, Haifa, Israel ABSTRACT We show that following two conditions are equivalent: 1) The existence of pseudorandom
More informationdivision of the residence time into two parts is motivated by a simple model of a ow system. The ow system is modeled as a transportation part without
Estimation of Residence Time in Continuous Flow Systems with Dynamics Torbjorn Andersson email: toand@isy.liu.se Department of Electrical Engineering Linkoping University S-581 83 Linkoping, Sweden Predrag
More informationR. Schaback. numerical method is proposed which rst minimizes each f j separately. and then applies a penalty strategy to gradually force the
A Multi{Parameter Method for Nonlinear Least{Squares Approximation R Schaback Abstract P For discrete nonlinear least-squares approximation problems f 2 (x)! min for m smooth functions f : IR n! IR a m
More informationFinal Suggested Solutions 1
Econ 2010 Mathematics for Economists 1 Final 2011 - Suggested Solutions 1 Question 1 See question 2 Homework 4 Question 2.- Dene the binary operations + 1 and + 2 on R 2 by x + 1 y = (x 1 + y 1, x 2 +
More informationy Ray of Half-line or ray through in the direction of y
Chapter LINEAR COMPLEMENTARITY PROBLEM, ITS GEOMETRY, AND APPLICATIONS. THE LINEAR COMPLEMENTARITY PROBLEM AND ITS GEOMETRY The Linear Complementarity Problem (abbreviated as LCP) is a general problem
More information1/sqrt(B) convergence 1/B convergence B
The Error Coding Method and PICTs Gareth James and Trevor Hastie Department of Statistics, Stanford University March 29, 1998 Abstract A new family of plug-in classication techniques has recently been
More informationSimultaneous Choice Models: The Sandwich Approach to Nonparametric Analysis
Simultaneous Choice Models: The Sandwich Approach to Nonparametric Analysis Natalia Lazzati y November 09, 2013 Abstract We study collective choice models from a revealed preference approach given limited
More informationApplied Mathematics &Optimization
Appl Math Optim 29: 211-222 (1994) Applied Mathematics &Optimization c 1994 Springer-Verlag New Yor Inc. An Algorithm for Finding the Chebyshev Center of a Convex Polyhedron 1 N.D.Botin and V.L.Turova-Botina
More informationMidterm 1. Every element of the set of functions is continuous
Econ 200 Mathematics for Economists Midterm Question.- Consider the set of functions F C(0, ) dened by { } F = f C(0, ) f(x) = ax b, a A R and b B R That is, F is a subset of the set of continuous functions
More informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
More informationA Stable Block Model Predictive Control with Variable Implementation Horizon
American Control Conference June 8-,. Portland, OR, USA WeB9. A Stable Block Model Predictive Control with Variable Implementation Horizon Jing Sun, Shuhao Chen, Ilya Kolmanovsky Abstract In this paper,
More information1 INTRODUCTION This paper will be concerned with a class of methods, which we shall call Illinois-type methods, for the solution of the single nonline
IMPROVED ALGORITMS OF ILLINOIS-TYPE FOR TE NUMERICAL SOLUTION OF NONLINEAR EQUATIONS J.A. FORD Department of Computer Science, University of Essex, Wivenhoe Park, Colchester, Essex, United Kingdom ABSTRACT
More information1. Introduction Let the least value of an objective function F (x), x2r n, be required, where F (x) can be calculated for any vector of variables x2r
DAMTP 2002/NA08 Least Frobenius norm updating of quadratic models that satisfy interpolation conditions 1 M.J.D. Powell Abstract: Quadratic models of objective functions are highly useful in many optimization
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationON THE COMPLEXITY OF SOLVING THE GENERALIZED SET PACKING PROBLEM APPROXIMATELY. Nimrod Megiddoy
ON THE COMPLEXITY OF SOLVING THE GENERALIZED SET PACKING PROBLEM APPROXIMATELY Nimrod Megiddoy Abstract. The generalized set packing problem (GSP ) is as follows. Given a family F of subsets of M = f mg
More informationA note on continuous behavior homomorphisms
Available online at www.sciencedirect.com Systems & Control Letters 49 (2003) 359 363 www.elsevier.com/locate/sysconle A note on continuous behavior homomorphisms P.A. Fuhrmann 1 Department of Mathematics,
More information[A + 1 ] + (1 ) v: : (b) Show: the derivative of T at v = v 0 < 0 is: = (v 0 ) (1 ) ; [A + 1 ]
Homework #2 Economics 4- Due Wednesday, October 5 Christiano. This question is designed to illustrate Blackwell's Theorem, Theorem 3.3 on page 54 of S-L. That theorem represents a set of conditions that
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
An Alternative Motivation for the Indirect Approach to Closed-loop Identication Lennart Ljung and Urban Forssell Department of Electrical Engineering Linkping University, S-581 83 Linkping, Sweden WWW:
More informationPlan of Class 4. Radial Basis Functions with moving centers. Projection Pursuit Regression and ridge. Principal Component Analysis: basic ideas
Plan of Class 4 Radial Basis Functions with moving centers Multilayer Perceptrons Projection Pursuit Regression and ridge functions approximation Principal Component Analysis: basic ideas Radial Basis
More information2 Tikhonov Regularization and ERM
Introduction Here we discusses how a class of regularization methods originally designed to solve ill-posed inverse problems give rise to regularized learning algorithms. These algorithms are kernel methods
More informationRichard DiSalvo. Dr. Elmer. Mathematical Foundations of Economics. Fall/Spring,
The Finite Dimensional Normed Linear Space Theorem Richard DiSalvo Dr. Elmer Mathematical Foundations of Economics Fall/Spring, 20-202 The claim that follows, which I have called the nite-dimensional normed
More informationBranching Rules for Minimum Congestion Multi- Commodity Flow Problems
Clemson University TigerPrints All Theses Theses 8-2012 Branching Rules for Minimum Congestion Multi- Commodity Flow Problems Cameron Megaw Clemson University, cmegaw@clemson.edu Follow this and additional
More informationand are based on the precise formulation of the (vague) concept of closeness. Traditionally,
LOCAL TOPOLOGY AND A SPECTRAL THEOREM Thomas Jech 1 1. Introduction. The concepts of continuity and convergence pervade the study of functional analysis, and are based on the precise formulation of the
More informationRICE UNIVERSITY. System Identication for Robust Control. Huipin Zhang. A Thesis Submitted. in Partial Fulfillment of the. Requirements for the Degree
RICE UNIVERSITY System Identication for Robust Control by Huipin Zhang A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree Master of Science Approved, Thesis Committee: Athanasios
More informationAlgorithms for Nonsmooth Optimization
Algorithms for Nonsmooth Optimization Frank E. Curtis, Lehigh University presented at Center for Optimization and Statistical Learning, Northwestern University 2 March 2018 Algorithms for Nonsmooth Optimization
More informationThese outputs can be written in a more convenient form: with y(i) = Hc m (i) n(i) y(i) = (y(i); ; y K (i)) T ; c m (i) = (c m (i); ; c m K(i)) T and n
Binary Codes for synchronous DS-CDMA Stefan Bruck, Ulrich Sorger Institute for Network- and Signal Theory Darmstadt University of Technology Merckstr. 25, 6428 Darmstadt, Germany Tel.: 49 65 629, Fax:
More information1 Introduction 198; Dugard et al, 198; Dugard et al, 198) A delay matrix in such a lower triangular form is called an interactor matrix, and almost co
Multivariable Receding-Horizon Predictive Control for Adaptive Applications Tae-Woong Yoon and C M Chow y Department of Electrical Engineering, Korea University 1, -a, Anam-dong, Sungbu-u, Seoul 1-1, Korea
More informationAUTOMATIC CONTROL COMMUNICATION SYSTEMS LINKÖPINGS UNIVERSITET. Questions AUTOMATIC CONTROL COMMUNICATION SYSTEMS LINKÖPINGS UNIVERSITET
The Problem Identification of Linear and onlinear Dynamical Systems Theme : Curve Fitting Division of Automatic Control Linköping University Sweden Data from Gripen Questions How do the control surface
More informationError Empirical error. Generalization error. Time (number of iteration)
Submitted to Neural Networks. Dynamics of Batch Learning in Multilayer Networks { Overrealizability and Overtraining { Kenji Fukumizu The Institute of Physical and Chemical Research (RIKEN) E-mail: fuku@brain.riken.go.jp
More informationExplicit Approximate Model Predictive Control of Constrained Nonlinear Systems with Quantized Input
Explicit Approximate Model Predictive Control of Constrained Nonlinear Systems with Quantized Input Alexandra Grancharova and Tor A. Johansen Institute of Control and System Research, Bulgarian Academy
More informationLearning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System
Learning Model Predictive Control for Iterative Tasks: A Computationally Efficient Approach for Linear System Ugo Rosolia Francesco Borrelli University of California at Berkeley, Berkeley, CA 94701, USA
More informationR. Balan. Splaiul Independentei 313, Bucharest, ROMANIA D. Aur
An On-line Robust Stabilizer R. Balan University "Politehnica" of Bucharest, Department of Automatic Control and Computers, Splaiul Independentei 313, 77206 Bucharest, ROMANIA radu@karla.indinf.pub.ro
More informationOPTIMAL CONSTRAINED CONTROL ALLOCATION IN MARINE SURFACE VESSELS WITH RUDDERS. Tor A. Johansen Λ Thomas P. Fuglseth Λ Petter Tøndel Λ Thor I.
OPTIMAL CONSTRAINED CONTROL ALLOCATION IN MARINE SURFACE VESSELS WITH RUDDERS Tor A. Johansen Λ Thomas P. Fuglseth Λ Petter Tøndel Λ Thor I. Fossen Λ Λ Department of Engineering Cybernetics, Norwegian
More information1. Introduction The nonlinear complementarity problem (NCP) is to nd a point x 2 IR n such that hx; F (x)i = ; x 2 IR n + ; F (x) 2 IRn + ; where F is
New NCP-Functions and Their Properties 3 by Christian Kanzow y, Nobuo Yamashita z and Masao Fukushima z y University of Hamburg, Institute of Applied Mathematics, Bundesstrasse 55, D-2146 Hamburg, Germany,
More informationFractional Roman Domination
Chapter 6 Fractional Roman Domination It is important to discuss minimality of Roman domination functions before we get into the details of fractional version of Roman domination. Minimality of domination
More informationg(.) 1/ N 1/ N Decision Decision Device u u u u CP
Distributed Weak Signal Detection and Asymptotic Relative Eciency in Dependent Noise Hakan Delic Signal and Image Processing Laboratory (BUSI) Department of Electrical and Electronics Engineering Bogazici
More informationComments on \Wavelets in Statistics: A Review" by. A. Antoniadis. Jianqing Fan. University of North Carolina, Chapel Hill
Comments on \Wavelets in Statistics: A Review" by A. Antoniadis Jianqing Fan University of North Carolina, Chapel Hill and University of California, Los Angeles I would like to congratulate Professor Antoniadis
More informationAverage Reward Parameters
Simulation-Based Optimization of Markov Reward Processes: Implementation Issues Peter Marbach 2 John N. Tsitsiklis 3 Abstract We consider discrete time, nite state space Markov reward processes which depend
More informationChapter 7 Interconnected Systems and Feedback: Well-Posedness, Stability, and Performance 7. Introduction Feedback control is a powerful approach to o
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 7 Interconnected
More informationAdaptive Track Fusion in a Multisensor Environment. in this work. It is therefore assumed that the local
Adaptive Track Fusion in a Multisensor Environment Celine Beugnon Graduate Student Mechanical & Aerospace Engineering SUNY at Bualo Bualo, NY 14260, U.S.A. beugnon@eng.bualo.edu James Llinas Center for
More informationOptimal blocking of two-level fractional factorial designs
Journal of Statistical Planning and Inference 91 (2000) 107 121 www.elsevier.com/locate/jspi Optimal blocking of two-level fractional factorial designs Runchu Zhang a;, DongKwon Park b a Department of
More informationCalculus and linear algebra for biomedical engineering Week 3: Matrices, linear systems of equations, and the Gauss algorithm
Calculus and linear algebra for biomedical engineering Week 3: Matrices, linear systems of equations, and the Gauss algorithm Hartmut Führ fuehr@matha.rwth-aachen.de Lehrstuhl A für Mathematik, RWTH Aachen
More informationDegradable Agreement in the Presence of. Byzantine Faults. Nitin H. Vaidya. Technical Report #
Degradable Agreement in the Presence of Byzantine Faults Nitin H. Vaidya Technical Report # 92-020 Abstract Consider a system consisting of a sender that wants to send a value to certain receivers. Byzantine
More informationIntrinsic diculties in using the. control theory. 1. Abstract. We point out that the natural denitions of stability and
Intrinsic diculties in using the doubly-innite time axis for input-output control theory. Tryphon T. Georgiou 2 and Malcolm C. Smith 3 Abstract. We point out that the natural denitions of stability and
More informationNear convexity, metric convexity, and convexity
Near convexity, metric convexity, and convexity Fred Richman Florida Atlantic University Boca Raton, FL 33431 28 February 2005 Abstract It is shown that a subset of a uniformly convex normed space is nearly
More informationA Geometric Framework for Nonconvex Optimization Duality using Augmented Lagrangian Functions
A Geometric Framework for Nonconvex Optimization Duality using Augmented Lagrangian Functions Angelia Nedić and Asuman Ozdaglar April 15, 2006 Abstract We provide a unifying geometric framework for the
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationCONTINUOUS TIME D=0 ZOH D 0 D=0 FOH D 0
IDENTIFICATION ASPECTS OF INTER- SAMPLE INPUT BEHAVIOR T. ANDERSSON, P. PUCAR and L. LJUNG University of Linkoping, Department of Electrical Engineering, S-581 83 Linkoping, Sweden Abstract. In this contribution
More informationMath 1180, Notes, 14 1 C. v 1 v n v 2. C A ; w n. A and w = v i w i : v w = i=1
Math 8, 9 Notes, 4 Orthogonality We now start using the dot product a lot. v v = v v n then by Recall that if w w ; w n and w = v w = nx v i w i : Using this denition, we dene the \norm", or length, of
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
More informationPARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation
PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In
More informationNonlinearity assessment of chemical processes
Scientia Iranica C (2015) 22(3), 967{980 Sharif University of Technology Scientia Iranica Transactions C: Chemistry and Chemical Engineering www.scientiairanica.com Nonlinearity assessment of chemical
More informationTechinical Proofs for Nonlinear Learning using Local Coordinate Coding
Techinical Proofs for Nonlinear Learning using Local Coordinate Coding 1 Notations and Main Results Denition 1.1 (Lipschitz Smoothness) A function f(x) on R d is (α, β, p)-lipschitz smooth with respect
More informationConvex relaxation. In example below, we have N = 6, and the cut we are considering
Convex relaxation The art and science of convex relaxation revolves around taking a non-convex problem that you want to solve, and replacing it with a convex problem which you can actually solve the solution
More information2 FRED J. HICKERNELL the sample mean of the y (i) : (2) ^ 1 N The mean square error of this estimate may be written as a sum of two parts, a bias term
GOODNESS OF FIT STATISTICS, DISCREPANCIES AND ROBUST DESIGNS FRED J. HICKERNELL Abstract. The Cramer{Von Mises goodness-of-t statistic, also known as the L 2 -star discrepancy, is the optimality criterion
More informationWerner Romisch. Humboldt University Berlin. Abstract. Perturbations of convex chance constrained stochastic programs are considered the underlying
Stability of solutions to chance constrained stochastic programs Rene Henrion Weierstrass Institute for Applied Analysis and Stochastics D-7 Berlin, Germany and Werner Romisch Humboldt University Berlin
More informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
More informationM 2 M 3. Robot M (O)
R O M A TRE DIA Universita degli Studi di Roma Tre Dipartimento di Informatica e Automazione Via della Vasca Navale, 79 { 00146 Roma, Italy Part Sequencing in Three Machine No-Wait Robotic Cells Alessandro
More informationAn Optimization-based Approach to Decentralized Assignability
2016 American Control Conference (ACC) Boston Marriott Copley Place July 6-8, 2016 Boston, MA, USA An Optimization-based Approach to Decentralized Assignability Alborz Alavian and Michael Rotkowitz Abstract
More informationON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL
ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,
More information2 IM Preprint series A, No. 1/2003 If vector x(k) = (x 1 (k) x 2 (k) ::: x n (k)) denotes the time instants in which all jobs started for the k th tim
Eigenvectors of interval matrices over max-plus algebra Katarna Cechlarova Institute of mathematics, Faculty ofscience, P.J. Safarik University, Jesenna 5, 041 54 Kosice, Slovakia, e-mail: cechlarova@science.upjs.sk
More informationCV-NP BAYESIANISM BY MCMC. Cross Validated Non Parametric Bayesianism by Markov Chain Monte Carlo CARLOS C. RODRIGUEZ
CV-NP BAYESIANISM BY MCMC Cross Validated Non Parametric Bayesianism by Markov Chain Monte Carlo CARLOS C. RODRIGUE Department of Mathematics and Statistics University at Albany, SUNY Albany NY 1, USA
More informationNew concepts: Span of a vector set, matrix column space (range) Linearly dependent set of vectors Matrix null space
Lesson 6: Linear independence, matrix column space and null space New concepts: Span of a vector set, matrix column space (range) Linearly dependent set of vectors Matrix null space Two linear systems:
More informationSupervisory Control of Petri Nets with. Uncontrollable/Unobservable Transitions. John O. Moody and Panos J. Antsaklis
Supervisory Control of Petri Nets with Uncontrollable/Unobservable Transitions John O. Moody and Panos J. Antsaklis Department of Electrical Engineering University of Notre Dame, Notre Dame, IN 46556 USA
More informationStagnation proofness and individually monotonic bargaining solutions. Jaume García-Segarra Miguel Ginés-Vilar 2013 / 04
Stagnation proofness and individually monotonic bargaining solutions Jaume García-Segarra Miguel Ginés-Vilar 2013 / 04 Stagnation proofness and individually monotonic bargaining solutions Jaume García-Segarra
More information