Generation of Linear Models using Simulation Results
|
|
- Earl Chandler
- 5 years ago
- Views:
Transcription
1 4. IMACS-Symosium MATHMOD, Wien, , Generation of Linear Models using Simulation Results Georg Otte, Sven Reitz, Joachim Haase Fraunhofer Institute for Integrated Circuits, Branch Lab Design Automation Dresden Zeunerstraße 38, D Dresden, Germany Phone: +49 (0) , Fax +49 (0) Abstract: A two ste aroach for the generation of linear models using imulse resonse values comuted by simulation will be described. In the first ste a model of a "middle" order is determined using a realization based method. In the second ste the order of this model will be reduced using nown model order reduction methods. The objective of this aer is to show, that this ind of modeling can be realized by a combination of well nown methods and leads to accurate models of very low order. 1. Introduction In many alications very fast models of as low order as ossible are required. In the last time very owerful methods have been nown which are used for the order reduction of high dimensional models [1]. Basing on the nowledge of the matrices Σ= ( A, BCD,, ) of the high dimensional model these methods determine a suitable subsace and roject transient rocesses of the system into this subsace. On the owerful methods the subsace is determined using balanced truncation or as a Krylov-subsace, resectively. If the matrices Σ= ( A, BCD,, ) are not available but only the oututs of the system, e.g. comuted via simulation, these methods cannot be used. The aim of this aer is the generation of linear time discrete models from simulation results. Given the inut and outut values of a rocess the determination of a model is called an identification tas. But the existing identification methods do not use all advantages of simulated data. For rocess identification there are established solutions, nown as for examle the "rediction error method PEM" or the "instrumental variable method IVM", resectively. These identification methods are develoed for the determination of models from exerimental, i.e. disturbed measurements. They are secialized for the handling of disturbances in the rocess outut which occur in all real measurements. But if there is a high number of inuts and oututs or a high systems order their numerical roerties are very oor []. In the case of nearly undisturbed data lie simulation results these methods are without advantage. Recently develoed subsace-methods [], [3], [4] have better numerical roerties. Otimization stes are not needed. The multi-inut multi-outut case can be handled in the same way as the single-inut single-outut one. Initial values are not needed. An imortant advantage in ractice is the simle determination of the model structure by only one arameter, the order of the model [6]. If it is ossible to determine the ste-resonse or imulse-resonse directly, simlified (basic-) versions of these methods can be used. These are much more robust and their alication area can be extended to large systems. Our aroach to determine time discrete state-sace-models directly from imulse resonse values follows this idea. Methods used here are closely related to the realization theory []. Partial realizations that are realizations of lower order aroximate the given rocess data. The aroximation is excellent at the beginning of a transient rocess, but it is not so good at the end. Often the artial realization has to be of relative high ("middle") order to obtain models with small static errors. Therefore it will be suggested in such cases to reduce the order of the analytic "middle" order model in a second ste with nown standard methods lie modal or balanced truncation, with or without singular erturbation, resectively [5]. This aer is organized lie this: In section the basic relations between system matrices, observability, controllability, Hanel matrices and the imulse resonse are resented. The construction of a suitable Hanel matrix from given imulse resonse values and the determination of the model from these Hanel matrix will be discussed. In section 3 methods for model reduction as second ste of modeling are mentioned. Section 4 contains examles and section 5 concludes the aer. 436
2 . First ste: Realization based determination of a model Suose Σ= ( A, BCD,, ) is a linear time discrete system x = 1 Ax + + Bu y = Cx + Du Here x is an N-dimensional state vector, u is a m-dimensional inut vector and y is a q-dimensional outut vector. The matrix A is N N, B is N m, C is q N and D is q m. To this system belong an observability matrix O, a controllability matrix R and a Hanel matrix H. 1 C CB CAB CA B CA O 1 CAB CA B CA B =, R = B AB A B and R H = O = 1 CA 1 CA CA B CA B If is full controllable and observable and if N, all these matrices are of the same ran N. 1 The elements of the Hanel matrix are m matrices g = CA B, 1, g1 g g3 g g3 g 4 H = g3 g4 g5 which are equal to the imulse resonse values of for 1. If the g are nown, e.g. as a simulation result, matrices H can be constructed without exlicit nowledge of the matrices A,B,C,D and the Order N. Given the g the first tas is to assemble a Hanel matrix H having a sufficient high numeric ran N n. The second tas is to determine the matrices Σ= ( A, BCD,, ), n N of an aroximate system. n.1. Construction of a Hanel matrix and choosing the model order If the system matrices and the system order are nown, the Hanel matrix of correct order N can be constructed easily. Lie in all identification rocesses, these information is not available here. In our case the matrix H is to assemble by a sufficient number of bloc-columns and bloc-rows of the g by numeric exeriments. It will be assumed that there exists a large number K of instants of the imulse resonse{ g, = 1, K}. Thus the numeric ran of H is not limited by the columns number but it can be limited by the rows number. To find a sufficient number of rows, these can be increased ste by ste until the ran does not more grow or it grows only slowly, Fig 1 numeric ran of Hanel matrices of different rows (examle) Fig singular values of Hanel matrices of a system with a high number of oututs (examle) caused by small numeric disturbances. In this way a Hanel matrix of sufficient numeric ran can be established. Fig. 1 shows the ran values versus row bloc numbers of an examle. In such cases, the order of the model looed for can be chosen as high as the ran. 437
3 In other alications, if there exits a very high number of oututs, may be, many thousands, one bloc row can be sufficient. The numeric ran of the Hanel matrix can be as high as the outut number. Here the model order cannot be chosen as high as the ran. Fig. shows the singular values of such an Hanel matrix of an examle with aroximately 100 oututs. As can be seen, the use of more bloc rows does not really change the relations. The large number of singular values which are very small and do not decrease further on is caused by numeric disturbances. In such cases the model order should be chosen near the bottom of the left art of the grah, before the bli. In other cases, if the ran of the Hanel matrix in the range of interest increases ermanently but slowly and the singular values decrease regularly slow, too, the sufficient ran and the system order are to determine by numerical exeriments. The choice of a relative high model order is not critical. Suerfluous states will be removed in the next ste... Determination of the model from Hanel matrix From H it is ossible to comute two matrices Oˆ and R ˆ with any full ran factorization, so that H = OR ˆ ˆ. Different factorizations lead to different realizations of the same system. A well suitable factorization method is the singular value decomosition (svd) which calculates matrices U,S,V with H = USV. Furthermore, the matrices Oˆ and Rˆ calculated by svd are well suitable to carry out an aroximation of the Hanel matrix if the numeric ran is too large. In S 0 V H = ( U Uneg ) USV U = + negsnegvneg 0 S. neg V neg the comonents indexed by "neg" are connected to small singular values si, n< i Nn. These comonents are negligible. This can be done, by setting si = 0, n< i Nn. From the aroximation H of ran n of the Hanel matrix H H =USV aroriate observability and controllability matrices can be found, e.g.: O = U S and R = SV. The aroximations O and R contain on the to the matrices B and C. In MATLAB notation B = R (1: n,1: m) (m: number of inuts) C = O (1: q,1: n) (q: number of oututs) Because of the shift roerty of O one gets O(1: end q,1: n) A = O( q + 1: end,1: n) A + = O(1: end q,1: n) O( q + 1: end,1: n) where ( ) + denotes the Moore-Penrose seudo inverse. These inverse is to build for a matrix of the low order n. Directly from given imulse resonse there follows D = g0. So the system Σ is determined. 3. Second ste: Order reduction Using the shortly described method, the demand of a small static error leads to models of relatively high ("middle") order n. In most alications these order can be reduced further on without significant reducing the accuracy of the model. The model matrices are nown now. This enables the alication of well established order reduction methods lie truncated modal or balanced realization, with or without singular erturbation, resectively. The main roerties of these methods are summarized in [5]. For a truncation of the model it is to transform in a suitable form. The state sace vector has to be divided in a x1 retained ν-vector x1 and a (n-ν)-vector x that can be discarded. x = x artitioned according these reordering A11 A1 B1 A = A1 A B = B By omitting the states x the following lower order system is obtained. C = [ C1 C].The matrices A, BC, have to be 438
4 ξ + 1 = A11ξ + Bu 1 η = C1ξ + Du These elimination of x can cause a static error. Therefore in the second way, the singular erturbation aroximation will be used, to correct the static behavior. Only the "dynamic" of the omitted states will be neglected. Setting x =, 1 x and eliminating +, x from the remaining equations yields (if ( E A) is nonsingular) 1 1 ξ + 1 = ( A11+ A1( E A) A1) ξ + ( B1+ ( E A) B) u 1 1 η = ( C1+ C( E A) A1) ξ + ( D+ C( E A) B) u The truncation of modal realizations is common in engineering ractice. In this case the eigenvalues of A have to be comuted and ordered by asects of ractice (e.g. ordered so that λ i, Re( λ i ) or Im( λ i ), resectively, is nondecreasing). For the balancing truncation a state transformation is needed, which arranges the (transformed) states according to their amount of energy transformation from inut to outut, i.e. according to their transfer significance. So, the states are omitted, which are most nonsignificant. It is very difficult to redict, which method in a secial alication gives the better result. This question is to answer by numerical exeriments. Thus, the order of the model which is determined in the first ste can be reduced in a second ste. Finally one gets a model of the order ν with v< n< N. 4. Examles By three examles it will be shown that the two-ste method yields very comact models and wors reliable. 4.1 An examle from MATLAB The first examle "gasf" was taen from the MATLAB Control System Toolbox, described in get_start.df,.6. The original system contains 4 oututs, 6 inuts, and 5 states. It is characterized by values with great differences (static factors from 10-5 to ) between amlifications of different inut-outut branches. The imulse resonse instants were obtained using the ste und diff functions. The start of the first ste is the construction of the Hanel matrix, the svd and the choice of order n=17. This setting is motivated by the grah of the singular values of the Hanel matrix, see the vertical blac bar in fig rel. Singulärwerte versch. Hanelmatrizen über der Aroximationsordnung 1 Bloczeile Bloczeilen 4 Bloczeilen 7 Bloczeilen Bloczeilen 0 Bloczeilen Ordnung Fig.3 Relative singular values of the Hanel matrices, constructed with a different number of bloc rows In the following only the ste resonses with the largest and with the smallest magnitude are considered. From the figures 4 and 5 you can see the good corresondence between the original and the both aroximations. Only at time ste =0 the second aroximation of order 1 doesnt match the green circle (the original). For the order reduction in the second ste modal truncation (MATLAB function modreal) was used. By choosing the order ν=1 already for the model 1 the corresondence of the ste resonses in the finish art will be very oor, comare fig
5 3 x 105 Teilsystem "MBs51" - ElemNr.: 1.5 Original reduz. Modell 1 Ordnung 17 reduz. Modell Ordnung 1 rel Fehlernorm: ? O ut m d i dimt? Fig. 4 Ste resonse from inut #1 to outut #1, ste resonse with maximal magnitude, the left diagram shows original und both aroximations at the start of ste resonse 0 x 10-5 Teilsystem "MBs56" - ElemNr.: 4-1 -? O ut m d i Original reduz. Modell 1 Ordnung 17 reduz. Modell Ordnung 1 rel Fehlernorm: dimt? Fig. 5 Ste resonse from inut #6 to outut #4, ste resonse with minimal magnitude, the left diagram shows original und both aroximations at the start of ste resonse x Teilsystem "MBs56" - ElemNr.: 4 Original reduziertes Modell rel. stat. Fehler: ? O ut m d i dimt? Fig. 6 Ste resonse from inut #6 to outut #4, original and aroximation of order 1 already in the first ste This examle demonstrates the robustness of the first modeling ste. It is caable to wor on systems with extreme different static branches. the advantage of the second ste for reduction the model order further on. 440
6 4. Temerature distribution of an analog Amlifier-IC Fig. 7 shows the locations of critical oints, where heat is roduced (blac arrows). The temerature distribution is of interest in 10 oints, mared in fig. 8 on the surface of the IC. The arrows show the locations of the temerature sensible inut stages (in our model the oututs) numbered by 36 and 41. inut 1 inut outut #36 outut 41 Fig. 7 Energy entry oints (ower transistors) Fig.8 10 outut-oints on the surface (magenta color) The ste resonses between the (normalized) two heat ower inuts of the ower stages and the 10 mared nodes were calculated by an FEM model, consisting of nodes. Fig.9 Relative singular values of the Hanel matrices, constructed with a different number of bloc rows. The bloc row consists of 10 rows. Outut No. #36 Outut No. #41 I n u t 1 I n u t Fig. 10 Temerature ste resonses of oututs #36 and #41 caused by temerature stes of the ower stages. 441
7 The order of the first model was chosen as n=1, comare fig. 9. For the order reduction in the second ste balanced truncation (MATLAB functions balreal and modred) was used. The order was reduced to ν=15. The ste resonses of the second model agree with the original and the first model aart from the first three / four time instants. In this examle already the first model is of comarative small order. Nevertheless these order can be reduced further on. The corresondence of the model and the original remains aart from the first time oints. 4.3 A micro system acceleration sensor Fig. 11 shows the scheme of a micro system acceleration sensor. The FEM model describing the outut signal as a function of the acceleration consists of 5365 nodes with a degree of freedom of 3. Therefore the system is of order Fig. 11 Scheme of a micro system acceleration sensor Fig. 1 Relative singular values of the Hanel matrices, constructed with a different number of bloc rows The singular values of the Hanel matrix in fig. 1 suggest an order n=8. In the second ste the order can be reduced to ν= without loss of accuracy. The ste resonse oututs of the original and of the two models are shown in fig. 13. Fig. 14 shows zoomed the starting sector of the ste resonse. Both of the models match the oints given. The order reduction in the second ste was carried out by modal truncation (MATLAB function modreal) Fig. 13 Ste resonses of original and models Fig. 14 Origin sector of the Ste resonses A difference between the models determined is visible in the Bode-lot, see fig. 15. As can be seen, they differ on very high frequencies at the end of the lot. 44
8 Fig. 15 Bode lot of the two accelerator models of order n=8 and order ν= 5. Conclusions The combination of well nown aroaches to calculate a artial realization via factorization of the Hanel matrix using svd and model order reduction by modal or balanced truncation yields a owerful and fast tool to generate time discrete state sace models from given imulse resonse values. The two-ste aroach yields in most alications significantly better results than the modeling of order ν already in the first ste. The method wors well on large systems too. The models are accurate and of very low order. 6. References 1. Antoulas, A.C. and Sorensen, D. C., Aroximation of large-scale dynamical systems: An overview Technical Reort, February 001, htt://www-ece.rice.edu/~aca/mtns00.df. Viberg, M., Subsace-based Methods for the Identification of Linear Time-invariant Systems Automatica, (1995) Vol. 31 No van Overschee, P. and De Moor, B., A Unifying Theorem for three Subsace System Identification Algorithms Automatica, (1995) Vol. 31 No Ljung, L., System Identification, Theory for the User, nd ed. Prentice Hall PTR, Uer Saddle River, NJ, Green, M. and Limebeer, D. J. N., Linear Robust Control Prentice Hall, Englewood Cliffs, New Jersey 0763, Varga, A., Subsace Identifiation htt:// Jan
Radial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationA Recursive Block Incomplete Factorization. Preconditioner for Adaptive Filtering Problem
Alied Mathematical Sciences, Vol. 7, 03, no. 63, 3-3 HIKARI Ltd, www.m-hiari.com A Recursive Bloc Incomlete Factorization Preconditioner for Adative Filtering Problem Shazia Javed School of Mathematical
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More information2-D Analysis for Iterative Learning Controller for Discrete-Time Systems With Variable Initial Conditions Yong FANG 1, and Tommy W. S.
-D Analysis for Iterative Learning Controller for Discrete-ime Systems With Variable Initial Conditions Yong FANG, and ommy W. S. Chow Abstract In this aer, an iterative learning controller alying to linear
More informationPosition Control of Induction Motors by Exact Feedback Linearization *
BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 8 No Sofia 008 Position Control of Induction Motors by Exact Feedback Linearization * Kostadin Kostov Stanislav Enev Farhat
More informationMODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS
MODEL-BASED MULIPLE FAUL DEECION AND ISOLAION FOR NONLINEAR SYSEMS Ivan Castillo, and homas F. Edgar he University of exas at Austin Austin, X 78712 David Hill Chemstations Houston, X 77009 Abstract A
More informationFactors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier
Australian Journal of Basic and Alied Sciences, 5(12): 2010-2020, 2011 ISSN 1991-8178 Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doed Fiber Amlifier
More informationLINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL
LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL Mohammad Bozorg Deatment of Mechanical Engineering University of Yazd P. O. Box 89195-741 Yazd Iran Fax: +98-351-750110
More informationESTIMATION OF THE OUTPUT DEVIATION NORM FOR UNCERTAIN, DISCRETE-TIME NONLINEAR SYSTEMS IN A STATE DEPENDENT FORM
Int. J. Al. Math. Comut. Sci. 2007 Vol. 17 No. 4 505 513 DOI: 10.2478/v10006-007-0042-z ESTIMATION OF THE OUTPUT DEVIATION NORM FOR UNCERTAIN DISCRETE-TIME NONLINEAR SYSTEMS IN A STATE DEPENDENT FORM PRZEMYSŁAW
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More informationIntroduction to MVC. least common denominator of all non-identical-zero minors of all order of G(s). Example: The minor of order 2: 1 2 ( s 1)
Introduction to MVC Definition---Proerness and strictly roerness A system G(s) is roer if all its elements { gij ( s)} are roer, and strictly roer if all its elements are strictly roer. Definition---Causal
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationFigure : An 8 bridge design grid. (a) Run this model using LOQO. What is the otimal comliance? What is the running time?
5.094/SMA53 Systems Otimization: Models and Comutation Assignment 5 (00 o i n ts) Due Aril 7, 004 Some Convex Analysis (0 o i n ts) (a) Given ositive scalars L and E, consider the following set in three-dimensional
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationOn Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.**
On Fractional Predictive PID Controller Design Method Emmanuel Edet*. Reza Katebi.** * echnology and Innovation Centre, Level 4, Deartment of Electronic and Electrical Engineering, University of Strathclyde,
More informationSpecialized and hybrid Newton schemes for the matrix pth root
Secialized and hybrid Newton schemes for the matrix th root Braulio De Abreu Marlliny Monsalve Marcos Raydan June 15, 2006 Abstract We discuss different variants of Newton s method for comuting the th
More informationWolfgang POESSNECKER and Ulrich GROSS*
Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND
More informationChater Matrix Norms and Singular Value Decomosition Introduction In this lecture, we introduce the notion of a norm for matrices The singular value de
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A Dahleh George Verghese Deartment of Electrical Engineering and Comuter Science Massachuasetts Institute of Technology c Chater Matrix Norms
More informationGradient Based Iterative Algorithms for Solving a Class of Matrix Equations
1216 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 50, NO. 8, AUGUST 2005 Exanding (A3)gives M 0 8 T +1 T +18 =( T +1 +1 ) F = I (A4) g = 0F 8 T +1 =( T +1 +1 ) f = T +1 +1 + T +18 F 8 T +1 =( T +1 +1 )
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More information1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum
Estimation of fluid mobility from frequency deendent azimuthal AVO a synthetic model study Yingrui Ren 1*, Xiaoyang Wu 2, Mark Chaman 1 and Xiangyang Li 2,3 1 University of Edinburgh, 2 British Geological
More informationDETERMINISTIC STOCHASTIC SUBSPACE IDENTIFICATION FOR BRIDGES
DEERMINISIC SOCHASIC SBSPACE IDENIFICAION FOR BRIDGES H. hai, V. DeBrunner, L. S. DeBrunner, J. P. Havlicek 2, K. Mish 3, K. Ford 2, A. Medda Deartment o Electrical and Comuter Engineering Florida State
More informationObserver/Kalman Filter Time Varying System Identification
Observer/Kalman Filter Time Varying System Identification Manoranjan Majji Texas A&M University, College Station, Texas, USA Jer-Nan Juang 2 National Cheng Kung University, Tainan, Taiwan and John L. Junins
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More information216 S. Chandrasearan and I.C.F. Isen Our results dier from those of Sun [14] in two asects: we assume that comuted eigenvalues or singular values are
Numer. Math. 68: 215{223 (1994) Numerische Mathemati c Sringer-Verlag 1994 Electronic Edition Bacward errors for eigenvalue and singular value decomositions? S. Chandrasearan??, I.C.F. Isen??? Deartment
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationFrequency-Weighted Robust Fault Reconstruction Using a Sliding Mode Observer
Frequency-Weighted Robust Fault Reconstruction Using a Sliding Mode Observer C.P. an + F. Crusca # M. Aldeen * + School of Engineering, Monash University Malaysia, 2 Jalan Kolej, Bandar Sunway, 4650 Petaling,
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationDetection Algorithm of Particle Contamination in Reticle Images with Continuous Wavelet Transform
Detection Algorithm of Particle Contamination in Reticle Images with Continuous Wavelet Transform Chaoquan Chen and Guoing Qiu School of Comuter Science and IT Jubilee Camus, University of Nottingham Nottingham
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More informationRobust Predictive Control of Input Constraints and Interference Suppression for Semi-Trailer System
Vol.7, No.7 (4),.37-38 htt://dx.doi.org/.457/ica.4.7.7.3 Robust Predictive Control of Inut Constraints and Interference Suression for Semi-Trailer System Zhao, Yang Electronic and Information Technology
More informationy p 2 p 1 y p Flexible System p 2 y p2 p1 y p u, y 1 -u, y 2 Component Breakdown z 1 w 1 1 y p 1 P 1 P 2 w 2 z 2
ROBUSTNESS OF FLEXIBLE SYSTEMS WITH COMPONENT-LEVEL UNCERTAINTIES Peiman G. Maghami Λ NASA Langley Research Center, Hamton, VA 368 Robustness of flexible systems in the resence of model uncertainties at
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationPassive Identification is Non Stationary Objects With Closed Loop Control
IOP Conerence Series: Materials Science and Engineering PAPER OPEN ACCESS Passive Identiication is Non Stationary Obects With Closed Loo Control To cite this article: Valeriy F Dyadik et al 2016 IOP Con.
More informationHow to Estimate Expected Shortfall When Probabilities Are Known with Interval or Fuzzy Uncertainty
How to Estimate Exected Shortfall When Probabilities Are Known with Interval or Fuzzy Uncertainty Christian Servin Information Technology Deartment El Paso Community College El Paso, TX 7995, USA cservin@gmail.com
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
More informationSTABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS
STABILITY ANALYSIS TOOL FOR TUNING UNCONSTRAINED DECENTRALIZED MODEL PREDICTIVE CONTROLLERS Massimo Vaccarini Sauro Longhi M. Reza Katebi D.I.I.G.A., Università Politecnica delle Marche, Ancona, Italy
More informationMultivariable Generalized Predictive Scheme for Gas Turbine Control in Combined Cycle Power Plant
Multivariable Generalized Predictive Scheme for Gas urbine Control in Combined Cycle Power Plant L.X.Niu and X.J.Liu Deartment of Automation North China Electric Power University Beiing, China, 006 e-mail
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationUncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition
TNN-2007-P-0332.R1 1 Uncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition Haiing Lu, K.N. Plataniotis and A.N. Venetsanooulos The Edward S. Rogers
More informationA Toolbox for Validation of Nonlinear Finite Element Models
A oolbox for Validation of Nonlinear Finite Element Models Dr.MarkC.Anderson ACA Inc. 279 Skyark Drive, Suite 3 orrance, CA 955-5345 Phone: (3) 53-8 Email: anderson@actainc.com imothy K. Hasselman ACA
More informationThe extreme case of the anisothermal calorimeter when there is no heat exchange is the adiabatic calorimeter.
.4. Determination of the enthaly of solution of anhydrous and hydrous sodium acetate by anisothermal calorimeter, and the enthaly of melting of ice by isothermal heat flow calorimeter Theoretical background
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationStudy of the circulation theory of the cooling system in vertical evaporative cooling generator
358 Science in China: Series E Technological Sciences 006 Vol.49 No.3 358 364 DOI: 10.1007/s11431-006-0358-1 Study of the circulation theory of the cooling system in vertical evaorative cooling generator
More informationOil Temperature Control System PID Controller Algorithm Analysis Research on Sliding Gear Reducer
Key Engineering Materials Online: 2014-08-11 SSN: 1662-9795, Vol. 621, 357-364 doi:10.4028/www.scientific.net/kem.621.357 2014 rans ech Publications, Switzerland Oil emerature Control System PD Controller
More informationFuzzy Automata Induction using Construction Method
Journal of Mathematics and Statistics 2 (2): 395-4, 26 ISSN 1549-3644 26 Science Publications Fuzzy Automata Induction using Construction Method 1,2 Mo Zhi Wen and 2 Wan Min 1 Center of Intelligent Control
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationNonlinear Static Analysis of Cable Net Structures by Using Newton-Raphson Method
Nonlinear Static Analysis of Cable Net Structures by Using Newton-Rahson Method Sayed Mahdi Hazheer Deartment of Civil Engineering University Selangor (UNISEL) Selangor, Malaysia hazheer.ma@gmail.com Abstract
More informationVariable Selection and Model Building
LINEAR REGRESSION ANALYSIS MODULE XIII Lecture - 38 Variable Selection and Model Building Dr. Shalabh Deartment of Mathematics and Statistics Indian Institute of Technology Kanur Evaluation of subset regression
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
More informationMetrics Performance Evaluation: Application to Face Recognition
Metrics Performance Evaluation: Alication to Face Recognition Naser Zaeri, Abeer AlSadeq, and Abdallah Cherri Electrical Engineering Det., Kuwait University, P.O. Box 5969, Safat 6, Kuwait {zaery, abeer,
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationUncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition
Uncorrelated Multilinear Discriminant Analysis with Regularization and Aggregation for Tensor Object Recognition Haiing Lu, K.N. Plataniotis and A.N. Venetsanooulos The Edward S. Rogers Sr. Deartment of
More informationConvergence Analysis of Terminal ILC in the z Domain
25 American Control Conference June 8-, 25 Portlan, OR, USA WeA63 Convergence Analysis of erminal LC in the Domain Guy Gauthier, an Benoit Boulet, Member, EEE Abstract his aer shows how we can aly -transform
More informationKeywords: pile, liquefaction, lateral spreading, analysis ABSTRACT
Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,
More informationA Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split
A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com
More informationModule 4. Analysis of Statically Indeterminate Structures by the Direct Stiffness Method. Version 2 CE IIT, Kharagpur
Module Analysis of Statically Indeterminate Structures by the Direct Stiffness Method Version CE IIT, Kharagur Lesson 9 The Direct Stiffness Method: Beams (Continued) Version CE IIT, Kharagur Instructional
More informationAR PROCESSES AND SOURCES CAN BE RECONSTRUCTED FROM. Radu Balan, Alexander Jourjine, Justinian Rosca. Siemens Corporation Research
AR PROCESSES AND SOURCES CAN BE RECONSTRUCTED FROM DEGENERATE MIXTURES Radu Balan, Alexander Jourjine, Justinian Rosca Siemens Cororation Research 7 College Road East Princeton, NJ 8 fradu,jourjine,roscag@scr.siemens.com
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationFocal waveforms for various source waveforms driving a prolate-spheroidal impulse radiating antenna (IRA)
RADIO SCIENCE, VOL. 43,, doi:10.1029/2007rs003775, 2008 Focal waveforms for various source waveforms driving a rolate-sheroidal imulse radiating antenna (IRA) Serhat Altunc, 1 Carl E. Baum, 1 Christos
More informationGeneralized Coiflets: A New Family of Orthonormal Wavelets
Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University
More informationThe Recursive Fitting of Multivariate. Complex Subset ARX Models
lied Mathematical Sciences, Vol. 1, 2007, no. 23, 1129-1143 The Recursive Fitting of Multivariate Comlex Subset RX Models Jack Penm School of Finance and lied Statistics NU College of Business & conomics
More informationDETC2003/DAC AN EFFICIENT ALGORITHM FOR CONSTRUCTING OPTIMAL DESIGN OF COMPUTER EXPERIMENTS
Proceedings of DETC 03 ASME 003 Design Engineering Technical Conferences and Comuters and Information in Engineering Conference Chicago, Illinois USA, Setember -6, 003 DETC003/DAC-48760 AN EFFICIENT ALGORITHM
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationI Poles & zeros. I First-order systems. I Second-order systems. I E ect of additional poles. I E ect of zeros. I E ect of nonlinearities
EE C28 / ME C34 Lecture Chater 4 Time Resonse Alexandre Bayen Deartment of Electrical Engineering & Comuter Science University of California Berkeley Lecture abstract Toics covered in this resentation
More informationGeneral Linear Model Introduction, Classes of Linear models and Estimation
Stat 740 General Linear Model Introduction, Classes of Linear models and Estimation An aim of scientific enquiry: To describe or to discover relationshis among events (variables) in the controlled (laboratory)
More informationResearch of PMU Optimal Placement in Power Systems
Proceedings of the 5th WSEAS/IASME Int. Conf. on SYSTEMS THEORY and SCIENTIFIC COMPUTATION, Malta, Setember 15-17, 2005 (38-43) Research of PMU Otimal Placement in Power Systems TIAN-TIAN CAI, QIAN AI
More informationDynamic System Eigenvalue Extraction using a Linear Echo State Network for Small-Signal Stability Analysis a Novel Application
Dynamic System Eigenvalue Extraction using a Linear Echo State Network for Small-Signal Stability Analysis a Novel Alication Jiaqi Liang, Jing Dai, Ganesh K. Venayagamoorthy, and Ronald G. Harley Abstract
More informationRobust Performance Design of PID Controllers with Inverse Multiplicative Uncertainty
American Control Conference on O'Farrell Street San Francisco CA USA June 9 - July Robust Performance Design of PID Controllers with Inverse Multilicative Uncertainty Tooran Emami John M Watkins Senior
More informationTHE 3-DOF helicopter system is a benchmark laboratory
Vol:8, No:8, 14 LQR Based PID Controller Design for 3-DOF Helicoter System Santosh Kr. Choudhary International Science Index, Electrical and Information Engineering Vol:8, No:8, 14 waset.org/publication/9999411
More informationRUN-TO-RUN CONTROL AND PERFORMANCE MONITORING OF OVERLAY IN SEMICONDUCTOR MANUFACTURING. 3 Department of Chemical Engineering
Coyright 2002 IFAC 15th Triennial World Congress, Barcelona, Sain RUN-TO-RUN CONTROL AND PERFORMANCE MONITORING OF OVERLAY IN SEMICONDUCTOR MANUFACTURING C.A. Bode 1, B.S. Ko 2, and T.F. Edgar 3 1 Advanced
More informationF. Augustin, P. Rentrop Technische Universität München, Centre for Mathematical Sciences
NUMERICS OF THE VAN-DER-POL EQUATION WITH RANDOM PARAMETER F. Augustin, P. Rentro Technische Universität München, Centre for Mathematical Sciences Abstract. In this article, the roblem of long-term integration
More informationAnalysis of some entrance probabilities for killed birth-death processes
Analysis of some entrance robabilities for killed birth-death rocesses Master s Thesis O.J.G. van der Velde Suervisor: Dr. F.M. Sieksma July 5, 207 Mathematical Institute, Leiden University Contents Introduction
More informationCalculation of MTTF values with Markov Models for Safety Instrumented Systems
7th WEA International Conference on APPLIE COMPUTE CIENCE, Venice, Italy, November -3, 7 3 Calculation of MTTF values with Markov Models for afety Instrumented ystems BÖCÖK J., UGLJEA E., MACHMU. University
More informationNumerical Methods for Particle Tracing in Vector Fields
On-Line Visualization Notes Numerical Methods for Particle Tracing in Vector Fields Kenneth I. Joy Visualization and Grahics Research Laboratory Deartment of Comuter Science University of California, Davis
More informationAI*IA 2003 Fusion of Multiple Pattern Classifiers PART III
AI*IA 23 Fusion of Multile Pattern Classifiers PART III AI*IA 23 Tutorial on Fusion of Multile Pattern Classifiers by F. Roli 49 Methods for fusing multile classifiers Methods for fusing multile classifiers
More informationA Simple Weight Decay Can Improve. Abstract. It has been observed in numerical simulations that a weight decay can improve
In Advances in Neural Information Processing Systems 4, J.E. Moody, S.J. Hanson and R.P. Limann, eds. Morgan Kaumann Publishers, San Mateo CA, 1995,. 950{957. A Simle Weight Decay Can Imrove Generalization
More informationGeneration of Linear Models using Simulation Results
4. IMACS-Symosum MATHMOD, We, 5..003,. 436-443 Geerato of Lear Models usg Smulato Results Georg Otte, Sve Retz, Joachm Haase Frauhofer Isttute for Itegrated Crcuts, Brach Lab Desg Automato Dresde Zeuerstraße
More informationTIME VARYING SYSTEM IDENTIFICATION
TIME VARYING SYSTEM IDENTIFICATION FOR GUIDANCE, NAVIGATION AND CONTROL APPLICATIONS 1 Manoranjan Majji Research Associate, Texas A&M University, College Station, TX I warmly acnowledge my colleagues:
More informationGeneralized analysis method of engine suspensions based on bond graph modeling and feedback control theory
Generalized analysis method of ine susensions based on bond grah modeling and feedback ntrol theory P.Y. RICHARD *, C. RAMU-ERMENT **, J. BUION *, X. MOREAU **, M. LE FOL *** *Equie Automatique des ystèmes
More informationA Parallel Algorithm for Minimization of Finite Automata
A Parallel Algorithm for Minimization of Finite Automata B. Ravikumar X. Xiong Deartment of Comuter Science University of Rhode Island Kingston, RI 02881 E-mail: fravi,xiongg@cs.uri.edu Abstract In this
More informationSIMULATED ANNEALING AND JOINT MANUFACTURING BATCH-SIZING. Ruhul SARKER. Xin YAO
Yugoslav Journal of Oerations Research 13 (003), Number, 45-59 SIMULATED ANNEALING AND JOINT MANUFACTURING BATCH-SIZING Ruhul SARKER School of Comuter Science, The University of New South Wales, ADFA,
More informationEfficient algorithms for the smallest enclosing ball problem
Efficient algorithms for the smallest enclosing ball roblem Guanglu Zhou, Kim-Chuan Toh, Jie Sun November 27, 2002; Revised August 4, 2003 Abstract. Consider the roblem of comuting the smallest enclosing
More informationDeformation Effect Simulation and Optimization for Double Front Axle Steering Mechanism
0 4th International Conference on Comuter Modeling and Simulation (ICCMS 0) IPCSIT vol. (0) (0) IACSIT Press, Singaore Deformation Effect Simulation and Otimization for Double Front Axle Steering Mechanism
More informationAN APPROACH FOR THE MODEL BASED MONITORING OF PIEZOELECTRIC ACTUATORS
II ECCOMAS THEMATIC CONFERENCE ON SMART STRUCTURES AND MATERIALS C.A. Mota Soares et al. (Eds.) Lisbon, Portugal, July 18-1, 005 AN APPROACH FOR THE MODEL BASED MONITORING OF PIEZOELECTRIC ACTUATORS Dirk
More informationCryptanalysis of Pseudorandom Generators
CSE 206A: Lattice Algorithms and Alications Fall 2017 Crytanalysis of Pseudorandom Generators Instructor: Daniele Micciancio UCSD CSE As a motivating alication for the study of lattice in crytograhy we
More informationModelling a Partly Filled Road Tanker during an Emergency Braking
Proceedings of the World Congress on Engineering and Comuter Science 217 Vol II, October 25-27, 217, San Francisco, USA Modelling a Partly Filled Road Tanker during an Emergency Braking Frank Otremba,
More informationA New GP-evolved Formulation for the Relative Permittivity of Water and Steam
ew GP-evolved Formulation for the Relative Permittivity of Water and Steam S. V. Fogelson and W. D. Potter rtificial Intelligence Center he University of Georgia, US Contact Email ddress: sergeyf1@uga.edu
More informationFinite Mixture EFA in Mplus
Finite Mixture EFA in Mlus November 16, 2007 In this document we describe the Mixture EFA model estimated in Mlus. Four tyes of deendent variables are ossible in this model: normally distributed, ordered
More informationDesign of NARMA L-2 Control of Nonlinear Inverted Pendulum
International Research Journal of Alied and Basic Sciences 016 Available online at www.irjabs.com ISSN 51-838X / Vol, 10 (6): 679-684 Science Exlorer Publications Design of NARMA L- Control of Nonlinear
More informationChapter 1 Fundamentals
Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors
More informationIterative Methods for Designing Orthogonal and Biorthogonal Two-channel FIR Filter Banks with Regularities
R. Bregović and T. Saramäi, Iterative methods for designing orthogonal and biorthogonal two-channel FIR filter bans with regularities, Proc. Of Int. Worsho on Sectral Transforms and Logic Design for Future
More informationA MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST BASED ON THE WEIBULL DISTRIBUTION
O P E R A T I O N S R E S E A R C H A N D D E C I S I O N S No. 27 DOI:.5277/ord73 Nasrullah KHAN Muhammad ASLAM 2 Kyung-Jun KIM 3 Chi-Hyuck JUN 4 A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST
More informationRoots Blower with Gradually Expanding Outlet Gap: Mathematical Modelling and Performance Simulation Yingjie Cai 1, a, Ligang Yao 2, b
2nd International Conference on Advances in Mechanical Engineering and Industrial Informatics (AMEII 216) Roots Bloer ith Gradually Exanding Outlet Ga: Mathematical Modelling and Performance Simulation
More informationThe Knuth-Yao Quadrangle-Inequality Speedup is a Consequence of Total-Monotonicity
The Knuth-Yao Quadrangle-Ineuality Seedu is a Conseuence of Total-Monotonicity Wolfgang W. Bein Mordecai J. Golin Lawrence L. Larmore Yan Zhang Abstract There exist several general techniues in the literature
More informationKEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS
4 th International Conference on Earthquake Geotechnical Engineering June 2-28, 27 KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS Misko CUBRINOVSKI 1, Hayden BOWEN 1 ABSTRACT Two methods for analysis
More information