Lecture abstract. EE C128 / ME C134 Feedback Control Systems. Chapter outline. Intro, [1, p. 664] Lecture Chapter 12 Design via State Space
|
|
- Myron Whitehead
- 6 years ago
- Views:
Transcription
1 EE C18 / ME C134 Feedback Control Systems Lecture Chapter 1 Design via State Space Lecture abstract Alexandre Bayen Department of Electrical Engineering & Computer Science University of California Berkeley Topics covered in this presentation Controller design Controllability Observer design Observability September 10, 013 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, 013 / 58 Chapter outline 1 Design via state space 1.1 ntro Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 ntro, [1, p. 664] 1 Design via state space 1.1 ntro 1 Design via state space Design via state space Can be applied to a wider class of systems than transform methods Systems with nonlinearities Multiple-input, multiple output (MMO) systems Systems of higher order than We will focus on the application to linear systems Specify all CL poles Parameters for each CL pole Technique for finding these parameter values Cannot specify CL zero locations Sensitive to parameter changes Wide range of computational support Loss of graphic insight into a design problem Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
2 1 Design via state space State-variable FB control, [1, p. 665] 1 Design via state space State-variable FB control, [1, p. 666] n th -order CL characteristic equation (CE) There are n coe det(s A) =s n + a n 1 s n a 1 s + a 0 =0 cients whose values determine the n CL poles ntroduce n adjustable parameters into the system and relate them to the n coe cients, so that we can place the n CL poles Before, output-variable FB, now, state-variable FB Each state variable is fed back to the control, u, through a gain, k i State-variable FB gain: K CL system is plant with state-variable FB ẋ = Ax + Bu = Ax + B( Kx + r) =(A BK)x + Br y = Cx Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 1 Design via state space State-variable FB control, [1, p. 666] 1 Design via state space State-variable FB control, [1, p. 666] Figure: Plant with state-variable FB Figure: State-space representation of a plant Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Methodology, [1, p. 668] 1 Design via state space 1 Design via state space State-variable FB control in PV form, [1, p. 668] Pole placement for plants in phase-variable (PV) form 1. Represent the plant in PV form. FB each PV to the input of the plant through a gain, k i 3. Find the CE for the CL system 4. Decide upon all CL pole locations and determine equivalent CE 5. Equate like coe cients of the CE and solve for k i Plant Plant CE A = ; B = ; a 0 a 1 a... a n 1 1 C = c 1 c... c n det(s A) = s n + a n 1 s n a 1 s + a 0 =0 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
3 1 Design via state space State-variable FB control in PV form, [1, p. 668] State-variable FB CL system u = Kx; K = k 1 k... k n A BK = (a 0 + k 1 ) (a 1 + k ) (a + k 3 )... (a n 1 + K n ) CL system CE det(s (A BK)) = s n +(a n 1 + k n )s n (a 1 + k )s +(a 0 + k 1 )=0 1 Design via state space State-variable FB control in PV form, [1, p. 669] Desired CL system CE det(s (A BK)) = s n + d n 1 s n d 1 s + d 0 =0 d i = a i + k i+1 ; i =0, 1,,...,n 1 CL system TF Denominator polynomial: the CE Numerator polynomial: formed from the coe cients of the output coupling matrix, C, forsystemswritteninpvform Same for plant and CL system Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Example, [1, p. 669] 1 Design via state space 1 Design via state space Example (Controller design for PV form) Problem: Design the PV FB gains to yield %OS =9.5% Ts =0.74 seconds Solution: On the board Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Definitions, [1, p. 67] 1 Design via state space 1 Design via state space Controllability by inspection, [1, p. 673] Controllability: f an input to a system can be found that takes every state variable from a desired initial state to a desired final state, the system is said to be controllable; otherwise the system is uncontrollable. Control variable, u, can be used to control the behavior of each state variable in x Poles of the control system can be placed where we desire Determine, a priori, whether pole placement is a viable design technique for a controller When the system matrix, A, is in diagonal or parallel form, it is apparent whether or not the system is controllable A system with distinct (no repeat) eigenvalues and a diagonal system matrix, A, is controllable if the input coupling matrix, B, does not have any rows that are zero Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
4 1 Design via state space The controllability matrix, [1, p. 674] Example, [1, p. 675] 1 Design via state space n other forms, the existence of paths from the input to the state variables is not a criterion for controllability since the equations are not decoupled n th -order plant ẋ = Ax + Bu is completely controllable if the matrix C M = B AB... A n 1 B is of rank n, wherec M is called the controllability matrix Example (Controllability via the controllability matrix) Problem: Determine if the system is controllable A = ; B = Solution: On the board Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 1 Design via state space 1 Design via state space Alternative approaches to controller design, [1, p. 676] For systems not represented in PV form approaches Matching coe cients: Matching the coe cients of det(s (A BK)) with coe cients of the desired CE Same method used for systems in PV representation Transformation: Transforming the system to PV form, designing the control FB gain, & transforming the designed system back to its original state-variable representation Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, 013 / 58 1 Design via state space Approach matching coe cients, [1, p. 677] Example, [1, p. 677] 1 Design via state space Matching the coe with coe cients of cients of the desired CE det(s (A BK)) Leads to di cult calculations of the control gains, especially for higher-order systems not represented with PVs Example (Controller design by matching coe cients) Problem: Design state-variable control FB gain for the plant to yield %OS = 15% Ts =0.5 second Solution: On the board G(s) = 10 (s + 1)(s + ) Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
5 1 Design via state space Approach transformation, [1, p. 678] 1 Design via state space Approach transformation, [1, p. 678]. Design the state-variable control FB gain 3. Transform the system in PV representation back to the original representation Plant not in PV representation ż = Az + Bu y = Cz with controllability matrix C Mz = B AB A B... A n 1 B Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 1 Design via state space Approach transformation, [1, p. 678] 1 Design via state space Approach transformation, [1, p. 678] Assume the system can be transformed into the PV representation with the transformation z = Px Transformed plant. Design the control FB nclude both state-variable control FB and input u = K x x + r with controllability matrix ẋ = P 1 AP x + P 1 Bu y = CPx Transformed plant with state-variable control FB ẋ =(P 1 AP P 1 BK x )x + P 1 Br y = CPx C Mx = P 1 C Mz Zeros of this CL system are determined from the polynomial formed from the elements of CP Solving for P P = C MzC 1 Mx Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 1 Design via state space Approach transformation, [1, p. 679] Example, [1, p. 679] 1 Design via state space 3. Transform the system in PV representation back to the original representation Plant not in PV representation with state-variable control FB ż =(A BK x P 1 )z + Br y = Cz State-variable control FB gain K z = K x P 1 Example (Controller design by transformation) Problem: Design state-variable control FB gain for the plant to yield %OS = 0.8% Ts =4seconds Solution: On the board G(s) = s +4 (s + 1)(s + )(s + 5) Zeros of the CL TF are the same as the zeros of the uncompensated plant Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
6 1 Design via state space 1 Design via state space Observer motivation, [1, p. 68] Controller design relies upon access to the state variables for FB through adjustable gains Estimate states can be fed to the controller Observer: Estimator used to calculate state variables that are not accessible from the plant Plant Observer ẋ = Ax + Bu y = Cx ˆx = Aˆx + Bu ŷ = C ˆx Observer error e x =ẋ ˆx = A(x ˆx) y ŷ = C(x ˆx) = Ce x Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 1 Design via state space Observer motivation, [1, p. 683] 1 Design via state space Observer motivation, [1, p. 683] Figure: State-FB design using an OL observer Figure: State-FB design using a CL observer Dynamics of the di erence between the actual & estimated states is unforced, & if the plant is stable, this di erence, due to the di erences in initial state vectors,! 0 Speed of convergence between the actual state & the estimated state is the same as the TR of the plant since the CE of the observer error is the same as the plant Convergence is too slow,! speed up the observer and make its response time much faster than that of the controlled CL system,! the controller will receive the estimated states instantaneously Error between the outputs of the plant and the observer is fed back to the derivatives of the observer s states The system corrects to drive this error! 0 Design a desired TR into the observer that is much quicker than that of the plant or controlled CL system Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 1 Design via state space Observer motivation, [1, p. 683] 1 Design via state space Observer motivation, [1, p. 683] Observer canonical form yields an easy solution for the observer FB gain Observer FB gain, L: Ensures the TR of the observer is faster than the response of the controlled loop in order to yield a rapidly updated estimate of the state vector Figure: State-FB design using a CL observer exploded view showing FB arrangement to reduce state-variable estimation error Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
7 1 Design via state space Design methodology, [1, p. 684] 1 Design via state space 1. Find error system, i.e., state equations for error between actual state vector & estimated state vector, x ˆx. Find CE for error system ˆx = Aˆx + Bu + L(y ŷ) ŷ = C ˆx 3. Evaluate required observer FB gain, L, to meet rapid TR for observer 4. Select eigenvalues of observer to yield stability & desired TR that is faster than controlled CL response Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Definitions, [1, p. 690] 1 Design via state space 1 Design via state space Observability by inspection, [1, p. 690] Observability: f the initial-state vector, x(t 0 ), can be found from inputs, u(t), andmeasurements, y(t), over a finite interval of time from t 0, the system is said to be observable; otherwise the system is said to be unobservable. Knowledge of measured output variables, y, and control inputs, u(t), can be used to observe the behavior of each state variable in x Poles of the observer system can be placed where we desire Determine, a priori, whether pole placement is a viable design technique for an observer When the system matrix, A, is in diagonal or parallel form, it is apparent whether or not the system is observable A system with distinct (no repeat) eigenvalues and a diagonal system matrix, A, is observable if the output coupling matrix, C, does not have any columns that are zero Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 1 Design via state space The observability matrix, [1, p. 691] Example, [1, p. 691] 1 Design via state space n other forms, the existence of paths from the output to the state variables is not a criterion for observability since the equations are not decoupled n th -order plant ẋ = Ax + Bu y = Cx is completely observable if the matrix C 3 CA O M = CA n 1 Example (Observability via the observability matrix) Problem: Determine if the system is observable A = ; B = 405 ; C = Solution: On the board is of rank n, whereo M is called the observability matrix Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
8 Example, [1, p. 69] 1 Design via state space Example (Unobservability via the observability matrix) Problem: Determine if the system is observable apple 0 1 A = Solution: On the board apple 0 ; B = 1 ; C = 5 4 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Alternative approaches to controller design, [1, p. 676] Approach - matching coe cients, [1, p. 693] For systems not represented in observer canonical form approaches Matching coe cients: Matching the coe cients of det(s (A LC)) with coe cients of the desired CE Same method used for systems in PV representation Transformation: Transforming the system to observer canonical form, designing the observer FB gain, & transforming the designed system back to its original state-variable representation Matching the coe with coe cients of cients of the desired CE det(s (A LC)) Leads to di cult calculations of the observer FB gain, especially for higher-order systems not in PV representation Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Example, [1, p. 698] Approach transformation, [1, p. 695] Example (Observer design by matching coe cients) Problem: Design an observer FB gain for the system in PV representation with a TR described by =0.7!n = (s ) G(s) = (s +1.7)(s +.69) Solution: On the board. Design the observer FB gain 3. Transform the system in PV representation back to the original representation Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
9 Approach transformation, [1, p. 695] Approach transformation, [1, p. 695] Plant not in PV representation with observability matrix ż = Az + Bu y = Cz O Mz = 6 4 C CA. CA n Assume the system can be transformed into the PV, x, representation with the transformation z = Px Transformed plant with observability matrix ẋ = P 1 AP x + P 1 Bu y = CPx O Mx = O Mz P Solving for P P = O 1 Mz O Mx Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Approach transformation, [1, p. 695] Approach transformation, [1, p. 695]. Design the observer FB Transformed plant with FB gains ė x =(P 1 AP L x CP)e x y ŷ = CPe x 3. Transform the system in PV representation back to the original representation Plant not in PV representation with observer FB gain y ė z =(A ŷ = Ce z PL x C)e z Observer FB gain L z = PL x Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Example, [1, p. 695] 1 Design via state space Example (Observer design by transformation) Problem: Design an observer in cascade form Solution: On the board G(s) = 1 (s + 1)(s + )(s + 5) Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
10 s, [1, p. 700] 1 Design via state space Design systems in state-space representation for steady-state error Error is fed forward to the controlled plant via an integrator Additional state variable Plant Control FB u = ẋ N = r Cx ẋ = Ax + Bu ẋ N = y = Cx Cx + r Kx + K e x N = K K e apple x x N s, [1, p. 701] 1 Design via state space Error is fed forward to the controlled plant via an integrator Augmented representation apple apple ẋ A BK BKe = ẋ N C 0 y = C 0 apple x x N apple apple x B + x N 0 K and K e can be selected to yield the desired TR We have an additional pole to place Keep an eye on the CL zeros and their e ect on TR u + apple 0 1 r Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Example, [1, p. 701] 1 Design via state space Bibliography 1 Design via state space Example (Design of integral control) Problem: a. Design a controller without integral control to yield %OS =10% T s =0.5 second Evaluate the steady-state error for a unit step b. Repeat the design using integral control. Evaluate the steady-state error for a unit step input. Solution: On the board Norman S. Nise. Control Systems Engineering, 011. Bayen (EECS, UCB) Feedback Control Systems September 10, / 58 Bayen (EECS, UCB) Feedback Control Systems September 10, / 58
Moving Beyond Linearity
Moving Beyond Linearity The truth is never linear! 1/23 Moving Beyond Linearity The truth is never linear! r almost never! 1/23 Moving Beyond Linearity The truth is never linear! r almost never! But often
More informationWorkshop Matlab/Simulink in Drives and Power electronics Lecture 3
Workshop Matlab/Simulink in Drives and Power electronics Lecture 3 : DC-Motor Control design Ghislain REMY Jean DEPREZ 1 / 29 Workshop Program 8 lectures will be presented based on Matlab/Simulink : 1
More informationSimulation of sensor failure accommodation in flight control system of transport aircraft: a modular approach
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 11 (2015) No. 1, pp. 55-68 Simulation of sensor failure accommodation in flight control system of transport aircraft: a modular
More informationSMF Transient Voltage Suppressor Diode Series
SMF Transient Voltage Suppressor Diode Series General Information The SMF series is designed specifically to protect sensitive electronic equipment from voltage transients induced by lightning and other
More informationAM205: lecture 2. 1 These have been shifted to MD 323 for the rest of the semester.
AM205: lecture 2 Luna and Gary will hold a Python tutorial on Wednesday in 60 Oxford Street, Room 330 Assignment 1 will be posted this week Chris will hold office hours on Thursday (1:30pm 3:30pm, Pierce
More information3.5 Equations of Lines and Planes
3.5 Equations of Lines and Planes Objectives Iknowhowtodefinealineinthree-dimensionalspace. I can write a line as a parametric equation, a symmetric equation, and a vector equation. I can define a plane
More informationMATLAB Control Software Bharat Balagopal, Bharathram Balasubramanian, and Eric Stratton Green
ECE-536 DIGITAL CONTROL SYSTEMS Page 1 of 12 MATLAB Control Software Bharat Balagopal, Bharathram Balasubramanian, and Eric Stratton Green Electrical and Computer Engineering Department, North Carolina
More informationIntroduction to Mathematical Programming IE406. Lecture 4. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 4 Dr. Ted Ralphs IE406 Lecture 4 1 Reading for This Lecture Bertsimas 2.2-2.4 IE406 Lecture 4 2 The Two Crude Petroleum Example Revisited Recall the
More informationSolution to ECE680 HW3
ECE 680 Solution to HW 3 April 30, 2008 1 Solution to ECE680 HW3 In this homework assignment, we consider the one-link robot manipulator described in Example 5.13 on pp. 256 262 of the text. This system
More informationMath 5593 Linear Programming Lecture Notes
Math 5593 Linear Programming Lecture Notes Unit II: Theory & Foundations (Convex Analysis) University of Colorado Denver, Fall 2013 Topics 1 Convex Sets 1 1.1 Basic Properties (Luenberger-Ye Appendix B.1).........................
More informationRobust Pole Placement using Linear Quadratic Regulator Weight Selection Algorithm
329 Robust Pole Placement using Linear Quadratic Regulator Weight Selection Algorithm Vishwa Nath 1, R. Mitra 2 1,2 Department of Electronics and Communication Engineering, Indian Institute of Technology,
More informationSPAREPARTSCATALOG: CONNECTORS SPARE CONNECTORS KTM ART.-NR.: 3CM EN
SPAREPARTSCATALOG: CONNECTORS ART.-NR.: 3CM3208201EN CONTENT SPARE CONNECTORS AA-AN SPARE CONNECTORS AO-BC SPARE CONNECTORS BD-BQ SPARE CONNECTORS BR-CD 3 4 5 6 SPARE CONNECTORS CE-CR SPARE CONNECTORS
More information4.4 Example of [ Z Bus ] matrix formulation in the presence of mutual impedances
4.4 Example of [ Z matrix formulation in the presence of mutual impedances Consider the network shown in Fig. 4.32. Figure 4.32: The power system for [ Z example A tree for the network is shown in Fig.
More informationChapter 3. Quadric hypersurfaces. 3.1 Quadric hypersurfaces Denition.
Chapter 3 Quadric hypersurfaces 3.1 Quadric hypersurfaces. 3.1.1 Denition. Denition 1. In an n-dimensional ane space A; given an ane frame fo;! e i g: A quadric hypersurface in A is a set S consisting
More informationSPARE CONNECTORS KTM 2014
SPAREPARTSCATALOG: // ENGINE ART.-NR.: 3208201EN CONTENT CONNECTORS FOR WIRING HARNESS AA-AN CONNECTORS FOR WIRING HARNESS AO-BC CONNECTORS FOR WIRING HARNESS BD-BQ CONNECTORS FOR WIRING HARNESS BR-CD
More informationComputer Vision. Geometric Camera Calibration. Samer M Abdallah, PhD
Computer Vision Samer M Abdallah, PhD Faculty of Engineering and Architecture American University of Beirut Beirut, Lebanon Geometric Camera Calibration September 2, 2004 1 Computer Vision Geometric Camera
More informationIntroduction to Optimization Problems and Methods
Introduction to Optimization Problems and Methods wjch@umich.edu December 10, 2009 Outline 1 Linear Optimization Problem Simplex Method 2 3 Cutting Plane Method 4 Discrete Dynamic Programming Problem Simplex
More informationCS231A Course Notes 4: Stereo Systems and Structure from Motion
CS231A Course Notes 4: Stereo Systems and Structure from Motion Kenji Hata and Silvio Savarese 1 Introduction In the previous notes, we covered how adding additional viewpoints of a scene can greatly enhance
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-4: Constrained optimization Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428 June
More information(b) Linking and dynamic graph t=
1 (a) (b) (c) 2 2 2 1 1 1 6 3 4 5 6 3 4 5 6 3 4 5 7 7 7 Supplementary Figure 1: Controlling a directed tree of seven nodes. To control the whole network we need at least 3 driver nodes, which can be either
More informationInvestigating an iterative method for tuning decentralised PI-controllers in TITO systems
u 1 y 1 step G 2x2 C 2 u 2 y 2 C 1 u 1 y 1 G 2x2 u 2 y 2 step Investigating an iterative method for tuning decentralised PI-controllers in TITO systems Master s thesis in Systems, Control and Mechatronics
More informationx n x n stepnumber k order r error constant C r+1 1/2 5/12 3/8 251/720 abs. stab. interval (α,0) /11-3/10
MATH 573 LECTURE NOTES 77 13.8. Predictor-corrector methods. We consider the Adams methods, obtained from the formula xn+1 xn+1 y(x n+1 y(x n ) = y (x)dx = f(x,y(x))dx x n x n by replacing f by an interpolating
More information22ND CENTURY_J1.xls Government Site Hourly Rate
Escalation rate 000 AA0 Administrative Assistant Level I 000 AA0 Administrative Assistant Level II 000 AB0 Application Engineer Level I 000 AB0 Application Engineer Level II 000 AC0 Application Programmer
More informationtutorial i: Using matcont for numerical integration of ODEs
tutorial i: Using matcont for numerical integration of ODEs u.a. Kuznetsov Department of Mathematics Utrecht University Budapestlaan TA, Utrecht September, This session illustrates how to input a system
More informationMatrix Transformations The position of the corners of this triangle are described by the vectors: 0 1 ] 0 1 ] Transformation:
Matrix Transformations The position of the corners of this triangle are described by the vectors: [ 2 1 ] & [4 1 ] & [3 3 ] Use each of the matrices below to transform these corners. In each case, draw
More informationMultiresponse Sparse Regression with Application to Multidimensional Scaling
Multiresponse Sparse Regression with Application to Multidimensional Scaling Timo Similä and Jarkko Tikka Helsinki University of Technology, Laboratory of Computer and Information Science P.O. Box 54,
More informationUC Berkeley College of Engineering, EECS Department CS61C: Combinational Logic Blocks
2 Wawrzynek, Garcia 2004 c UCB UC Berkeley College of Engineering, EECS Department CS61C: Combinational Logic Blocks 1 Introduction Original document by J. Wawrzynek (2003-11-15) Revised by Chris Sears
More informationUC Berkeley College of Engineering, EECS Department CS61C: Combinational Logic Blocks
UC Berkeley College of Engineering, EECS Department CS61C: Combinational Logic Blocks Original document by J. Wawrzynek (2003-11-15) Revised by Chris Sears and Dan Garcia (2004-04-26) 1 Introduction Last
More informationLinear quadratic regulator
Chapter 5 Linear quadratic regulator 5.1 Introduction The linear quadratic regulator is one of the most simple optimal controller that can be developed. It is a stable and explicit controller. Its key
More informationComputational Optimization. Constrained Optimization Algorithms
Computational Optimization Constrained Optimization Algorithms Same basic algorithms Repeat Determine descent direction Determine step size Take a step Until Optimal But now must consider feasibility,
More informationQuasilinear First-Order PDEs
MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 16 Lecture 3 Quasilinear First-Order PDEs A first order quasilinear PDE is of the form a(x, y, z) + b(x, y, z) x y = c(x, y, z). (1) Such equations
More informationEE123 Digital Signal Processing
EE23 Digital Signal Processing Lecture 8 FFT II Lab based on slides by JM Kahn M Lustig, EECS UC Berkeley Announcements Last time: Started FFT Today Lab Finish FFT Read Ch 002 Midterm : Feb 22nd M Lustig,
More informationProblem. Prove that the square of any whole number n is a multiple of 4 or one more than a multiple of 4.
CHAPTER 8 Integers Problem. Prove that the square of any whole number n is a multiple of 4 or one more than a multiple of 4. Strategy 13 Use cases. This strategy may be appropriate when A problem can be
More informationAMSC/CMSC 460 Final Exam, Fall 2007
AMSC/CMSC 460 Final Exam, Fall 2007 Show all work. You may leave arithmetic expressions in any form that a calculator could evaluate. By putting your name on this paper, you agree to abide by the university
More informationDubna 2018: lines on cubic surfaces
Dubna 2018: lines on cubic surfaces Ivan Cheltsov 20th July 2018 Lecture 1: projective plane Complex plane Definition A line in C 2 is a subset that is given by ax + by + c = 0 for some complex numbers
More informationIf the center of the sphere is the origin the the equation is. x y z 2ux 2vy 2wz d 0 -(2)
Sphere Definition: A sphere is the locus of a point which remains at a constant distance from a fixed point. The fixed point is called the centre and the constant distance is the radius of the sphere.
More informationSolving Sparse Linear Systems. Forward and backward substitution for solving lower or upper triangular systems
AMSC 6 /CMSC 76 Advanced Linear Numerical Analysis Fall 7 Direct Solution of Sparse Linear Systems and Eigenproblems Dianne P. O Leary c 7 Solving Sparse Linear Systems Assumed background: Gauss elimination
More information10/26/ Solving Systems of Linear Equations Using Matrices. Objectives. Matrices
6.1 Solving Systems of Linear Equations Using Matrices Objectives Write the augmented matrix for a linear system. Perform matrix row operations. Use matrices and Gaussian elimination to solve systems.
More informationAppendix 5-1: Attachment J.1 Pricing Table -1: IMS Ceiling Loaded Rates at Contractor Site
Appendix 5-1: Attachment J.1 Pricing Table -1: IMS Ceiling Loaded Rates at Contractor Site Escalation rate 4.6% 4.6% 4.6% 4.6% 4.6% 4.6% 4.6% 4.6% 4.6% 0001 AA01 Administrative Assistant Level I $51.00
More informationAH Matrices.notebook November 28, 2016
Matrices Numbers are put into arrays to help with multiplication, division etc. A Matrix (matrices pl.) is a rectangular array of numbers arranged in rows and columns. Matrices If there are m rows and
More informationRobust Signal-Structure Reconstruction
Robust Signal-Structure Reconstruction V. Chetty 1, D. Hayden 2, J. Gonçalves 2, and S. Warnick 1 1 Information and Decision Algorithms Laboratories, Brigham Young University 2 Control Group, Department
More informationMATLAB CONTROL SYSTEM TOOLBOX IN LTI SYSTEM MODEL ANALYSIS
MATLAB CONTROL SYSTEM TOOLBOX IN LTI SYSTEM MODEL ANALYSIS Asist.univ. Luminiţa Giurgiu Abstract The MATLAB environment has important numerical tools. One of them provides a reliable foundation for control
More informationAmerican International University- Bangladesh Faculty of Engineering (EEE) Control Systems Laboratory
Experiment 1 American International University- Bangladesh Faculty of Engineering (EEE) Control Systems Laboratory Title: Introduction to design and simulation of open loop and close loop control systems
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 6:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_1.htm 1 of 1 6/20/2012 12:24 PM The Lecture deals with: ADI Method file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_2.htm 1 of 2 6/20/2012
More informationLaboratory 1 Introduction to MATLAB for Signals and Systems
Laboratory 1 Introduction to MATLAB for Signals and Systems INTRODUCTION to MATLAB MATLAB is a powerful computing environment for numeric computation and visualization. MATLAB is designed for ease of use
More informationOn the maximum rank of completions of entry pattern matrices
Linear Algebra and its Applications 525 (2017) 1 19 Contents lists available at ScienceDirect Linear Algebra and its Applications wwwelseviercom/locate/laa On the maximum rank of completions of entry pattern
More informationScicos is a Scilab toolbox included in the ScicosLab package. The Scicos editor can be opened by the scicos command
7 Getting Started 7.1 Construction of a Simple Diagram Scicos contains a graphical editor that can be used to construct block diagram models of dynamical systems. The blocks can come from various palettes
More informationCore Mathematics 1 Graphs of Functions
Regent College Maths Department Core Mathematics 1 Graphs of Functions Graphs of Functions September 2011 C1 Note Graphs of functions; sketching curves defined by simple equations. Here are some curves
More informationx = 12 x = 12 1x = 16
2.2 - The Inverse of a Matrix We've seen how to add matrices, multiply them by scalars, subtract them, and multiply one matrix by another. The question naturally arises: Can we divide one matrix by another?
More informationMATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix.
MATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix. Row echelon form A matrix is said to be in the row echelon form if the leading entries shift to the
More informationInterior Point I. Lab 21. Introduction
Lab 21 Interior Point I Lab Objective: For decades after its invention, the Simplex algorithm was the only competitive method for linear programming. The past 30 years, however, have seen the discovery
More informationNumerical Methods for PDEs : Video 9: 2D Finite Difference February 14, Equations / 29
22.520 Numerical Methods for PDEs Video 9 2D Finite Difference Equations February 4, 205 22.520 Numerical Methods for PDEs Video 9 2D Finite Difference February 4, Equations 205 / 29 Thought Experiment
More informationNumerical Analysis I - Final Exam Matrikelnummer:
Dr. Behrens Center for Mathematical Sciences Technische Universität München Winter Term 2005/2006 Name: Numerical Analysis I - Final Exam Matrikelnummer: I agree to the publication of the results of this
More informationReducing Communication Costs Associated with Parallel Algebraic Multigrid
Reducing Communication Costs Associated with Parallel Algebraic Multigrid Amanda Bienz, Luke Olson (Advisor) University of Illinois at Urbana-Champaign Urbana, IL 11 I. PROBLEM AND MOTIVATION Algebraic
More informationand n is an even positive integer, then A n is a
1) If A is a skew symmetric matrix and n is an even positive integer, then A n is a a) Symmetric Matrix b) Skew Symmetric Matrix c) Diagonal Matrix d) Scalar a Matrix Correct answer is (a) a) 5 b) 3 c)
More informationP2P VoD Systems: Modelling and Performance
P2P VoD Systems: Modelling and Performance Samuli Aalto, Aalto University, Finland Pasi Lassila, Aalto University, Finland Niklas Raatikainen, HIIT, Finland Petri Savolainen, HIIT, Finland Sasu Tarkoma,
More informationCS242: Probabilistic Graphical Models Lecture 2B: Loopy Belief Propagation & Junction Trees
CS242: Probabilistic Graphical Models Lecture 2B: Loopy Belief Propagation & Junction Trees Professor Erik Sudderth Brown University Computer Science September 22, 2016 Some figures and materials courtesy
More informationLecture 19: November 5
0-725/36-725: Convex Optimization Fall 205 Lecturer: Ryan Tibshirani Lecture 9: November 5 Scribes: Hyun Ah Song Note: LaTeX template courtesy of UC Berkeley EECS dept. Disclaimer: These notes have not
More informationMultiple View Geometry in Computer Vision
Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 Structure Computation Lecture 18 March 22, 2005 2 3D Reconstruction The goal of 3D reconstruction
More informationDM545 Linear and Integer Programming. Lecture 2. The Simplex Method. Marco Chiarandini
DM545 Linear and Integer Programming Lecture 2 The Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 3. 4. Standard Form Basic Feasible Solutions
More informationNote 10 Introduction to MATLAB & SIMULINK
Note 10 Introduction to MATLAB & SIMULINK Department of Mechanical Engineering, University Of Saskatchewan, 57 Campus Drive, Saskatoon, SK S7N 5A9, Canada 1 1 Introduction to MATLAB MATLAB stands for the
More informationParametric Curves, Lines in Space
Parametric Curves, Lines in Space Calculus III Josh Engwer TTU 02 September 2014 Josh Engwer (TTU) Parametric Curves, Lines in Space 02 September 2014 1 / 37 PART I PART I: 2D PARAMETRIC CURVES Josh Engwer
More informationEECS150 - Digital Design Lecture 14 FIFO 2 and SIFT. Recap and Outline
EECS150 - Digital Design Lecture 14 FIFO 2 and SIFT Oct. 15, 2013 Prof. Ronald Fearing Electrical Engineering and Computer Sciences University of California, Berkeley (slides courtesy of Prof. John Wawrzynek)
More informationOPTIMAL CONTROL OF A ROBOTIC SYSTEM WITH TWO DEGREE OF FREEDOM
International Journal of Electrical and Electronics Engineering (IJEEE) ISSN(P): 2278-9944; ISSN(E): 2278-9952 Vol. 4, Issue 6, Oct Nov 2015, 1-10 IASET OPTIMAL CONTROL OF A ROBOTIC SYSTEM WITH TWO DEGREE
More information.(3, 2) Co-ordinate Geometry Co-ordinates. Every point has two co-ordinates. Plot the following points on the plane. A (4, 1) D (2, 5) G (6, 3)
Co-ordinate Geometry Co-ordinates Every point has two co-ordinates. (3, 2) x co-ordinate y co-ordinate Plot the following points on the plane..(3, 2) A (4, 1) D (2, 5) G (6, 3) B (3, 3) E ( 4, 4) H (6,
More informationReal-time Model Predictive Control
Real-time Model Predictive Control MARIÁN MROSKO, EVA MIKLOVIČOVÁ Institute of Control and Industrial Informatics Faculty of Electrical Engineering and Information Technology Slovak University of Technology
More informationA Study of Numerical Methods for Simultaneous Equations
A Study of Numerical Methods for Simultaneous Equations Er. Chandan Krishna Mukherjee B.Sc.Engg., ME, MBA Asstt. Prof. ( Mechanical ), SSBT s College of Engg. & Tech., Jalgaon, Maharashtra Abstract: -
More informationModule 4F12: Computer Vision and Robotics Solutions to Examples Paper 2
Engineering Tripos Part IIB FOURTH YEAR Module 4F2: Computer Vision and Robotics Solutions to Examples Paper 2. Perspective projection and vanishing points (a) Consider a line in 3D space, defined in camera-centered
More informationInnite Impulse Response Filters Lecture by Prof Brian L Evans Slides by Niranjan Damera-Venkata Embedded Signal Processing Laboratory Dept of Electrical and Computer Engineering The University of Texas
More informationTHREE DIMENSIONAL GEOMETRY
For more important questions visit : www4onocom CHAPTER 11 THREE DIMENSIONAL GEOMETRY POINTS TO REMEMBER Distance between points P(x 1 ) and Q(x, y, z ) is PQ x x y y z z 1 1 1 (i) The coordinates of point
More informationGraded Assignment 2 Maple plots
Graded Assignment 2 Maple plots The Maple part of the assignment is to plot the graphs corresponding to the following problems. I ll note some syntax here to get you started see tutorials for more. Problem
More informationLinear, Quadratic, Exponential, and Absolute Value Functions
Linear, Quadratic, Exponential, and Absolute Value Functions Linear Quadratic Exponential Absolute Value Y = mx + b y = ax 2 + bx + c y = a b x y = x 1 What type of graph am I? 2 What can you tell me about
More informationHONOM EUROPEAN WORKSHOP on High Order Nonlinear Numerical Methods for Evolutionary PDEs, Bordeaux, March 18-22, 2013
HONOM 2013 EUROPEAN WORKSHOP on High Order Nonlinear Numerical Methods for Evolutionary PDEs, Bordeaux, March 18-22, 2013 A priori-based mesh adaptation for third-order accurate Euler simulation Alexandre
More informationSMBJ SERIES SMB SURFACE MOUNT UNIDIRECTIONAL AND BIDIRECTIONAL TRANSIENT VOLTAGE SUPPRESSORS
LITE-ON SURFACE MOUNT UNIDIRECTIONAL AND BIDIRECTIONAL TRANSIENT VOLTAGE SUPPRESSORS STAND-OFF VOLTAGE - 5.0 to 350 Volts POWER DISSIPATION - 600 WATTS FEATURES For surface mounted applications Reliable
More informationCentre for Autonomous Systems
Robot Henrik I Centre for Autonomous Systems Kungl Tekniska Högskolan hic@kth.se 27th April 2005 Outline 1 duction 2 Kinematic and Constraints 3 Mobile Robot 4 Mobile Robot 5 Beyond Basic 6 Kinematic 7
More informationThe Pre-Image Problem in Kernel Methods
The Pre-Image Problem in Kernel Methods James Kwok Ivor Tsang Department of Computer Science Hong Kong University of Science and Technology Hong Kong The Pre-Image Problem in Kernel Methods ICML-2003 1
More informationFreeMat Tutorial. 3x + 4y 2z = 5 2x 5y + z = 8 x x + 3y = -1 xx
1 of 9 FreeMat Tutorial FreeMat is a general purpose matrix calculator. It allows you to enter matrices and then perform operations on them in the same way you would write the operations on paper. This
More informationFast-Lipschitz Optimization
Fast-Lipschitz Optimization DREAM Seminar Series University of California at Berkeley September 11, 2012 Carlo Fischione ACCESS Linnaeus Center, Electrical Engineering KTH Royal Institute of Technology
More informationMATLAB Premier. Asst. Prof. Dr. Melik DÖLEN. Middle East Technical University Department of Mechanical Engineering 10/30/04 ME 304 1
MATLAB Premier Asst. Prof. Dr. Melik DÖLEN Middle East Technical University Department of Mechanical Engineering 0/0/04 ME 04 Outline! Introduction! Basic Features of MATLAB! Prompt Level and Basic Aritmetic
More informationModule 04 Block Diagrams and Graphical Representations of Intertwined Dynamic Systems
Module 04 Block Diagrams and Graphical Representations of Intertwined Dynamic Systems Ahmad F. Taha EE 3413: Analysis and Desgin of Control Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/
More informationSimultaneous equations 11 Row echelon form
1 Simultaneous equations 11 Row echelon form J A Rossiter For a neat organisation of all videos and resources http://controleducation.group.shef.ac.uk/indexwebbook.html Introduction 2 The previous videos
More informationUSING THE ORTHOGONAL PROJECTION FOR PARAMETER INITIALIZATION IN THE 3D RECONSTRUCTION OF DISTANT OBJECTS INTRODUCTION
USING THE ORTHOGONAL PROJECTION FOR PARAMETER INITIALIZATION IN THE 3D RECONSTRUCTION OF DISTANT OBJECTS Keith F. Blonquist, Junior Project Engineer Lidar Pacific Corporation Logan, UT 84341 kfblonquist@lidarpacific.com
More informationModeling, analysis and control of fuzzy systems
Electrical and Computer Engineering Faculty Publications Electrical & Computer Engineering 1995 Modeling, analysis and control of fuzzy systems Pushkin Kachroo University of Nevada, Las Vegas, pushkin@unlv.edu
More informationAdvanced Computer Graphics
G22.2274 001, Fall 2009 Advanced Computer Graphics Project details and tools 1 Project Topics Computer Animation Geometric Modeling Computational Photography Image processing 2 Optimization All projects
More informationFinite difference methods
Finite difference methods Siltanen/Railo/Kaarnioja Spring 8 Applications of matrix computations Applications of matrix computations Finite difference methods Spring 8 / Introduction Finite difference methods
More informationComputer Packet 1 Row Operations + Freemat
Computer Packet 1 Row Operations + Freemat For this packet, you will use a website to do row operations, and then learn to use a general purpose matrix calculator called FreeMat. To reach the row operations
More informationAn Improved Measurement Placement Algorithm for Network Observability
IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 16, NO. 4, NOVEMBER 2001 819 An Improved Measurement Placement Algorithm for Network Observability Bei Gou and Ali Abur, Senior Member, IEEE Abstract This paper
More informationMotion Tracking and Event Understanding in Video Sequences
Motion Tracking and Event Understanding in Video Sequences Isaac Cohen Elaine Kang, Jinman Kang Institute for Robotics and Intelligent Systems University of Southern California Los Angeles, CA Objectives!
More informationREPRESENTATION OF CONTROL SYSTEMS USING JAVA JAMES NEAL CONN, B.S. A THESIS MATHEMATICS
REPRESENTATION OF CONTROL SYSTEMS USING JAVA by JAMES NEAL CONN, B.S. A THESIS IN MATHEMATICS Submitted to the Graduate Faculty of Texas Tech University in Partial Fulfillment of the Requirements for the
More informationBasics of Multivariate Modelling and Data Analysis
Basics of Multivariate Modelling and Data Analysis Kurt-Erik Häggblom 9. Linear regression with latent variables 9.1 Principal component regression (PCR) 9.2 Partial least-squares regression (PLS) [ mostly
More informationInformation Retrieval. Lecture 11 - Link analysis
Information Retrieval Lecture 11 - Link analysis Seminar für Sprachwissenschaft International Studies in Computational Linguistics Wintersemester 2007 1/ 35 Introduction Link analysis: using hyperlinks
More informationAdvanced Computer Graphics
G22.2274 001, Fall 2010 Advanced Computer Graphics Project details and tools 1 Projects Details of each project are on the website under Projects Please review all the projects and come see me if you would
More informationFuzzy Neural Network-based Time Delay Prediction for Networked Control Systems
Appl. Math. Inf. Sci. 8, No. 1, 407-413 (2014) 407 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/080151 Fuzzy Neural Network-based Time Delay Prediction
More informationEstimation of Unknown Disturbances in Gimbal Systems
Estimation of Unknown Disturbances in Gimbal Systems Burak KÜRKÇÜ 1, a, Coşku KASNAKOĞLU 2, b 1 ASELSAN Inc., Ankara, Turkey 2 TOBB University of Economics and Technology, Ankara, Turkey a bkurkcu@aselsan.com.tr,
More informationMATLAB Premier. Middle East Technical University Department of Mechanical Engineering ME 304 1/50
MATLAB Premier Middle East Technical University Department of Mechanical Engineering ME 304 1/50 Outline Introduction Basic Features of MATLAB Prompt Level and Basic Arithmetic Operations Scalars, Vectors,
More informationSampling Approaches to Metrology in Semiconductor Manufacturing
Sampling Approaches to Metrology in Semiconductor Manufacturing Tyrone Vincent 1 and Broc Stirton 2, Kameshwar Poolla 3 1 Colorado School of Mines, Golden CO 2 GLOBALFOUNDRIES, Austin TX 3 University of
More informationFor example, the system. 22 may be represented by the augmented matrix
Matrix Solutions to Linear Systems A matrix is a rectangular array of elements. o An array is a systematic arrangement of numbers or symbols in rows and columns. Matrices (the plural of matrix) may be
More information15. Cutting plane and ellipsoid methods
EE 546, Univ of Washington, Spring 2012 15. Cutting plane and ellipsoid methods localization methods cutting-plane oracle examples of cutting plane methods ellipsoid method convergence proof inequality
More informationFinite Element Analysis Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology Madras. Module - 01 Lecture - 15
Finite Element Analysis Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology Madras Module - 01 Lecture - 15 In the last class we were looking at this 3-D space frames; let me summarize
More informationPick and Place Robot Simulation
Pick and Place Robot Simulation James Beukers Jordan Jacobson ECE 63 Fall 4 December 6, 4 Contents Introduction System Overview 3 3 State Space Model 3 4 Controller Design 6 5 Simulation and Results 7
More information