Advanced Control Systems Problem Sheet for Part B: Multivariable Systems

Similar documents
4.2 Continuous-Time Systems and Processes Problem Definition Let the state variable representation of a linear system be

A Simulative Study on Active Disturbance Rejection Control (ADRC) as a Control Tool for Practitioners

Introduction to Bayesian Estimation. McGill COMP 765 Sept 12 th, 2017

Probabilistic Robotics Sebastian Thrun-- Stanford

first-order circuit Complete response can be regarded as the superposition of zero-input response and zero-state response.

Localization and Map Making

Uncertainty & Localization I

ME 321: FLUID MECHANICS-I

DESIGN OF TENSION MEMBERS

The Research of Active Disturbance Rejection Control on Shunt Hybrid Active Power Filter

Lecture 8 Backlash and Quantization. Material. Linear and Angular Backlash. Example: Parallel Kinematic Robot. Backlash.

Scalar Conservation Laws

Optimal Control and Online Game Solutions Using Ui ADP.

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2

Vector autoregression VAR. Case 1

Redundancy Allocation of Partitioned Linear Block Codes

Earthquake, Volcano and Tsunami

Math 2214 Sol Test 2B Spring 2015

Mat 267 Engineering Calculus III Updated on 04/30/ x 4y 4z 8x 16y / 4 0. x y z x y. 4x 4y 4z 24x 16y 8z.

Study and Evaluation of Performances of the Digital Multimeter

A Mathematical model to Solve Reaction Diffusion Equation using Differential Transformation Method

I Let E(v! v 0 ) denote the event that v 0 is selected instead of v I The block error probability is the union of such events

1 First Order Partial Differential Equations

PH2130 Mathematical Methods Lab 3. z x

, u denotes uxt (,) and u. mean first partial derivatives of u with respect to x and t, respectively. Equation (1.1) can be simply written as

The co-movement of inflation and the real growth of output

Riemann Function and Methods of Group Analysis

CS 4495 Computer Vision Tracking 1- Kalman,Gaussian

10. State Space Methods

Data Fusion using Kalman Filter. Ioannis Rekleitis

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

The Euler-Lagrange Approach for Steady and Unsteady Flows. M. Sommerfeld. www-mvt.iw.uni-halle.de. Title. Zentrum für Ingenieurwissenschaften

Space truss bridge optimization by dynamic programming and linear programming

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES

Computer Vision. Motion Extraction

Homotopy Perturbation Method for Solving Partial Differential Equations

DISPLACEMENT ESTIMATION FOR IMAGE PREDICTIVE CODING AND FRAME MOTION-ADAPTIVE INTERPOLATION

ME 425: Aerodynamics

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Tracking. Many slides adapted from Kristen Grauman, Deva Ramanan

Srednicki Chapter 20

Description of the MS-Regress R package (Rmetrics)

Mixing of two different temperature water flows as feedback controlled system mathematically modeled and simulated in MS Excel spreadsheet

APPLIED ECONOMETRIC TIME SERIES (2nd edition)

Learning from a Golf Ball

A Direct Method for Solving Nonlinear PDEs and. New Exact Solutions for Some Examples

Computing with diode model

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Dirac s hole theory and the Pauli principle: clearing up the confusion.

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Miscellanea Miscellanea

UWE has obtained warranties from all depositors as to their title in the material deposited and as to their right to deposit such material.

Centre for Computational Finance and Economic Agents WP Working Paper Series. Hengxu Wang, John G. O Hara, Nick Constantinou

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

Numerical Simulation of Reaction-Diffusion Systems of Turing Pattern Formation

Nonlinear Network Structures for Optimal Control

Atmospheric Dynamics 11:670:324. Class Time: Tuesdays and Fridays 9:15-10:35

Optimal Time Varying Fuzzy Boundary Layer Thickness for Decouple Sliding-mode Control Base on Moving Switching Surface Using Pareto set.

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

Presentation Overview

A Review of Gradient Algorithms for Numerical Computation of Optimal Trajectories

