into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

Size: px
Start display at page:

Download "into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get"

Transcription

1 Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }} Recall that the table of Laplace/-tranform i contructed by (i) electing h(t) to get h(k), and then taking the -tranform of h(k) to arrive at H () In the cae where, (ii) ampling i a ytem tranfer function, then h(t) i the ytem impule repone In thi cae, h(k) i the ampled impule repone In word, we are ampling in uch a manner that the ampled impule repone matche the original impule repone at the ample time here i no variability between the two at the ampling time For thi reaon, we have the following Fact : he Laplace/-tranform table i contructed via impule-invariant ampling Per the Matlab decription below, thi method i alo called the impule-invariant mapping method Chooing a Converion Method he c2d command dicretie continuou-time model Converely, d2c convert dicrete-time model to continuou time Both command upport everal dicretiation and interpolation method, a hown in the following table Dicretiation Method Ue when: Zero-Order Hold You want an exact dicretiation in the time domain for taircae input Firt-Order Hold You want an exact dicretiation in the time domain for piecewie linear input Impule-Invariant Mapping (c2d only) You want an exact dicretiation in the time domain for impule train input utin Approximation Zero-Pole Matching Equivalent You want good matching in the frequency domain between the continuou- and dicrete-time model Your model ha important dynamic at ome particular frequency You have a SISO model, and you want good matching in the frequency domain between the continuou- and dicrete-time model We will now addre the Zero-Order Hold (ZOH) method in the above table he following i Matlab decription Zero-Order Hold he Zero-Order Hold (ZOH) method provide an exact match between the continuou- and dicrete-time ytem in the time domain for taircae input he following block diagram illutrate the ero-order-hold dicretiation H d() of a continuou-time linear model he ZOH block generate the continuoutime input ignal u(t) by holding each ample value u(k) contant over one ample period he ignal u(t) i the input to the continuou ytem he output y[k] reult from ampling y(t) every econd Converely, given a dicrete ytem H d(), d2c produce a continuou ytem he ZOH dicretiation of coincide with H d() Clearly, the ZOH method i different from the impule-invariant method We will now elaborate on the above decription via the following example Example Conider the firt order ytem Y ( / U ( For a unit tep, Y ( Wherea the unit ( ) impule i not a taircae input, the unit tep i a trivial example of one (ie it ha only one tep) From a table of Laplace t tranform, we have y( t) e he above decription tate that, uing the ZOH method, we will have y k ( k) e for the ampled input u ( k) Uing the Matlab command

2 c2d with the oh flag (which, by the way, i the default flag), we have: H=tf(,[ ]); =005; H=c2d(H,,'oh') H =004877/(-0952) he Figure at right verifie the fact that the ampled ytem tep repone Matche that of the continuou time ytem For thi reaon, one could Decribe the ZOH method a tep-invariant ampling We will now develop the reult: H ( ) o thi end, recall that 0952 Figure Step repone for ZO, giving ( e ) ) ( )( e ZOH( From a table of Laplace/-tranform, we have ) ) Hence, the ampled verion of ZO i ) e H( ) he Bode plot for and are e 0952 hown at right A expected, we ee amplitude and phae aliaing near the Nyquit frequency / 6283r N / H () and H () he atute tudent might wonder why we did not include the factor of in the Laplace/-tranform pair above Recall that thi factor wa included in the impule-invariant ampling method, where the impule i Figure Bode plot for H ( and H () approximated by a pule with width and height / (hence, having area equal to ) In thi example the input i not an impule It i a more well-behaved function of time (ie a unit tep, which i a well-defined function) Hence, Matlab doe not include it in the c2d command uing the oh flag Uing the imp flag give: (t) H=c2d(H,,'imp') H =( e-8)/( ) = 005/(-0952) From the table of Laplace/-tranform we have: of i included ) e Since 0 05, it i clear that for thi flag, the factor At thi point, it i reaonable to ak: When hould I ue the imp flag, and when hould I ue the oh flag? he anwer to thi quetion i given in Matlab above block diagram When you have a dicrete-time input u (k) from a digital controller that you want to power a plant H ( with, you mut ue a D/A circuit, uch a the ZOH, to convert number into voltage Even though the output of H (, y (t), i a continuou-time output, we cannot eaily analye it in continuou time hi i becaue the digital input i only defined at dicrete time So, the imple way to conduct an analyi i to ample the output of at the ample time uing the oh flag On the other hand, if you want a dicrete approximation for H (, then you hould ue the imp flag Before leaving thi example, we will cover one more point Specifically, conider the three tranfer function: (i) (ii) H '( ZO (iii) ) H ( ) 0952

3 Recall that the ZOH tranfer function ZO ( ) / wa obtained by computing the continuou time repone to a continuou time unit impule, (t) Clearly, (i) and (ii) are different However, they can both be viewed a continuou time tranfer function We will now proceed to compare their FRF o thi end, we will firt cale (ii) to have unity tatic gain For mall value of, we have e Hence, we will cale (ii) by dividing by hi cale factor change the output of the ZOH form a pule of width and height 0 to a pule of height / And o, the input to i now an approximate unit impule Unlike (t), which ha equal energy at all frequencie, thi pule ha Laplace tranform ZO / ( ) / It FRF i hown at right It energy i contant at low frequencie However, it decay rapidly at frequencie above the Nyquit frequency hi obervation, in itelf, i an important one in relation to tet engineering A mall value of correpond to a hammer Figure 3 FRF for a pule of width =05 with a hard tip he ofter the tip, the le excitation energy at higher frequencie will be imparted to the tet object Hence, the high frequency dynamic of the object, not being excited, will remain hidden We ee thi effect in relation to he FRF for (ii) i hown at right he energy above the Nyquit frequency i notably le than it would be if were excited by In other word, uing the ZOH i akin to (t) applying a low pa filter to H ( he bandwidth of thi filter extend to the Nyquit frequency Figure 4 FRF for the tranfer function (ii) Viewing D/A Converion a Numerical Integration Recall that ZO / ( ) / wa arrived at by contructing it impule repone a the uperpoition of a unit tep followed by a negative unit tep delayed by, a hown at right Since the impule repone i h( t) u ( t) u ( t ), it tranfer function i e H( o view thi a an integration problem, write ( ) H( H 2( t From the block diagram at right, we have v( t) u( ) d and y( t) v( t) v( t ) u( ) d Hence, for a general input, u (t), the output, y (t), i imply thi integral Written in a recurive fahion, we have: 0 k y( k) y[( k ) u( ) d () ( k) t t u H ( v H 2 ( y

4 he ZOH A/D circuit approximate the integral in () a area of a rectangle having width and height thi cae, () become he -tranform of (2a) give : utin method doe not hold the value the average of u[( k ) and u(k) u[( k ) k u ( k) ( ) d u[( k ) In word, it approximate it a the ; that i, it hold the value u[( k ) for one ample period In y( k) y[( k ) u[( k ) (2a) u[( k ) p602 of the book) Uing thi integration method, (2a) change to: Y ( ) (2b) U ( ) Intead, it hold the value { u[( k ) u( k)}/ 2 In word, it hold hi method of integration i known a trapeoidal integration (ee Figure 87 on he -tranform of (3a) give: u[( k ) u[( k) y( k) y[( k ) (3a) 2 Y( ) U( ) 2 (3b) Conider now the tranfer function a From the table of Laplace/-tranform we have a QUESION: What value of hould we ubtitute into H ( to obtain H ()? H a e ( ) a ANSWER: Setting H( H( ) and olving for give: a ( a) e Recall that the table of Laplace/ tranform wa created uing the Riemann um method of numerical integration QUESION: What value of hould we ubtitute into H ( to obtain H () aociated with trapeoidal integration? ANSWER: 2 hi anwer i not eay to prove, even for the imple a However, it can be hown to hold for any H ( a Example 2 Suppoe that we have a command ytem that include a forward loop continuou time lead compenator with tranfer function / 2 G c ( 0 U( / E(, and that we want to replace it by a digital one /0 (including the D/A circuit) Compare the Bode plot of the digital compenator uing a oh and a tutin A/D for 0 0 Solution: he figure at right how that thee method are comparable Figure 2 Digital controller Bode plot

5 Summary he goal of thee note wa to introduce you to ome of the Matlab c2d command flag he imp flag can be ued a an e-table of Laplace/-tranform pair he oh and tutin flag are ued to incorporate an A/D converter Even though the output of an A/D converter i analog, the dicrete tranfer function allow one to evaluate thi analog output at the ampling time

Digital Control System

Digital Control System Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)

More information

Design By Emulation (Indirect Method)

Design By Emulation (Indirect Method) Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal

More information

Digital Control System

Digital Control System Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital

More information

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -

More information

( 1) EE 313 Linear Signals & Systems (Fall 2018) Solution Set for Homework #10 on Laplace Transforms

( 1) EE 313 Linear Signals & Systems (Fall 2018) Solution Set for Homework #10 on Laplace Transforms EE 33 Linear Signal & Sytem (Fall 08) Solution Set for Homework #0 on Laplace Tranform By: Mr. Houhang Salimian & Prof. Brian L. Evan Problem. a) xt () = ut () ut ( ) From lecture Lut { ()} = and { } t

More information

Design of Digital Filters

Design of Digital Filters Deign of Digital Filter Paley-Wiener Theorem [ ] ( ) If h n i a caual energy ignal, then ln H e dω< B where B i a finite upper bound. One implication of the Paley-Wiener theorem i that a tranfer function

More information

Sampling and the Discrete Fourier Transform

Sampling and the Discrete Fourier Transform Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at

More information

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou

More information

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder R. W. Erickon Department of Electrical, Computer, and Energy Engineering Univerity of Colorado, Boulder Cloed-loop buck converter example: Section 9.5.4 In ECEN 5797, we ued the CCM mall ignal model to

More information

Properties of Z-transform Transform 1 Linearity a

Properties of Z-transform Transform 1 Linearity a Midterm 3 (Fall 6 of EEG:. Thi midterm conit of eight ingle-ided page. The firt three page contain variou table followed by FOUR eam quetion and one etra workheet. You can tear out any page but make ure

More information

Chapter 2: Problem Solutions

Chapter 2: Problem Solutions Chapter 2: Solution Dicrete Time Proceing of Continuou Time Signal Sampling à 2.. : Conider a inuoidal ignal and let u ample it at a frequency F 2kHz. xt 3co000t 0. a) Determine and expreion for the ampled

More information

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

Module 4: Time Response of discrete time systems Lecture Note 1 Digital Control Module 4 Lecture Module 4: ime Repone of dicrete time ytem Lecture Note ime Repone of dicrete time ytem Abolute tability i a baic requirement of all control ytem. Apart from that, good

More information

EE Control Systems LECTURE 6

EE Control Systems LECTURE 6 Copyright FL Lewi 999 All right reerved EE - Control Sytem LECTURE 6 Updated: Sunday, February, 999 BLOCK DIAGRAM AND MASON'S FORMULA A linear time-invariant (LTI) ytem can be repreented in many way, including:

More information

March 18, 2014 Academic Year 2013/14

March 18, 2014 Academic Year 2013/14 POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of

More information

Lecture 10 Filtering: Applied Concepts

Lecture 10 Filtering: Applied Concepts Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering

More information

Lecture #9 Continuous time filter

Lecture #9 Continuous time filter Lecture #9 Continuou time filter Oliver Faut December 5, 2006 Content Review. Motivation......................................... 2 2 Filter pecification 2 2. Low pa..........................................

More information

5.5 Application of Frequency Response: Signal Filters

5.5 Application of Frequency Response: Signal Filters 44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order

More information

EE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis

EE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis EE/ME/AE34: Dynamical Sytem Chapter 8: Tranfer Function Analyi The Sytem Tranfer Function Conider the ytem decribed by the nth-order I/O eqn.: ( n) ( n 1) ( m) y + a y + + a y = b u + + bu n 1 0 m 0 Taking

More information

Wolfgang Hofle. CERN CAS Darmstadt, October W. Hofle feedback systems

Wolfgang Hofle. CERN CAS Darmstadt, October W. Hofle feedback systems Wolfgang Hofle Wolfgang.Hofle@cern.ch CERN CAS Darmtadt, October 9 Feedback i a mechanim that influence a ytem by looping back an output to the input a concept which i found in abundance in nature and

More information

Mathematical modeling of control systems. Laith Batarseh. Mathematical modeling of control systems

Mathematical modeling of control systems. Laith Batarseh. Mathematical modeling of control systems Chapter two Laith Batareh Mathematical modeling The dynamic of many ytem, whether they are mechanical, electrical, thermal, economic, biological, and o on, may be decribed in term of differential equation

More information

Chapter 4: Applications of Fourier Representations. Chih-Wei Liu

Chapter 4: Applications of Fourier Representations. Chih-Wei Liu Chapter 4: Application of Fourier Repreentation Chih-Wei Liu Outline Introduction Fourier ranform of Periodic Signal Convolution/Multiplication with Non-Periodic Signal Fourier ranform of Dicrete-ime Signal

More information

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002 Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in

More information

Roadmap for Discrete-Time Signal Processing

Roadmap for Discrete-Time Signal Processing EE 4G Note: Chapter 8 Continuou-time Signal co(πf Roadmap for Dicrete-ime Signal Proceing.5 -.5 -..4.6.8..4.6.8 Dicrete-time Signal (Section 8.).5 -.5 -..4.6.8..4.6.8 Sampling Period econd (or ampling

More information

Bogoliubov Transformation in Classical Mechanics

Bogoliubov Transformation in Classical Mechanics Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How

More information

DYNAMIC MODELS FOR CONTROLLER DESIGN

DYNAMIC MODELS FOR CONTROLLER DESIGN DYNAMIC MODELS FOR CONTROLLER DESIGN M.T. Tham (996,999) Dept. of Chemical and Proce Engineering Newcatle upon Tyne, NE 7RU, UK.. INTRODUCTION The problem of deigning a good control ytem i baically that

More information

ME 375 FINAL EXAM Wednesday, May 6, 2009

ME 375 FINAL EXAM Wednesday, May 6, 2009 ME 375 FINAL EXAM Wedneday, May 6, 9 Diviion Meckl :3 / Adam :3 (circle one) Name_ Intruction () Thi i a cloed book examination, but you are allowed three ingle-ided 8.5 crib heet. A calculator i NOT allowed.

More information

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3

More information

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning

More information

Question 1 Equivalent Circuits

Question 1 Equivalent Circuits MAE 40 inear ircuit Fall 2007 Final Intruction ) Thi exam i open book You may ue whatever written material you chooe, including your cla note and textbook You may ue a hand calculator with no communication

More information

EE Control Systems LECTURE 14

EE Control Systems LECTURE 14 Updated: Tueday, March 3, 999 EE 434 - Control Sytem LECTURE 4 Copyright FL Lewi 999 All right reerved ROOT LOCUS DESIGN TECHNIQUE Suppoe the cloed-loop tranfer function depend on a deign parameter k We

More information

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. R 4 := 100 kohm

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. R 4 := 100 kohm SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2003. Cacaded Op Amp [DC&L, problem 4.29] An ideal op amp ha an output impedance of zero,

More information

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL 98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i

More information

EE 477 Digital Signal Processing. 4 Sampling; Discrete-Time

EE 477 Digital Signal Processing. 4 Sampling; Discrete-Time EE 477 Digital Signal Proceing 4 Sampling; Dicrete-Time Sampling a Continuou Signal Obtain a equence of ignal ample uing a periodic intantaneou ampler: x [ n] = x( nt ) Often plot dicrete ignal a dot or

More information

Spring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material

Spring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material Spring 4 EE 445S Real-Time Digital Signal Proceing Laboratory Prof. Evan Homework # Solution on Review of Signal and Sytem Material Problem.. Continuou-Time Sinuoidal Generation. In practice, we cannot

More information

SAMPLING. Sampling is the acquisition of a continuous signal at discrete time intervals and is a fundamental concept in real-time signal processing.

SAMPLING. Sampling is the acquisition of a continuous signal at discrete time intervals and is a fundamental concept in real-time signal processing. SAMPLING Sampling i the acquiition of a continuou ignal at dicrete time interval and i a fundamental concept in real-time ignal proceing. he actual ampling operation can alo be defined by the figure belo

More information

ME2142/ME2142E Feedback Control Systems

ME2142/ME2142E Feedback Control Systems Root Locu Analyi Root Locu Analyi Conider the cloed-loop ytem R + E - G C B H The tranient repone, and tability, of the cloed-loop ytem i determined by the value of the root of the characteritic equation

More information

MAE140 Linear Circuits Fall 2012 Final, December 13th

MAE140 Linear Circuits Fall 2012 Final, December 13th MAE40 Linear Circuit Fall 202 Final, December 3th Intruction. Thi exam i open book. You may ue whatever written material you chooe, including your cla note and textbook. You may ue a hand calculator with

More information

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical

More information

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281 72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition

More information

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax:

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax: Control Sytem Engineering ( Chapter 7. Steady-State Error Prof. Kwang-Chun Ho kwangho@hanung.ac.kr Tel: 0-760-453 Fax:0-760-4435 Introduction In thi leon, you will learn the following : How to find the

More information

Midterm Test Nov 10, 2010 Student Number:

Midterm Test Nov 10, 2010 Student Number: Mathematic 265 Section: 03 Verion A Full Name: Midterm Tet Nov 0, 200 Student Number: Intruction: There are 6 page in thi tet (including thi cover page).. Caution: There may (or may not) be more than one

More information

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005.

SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuits II. Solutions to Assignment 3 February 2005. SIMON FRASER UNIVERSITY School of Engineering Science ENSC 320 Electric Circuit II Solution to Aignment 3 February 2005. Initial Condition Source 0 V battery witch flip at t 0 find i 3 (t) Component value:

More information

EE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505)

EE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505) EE 4443/5329 LAB 3: Control of Indutrial Sytem Simulation and Hardware Control (PID Deign) The Inverted Pendulum (ECP Sytem-Model: 505) Compiled by: Nitin Swamy Email: nwamy@lakehore.uta.edu Email: okuljaca@lakehore.uta.edu

More information

Lecture 8: Period Finding: Simon s Problem over Z N

Lecture 8: Period Finding: Simon s Problem over Z N Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing

More information

Root Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0

Root Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0 Objective Root Locu Diagram Upon completion of thi chapter you will be able to: Plot the Root Locu for a given Tranfer Function by varying gain of the ytem, Analye the tability of the ytem from the root

More information

Linear System Fundamentals

Linear System Fundamentals Linear Sytem Fundamental MEM 355 Performance Enhancement of Dynamical Sytem Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Content Sytem Repreentation Stability Concept

More information

ME 375 FINAL EXAM SOLUTIONS Friday December 17, 2004

ME 375 FINAL EXAM SOLUTIONS Friday December 17, 2004 ME 375 FINAL EXAM SOLUTIONS Friday December 7, 004 Diviion Adam 0:30 / Yao :30 (circle one) Name Intruction () Thi i a cloed book eamination, but you are allowed three 8.5 crib heet. () You have two hour

More information

ECE 3510 Root Locus Design Examples. PI To eliminate steady-state error (for constant inputs) & perfect rejection of constant disturbances

ECE 3510 Root Locus Design Examples. PI To eliminate steady-state error (for constant inputs) & perfect rejection of constant disturbances ECE 350 Root Locu Deign Example Recall the imple crude ervo from lab G( ) 0 6.64 53.78 σ = = 3 23.473 PI To eliminate teady-tate error (for contant input) & perfect reection of contant diturbance Note:

More information

Chapter 13. Root Locus Introduction

Chapter 13. Root Locus Introduction Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will

More information

ECE382/ME482 Spring 2004 Homework 4 Solution November 14,

ECE382/ME482 Spring 2004 Homework 4 Solution November 14, ECE382/ME482 Spring 2004 Homework 4 Solution November 14, 2005 1 Solution to HW4 AP4.3 Intead of a contant or tep reference input, we are given, in thi problem, a more complicated reference path, r(t)

More information

A Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking

A Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking A Simple Approach to Syntheizing Naïve Quantized Control for Reference Tracking SHIANG-HUA YU Department of Electrical Engineering National Sun Yat-Sen Univerity 70 Lien-Hai Road, Kaohiung 804 TAIAN Abtract:

More information

Lecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004

Lecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004 METR4200 Advanced Control Lecture 4 Chapter Nie Controller Deign via Frequency Repone G. Hovland 2004 Deign Goal Tranient repone via imple gain adjutment Cacade compenator to improve teady-tate error Cacade

More information

MM1: Basic Concept (I): System and its Variables

MM1: Basic Concept (I): System and its Variables MM1: Baic Concept (I): Sytem and it Variable A ytem i a collection of component which are coordinated together to perform a function Sytem interact with their environment. The interaction i defined in

More information

376 CHAPTER 6. THE FREQUENCY-RESPONSE DESIGN METHOD. D(s) = we get the compensated system with :

376 CHAPTER 6. THE FREQUENCY-RESPONSE DESIGN METHOD. D(s) = we get the compensated system with : 376 CHAPTER 6. THE FREQUENCY-RESPONSE DESIGN METHOD Therefore by applying the lead compenator with ome gain adjutment : D() =.12 4.5 +1 9 +1 we get the compenated ytem with : PM =65, ω c = 22 rad/ec, o

More information

ME 375 EXAM #1 Tuesday February 21, 2006

ME 375 EXAM #1 Tuesday February 21, 2006 ME 375 EXAM #1 Tueday February 1, 006 Diviion Adam 11:30 / Savran :30 (circle one) Name Intruction (1) Thi i a cloed book examination, but you are allowed one 8.5x11 crib heet. () You have one hour to

More information

Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:

Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web:     Ph: Serial :. PT_EE_A+C_Control Sytem_798 Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubanewar olkata Patna Web: E-mail: info@madeeay.in Ph: -4546 CLASS TEST 8-9 ELECTRICAL ENGINEERING Subject

More information

ECE-320 Linear Control Systems. Spring 2014, Exam 1. No calculators or computers allowed, you may leave your answers as fractions.

ECE-320 Linear Control Systems. Spring 2014, Exam 1. No calculators or computers allowed, you may leave your answers as fractions. ECE-0 Linear Control Sytem Spring 04, Exam No calculator or computer allowed, you may leave your anwer a fraction. All problem are worth point unle noted otherwie. Total /00 Problem - refer to the unit

More information

Linearteam tech paper. The analysis of fourth-order state variable filter and it s application to Linkwitz- Riley filters

Linearteam tech paper. The analysis of fourth-order state variable filter and it s application to Linkwitz- Riley filters Linearteam tech paper The analyi of fourth-order tate variable filter and it application to Linkwitz- iley filter Janne honen 5.. TBLE OF CONTENTS. NTOCTON.... FOTH-OE LNWTZ-LEY (L TNSFE FNCTON.... TNSFE

More information

55:041 Electronic Circuits

55:041 Electronic Circuits 55:04 Electronic ircuit Frequency epone hapter 7 A. Kruger Frequency epone- ee page 4-5 of the Prologue in the text Important eview co Thi lead to the concept of phaor we encountered in ircuit In Linear

More information

FRTN10 Exercise 3. Specifications and Disturbance Models

FRTN10 Exercise 3. Specifications and Disturbance Models FRTN0 Exercie 3. Specification and Diturbance Model 3. A feedback ytem i hown in Figure 3., in which a firt-order proce if controlled by an I controller. d v r u 2 z C() P() y n Figure 3. Sytem in Problem

More information

EE C128 / ME C134 Problem Set 1 Solution (Fall 2010) Wenjie Chen and Jansen Sheng, UC Berkeley

EE C128 / ME C134 Problem Set 1 Solution (Fall 2010) Wenjie Chen and Jansen Sheng, UC Berkeley EE C28 / ME C34 Problem Set Solution (Fall 200) Wenjie Chen and Janen Sheng, UC Berkeley. (0 pt) BIBO tability The ytem h(t) = co(t)u(t) i not BIBO table. What i the region of convergence for H()? A bounded

More information

Introduction to Laplace Transform Techniques in Circuit Analysis

Introduction to Laplace Transform Techniques in Circuit Analysis Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found

More information

Root Locus Contents. Root locus, sketching algorithm. Root locus, examples. Root locus, proofs. Root locus, control examples

Root Locus Contents. Root locus, sketching algorithm. Root locus, examples. Root locus, proofs. Root locus, control examples Root Locu Content Root locu, ketching algorithm Root locu, example Root locu, proof Root locu, control example Root locu, influence of zero and pole Root locu, lead lag controller deign 9 Spring ME45 -

More information

Department of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002

Department of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002 Department of Mechanical Engineering Maachuett Intitute of Technology 2.010 Modeling, Dynamic and Control III Spring 2002 SOLUTIONS: Problem Set # 10 Problem 1 Etimating tranfer function from Bode Plot.

More information

Lecture 8 - SISO Loop Design

Lecture 8 - SISO Loop Design Lecture 8 - SISO Loop Deign Deign approache, given pec Loophaping: in-band and out-of-band pec Fundamental deign limitation for the loop Gorinevky Control Engineering 8-1 Modern Control Theory Appy reult

More information

Main Topics: The Past, H(s): Poles, zeros, s-plane, and stability; Decomposition of the complete response.

Main Topics: The Past, H(s): Poles, zeros, s-plane, and stability; Decomposition of the complete response. EE202 HOMEWORK PROBLEMS SPRING 18 TO THE STUDENT: ALWAYS CHECK THE ERRATA on the web. Quote for your Parent' Partie: 1. Only with nodal analyi i the ret of the emeter a poibility. Ray DeCarlo 2. (The need

More information

What lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine?

What lies between Δx E, which represents the steam valve, and ΔP M, which is the mechanical power into the synchronous machine? A 2.0 Introduction In the lat et of note, we developed a model of the peed governing mechanim, which i given below: xˆ K ( Pˆ ˆ) E () In thee note, we want to extend thi model o that it relate the actual

More information

Lecture 8. PID control. Industrial process control ( today) PID control. Insights about PID actions

Lecture 8. PID control. Industrial process control ( today) PID control. Insights about PID actions Lecture 8. PID control. The role of P, I, and D action 2. PID tuning Indutrial proce control (92... today) Feedback control i ued to improve the proce performance: tatic performance: for contant reference,

More information

Function and Impulse Response

Function and Impulse Response Tranfer Function and Impule Repone Solution of Selected Unolved Example. Tranfer Function Q.8 Solution : The -domain network i hown in the Fig... Applying VL to the two loop, R R R I () I () L I () L V()

More information

Massachusetts Institute of Technology Dynamics and Control II

Massachusetts Institute of Technology Dynamics and Control II I E Maachuett Intitute of Technology Department of Mechanical Engineering 2.004 Dynamic and Control II Laboratory Seion 5: Elimination of Steady-State Error Uing Integral Control Action 1 Laboratory Objective:

More information

Figure 1 Siemens PSSE Web Site

Figure 1 Siemens PSSE Web Site Stability Analyi of Dynamic Sytem. In the lat few lecture we have een how mall ignal Lalace domain model may be contructed of the dynamic erformance of ower ytem. The tability of uch ytem i a matter of

More information

Control Systems Analysis and Design by the Root-Locus Method

Control Systems Analysis and Design by the Root-Locus Method 6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If

More information

Solving Differential Equations by the Laplace Transform and by Numerical Methods

Solving Differential Equations by the Laplace Transform and by Numerical Methods 36CH_PHCalter_TechMath_95099 3//007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the

More information

NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE

NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE POLITONG SHANGHAI BASIC AUTOMATIC CONTROL June Academic Year / Exam grade NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE Ue only thee page (including the bac) for anwer. Do not ue additional

More information

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder R. W. Erickon Department of Electrical, Computer, and Energy Engineering Univerity of Colorado, Boulder ZOH: Sampled Data Sytem Example v T Sampler v* H Zero-order hold H v o e = 1 T 1 v *( ) = v( jkω

More information

MODERN CONTROL SYSTEMS

MODERN CONTROL SYSTEMS MODERN CONTROL SYSTEMS Lecture 1 Root Locu Emam Fathy Department of Electrical and Control Engineering email: emfmz@aat.edu http://www.aat.edu/cv.php?dip_unit=346&er=68525 1 Introduction What i root locu?

More information

Chapter 7. Root Locus Analysis

Chapter 7. Root Locus Analysis Chapter 7 Root Locu Analyi jw + KGH ( ) GH ( ) - K 0 z O 4 p 2 p 3 p Root Locu Analyi The root of the cloed-loop characteritic equation define the ytem characteritic repone. Their location in the complex

More information

Chapter #4 EEE8013. Linear Controller Design and State Space Analysis. Design of control system in state space using Matlab

Chapter #4 EEE8013. Linear Controller Design and State Space Analysis. Design of control system in state space using Matlab EEE83 hapter #4 EEE83 Linear ontroller Deign and State Space nalyi Deign of control ytem in tate pace uing Matlab. ontrollabilty and Obervability.... State Feedback ontrol... 5 3. Linear Quadratic Regulator

More information

MEM 355 Performance Enhancement of Dynamical Systems Root Locus Analysis

MEM 355 Performance Enhancement of Dynamical Systems Root Locus Analysis MEM 355 Performance Enhancement of Dynamical Sytem Root Locu Analyi Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Outline The root locu method wa introduced by Evan in

More information

e st t u(t 2) dt = lim t dt = T 2 2 e st = T e st lim + e st

e st t u(t 2) dt = lim t dt = T 2 2 e st = T e st lim + e st Math 46, Profeor David Levermore Anwer to Quetion for Dicuion Friday, 7 October 7 Firt Set of Quetion ( Ue the definition of the Laplace tranform to compute Lf]( for the function f(t = u(t t, where u i

More information

G(s) = 1 s by hand for! = 1, 2, 5, 10, 20, 50, and 100 rad/sec.

G(s) = 1 s by hand for! = 1, 2, 5, 10, 20, 50, and 100 rad/sec. 6003 where A = jg(j!)j ; = tan Im [G(j!)] Re [G(j!)] = \G(j!) 2. (a) Calculate the magnitude and phae of G() = + 0 by hand for! =, 2, 5, 0, 20, 50, and 00 rad/ec. (b) ketch the aymptote for G() according

More information

Chapter 4. The Laplace Transform Method

Chapter 4. The Laplace Transform Method Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination

More information

Social Studies 201 Notes for March 18, 2005

Social Studies 201 Notes for March 18, 2005 1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,

More information

LTV System Modelling

LTV System Modelling Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction

More information

Laplace Transformation

Laplace Transformation Univerity of Technology Electromechanical Department Energy Branch Advance Mathematic Laplace Tranformation nd Cla Lecture 6 Page of 7 Laplace Tranformation Definition Suppoe that f(t) i a piecewie continuou

More information

Lecture 2: The z-transform

Lecture 2: The z-transform 5-59- Control Sytem II FS 28 Lecture 2: The -Tranform From the Laplace Tranform to the tranform The Laplace tranform i an integral tranform which take a function of a real variable t to a function of a

More information

What is automatic control all about? Welcome to Automatic Control III!! "Automatic control is the art of getting things to behave as you want.

What is automatic control all about? Welcome to Automatic Control III!! Automatic control is the art of getting things to behave as you want. What i automatic control all about? Welcome to Automatic Control III!! Lecture Introduction and linear ytem theory Thoma Schön Diviion of Sytem and Control Department of Information Technology Uppala Univerity.

More information

Moment of Inertia of an Equilateral Triangle with Pivot at one Vertex

Moment of Inertia of an Equilateral Triangle with Pivot at one Vertex oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.

More information

POWER SYSTEM SMALL SIGNAL STABILITY ANALYSIS BASED ON TEST SIGNAL

POWER SYSTEM SMALL SIGNAL STABILITY ANALYSIS BASED ON TEST SIGNAL POWE YEM MALL INAL ABILIY ANALYI BAE ON E INAL Zheng Xu, Wei hao, Changchun Zhou Zheang Univerity, Hangzhou, 37 PChina Email: hvdc@ceezueducn Abtract - In thi paper, a method baed on ome tet ignal (et

More information

SKEE 3143 CONTROL SYSTEM DESIGN. CHAPTER 3 Compensator Design Using the Bode Plot

SKEE 3143 CONTROL SYSTEM DESIGN. CHAPTER 3 Compensator Design Using the Bode Plot SKEE 3143 CONTROL SYSTEM DESIGN CHAPTER 3 Compenator Deign Uing the Bode Plot 1 Chapter Outline 3.1 Introduc4on Re- viit to Frequency Repone, ploang frequency repone, bode plot tability analyi. 3.2 Gain

More information

The Root Locus Method

The Root Locus Method The Root Locu Method MEM 355 Performance Enhancement of Dynamical Sytem Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Outline The root locu method wa introduced by Evan

More information

MATHEMATICAL MODELS OF PHYSICAL SYSTEMS

MATHEMATICAL MODELS OF PHYSICAL SYSTEMS hapter MATHEMATIAL MODELS OF PHYSIAL SYSTEMS.. INTODUTION For the analyi and deign of control ytem, we need to formulate a mathematical decription of the ytem. The proce of obtaining the deired mathematical

More information

TMA4125 Matematikk 4N Spring 2016

TMA4125 Matematikk 4N Spring 2016 Norwegian Univerity of Science and Technology Department of Mathematical Science TMA45 Matematikk 4N Spring 6 Solution to problem et 6 In general, unle ele i noted, if f i a function, then F = L(f denote

More information

MAHALAKSHMI ENGINEERING COLLEGE-TRICHY

MAHALAKSHMI ENGINEERING COLLEGE-TRICHY DIGITAL SIGNAL PROCESSING DEPT./SEM.: CSE /VII DIGITAL FILTER DESIGN-IIR & FIR FILTER DESIGN PART-A. Lit the different type of tructure for realiation of IIR ytem? AUC APR 09 The different type of tructure

More information

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems Control and Dynamical Sytem CDS 0 Problem Set #5 Iued: 3 Nov 08 Due: 0 Nov 08 Note: In the upper left hand corner of the econd page of your homework et, pleae put the number of hour that you pent on thi

More information

Filter Dispersion. with respect to frequency otherwise signal dispersion (and thus signal distortion) will result. Right?

Filter Dispersion. with respect to frequency otherwise signal dispersion (and thus signal distortion) will result. Right? 3/1/005 Filter Diperion.doc 1/6 Filter Diperion Any ignal that carrie ignificant information mut ha ome non-zero bandwidth. In other word, the ignal energy (a well a the information it carrie) i pread

More information

The Laplace Transform , Haynes Miller and Jeremy Orloff

The Laplace Transform , Haynes Miller and Jeremy Orloff The Laplace Tranform 8.3, Hayne Miller and Jeremy Orloff Laplace tranform baic: introduction An operator take a function a input and output another function. A tranform doe the ame thing with the added

More information

Factor Analysis with Poisson Output

Factor Analysis with Poisson Output Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu

More information

1 Routh Array: 15 points

1 Routh Array: 15 points EE C28 / ME34 Problem Set 3 Solution Fall 2 Routh Array: 5 point Conider the ytem below, with D() k(+), w(t), G() +2, and H y() 2 ++2 2(+). Find the cloed loop tranfer function Y () R(), and range of k

More information