Discrete time systems (sampling data systems) 1 Theoretical aspects. 2 Aims. 3 Establish the sampling period (Ts) Lab. no. 9

Size: px
Start display at page:

Download "Discrete time systems (sampling data systems) 1 Theoretical aspects. 2 Aims. 3 Establish the sampling period (Ts) Lab. no. 9"

Transcription

1 Amplitudine Discrete time systems (sampling data systems) Theoretical aspects Sampling theorem (Shannon s theorem) : Discretizing continuous time systems: Zero Order Hold element (zoh) Aims Obtain sampling data systems starting from the continuous time systems; Discretizing continuous transfer functions using Matlab tools; Simulate the time response of discrete time systems: step response, impulse response, ramp response and the time response to a sine wave of amplitude. 3 Establish the sampling period (Ts) When Shannon theorem is not respected, it will appear false frequencies inducing aliasing and folding effects. In the next graphic it can be seen the aliasing effect. x(t)=cos( ft) f=hz f=9hz f=hz timp (sec.) Figura Fenomenul de aliasing When dealing with sampling data control systems, the sampling period must be enough higher to assure not only data sampling and reconstruction but other different task like arithmetic calculations, memory allocation, etc. The sampling period will appear in the discrete (pulse) transfer function as a parameter.

2 Amplitudine To obtain correctly the behaviors of a discrete system (in comparison to the initial continuous time one, the sampling period must be connected to the time constants of the continuous transfer function. Having the continuous transfer function as a ratio of polynomials expressed by products of firs and second order terms: The sampling period value must be strictly smaller than the highest time constant in the considered transfer function: In practice, the sampling period is established times smaller than the highest time constant, or the smallest period of the input signal. When the sampling theorem is not respected, bad occurrence the real behavior is detected (figure )..8.6 Raspuns la treapta el. de ord. II el de ord. I timp (sec.) (seconds) Figure Incorrect sampling vs. the continuous time response 4 Discrete time signal representation in Matlab To graphically represent sampled signals the suggestion is to use the stem function in Matlab. To put in evidence the sampled and reconstructed signal, use the Matlab function stairs; Became familiar with the use of the two functions (access Matlab help) and elaborate a Matlab script in order to obtain the next graphical representations in Figure 3.

3 Amplitudine stairs vs. stem stairs stem timp (sec.) esantioane de esantioane de esantioane Figura 3 Funcţia stem şi stairs din Matlab %% figura 3 f=; % frecventa semnalului n=; % numar de esantioane Te=/f/n ; % perioada de esantionare t=:te:; % timpul de simulare esantionat cu Te x=cos(*pi*f*t); % functia cosinus de frecventa f subplot(3,,[,3,5]); stem(t,x);hold;stairs(t,x,'r');hold;% comparatie stem vs stairs legend('stem','stairs'); title('stem vs. stairs'); ylabel('cos(\pift)');xlabel('timp (sec.)') Te=/f/;% de esantioane Te3=/f/;% de esantioane t=:te:;t3=:te3:; x=cos(*pi*f*t);% semnal cu de esantioane pe o perioada x3=cos(*pi*f*t3);% semnal cu de esantioane pe o perioada % se reprezinta grafic semnalele subplot(3,,);stem(t,x);title(' esantioane'); subplot(3,,4);stem(t,x);title(' de esantioane'); subplot(3,,6);stem(t3,x3);title(' de esantioane'); 3

4 5 Discretizing the transfer function in Matlab There are two functions available to discretize transfer functions in Matlab: a) cdm (continuous to discrete method, the old version) The function syntax is as follows: [nz,dz]=cdm(nc,dc,te, metoda ) - nz and dz (output parameters) the nominator and denominator of the discrete transfer function; - nc and dc are the nominator and denominator of the continuous transfer function; - Te is the sampling period - metoda represents the discretizing method with the following values 'zoh' Convert to discrete time assuming a zero order hold on the inputs. foh' Convert to discrete time assuming a first order hold on the inputs. 'tustin' Convert to discrete time using the bilinear (Tustin) approximation to the derivative. 'prewarp' Convert to discrete time using the bilinear (Tustin) approximation with frequency prewarping. Specify the critical frequency with an additional argument, i.e cdm(a,b,c,d,ts,'prewarp',wc) 'matched' Convert the SISO system to discrete time using the matched pole-zero method. b) cd (continuous to discrete) This alternative of discretizing offers the possibility of working directly with transfer function objects. The function syntax is quite similar to the old one: sist_discret=cd(sist_continuu,te, metoda ) - in addition to the methods presented at the previous function, in this case it is available also: o imp (based on impulse response invariance) 6 Dynamic behaviors for discrete time systems After mapping the left half plane s into the z plane through the transformation, the theoretical results can be extended to the step response type (overdamped, critically damped, underdamped, undamped) depending on the pole location - the real and negative axis (from s plane) is mapped into real, positive and under unity axis in the z plane (having only positive and under unity real poles, the overdamped case will occur); - the content of the primary strip is mapped inside the unity circle in the z plane ( having at least one pair of complex conjugated poles inside the unity circle, the underdamped case will occur); - the imaginary axis from the s plane is mapped into the z plane exactly on the contour of the unity circle (having at least one pole situated on the unity circle, the system is only simple stable) for one pole in =, step response will be a ramp for a pair of complex conjugated poles on the unity circle, the undamped case will occur; - for the case when at least one pole of the pulse transfer function in situated outside the unity circle, the discrete system is unstable. 4