Key Chemistry 102 Discussion #4, Chapter 11 and 12 Student name TA name Section. ; u= M. and T(red)=2*T(yellow) ; t(yellow)=4*t(red) or

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

Problemas das Aulas Práticas

Deterministic and Stochastic Optimal Control Analysis of an SIR Epidemic model

Relay Feedback Based Time Domain Modeling of Linear 3-by-3 MIMO System

ON JENSEN S INEQUALITY FOR g-expectation

Homework 6 AERE331 Spring 2019 Due 4/24(W) Name Sec. 1 / 2

Notes on Kalman Filtering

Applied Econometrics GARCH Models - Extensions. Roman Horvath Lecture 2

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

Chapter 3 Common Families of Distributions

Math Spring 2015 PRACTICE FINAL EXAM (modified from Math 2280 final exam, April 29, 2011)

Module 4: Time Response of discrete time systems Lecture Note 2

SYNTHESIS OF DYNAMICAL STRUSTURES FOR DESIGN OF VIBRO-SAFE HAND-HELD PERCUSSION MACHINES. V.I. Babitsky, I.J.Sokolov, N.A.

Compatible Versus Regular Well-Posed Linear Systems

Applied Mathematics Letters. Oscillation results for fourth-order nonlinear dynamic equations

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Vector autoregression VAR

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum

SOLVING AN OPTIMAL CONTROL PROBLEM WITH MATLAB

Linear Control System EE 711. Design. Lecture 8 Dr. Mostafa Abdel-geliel

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

ELEMENTS OF ACOUSTIC WAVES IN POROUS MEDIA

Module 4. Analysis of Statically Indeterminate Structures by the Direct Stiffness Method. Version 2 CE IIT, Kharagpur

Linear Quadratic Regulator (LQR) - State Feedback Design

An impact of noise on invariant manifolds in nonlinear dynamical systems

Basilio Bona ROBOTICA 03CFIOR 1

ON THE ANALYSIS OF SYMBOLIC DATA

Applications in Industry (Extended) Kalman Filter. Week Date Lecture Title

IB Physics Kinematics Worksheet

Optimal Stochastic Control in Continuous Time with Wiener Processes: General Results and Applications to Optimal Wildlife Management

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

Tracking. Many slides adapted from Kristen Grauman, Deva Ramanan

STAD57 Time Series Analysis. Lecture 14

Transcription:

436-45 Advanced Conrol Ssems Problem Shee for Par B: Mlivariable Ssems Qesion B 998 Given a lan o be conrolled, which is described b a sae-sace model A B C Oline he rocess b which o wold design a discree LQG reglaor incororaing a Kalman sae esimaor. Draw a simlaion diagram of he closed-loo ssem ha dislas he srcre of he reglaor. Wrie down he eqaions describing he Kalman esimaor, and describe how he wold be imlemened in a conrol comer. Qesion B 998 Under wha circmsances wold i be aroriae o se a fll-order esimaor, even hogh some of he lan saes are available for direc measremen? Derive he eqaions describing a redced-order esimaor. Show ha for an ideal redced-order esimaor erfec lan model, no rocess or measremen disrbances, he errors in he esimaed saes will aroach zero asmoicall. Qesion B3 998 Show how he reference ins ma be inrodced ino an LQG-reglaed ssem so ha i will have zero sead-sae errors in resonse o se-oin changes. Show how inegral conrol ma be inrodced ino an LQG-reglaed ssem so ha i will have zero sead-sae errors in resonse o se-oin changes and consan lan load disrbances. Qesion B4 998 Consider a lan whose ransfer fncion is /s and which has a ZOH on is in. If he o is samled wih a samling eriod s, he discree sae eqaions are.5 [ ] Wha is he discree ransfer fncion from o? Wrie down he smmeric roo locs eqaion for his ssem. Wha relaion does his have o i oimal conrol of he lan, and ii oimal esimaion of he saes? he SRL is drawn in figre. If oles are chosen for a sae-feedbac conroller which have roeries eqivalen o s-lane oles wih ω n.π rad/s and ζ.7, wha is he weigh laced on o ecrsions relaive o conrol effor in he erformance inde? Problem_A.doc

