Servomechanism Design

Size: px
Start display at page:

Download "Servomechanism Design"

Transcription

1 Sevomechanism Design Sevomechanism (sevo-sysem) is a conol sysem in which he efeence () (age, Se poin) changes as ime passes. Design mehods PID Conol u () Ke P () + K I ed () + KDe () Sae Feedback u() Kx() + S K x ( ) () Gain K is obained based on Pole Assignmen Opimal egulao design Sevo Seve / Slave Refeence sysem () Ax (), x () x y () Cx()

2 Refeence Sysem x () A x () x () x () c () Linea sysem which geneaes he efeence signal y () 3 Ex. Sep signal (Consan) () x( ) 3 () x (), x() 3 y( ) [ ] 3 y [ ] () x() { x2, x() 2, x2() Ex.2 Ramp signal (Polynomial of ) () x (), x() x() d + + 2( ) x y [ ] + () x() y() [ ] + Ex.3 igonomeic funcions x2, 2 x () (), () x2 x2, x2() x x x 2() cos, x() sin sin y [ ] () x() y () [ ] cos cos y x,

3 Sevo Poblem Obain a conol inpu u() so ha he sysem oupu y() coincides wih he efeence y (). () () () A x, x x, x R y () () () () Ax + bu, x x, x R y c () x() c () x () Deemine a conol inpu u*() so ha he sysem oupu y() y (),. y* () y () u*(): ideal inpu, x*(): ideal sae, y*(): ideal oupu *() Ax* ()+ bu* (), x* () x () A x(), x() x x* () c x () c n n

4 Goal y*( ) () Inenal Model Pinciple x () x() x* () x () y c * c cx( ) c x( ) x * x () () A * A () x + * () bu () A () + bu* () bu *( ) A A x () * ( ) ( ) c bs x () A A x () Inenal Model Pinciple A c A c b s x c Linea elaion u () s x () s Linea elaion b A A and s obained fom Inenal Model Pinciple give he ideal conol inpu u* ha achieves he Goal y*()y ()

5 Sevomechanism design based on Sae Feedback Conol x() x () : impossible in geneal x () x () x() x () x() * * ( ) u () * ε u () u() : diffeence beween ideal conol and eal conol x () * ε x () x () : diffeence beween ideal sae and eal sae Eo sysem () Ax () + bu (), x () x x ε ε ε ε Deemine uε so ha xε and uε conveges o zeo as fas as possible. Design Sae Feedback uεk xε fo Pole Assignmen fo A+b k Possible if (A,b) is conollable

6 Sevomechanism design based on Sae Feedback Conol u () k x () : Sae feedback law fo eo sysem ε ε uε u u * () () (), xε x x ( ) * * u u k x x () () () () * u s () x (), x * * () () () () x () coninued ( ) u () k x() + s k x () Feed fowad of efeence sysem sae Sae feedback fo oiginal sysem Obseve design

7 Block Diagam of Sae-feedback-based Sevo Mechanism ( s ) x u () k x() + k () Esimaed by Sae Obseve

8 Design Example Design a sevomechanism o follow he efeence wih convegence speed specified wih poles {-2,-3} fo he linea sysem: Seps fo soluion () () (), () 3 2 x + u x y () x() [ ] y () cos. Make a efeence sysem o geneae he efeence signal y () 2. Calculae and s fo he esulan sevo-sysem o saisfy he Inenal Model Pinciple A A bs, c c 3. Calculae he sae feedback gain k o ealize he specified convegence speed by pole assignmen 4. Simulae he behavio of he esulan sevomechanism (Make an augmened sysem fo () x() x () ) [ ]

9 Design example - sep () x (), x() y () x () [ ] x( ) sin x () cos 2 (Alenaive) x () y () cos, x () () sin 2 x () sin () () x, x 2() cos x2 () x2() () x (), x() y () x () [ ] x x 2 ( ) () cos sin

10 Design example - sep 2/3 xy, v zu w s A A, Sep 2 Pu s b c c ha is, And solve he pinciple: x y x y x y z u 3 2 z u z u 2, s 2 Sep 3 Oiginal chaaceisic polynomial is s s s s s I [ v w], [ ] [ ] Fom he desiable poles {-2,-3}, he ideal chaaceisic polynomial is 2 ( s+ 2)( s+ 3) s + 5s+ 6 Hence he sae feedback gain k is obained as k [ k k ] [ 3 2 ] [ 3 3] 6 5

11 Design example - sep 4 ( ) Ax( ) + bu (), x() x, x R, y( ) x( ) A x x x x R () n () (), (),, y() c x ) k x( ) + x( ) u ( ( s k ) k : deemined so ha x ε x*-x quickly conveges o () Ax() + b k x() + k x () { ( s ) } ( ) ( ) ( A bk x + b s k ) x () + ( ) Ox() + A x ( ) ( s ) x () A+ bk b k x() x(), () () () O A x x x y () c x() y() () c x n () Ax (), (), x R 2 y () Cx (), y R c n+ n

12 Design example Design/Simulaion pogam /3 H Pocess L n 2; n 2; m ; A 88, <, 8 3, 2<<; b 88<, 8<<; c 88, <<; x 8, <; AJ 3 2 N bj N c^h L x_j N A_J N c_^h L x_j N H Refeence geneao L A 88, <, 8, <<; c 88, <<; x 8, <; Pin@"A", A êê MaixFom, " b", b êê MaixFom, " c^", c êê MaixFom, " x_", x êê MaixFomD Pin@"A_", A êê MaixFom, " c_^", c êê MaixFom, " x_", x êê MaixFomD

13 Design example Design/Simulaion pogam 2/3 H Inenal Model Pinciple L 88, 2<, 82, 22<<; s 88s, s2<<; Eq 8.A A. b.s, c. c<; Sol Flaen@Solve@Flaen@EqD, Flaen@8, s<ddd; ê. Sol; J s s ê. Sol; Pin@"", êê MaixFom, " s^", s êê MaixFomD N s^h 2 2L H Pole Assignmen fo Eo Sysem L ChCoef Flaen@ake@A, DD; H Coefficiens of chaceisic equaion L OpPole 8 2, 3<; OpChPoly ; H Opimal chaaceisic polynomial L Do@OpChPoly OpChPoly Hs OpPole@@iDDL, 8i,, n<d OpCoef 8<; H Opimal coefficiens L Do@OpCoef Append@OpCoef, Coefficien@OpChPoly, s, idd, 8i,, n <D Feedback gain fo σ^ 8 2, 3< is k^h 3 k 8ChCoef OpCoef<; Pin@"Feedback gain fo σ^ ", Ldes, " is k^", k êê MaixFomD 3 L

14 Design example Design/Simulaion pogam 3/3 H Make Simulaion Sysem HAugmened SysemL L M Join@A + b.k, able@, 8i, 2<, 8j, 2<DD; M2 Join@b.Hs k.l, AD; AA anspose@join@anspose@md, anspose@m2ddd; Pin@"AA", AA êê MaixFom, " CC", CC êê MaixFomD Pin@"Eigenvalues of AA", Eigenvalues@AADD M Join@c, able@, 8i, <, 8j, 2<DD; M2 Join@able@, 8i, <, 8j, 2<D, cd; CC anspose@join@anspose@md, anspose@m2ddd; H Rewie fom SaeEquaion L i n 4; vx able@x@id@d, 8i, n<d; AA vx Join@x, xd; j k f 3 Pi; Equ HAA.vxLPi, 8i,, n<d; IC able@x@id@d vxpi, 8i,, n<d; DE Join@Equ, ICD; Sol NDSolve@DE, vx, 8,, f<d; y z { CCJ N Plo@Evaluae@CC.vx ê. SolD, 8,, f<, PloRange All, Fame ue, FameLabel "yhl & yhl"<, DefaulFon > 8"Helveica", 5<, GidLines Auomaic, PloSyle 88RGBColo@,, D, hickness@.3d<, 8RGBColo@,, D, hickness@.7d<<d

15 Refeence & Oupu.5 y HL &yhl

16 Saes of efeence sysem & plan.5 x HL &xhl

17 Response fo anohe iniial sae AJ 3 2 N bj N c^h L x_j N y HL &yhl x HL &xhl

18 Pole assignmen fo fase convegence Feedback gain fo σ^ 8 5, 6< is k^h 27 i y 3 29 A j z k { CCJ N 9 L y HL &yhl x HL &xhl

19 Pole assignmen fo slowe convegence Feedback gain fo σ^ 8.5,.6< is k^h L i.... y AA j z k { CCJ N y HL &yhl x HL &xhl

20 Response fo fase efeence A_J 2 2 N c_^h L x_j i 2 j N y z k s^i 2 { 2 M y HL &yhl x HL &xhl

21 Design Pacice Design a sevomechanism o follow he efeence wih convegence speed specified wih poles {-4±2j} fo he linea sysem: () () (), () 3 2 x + u x y () x() [ ] y () sin2

On Control Problem Described by Infinite System of First-Order Differential Equations

On Control Problem Described by Infinite System of First-Order Differential Equations Ausalian Jounal of Basic and Applied Sciences 5(): 736-74 ISS 99-878 On Conol Poblem Descibed by Infinie Sysem of Fis-Ode Diffeenial Equaions Gafujan Ibagimov and Abbas Badaaya J'afau Insiue fo Mahemaical

More information

LAPLACE TRANSFORM AND TRANSFER FUNCTION

LAPLACE TRANSFORM AND TRANSFER FUNCTION CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions

More information

Intermediate Differential Equations Review and Basic Ideas

Intermediate Differential Equations Review and Basic Ideas Inermediae Differenial Equaions Review and Basic Ideas John A. Burns Cener for Opimal Design And Conrol Inerdisciplinary Cener forappliedmahemaics Virginia Polyechnic Insiue and Sae Universiy Blacksburg,

More information

Risk tolerance and optimal portfolio choice

Risk tolerance and optimal portfolio choice Risk oleance and opimal pofolio choice Maek Musiela BNP Paibas London Copoae and Invesmen Join wok wih T. Zaiphopoulou (UT usin) Invesmens and fowad uiliies Pepin 6 Backwad and fowad dynamic uiliies and

More information

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING MEEN 67 Handou # MODAL ANALYSIS OF MDOF Sysems wih VISCOS DAMPING ^ Symmeic Moion of a n-dof linea sysem is descibed by he second ode diffeenial equaions M+C+K=F whee () and F () ae n ows vecos of displacemens

More information

Lecture 1: Contents of the course. Advanced Digital Control. IT tools CCSDEMO

Lecture 1: Contents of the course. Advanced Digital Control. IT tools CCSDEMO Goals of he course Lecure : Advanced Digial Conrol To beer undersand discree-ime sysems To beer undersand compuer-conrolled sysems u k u( ) u( ) Hold u k D-A Process Compuer y( ) A-D y ( ) Sampler y k

More information

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

Q1) [20 points] answer for the following questions (ON THIS SHEET): Dr. Anas Al Tarabsheh The Hashemie Universiy Elecrical and Compuer Engineering Deparmen (Makeup Exam) Signals and Sysems Firs Semeser 011/01 Final Exam Dae: 1/06/01 Exam Duraion: hours Noe: means convoluion

More information

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

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information

Linear Quadratic Regulator (LQR) - State Feedback Design

Linear Quadratic Regulator (LQR) - State Feedback Design Linear Quadrai Regulaor (LQR) - Sae Feedbak Design A sysem is expressed in sae variable form as x = Ax + Bu n m wih x( ) R, u( ) R and he iniial ondiion x() = x A he sabilizaion problem using sae variable

More information

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,

More information

EE363 homework 1 solutions

EE363 homework 1 solutions EE363 Prof. S. Boyd EE363 homework 1 soluions 1. LQR for a riple accumulaor. We consider he sysem x +1 = Ax + Bu, y = Cx, wih 1 1 A = 1 1, B =, C = [ 1 ]. 1 1 This sysem has ransfer funcion H(z) = (z 1)

More information

System of Linear Differential Equations

System of Linear Differential Equations Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y

More information

Chapter 1 Fundamental Concepts

Chapter 1 Fundamental Concepts Chaper 1 Fundamenal Conceps 1 Signals A signal is a paern of variaion of a physical quaniy, ofen as a funcion of ime (bu also space, disance, posiion, ec). These quaniies are usually he independen variables

More information

Non-sinusoidal Signal Generators

Non-sinusoidal Signal Generators Non-sinusoidal Signal Geneaos ecangle, iangle, saw ooh, pulse, ec. Muliibao cicuis: asable no sable saes (wo quasi-sable saes; i emains in each sae fo pedeemined imes) monosable one sable sae, one unsable

More information

non-linear oscillators

non-linear oscillators non-linear oscillaors The invering comparaor operaion can be summarized as When he inpu is low, he oupu is high. When he inpu is high, he oupu is low. R b V REF R a and are given by he expressions derived

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

Differential Equations

Differential Equations Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding

More information

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions Inenaional Mahemaical Foum, Vol 8, 03, no 0, 463-47 HIKARI Ld, wwwm-hikaicom Combinaoial Appoach o M/M/ Queues Using Hypegeomeic Funcions Jagdish Saan and Kamal Nain Depamen of Saisics, Univesiy of Delhi,

More information

156 There are 9 books stacked on a shelf. The thickness of each book is either 1 inch or 2

156 There are 9 books stacked on a shelf. The thickness of each book is either 1 inch or 2 156 Thee ae 9 books sacked on a shelf. The hickness of each book is eihe 1 inch o 2 F inches. The heigh of he sack of 9 books is 14 inches. Which sysem of equaions can be used o deemine x, he numbe of

More information

KINEMATICS OF RIGID BODIES

KINEMATICS OF RIGID BODIES KINEMTICS OF RIGID ODIES In igid body kinemaics, we use he elaionships govening he displacemen, velociy and acceleaion, bu mus also accoun fo he oaional moion of he body. Descipion of he moion of igid

More information

DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2008

DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2008 [E5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 008 EEE/ISE PART II MEng BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: :00 hours There are FOUR quesions

More information

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004 Augmened Realiy II Kalman Filers Gudrun Klinker May 25, 2004 Ouline Moivaion Discree Kalman Filer Modeled Process Compuing Model Parameers Algorihm Exended Kalman Filer Kalman Filer for Sensor Fusion Lieraure

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay) Secions 3.1 and 3.4 Eponenial Funcions (Gowh and Decay) Chape 3. Secions 1 and 4 Page 1 of 5 Wha Would You Rahe Have... $1million, o double you money evey day fo 31 days saing wih 1cen? Day Cens Day Cens

More information

r P + '% 2 r v(r) End pressures P 1 (high) and P 2 (low) P 1 , which must be independent of z, so # dz dz = P 2 " P 1 = " #P L L,

r P + '% 2 r v(r) End pressures P 1 (high) and P 2 (low) P 1 , which must be independent of z, so # dz dz = P 2  P 1 =  #P L L, Lecue 36 Pipe Flow and Low-eynolds numbe hydodynamics 36.1 eading fo Lecues 34-35: PKT Chape 12. Will y fo Monday?: new daa shee and daf fomula shee fo final exam. Ou saing poin fo hydodynamics ae wo equaions:

More information

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

e 2t u(t) e 2t u(t) =? EE : Signals, Sysems, and Transforms Fall 7. Skech he convoluion of he following wo signals. Tes No noes, closed book. f() Show your work. Simplify your answers. g(). Using he convoluion inegral, find

More information

Reinforcement learning

Reinforcement learning Lecue 3 Reinfocemen leaning Milos Hauskech milos@cs.pi.edu 539 Senno Squae Reinfocemen leaning We wan o lean he conol policy: : X A We see examples of x (bu oupus a ae no given) Insead of a we ge a feedback

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 20: Riccati Equations and Least Squares Feedback Control 34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he

More information

Signal and System (Chapter 3. Continuous-Time Systems)

Signal and System (Chapter 3. Continuous-Time Systems) Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b

More information

Chapter 3 Boundary Value Problem

Chapter 3 Boundary Value Problem Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le

More information

Control Volume Derivation

Control Volume Derivation School of eospace Engineeing Conol Volume -1 Copyigh 1 by Jey M. Seizman. ll ighs esee. Conol Volume Deiaion How o cone ou elaionships fo a close sysem (conol mass) o an open sysem (conol olume) Fo mass

More information

Differential Geometry: Numerical Integration and Surface Flow

Differential Geometry: Numerical Integration and Surface Flow Differenial Geomery: Numerical Inegraion and Surface Flow [Implici Fairing of Irregular Meshes using Diffusion and Curaure Flow. Desbrun e al., 1999] Energy Minimizaion Recall: We hae been considering

More information

ECE Spring Prof. David R. Jackson ECE Dept. Notes 20

ECE Spring Prof. David R. Jackson ECE Dept. Notes 20 ECE 6345 Spring 5 Prof. David R. Jacson ECE Dep. Noes Overview In his se of noes we apply he SDI mehod o invesigae he fields produced by a pach curren. We calculae he field due o a recangular pach on op

More information

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

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Lecture 22 Electromagnetic Waves

Lecture 22 Electromagnetic Waves Lecue Elecomagneic Waves Pogam: 1. Enegy caied by he wave (Poyning veco).. Maxwell s equaions and Bounday condiions a inefaces. 3. Maeials boundaies: eflecion and efacion. Snell s Law. Quesions you should

More information

ECE 2100 Circuit Analysis

ECE 2100 Circuit Analysis ECE 1 Circui Analysis Lesson 35 Chaper 8: Second Order Circuis Daniel M. Liynski, Ph.D. ECE 1 Circui Analysis Lesson 3-34 Chaper 7: Firs Order Circuis (Naural response RC & RL circuis, Singulariy funcions,

More information

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s MÜHENDİSLİK MEKANİĞİ. HAFTA İMPULS- MMENTUM-ÇARPIŞMA Linea oenu of a paicle: The sybol L denoes he linea oenu and is defined as he ass ies he elociy of a paicle. L ÖRNEK : THE LINEAR IMPULSE-MMENTUM RELATIN

More information

ENGI 9420 Engineering Analysis Assignment 2 Solutions

ENGI 9420 Engineering Analysis Assignment 2 Solutions ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion

More information

EECS20n, Solution to Midterm 2, 11/17/00

EECS20n, Solution to Midterm 2, 11/17/00 EECS20n, Soluion o Miderm 2, /7/00. 0 poins Wrie he following in Caresian coordinaes (i.e. in he form x + jy) (a) 2 poins j 3 j 2 + j += j ++j +=2 (b) 2 poins ( j)/( + j) = j (c) 2 poins cos π/4+jsin π/4

More information

The Global Trade and Environment Model: GTEM

The Global Trade and Environment Model: GTEM The Global Tade and Envionmen Model: A pojecion of non-seady sae daa using Ineempoal GTEM Hom Pan, Vivek Tulpulé and Bian S. Fishe Ausalian Bueau of Agiculual and Resouce Economics OBJECTIVES Deive an

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 31 Signals & Sysems Prof. Mark Fowler Noe Se #1 C-T Sysems: Convoluion Represenaion Reading Assignmen: Secion 2.6 of Kamen and Heck 1/11 Course Flow Diagram The arrows here show concepual flow beween

More information

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations Today - Lecue 13 Today s lecue coninue wih oaions, oque, Noe ha chapes 11, 1, 13 all inole oaions slide 1 eiew Roaions Chapes 11 & 1 Viewed fom aboe (+z) Roaional, o angula elociy, gies angenial elociy

More information

Low-complexity Algorithms for MIMO Multiplexing Systems

Low-complexity Algorithms for MIMO Multiplexing Systems Low-complexiy Algoihms fo MIMO Muliplexing Sysems Ouline Inoducion QRD-M M algoihm Algoihm I: : o educe he numbe of suviving pahs. Algoihm II: : o educe he numbe of candidaes fo each ansmied signal. :

More information

Homework 4 SOLUTION EE235, Summer 2012

Homework 4 SOLUTION EE235, Summer 2012 Homework 4 SOLUTION EE235, Summer 202. Causal and Sable. These are impulse responses for LTI sysems. Which of hese LTI sysem impulse responses represen BIBO sable sysems? Which sysems are causal? (a) h()

More information

Stochastic Model Predictive Control for Gust Alleviation during Aircraft Carrier Landing

Stochastic Model Predictive Control for Gust Alleviation during Aircraft Carrier Landing Ouline Moivaion Sochasic formulaion Sochasic Model Predicive Conrol for Gus Alleviaion during Aircraf Carrier Landing Aircraf and gus modeling Resuls Deparmen of Mechanical and Aerospace Engineering Rugers,

More information

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

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by

More information

Module MA1132 (Frolov), Advanced Calculus Tutorial Sheet 1. To be solved during the tutorial session Thursday/Friday, 21/22 January 2016

Module MA1132 (Frolov), Advanced Calculus Tutorial Sheet 1. To be solved during the tutorial session Thursday/Friday, 21/22 January 2016 Module MA113 (Frolov), Advanced Calculus Tuorial Shee 1 To be solved during he uorial session Thursday/Friday, 1/ January 16 A curve C in he xy-plane is represened by he equaion Ax + Bxy + Cy + Dx + Ey

More information

Lecture 17: Kinetics of Phase Growth in a Two-component System:

Lecture 17: Kinetics of Phase Growth in a Two-component System: Lecue 17: Kineics of Phase Gowh in a Two-componen Sysem: descipion of diffusion flux acoss he α/ ineface Today s opics Majo asks of oday s Lecue: how o deive he diffusion flux of aoms. Once an incipien

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

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.

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. . Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.

More information

Sliding Mode Controller for Unstable Systems

Sliding Mode Controller for Unstable Systems S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.

More information

Solutions to Assignment 1

Solutions to Assignment 1 MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Anti-Disturbance Control for Multiple Disturbances

Anti-Disturbance Control for Multiple Disturbances Workshop a 3 ACC Ani-Disurbance Conrol for Muliple Disurbances Lei Guo (lguo@buaa.edu.cn) Naional Key Laboraory on Science and Technology on Aircraf Conrol, Beihang Universiy, Beijing, 9, P.R. China. Presened

More information

6.003 Homework #9 Solutions

6.003 Homework #9 Solutions 6.003 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 3 0 a 0 5 a k a k 0 πk j3 e 0 e j πk 0 jπk πk e 0

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

7 Wave Equation in Higher Dimensions

7 Wave Equation in Higher Dimensions 7 Wave Equaion in Highe Dimensions We now conside he iniial-value poblem fo he wave equaion in n dimensions, u c u x R n u(x, φ(x u (x, ψ(x whee u n i u x i x i. (7. 7. Mehod of Spheical Means Ref: Evans,

More information

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

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB EECS 6 Winer 7 Laboraory 3 Fourier series, Fourier ransform and Bode Plos in MATLAB Inroducion: The objecives of his lab are o use MATLAB:. To plo periodic signals wih Fourier series represenaion. To obain

More information

6.003 Homework #8 Solutions

6.003 Homework #8 Solutions 6.003 Homework #8 Soluions Problems. Fourier Series Deermine he Fourier series coefficiens a k for x () shown below. x ()= x ( + 0) 0 a 0 = 0 a k = e /0 sin(/0) for k 0 a k = π x()e k d = 0 0 π e 0 k d

More information

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP). Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions

Signals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions 8-90 Signals and Sysems Profs. Byron Yu and Pulki Grover Fall 07 Miderm Soluions Name: Andrew ID: Problem Score Max 0 8 4 6 5 0 6 0 7 8 9 0 6 Toal 00 Miderm Soluions. (0 poins) Deermine wheher he following

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Chapter 6. Systems of First Order Linear Differential Equations

Chapter 6. Systems of First Order Linear Differential Equations Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh

More information

Two-dimensional Effects on the CSR Interaction Forces for an Energy-Chirped Bunch. Rui Li, J. Bisognano, R. Legg, and R. Bosch

Two-dimensional Effects on the CSR Interaction Forces for an Energy-Chirped Bunch. Rui Li, J. Bisognano, R. Legg, and R. Bosch Two-dimensional Effecs on he CS Ineacion Foces fo an Enegy-Chiped Bunch ui Li, J. Bisognano,. Legg, and. Bosch Ouline 1. Inoducion 2. Pevious 1D and 2D esuls fo Effecive CS Foce 3. Bunch Disibuion Vaiaion

More information

EE 330 Lecture 23. Small Signal Analysis Small Signal Modelling

EE 330 Lecture 23. Small Signal Analysis Small Signal Modelling EE 330 Lecure 23 Small Signal Analysis Small Signal Modelling Exam 2 Friday March 9 Exam 3 Friday April 13 Review Session for Exam 2: 6:00 p.m. on Thursday March 8 in Room Sweeney 1116 Review from Las

More information

Lecture 13 RC/RL Circuits, Time Dependent Op Amp Circuits

Lecture 13 RC/RL Circuits, Time Dependent Op Amp Circuits Lecure 13 RC/RL Circuis, Time Dependen Op Amp Circuis RL Circuis The seps involved in solving simple circuis conaining dc sources, resisances, and one energy-sorage elemen (inducance or capaciance) are:

More information

The sudden release of a large amount of energy E into a background fluid of density

The sudden release of a large amount of energy E into a background fluid of density 10 Poin explosion The sudden elease of a lage amoun of enegy E ino a backgound fluid of densiy ceaes a song explosion, chaaceized by a song shock wave (a blas wave ) emanaing fom he poin whee he enegy

More information

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC Ch.. Group Work Unis Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Uni 2 Jusify wheher he following saemens are rue or false: a) Two sreamlines, corresponding o a same insan of ime, can never inersec

More information

Problemas das Aulas Práticas

Problemas das Aulas Práticas Mesrado Inegrado em Engenharia Elecroécnica e de Compuadores Conrolo em Espaço de Esados Problemas das Aulas Práicas J. Miranda Lemos Fevereiro de 3 Translaed o English by José Gaspar, 6 J. M. Lemos, IST

More information

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

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006 2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 3 Signals & Sysems Prof. Mark Fowler Noe Se #2 Wha are Coninuous-Time Signals??? Reading Assignmen: Secion. of Kamen and Heck /22 Course Flow Diagram The arrows here show concepual flow beween ideas.

More information

5. Response of Linear Time-Invariant Systems to Random Inputs

5. Response of Linear Time-Invariant Systems to Random Inputs Sysem: 5. Response of inear ime-invarian Sysems o Random Inpus 5.. Discree-ime linear ime-invarian (IV) sysems 5... Discree-ime IV sysem IV sysem xn ( ) yn ( ) [ xn ( )] Inpu Signal Sysem S Oupu Signal

More information

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security 1 Geneal Non-Abiage Model I. Paial Diffeenial Equaion fo Picing A. aded Undelying Secuiy 1. Dynamics of he Asse Given by: a. ds = µ (S, )d + σ (S, )dz b. he asse can be eihe a sock, o a cuency, an index,

More information

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I Open loop vs Closed Loop Advanced I Moor Command Movemen Overview Open Loop vs Closed Loop Some examples Useful Open Loop lers Dynamical sysems CPG (biologically inspired ), Force Fields Feedback conrol

More information

ME 391 Mechanical Engineering Analysis

ME 391 Mechanical Engineering Analysis Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of

More information

Chapter 7 Response of First-order RL and RC Circuits

Chapter 7 Response of First-order RL and RC Circuits Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial

More information

AN EVOLUTIONARY APPROACH FOR SOLVING DIFFERENTIAL EQUATIONS

AN EVOLUTIONARY APPROACH FOR SOLVING DIFFERENTIAL EQUATIONS AN EVOLUTIONARY APPROACH FOR SOLVING DIFFERENTIAL EQUATIONS M. KAMESWAR RAO AND K.P. RAVINDRAN Depamen of Mechanical Engineeing, Calicu Regional Engineeing College, Keala-67 6, INDIA. Absac:- We eploe

More information

Tracking Adversarial Targets

Tracking Adversarial Targets A. Proofs Proof of Lemma 3. Consider he Bellman equaion λ + V π,l x, a lx, a + V π,l Ax + Ba, πax + Ba. We prove he lemma by showing ha he given quadraic form is he unique soluion of he Bellman equaion.

More information

4.5 Constant Acceleration

4.5 Constant Acceleration 4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),

More information

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH Fundamenal Jounal of Mahemaical Phsics Vol 3 Issue 013 Pages 55-6 Published online a hp://wwwfdincom/ MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH Univesias

More information

Chapter 5 Digital PID control algorithm. Hesheng Wang Department of Automation,SJTU 2016,03

Chapter 5 Digital PID control algorithm. Hesheng Wang Department of Automation,SJTU 2016,03 Chaper 5 Digial PID conrol algorihm Hesheng Wang Deparmen of Auomaion,SJTU 216,3 Ouline Absrac Quasi-coninuous PID conrol algorihm Improvemen of sandard PID algorihm Choosing parameer of PID regulaor Brief

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

Pseudosteady-State Flow Relations for a Radial System from Department of Petroleum Engineering Course Notes (1997)

Pseudosteady-State Flow Relations for a Radial System from Department of Petroleum Engineering Course Notes (1997) Pseudoseady-Sae Flow Relaions fo a Radial Sysem fom Deamen of Peoleum Engineeing Couse Noes (1997) (Deivaion of he Pseudoseady-Sae Flow Relaions fo a Radial Sysem) (Deivaion of he Pseudoseady-Sae Flow

More information

Homework-8(1) P8.3-1, 3, 8, 10, 17, 21, 24, 28,29 P8.4-1, 2, 5

Homework-8(1) P8.3-1, 3, 8, 10, 17, 21, 24, 28,29 P8.4-1, 2, 5 Homework-8() P8.3-, 3, 8, 0, 7, 2, 24, 28,29 P8.4-, 2, 5 Secion 8.3: The Response of a Firs Order Circui o a Consan Inpu P 8.3- The circui shown in Figure P 8.3- is a seady sae before he swich closes a

More information

Solutionbank Edexcel AS and A Level Modular Mathematics

Solutionbank Edexcel AS and A Level Modular Mathematics Page of 4 Soluionbank Edexcel AS and A Level Modular Mahemaics Exercise A, Quesion Quesion: Skech he graphs of (a) y = e x + (b) y = 4e x (c) y = e x 3 (d) y = 4 e x (e) y = 6 + 0e x (f) y = 00e x + 0

More information

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation Lecue 8: Kineics of Phase Gowh in a Two-componen Sysem: geneal kineics analysis based on he dilue-soluion appoximaion Today s opics: In he las Lecues, we leaned hee diffeen ways o descibe he diffusion

More information

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff

Laplace transfom: t-translation rule , Haynes Miller and Jeremy Orloff Laplace ransfom: -ranslaion rule 8.03, Haynes Miller and Jeremy Orloff Inroducory example Consider he sysem ẋ + 3x = f(, where f is he inpu and x he response. We know is uni impulse response is 0 for

More information

CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang

CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS Professor Dae Ryook Yang Fall 200 Dep. of Chemical and Biological Engineering Korea Universiy CHE302 Process Dynamics and Conrol Korea Universiy

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

ECE 2100 Circuit Analysis

ECE 2100 Circuit Analysis ECE 1 Circui Analysis Lesson 37 Chaper 8: Second Order Circuis Discuss Exam Daniel M. Liynski, Ph.D. Exam CH 1-4: On Exam 1; Basis for work CH 5: Operaional Amplifiers CH 6: Capaciors and Inducor CH 7-8:

More information

3, so θ = arccos

3, so θ = arccos Mahemaics 210 Professor Alan H Sein Monday, Ocober 1, 2007 SOLUTIONS This problem se is worh 50 poins 1 Find he angle beween he vecors (2, 7, 3) and (5, 2, 4) Soluion: Le θ be he angle (2, 7, 3) (5, 2,

More information

Math Final Exam Solutions

Math Final Exam Solutions Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,

More information

Analytic Model and Bilateral Approximation for Clocked Comparator

Analytic Model and Bilateral Approximation for Clocked Comparator Analyic Model and Bilaeral Approximaion for Clocked Comparaor M. Greians, E. Hermanis, G.Supols Insiue of, Riga, Lavia, e-mail: gais.supols@edi.lv Research is suppored by: 1) ESF projec Nr.1DP/1.1.1.2.0/09/APIA/VIAA/020,

More information

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

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

More information

CHE302 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES. Professor Dae Ryook Yang

CHE302 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES. Professor Dae Ryook Yang CHE30 LECTURE VI DYNAMIC BEHAVIORS OF REPRESENTATIVE PROCESSES Professor Dae Ryook Yang Fall 00 Dep. of Chemical and Biological Engineering orea Universiy CHE30 Process Dynamics and Conrol orea Universiy

More information

Structural Dynamics and Earthquake Engineering

Structural Dynamics and Earthquake Engineering Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/

More information

( ) ( ) ( ) () () Signals And Systems Exam#1. 1. Given x(t) and y(t) below: x(t) y(t) (A) Give the expression of x(t) in terms of step functions.

( ) ( ) ( ) () () Signals And Systems Exam#1. 1. Given x(t) and y(t) below: x(t) y(t) (A) Give the expression of x(t) in terms of step functions. Signals And Sysems Exam#. Given x() and y() below: x() y() 4 4 (A) Give he expression of x() in erms of sep funcions. (%) x () = q() q( ) + q( 4) (B) Plo x(.5). (%) x() g() = x( ) h() = g(. 5) = x(. 5)

More information

Math 4600: Homework 11 Solutions

Math 4600: Homework 11 Solutions Mah 46: Homework Soluions Gregory Handy [.] One of he well-known phenomenological (capuring he phenomena, bu no necessarily he mechanisms) models of cancer is represened by Gomperz equaion dn d = bn ln(n/k)

More information