5 7 Solved problems For the continuous time system with the transfer function. Obtain the pulse transfer function using the zero order hold element. Hc=tf(5*[ 4],conv([ ],[ 3])); Te=max([/,/3])/;% establish sampling period Hd=cd(Hc,Te,'zoh'). Obtain the pulse transfer function as a ratio of polynomials expressed as products of first and second order; then the poles and zeros of the transfer function can be observed; for this case, the transfer function singularities (all real and negative) results after discretizing on the real and under unity axis. zpk(hd) 3. Modify the number representation format from default (4 decimals) into the long format (6 decimals) and reconsider the values of the pulse transfer function coefficients. format long [num,den]=tfdata(hd,'v') 4. Make the comparison between the step response from continuous and discrete time systems. Notice the stairs of the discrete time response that can be increased by increasing the sampling period. Step(Hc,Hd); legend('hc','hd') 5

6 3.5 Step Response Hc Hd Represent the step response of discrete time systems by discrete values and not by constant stairs. The syntax for the step function is the same as the one for the continuous case; when calling the step function with output parameters, there will no longer be a graphical representation but a memory allocation in order to make additional computations on the numerical response. [y,t]=step(hd); step(hc);hold;stem(t,y,'r');hold legend('hc','hd'); 3.5 Step Response Hc Hd Compare the step response for the cases when the transfer function coefficients are rounded to the first 4 decimals, respectively to the first decimals. 6

7 Hd4=tf([ 3 -.9],[ ],Te); Hd3=tf([ 3 -.9],[ ],Te); Hd=tf([ 3 -.9],[ ],Te); step(hd,hd4,hd3,hd);shg legend('5 zecimale','4 zecimale','3 zecimale',' zecimale') Step Response 5 zecimale 4 zecimale 3 zecimale zecimale Make a function script in order to represent graphically the poles and zeros of the pulse transfer function together with the unity circle; call the function script with one parameter represent(hd), where the parameter is the pulse transfer function. function []=reprezentare(hd) x=:.:*pi; plot(sin(x),cos(x),'r');% cercul de raza unitate [num,den]=tfdata(hd,'v'); poli=roots(den); zerouri=roots(num); hold; plot([ ],[-..],'--k',[-..],[ ],'--k') plot(real(poli),imag(poli),'x','markersize',6); plot(real(zerouri),imag(zerouri),'o','markersize',6); title('') axis([ ]) hold;shg 7

8 ^z =.9484 ^z = 5866 z _ = Compare the pulse transfer functions obtaing after discretizing the next continuous transfer functions with the specified sampling period: sampled with and sampled with clear all;clc Hc=tf(5*[ 4],conv([ ],[ 3]));Te=/; Hc=tf(5*[ 4],conv([ ],[ 3]));Te=/; Hd=cd(Hc,Te,'zoh') Hd=cd(Hc,Te,'zoh') Notice that both pulse transfer functions have the same coefficients; the difference is made only by the sampling period. 8

9 8 Problems to be solved. Reconsider the solved problems for the case of a second order transfer function with complex conjugated poles:. Make a matlab script divided in sections in order to understand the connection between the poles locations and the behavior of the system (for the step response) (z-.95) (z-.75) Regim aperiodic amortizat ^z =.95 ^z =.75 Sample time:. seconds Discrete-time zero/pole/gain model H = z^ -.9 z Regim oscilant amortizat ^z =.95+.5i ^z =.95-.5i Sample time:. seconds 3 - roots([ ]) ans = i i

10 H = z^ -.86 z + 8 Regim oscilant intretinut ^z = i ^z = i Sample time:. seconds >> roots([ -.86 ]) ans = i.93-8i z - 9 Limita de stabilitate Sample time:. seconds 8 7 ^z = z^ z +.64 Sample time:. seconds roots([ ]) ans = 3 - Regim oscilant neamortizat - - ^z = i ^z = i i i

11 9 References []. Hăngănuţ Marius. Teoria Sistemelor. vol. II, Ed. Univ. Tehnice din Cluj-Napoca (996). []. Kuo, Benjamin C. Automatic control systems. Prentice Hall PTR, 98. [3]. Leonard, Naomi Ehrich, and William S. Levine. Using MATLAB to analyze and design Control Systems. Benjamin-Cummings Publishing Co., Inc., 995. [4]. Nise, Norman S. Control systems engineering, (With CD). John Wiley & Sons, 7. [5]. Ogata, Katsuhiko, and Yanjuan Yang. Modern control engineering. Prentice Hall, 97.

16.30 Learning Objectives and Practice Problems - - Lectures 16 through 20

16.30 Learning Objectives and Practice Problems - - Lectures 16 through 20 16.30 Learning Objectives and Practice Problems - - Lectures 16 through 20 IV. Lectures 16-20 IVA : Sampling, Aliasing, and Reconstruction JVV 9.5, Lecture Notes on Shannon - Understand the mathematical

More information

Lecture 18 Stability of Feedback Control Systems

Lecture 18 Stability of Feedback Control Systems 16.002 Lecture 18 Stability of Feedback Control Systems May 9, 2008 Today s Topics Stabilizing an unstable system Stability evaluation using frequency responses Take Away Feedback systems stability can

More information

Introduction to Signals and Systems Lecture #9 - Frequency Response. Guillaume Drion Academic year

Introduction to Signals and Systems Lecture #9 - Frequency Response. Guillaume Drion Academic year Introduction to Signals and Systems Lecture #9 - Frequency Response Guillaume Drion Academic year 2017-2018 1 Transmission of complex exponentials through LTI systems Continuous case: LTI system where

More information

Concordia University. Discrete-Time Signal Processing. Lab Manual (ELEC442) Dr. Wei-Ping Zhu

