F This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.

Similar documents
23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Lecture 20: Riccati Equations and Least Squares Feedback Control

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

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

The general Solow model

2 int T. is the Fourier transform of f(t) which is the inverse Fourier transform of f. i t e

4. The multiple use forestry maximum principle. This principle will be derived as in Heaps (1984) by considering perturbations

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

THE SINE INTEGRAL. x dt t

Sterilization D Values

t dt t SCLP Bellman (1953) CLP (Dantzig, Tyndall, Grinold, Perold, Anstreicher 60's-80's) Anderson (1978) SCLP

Sliding Mode Controller for Unstable Systems

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

Lectures 29 and 30 BIQUADRATICS AND STATE SPACE OP AMP REALIZATIONS. I. Introduction

Chapter 2. First Order Scalar Equations

Differential Equations

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

Problem Set #3: AK models

Solutions to Assignment 1

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Chapter 4. Truncation Errors

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

Testing the Random Walk Model. i.i.d. ( ) r

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

10. State Space Methods

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

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

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Chapter 8 The Complete Response of RL and RC Circuits

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Mathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

Hamilton Jacobi equations

Linear Response Theory: The connection between QFT and experiments

Continuous Time Linear Time Invariant (LTI) Systems. Dr. Ali Hussein Muqaibel. Introduction

ME 391 Mechanical Engineering Analysis

on the interval (x + 1) 0! x < ", where x represents feet from the first fence post. How many square feet of fence had to be painted?

Linear Time-invariant systems, Convolution, and Cross-correlation

6.003 Homework #9 Solutions

Notes 04 largely plagiarized by %khc

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Two Coupled Oscillators / Normal Modes

Math 334 Fall 2011 Homework 11 Solutions

Convolution. Lecture #6 2CT.3 8. BME 333 Biomedical Signals and Systems - J.Schesser

Circuit Variables. AP 1.1 Use a product of ratios to convert two-thirds the speed of light from meters per second to miles per second: 1 ft 12 in

Notes for Lecture 17-18

ME 452 Fourier Series and Fourier Transform

i L = VT L (16.34) 918a i D v OUT i L v C V - S 1 FIGURE A switched power supply circuit with diode and a switch.

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

System of Linear Differential Equations

Second Law. first draft 9/23/04, second Sept Oct 2005 minor changes 2006, used spell check, expanded example

6.003 Homework #9 Solutions

Math 106: Review for Final Exam, Part II. (x x 0 ) 2 = !

Expert Advice for Amateurs

Chapter 3 Boundary Value Problem

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Lecture 2: Optics / C2: Quantum Information and Laser Science

Math From Scratch Lesson 34: Isolating Variables

Q1) [20 points] answer for the following questions (ON THIS SHEET):

1 Evaluating Chromatograms

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

SZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1

The Laplace Transform

6.2 Transforms of Derivatives and Integrals.

6.003 Homework #8 Solutions

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

EE 301 Lab 2 Convolution

BU Macro BU Macro Fall 2008, Lecture 4

Chapter 6. Systems of First Order Linear Differential Equations

Ordinary dierential equations

