MATLAB Introduction and Review
|
|
- Sharleen Mills
- 6 years ago
- Views:
Transcription
1 MATLAB Introduction and Review Alireza Abouhossein Ph.D. 1
2 Introduction MATLAB: MATrix LABoratory System Analysis: SISO (single input -single output) MIMO (multiple input multiple output) Linear and non-linear Time domain and frequency domain (Bode, Nyquist) Stability (Lyapunov) Controller development PID Pole placement Optimal Control regulators (LQR LQG) Control System Toolbox: Provides tools for modeling, analysis of a system 2
3 Matlab arithmetic i. Definition of variables and evaluation of numerical expressions >> a= 3; b=5; >> a*b Ans = 15 ii. Vectors: >>A = [ ] >> a = >> a= [1:2:10] a = iii. Matrices >> A = [1 2; 3 4] A = ; operator supresses display of results on the monitor 3
4 >>A=[ ]*[4;3;2;1;0] A= 23 >> A=[ ].*[ ] A = Matrix arithmetic
5 Built-in functions Matlab has series of predefined functions >> sin(pi/3) ans = Functions "Element-wise" on Vectors >> a=[0:pi/3:pi]; sin(a) ans= >>t= 0:0.1:1 t =
6 Matrices and Vectors operators Vectors Simple arthmatic operators are multiplication, sumation and subtraction: a*b, a+b, a-b Vector transpose: transpose(a) or a Dimension: length(a) Identifying part of the vector: a(1), a(2:2:length(a)) Accessing matrix elements: A (:, 2) the second column A (1,2) element of the first row / second column matrix operators: A+A, A-A, A*A, A^2 matrix traspose (conjucate): A matrxi inverse: inv(a) determinante: det(a) eigenvalues: eig(a) rank(a) trace(a) norm(a) Finding dimention of a matrix: size(a) 6
7 Linear Control (LTI) systems Specify Linear time-invariant (LTI) systems as Transfer function sys=tf(num,den) Zero/pole/gain models sys=zpk(z,p,k) State/space models sys=ss(a,b,c,d) e.g. h = tf(1,[1 1]) % creates transfer function 1/(s+1) Type 1 + h (s+2)/s+1 To create discrete-time systems, append the sample time Ts to the previous calling sequences. sys=tf(num,den,ts) e.g. sys= zpk(0.5,[ ],1,0.05) (z-0.5)/(z+0.1)(z-0.3) : sampling time
8 System transfer functions G = tf(numg, deng) Input the transfer function in the form of numerator and denomerator. Or you can enter a system function based on zero, pole and gain. G = zpk(z, P, K) Create the continuous-time SISO transfer function: H(s)=-2s/(s-1+j)(s-1-j)(s-2) Create h(s) as a zpk object using: h = zpk(0, [1-i 1+i 2], -2); 8
9 More matrix commands Matrices: Identity matrix nxn: eye(n) Zero matrix nxm: zeros(n,m) Create array of all ones nxm: ones(n,m) Random matrix nxm random values between 0 and 1: rand(n,m) Polynomial characteristic = poly(a) where A is an n-by-n matrix returns an n+1 element row vector whose elements are the coefficients of the characteristic polynomial, (λi-a): poly(a) Clear all: removes all saved variables from the working space whos: lists the variables defined and specifies the type and size 9
10 Polynomials Assume the polynomial: s 4 +3s 3-15s 2-2s+9=0 >> pol=[ ] pol = Finding the roots: roots(pol) ans =
11 Plots Plot >> t=0:0.25:7 >> y =sin(t) >> plot(t,y) A few important commands axis[xminxmax YminYmax] grid on / grid off title, xlabel/ ylabel, Legend hold on, hold off print -depsc nomefile.eps set(gca, )
12 M-file Using Matlab editor you can create.m file Many commands and arithmetic can be written in the file. Variables can be passed between different m files. Input-out can be defined Similar to basic programming high level languages you can use if/then/else, while,end, for/end, etc. New functions may be added to MATLAB's vocabulary if they are expressed in terms of other existing functions. The commands and functions that comprise the new function must be put in a file whose name defines the name of the new function, with a filename extension of '.m'. e.g. function [mean,stdev] = stat(x) %STAT Interesting statistics. n = length(x); mean = sum(x) / n; stdev = sqrt(sum((x - mean).^2)/n); 12
13 Example: mechanical system modeling F(t) K 1 M2 M1 x 1 K 2 D 1 D 2 x 2 Remember x 1, x 2 are function of time x 1 (t) Input to the system: F(t) we show by U(s) by Lapalce transformation to frequency domain Output: displacement of the system/ response of the system to the disturbance force shown by Y(s) in this case it is desired to find x 1? G(s) =Y(s)/U(s) 13
14 Mechanical Example con s F(t) K 2 FBD M1 M2 Force balance K 1 b 1 b 2 Understand the input(s) and output(s) x 1 x 2 m1 = 10 m2 = 5 k1 = 100 Datum k2 = 50 FBD b1 = 10 b2 = 20 F(t) K 1 X 1 M1 b 1 (dx 2 /dt-dx 1 /dt) -b 1 (dx 2 /dt-dx 1 /dt) M2 K 2 ( X 2 ) b 1 (dx 2 /dt) *Don t forget inertia forces 14
15 Mechanical example cont s After combining two equations we have: (1) (2) F(s) = + 15
16 Mechanical example cont s Transfer Function: >>G= tf(num,den) >>numg=[m 2 (b 1 +b 2 ) k 2 ] >>DenG= [m 1 m 2 (m 1 (b 1 +b 2 )+b 1 m 2 ) (b 1 b 2 +k 2 m 1 +k 1 m 2 ) (k 2 b 1 +k 1 b 1 +k 1 b 2 ) k 1 k 2 ] Response of the system to step input: >>impulse(g, t) 16
17 System response to impulse input >>Figure (1) >>t = 10; >>impulse(g, t) >>grid on You can assign the results to matrices. >>[yi,t]=impulse(g, t) >> plot (t,yi) Amplitude Impulse Response Time (seconds) 17
18 System response to step input System response to step input >> figure(2) >> t = 10; % sec >> step(g, t) >> grid on Amplitude Step Response Time (seconds) 18
19 Response to a known signal Response to an input known signal: f(t) >> figure(3) >> t = 0:0.01:50; % signal for the given time interval (sec) >> f = 50*[ones(1,2500), zeros(1,2500) ones(1,1)]; % input signal >> lsim(g, f, t) >> grid on >> axis([ ]) Linear Simulation Results Time (seconds) 19
20 State-Space representation Let and So the state variables are X= [ ] T (1) (2) Isolate acceleration variable and then rewrite: The equations are in the form of, = Ax(t)+Bu(t) y = Cx(t) 20
21 Mechanical example: SS 21
22 Response to the step input Use lsim to find SS equation response to a step input >>A=[ ; ;-k1/m1 0 -b1/m1 b1/m1;0 -k2/m2 b1/m2 -(b1+b2)/m2] >>B=[0; 0; 1/m1; 0] >>C=[ ; ] >>D=0 >>T=0:0.05:10; >>U=0.2*ones(size(T)); >>[Y,X]=lsim(A,B,C,D,U,T); >>plot(t,y) 22
23 Eigenvalues/Eigenvectors [A,B,C,D]=tf2ss(b,a) [EV, EV] = eig (A) Use ss to make state-space model from the transfer function Is there any differences between using ss command and the ss you calculated by hand? Any differences in eigenvalues and eigenvectors of two state-space models? 23
24 Discrete-time/ Continuous-time Converting continuous time to discrete time >> G(z)= c2d (G, sampling_time, method ) Converting discrete time to continuous time: >> G = d2c(g_discrete, method ) Methods: 'zoh : Zero-order hold on the inputs. The control inputs are assumed piecewise constant over the sampling period. tustin : Bilinear (Tustin) approximation to the derivative. 'matched : Zero-pole matching method for SISO systems only 24
25 Bode plot is the representation of the magnitude and phase of G(j*w) bode(g), grid on G(s)= ω c 5 Frequency analysis-bode Magnitude (db) Diagram:c2d Bode Diagram System: G Frequency (rad/s): 4.98 Magnitude (db): Phase (deg) Frequency (rad/s) Note: a decibel is defined as 20*log10 ( G(j*w ) 25
26 Discretization Shannon practice αω c ω s 10αω c in this case: 5 α 10; ω s =2 /T ω c = 5, α = 5 25 ω s /T 40 Converting the continuous-time to discrete-time: >> T = 1/ ; %sec. ω s 100 rad/sec 2 / >> G_T1 =c2d(g,t, tustin ); >> G_M1 = c2d(g,t, matched ); >> Bodemag(G, b,g_t1, r,g_m1, g-- ) >> T = 1/3 0.33; %sec. ω s 20 rad/sec 2 / >> G_T1 =c2d(g,t, tustin ); >> G_M1 = c2d(g,t, matched ); >> Bodemag(G, b,g_t1, r,g_m1, g-- ) 26
27 Discretization 0 Bode Diagram 0 Bode Diagram G(s) Tustin Matched G(s) Tustin Matched Magnitude (db) Magnitude (db) Frequency (rad/s) Frequency (rad/s) T = 1/ ; %sec. ω s 100 rad/sec T = 1/3 0.33; %sec. ω s 20 rad/sec 27
28 Discretization T = 1/15 ω s 100 rad/sec T = 1/ ω s 20 rad/sec 28
29 Hard-disk read/write head controller θ J: inertia of the head assembly : 0.01 kgm^2 D: viscous damping in the bearings : Nm.(rad/sec) K: spring back constant: 10 Nm/rad K i =motor torque constant :0.05 Nm/rad Θ: angular position of the head i: input current l Laplace transform from ito θ, Recall, H=tf(num,den) 29
30 Hard disk cont s Design a digital controller that provides accurate positioning of the read/write head. The design is preformed in digital domain. First, discretize the continues plant we have plant equipped with a D/A converter (ZOH) connected to its input Ts= 0.005; %sampling period Hd= c2d(h,ts, ZOH ) Draw bode plot of both systems Hd and H bode(h,'-',hd,'--') To analyze the discrete system, plot its step response, Step(Hd) 30
31 Hard disk cont s The system oscillates quite a bit. This is probably due to very light damping. You can check this by calculating open loop poles. openloop poles of discrete model: Damp(Hd) The poles are very light equivalent damping and near the unit circle. You need to design a compensator that increases the Damping of these poles. Amplitude Step Response Time (seconds) 31
32 Hard disk cont s The simplest compensator is gain. Try root locus to select an appropriate feedback gain. >>rlocus(hd) The poles quickly leave the unit circle and go unstable Root Locus 2 >> rlocfind(oloop) 1.5 To pick a set of poles, pick a point inside the circle Imaginary Axis /T 1 /T 1 /T 0.9 /T 0.8 /T 0.7 /T 0.5 /T 0.6 /T 0.4 /T 0.3 /T /T /T /T /T 0.2 /T /T 0.6 /T 0.4 /T 0.5 /T 0.3 /T Real Axis
33 Hard disk cont s >> ddamp(poles,ts) To analyze this design, form the closed-loop system and plot the closed loop step response: >> cloop=feedback(oloop,k); >>step(cloop) 33
34 Hard disk cont s The response is fast and settles in 14 samples that is 14 * T s = 0.07 sec. You disk drive has seek time of x 10-4 Step Response Amplitude Time (seconds) 34
35 Hard disk cont s You need to introduce some lead or a compensator with some zeros D(z)=(z+a)/(z+a) With a=-0.85 and b=0 D=zpk(o.85,0,1,Ts) Oloop =Hd*D Bode(Hd, --,oloop, - ) The plot shows the compensator has shifted up the phase plot (added lead) in the fraquncy range of ω >10 rad/sec. Magnitude (db) Phase (deg) Bode Diagram Frequency (rad/s) 35
36 Let s try root locus again: >>rlocus(oloop) >>zgrid The poles stay within the unit circle Root Locus Root Locus Imaginary Axis Imaginary Axis Real Axis Without lead compensator Real Axis With lead compensator 36
37 The Singular Value Decomposition >> A=[1 2;3 4]; >> [U,S,V] = svd(a) U = S = V = Where Would you use SVD? The Singular Value Decomposition (SVD) is a widely used technique to decompose a matrix into several component matrices. E.g. when the matrix is not full rank: that means the row and columns and the matrix are linearly dependent. 37
38 Using SIMULINK First bring up the Simulink Go to continuous icon Click on Transfer function You can either enter the coeff. Of Num/Denum or enter Num/Denum as variable but values of the available should be loaded in the working space. 38
39 Simulink Cont s 39
40 Simulink Cont s 40
41 Simulink Solver Solver type Step of integration: variable or fixed step 41
42 MIMO Creating multi-input/multi-output (MIMO) cell arrays provide an ideal means to specify the resulting arrays of numerators and denominators num = {0.5,[1 1]} %1-by-2 cell array of numerators den = {[1 0],[1 2]} %1-by-2 cell array of denominators h = tf (num,den) h11 = tf(0.5,[1 0]) % 0.5/s h12 = tf ([1 1], [1 2]) % (s+1)/(s+2) H = [h11,h12] 42
43 MIMO Computer low frequncy (DC) gain of linear time invariant system K=dcgain(sys) The continous-time DC gain is the transfer function value at the frequncy s = 0. For the state space models matrices (A,B,C,D), this value is K= D CA -1 B 43
44 MIMO The discrete-time DC gain is the transfer function value at z= 1. For SS models with matrices (A,B,C,D), this value is K= D+C(I-A) -1 B The DC gain is infinite for system with integrators. 44
45 Example: MIMO H=[1 tf([1-1],[1 1 3]); tf(1,[1 1]) tf([1 2], [1-3])] Dcgain(H)? 45
46 ACTIVE Suspension: Homework Active suspension: This example concerns active suspension of a vehicle for passenger comfort. To drive the equation of motions We draw FBD of the entire system m1 = 250kg m2 = 40kg k1 = 1000N/m k2 = 20000N/m D 0 = 500Ns/m r = 0.5 m Road Surface K1 M1 M2 K2 D 0 Road radius r X 1 M1: Mass of chassis and passengers M2: Mass of Wheel carriage X 2 46
47 Active Suspension M 1 = u K 1 (x 1 -x 2 )-M 1 g- D( 1- M1 K 1 (x 1 x 2 ) - u M 1 g M 2 = K 1 (x 1 -x 2 ) + - u M 2 g K 2 (x 2 -r) K 1 (x 1 x 2 ) u M2 M 2 g K 2 (x 2 r) 47
48 Active suspension Define x 3 = 1, x 4 = can be assembled as. Then the equations = Ax + B 1 u + B 2 r + c 1 48
49 Automatic control and system theory Let s download this toolbox Books and Matlab Toolbox of Prof. Marro available at: nnimarro/gm_books.htm 49
50 Solution to the active suspension system The system can be written in form of = Ax + B 1 u + B 2 r + c 1 Define x 3 = 1, x 4 =. Then, the SS can be written as: 50
Outline. Classical Control. Lecture 1
Outline Outline Outline 1 Introduction 2 Prerequisites Block diagram for system modeling Modeling Mechanical Electrical Outline Introduction Background Basic Systems Models/Transfers functions 1 Introduction
More informationDepartment of Electrical and Computer Engineering ECED4601 Digital Control System Lab3 Digital State Space Model
Department of Electrical and Computer Engineering ECED46 Digital Control System Lab3 Digital State Space Model Objectives. To learn some MATLAB commands that deals with the discrete time systems.. To give
More informationAPPLICATIONS FOR ROBOTICS
Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table
More informationECE 3793 Matlab Project 3 Solution
ECE 3793 Matlab Project 3 Solution Spring 27 Dr. Havlicek. (a) In text problem 9.22(d), we are given X(s) = s + 2 s 2 + 7s + 2 4 < Re {s} < 3. The following Matlab statements determine the partial fraction
More informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
More informationIntroduction to Feedback Control
Introduction to Feedback Control Control System Design Why Control? Open-Loop vs Closed-Loop (Feedback) Why Use Feedback Control? Closed-Loop Control System Structure Elements of a Feedback Control System
More informationBangladesh University of Engineering and Technology. EEE 402: Control System I Laboratory
Bangladesh University of Engineering and Technology Electrical and Electronic Engineering Department EEE 402: Control System I Laboratory Experiment No. 4 a) Effect of input waveform, loop gain, and system
More information1 x(k +1)=(Φ LH) x(k) = T 1 x 2 (k) x1 (0) 1 T x 2(0) T x 1 (0) x 2 (0) x(1) = x(2) = x(3) =
567 This is often referred to as Þnite settling time or deadbeat design because the dynamics will settle in a Þnite number of sample periods. This estimator always drives the error to zero in time 2T or
More information] [ 200. ] 3 [ 10 4 s. [ ] s + 10 [ P = s [ 10 8 ] 3. s s (s 1)(s 2) series compensator ] 2. s command pre-filter [ 0.
EEE480 Exam 2, Spring 204 A.A. Rodriguez Rules: Calculators permitted, One 8.5 sheet, closed notes/books, open minds GWC 352, 965-372 Problem (Analysis of a Feedback System) Consider the feedback system
More informationOKLAHOMA STATE UNIVERSITY
OKLAHOMA STATE UNIVERSITY ECEN 4413 - Automatic Control Systems Matlab Lecture 1 Introduction and Control Basics Presented by Moayed Daneshyari 1 What is Matlab? Invented by Cleve Moler in late 1970s to
More informationECE382/ME482 Spring 2005 Homework 7 Solution April 17, K(s + 0.2) s 2 (s + 2)(s + 5) G(s) =
ECE382/ME482 Spring 25 Homework 7 Solution April 17, 25 1 Solution to HW7 AP9.5 We are given a system with open loop transfer function G(s) = K(s +.2) s 2 (s + 2)(s + 5) (1) and unity negative feedback.
More informationPositioning Servo Design Example
Positioning Servo Design Example 1 Goal. The goal in this design example is to design a control system that will be used in a pick-and-place robot to move the link of a robot between two positions. Usually
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationMAE 143B - Homework 7
MAE 143B - Homework 7 6.7 Multiplying the first ODE by m u and subtracting the product of the second ODE with m s, we get m s m u (ẍ s ẍ i ) + m u b s (ẋ s ẋ u ) + m u k s (x s x u ) + m s b s (ẋ s ẋ u
More informationLABORATORY INSTRUCTION MANUAL CONTROL SYSTEM II LAB EE 693
LABORATORY INSTRUCTION MANUAL CONTROL SYSTEM II LAB EE 693 ELECTRICAL ENGINEERING DEPARTMENT JIS COLLEGE OF ENGINEERING (AN AUTONOMOUS INSTITUTE) KALYANI, NADIA EXPERIMENT NO : CS II/ TITLE : FAMILIARIZATION
More informationMAE 143B - Homework 9
MAE 143B - Homework 9 7.1 a) We have stable first-order poles at p 1 = 1 and p 2 = 1. For small values of ω, we recover the DC gain K = lim ω G(jω) = 1 1 = 2dB. Having this finite limit, our straight-line
More informationAMME3500: System Dynamics & Control
Stefan B. Williams May, 211 AMME35: System Dynamics & Control Assignment 4 Note: This assignment contributes 15% towards your final mark. This assignment is due at 4pm on Monday, May 3 th during Week 13
More informationLecture 14 - Using the MATLAB Control System Toolbox and Simulink Friday, February 8, 2013
Today s Objectives ENGR 105: Feedback Control Design Winter 2013 Lecture 14 - Using the MATLAB Control System Toolbox and Simulink Friday, February 8, 2013 1. introduce the MATLAB Control System Toolbox
More informationDr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review
Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics
More informationECSE 4962 Control Systems Design. A Brief Tutorial on Control Design
ECSE 4962 Control Systems Design A Brief Tutorial on Control Design Instructor: Professor John T. Wen TA: Ben Potsaid http://www.cat.rpi.edu/~wen/ecse4962s04/ Don t Wait Until The Last Minute! You got
More informationECE 3793 Matlab Project 3
ECE 3793 Matlab Project 3 Spring 2017 Dr. Havlicek DUE: 04/25/2017, 11:59 PM What to Turn In: Make one file that contains your solution for this assignment. It can be an MS WORD file or a PDF file. Make
More informationMatlab for Review. NDSU Matlab Review pg 1
NDSU Matlab Review pg 1 Becoming familiar with MATLAB The console The editor The graphics windows The help menu Saving your data (diary) General environment and the console Matlab for Review Simple numerical
More informationAPPENDIX 1 MATLAB AND ANSYS PROGRAMS
APPENDIX 1 MATLAB AND ANSYS PROGRAMS This appendix lists all the MATLAB and ANSYS codes used in each chapter, along with a short description of the purpose of each. MATLAB codes have the suffix.m and the
More informationUniversity of Alberta ENGM 541: Modeling and Simulation of Engineering Systems Laboratory #7. M.G. Lipsett & M. Mashkournia 2011
ENG M 54 Laboratory #7 University of Alberta ENGM 54: Modeling and Simulation of Engineering Systems Laboratory #7 M.G. Lipsett & M. Mashkournia 2 Mixed Systems Modeling with MATLAB & SIMULINK Mixed systems
More informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More informationLABORATORY INSTRUCTION MANUAL CONTROL SYSTEM I LAB EE 593
LABORATORY INSTRUCTION MANUAL CONTROL SYSTEM I LAB EE 593 ELECTRICAL ENGINEERING DEPARTMENT JIS COLLEGE OF ENGINEERING (AN AUTONOMOUS INSTITUTE) KALYANI, NADIA CONTROL SYSTEM I LAB. MANUAL EE 593 EXPERIMENT
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING
KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK SUB.NAME : CONTROL SYSTEMS BRANCH : ECE YEAR : II SEMESTER: IV 1. What is control system? 2. Define open
More informationSTABILITY ANALYSIS. Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated using cones: Stable Neutral Unstable
ECE4510/5510: Feedback Control Systems. 5 1 STABILITY ANALYSIS 5.1: Bounded-input bounded-output (BIBO) stability Asystemmaybe stable, neutrallyormarginallystable, or unstable. This can be illustrated
More informationOutline. Classical Control. Lecture 5
Outline Outline Outline 1 What is 2 Outline What is Why use? Sketching a 1 What is Why use? Sketching a 2 Gain Controller Lead Compensation Lag Compensation What is Properties of a General System Why use?
More informationTutorial 4 (Week 11): Matlab - Digital Control Systems
ELEC 3004 Systems: Signals & Controls Tutorial 4 (Week ): Matlab - Digital Control Systems The process of designing and analysing sampled-data systems is enhanced by the use of interactive computer tools
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
More information6.1 Sketch the z-domain root locus and find the critical gain for the following systems K., the closed-loop characteristic equation is K + z 0.
6. Sketch the z-domain root locus and find the critical gain for the following systems K (i) Gz () z 4. (ii) Gz K () ( z+ 9. )( z 9. ) (iii) Gz () Kz ( z. )( z ) (iv) Gz () Kz ( + 9. ) ( z. )( z 8. ) (i)
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Guidance and Control Introduction and PID Loops Dr. Kostas Alexis (CSE) Autonomous Robot Challenges How do I control where to go? Autonomous Mobile Robot Design Topic:
More informationMechanical Systems Part A: State-Space Systems Lecture AL12
AL: 436-433 Mechanical Systems Part A: State-Space Systems Lecture AL Case study Case study AL: Design of a satellite attitude control system see Franklin, Powell & Emami-Naeini, Ch. 9. Requirements: accurate
More informationDigital Control Systems
Digital Control Systems Lecture Summary #4 This summary discussed some graphical methods their use to determine the stability the stability margins of closed loop systems. A. Nyquist criterion Nyquist
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationDESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD
206 Spring Semester ELEC733 Digital Control System LECTURE 7: DESIGN USING TRANSFORMATION TECHNIQUE CLASSICAL METHOD For a unit ramp input Tz Ez ( ) 2 ( z ) D( z) G( z) Tz e( ) lim( z) z 2 ( z ) D( z)
More informationSchool of Mechanical Engineering Purdue University. ME375 Feedback Control - 1
Introduction to Feedback Control Control System Design Why Control? Open-Loop vs Closed-Loop (Feedback) Why Use Feedback Control? Closed-Loop Control System Structure Elements of a Feedback Control System
More informationLaboratory 11 Control Systems Laboratory ECE3557. State Feedback Controller for Position Control of a Flexible Joint
Laboratory 11 State Feedback Controller for Position Control of a Flexible Joint 11.1 Objective The objective of this laboratory is to design a full state feedback controller for endpoint position control
More informationR a) Compare open loop and closed loop control systems. b) Clearly bring out, from basics, Force-current and Force-Voltage analogies.
SET - 1 II B. Tech II Semester Supplementary Examinations Dec 01 1. a) Compare open loop and closed loop control systems. b) Clearly bring out, from basics, Force-current and Force-Voltage analogies..
More informationRoot Locus. Motivation Sketching Root Locus Examples. School of Mechanical Engineering Purdue University. ME375 Root Locus - 1
Root Locus Motivation Sketching Root Locus Examples ME375 Root Locus - 1 Servo Table Example DC Motor Position Control The block diagram for position control of the servo table is given by: D 0.09 Position
More informationInverted Pendulum: State-Space Methods for Controller Design
1 de 12 18/10/2015 22:45 Tips Effects TIPS ABOUT BASICS HARDWARE INDEX NEXT INTRODUCTION CRUISE CONTROL MOTOR SPEED MOTOR POSITION SYSTEM MODELING ANALYSIS Inverted Pendulum: State-Space Methods for Controller
More informationMAS107 Control Theory Exam Solutions 2008
MAS07 CONTROL THEORY. HOVLAND: EXAM SOLUTION 2008 MAS07 Control Theory Exam Solutions 2008 Geir Hovland, Mechatronics Group, Grimstad, Norway June 30, 2008 C. Repeat question B, but plot the phase curve
More informationMassachusetts Institute of Technology Department of Mechanical Engineering Dynamics and Control II Design Project
Massachusetts Institute of Technology Department of Mechanical Engineering.4 Dynamics and Control II Design Project ACTIVE DAMPING OF TALL BUILDING VIBRATIONS: CONTINUED Franz Hover, 5 November 7 Review
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using the
More informationChapter 4: The State Space and Numerical Simulation
Chapter 4: The State Space and Numerical Simulation Samantha Ramirez Preview Questions What mathematical formulation might we apply to the models we derived in the previous chapter in order to facilitate
More informationCONTROL OF DIGITAL SYSTEMS
AUTOMATIC CONTROL AND SYSTEM THEORY CONTROL OF DIGITAL SYSTEMS Gianluca Palli Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI) Università di Bologna Email: gianluca.palli@unibo.it
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More information16.31 Homework 2 Solution
16.31 Homework Solution Prof. S. R. Hall Issued: September, 6 Due: September 9, 6 Problem 1. (Dominant Pole Locations) [FPE 3.36 (a),(c),(d), page 161]. Consider the second order system ωn H(s) = (s/p
More information1 An Overview and Brief History of Feedback Control 1. 2 Dynamic Models 23. Contents. Preface. xiii
Contents 1 An Overview and Brief History of Feedback Control 1 A Perspective on Feedback Control 1 Chapter Overview 2 1.1 A Simple Feedback System 3 1.2 A First Analysis of Feedback 6 1.3 Feedback System
More informationCYBER EXPLORATION LABORATORY EXPERIMENTS
CYBER EXPLORATION LABORATORY EXPERIMENTS 1 2 Cyber Exploration oratory Experiments Chapter 2 Experiment 1 Objectives To learn to use MATLAB to: (1) generate polynomial, (2) manipulate polynomials, (3)
More informationVALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur 603 203. DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING SUBJECT QUESTION BANK : EC6405 CONTROL SYSTEM ENGINEERING SEM / YEAR: IV / II year
More informationAutomatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: Student ID number... Signature...
Automatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: 29..23 Given and family names......................solutions...................... Student ID number..........................
More informationFrequency domain analysis
Automatic Control 2 Frequency domain analysis Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationLaboratory handouts, ME 340
Laboratory handouts, ME 340 This document contains summary theory, solved exercises, prelab assignments, lab instructions, and report assignments for Lab 4. 2014-2016 Harry Dankowicz, unless otherwise
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationUniversity of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.
University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is
More informationMatlab Controller Design. 1. Control system toolbox 2. Functions for model analysis 3. Linear system simulation 4. Biochemical reactor linearization
Matlab Controller Design. Control system toolbox 2. Functions for model analysis 3. Linear system simulation 4. Biochemical reactor linearization Control System Toolbox Provides algorithms and tools for
More informationQuanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control. DC Motor Control Trainer (DCMCT) Student Manual
Quanser NI-ELVIS Trainer (QNET) Series: QNET Experiment #02: DC Motor Position Control DC Motor Control Trainer (DCMCT) Student Manual Table of Contents 1 Laboratory Objectives1 2 References1 3 DCMCT Plant
More informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 9. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationControl Systems Engineering ( Chapter 8. Root Locus Techniques ) Prof. Kwang-Chun Ho Tel: Fax:
Control Systems Engineering ( Chapter 8. Root Locus Techniques ) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 02-760-4253 Fax:02-760-4435 Introduction In this lesson, you will learn the following : The
More informationLinear State Feedback Controller Design
Assignment For EE5101 - Linear Systems Sem I AY2010/2011 Linear State Feedback Controller Design Phang Swee King A0033585A Email: king@nus.edu.sg NGS/ECE Dept. Faculty of Engineering National University
More informationCDS 101/110: Lecture 3.1 Linear Systems
CDS /: Lecture 3. Linear Systems Goals for Today: Describe and motivate linear system models: Summarize properties, examples, and tools Joel Burdick (substituting for Richard Murray) jwb@robotics.caltech.edu,
More informationTopic # Feedback Control. State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback
Topic #17 16.31 Feedback Control State-Space Systems Closed-loop control using estimators and regulators. Dynamics output feedback Back to reality Copyright 21 by Jonathan How. All Rights reserved 1 Fall
More informationDr. 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
More informationEEL2216 Control Theory CT1: PID Controller Design
EEL6 Control Theory CT: PID Controller Design. Objectives (i) To design proportional-integral-derivative (PID) controller for closed loop control. (ii) To evaluate the performance of different controllers
More information2.010 Fall 2000 Solution of Homework Assignment 1
2. Fall 2 Solution of Homework Assignment. Compact Disk Player. This is essentially a reprise of Problems and 2 from the Fall 999 2.3 Homework Assignment 7. t is included here to encourage you to review
More informationEE Experiment 11 The Laplace Transform and Control System Characteristics
EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal
More informationLaboratory handout 5 Mode shapes and resonance
laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.04A Systems and Controls Spring 2013
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.04A Systems and Controls Spring 2013 Problem Set #4 Posted: Thursday, Mar. 7, 13 Due: Thursday, Mar. 14, 13 1. Sketch the Root
More informationAppendix 3B MATLAB Functions for Modeling and Time-domain analysis
Appendix 3B MATLAB Functions for Modeling and Time-domain analysis MATLAB control system Toolbox contain the following functions for the time-domain response step impulse initial lsim gensig damp ltiview
More informationEigenvalues and eigenvectors System Theory: electricmotor
Eigenvalues and eigenvectors System Theory: electricmotor Alireza Abouhossein Ph.D. Alireza.abouhossein@unibo.it 1 Example-1 Can you calculate A^2, A -1 -I? without multiplying A or finding the inverse
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationCHAPTER 2 TRANSFER FUNCTION ANALYSIS
. Introduction CHAPTER TRANSFER FUNCTION ANALYSIS The purpose of this chapter is to illustrate how to derive equations of motion for Multi Degree of Freedom (mdof) systems and how to solve for their transfer
More informationLinear System Theory
Linear System Theory - Laplace Transform Prof. Robert X. Gao Department of Mechanical Engineering University of Connecticut Storrs, CT 06269 Outline What we ve learned so far: Setting up Modeling Equations
More informationCONTROL SYSTEMS LABORATORY ECE311 LAB 1: The Magnetic Ball Suspension System: Modelling and Simulation Using Matlab
CONTROL SYSTEMS LABORATORY ECE311 LAB 1: The Magnetic Ball Suspension System: Modelling and Simulation Using Matlab 1 Introduction and Purpose The purpose of this experiment is to familiarize you with
More informationControl Systems I Lecture 10: System Specifications
Control Systems I Lecture 10: System Specifications Readings: Guzzella, Chapter 10 Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich November 24, 2017 E. Frazzoli (ETH) Lecture
More informationControl of Electromechanical Systems
Control of Electromechanical Systems November 3, 27 Exercise Consider the feedback control scheme of the motor speed ω in Fig., where the torque actuation includes a time constant τ A =. s and a disturbance
More informationINTRODUCTION TO DIGITAL CONTROL
ECE4540/5540: Digital Control Systems INTRODUCTION TO DIGITAL CONTROL.: Introduction In ECE450/ECE550 Feedback Control Systems, welearnedhow to make an analog controller D(s) to control a linear-time-invariant
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationImplementation Issues for the Virtual Spring
Implementation Issues for the Virtual Spring J. S. Freudenberg EECS 461 Embedded Control Systems 1 Introduction One of the tasks in Lab 4 is to attach the haptic wheel to a virtual reference position with
More informationH(s) = s. a 2. H eq (z) = z z. G(s) a 2. G(s) A B. s 2 s(s + a) 2 s(s a) G(s) 1 a 1 a. } = (z s 1)( z. e ) ) (z. (z 1)(z e at )(z e at )
.7 Quiz Solutions Problem : a H(s) = s a a) Calculate the zero order hold equivalent H eq (z). H eq (z) = z z G(s) Z{ } s G(s) a Z{ } = Z{ s s(s a ) } G(s) A B Z{ } = Z{ + } s s(s + a) s(s a) G(s) a a
More informationK(s +2) s +20 K (s + 10)(s +1) 2. (c) KG(s) = K(s + 10)(s +1) (s + 100)(s +5) 3. Solution : (a) KG(s) = s +20 = K s s
321 16. Determine the range of K for which each of the following systems is stable by making a Bode plot for K = 1 and imagining the magnitude plot sliding up or down until instability results. Verify
More informationPower System Control
Power System Control Basic Control Engineering Prof. Wonhee Kim School of Energy Systems Engineering, Chung-Ang University 2 Contents Why feedback? System Modeling in Frequency Domain System Modeling in
More informationAppendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2)
Appendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2) For all calculations in this book, you can use the MathCad software or any other mathematical software that you are familiar
More informationChapter 6 - Solved Problems
Chapter 6 - Solved Problems Solved Problem 6.. Contributed by - James Welsh, University of Newcastle, Australia. Find suitable values for the PID parameters using the Z-N tuning strategy for the nominal
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Illinois Institute of Technology Lecture 23: Drawing The Nyquist Plot Overview In this Lecture, you will learn: Review of Nyquist Drawing the Nyquist Plot Using
More informationEECS C128/ ME C134 Final Wed. Dec. 14, am. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C128/ ME C134 Final Wed. Dec. 14, 211 81-11 am Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth 1 points total. Problem Points Score 1 16 2 12
More informationExample: DC Motor Speed Modeling
Page 1 of 5 Example: DC Motor Speed Modeling Physical setup and system equations Design requirements MATLAB representation and open-loop response Physical setup and system equations A common actuator in
More information(a) Torsional spring-mass system. (b) Spring element.
m v s T s v a (a) T a (b) T a FIGURE 2.1 (a) Torsional spring-mass system. (b) Spring element. by ky Wall friction, b Mass M k y M y r(t) Force r(t) (a) (b) FIGURE 2.2 (a) Spring-mass-damper system. (b)
More informationLinear Algebra and Matrices
Linear Algebra and Matrices 4 Overview In this chapter we studying true matrix operations, not element operations as was done in earlier chapters. Working with MAT- LAB functions should now be fairly routine.
More informationIntroduction to Matlab
History of Matlab Starting Matlab Matrix operation Introduction to Matlab Useful commands in linear algebra Scripts-M file Use Matlab to explore the notion of span and the geometry of eigenvalues and eigenvectors.
More informationCDS 101/110: Lecture 3.1 Linear Systems
CDS /: Lecture 3. Linear Systems Goals for Today: Revist and motivate linear time-invariant system models: Summarize properties, examples, and tools Convolution equation describing solution in response
More informationChapter 7. Digital Control Systems
Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,
More informationECE382/ME482 Spring 2005 Homework 6 Solution April 17, (s/2 + 1) s(2s + 1)[(s/8) 2 + (s/20) + 1]
ECE382/ME482 Spring 25 Homework 6 Solution April 17, 25 1 Solution to HW6 P8.17 We are given a system with open loop transfer function G(s) = 4(s/2 + 1) s(2s + 1)[(s/8) 2 + (s/2) + 1] (1) and unity negative
More informationRecursive, Infinite Impulse Response (IIR) Digital Filters:
Recursive, Infinite Impulse Response (IIR) Digital Filters: Filters defined by Laplace Domain transfer functions (analog devices) can be easily converted to Z domain transfer functions (digital, sampled
More informationFrequency (rad/s)
. The frequency response of the plant in a unity feedback control systems is shown in Figure. a) What is the static velocity error coefficient K v for the system? b) A lead compensator with a transfer
More informationCourse roadmap. Step response for 2nd-order system. Step response for 2nd-order system
ME45: Control Systems Lecture Time response of nd-order systems Prof. Clar Radcliffe and Prof. Jongeun Choi Department of Mechanical Engineering Michigan State University Modeling Laplace transform Transfer
More informationDepartment of Electrical and Computer Engineering. EE461: Digital Control - Lab Manual
Department of Electrical and Computer Engineering EE461: Digital Control - Lab Manual Winter 2011 EE 461 Experiment #1 Digital Control of DC Servomotor 1 Objectives The objective of this lab is to introduce
More informationDepartment of Electronics and Instrumentation Engineering M. E- CONTROL AND INSTRUMENTATION ENGINEERING CL7101 CONTROL SYSTEM DESIGN Unit I- BASICS AND ROOT-LOCUS DESIGN PART-A (2 marks) 1. What are the
More information