Concordia University. Discrete-Time Signal Processing. Lab Manual (ELEC442) Dr. Wei-Ping Zhu Concordia University Discrete-Time Signal Processing Lab Manual (ELEC442) Course Instructor: Dr. Wei-Ping Zhu Fall 2012 Lab 1: Linear Constant Coefficient Difference Equations (LCCDE) Objective In this

More information

Basic Signals and Systems

Basic Signals and Systems Chapter 2 Basic Signals and Systems A large part of this chapter is taken from: C.S. Burrus, J.H. McClellan, A.V. Oppenheim, T.W. Parks, R.W. Schafer, and H. W. Schüssler: Computer-based exercises for

More information

Digital Video and Audio Processing. Winter term 2002/ 2003 Computer-based exercises

Digital Video and Audio Processing. Winter term 2002/ 2003 Computer-based exercises Digital Video and Audio Processing Winter term 2002/ 2003 Computer-based exercises Rudolf Mester Institut für Angewandte Physik Johann Wolfgang Goethe-Universität Frankfurt am Main 6th November 2002 Chapter

More information

Outline. Digital Control. Lecture 3

Outline. Digital Control. Lecture 3 Outline Outline Outline 1 ler Design 2 What have we talked about in MM2? Sampling rate selection Equivalents between continuous & digital Systems Outline ler Design Emulation Method for 1 ler Design

More information

Poles and Zeros of H(s), Analog Computers and Active Filters

Poles and Zeros of H(s), Analog Computers and Active Filters Poles and Zeros of H(s), Analog Computers and Active Filters Physics116A, Draft10/28/09 D. Pellett LRC Filter Poles and Zeros Pole structure same for all three functions (two poles) HR has two poles and

More information

Infinite Impulse Response Filters