d e f Design a digial sae-feedbac conrol law ha will ield hese closed-loo oles. Design a fll-order esimaor wih dnamics five imes faser han he conroller. If he esimaor oles are chosen from he SRL, wha do he chosen oles iml abo he relaive levels of he assmed rocess and measremen noise? Draw a simlaion diagram of he closed-loo ssem, srcred so ha i will resond o se in commands wih zero sead-sae error..5 Smmeric roo locs.5 Im ag A is -.5 - Qesion A 999 -.5 -.5 - -.5.5.5 Real Ais Figre Given a lan o be conrolled, which is described b a sae-sace model A B C Oline he rocess b which o wold design a discree LQ reglaor incororaing a crren sae esimaor. Draw a simlaion diagram of he closed-loo ssem ha dislas he srcre of he reglaor. Wrie down he eqaions describing he crren esimaor, and describe how he wold be imlemened in a conrol comer. Problem_A.doc

Qesion A Describe he rocess b which disrbances of nown form, b nnown magnide, can be esimaed and comensaed for in a sae-sace conrol ssem design. A endlm wih a orqe moor a is hinge is described b [ ] b w where [ ] [ ] θ θ, θ being he angle of he endlm from verical, is a conrol orqe, w is a orqe bias, and b is a measremen bias. he o is samled wih a samle eriod. s, and he oal orqe w is alied hrogh a zero-order hold. Wih his discreisaion, he ssem and in marices become ξ ξ.998.5,.995.998.998.995 d e e B A A i Wih no measremen bias b, agmen he discree ssem model so ha he orqe bias is an elemen of he sae vecor. Is his ssem observable? ii Wih no orqe bias w, agmen he discree ssem model so ha he measremen bias is an elemen of he sae vecor. Is his ssem observable? Problem_A.doc 3

Qesion A3 Consider a doble-inegraor lan: s s U s Y Sech a smmeric roo locs SRL for his lan. Elain how his SRL cold be sed in he design of a feedbac reglaor. Wha is he heoreical basis for sch ses? A conroller canonical realisaion of he lan is [ ] Show ha he oimal conrol law which will minimise he erformance inde [ ]d J ρ is ρ ρ Hin: he erformance inde [ ]d J R Q is minimised b he conrol law K, where P B R K and P is a smmeric, osiive-definie solion of he coninos algebraic Riccai eqaion Q B B PBR PA P A For his simle lan i is sraighforward o manall solve he mari Riccai eqaion for P. he mari 3 P is osiive-definie if and onl if and >. 3 > Problem_A.doc 4

Qesion Simlaion diagrams for reglaed ssems incororaing i a discree redicion sae esimaor, and ii a discree crren sae esimaor, are shown in figres 4 and 5 resecivel. Show ha he esimaion error dnamics for he redicion esimaor are [ C] ~ ~ L where he redicion esimae error is defined as ~ ˆ. For a mli-in, mli-o MIMO ssem, discss some sraegies for deermining he esimaor gain L which will lace he esimaor oles a desired locaions in he z-lane. In a ole-lacemen design of a redicion esimaor for a MIMO ssem, he esimaor gain L can be fond wih he MALAB Conrol Ssems oolbo command L lacephi, C, de where he desired esimaor oles are defined in he vecor de. If he esimaion error of he crren esimaor is defined as ~ ˆ c derive an eqaion for he esimaion error dnamics, and hence show how he lace command can be sed o calclae he crren esimaor gain L c which will resl in he same esimaor oles, de. Qesion 3 Describe he fndamenal objecives which nderlie he mahemaical formlaion of an oimal conroller and an oimal esimaor. Describe he rocess of designing an LQG reglaor for a MIMO ssem. m q C m q C Plan Plan K ˆ q L ˆ q C ˆ Feedbac comensaor C ˆ Esimaor K Feedbac comensaor ˆ L c Esimaor Figre 4: Predicion esimaor Figre 5: Crren esimaor Problem_A.doc 5

