Final exam: Automatic Control II (Reglerteknik II, 1TT495)
|
|
- Sophie Wood
- 5 years ago
- Views:
Transcription
1 Uppsala University Department of Information Technology Systems and Control Professor Torsten Söderström Final exam: Automatic Control II (Reglerteknik II, TT495) Date: October 22, 2 Responsible examiner: Torsten Söderström Preliminary grades: p, p, p. Instructions The solutions to the problems can be given in Swedish or in English, except for Problem 2 that must be solved in English. Problem is an alternative to the homework assignment. (In case you choose to hand in a solution to Problem you will be accounted for the best performance of the homework assignments and Problem.) Solve each problem on a separate page. Write your code on every page. Provide motivations for your solutions. Vague or lacking motivations may lead to a reduced number of points. Aiding material: Textbooks in automatic control (such as Reglerteori flervariabla och olinjära metoder, Reglerteknik Grundläggande teori, and others), mathematical handbooks, collection of formulas (formelsamlingar), textbooks in mathematics, calculators, copies of OH transparencies. Note that the following is not allowed: Exempelsamling med lösningar. Good luck!
2 Problem Consider the system Ü(Ø + ) Ý(Ø) Ü(Ø) Ü(Ø) Ù(Ø) (a) Is the system controllable? (b) Is the system observable? (c) Is the system controllable from the input Ù? Problem 2 Consider the system ( ) controlled with the regulator À ( ) according to the block diagram. ¹ Σ ¹ À ¹ In this case, the system is a simple DC motor. Set ( ) ( + ) À ( ) à (a) For which values of à is the closed-loop system stable? (b) Assume the system is computer-controlled, so that the closed loop system actually can be represented as in the block diagram below. Let the sampling interval be equal to. ¹ Σ ¹ À ¹ ZOH ¹ 2
3 Determine the discrete-time transfer function (Þ) that corresponds to sampling ( ). Set «. 4 points Hint. It may be an advantage to first split ( ) as a sum of two first order systems: ( ) ( + ). (c) Assume that the discrete-time regulator is a proportional regulator, so that that À (Õ) Ã. Show that the closed-loop system can be unstable for some values of Ã. 3 points Problem 3 Consider a random process in discrete-time Ý(Ø) with covariance function Ö() 5 ( ) (a) Assume that we want to make prediction steps ahead of the process. A very simple estimator of the Ý(Ø + ) from data available at time Ø would be to take ˆÝ(Ø + ) Ý(Ø) that is to take the last available measurement. How large will the prediction error variance Î [Ý(Ø + ) ˆÝ(Ø + )] 2 become? (b) A somewhat more sophisticated predictor would be to multiply the last available measurement with a constant, say ˆÝ(Ø + ) «Ý(Ø) Determine how the prediction error variance depends on the parameter «, and how it can be minimized with respect to «. 3 points (c) Determine the spectrum of the process. 4 points (d) Use the spectrum to derive a time-domain model of the process on state-space form Ü(Ø + ) Ü(Ø) + Ú(Ø) Ý(Ø) Ü(Ø) + (Ø) and specify also the variances of Ú(Ø) and (Ø). 3 points Problem 4 Consider a first order system given on polynomial form ( + Õ )Ý(Ø) Õ Ù(Ø) + ( + Õ )(Ø) where (Ø) is white noise, with zero mean and variance 2. Assume that holds. 3
4 (a) Show that the system can be written in standard state-space form as follows Ü(Ø + ) Ü(Ø) + Ù(Ø) + Ú(Ø) Ý(Ø) Ü(Ø) How are Ú and related? (b) Assume both states can be measured. Determine the optimal feedback Ù(Ø) ÄÜ(Ø) that minimizes È Ø Ý2 (Ø). 4 points (c) Determine the Kalman filter (the predictor form) giving ˆÜ(ØØ ). 5 points Problem 5 Consider a first order process Ü(Ø) observed with some noise (Ø). This may be modelled as Ü Ý Ü + Ú Ü + where Ú and are white noise processes, with intensities Ö and Ö 2, respectively. Also, assume that. (a) Determine the spectrum and the variance of the signal Ü(Ø). (b) As the signal Ü(Ø) is not measured directly, consider estimating it with an observer ˆÜ ˆÜ + à (Ý ˆÜ) Show that this can be written as ˆÜ(Ø) (Ô)Ý(Ø) and determine the transfer function operator (Ô), and sketch the character of its Bode plot, that is sketch how () varies with. Let à be a fixed parameter. (c) Examine the estimation error Ü Ü the error can be written in the form ˆÜ as a function of Ú and. Show that Ü(Ø) À (Ô)Ú(Ø) + À 2 (Ô)(Ø) and determine the two transfer function operators À (Ô) and À 2 (Ô). (d) Determine the variances of the two error terms in part (c), that is the variances of À (Ô)Ú(Ø) and of À 2 (Ô)(Ø). 4
5 (e) Determine the observer gain à that minimizes the total error variance Î (Ã) [À (Ô)Ú(Ø)] 2 + [À 2 (Ô)(Ø)] 2 (f) Determine the Kalman filter related to the state space model set up in this problem. What is the variance Ü 2 (Ø) when the Kalman filter is used? 5
6 Uppsala University Department of Information Technology Systems and Control Prof Torsten Söderström Automatic control II, October 22, 2 Answers and brief solutions Problem (a) As is nonsingular, the controllability matrix has full rank. The system is controllable. (b) As is nonsingular, the observability matrix has full rank. The system is observable. (c) As the element and the system is in diagonal form, the state Ü is not controllable. The system is not controllable using Ù. Problem 2 (a) The closed loop characteristic equation becomes + ( )À ( ) µ à (b) The poles will lie in the left half plane for all positive values of Ã. ( ) ( + ) + The sampled-data system can be written as follows. Using «, (Þ) Þ «Þ «Þ «(Þ )( «) (Þ )(Þ «) (c) The closed loop characteristic equation becomes leading to or yet + (Þ)À (Þ) (Þ )(Þ «) + à [Þ( + «) + ( ««)] Þ 2 + Þ ( «+ Ã( + «)) ßÞ Ð A necessary condition for stability is «+ Ã( ««) ßÞ Ð 2 Both these two inequalities are violated if à is large enough. 6 Þ( + «) + ( ««) (Þ )(Þ «)
7 Problem 3 (a) In this case (b) In this case Î [Ý(Ø + ) ˆÝ(Ø + )] 2 [Ý(Ø + ) Ý(Ø)] 2 Ý 2 (Ø + ) + Ý 2 (Ø) 2Ý(Ø + )Ý(Ø) Ö() + Ö() 2Ö() ( ) Î («) [Ý(Ø + ) ˆÝ(Ø + )] 2 [Ý(Ø + ) «Ý(Ø)] 2 Ö() + «2 Ö() 2«Ö() Minimizing Î («) with respect to «gives that the best value of «is «Ö()Ö(), and min Î («) Ö() Ö 2 ()Ö() ( 2 ) «which is always smaller than ( ). (c) The spectrum becomes directly from the definition () ½ Ö() ½ ½ Ö() + Ö() ½ ½ 5( ) + ½ 5( ) 5 Ö() ( )( ) 5( 2 ) cos() (d) The process is apparently an AR() process, and can be represented as a first order system with Ê 5( 2 ) Ê 2 7
8 Problem 4 (a) The transfer function operators from Ù(Ø) and Ú(Ø) to Ý(Ø) appear from the calculations Ý(Ø) Õ + Õ(Õ + ) Õ Õ(Õ + ) Õ Õ Õ + (Õ + ) Ù(Ø) + Õ + Õ(Õ + ) Ú(Ø) Ù(Ø) + Ù(Ø) + Ù(Ø) + Ú(Ø) This describes precisely the given system if Ú(Ø) (Ø + ). Ú(Ø) (b) In this case É Ì É 2. The Riccati equation becomes Ë Ì Ë + Ì Ì Ë Ë Ì Ì Ë Set Ë Spelling out the Riccati equation elementwise leads to The feedback gain becomes Ú(Ø) Ä ( Ì Ë) Ì Ë 2 (c) Set È Ô Ô 2 Ô 2 Ô 22 The Riccati equation is in this case È È Ì + Æ 2 Æ Ì È Ì È Ì È Ì È È + 2 Ô È 8
9 2 + È Ô Ô 2 Ô Ô Ô Ô 22 Ô 2 2 Ô Ô 22 Ô 2 2 Ô Hence È 2 + where is the remaining unknown to be determined. It holds leading to Ô 22 Ô 2 2 Ô µ 2 2 ( 2 ) µ with the two solutions 2 2 ( 2 ) and hence must be chosen. Therefore, È The Kalman gain becomes Ã È Ì È Ì 2 2 ( ) The Kalman filter will be ˆÜ(Ø + Ø) ( Ã)ˆÜ(ØØ ) + Ù(Ø) + ÃÝ(Ø) ˆÜ(ØØ ) + Ù(Ø) + Ý(Ø) It even follows that ˆÜ 2 (Ø + Ø). As this applies at all times, we also have ˆÜ (Ø + Ø) ˆÜ (ØØ ) + Ù(Ø) + ( )Ý(Ø) 9
10 Problem 5 (a) Ü(Ø) Ô + Ú(Ø) µ Ü() Ö The variance of Ü(Ø) can be obtained by integrating the spectrum, but even easier by solving a Lyapunov equation, (b) One gets directly 2È + Ö µ È Ü 2 (Ø) Ö 2 (Ô) Ã Ô + + à which is a first order, lowpass, filter. The static gain is Ã( + Ã). (c) Set Ü Ü ˆÜ. Then Ü ( Ü + Ú) ( ˆÜ + ÃÝ ÃˆÜ) ( Ã) Ü + Ú Ã It follows that À (Ô) Ô + + à À 2(Ô) Ã Ô + + à (d) Using part (a), [À (Ô)Ú(Ø)] 2 Ö 2( + Ã) [À 2(Ô)Ú(Ø)] 2 à 2 Ö 2 2( + Ã) (e) Î (Ã) Ö + à 2 Ö 2 2( + Ã) Î ¼ (Ã) µ ( + Ã) 2ÃÖ 2 (Ö + à 2 Õ Ö 2 ) Õ µ à 2 Ö 2 + 2ÃÖ 2 Ö µ à 2 + Ö Ö 2 Ã Ö Ö 2 The positive sign is chosen as + à must hold to guarantee stability. (f) The Riccati equation gives in this case È È + Ö È 2 Ö 2 È 2 + 2Ö 2 È Ö Ö 2 Õ È Ö Ö2 2 + Ö Ö Õ 2 à ÈÖ Ö Ö 2 which is the same result as in part (e).
Final exam: Automatic Control II (Reglerteknik II, 1TT495)
Uppsala University Department of Information Technology Systems and Control Professor Torsten Söderström Final exam: Automatic Control II (Reglerteknik II, TT495) Date: October 6, Responsible examiner:
More informationFinal exam: Computer-controlled systems (Datorbaserad styrning, 1RT450, 1TS250)
Uppsala University Department of Information Technology Systems and Control Professor Torsten Söderström Final exam: Computer-controlled systems (Datorbaserad styrning, RT450, TS250) Date: December 9,
More informationAutomatic Control III (Reglerteknik III) fall Nonlinear systems, Part 3
Automatic Control III (Reglerteknik III) fall 20 4. Nonlinear systems, Part 3 (Chapter 4) Hans Norlander Systems and Control Department of Information Technology Uppsala University OSCILLATIONS AND DESCRIBING
More informationRäkneövningar Empirisk modellering
Räkneövningar Empirisk modellering Bengt Carlsson Systems and Control Dept of Information Technology, Uppsala University 5th February 009 Abstract Räkneuppgifter samt lite kompletterande teori. Contents
More informationAutomatic Control II: Summary and comments
Automatic Control II: Summary and comments Hints for what is essential to understand the course, and to perform well at the exam. You should be able to distinguish between continuous-time (c-t) and discrete-time
More information! " # $! % & '! , ) ( + - (. ) ( ) * + / 0 1 2 3 0 / 4 5 / 6 0 ; 8 7 < = 7 > 8 7 8 9 : Œ Š ž P P h ˆ Š ˆ Œ ˆ Š ˆ Ž Ž Ý Ü Ý Ü Ý Ž Ý ê ç è ± ¹ ¼ ¹ ä ± ¹ w ç ¹ è ¼ è Œ ¹ ± ¹ è ¹ è ä ç w ¹ ã ¼ ¹ ä ¹ ¼ ¹ ±
More informationExam in Automatic Control II Reglerteknik II 5hp (1RT495)
Exam in Automatic Control II Reglerteknik II 5hp (1RT495) Date: August 4, 018 Venue: Bergsbrunnagatan 15 sal Responsible teacher: Hans Rosth. Aiding material: Calculator, mathematical handbooks, textbooks
More informationF O R SOCI AL WORK RESE ARCH
7 TH EUROPE AN CONFERENCE F O R SOCI AL WORK RESE ARCH C h a l l e n g e s i n s o c i a l w o r k r e s e a r c h c o n f l i c t s, b a r r i e r s a n d p o s s i b i l i t i e s i n r e l a t i o n
More informationCONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN. Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren
CONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren CESAME, Université catholique de Louvain Bâtiment Euler, Avenue G. Lemaître 4-6 B-1348 Louvain-la-Neuve,
More informationFRTN 15 Predictive Control
Department of AUTOMATIC CONTROL FRTN 5 Predictive Control Final Exam March 4, 27, 8am - 3pm General Instructions This is an open book exam. You may use any book you want, including the slides from the
More informationExam in Systems Engineering/Process Control
Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 7-6- Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total number
More informationT i t l e o f t h e w o r k : L a M a r e a Y o k o h a m a. A r t i s t : M a r i a n o P e n s o t t i ( P l a y w r i g h t, D i r e c t o r )
v e r. E N G O u t l i n e T i t l e o f t h e w o r k : L a M a r e a Y o k o h a m a A r t i s t : M a r i a n o P e n s o t t i ( P l a y w r i g h t, D i r e c t o r ) C o n t e n t s : T h i s w o
More informationECE 388 Automatic Control
Controllability and State Feedback Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage:
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More information4.3 Laplace Transform in Linear System Analysis
4.3 Laplace Transform in Linear System Analysis The main goal in analysis of any dynamic system is to find its response to a given input. The system response in general has two components: zero-state response
More informationAN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS
AN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS First the X, then the AR, finally the MA Jan C. Willems, K.U. Leuven Workshop on Observation and Estimation Ben Gurion University, July 3, 2004 p./2 Joint
More informationFramework for functional tree simulation applied to 'golden delicious' apple trees
Purdue University Purdue e-pubs Open Access Theses Theses and Dissertations Spring 2015 Framework for functional tree simulation applied to 'golden delicious' apple trees Marek Fiser Purdue University
More informationUNIQUE FJORDS AND THE ROYAL CAPITALS UNIQUE FJORDS & THE NORTH CAPE & UNIQUE NORTHERN CAPITALS
Q J j,. Y j, q.. Q J & j,. & x x. Q x q. ø. 2019 :. q - j Q J & 11 Y j,.. j,, q j q. : 10 x. 3 x - 1..,,. 1-10 ( ). / 2-10. : 02-06.19-12.06.19 23.06.19-03.07.19 30.06.19-10.07.19 07.07.19-17.07.19 14.07.19-24.07.19
More informationPeriodic monopoles and difference modules
Periodic monopoles and difference modules Takuro Mochizuki RIMS, Kyoto University 2018 February Introduction In complex geometry it is interesting to obtain a correspondence between objects in differential
More informationCitation Osaka Journal of Mathematics. 43(2)
TitleIrreducible representations of the Author(s) Kosuda, Masashi Citation Osaka Journal of Mathematics. 43(2) Issue 2006-06 Date Text Version publisher URL http://hdl.handle.net/094/0396 DOI Rights Osaka
More informationSelf-Testing Polynomial Functions Efficiently and over Rational Domains
Chapter 1 Self-Testing Polynomial Functions Efficiently and over Rational Domains Ronitt Rubinfeld Madhu Sudan Ý Abstract In this paper we give the first self-testers and checkers for polynomials over
More informationLab 2 Iterative methods and eigenvalue problems. Introduction. Iterative solution of the soap film problem. Beräkningsvetenskap II/NV2, HT (6)
Beräkningsvetenskap II/NV2, HT 2008 1 (6) Institutionen för informationsteknologi Teknisk databehandling Besöksadress: MIC hus 2, Polacksbacken Lägerhyddsvägen 2 Postadress: Box 337 751 05 Uppsala Telefon:
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Basics of System Identification Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont (UBC) EECE574 - Basics of
More informationF(jω) = a(jω p 1 )(jω p 2 ) Û Ö p i = b± b 2 4ac. ω c = Y X (jω) = 1. 6R 2 C 2 (jω) 2 +7RCjω+1. 1 (6jωRC+1)(jωRC+1) RC, 1. RC = p 1, p
ÓÖ Ò ÊÄ Ò Ò Û Ò Ò Ö Ý ½¾ Ù Ö ÓÖ ÖÓÑ Ö ÓÒ Ò ÄÈ ÐØ Ö ½¾ ½¾ ½» ½½ ÓÖ Ò ÊÄ Ò Ò Û Ò Ò Ö Ý ¾ Á b 2 < 4ac Û ÒÒÓØ ÓÖ Þ Û Ö Ð Ó ÒØ Ó Û Ð Ú ÕÙ Ö º ËÓÑ Ñ ÐÐ ÕÙ Ö Ö ÓÒ Ò º Ù Ö ÓÖ ½¾ ÓÖ Ù Ö ÕÙ Ö ÓÖ Ò ØÖ Ò Ö ÙÒØ ÓÒ
More informationAutomatic Control II Computer exercise 3. LQG Design
Uppsala University Information Technology Systems and Control HN,FS,KN 2000-10 Last revised by HR August 16, 2017 Automatic Control II Computer exercise 3 LQG Design Preparations: Read Chapters 5 and 9
More informationBlock-tridiagonal matrices
Block-tridiagonal matrices. p.1/31 Block-tridiagonal matrices - where do these arise? - as a result of a particular mesh-point ordering - as a part of a factorization procedure, for example when we compute
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationFinding small factors of integers. Speed of the number-field sieve. D. J. Bernstein University of Illinois at Chicago
The number-field sieve Finding small factors of integers Speed of the number-field sieve D. J. Bernstein University of Illinois at Chicago Prelude: finding denominators 87366 22322444 in R. Easily compute
More informationGeneral Neoclassical Closure Theory: Diagonalizing the Drift Kinetic Operator
General Neoclassical Closure Theory: Diagonalizing the Drift Kinetic Operator E. D. Held eheld@cc.usu.edu Utah State University General Neoclassical Closure Theory:Diagonalizing the Drift Kinetic Operator
More informationEXAM IN MODELING AND SIMULATION (TSRT62)
EXAM IN MODELING AND SIMULATION (TSRT62) SAL: ISY:s datorsalar TID: Friday 5th January 2018, kl. 8.00 12.00 KURS: TSRT62 Modeling and Simulation PROVKOD: DAT1 INSTITUTION: ISY ANTAL UPPGIFTER: 5 ANTAL
More information: œ Ö: =? À =ß> real numbers. œ the previous plane with each point translated by : Ðfor example,! is translated to :)
â SpanÖ?ß@ œ Ö =? > @ À =ß> real numbers : SpanÖ?ß@ œ Ö: =? > @ À =ß> real numbers œ the previous plane with each point translated by : Ðfor example, is translated to :) á In general: Adding a vector :
More informationLA PRISE DE CALAIS. çoys, çoys, har - dis. çoys, dis. tons, mantz, tons, Gas. c est. à ce. C est à ce. coup, c est à ce
> ƒ? @ Z [ \ _ ' µ `. l 1 2 3 z Æ Ñ 6 = Ð l sl (~131 1606) rn % & +, l r s s, r 7 nr ss r r s s s, r s, r! " # $ s s ( ) r * s, / 0 s, r 4 r r 9;: < 10 r mnz, rz, r ns, 1 s ; j;k ns, q r s { } ~ l r mnz,
More informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationprobability of k samples out of J fall in R.
Nonparametric Techniques for Density Estimation (DHS Ch. 4) n Introduction n Estimation Procedure n Parzen Window Estimation n Parzen Window Example n K n -Nearest Neighbor Estimation Introduction Suppose
More informationThis document has been prepared by Sunder Kidambi with the blessings of
Ö À Ö Ñ Ø Ò Ñ ÒØ Ñ Ý Ò Ñ À Ö Ñ Ò Ú º Ò Ì ÝÊ À Å Ú Ø Å Ê ý Ú ÒØ º ÝÊ Ú Ý Ê Ñ º Å º ² ºÅ ý ý ý ý Ö Ð º Ñ ÒÜ Æ Å Ò Ñ Ú «Ä À ý ý This document has been prepared by Sunder Kidambi with the blessings of Ö º
More informationConnection equations with stream variables are generated in a model when using the # $ % () operator or the & ' %
7 9 9 7 The two basic variable types in a connector potential (or across) variable and flow (or through) variable are not sufficient to describe in a numerically sound way the bi-directional flow of matter
More informationEXAM IN MODELING AND SIMULATION (TSRT62)
EXAM IN MODELING AND SIMULATION (TSRT62) SAL: ISY:s datorsalar TID: Wednesday 4th January 2017, kl. 8.00 12.00 KURS: TSRT62 Modeling and Simulation PROVKOD: DAT1 INSTITUTION: ISY ANTAL UPPGIFTER: 5 ANTAL
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,
More informationLecture 11: Regression Methods I (Linear Regression)
Lecture 11: Regression Methods I (Linear Regression) Fall, 2017 1 / 40 Outline Linear Model Introduction 1 Regression: Supervised Learning with Continuous Responses 2 Linear Models and Multiple Linear
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More informationLecture 11: Regression Methods I (Linear Regression)
Lecture 11: Regression Methods I (Linear Regression) 1 / 43 Outline 1 Regression: Supervised Learning with Continuous Responses 2 Linear Models and Multiple Linear Regression Ordinary Least Squares Statistical
More informationOC330C. Wiring Diagram. Recommended PKH- P35 / P50 GALH PKA- RP35 / RP50. Remarks (Drawing No.) No. Parts No. Parts Name Specifications
G G " # $ % & " ' ( ) $ * " # $ % & " ( + ) $ * " # C % " ' ( ) $ * C " # C % " ( + ) $ * C D ; E @ F @ 9 = H I J ; @ = : @ A > B ; : K 9 L 9 M N O D K P D N O Q P D R S > T ; U V > = : W X Y J > E ; Z
More informationScandinavia SUMMER / GROUPS. & Beyond 2014
/ & 2014 25 f f Fx f 211 9 Öæ Höf æ å f 807 ø 19 øø ä 2111 1 Z F ø 1328 H f fö F H å fö ö 149 H 1 ö f Hø ø Hf 1191 2089 ä ø å F ä 1907 ä 1599 H 1796 F ø ä Jä J ( ) ø F ø 19 ö ø 15 á Å f 2286 æ f ó ä H
More informationApplications of Discrete Mathematics to the Analysis of Algorithms
Applications of Discrete Mathematics to the Analysis of Algorithms Conrado Martínez Univ. Politècnica de Catalunya, Spain May 2007 Goal Given some algorithm taking inputs from some set Á, we would like
More informationContents. PART I METHODS AND CONCEPTS 2. Transfer Function Approach Frequency Domain Representations... 42
Contents Preface.............................................. xiii 1. Introduction......................................... 1 1.1 Continuous and Discrete Control Systems................. 4 1.2 Open-Loop
More informationProblem Weight Score Total 100
EE 350 EXAM IV 15 December 2010 Last Name (Print): First Name (Print): ID number (Last 4 digits): Section: DO NOT TURN THIS PAGE UNTIL YOU ARE TOLD TO DO SO Problem Weight Score 1 25 2 25 3 25 4 25 Total
More informationDerivation of the Kalman Filter
Derivation of the Kalman Filter Kai Borre Danish GPS Center, Denmark Block Matrix Identities The key formulas give the inverse of a 2 by 2 block matrix, assuming T is invertible: T U 1 L M. (1) V W N P
More informationAn Example file... log.txt
# ' ' Start of fie & %$ " 1 - : 5? ;., B - ( * * B - ( * * F I / 0. )- +, * ( ) 8 8 7 /. 6 )- +, 5 5 3 2( 7 7 +, 6 6 9( 3 5( ) 7-0 +, => - +< ( ) )- +, 7 / +, 5 9 (. 6 )- 0 * D>. C )- +, (A :, C 0 )- +,
More informationExamination paper for TFY4240 Electromagnetic theory
Department of Physics Examination paper for TFY4240 Electromagnetic theory Academic contact during examination: Associate Professor John Ove Fjærestad Phone: 97 94 00 36 Examination date: 16 December 2015
More informationLinear and Nonlinear Regression with Application to Unbalance Estimation
Linear and Nonlinear Regression with Application to Unbalance Estimation Peter Nauclér Torsten Söderström Abstract This paper considers estimation of parameters that enters nonlinearly in a regression
More informationFrequency domain representation and singular value decomposition
EOLSS Contribution 643134 Frequency domain representation and singular value decomposition AC Antoulas Department of Electrical and Computer Engineering Rice University Houston, Texas 77251-1892, USA e-mail:
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationI. D. Landau, A. Karimi: A Course on Adaptive Control Adaptive Control. Part 9: Adaptive Control with Multiple Models and Switching
I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 1 Adaptive Control Part 9: Adaptive Control with Multiple Models and Switching I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 2 Outline
More informationM30-1: Polynomial, Radical and Rational Functions, Graphs and Equations Exam /20
Class: Date: ID: A M30-1: Polynomial, Radical and Rational Functions, Graphs and Equations Exam /20 Multiple Choice Identify the choice that best completes the statement or answers the question. 1. Which
More informationLecture 15: SQP methods for equality constrained optimization
Lecture 15: SQP methods for equality constrained optimization Coralia Cartis, Mathematical Institute, University of Oxford C6.2/B2: Continuous Optimization Lecture 15: SQP methods for equality constrained
More informationMutually orthogonal latin squares (MOLS) and Orthogonal arrays (OA)
and Orthogonal arrays (OA) Bimal Roy Indian Statistical Institute, Kolkata. Bimal Roy, Indian Statistical Institute, Kolkata. and Orthogonal arrays (O Outline of the talk 1 Latin squares 2 3 Bimal Roy,
More information2 Hallén s integral equation for the thin wire dipole antenna
Ú Ð Ð ÓÒÐ Ò Ø ØØÔ»» Ѻ Ö Ùº º Ö ÁÒغ º ÁÒ Ù ØÖ Ð Å Ø Ñ Ø ÎÓк ÆÓº ¾ ¾¼½½µ ½ ¹½ ¾ ÆÙÑ Ö Ð Ñ Ø Ó ÓÖ Ò ÐÝ Ó Ö Ø ÓÒ ÖÓÑ Ø Ò Û Ö ÔÓÐ ÒØ ÒÒ Ëº À Ø ÑÞ ¹Î ÖÑ ÞÝ Ö Åº Æ Ö¹ÅÓ Êº Ë Þ ¹Ë Ò µ Ô ÖØÑ ÒØ Ó Ð ØÖ Ð Ò Ò
More informationRECURSIVE ESTIMATION AND KALMAN FILTERING
Chapter 3 RECURSIVE ESTIMATION AND KALMAN FILTERING 3. The Discrete Time Kalman Filter Consider the following estimation problem. Given the stochastic system with x k+ = Ax k + Gw k (3.) y k = Cx k + Hv
More informationMATH 174: Numerical Analysis I. Math Division, IMSP, UPLB 1 st Sem AY
MATH 74: Numerical Analysis I Math Division, IMSP, UPLB st Sem AY 0809 Eample : Prepare a table or the unction e or in [0,]. The dierence between adjacent abscissas is h step size. What should be the step
More informationLecture 16: Modern Classification (I) - Separating Hyperplanes
Lecture 16: Modern Classification (I) - Separating Hyperplanes Outline 1 2 Separating Hyperplane Binary SVM for Separable Case Bayes Rule for Binary Problems Consider the simplest case: two classes are
More informationFull file at 2CT IMPULSE, IMPULSE RESPONSE, AND CONVOLUTION CHAPTER 2CT
2CT.2 UNIT IMPULSE 2CT.2. CHAPTER 2CT 2CT.2.2 (d) We can see from plot (c) that the limit?7 Ä!, B :?7 9 œb 9$ 9. This result is consistent with the sampling property: B $ œb $ 9 9 9 9 2CT.2.3 # # $ sin
More informationExam in Systems Engineering/Process Control
Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total
More informationETIKA V PROFESII PSYCHOLÓGA
P r a ž s k á v y s o k á š k o l a p s y c h o s o c i á l n í c h s t u d i í ETIKA V PROFESII PSYCHOLÓGA N a t á l i a S l o b o d n í k o v á v e d ú c i p r á c e : P h D r. M a r t i n S t r o u
More informationLund Institute of Technology Centre for Mathematical Sciences Mathematical Statistics
Lund Institute of Technology Centre for Mathematical Sciences Mathematical Statistics STATISTICAL METHODS FOR SAFETY ANALYSIS FMS065 ÓÑÔÙØ Ö Ü Ö Ì ÓÓØ ØÖ Ô Ð ÓÖ Ø Ñ Ò Ý Ò Ò ÐÝ In this exercise we will
More informationState Observers and the Kalman filter
Modelling and Control of Dynamic Systems State Observers and the Kalman filter Prof. Oreste S. Bursi University of Trento Page 1 Feedback System State variable feedback system: Control feedback law:u =
More informationOptimal Control of PDEs
Optimal Control of PDEs Suzanne Lenhart University of Tennessee, Knoville Department of Mathematics Lecture1 p.1/36 Outline 1. Idea of diffusion PDE 2. Motivating Eample 3. Big picture of optimal control
More informationHere are proofs for some of the results about diagonalization that were presented without proof in class.
Suppose E is an 8 8 matrix. In what follows, terms like eigenvectors, eigenvalues, and eigenspaces all refer to the matrix E. Here are proofs for some of the results about diagonalization that were presented
More informationTime Series Examples Sheet
Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationSKMM 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SKMM 3023 Applied Numerical Methods Ordinary Differential Equations ibn Abdullah Faculty of Mechanical Engineering Òº ÙÐÐ ÚºÒÙÐÐ ¾¼½ SKMM 3023 Applied Numerical Methods Ordinary
More informationVectors. Teaching Learning Point. Ç, where OP. l m n
Vectors 9 Teaching Learning Point l A quantity that has magnitude as well as direction is called is called a vector. l A directed line segment represents a vector and is denoted y AB Å or a Æ. l Position
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationB œ c " " ã B œ c 8 8. such that substituting these values for the B 3 's will make all the equations true
System of Linear Equations variables Ð unknowns Ñ B" ß B# ß ÞÞÞ ß B8 Æ Æ Æ + B + B ÞÞÞ + B œ, "" " "# # "8 8 " + B + B ÞÞÞ + B œ, #" " ## # #8 8 # ã + B + B ÞÞÞ + B œ, 3" " 3# # 38 8 3 ã + 7" B" + 7# B#
More informationA Tutorial on Data Reduction. Principal Component Analysis Theoretical Discussion. By Shireen Elhabian and Aly Farag
A Tutorial on Data Reduction Principal Component Analysis Theoretical Discussion By Shireen Elhabian and Aly Farag University of Louisville, CVIP Lab November 2008 PCA PCA is A backbone of modern data
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationPH Nuclear Physics Laboratory Gamma spectroscopy (NP3)
Physics Department Royal Holloway University of London PH2510 - Nuclear Physics Laboratory Gamma spectroscopy (NP3) 1 Objectives The aim of this experiment is to demonstrate how γ-ray energy spectra may
More informationRead Section 1.1, Examples of time series, on pages 1-8. These example introduce the book; you are not tested on them.
TS Module 1 Time series overview (The attached PDF file has better formatting.)! Model building! Time series plots Read Section 1.1, Examples of time series, on pages 1-8. These example introduce the book;
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 17, pp. 76-92, 2004. Copyright 2004,. ISSN 1068-9613. ETNA STRONG RANK REVEALING CHOLESKY FACTORIZATION M. GU AND L. MIRANIAN Ý Abstract. For any symmetric
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
More informationExercises Automatic Control III 2015
Exercises Automatic Control III 205 Foreword This exercise manual is designed for the course "Automatic Control III", given by the Division of Systems and Control. The numbering of the chapters follows
More informationMargin Maximizing Loss Functions
Margin Maximizing Loss Functions Saharon Rosset, Ji Zhu and Trevor Hastie Department of Statistics Stanford University Stanford, CA, 94305 saharon, jzhu, hastie@stat.stanford.edu Abstract Margin maximizing
More information3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 REDCLIFF MUNICIPAL PLANNING COMMISSION FOR COMMENT/DISCUSSION DATE: TOPIC: April 27 th, 2018 Bylaw 1860/2018, proposed amendments to the Land Use Bylaw regarding cannabis
More informationStochastic invariances and Lamperti transformations for Stochastic Processes
Stochastic invariances and Lamperti transformations for Stochastic Processes Pierre Borgnat, Pierre-Olivier Amblard, Patrick Flandrin To cite this version: Pierre Borgnat, Pierre-Olivier Amblard, Patrick
More informationhal , version 1-27 Mar 2014
Author manuscript, published in "2nd Multidisciplinary International Conference on Scheduling : Theory and Applications (MISTA 2005), New York, NY. : United States (2005)" 2 More formally, we denote by
More informationMax. Input Power (W) Input Current (Arms) Dimming. Enclosure
Product Overview XI025100V036NM1M Input Voltage (Vac) Output Power (W) Output Voltage Range (V) Output urrent (A) Efficiency@ Max Load and 70 ase Max ase Temp. ( ) Input urrent (Arms) Max. Input Power
More information13th IFAC Symposium on System Identification
ESTIMATING PARAMETERS IN A LUMPED PARAMETER SYSTEM WITH FIRST PRINCIPLE MODELING AND DYNAMIC EXPERIMENTS RAdeCallafon UniversityofCalifornia,SanDiego Dept of Mechanical and Aerospace Engineering 9500 Gilman
More informationSME 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SME 3023 Applied Numerical Methods Ordinary Differential Equations Abu Hasan Abdullah Faculty of Mechanical Engineering Sept 2012 Abu Hasan Abdullah (FME) SME 3023 Applied
More informationOn the Performance of Random Vector Quantization Limited Feedback Beamforming in a MISO System
1 On the Performance of Random Vector Quantization Limited Feedback Beamforming in a MISO System Chun Kin Au-Yeung, Student Member, IEEE, and David J. Love, Member, IEEE Abstract In multiple antenna wireless
More informationDistributed Real-Time Control Systems
Distributed Real-Time Control Systems Chapter 9 Discrete PID Control 1 Computer Control 2 Approximation of Continuous Time Controllers Design Strategy: Design a continuous time controller C c (s) and then
More informationMotor Controller. A block diagram for the motor with a feedback controller is shown below
Motor Controller A block diagram for the motor with a feedback controller is shown below A few things to note 1. In this modeling problem, there is no established method or set of criteria for selecting
More informationFinal Exam Solutions
EE55: Linear Systems Final Exam SIST, ShanghaiTech Final Exam Solutions Course: Linear Systems Teacher: Prof. Boris Houska Duration: 85min YOUR NAME: (type in English letters) I Introduction This exam
More informationEffect of isotope substitution on the stereodynamics for O+H(D)Br OH(D)+Br reactions. 1 Introduction
J. At. Mol. Sci. doi: 10.4208/jams.052411.070811a Vol. 3, No. 2, pp. 114-121 May 2012 Effect of isotope substitution on the stereodynamics for O+H(D)Br OH(D)+Br reactions Hong Li, Bin Zheng, Ji-Qing Yin,
More informationLoop parallelization using compiler analysis
Loop parallelization using compiler analysis Which of these loops is parallel? How can we determine this automatically using compiler analysis? Organization of a Modern Compiler Source Program Front-end
More informationâ, Đ (Very Long Baseline Interferometry, VLBI)
½ 55 ½ 5 Í Vol.55 No.5 2014 9 ACTA ASTRONOMICA SINICA Sep., 2014» Á Çý è 1,2 1,2 å 1,2 Ü ô 1,2 ï 1,2 ï 1,2 à 1,3 Æ Ö 3 ý (1 Á Í 200030) (2 Á Í û À 210008) (3 541004) ÇÅè 1.5 GHz Á è, î Í, û ÓÆ Å ò ½Ò ¼ï.
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationT T V e g em D e j ) a S D } a o "m ek j g ed b m "d mq m [ d, )
. ) 6 3 ; 6 ;, G E E W T S W X D ^ L J R Y [ _ ` E ) '" " " -, 7 4-4 4-4 ; ; 7 4 4 4 4 4 ;= : " B C CA BA " ) 3D H E V U T T V e g em D e j ) a S D } a o "m ek j g ed b m "d mq m [ d, ) W X 6 G.. 6 [ X
More informationEKF and SLAM. McGill COMP 765 Sept 18 th, 2017
EKF and SLAM McGill COMP 765 Sept 18 th, 2017 Outline News and information Instructions for paper presentations Continue on Kalman filter: EKF and extension to mapping Example of a real mapping system:
More informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
More information