Infinite Impulse Response Filters 6 Infinite Impulse Response Filters Ren Zhou In this chapter we introduce the analysis and design of infinite impulse response (IIR) digital filters that have the potential of sharp rolloffs (Tompkins

More information

Lab 1: Simulating Control Systems with Simulink and MATLAB

Lab 1: Simulating Control Systems with Simulink and MATLAB Lab 1: Simulating Control Systems with Simulink and MATLAB EE128: Feedback Control Systems Fall, 2006 1 Simulink Basics Simulink is a graphical tool that allows us to simulate feedback control systems.

More information

NH 67, Karur Trichy Highways, Puliyur C.F, Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3

NH 67, Karur Trichy Highways, Puliyur C.F, Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3 NH 67, Karur Trichy Highways, Puliyur C.F, 639 114 Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3 IIR FILTER DESIGN Structure of IIR System design of Discrete time

More information

LECTURE FOUR Time Domain Analysis Transient and Steady-State Response Analysis

LECTURE FOUR Time Domain Analysis Transient and Steady-State Response Analysis LECTURE FOUR Time Domain Analysis Transient and Steady-State Response Analysis 4.1 Transient Response and Steady-State Response The time response of a control system consists of two parts: the transient

More information

ELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet

ELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet ELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet Lecture 10: Summary Taneli Riihonen 16.05.2016 Lecture 10 in Course Book Sanjit K. Mitra, Digital Signal Processing: A Computer-Based Approach, 4th

More information

Modeling and Analysis of Systems Lecture #9 - Frequency Response. Guillaume Drion Academic year

Modeling and Analysis of Systems Lecture #9 - Frequency Response. Guillaume Drion Academic year Modeling and Analysis of Systems Lecture #9 - Frequency Response Guillaume Drion Academic year 2015-2016 1 Outline Frequency response of LTI systems Bode plots Bandwidth and time-constant 1st order and

More information

Understanding Digital Signal Processing

Understanding Digital Signal Processing Understanding Digital Signal Processing Richard G. Lyons PRENTICE HALL PTR PRENTICE HALL Professional Technical Reference Upper Saddle River, New Jersey 07458 www.photr,com Contents Preface xi 1 DISCRETE

More information

[ á{tå TÄàt. Chapter Four. Time Domain Analysis of control system

[ á{tå TÄàt. Chapter Four. Time Domain Analysis of control system Chapter Four Time Domain Analysis of control system The time response of a control system consists of two parts: the transient response and the steady-state response. By transient response, we mean that

More information

Signals and Systems Using MATLAB

Signals and Systems Using MATLAB Signals and Systems Using MATLAB Second Edition Luis F. Chaparro Department of Electrical and Computer Engineering University of Pittsburgh Pittsburgh, PA, USA AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK

More information

Implementation and Simulation of Digital Control Compensators from Continuous Compensators Using MATLAB Software

Implementation and Simulation of Digital Control Compensators from Continuous Compensators Using MATLAB Software Implementation and Simulation of Digital Control Compensators from Continuous Compensators Using MATLAB Software MAHMOUD M. EL -FANDI Electrical and Electronic Dept. University of Tripoli/Libya m_elfandi@hotmail.com

More information

System analysis and signal processing

System analysis and signal processing System analysis and signal processing with emphasis on the use of MATLAB PHILIP DENBIGH University of Sussex ADDISON-WESLEY Harlow, England Reading, Massachusetts Menlow Park, California New York Don Mills,

More information

Experimental drawing of Nyquist and Bode diagrams

Experimental drawing of Nyquist and Bode diagrams Experimental drawing of Nyquist and Bode diagrams 1 Theoretical aspects Frequency response i.e. can be plotted in several ways; two of them are (1) as a polar plot, where the phasor length is the magnitude

More information

Signals. Continuous valued or discrete valued Can the signal take any value or only discrete values?

Signals. Continuous valued or discrete valued Can the signal take any value or only discrete values? Signals Continuous time or discrete time Is the signal continuous or sampled in time? Continuous valued or discrete valued Can the signal take any value or only discrete values? Deterministic versus random

More information

UNIT II IIR FILTER DESIGN

UNIT II IIR FILTER DESIGN UNIT II IIR FILTER DESIGN Structures of IIR Analog filter design Discrete time IIR filter from analog filter IIR filter design by Impulse Invariance, Bilinear transformation Approximation of derivatives

More information

Massachusetts Institute of Technology Department of Electrical Engineering & Computer Science 6.341: Discrete-Time Signal Processing Fall 2005

Massachusetts Institute of Technology Department of Electrical Engineering & Computer Science 6.341: Discrete-Time Signal Processing Fall 2005 Massachusetts Institute of Technology Department of Electrical Engineering & Computer Science 6.341: Discrete-Time Signal Processing Fall 2005 Project Assignment Issued: Sept. 27, 2005 Project I due: Nov.

More information

Chapter 7 Filter Design Techniques. Filter Design Techniques

Chapter 7 Filter Design Techniques. Filter Design Techniques Chapter 7 Filter Design Techniques Page 1 Outline 7.0 Introduction 7.1 Design of Discrete Time IIR Filters 7.2 Design of FIR Filters Page 2 7.0 Introduction Definition of Filter Filter is a system that

More information

ECE503: Digital Filter Design Lecture 9

ECE503: Digital Filter Design Lecture 9 ECE503: Digital Filter Design Lecture 9 D. Richard Brown III WPI 26-March-2012 WPI D. Richard Brown III 26-March-2012 1 / 33 Lecture 9 Topics Within the broad topic of digital filter design, we are going

More information

Digital Signal Processing

Digital Signal Processing COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #25 Wednesday, November 5, 23 Aliasing in the impulse invariance method: The impulse invariance method is only suitable for filters with a bandlimited

More information

Lecture 17 z-transforms 2

Lecture 17 z-transforms 2 Lecture 17 z-transforms 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/3 1 Factoring z-polynomials We can also factor z-transform polynomials to break down a large system into

More information

EES42042 Fundamental of Control Systems Bode Plots

EES42042 Fundamental of Control Systems Bode Plots EES42042 Fundamental of Control Systems Bode Plots DR. Ir. Wahidin Wahab M.Sc. Ir. Aries Subiantoro M.Sc. 2 Bode Plots Plot of db Gain and phase vs frequency It is assumed you know how to construct Bode

More information

Classical Control Design Guidelines & Tools (L10.2) Transfer Functions

Classical Control Design Guidelines & Tools (L10.2) Transfer Functions Classical Control Design Guidelines & Tools (L10.2) Douglas G. MacMartin Summarize frequency domain control design guidelines and approach Dec 4, 2013 D. G. MacMartin CDS 110a, 2013 1 Transfer Functions

More information

DIGITAL FILTERS. !! Finite Impulse Response (FIR) !! Infinite Impulse Response (IIR) !! Background. !! Matlab functions AGC DSP AGC DSP

DIGITAL FILTERS. !! Finite Impulse Response (FIR) !! Infinite Impulse Response (IIR) !! Background. !! Matlab functions AGC DSP AGC DSP DIGITAL FILTERS!! Finite Impulse Response (FIR)!! Infinite Impulse Response (IIR)!! Background!! Matlab functions 1!! Only the magnitude approximation problem!! Four basic types of ideal filters with magnitude

More information

1 PeZ: Introduction. 1.1 Controls for PeZ using pezdemo. Lab 15b: FIR Filter Design and PeZ: The z, n, and O! Domains

1 PeZ: Introduction. 1.1 Controls for PeZ using pezdemo. Lab 15b: FIR Filter Design and PeZ: The z, n, and O! Domains DSP First, 2e Signal Processing First Lab 5b: FIR Filter Design and PeZ: The z, n, and O! Domains The lab report/verification will be done by filling in the last page of this handout which addresses a

More information

Answers to Problems of Chapter 4

Answers to Problems of Chapter 4 Answers to Problems of Chapter 4 The answers to the problems of this chapter are based on the use of MATLAB. Thus, if the readers have some prior elementary knowledge on it, it will be easier for them

More information

Instruction Manual for Concept Simulators. Signals and Systems. M. J. Roberts

Instruction Manual for Concept Simulators. Signals and Systems. M. J. Roberts Instruction Manual for Concept Simulators that accompany the book Signals and Systems by M. J. Roberts March 2004 - All Rights Reserved Table of Contents I. Loading and Running the Simulators II. Continuous-Time

More information

ANNA UNIVERSITY :: CHENNAI MODEL QUESTION PAPER(V-SEMESTER) B.E. ELECTRONICS AND COMMUNICATION ENGINEERING EC334 - CONTROL SYSTEMS

ANNA UNIVERSITY :: CHENNAI MODEL QUESTION PAPER(V-SEMESTER) B.E. ELECTRONICS AND COMMUNICATION ENGINEERING EC334 - CONTROL SYSTEMS ANNA UNIVERSITY :: CHENNAI - 600 025 MODEL QUESTION PAPER(V-SEMESTER) B.E. ELECTRONICS AND COMMUNICATION ENGINEERING EC334 - CONTROL SYSTEMS Time: 3hrs Max Marks: 100 Answer all Questions PART - A (10

More information

Electrical & Computer Engineering Technology

Electrical & Computer Engineering Technology Electrical & Computer Engineering Technology EET 419C Digital Signal Processing Laboratory Experiments by Masood Ejaz Experiment # 1 Quantization of Analog Signals and Calculation of Quantized noise Objective:

More information

Laboration Exercises in Digital Signal Processing

Laboration Exercises in Digital Signal Processing Laboration Exercises in Digital Signal Processing Mikael Swartling Department of Electrical and Information Technology Lund Institute of Technology revision 215 Introduction Introduction The traditional

More information

ECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015

ECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015 Purdue University: ECE438 - Digital Signal Processing with Applications 1 ECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015 1 Introduction

More information

Design IIR Band-Reject Filters

Design IIR Band-Reject Filters db Design IIR Band-Reject Filters In this post, I show how to design IIR Butterworth band-reject filters, and provide two Matlab functions for band-reject filter synthesis. Earlier posts covered IIR Butterworth

More information

EE 422G - Signals and Systems Laboratory

EE 422G - Signals and Systems Laboratory EE 422G - Signals and Systems Laboratory Lab 3 FIR Filters Written by Kevin D. Donohue Department of Electrical and Computer Engineering University of Kentucky Lexington, KY 40506 September 19, 2015 Objectives:

More information

Platypus (/platypus/)

Platypus (/platypus/) Viewing as s4251143 Platypus (/ Page 1 assignments/ studentlist/ markers/ / Problem Set 2 Problem Set 2 Question 10 Question: Q1. Graphical Convolution Consider an LTI system with impulse response h(t

More information

Project I: Phase Tracking and Baud Timing Correction Systems

Project I: Phase Tracking and Baud Timing Correction Systems Project I: Phase Tracking and Baud Timing Correction Systems ECES 631, Prof. John MacLaren Walsh, Ph. D. 1 Purpose In this lab you will encounter the utility of the fundamental Fourier and z-transform

More information

MATLAB Assignment. The Fourier Series

MATLAB Assignment. The Fourier Series MATLAB Assignment The Fourier Series Read this carefully! Submit paper copy only. This project could be long if you are not very familiar with Matlab! Start as early as possible. This is an individual

More information

Experiments #6. Convolution and Linear Time Invariant Systems

Experiments #6. Convolution and Linear Time Invariant Systems Experiments #6 Convolution and Linear Time Invariant Systems 1) Introduction: In this lab we will explain how to use computer programs to perform a convolution operation on continuous time systems and

More information

IIR Filter Design Chapter Intended Learning Outcomes: (i) Ability to design analog Butterworth filters

IIR Filter Design Chapter Intended Learning Outcomes: (i) Ability to design analog Butterworth filters IIR Filter Design Chapter Intended Learning Outcomes: (i) Ability to design analog Butterworth filters (ii) Ability to design lowpass IIR filters according to predefined specifications based on analog

More information

Pulse Code Modulation (PCM)

Pulse Code Modulation (PCM) Project Title: e-laboratories for Physics and Engineering Education Tempus Project: contract # 517102-TEMPUS-1-2011-1-SE-TEMPUS-JPCR 1. Experiment Category: Electrical Engineering >> Communications 2.

More information

SMS045 - DSP Systems in Practice. Lab 1 - Filter Design and Evaluation in MATLAB Due date: Thursday Nov 13, 2003

SMS045 - DSP Systems in Practice. Lab 1 - Filter Design and Evaluation in MATLAB Due date: Thursday Nov 13, 2003 SMS045 - DSP Systems in Practice Lab 1 - Filter Design and Evaluation in MATLAB Due date: Thursday Nov 13, 2003 Lab Purpose This lab will introduce MATLAB as a tool for designing and evaluating digital

More information

BSNL TTA Question Paper Control Systems Specialization 2007

BSNL TTA Question Paper Control Systems Specialization 2007 BSNL TTA Question Paper Control Systems Specialization 2007 1. An open loop control system has its (a) control action independent of the output or desired quantity (b) controlling action, depending upon

More information

EE 435. Lecture 16. Compensation Systematic Two-Stage Op Amp Design

EE 435. Lecture 16. Compensation Systematic Two-Stage Op Amp Design EE 435 Lecture 16 Compensation Systematic Two-Stage Op Amp Design Review from last lecture Review of Basic Concepts Pole Locations and Stability Theorem: A system is stable iff all closed-loop poles lie

More information

Lecture 3, Multirate Signal Processing

Lecture 3, Multirate Signal Processing Lecture 3, Multirate Signal Processing Frequency Response If we have coefficients of an Finite Impulse Response (FIR) filter h, or in general the impulse response, its frequency response becomes (using

More information

ECE 203 LAB 2 PRACTICAL FILTER DESIGN & IMPLEMENTATION

ECE 203 LAB 2 PRACTICAL FILTER DESIGN & IMPLEMENTATION Version 1. 1 of 7 ECE 03 LAB PRACTICAL FILTER DESIGN & IMPLEMENTATION BEFORE YOU BEGIN PREREQUISITE LABS ECE 01 Labs ECE 0 Advanced MATLAB ECE 03 MATLAB Signals & Systems EXPECTED KNOWLEDGE Understanding

More information

Sampling and Signal Processing

Sampling and Signal Processing Sampling and Signal Processing Sampling Methods Sampling is most commonly done with two devices, the sample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquires a continuous-time signal

More information

Discrete-Time Signal Processing (DTSP) v14

Discrete-Time Signal Processing (DTSP) v14 EE 392 Laboratory 5-1 Discrete-Time Signal Processing (DTSP) v14 Safety - Voltages used here are less than 15 V and normally do not present a risk of shock. Objective: To study impulse response and the

More information

Bibliography. Practical Signal Processing and Its Applications Downloaded from

Bibliography. Practical Signal Processing and Its Applications Downloaded from Bibliography Practical Signal Processing and Its Applications Downloaded from www.worldscientific.com Abramowitz, Milton, and Irene A. Stegun. Handbook of mathematical functions: with formulas, graphs,

More information

Review of Filter Types

Review of Filter Types ECE 440 FILTERS Review of Filters Filters are systems with amplitude and phase response that depends on frequency. Filters named by amplitude attenuation with relation to a transition or cutoff frequency.

More information

Comparison Methods for Converting a Spindle Plant to Discrete System

Comparison Methods for Converting a Spindle Plant to Discrete System Indonesian Journal of Electrical Engineering and Computer Science Vol., No., March 6, pp. 575 ~ 58 DOI:.59/ijeecs.v.i.pp575-58 575 Comparison Methods for Converting a Spindle Plant to Discrete System M

More information

Signal Processing Toolbox

Signal Processing Toolbox Signal Processing Toolbox Perform signal processing, analysis, and algorithm development Signal Processing Toolbox provides industry-standard algorithms for analog and digital signal processing (DSP).

More information

Module 3 : Sampling and Reconstruction Problem Set 3

Module 3 : Sampling and Reconstruction Problem Set 3 Module 3 : Sampling and Reconstruction Problem Set 3 Problem 1 Shown in figure below is a system in which the sampling signal is an impulse train with alternating sign. The sampling signal p(t), the Fourier

More information

Lecture 3 Complex Exponential Signals

Lecture 3 Complex Exponential Signals Lecture 3 Complex Exponential Signals Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/3/1 1 Review of Complex Numbers Using Euler s famous formula for the complex exponential The

More information

The period is the time required for one complete oscillation of the function.

The period is the time required for one complete oscillation of the function. Trigonometric Curves with Sines & Cosines + Envelopes Terminology: AMPLITUDE the maximum height of the curve For any periodic function, the amplitude is defined as M m /2 where M is the maximum value and

More information

Basics of Digital Filtering

Basics of Digital Filtering 4 Basics of Digital Filtering Willis J. Tompkins and Pradeep Tagare In this chapter we introduce the concept of digital filtering and look at the advantages, disadvantages, and differences between analog

More information

SYLLABUS. For B.TECH. PROGRAMME ELECTRONICS & COMMUNICATION ENGINEERING

SYLLABUS. For B.TECH. PROGRAMME ELECTRONICS & COMMUNICATION ENGINEERING SYLLABUS For B.TECH. PROGRAMME In ELECTRONICS & COMMUNICATION ENGINEERING INSTITUTE OF TECHNOLOGY UNIVERSITY OF KASHMIR ZAKURA CAMPUS SRINAGAR, J&K, 190006 Course No. Lect Tut Prac ECE5117B Digital Signal

More information

Lab 11. Speed Control of a D.C. motor. Motor Characterization

Lab 11. Speed Control of a D.C. motor. Motor Characterization Lab 11. Speed Control of a D.C. motor Motor Characterization Motor Speed Control Project 1. Generate PWM waveform 2. Amplify the waveform to drive the motor 3. Measure motor speed 4. Estimate motor parameters

More information

Pole, zero and Bode plot

Pole, zero and Bode plot Pole, zero and Bode plot EC04 305 Lecture notes YESAREKEY December 12, 2007 Authored by: Ramesh.K Pole, zero and Bode plot EC04 305 Lecture notes A rational transfer function H (S) can be expressed as

More information

Design IIR Filters Using Cascaded Biquads

Design IIR Filters Using Cascaded Biquads Design IIR Filters Using Cascaded Biquads This article shows how to implement a Butterworth IIR lowpass filter as a cascade of second-order IIR filters, or biquads. We ll derive how to calculate the coefficients

More information

Electrical and Telecommunication Engineering Technology NEW YORK CITY COLLEGE OF TECHNOLOGY THE CITY UNIVERSITY OF NEW YORK

Electrical and Telecommunication Engineering Technology NEW YORK CITY COLLEGE OF TECHNOLOGY THE CITY UNIVERSITY OF NEW YORK NEW YORK CITY COLLEGE OF TECHNOLOGY THE CITY UNIVERSITY OF NEW YORK DEPARTMENT: Electrical and Telecommunication Engineering Technology SUBJECT CODE AND TITLE: DESCRIPTION: REQUIRED TCET 4202 Advanced

More information

SECTION 7: FREQUENCY DOMAIN ANALYSIS. MAE 3401 Modeling and Simulation

SECTION 7: FREQUENCY DOMAIN ANALYSIS. MAE 3401 Modeling and Simulation SECTION 7: FREQUENCY DOMAIN ANALYSIS MAE 3401 Modeling and Simulation 2 Response to Sinusoidal Inputs Frequency Domain Analysis Introduction 3 We ve looked at system impulse and step responses Also interested

More information

Dr Ian R. Manchester

Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus 7 Root Locus 2 Assign

More information

Experiment 8: Sampling

Experiment 8: Sampling Prepared By: 1 Experiment 8: Sampling Objective The objective of this Lab is to understand concepts and observe the effects of periodically sampling a continuous signal at different sampling rates, changing

More information

Outline. Discrete time signals. Impulse sampling z-transform Frequency response Stability INF4420. Jørgen Andreas Michaelsen Spring / 37 2 / 37

Outline. Discrete time signals. Impulse sampling z-transform Frequency response Stability INF4420. Jørgen Andreas Michaelsen Spring / 37 2 / 37 INF4420 Discrete time signals Jørgen Andreas Michaelsen Spring 2013 1 / 37 Outline Impulse sampling z-transform Frequency response Stability Spring 2013 Discrete time signals 2 2 / 37 Introduction More

More information

Fall Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class

Fall Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class Fall 2018 2019 Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class Theory Problems 1. 15 pts) [Sinusoids] Define xt) as xt) = 2sin

More information

SIGNALS AND SYSTEMS: 3C1 LABORATORY 1. 1 Dr. David Corrigan Electronic and Electrical Engineering Dept.

SIGNALS AND SYSTEMS: 3C1 LABORATORY 1. 1 Dr. David Corrigan Electronic and Electrical Engineering Dept. 2012 Signals and Systems: Laboratory 1 1 SIGNALS AND SYSTEMS: 3C1 LABORATORY 1. 1 Dr. David Corrigan Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.mee.tcd.ie/ corrigad The aims of this

More information

Lab S-9: Interference Removal from Electro-Cardiogram (ECG) Signals

Lab S-9: Interference Removal from Electro-Cardiogram (ECG) Signals DSP First, 2e Signal Processing First Lab S-9: Interference Removal from Electro-Cardiogram (ECG) Signals Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab.

More information

EE422G Solution to Homework #8

EE422G Solution to Homework #8 EE4G Solution to Homework #8. MATLAB >> H = tf([ 4],[ 6 6]); >> H = tf([ ],[ - 5 5 4]); >> step(h).7 Step Response.6.5 Amplitude.4... 4 5 6 >> step(h) Time (sec).5 Step Response.5 Amplitude.5.5.5..5..5..5.4.45

More information

The Discussion of this exercise covers the following points: Angular position control block diagram and fundamentals. Power amplifier 0.

The Discussion of this exercise covers the following points: Angular position control block diagram and fundamentals. Power amplifier 0. Exercise 6 Motor Shaft Angular Position Control EXERCISE OBJECTIVE When you have completed this exercise, you will be able to associate the pulses generated by a position sensing incremental encoder with

More information

DEPARTMENT OF ELECTRONIC ENGINEERING PRACTICAL MANUAL CONTROL SYSTEMS 3 CSYS 302

DEPARTMENT OF ELECTRONIC ENGINEERING PRACTICAL MANUAL CONTROL SYSTEMS 3 CSYS 302 Name: Student number: Mark: DEPARTMENT OF ELECTRONIC ENGINEERING PRACTICAL MANUAL CONTROL SYSTEMS 3 (Process Instrumentation and Mechatronics) CSYS 30 Latest Revision: Semester 1-016 1 INTRODUCTION The

More information

Design IIR Filter using MATLAB

Design IIR Filter using MATLAB International Journal of Science, Engineering and Technology Research (IJSETR), Volume 4, Issue 2, December 25 Design IIR Filter using MATLAB RainuArya Abstract in Digital Signal Processing (DSP), most

More information

Open Access Pulse-Width Modulated Amplifier for DC Servo System and Its Matlab Simulation

Open Access Pulse-Width Modulated Amplifier for DC Servo System and Its Matlab Simulation Send Orders for Reprints to reprints@benthamscience.ae The Open Electrical & Electronic Engineering Journal, 25, 9, 625-63 625 Open Access Pulse-Width Modulated Amplifier for DC Servo System and Its Matlab

More information

System Identification & Parameter Estimation

System Identification & Parameter Estimation System Identification & Parameter Estimation Wb2301: SIPE lecture 4 Perturbation signal design Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE 3/9/2010 Delft University of Technology

More information

Digital Signal Processing

Digital Signal Processing Digital Signal Processing System Analysis and Design Paulo S. R. Diniz Eduardo A. B. da Silva and Sergio L. Netto Federal University of Rio de Janeiro CAMBRIDGE UNIVERSITY PRESS Preface page xv Introduction

More information

Contents. An introduction to MATLAB for new and advanced users

Contents. An introduction to MATLAB for new and advanced users An introduction to MATLAB for new and advanced users (Using Two-Dimensional Plots) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional

More information

1. In the command window, type "help conv" and press [enter]. Read the information displayed.

1. In the command window, type help conv and press [enter]. Read the information displayed. ECE 317 Experiment 0 The purpose of this experiment is to understand how to represent signals in MATLAB, perform the convolution of signals, and study some simple LTI systems. Please answer all questions

More information

ELT Receiver Architectures and Signal Processing Fall Mandatory homework exercises

ELT Receiver Architectures and Signal Processing Fall Mandatory homework exercises ELT-44006 Receiver Architectures and Signal Processing Fall 2014 1 Mandatory homework exercises - Individual solutions to be returned to Markku Renfors by email or in paper format. - Solutions are expected

More information

EC CONTROL SYSTEMS ENGINEERING

EC CONTROL SYSTEMS ENGINEERING 1 YEAR / SEM: II / IV EC 1256. CONTROL SYSTEMS ENGINEERING UNIT I CONTROL SYSTEM MODELING PART-A 1. Define open loop and closed loop systems. 2. Define signal flow graph. 3. List the force-voltage analogous

More information

Digital Filtering: Realization

Digital Filtering: Realization Digital Filtering: Realization Digital Filtering: Matlab Implementation: 3-tap (2 nd order) IIR filter 1 Transfer Function Differential Equation: z- Transform: Transfer Function: 2 Example: Transfer Function

More information

Produced in cooperation with. Revision: May 26, Overview

Produced in cooperation with. Revision: May 26, Overview Lab Aignment 6: Tranfer Function Analyi Reviion: May 6, 007 Produced in cooperation with www.digilentinc.com Overview In thi lab, we will employ tranfer function to determine the frequency repone and tranient

More information

THE HONG KONG POLYTECHNIC UNIVERSITY Department of Electronic and Information Engineering. EIE2106 Signal and System Analysis Lab 2 Fourier series

THE HONG KONG POLYTECHNIC UNIVERSITY Department of Electronic and Information Engineering. EIE2106 Signal and System Analysis Lab 2 Fourier series THE HONG KONG POLYTECHNIC UNIVERSITY Department of Electronic and Information Engineering EIE2106 Signal and System Analysis Lab 2 Fourier series 1. Objective The goal of this laboratory exercise is to

More information

LAB Week 7: Data Acquisition

LAB Week 7: Data Acquisition LAB Week 7: Data Acquisition Wright State University: Mechanical Engineering ME 3600L Section 01 Report and experiment by: Nicholas Smith Experiment performed on February 23, 2015 Due: March 16, 2015 Instructor:

More information

AC : FIR FILTERS FOR TECHNOLOGISTS, SCIENTISTS, AND OTHER NON-PH.D.S

AC : FIR FILTERS FOR TECHNOLOGISTS, SCIENTISTS, AND OTHER NON-PH.D.S AC 29-125: FIR FILTERS FOR TECHNOLOGISTS, SCIENTISTS, AND OTHER NON-PH.D.S William Blanton, East Tennessee State University Dr. Blanton is an associate professor and coordinator of the Biomedical Engineering

More information

Bode plot, named after Hendrik Wade Bode, is usually a combination of a Bode magnitude plot and Bode phase plot:

Bode plot, named after Hendrik Wade Bode, is usually a combination of a Bode magnitude plot and Bode phase plot: Bode plot From Wikipedia, the free encyclopedia A The Bode plot for a first-order (one-pole) lowpass filter Bode plot, named after Hendrik Wade Bode, is usually a combination of a Bode magnitude plot and

More information

THE SINUSOIDAL WAVEFORM

THE SINUSOIDAL WAVEFORM Chapter 11 THE SINUSOIDAL WAVEFORM The sinusoidal waveform or sine wave is the fundamental type of alternating current (ac) and alternating voltage. It is also referred to as a sinusoidal wave or, simply,

More information

JUNE 2014 Solved Question Paper

JUNE 2014 Solved Question Paper JUNE 2014 Solved Question Paper 1 a: Explain with examples open loop and closed loop control systems. List merits and demerits of both. Jun. 2014, 10 Marks Open & Closed Loop System - Advantages & Disadvantages

More information

Positive Feedback and Oscillators

Positive Feedback and Oscillators Physics 3330 Experiment #5 Fall 2011 Positive Feedback and Oscillators Purpose In this experiment we will study how spontaneous oscillations may be caused by positive feedback. You will construct an active

More information

Lecture 10. Lab next week: Agenda: Control design fundamentals. Proportional Control Proportional-Integral Control

Lecture 10. Lab next week: Agenda: Control design fundamentals. Proportional Control Proportional-Integral Control 264 Lab next week: Lecture 10 Lab 17: Proportional Control Lab 18: Proportional-Integral Control (1/2) Agenda: Control design fundamentals Objectives (Tracking, disturbance/noise rejection, robustness)

More information

ENGR 210 Lab 12: Sampling and Aliasing

ENGR 210 Lab 12: Sampling and Aliasing ENGR 21 Lab 12: Sampling and Aliasing In the previous lab you examined how A/D converters actually work. In this lab we will consider some of the consequences of how fast you sample and of the signal processing

More information

Brief Introduction to Signals & Systems. Phani Chavali

Brief Introduction to Signals & Systems. Phani Chavali Brief Introduction to Signals & Systems Phani Chavali Outline Signals & Systems Continuous and discrete time signals Properties of Systems Input- Output relation : Convolution Frequency domain representation

More information

EEM478-DSPHARDWARE. WEEK12:FIR & IIR Filter Design

EEM478-DSPHARDWARE. WEEK12:FIR & IIR Filter Design EEM478-DSPHARDWARE WEEK12:FIR & IIR Filter Design PART-I : Filter Design/Realization Step-1 : define filter specs (pass-band, stop-band, optimization criterion, ) Step-2 : derive optimal transfer function

More information

Recall. Sampling. Why discrete time? Why discrete time? Many signals are continuous-time signals Light Object wave CCD

Recall. Sampling. Why discrete time? Why discrete time? Many signals are continuous-time signals Light Object wave CCD Recall Many signals are continuous-time signals Light Object wave CCD Sampling mic Lens change of voltage change of voltage 2 Why discrete time? With the advance of computer technology, we want to process

More information

CHAPTER 9 FEEDBACK. NTUEE Electronics L.H. Lu 9-1

CHAPTER 9 FEEDBACK. NTUEE Electronics L.H. Lu 9-1 CHAPTER 9 FEEDBACK Chapter Outline 9.1 The General Feedback Structure 9.2 Some Properties of Negative Feedback 9.3 The Four Basic Feedback Topologies 9.4 The Feedback Voltage Amplifier (Series-Shunt) 9.5

More information

Andrea Zanchettin Automatic Control 1 AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear control systems design

Andrea Zanchettin Automatic Control 1 AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Linear control systems design Andrea Zanchettin Automatic Control 1 AUTOMATIC CONTROL Andrea M. Zanchettin, PhD Spring Semester, 2018 Linear control systems design Andrea Zanchettin Automatic Control 2 The control problem Let s introduce

More information

EE Experiment 8 Bode Plots of Frequency Response

EE Experiment 8 Bode Plots of Frequency Response EE16:Exp8-1 EE 16 - Experiment 8 Bode Plots of Frequency Response Objectives: To illustrate the relationship between a system frequency response and the frequency response break frequencies, factor powers,

More information