Problem 7.58 [adaed from Franlin e al. ] Consider he endlm roblem wih conrol orqe c and disrbance orqe d : θ θ 4 c d here g/l 4. Assme here is a oeniomeer a he in ha measres he o angle θ, b wih a consan nnown bias b. hs he measremen eqaion is θ b. a b c d e f ae he agmened sae vecor o be θ where w is he in-eqivalen bias. Wrie he ssem eqaions in sae-sace form. Give vales for he marices A, B, and C. Using sae-variable mehods, show ha he characerisic eqaion of he model is ss 4. Show ha w is observable if we assme θ, and wrie he esimaor eqaions for ˆ ˆ θ θ wˆ. Pic esimaor gains [L L L 3 ] o lace all he roos of he esimaor-error characerisic eqaion a. Using fll sae feedbac of he esimaed conrollable sae-variables, derive a conrol law o lace he closed-loo oles a ± j. Draw a simlaion diagram of he comlee closed-loo ssem esimaor, lan, and conroller sing inegraor blocs. Inrodce he esimaed bias ino he conrol so as o ield zero sead-sae error o he o bias b. Demonsrae he erformance of or design b loing he resonse of he ssem o a se change in b; ha is, b changes from o some consan vale. θ w Problem_A.doc 6

Problem Problems reqiring MALAB An aroimae linear model for he laeral dnamics of an aircraf, wrien in erms of small errbaions from a rim condiion, is r β φ.7 9.7 r β φ.8 3.3 δ δ Sose a malfncion revens manilaion of he rdder in δ r. Is i ossible o conrol he aircraf sing onl he aileron in δ a? If no, wha aircraf moions can and canno be conrolled? Is he aircraf conrollable wih js he rdder in δ r? d Verif ha he aircraf is conrollable wih boh ins oerable. If he onl o is a measremen of he roll rae rovided b a rae gro, is he ssem observable? e Is he ssem observable from ban indicaor measremens of φ? f g h Design a discree sae-feedbac conroller, wih a samle eriod of s. s, wih closed-loo oles which have roeries corresonding o s-lane oles a s, 5, 8,. Show how reference ins for aw rae and roll angle can be inrodced sch ha he resonse o se commands will ehibi zero sead-sae errors. Design a discree LQ oimal reglaor, wih a samle eriod of s. s. Sar wih he following rial weighing marices: Q diag.,.,.,., R I. Se he closed loo ssem o rac aw and roll commands. Problem An nsable discree-ime ssem is described b [ ] 3 a r Find he infinie horizon feedbac gain mari for he oimal reglaor roblem, assming all saes are available for se. Also deermine he resling closed-loo eigenvales. Use he following weighs: i Q I, R ii Q I, R iii Q I, R iv Q I, R Problem_A.doc 7

he ssem is now assmed o have addiive rocess noise w and measremen noise v. he noise covariance marices are: 9 Q and 6, resecivel 4 R A Kalman filer is o be sed raher han a deerminisic observer. Find he seadsae Kalman gain K and he oles of he filer s dnamics. he above filer is sed wih he LQ conroller gains fond above. How do he resence of he filer and he se of esimaed saes raher han acal saes affec he overall ssem erformance in each case? he searaion heorem garanees ha his aroach is he bes ha can be done. However, i does no garanee ha he resls will be wha he conroller ased for, or ha he ssem will even wor well. d Reea wih he measremen noise covariance redced o R. Problem 3 For a lan described b:. 3..5......5 8.9.5 4. [...] design a linear/qadraic conroller ha minimises he cos fncion: q r J d Perform he design for differen vales of q and r. Chec he sabili of he closed-loo ssem wih a linear/qadraic conroller for all design cases. Obain he oen-loo ransfer fncion o calclae he gain and hase margins. Comare he resls for differen vales of q and r. Comare he sensiivi of ole and zero locaions and se resonses of he closedloo ssem, for changing vales of q and r. Commen on he robsness of LQ-designed conrollers Problem_A.doc 8