Let us start with a two dimensional case. We consider a vector ( x,

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Numerical Dispersion

EECE 301 Signals & Systems Prof. Mark Fowler

6.302 Feedback Systems Recitation 4: Complex Variables and the s-plane Prof. Joel L. Dawson

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

Math Final Exam Solutions

Vehicle Arrival Models : Headway

Economics 6130 Cornell University Fall 2016 Macroeconomics, I - Part 2

Theory of! Partial Differential Equations!

DIFFERENTIAL GEOMETRY HW 5

RC, RL and RLC circuits

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Physical Limitations of Logic Gates Week 10a

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

Power of Random Processes 1/40

SOLUTIONS TO ECE 3084

-e x ( 0!x+1! ) -e x 0!x 2 +1!x+2! e t dt, the following expressions hold. t

Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis

Final Spring 2007

e 2t u(t) e 2t u(t) =?

Transcription:

Lecure 8 Las ime: Semi-free configuraion design This is equivalen o: Noe ns, ener he sysem a he same place. is fixed. We design C (and perhaps B. We mus sabilize if i is given as unsable. Cs ( H( s = + CssBs ( ( ( so ha having he opimum H, we deermine C from H( s Cs ( = H ( s ( s B ( s We do no collec H and ogeher because if is non-minimum phase, we would no wish o define H by ( H H = op This leads o an unsable mode which is no observable a he oupu hus canno be conrolled by feeding back. Associae weighing funcions wih he given ransfer funcions. H( s wh ( s ( w ( Ds ( w( D Page of 5

If s ( is unsable, pu a sabilizing feedback around i, laer associae i wih he res of he sysem. Error Analysis We require he mean squared error. c ( = w ( ( i d H o ( = w ( ( c d = d w ( d w ( i( H d ( = w( ( s d D e ( = o ( d ( 3 3 3 e ( = o ( od ( ( + d ( o( = d w( dwh( ( i d4w( 4 d3wh( 3( i 3 4 = d w ( d w ( d w ( d w ( i( i( H 3 H 3 4 4 3 4 = d w ( d w ( d w ( d4w( 4 Rii( + 3 4 H 3 H 3 od ( ( = d w( dwh( ( i d3wd( 3( s 3 = d w ( d w ( d w ( i( s( H 3 D 3 3 = d w ( d w ( d w ( R ( + H 3 D 3 is 3 We shall no require d ( in inegral form. Page of 5

The problem now is o choose wh ( so as o minimize his use variaional calculus. Le: w ( = w ( + δ w( H e (, for which we where w ( is he opimum weighing funcion (o be deermined and δ w ( is an arbirary variaion arbirary excep ha i mus be physically realizable. Calculae he opimum e and is firs and second variaions. e = e + δe + δ e e = o( + o( d( + d( The opimum e ( e for δ w ( = : = 3 3 4 4 ii + 3 4 dw( dw( d3wd( 3 Ris( + 3 + d( e( d w ( d w ( d w ( d w ( R ( The firs variaion in e ( is ( = ( ( 3 ( 3 4 ( 4 ii( + 3 4 δ e dδw d w d w d w R + d w ( d w ( d δw( d w ( R ( + 3 3 4 4 ii 3 4 dδw( d w ( d w ( R ( + 3 D 3 is 3 In he second erm, le: = 3 = 4 3 = 4 = and inerchange he order of inegraion. nd erm = d δw( d w( d 3w ( 3 d 4w( 4 Rii( 3+ 4 Page 3 of 5

bu since Rii( 3+ 4 = Rii( + 3 4 we see ha he second erm is exacly equal o he firs erm. Collecing hese erms and separaing ou he common inegral wih respec o gives δ e( = dδ w( dw( d3w( 3 d4w( 4 Rii( + 3 4 dw( d3wd( 3 Ris( + 3 The second variaion of e ( is δ e = dδ w dw dδ 3 w 3 d4w 4 Rii + 3 4 ( ( ( ( ( ( By comparison wih he expression for oupu of he sysem o (, his is seen o be he mean squared ( δ e oupu = ( >, non-zero inpu This second variaion mus be greaer han zero, so he saionary poin defined by he vanishing of he firs variaion is shown o be a minimum. In he expression for he firs variaion, δ w( = for < by he requiremen ha he variaion be physically realizable. Bu δ w( is arbirary for, so we can be assured of he vanishing of δ e ( only if he { } erm vanishes almos everywhere for. The condiion which defines he minimum in e ( is hen d w ( d w ( d w ( R ( + 3 3 4 4 ii 3 4 d w ( d w ( R ( + = 3 D 3 is 3 for all, non-real-ime. Using his condiion in he expression for for < gives he resul e ( and remembering ha w ( = Page 4 of 5

e ( = d ( o ( which is convenien for he calculaion of e (. Also since o ( = d ( e (, his says he opimum mean squared oupu is always less han he mean squared desired oupu. Auocorrelaion uncions We have arrived a an exended form of he Wiener-Kopf equaion which defines he opimum linear sysem under he ground rules saed before. Recall ha: Rii( = Rss ( + Rsn( + Rns ( + Rnn ( Ris( = Rss ( + Rns ( since i = s+ n. The free configuraion problem is a specializaion of he semi-free configuraion. In his expression we would ake s= (, or w ( = δ (. In ha case we have d d( d w ( d δ( R ( + 3 3 4 4 ii 3 4 dδ ( dw ( R ( + = 3 D 3 is 3 w ( R ( d w ( R ( d = for 3 ii 3 3 D 3 is 3 3 Page 5 of 5