Session outline. Introduction to Feedback Control. The Idea of Feedback. Automatic control. Basic setting. The feedback principle
|
|
- Lynne Cummings
- 6 years ago
- Views:
Transcription
1 Session otline Intodction to Feedback Contol Kal-Eik Åzen, Anton Cevin Feedback and feedfowad PID Contol State-space models Tansfe fnction models Contol design sing pole placement State feedback and obseves 2 The Idea of Feedback Atomatic contol Feedback: Compae the actal eslt with the desied eslt. Take actions based on the diffeence. A seemingl simple idea that is temendosl powefl. Use of feedback has often been evoltiona. Feedback is also called closed loop contol. The opposite is feedfowad o open loop contol: make a plan and execte it. Feedback and feedfowad ae ke ideas in the discipline of contol. Use of models and feedback Activities: Modeling Analsis and simlation Contol design Implementation Inpt Distbance Otpt Pocess Contolle Refeence 3 4 Basic setting The feedback pinciple Model Contol Pocess Distbance A ve powefl idea, that often leads to evoltiona changes in the wa sstems ae designed. The pima paadigm in atomatic contol. Ref. signal e Feedback Σ Pocess Contolle Mst handle two tasks: Follow efeence signals, Compensate fo distbances How to do seveal things with the contol signal 5 Base coective action on an eo that has occed Closed loop 6
2 Popeties of feedback + Redces inflence of distbances + Redces effect of pocess vaiations + Does not eqie exact models Feeds senso noise into the sstem Ma lead to instabilit, e.g.: if the contolle has too high gain if the feedback loop contains too lage time delas Ref. signal The feedfowad pinciple Feedfowad contolle Measable Distbance Pocess Take coective action befoe an eo has occed Mease the distbance and compensate fo it Use the fact that the efeence signal is known and adjst the contol signal to the efeence signal Open loop 7 8 Popeties of feedfowad Example: Cise contol sing feedfowad + Redces effect of distbances that cannot be edced b feedback + Allows faste set-point changes, withot intodcing contol eos Reqies good models Reqies stable sstems Open loop Poblems? Desied speed Table Ca Meased speed 9 Example: Cise contol sing feedback Exempel: Cise contol sing feedback and feedfowad Desied speed Eo Σ Contolle Thottle Ca Meased speed Table Desied speed Σ Contolle Σ Ca Meased speed Closed loop Simple contolle: Eo > : incease thottle Eo < : decease thottle Both poactive and eactive 2
3 Example: Segwa The sevo poblem Focs on efeence vale changes: Distbance Model Contol Pocess Tpical design citeia: Rise time, T Oveshoot, M. M 2p e Wh is this pocess moe difficlt to contol? Unstable dnamics 3 Settling time, T s Stead-state eo, e... T T s 4 t The eglato poblem Session otline Focs on pocess distbances: Model Contol Tpical design citeia: Otpt vaiance Contol signal vaiance Pobabilit densit.5 Test limit Pocess Distbance Set point fo eglato with low vaiance Set point fo eglato with high vaiance Feedback and feedfowad PID Contol State-space models Tansfe fnction models Contol design sing pole placement State feedback and obseves Pocess otpt 5 6 Example: Oven On/off contol = 2 C Σ e Contolle Oven (t) = { min, e(t) < max, e(t) > actal tempeate desied tempeate heating element powe ( ) Oscillations 8
4 Popotional contol Popotional contol P-contolle: (t) =Ke(t) (K gain) Inceased gain K: Stationa eo Smalle stationa eo Lage oscillations 2 Popotional Integal contol ( ) PI-contolle: (t) =K e(t)+ t e(s)ds T i (T i integal time) Popotional Integal contol Smalle integal time T i (lage integal action): No stationa eo Lage oscillations Popotional Integal Deivative contol PID-contolle: ( ) (t) =K e(t)+ t T i d de(t) (T d deivative time) PID: Pesent, past, and fte e D I P t P-pat: needed fo fast esponse I-pat: needed to emove stationa eo D-pat: ma be needed to stabilize the pocess The deivative pat edces oscillations 23 24
5 Session otline Feedback and feedfowad PID Contol State-space models Tansfe fnction models Contol design sing pole placement State feedback and obseves Dnamical sstems Sstem Static sstem: (t) = f ( (t) ) (The otpt at time t onl depends on the inpt at time t.) Dnamical sstem: (t) = f ( x(), [, t] ) 25 (The otpt at time t depends on the initial state x() and the inpt fom time to t.) 26 Linea sstems We will mainl deal with linea, time-invaiant (LTI) sstems Fo linea sstems, the pinciple of speposition holds: Sstem 2 2 Sstem Nonlinea sstems Almost all eal sstems ae nonlinea limited inpt and otpt signals nonlinea pocess geomet fiction, tblence,... Can be lineaized aond an opeating point If thee is feedback, a simple linea model is often enogh Bt, alwas emembe the limitations of the model! α + β 2 α + β 2 Sstem State space fom Standad sstem foms A nmbe of fist-ode diffeential eqations Descibes what happens inside the sstem and how inpts and otpt ae connected to this Nmeicall speio The heitage of mechanics Tansfe fnction fom The tansfom of a highe-ode linea diffeential eqation Descibes the elationship between the inpt and the otpt The sstem is a black box Compact notation, convenient fo hand calclations The heitage of electical engineeing 29 Nonlinea state-space model: State Space Models x dx = f (x,...,x n, ). dx n = f n (x,...,x n, ) = (x,...,x n, ) Sstem Linea state-space model: dx = a x a n x n + b. dx n = a n x a nn x n + b n = c x c n x n + d 3
6 State Space Models Intodce vectos and matices fo compact notation: x x =. x n Example: Pendlm n sstem ode Nonlinea state-space model: dx = f (x, ) = (x, ) Linea state-space model: dx = Ax + B = Cx + D 3 Nonlinea state-space model (x = angle, x 2 = angla velocit): ẋ = x 2 ẋ 2 = ω 2 sin x + kcos x = x whee ω = l is the natal feqenc of the pendlm. 32 Lineaization A nonlinea sstem can be lineaized aond an eqilibim point, whee it holds dx = f (x, )= Make fist-ode Talo appoximations of f and aond (x, ): f (x, ) f (x, ) + f }{{} x = (x, ) (x, ) }{{} x = + (x, ) (x, ) (x x )+ f (x, )( ) (x x )+ (x, )( ) Lineaization Intodce new vaiables Δx = x x, Δ = och Δ = The sstem can now be witten as dδx = dx f = f (x, ) x Δ = (x, ) x In matix fom: (x, ) (x, ) dδx = AΔx + BΔ Δ = CΔx + DΔ Δx + f (x, )Δ Δx + (x, )Δ Example Pendlm Lineaize ẋ = x 2 = f (x, x 2, ) ẋ 2 = ω 2 sin x + kcos x = f 2 (x, x 2, ) = x = (x, x 2, ) aond the ppe (nstable) eqilibim x = π, x 2 =, =. The lineaized sstem is given b dδx = AΔx + BΔ Δ = CΔx + DΔ whee Δx = x x, Δ =, Δ = and A = f f x x 2 f 2 f 2 x x 2 = ω 2 B = C = f f 2 x (x, ) (x, ) = = k cos x x 2 = (x, ) ω 2 cos x k sin x (x, ) (x, ) = k 35 D = = 36
7 Solving the sstem eqation The soltion to the sstem eqation {ẋ = Ax + B Stabilit concepts Stable = Cx + D is given b x(t) =e At x()+ (t) =Ce At x()+ t t e A(t τ ) B(τ )dτ Ce A(t τ ) B(τ )dτ + D(t) Unstable Asmptoticall stable Stabilit definitions Assme ẋ = Ax, x() =x The sstem is stable if x(t) is limited fo all x. The sstem is asmptoticall stable if x(t) fo all x. The sstem is nstable if x(t) is nlimited fo some x. Stabilit citeia {ẋ = Ax x() =x x(t) =x e At The behavio of the soltion depends on the eigenvales of A All eigenvales have negative eal pat: As. stab. Some eigenvale has positive eal pat: Unstable No eigenvales with positive eal pat and no mltiple eigenvales on the imagina axis: Stable 39 4 Session otline Feedback and feedfowad PID Contol State-space models Tansfe fnction models Contol design sing pole placement State feedback and obseves Tansfe fnction fom Std the sstem in the (complex) feqenc domain: U(s) Y(s) G(s) U(s) Laplace tansfom of (t) Y(s) Laplace tansfom of (t) G(s) tansfe fnction Y(s) =G(s)U(s) 4 (if the initial state is assmed to be zeo) 42
8 Some opeatos/signals and thei Laplace tansfoms Definition: L f = F(s) = ( ) df Deivative: L = sf(s) ( ) Integal: L f = s F(s) Diac implse: Lδ = Step fnction: Lθ = s e st f (t) Ramp fnction: L(tθ) = s 2 Exponential fnction: L(e at θ)= s a 43 Fom tansfe fnction to state space fom {ẋ = Ax + B x() = = Cx + D { sx(s) =AX(s)+BU(s) Y(s) =CX(s)+DU(s) Y(s) = [ C(sI A) B + D ] U(s) G(s) =C(sI A) B + D = p(s) q(s) q(s) = det(si A) is called chaacteistic polnomial 44 Often, Poles and zeos The oots of p(s) ae called zeos The oots of q(s) ae called poles Note that G(s) = p(s) q(s) Poles of G(s) Eigenvales of A Calclating sstem esponses. Find the tansfe fnction G(s) of the sstem 2. Find the Laplace tansfom U(s) of the inpt (t) 3. Y(s) =G(s)U(s) 4. Use invese Laplace tansfom to find (t) Calclating sstem esponses Example: Compte the step esponse of G(s) = s+ Inpt: U(s) =L {θ(t)} = s Otpt: Y(s) =G(s)U(s) = s(s+) Otpt in the time domain: { } (t) =L = e t s(s + ) G(s) = s + a Step esponse of fist-ode sstems G(s) = s + a = T + st Time constant: T = a Static gain: G() =/a step esponse (t) = a ( e at ) 47 48
9 Step esponse of second-ode sstems Block diagams Real poles: G(s) = (s + a)(s + b) step esponse: (t) = ae bt be at ab ab(b a) U G G 2 Y Y = G 2 G U Complex poles: ω 2 G(s) = s 2 + 2ζωs + ω 2 U G Σ Y Y = ( G + G 2 ) U step esponse: (t) = ζ 2 sin(ω ζ 2 t + φ) G 2 φ = accosζ ω = ndamped feqenc (ω > ) ζ = elative damping ( < ζ < ) 49 5 U Σ G G 2 Y. Feqenc esponse Y = G 2 G ( U Y) Y( + G 2 G )=G 2 G U Y = G 2G + G 2 G U Given a stable sstem G(s), the inpt (t) = sinω t will, afte a tansient, give the otpt ( ) (t) = G(iω ) sin ω t + ag G(iω ) The stead-state otpt is also sinsoidal 5 52 Bode diagam Daw G(iω ) as a fnction of ω (in log-log scale) Amplitde/magnitde/gain diagam ag G(iω ) as a fnction of ω (in log-lin scale) Phase/angle diagam d(t) Example: low-pass filte + (t) =(t) G(s) = s + G(iω )= iω + G(iω ) = ω 2 + ag G(iω )= actan ω 53 54
10 Example: low-pass filte Bode Diagam Nqist Diagam Daw G(iω ) in a pola diagam when ω goes fom to Im G(iω) Magnitde (abs) Ultimate point 2 ϕ a Re G(iω) Phase (deg) 45 ω Feqenc (ad/sec) Im Re G(iω ) Example of Nqist Diagam G(s) = s + G(iω )= iω + = iω ω 2 + Session otline Feedback and feedfowad PID Contol State-space models Tansfe fnction models Contol design sing pole placement State feedback and obseves Small ω : G(iω ) Lage ω : G(iω ) ω i 2 ω Closed-loop contol Analsis of the standad feedback loop efeence Contolle contol Pocess distbances measement Σ e C(s) P(s) Pima goals of the contolle: Follow the efeence Reject distbances C(s): contolle P(s): pocess Closed-loop tansfe fnction (fom to ): Y = PC + PC R 59 Contol design: Choose C to get the desied behavio! 6
11 Example cise contol Assme that the elationship between the thottle and the speed is given b d = P(s) = 5 s +.2 Fist t to eglate the speed with a P-contolle: whee e(t) =(t) (t) (t) =Ke(t) The closed-loop tansfe fnction is given b 5 PC + PC = K s K = 5K s K s+.2 The gain K affects the pole of the closed-loop sstem the static gain of the closed-loop sstem 6 62 Simlation of the contol sstem with diffeent vales of K: Speed Thottle K= K= K=.3 K=.3 K=. Now t a PI-contolle: t ) (t) =K (e(t)+ Ti e(τ )dτ ( U(s) =K + ) st }{{ i } C(s) E(s) K= Time Stationa eo The closed-loop tansfe fnction is given b 5 PC + PC = K( ) + s+.2 st i 5K ( ) s + T + 5 K( ) = i + s+.2 st i s 2 +(5K +.2)s + 5K T i The poles of the closed-loop sstem depend on K and T i The static gain of the closed-loop sstem is alwas Simlation of the contol sstem with = 2, K =.3 and diffeent vales of T i : Speed Thottle Ti=.5 Ti=2 Ti=.5 Ti= Ti=2 Ti= Time 65 No stationa eo 66
12 Whee to place the poles? Stabilit nde Feedback Pole placement accoding to the chaacteistic polnomial q(s) = s 2 + 2ζω s + ω 2 : Im Σ G ϕ ω Re The closed loop sstem is asmptoticall stable if and onl if all the zeos of + G (s) Lage ω faste sstem esponse Smalle ϕ bette damping (elative damping ζ = cosϕ). (Common choice: ζ = cos 45 =.7) 67 lies in the left half plane. 68 The Nqist Citeion Example If G (s) is stable then the closed loop sstem [ + G (s)] is stable if and onl if the the Nqist cveg(iω ) does not encicle. G (s) = G P (s)g R (s), i.e. modif G P sch that the Nqist cve does not encicle 69 G (iω )= = = Σ K s(s+)(s+2) K iω ( + iω )(2 + iω ) Ki( iω )(2 iω ) ω ( + ω 2 )(4 + ω 2 ) = Ki(2 ω 2 3iω ) ω ( + ω 2 )(4 + ω 2 ) 3K ( + ω 2 )(4 + ω 2 ) + i K(ω 2 2) ω ( + ω 2 )(4 + ω 2 ) 7 Stabilit fo the closed loop sstem Amplitde and phase magins G (iω ) Amplitde magina m ag G(iω )=8, G(iω ) = A m Phase magin φ m G(iω c ) =, ag G(iω c )=φ m 8 Stable if K < 6. G (i 2)= 3K 3 6 = K 6 7 (Rles of Thmb: A m [2, 6], φ m [3,6 ]) 72
13 Magins in the Bode diagam Session otline Feedback and feedfowad PID Contol State-space models Tansfe fnction models Contol design sing pole placement State feedback and obseves Session otline Feedback and feedfowad PID Contol State-space models Tansfe fnction models Contol design sing pole placement State feedback and obseves State feedback Pocess: dx = Ax + B = Cx Assme that the fll state vecto x is measable. Contol law: = Lx + l l Σ P L x Closed-loop sstem: dx =(A BL)x + Bl = Cx The closed loop poles ae given b det(si A + BL)= Tning: L is chosen to give the desied poles l is chosen to give the static gain fom to Example - Inveted pendlm State vaiables x =, x 2 = ẏ dx = x + = x 77 78
14 Detemine a state feedback law (assme = ) = Lx = l l 2 x x 2 Simlation fom x() =[.75 ] T :.5 = x sch that the closed-loop chaacteistic polnomial becomes s 2 +.4s +..5 x 2 Closed-loop poles: s det(si A + BL)= + l s + l = s2 + l 2 s + l 2 A compaison with the desied polnomial gives l = 2 l 2 = Obseve It is most often not possible to mease the fll state vecto x. The state can then be estimated sing an obseve: Obseve: d ˆx = A ˆx + B + K( ŷ ) ŷ = C ˆx P Obseve ˆx Dnamics of the estimation eo x = x ˆx: d x = Ax + B A ˆx B KC(x ˆx) =(A KC) x Obseve poles: Tning: det(si A + KC)= 8 K is chosen to give the desied poles fast poles fast convegence ˆx x bt sensitive to noise 82 slow poles slow convegence bt obst The complete contolle (obseve + state feedback) is given b d ˆx = A ˆx + B + K( C ˆx) = L ˆx The tansfe fnction of the contolle is given b C(s) = L(sI A + BL+ KC) K State feedback fom estimated states: x 2 = x
f(y) signal norms system gain bounded input bounded output (BIBO) stability For what G(s) and f( ) is the closed-loop system stable?
Lecte 5 Inpt otpt stabilit Cose Otline o How to make a cicle ot of the point + i, and diffeent was to sta awa fom it... Lecte -3 Lecte 4-6 Modelling and basic phenomena (lineaization, phase plane, limit
More informationME 3600 Control Systems Frequency Domain Analysis
ME 3600 Contol Systems Fequency Domain Analysis The fequency esponse of a system is defined as the steady-state esponse of the system to a sinusoidal (hamonic) input. Fo linea systems, the esulting steady-state
More informationIntegral Control via Bias Estimation
1 Integal Contol via Bias stimation Consie the sstem ẋ = A + B +, R n, R p, R m = C +, R q whee is an nknown constant vecto. It is possible to view as a step istbance: (t) = 0 1(t). (If in fact (t) vaies
More informationHow to Obtain Desirable Transfer Functions in MIMO Systems Under Internal Stability Using Open and Closed Loop Control
How to Obtain Desiable ansfe Functions in MIMO Sstems Unde Intenal Stabilit Using Open and losed Loop ontol echnical Repot of the ISIS Goup at the Univesit of Note Dame ISIS-03-006 June, 03 Panos J. Antsaklis
More informationInternal Model Control Using Neural Networks.
Poceedings of the IEEE Intenational Smposim on Indstial Electonics, jne 17-2 1996, Wasaw. Intenal Contol Using Neal Netwoks. Isabelle Rivals and Léon Pesonnaz ESPCI, Laboatoie d'électoniqe, 1 e Vaqelin
More informationConventional Paper-I (a) Explain the concept of gradient. Determine the gradient of the given field: ( )
EE-Conventional Pape-I IES-013 www.gatefoum.com Conventional Pape-I-013 1. (a) Eplain the concept of gadient. Detemine the gadient of the given field: V ρzsin φ+ z cos φ+ρ What is polaization? In a dielectic
More informationChapter 12: Kinematics of a Particle 12.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS. u of the polar coordinate system are also shown in
ME 01 DYNAMICS Chapte 1: Kinematics of a Paticle Chapte 1 Kinematics of a Paticle A. Bazone 1.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS Pola Coodinates Pola coodinates ae paticlaly sitable fo solving
More informationMechatronic system design
Mechatonic system design Mechatonic system design wb2414 2013/2014 Couse pat 5 Motion contol Pof.i. R.H.Munnig Schmidt Mechatonic System Design 1 Lectue outline: What did you lean about PID motion contol
More informationMaterial. Lecture 8 Backlash and Quantization. Linear and Angular Backlash. Example: Parallel Kinematic Robot. Backlash.
Lectre 8 Backlash and Qantization Material Toda s Goal: To know models and compensation methods for backlash Lectre slides Be able to analze the effect of qantization errors Note: We are sing analsis methods
More informationFrom Errors to Uncertainties in Basic Course in Electrical Measurements. Vladimir Haasz, Milos Sedlacek
Fom Eos to ncetainties in asic Cose in Electical Measements Vladimi Haasz, Milos Sedlacek Czech Technical nivesity in Page, Faclty of Electical Engineeing, Technicka, CZ-667 Page, Czech epblic phone:40
More informationLec 6: State Feedback, Controllability, Integral Action
Lec 6: State Feedback, Controllability, Integral Action November 22, 2017 Lund University, Department of Automatic Control Controllability and Observability Example of Kalman decomposition 1 s 1 x 10 x
More informationChapter 9 Dynamic stability analysis III Lateral motion (Lectures 33 and 34)
Pof. E.G. Tulapukaa Stability and contol Chapte 9 Dynamic stability analysis Lateal motion (Lectues 33 and 34) Keywods : Lateal dynamic stability - state vaiable fom of equations, chaacteistic equation
More informationA reduced polytopic LPV synthesis for a sampling varying controller : experimentation with a T inverted pendulum
A edced polytopic LPV synthesis fo a sampling vaying contolle : expeimentation with a T inveted pendlm David Robet, Olivie Sename and Daniel Simon Abstact This pape deals with the adaptation of a ealtime
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationVector Control. Application to Induction Motor Control. DSP in Motion Control - Seminar
Vecto Contol Application to Induction Moto Contol Vecto Contol - Pinciple The Aim of Vecto Contol is to Oient the Flux Poducing Component of the Stato Cuent to some Suitable Flux Vecto unde all Opeating
More informationLinear System Theory (Fall 2011): Homework 1. Solutions
Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts
More informationQualifying Examination Electricity and Magnetism Solutions January 12, 2006
1 Qualifying Examination Electicity and Magnetism Solutions Januay 12, 2006 PROBLEM EA. a. Fist, we conside a unit length of cylinde to find the elationship between the total chage pe unit length λ and
More informationarxiv:gr-qc/ v1 17 Feb 2001
Scatteing poblem of scala wave in womhole geomety Sng-Won Kim Depatment of Science Edcation, Ewha Women s Univesity, Seol -75, Koea Jne 5, 7 In this pape, we stdy the scatteing poblem of the scala wave
More informationParticle Systems. University of Texas at Austin CS384G - Computer Graphics Fall 2010 Don Fussell
Paticle Systems Univesity of Texas at Austin CS384G - Compute Gaphics Fall 2010 Don Fussell Reading Requied: Witkin, Paticle System Dynamics, SIGGRAPH 97 couse notes on Physically Based Modeling. Witkin
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationSteady State and Transient Performance Analysis of Three Phase Induction Machine using MATLAB Simulations
Intenational Jounal of Recent Tends in Engineeing, Vol, No., May 9 Steady State and Tansient Pefomance Analysis of Thee Phase Induction Machine using MATAB Simulations Pof. Himanshu K. Patel Assistant
More informationElectric Charge and Field
lectic Chage and ield Chapte 6 (Giancoli) All sections ecept 6.0 (Gauss s law) Compaison between the lectic and the Gavitational foces Both have long ange, The electic chage of an object plas the same
More informationSkps Media
Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps Media Skps
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationChapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS
Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In Chaptes 2 and 4 we have studied kinematics, i.e., we descibed the motion of objects using paametes such as the position vecto, velocity, and acceleation without any insights
More informationIntro. Computer Control Systems: F9
Intro. Computer Control Systems: F9 State-feedback control and observers Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 21 dave.zachariah@it.uu.se F8: Quiz! 2 / 21 dave.zachariah@it.uu.se
More informationAbsolute Specifications: A typical absolute specification of a lowpass filter is shown in figure 1 where:
FIR FILTER DESIGN The design of an digital filte is caied out in thee steps: ) Specification: Befoe we can design a filte we must have some specifications. These ae detemined by the application. ) Appoximations
More informationFigure 1. We will begin by deriving a very general expression before returning to Equations 1 and 2 to determine the specifics.
Deivation of the Laplacian in Spheical Coodinates fom Fist Pinciples. Fist, let me state that the inspiation to do this came fom David Giffiths Intodction to Electodynamics textbook Chapte 1, Section 4.
More informationDiscrete LQ optimal control with integral action: A simple controller on incremental form for MIMO systems
Modeling, Identification and Contol, Vol., No., 1, pp. 5, ISSN 189 18 Discete LQ optimal contol with integal action: A simple contolle on incemental fom fo MIMO systems David Di Ruscio Telemak Univesity
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 41 Digital Signal Pocessing Pof. Mak Fowle Note Set #31 Linea Phase FIR Design Optimum Equiipple (Paks-McClellan) Reading: Sect. 1.2.4 1.2.6 of Poakis & Manolakis 1/2 Motivation The window method and
More informationPhysics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!
Physics 161 Fall 011 Exta Cedit Investigating Black Holes - olutions The Following is Woth 50 Points!!! This exta cedit assignment will investigate vaious popeties of black holes that we didn t have time
More informationFully developed free convective flow of a Williamson fluid in a vertical channel under the effect of a magnetic field
Available online at www.pelagiaeseachliba.com Advances in Applied Science Reseach,, 3 (4):49-499 ISSN: 976-86 CODEN (USA): AASRFC Fll developed fee convective flow of a Williamson flid in a vetical channel
More informationPhysics 221 Lecture 41 Nonlinear Absorption and Refraction
Physics 221 Lectue 41 Nonlinea Absoption and Refaction Refeences Meye-Aendt, pp. 97-98. Boyd, Nonlinea Optics, 1.4 Yaiv, Optical Waves in Cystals, p. 22 (Table of cystal symmeties) 1. Intoductoy Remaks.
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationChapter 5 Force and Motion
Chapte 5 Foce and Motion In chaptes 2 and 4 we have studied kinematics i.e. descibed the motion of objects using paametes such as the position vecto, velocity and acceleation without any insights as to
More informationAdvanced Problems of Lateral- Directional Dynamics!
Advanced Poblems of Lateal- Diectional Dynamics! Robet Stengel, Aicaft Flight Dynamics! MAE 331, 216 Leaning Objectives 4 th -ode dynamics! Steady-state esponse to contol! Tansfe functions! Fequency esponse!
More informationGravity and isostasy
Gavity and isostasy Reading: owle p60 74 Theoy of gavity Use two of Newton s laws: ) Univesal law of gavitation: Gmm = m m Univesal gavitational constant G=6.67 x 0 - Nm /kg ) Second law of motion: = ma
More informationFE FORMULATIONS FOR PLASTICITY
G These slides ae designed based on the book: Finite Elements in Plasticity Theoy and Pactice, D.R.J. Owen and E. Hinton, 970, Pineidge Pess Ltd., Swansea, UK. Couse Content: A INTRODUCTION AND OVERVIEW
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More information2E1252 Control Theory and Practice
2E1252 Control Theory and Practice Lectre 11: Actator satration and anti wind-p Learning aims After this lectre, yo shold nderstand how satration can case controller states to wind p know how to modify
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More informationPhysics 111 Lecture 5 (Walker: 3.3-6) Vectors & Vector Math Motion Vectors Sept. 11, 2009
Physics 111 Lectue 5 (Walke: 3.3-6) Vectos & Vecto Math Motion Vectos Sept. 11, 2009 Quiz Monday - Chap. 2 1 Resolving a vecto into x-component & y- component: Pola Coodinates Catesian Coodinates x y =
More informationCOUPLED MODELS OF ROLLING, SLIDING AND WHIRLING FRICTION
ENOC 008 Saint Petesbug Russia June 30-July 4 008 COUPLED MODELS OF ROLLING SLIDING AND WHIRLING FRICTION Alexey Kieenkov Ins ti tu te fo P ob le ms in Me ch an ic s Ru ss ia n Ac ad em y of Sc ie nc es
More informationSTABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR
HUNGARIAN JOURNAL OF INDUSTRY AND CHEMISTRY VESZPRÉM Vol. 42(2) pp. 109 113 (2014) STABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR ATTILA FODOR 1, ROLAND BÁLINT 1, ATTILA MAGYAR 1, AND
More informationA Power Method for Computing Square Roots of Complex Matrices
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 13, 39345 1997 ARTICLE NO. AY975517 A Powe Method fo Computing Squae Roots of Complex Matices Mohammed A. Hasan Depatment of Electical Engineeing, Coloado
More informationPHYS 110B - HW #7 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased
PHYS 0B - HW #7 Sping 2004, Solutions by David Pace Any efeenced euations ae fom Giffiths Poblem statements ae paaphased. Poblem 0.3 fom Giffiths A point chage,, moves in a loop of adius a. At time t 0
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 5
ECE 634 Sping 06 Pof. David R. Jacson ECE Dept. Notes 5 TM x Sface-Wave Soltion Poblem nde consideation: x h ε, µ z A TM x sface wave is popagating in the z diection (no y vaiation). TM x Sface-Wave Soltion
More informationFUNDAMENTAL PROPERTIES OF LINEAR SHIP STEERING DYNAMIC MODELS
Jounal of Maine Science and Technology, Vol. 7, No. 2, pp. 79-88 (1999) 79 FUNDAMENTAL PROPERTIES OF LINEAR SHIP STEERING DYNAMIC MODELS Ching-Yaw Tzeng* and Ju-Fen Chen** Keywods: Nomoto model, Contollability
More informationCHAPTER 3. Section 1. Modeling Population Growth
CHAPTER 3 Section 1. Modeling Population Gowth 1.1. The equation of the Malthusian model is Pt) = Ce t. Apply the initial condition P) = 1. Then 1 = Ce,oC = 1. Next apply the condition P1) = 3. Then 3
More informationc n ψ n (r)e ient/ h (2) where E n = 1 mc 2 α 2 Z 2 ψ(r) = c n ψ n (r) = c n = ψn(r)ψ(r)d 3 x e 2r/a0 1 πa e 3r/a0 r 2 dr c 1 2 = 2 9 /3 6 = 0.
Poblem {a} Fo t : Ψ(, t ψ(e iet/ h ( whee E mc α (α /7 ψ( e /a πa Hee we have used the gound state wavefunction fo Z. Fo t, Ψ(, t can be witten as a supeposition of Z hydogenic wavefunctions ψ n (: Ψ(,
More information1. State-Space Linear Systems 2. Block Diagrams 3. Exercises
LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space
More informationAN INTRODUCTION TO THE CONTROL THEORY
Open-Loop controller An Open-Loop (OL) controller is characterized by no direct connection between the output of the system and its input; therefore external disturbance, non-linear dynamics and parameter
More informationOn the reconstruction of the coronal magnetic field from coronal Hanle / Zeeman observations
On the econstction of the coonal magnetic field fom coonal Hanle / Zeeman obsevations M. Kama & B. Inheste Max-Planck Institte fo Sola System Reseach GERMANY Linda 005 Coonal Magnetic Field Magnetic field
More informationStep Response Analysis. Frequency Response, Relation Between Model Descriptions
Step Response Analysis. Frequency Response, Relation Between Model Descriptions Automatic Control, Basic Course, Lecture 3 November 9, 27 Lund University, Department of Automatic Control Content. Step
More informationInteraction of Feedforward and Feedback Streams in Visual Cortex in a Firing-Rate Model of Columnar Computations. ( r)
Supplementay mateial fo Inteaction of Feedfowad and Feedback Steams in Visual Cotex in a Fiing-Rate Model of Columna Computations Tobias Bosch and Heiko Neumann Institute fo Neual Infomation Pocessing
More information3. Electromagnetic Waves II
Lectue 3 - Electomagnetic Waves II 9 3. Electomagnetic Waves II Last time, we discussed the following. 1. The popagation of an EM wave though a macoscopic media: We discussed how the wave inteacts with
More informationHammerstein Model Identification Based On Instrumental Variable and Least Square Methods
Intenational Jounal of Emeging Tends & Technology in Compute Science (IJETTCS) Volume 2, Issue, Januay Febuay 23 ISSN 2278-6856 Hammestein Model Identification Based On Instumental Vaiable and Least Squae
More informationSpring 2001 Physics 2048 Test 3 solutions
Sping 001 Physics 048 Test 3 solutions Poblem 1. (Shot Answe: 15 points) a. 1 b. 3 c. 4* d. 9 e. 8 f. 9 *emembe that since KE = ½ mv, KE must be positive Poblem (Estimation Poblem: 15 points) Use momentum-impulse
More informationChapter 7 Singular Value Decomposition
EE448/58 Vesion. John Stensy Chapte 7 Singla Vale Decomposition This chapte coves the singla vale decomposition (SVD) of an m n matix. This algoithm has applications in many aeas inclding nmeical analysis
More informationResearch Design - - Topic 17 Multiple Regression & Multiple Correlation: Two Predictors 2009 R.C. Gardner, Ph.D.
Reseach Design - - Topic 7 Multiple Regession & Multiple Coelation: Two Pedictos 009 R.C. Gadne, Ph.D. Geneal Rationale and Basic Aithmetic fo two pedictos Patial and semipatial coelation Regession coefficients
More informationTheWaveandHelmholtzEquations
TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationNumerical Integration
MCEN 473/573 Chapte 0 Numeical Integation Fall, 2006 Textbook, 0.4 and 0.5 Isopaametic Fomula Numeical Integation [] e [ ] T k = h B [ D][ B] e B Jdsdt In pactice, the element stiffness is calculated numeically.
More information15 Solving the Laplace equation by Fourier method
5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the
More information{ } ( ) ( ) ( 1 ( ) π 2. Trapezoidal rule: a ( ) ( ) Simpson rule: 3/8 Simpson rule: Boole's rule a ( ) ( ) ( ) ( + ) 45 Weddle's Rule Hardy's rule:
L. Yaoslavsy. Selected Topics in Image Pocessing Pat. Fast tansfom methods fo image esampling Lect. 5. Pecise numeical integation and eentiation Conventional numeical integation methods: T T T = S S S
More informationLC transfer of energy between the driving source and the circuit will be a maximum.
The Q of oscillatos efeences: L.. Fotney Pinciples of Electonics: Analog and Digital, Hacout Bace Jovanovich 987, Chapte (AC Cicuits) H. J. Pain The Physics of Vibations and Waves, 5 th edition, Wiley
More informationAn Adaptive Neural-Network Model-Following Speed Control of PMSM Drives for Electric Vehicle Applications
Poceedings of the 9th WSEAS Intenational Confeence on Applied Mathematics, Istanbul, Tuey, May 27-29, 2006 (pp412-417) An Adaptive Neual-Netwo Model-Following Speed Contol of PMSM Dives fo Electic Vehicle
More informationData Flow Anomaly Analysis
Pof. D. Liggesmeye, 1 Contents Data flows an ata flow anomalies State machine fo ata flow anomaly analysis Example withot loops Example with loops Data Flow Anomaly Analysis Softwae Qality Assance Softwae
More information4/18/2005. Statistical Learning Theory
Statistical Leaning Theoy Statistical Leaning Theoy A model of supevised leaning consists of: a Envionment - Supplying a vecto x with a fixed but unknown pdf F x (x b Teache. It povides a desied esponse
More informationPAPER 39 STOCHASTIC NETWORKS
MATHEMATICAL TRIPOS Pat III Tuesday, 2 June, 2015 1:30 pm to 4:30 pm PAPER 39 STOCHASTIC NETWORKS Attempt no moe than FOUR questions. Thee ae FIVE questions in total. The questions cay equal weight. STATIONERY
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More informationDymore User s Manual Two- and three dimensional dynamic inflow models
Dymoe Use s Manual Two- and thee dimensional dynamic inflow models Contents 1 Two-dimensional finite-state genealized dynamic wake theoy 1 Thee-dimensional finite-state genealized dynamic wake theoy 1
More informationINPUT-OUTPUT APPROACH NUMERICAL EXAMPLES
INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES EXERCISE s consider the linear dnamical sstem of order 2 with transfer fnction with Determine the gain 2 (H) of the inpt-otpt operator H associated with this sstem.
More informationEXAM NMR (8N090) November , am
EXA NR (8N9) Novembe 5 9, 9. 1. am Remaks: 1. The exam consists of 8 questions, each with 3 pats.. Each question yields the same amount of points. 3. You ae allowed to use the fomula sheet which has been
More information4. Some Applications of first order linear differential
August 30, 2011 4-1 4. Some Applications of fist ode linea diffeential Equations The modeling poblem Thee ae seveal steps equied fo modeling scientific phenomena 1. Data collection (expeimentation) Given
More informationBasic Bridge Circuits
AN7 Datafoth Copoation Page of 6 DID YOU KNOW? Samuel Hunte Chistie (784-865) was bon in London the son of James Chistie, who founded Chistie's Fine At Auctionees. Samuel studied mathematics at Tinity
More informationPHYSICS 151 Notes for Online Lecture #20
PHYSICS 151 Notes fo Online Lectue #20 Toque: The whole eason that we want to woy about centes of mass is that we ae limited to looking at point masses unless we know how to deal with otations. Let s evisit
More informationMU+CU+KU=F MU+CU+KU=0
MEEN 67 Handout # MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING ^ Symmetic he motion of a n-dof linea system is descibed by the set of 2 nd ode diffeential equations MU+CU+KU=F t whee U (t) and F
More informationLab 10: Newton s Second Law in Rotation
Lab 10: Newton s Second Law in Rotation We can descibe the motion of objects that otate (i.e. spin on an axis, like a popelle o a doo) using the same definitions, adapted fo otational motion, that we have
More informationA New Approach to General Relativity
Apeion, Vol. 14, No. 3, July 7 7 A New Appoach to Geneal Relativity Ali Rıza Şahin Gaziosmanpaşa, Istanbul Tukey E-mail: aizasahin@gmail.com Hee we pesent a new point of view fo geneal elativity and/o
More informationn 1 Cov(X,Y)= ( X i- X )( Y i-y ). N-1 i=1 * If variable X and variable Y tend to increase together, then c(x,y) > 0
Covaiance and Peason Coelation Vatanian, SW 540 Both covaiance and coelation indicate the elationship between two (o moe) vaiables. Neithe the covaiance o coelation give the slope between the X and Y vaiable,
More informationVector d is a linear vector function of vector d when the following relationships hold:
Appendix 4 Dyadic Analysis DEFINITION ecto d is a linea vecto function of vecto d when the following elationships hold: d x = a xxd x + a xy d y + a xz d z d y = a yxd x + a yy d y + a yz d z d z = a zxd
More informationQuantum Mechanics II
Quantum Mechanics II Pof. Bois Altshule Apil 25, 2 Lectue 25 We have been dicussing the analytic popeties of the S-matix element. Remembe the adial wave function was u kl () = R kl () e ik iπl/2 S l (k)e
More informationRelated Rates - the Basics
Related Rates - the Basics In this section we exploe the way we can use deivatives to find the velocity at which things ae changing ove time. Up to now we have been finding the deivative to compae the
More informationSchool of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007
School of Electical and Compute Engineeing, Conell Univesity ECE 303: Electomagnetic Fields and Waves Fall 007 Homewok 8 Due on Oct. 19, 007 by 5:00 PM Reading Assignments: i) Review the lectue notes.
More informationdq 1 (5) q 1 where the previously mentioned limit has been taken.
1 Vecto Calculus And Continuum Consevation Equations In Cuvilinea Othogonal Coodinates Robet Maska: Novembe 25, 2008 In ode to ewite the consevation equations(continuit, momentum, eneg) to some cuvilinea
More informationR dt 2. INDIRECT FIELD ORIENTED CONTROL OF INDUCTION MOTOR
Indiect Field Oiented Contol of Induction Moto Using Fuzzy Logic Aun Kuma R PG Schola, PED, SELECT, VIT Univesity, Chennai Campus aunkumalmc@gmail.com Abstact: This pape pesents a fuzzy logic based speed
More informationTEMPERATURE RISE OF POWER TRANSFORMERS: COMMENTS AND PROPOSALS TO IEC :2011
Jonal of Enegy VOLUME 63 04 jonal homepage: http://jonalofenegy.com/ Zdenko Godec Konča - Electical Engineeing Institte Inc. zgodec@konca-institt.h Vjenceslav Kpešanin Konča - Electical Engineeing Institte
More informationClassical Mechanics Homework set 7, due Nov 8th: Solutions
Classical Mechanics Homewok set 7, due Nov 8th: Solutions 1. Do deivation 8.. It has been asked what effect does a total deivative as a function of q i, t have on the Hamiltonian. Thus, lets us begin with
More informationIs there a magnification paradox in gravitational lensing?
Is thee a magnification paadox in gavitational ing? Olaf Wucknitz wucknitz@asto.uni-bonn.de Astophysics semina/colloquium, Potsdam, 6 Novembe 7 Is thee a magnification paadox in gavitational ing? gavitational
More informationDOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS
DOING PHYIC WITH MTLB COMPUTTIONL OPTIC FOUNDTION OF CLR DIFFRCTION THEORY Ian Coope chool of Physics, Univesity of ydney ian.coope@sydney.edu.au DOWNLOD DIRECTORY FOR MTLB CRIPT View document: Numeical
More informationNon-Linear Dynamics Homework Solutions Week 2
Non-Linea Dynamics Homewok Solutions Week Chis Small Mach, 7 Please email me at smach9@evegeen.edu with any questions o concens eguading these solutions. Fo the ececises fom section., we sketch all qualitatively
More informationPhys101 Lectures 30, 31. Wave Motion
Phys0 Lectues 30, 3 Wave Motion Key points: Types of Waves: Tansvese and Longitudinal Mathematical Repesentation of a Taveling Wave The Pinciple of Supeposition Standing Waves; Resonance Ref: -7,8,9,0,,6,,3,6.
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect
More informationMultiple Experts with Binary Features
Multiple Expets with Binay Featues Ye Jin & Lingen Zhang Decembe 9, 2010 1 Intoduction Ou intuition fo the poect comes fom the pape Supevised Leaning fom Multiple Expets: Whom to tust when eveyone lies
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationEFFECTS OF FRINGING FIELDS ON SINGLE PARTICLE DYNAMICS. M. Bassetti and C. Biscari INFN-LNF, CP 13, Frascati (RM), Italy
Fascati Physics Seies Vol. X (998), pp. 47-54 4 th Advanced ICFA Beam Dynamics Wokshop, Fascati, Oct. -5, 997 EFFECTS OF FRININ FIELDS ON SINLE PARTICLE DYNAMICS M. Bassetti and C. Biscai INFN-LNF, CP
More informationB da = 0. Q E da = ε. E da = E dv
lectomagnetic Theo Pof Ruiz, UNC Asheville, doctophs on YouTube Chapte Notes The Maxwell quations in Diffeential Fom 1 The Maxwell quations in Diffeential Fom We will now tansfom the integal fom of the
More informationGauss s Law Simulation Activities
Gauss s Law Simulation Activities Name: Backgound: The electic field aound a point chage is found by: = kq/ 2 If thee ae multiple chages, the net field at any point is the vecto sum of the fields. Fo a
More information