2 Solving Ordinary Differential Equations Using MATLAB
|
|
- Martin Gaines
- 5 years ago
- Views:
Transcription
1 Penn State Erie, The Behrend College School of Engineering E E 383 Signals and Control Lab Spring 2008 Lab 3 System Responses January 31, 2008 Due: February 7, 2008 Number of Lab Periods: 1 1 Objective The objective of this laboratory exercise is to learn how to use MATLAB to solve ordinary differential equations (ODEs) and find the responses of LTI systems for impulse, step and general inputs. 2 Solving Ordinary Differential Equations Using MATLAB MATLAB has the capability to solve a wide variety of problems involving differential equations. In the following, we introduce how to use MATLAB to solve initial value ordinary differential equations. We consider the following initial value problem (IVP) which is a set of first-order differential equations in state space representation along with initial conditions. ẋ = f(x, t), x(t 0 )=x 0. (1) Here x(t) =[x 1 (t),,x n (t)] T is the state variable in R n space; x 0 =[x 10,,x n0 ] T is the initial state. Note: When an initial value problem is not specified as a set of first-order differential equations, it can be written as one by using the technique that you will learn in EE BD 410 (It will also be briefly explained during the lab). MATLAB offers seven IVP solvers: ode23, ode23s, ode23t, ode23tb, ode45, ode113 and ode15s. Among then, ode45 is typically the first solver to try on a new problem. We will use ode45 in the followings. An example of using the IVP solver to solve the classic van der Pol equation ÿ μ(1 y 2 )ẏ + y =0 (2) 1
2 is presented as follows. Here μ is a parameter greater than zero. If we choose x 1 = y, x2 =ẏ, then the van der Pol equation becomes ẋ 1 = x 2, (3) ẋ 2 = μ(1 x 2 1)x 2 x 1. (4) Now we begin solving the equations. Note that the related m files vdpol.m and solution.m can be downloaded from the course web page (or you can type them in). To solve the equations (3)-(4), they must be coded in a function m file as xdot = odefile(t,x) (see Figure 1 for the function when μ = 2). That is, the file must accept a time t and a solution x and return values for the derivatives. Note that the input arguments are t and x even when the function does not use t explicitly. Note also that the output xdot must be a column vector. function xdot=vdpol(t,x) %VDPOL van der Pol equation. % Xdot = VDPOL(t,X) % Xdot(1) = X(2) % Xdot(2) = mu*(1-x(1)^2)*x(2)-x(1) % mu = 2 mu = 2; xdot = [x(2); mu*(1-x(1)^2)*x(2)-x(1)]; Figure 1: The MATLAB function vdpol.m. Given the above ODE file, this set of ODEs is solved using the following script m file(seefigure 2) (we assume the span of t is [0, 20] and x(0) = [2 0] T ). %solution.m %Solving the van der Pol equation %by using vdpol.m. tspan = [0 20]; % time span to integrate over x0 = [2; 0]; % initial conditions (must be a column) [t,x] = ode45( vdpol,tspan,x0); % solving the equation by calling % the ODE file name and using default % settings of ode45.m Figure 2: The MATLAB function solution.m. After you have run solution.m, you can view the sizes of t and x by using size(t) and size(x) respectively. You can also plot the two states on the same plot by typing in plot(t,x(:,1),t,x(:,2), -- ). Note: In solution.m, you can specify the solution time points you desired by simply adding them to tspan, for example tspan = linspace(0,20,100);. If you are not satisfied with the default 2
3 tolerances and options given by ode45, you can use the function odeset to specify your own options. For example, in solution.m, you can use the following lines options=odeset( AbsTol,1e-8, RelTol, 1e-6); [t,x]=ode45( vdpol,tspan,x0,options); 3 Responses of LTI Systems for Impulse, Step and General Inputs Responses Using Transfer Functions From the systems point of view, an LTI system is often represented by its transfer function H(s). Note that the impulse response h(t) is simply the inverse Laplace transform of the system transfer function, H(s). As such, being able to determine and plot the system impulse response will allow one to completely characterize a system in the time domain. MATLAB has the ability to plot the system impulse response by using the impulse command. If one examines the syntax for the impulse command via the MATLAB help facility (i.e., help impulse), one sees that there are numerous ways that one can describe a system so that its impulse response can be found. In this laboratory exercise, the system will be described by the coefficients of the numerator and denominator of the transfer function entered in vector form. As an example, consider a system represented by its transfer function 3 s The m file shown in Figure 3 will then allow one to plot the impulse response in the time interval 0 t 10. Note that if the fourth line in the m file is changed to impulse(num,den,t), then the remaining lines can be eliminated since MATLAB will automatically generate its standard labeled impulse response plot if one uses impulse(num,den,t). (The complete technique in Figure 3 is only necessary if one wants to customize the way the data is presented.) t=0:0.01:10; num=[3]; den=[1 0.5]; imp_resp=impulse(num,den,t); plot(t,imp_resp) xlabel( Time, t ) ylabel( Amplitude, h(t) ) title( Example impulse response ) %define time vector %define numerator coefficients %define denominator coefficients %store impulse response data in imp_resp %plot impulse response data versus t %x axis label %y axis label %plot title Figure 3: MATLAB m file used to plot the system impulse response for 3 s+0.5. In addition to the system impulse response, MATLAB can also determine the system step response. This is important as one is often required to characterize the system step response as underdamped, overdamped or critically damped. The MATLAB command for the step response is 3
4 step and the syntax for using this command can be obtained by using help step. In the context of the impulse response m file shown in Figure 3, one only needs to replace all occurrences of the word impulse with the word step to obtain the properly labeled system step response. Finally, MATLAB can plot the response of a system for some general input. Some common general inputs, which can be directly generated by MATLAB, include sinusoids (see help sin or help cos), square waves (see help square) and various types of pulse trains (see help pulstran, help rectpuls or help tripuls). To plot the response for these general inputs, the lsim command is used (see help lsim). As an example, the m file shown Figure 4 illustrates how one can plot the response of the previously defined H(s) for a 1 Hz, 75% duty cycle pulse train input that has a peak amplitude of 1 over the time interval 0 t 10. t=0:0.01:10; num=[3]; den=[1 0.5]; in_signal=(square(2*pi*t,75)+1)/2; out_signal=lsim(num,den,in_signal,t); plot(t,out_signal, k-,t,in_signal, k: ) xlabel( Time, t ) ylabel( Amplitude, y(t) and u(t) ) title( Example pulse train response ) legend( y(t), u(t) ) %define time vector %define numerator coefficients %define denominator coefficients %define the 1Hz, 75% duty cycle pulse %train of amplitude=1 %find system output %plot system output & input versus t %x axis label %y axis label %plot title %put legend to indicate y(t) and u(t) Figure 4: MATLAB m file used to plot the response of 3 s+0.5 pulse train input with an amplitude of 1. for a 75% duty cycle, 1 Hz Responses Using Convolution Another method to obtain the response of an LTI system is by using convolution. In EE BD 326, you have learned that the response for an LTI system can be represented by the convolution y(t) = h(τ)u(t τ)dτ (5) where u(t) is the system input (do not be confused with the unit step function 1(t)) and h(t) is the impulse response of the system. The continuous-time convolution (5) can be approximated by a discrete-time convolution lim h(nδτ)u(t nδτ)δτ. (6) Δτ 0 n= Stated in words, one creates a discrete-time vector, evaluates the two functions of interest along the time vector, convolves the two vectors and then multiplies the convolution by Δτ. IfΔτ is chosen to be small enough, n= h(nδτ)u(t nδτ)δτ can provide us with a good approximation of the system response. In particular, consider the case of a causal system with input u(t) which satisfies 4
5 u(t) = 0 for any t<0. Assume we want to compute the system response y(t) for0 t T. We can divide the interval [0,T]intoN equal intervals, each with length Δτ (i.e., NΔτ = T ). The system response at the grid point t = kδτ can then be approximated by y(kδτ) = k h(nδτ)u(kδτ nδτ)δτ. (7) n=0 It can be observed that y(kδτ)asexpressedin(7)isexactlythe(k+1)th element in the convolution of the following two vectors h = [h(0),h(δτ),h(2δτ),,h(nδτ)], (8) u = [u(0),u(δτ),u(2δτ),,u(nδτ)] (9) multiplied by Δτ. Figure 5 shows a MATLAB m file which computes and plots the response (for t [0, 6]) of a system whose impulse response is h(t) =3e 0.5t 1(t) for the input signal u(t) =1(t) 1(t 2). deltatau=0.005; %define deltatau to be a small number t=0:deltatau:6; %define time vector N=length(t); %the total number of grid points h=3*exp(-0.5*t); %the impulse reponse for t in [0,6] %since we want to know the system %response for t in [0,6] and u(t)=0 %for t<0, knowing h(t), t in [0,6] suffices %our needs for the solution of this example u=zeros(1,n); %this and next lines define u(1:(n-1)/3)=ones(1,(n-1)/3); %the input u(t)=1(t)-1(t-2) %as a row vector z=conv(h,u); %z is the convolution of h and u y=z(1:n)*deltatau; %for our problem, we only need to know %the first N elements in z multiplied by %deltatau plot(t,y) %plot the response data versus t xlabel( Time, t ) %t axis label ylabel( Output, y(t) ) %y axis label title( Example system response ) %plot title Figure 5: MATLAB m file used to plot the system response using convolution. 4 Exercises All exercises are to be completed by creating, saving and running the necessary MATLAB script m files and/or function m files. 5
6 1. Determine the state solution of the following linear system [ ] [ ẋ = x, x(0) = for t [0, 6]. Plot the state trajectory (phase plot) by using the command plot(x(:,1),x(:,2)) and put appropriate labels and title. 17 points. 2. A spring-mass-damper system has the differential equation mÿ + bẏ + ky = u. By defining x 1 = y, x2 = ẏ, it can be written as a set of first order ordinary differential equations ], ẋ 1 = x 2, (10) ẋ 2 = k m x 1 b m x u. (11) m Determine x 1 (t) andx 2 (t) fort [0, 100]. Here we assume that m = 10, k =4,b =2and u(t) =0.5cos(0.5t), for t 0. The initial condition is y(0) = ẏ(0) = 0 (i.e., x 1 (0) = x 2 (0) = 0). Plot the input u(t) and the output y(t) versus time on the same plot. Observe what is the relationship between their frequencies (be sure to clearly mention this in your lab report). (Hint: you can put the expression of u(t) directly into the system state equation so that the right hand side of the system equation only depends on t and x). 17 points. 3. Plot the step response and the impulse response (use separate graphs) for 200 s 2 +30s over the time interval 0 t 1 with a step size of overdamped system.) 17 points. (This system is an example of 4. Repeat Problem 3 for 29 s 2 +4s +29 but use a time interval of 0 t 3 with a step size of (This system is an example of underdamped system.) 16 points. 5. For the system defined by s s s , use the lsim command to determine and plot the output (plot on the same graph with the input) if the input to the system is a 3 Hz square wave with a peak to peak amplitude of 2. Use a time interval of 0 t 2andastepsizeof0.01. (Assume zero initial condition for the system.) 17 points. 6
7 6. Repeat Problem 5 but use 100 cos(20t) as the input. Use a time interval of 0 t 3witha step size of (Assume zero initial condition for the system.) 16 points. 7. (Optional Problem) Use the convolution method to find and plot the response of the system (for 0 t 1) in Problem 3 when the input is a unit step input. Show your results for Δτ = 0.01 and Δτ = Compare your results with the one you obtained from Problem 3 by showing both results on the same plot and making comments. (Hint: you should first find out the system impulse response.) 10 points. Final Report Expectations for Lab 3 For this laboratory report, submit the information noted below. m files and/or items entered into MATLAB and the plots generated by MATLAB should be included in the laboratory report in computer generated form. Items that require analysis/interpretation on your part can be presented in handwritten form. Also, include a title page with the course name and number, the date the exercise was performed and your name along with the name of your lab partner(s). 1. For Problems 1-2 of the laboratory exercise, include the script m files, the function m files and the plots (labeled appropriately and titled, each curve clearly indicated). For Problem 2, be sure to also clearly indicate the relationship between the frequencies of the input and the output. 2. For Problems 3-6 of the laboratory exercise, include the script m files and the plots (labeled appropriately and titled, each curve clearly indicated). 3. For Problem 7 (if you choose to do it), include the script m files, the plots (labeled appropriately and titled, each curve clearly indicated), and your comments on the comparisons. 7
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationComputer lab for MAN460
Computer lab for MAN460 (version 20th April 2006, corrected 20 May) Prerequisites Matlab is rather user friendly, and to do the first exercises, it is enough to write an m-file consisting of two lines,
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 informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
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 informationAMS 27L LAB #8 Winter 2009
AMS 27L LAB #8 Winter 29 Solving ODE s in Matlab Objectives:. To use Matlab s ODE Solvers 2. To practice using functions and in-line functions Matlab s ODE Suite Matlab offers a suite of ODE solvers including:
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationThe roots are found with the following two statements. We have denoted the polynomial as p1, and the roots as roots_ p1.
Part II Lesson 10 Numerical Analysis Finding roots of a polynomial In MATLAB, a polynomial is expressed as a row vector of the form [an an 1 a2 a1 a0]. The elements ai of this vector are the coefficients
More informationProblem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013
EE 56: Digital Control Systems Problem Set 3: Solution Due on Mon 7 th Oct in class Fall 23 Problem For the causal LTI system described by the difference equation y k + 2 y k = x k, () (a) By first finding
More informationTherefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1
Module 2 : Signals in Frequency Domain Problem Set 2 Problem 1 Let be a periodic signal with fundamental period T and Fourier series coefficients. Derive the Fourier series coefficients of each of the
More informationMAE143A Signals & Systems - Homework 5, Winter 2013 due by the end of class Tuesday February 12, 2013.
MAE43A Signals & Systems - Homework 5, Winter 23 due by the end of class Tuesday February 2, 23. If left under my door, then straight to the recycling bin with it. This week s homework will be a refresher
More informationLesson 14: Van der Pol Circuit and ode23s
Lesson 4: Van der Pol Circuit and ode3s 4. Applied Problem. A series LRC circuit when coupled via mutual inductance with a triode circuit can generate a sequence of pulsing currents that have very rapid
More informationProblem Weight Total 100
EE 350 Problem Set 3 Cover Sheet Fall 2016 Last Name (Print): First Name (Print): ID number (Last 4 digits): Section: Submission deadlines: Turn in the written solutions by 4:00 pm on Tuesday September
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 informationEE 4314 Lab 1 Handout Control Systems Simulation with MATLAB and SIMULINK Spring Lab Information
EE 4314 Lab 1 Handout Control Systems Simulation with MATLAB and SIMULINK Spring 2013 1. Lab Information This is a take-home lab assignment. There is no experiment for this lab. You will study the tutorial
More informationSolving ODEs and PDEs in MATLAB. Sören Boettcher
16.02.2009 Introduction Quick introduction to syntax ODE in the form of Initial Value Problems (IVP) what equations can handle how to code into how to choose the right solver how to get the solver to do
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 informationHomework 1 Solutions
18-9 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 18 Homework 1 Solutions Part One 1. (8 points) Consider the DT signal given by the algorithm: x[] = 1 x[1] = x[n] = x[n 1] x[n ] (a) Plot
More informationEE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL
EE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL EE 213 Fall 2009 LABORATORY #1 INTRODUCTION TO MATLAB INTRODUCTION The purpose of this laboratory is to introduce you to Matlab and to illustrate some of its circuit
More informationMAT 275 Laboratory 7 Laplace Transform and the Symbolic Math Toolbox
Laplace Transform and the Symbolic Math Toolbox 1 MAT 275 Laboratory 7 Laplace Transform and the Symbolic Math Toolbox In this laboratory session we will learn how to 1. Use the Symbolic Math Toolbox 2.
More informationa. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules.
Lab #1 - Free Vibration Name: Date: Section / Group: Procedure Steps (from lab manual): a. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules. b. Locate the various springs and
More informationExperiment 1: Linear Regression
Experiment 1: Linear Regression August 27, 2018 1 Description This first exercise will give you practice with linear regression. These exercises have been extensively tested with Matlab, but they should
More informationHomework 5 Solutions
18-290 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 2018 Homework 5 Solutions. Part One 1. (12 points) Calculate the following convolutions: (a) x[n] δ[n n 0 ] (b) 2 n u[n] u[n] (c) 2 n u[n]
More informationSolving a RLC Circuit using Convolution with DERIVE for Windows
Solving a RLC Circuit using Convolution with DERIVE for Windows Michel Beaudin École de technologie supérieure, rue Notre-Dame Ouest Montréal (Québec) Canada, H3C K3 mbeaudin@seg.etsmtl.ca - Introduction
More informationNew Mexico State University Klipsch School of Electrical Engineering. EE312 - Signals and Systems I Spring 2018 Exam #1
New Mexico State University Klipsch School of Electrical Engineering EE312 - Signals and Systems I Spring 2018 Exam #1 Name: Prob. 1 Prob. 2 Prob. 3 Prob. 4 Total / 30 points / 20 points / 25 points /
More informationHomework 6 Solutions
8-290 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 208 Homework 6 Solutions. Part One. (2 points) Consider an LTI system with impulse response h(t) e αt u(t), (a) Compute the frequency response
More informationMAT 275 Laboratory 6 Forced Equations and Resonance
MAT 275 Laboratory 6 Forced Equations and Resonance In this laboratory we take a deeper look at second-order nonhomogeneous equations. We will concentrate on equations with a periodic harmonic forcing
More informationSolutions. Chapter 1. Problem 1.1. Solution. Simulate free response of damped harmonic oscillator
Chapter Solutions Problem. Simulate free response of damped harmonic oscillator ẍ + ζẋ + x = for different values of damping ration ζ and initial conditions. Plot the response x(t) and ẋ(t) versus t and
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 informationSystem Control Engineering 0
System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/
More informationVector Fields and Solutions to Ordinary Differential Equations using Octave
Vector Fields and Solutions to Ordinary Differential Equations using Andreas Stahel 6th December 29 Contents Vector fields. Vector field for the logistic equation...............................2 Solutions
More informationFigure 1 A linear, time-invariant circuit. It s important to us that the circuit is both linear and time-invariant. To see why, let s us the notation
Convolution In this section we consider the problem of determining the response of a linear, time-invariant circuit to an arbitrary input, x(t). This situation is illustrated in Figure 1 where x(t) is
More informationECE 3084 QUIZ 2 SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING GEORGIA INSTITUTE OF TECHNOLOGY APRIL 2, Name:
ECE 3084 QUIZ 2 SCHOOL OF ELECTRICAL AND COMPUTER ENGINEERING GEORGIA INSTITUTE OF TECHNOLOGY APRIL 2, 205 Name:. The quiz is closed book, except for one 2-sided sheet of handwritten notes. 2. Turn off
More informationNAME: 13 February 2013 EE301 Signals and Systems Exam 1 Cover Sheet
NAME: February EE Signals and Systems Exam Cover Sheet Test Duration: 75 minutes. Coverage: Chaps., Open Book but Closed Notes. One 8.5 in. x in. crib sheet Calculators NOT allowed. This test contains
More informationLab 1: Dynamic Simulation Using Simulink and Matlab
Lab 1: Dynamic Simulation Using Simulink and Matlab Objectives In this lab you will learn how to use a program called Simulink to simulate dynamic systems. Simulink runs under Matlab and uses block diagrams
More informationLAB 1: MATLAB - Introduction to Programming. Objective:
LAB 1: MATLAB - Introduction to Programming Objective: The objective of this laboratory is to review how to use MATLAB as a programming tool and to review a classic analytical solution to a steady-state
More informationProblem Weight Score Total 100
EE 350 EXAM IV 15 December 2010 Last Name (Print): First Name (Print): ID number (Last 4 digits): Section: DO NOT TURN THIS PAGE UNTIL YOU ARE TOLD TO DO SO Problem Weight Score 1 25 2 25 3 25 4 25 Total
More informationES205 Analysis and Design of Engineering Systems: Lab 1: An Introductory Tutorial: Getting Started with SIMULINK
ES205 Analysis and Design of Engineering Systems: Lab 1: An Introductory Tutorial: Getting Started with SIMULINK What is SIMULINK? SIMULINK is a software package for modeling, simulating, and analyzing
More informationLTI Systems (Continuous & Discrete) - Basics
LTI Systems (Continuous & Discrete) - Basics 1. A system with an input x(t) and output y(t) is described by the relation: y(t) = t. x(t). This system is (a) linear and time-invariant (b) linear and time-varying
More informationProblem Value
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-04 COURSE: ECE-2025 NAME: GT #: LAST, FIRST Recitation Section: Circle the date & time when your Recitation
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 informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationOne-Sided Laplace Transform and Differential Equations
One-Sided Laplace Transform and Differential Equations As in the dcrete-time case, the one-sided transform allows us to take initial conditions into account. Preliminaries The one-sided Laplace transform
More informationEE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL
EE -213 BASIC CIRCUIT ANALYSIS LAB MANUAL EE 213 Spring 2008 LABORATORY #1 INTRODUCTION TO MATLAB INTRODUCTION The purpose of this laboratory is to introduce you to Matlab and to illustrate some of its
More informationMECH : a Primer for Matlab s ode suite of functions
Objectives MECH 4-563: a Primer for Matlab s ode suite of functions. Review the fundamentals of initial value problems and why numerical integration methods are needed.. Introduce the ode suite of numerical
More informationGeorge Mason University ECE 201: Introduction to Signal Analysis Spring 2017
Assigned: March 20, 2017 Due Date: Week of April 03, 2017 George Mason University ECE 201: Introduction to Signal Analysis Spring 2017 Laboratory Project #6 Due Date Your lab report must be submitted on
More informationStatistical methods. Mean value and standard deviations Standard statistical distributions Linear systems Matrix algebra
Statistical methods Mean value and standard deviations Standard statistical distributions Linear systems Matrix algebra Statistical methods Generating random numbers MATLAB has many built-in functions
More informationSIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB
INTRODUCTION SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB Laplace transform pairs are very useful tools for solving ordinary differential equations. Most
More informationDEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010
[E2.5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010 EEE/ISE PART II MEng. BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: 2:00 hours There are FOUR
More informationCHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN
CHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN 6. Introduction Frequency Response This chapter will begin with the state space form of the equations of motion. We will use Laplace transforms to
More informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear and
More informationto have roots with negative real parts, the necessary and sufficient conditions are that:
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 543 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 7 Sebastian A. Nugroho November 6, 7 Due date of the homework is: Sunday, November 6th @ :59pm.. The following
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 informationU(s) = 0 + 0s 3 + 0s 2 + 0s + 125
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5143 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 2 Andres Rainiero Hernandez Coronado September 6, 2017 The objective of this homework is to test your understanding
More informationTemperature measurement
Luleå University of Technology Johan Carlson Last revision: July 22, 2009 Measurement Technology and Uncertainty Analysis - E7021E Lab 3 Temperature measurement Introduction In this lab you are given a
More information06/12/ rws/jMc- modif SuFY10 (MPF) - Textbook Section IX 1
IV. Continuous-Time Signals & LTI Systems [p. 3] Analog signal definition [p. 4] Periodic signal [p. 5] One-sided signal [p. 6] Finite length signal [p. 7] Impulse function [p. 9] Sampling property [p.11]
More informationOrdinary Differential Equations
Ordinary Differential Equations for Engineers and Scientists Gregg Waterman Oregon Institute of Technology c 2017 Gregg Waterman This work is licensed under the Creative Commons Attribution 4.0 International
More information2D Plotting with Matlab
GEEN 1300 Introduction to Engineering Computing Class Meeting #22 Monday, Nov. 9 th Engineering Computing and Problem Solving with Matlab 2-D plotting with Matlab Script files User-defined functions Matlab
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 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 informationEE nd Order Time and Frequency Response (Using MatLab)
EE 2302 2 nd Order Time and Frequency Response (Using MatLab) Objective: The student should become acquainted with time and frequency domain analysis of linear systems using Matlab. In addition, the student
More informationPhysics with Matlab and Mathematica Exercise #1 28 Aug 2012
Physics with Matlab and Mathematica Exercise #1 28 Aug 2012 You can work this exercise in either matlab or mathematica. Your choice. A simple harmonic oscillator is constructed from a mass m and a spring
More informationThe Laplace Transform
The Laplace Transform Generalizing the Fourier Transform The CTFT expresses a time-domain signal as a linear combination of complex sinusoids of the form e jωt. In the generalization of the CTFT to the
More informationECE382/ME482 Spring 2005 Homework 1 Solution February 10,
ECE382/ME482 Spring 25 Homework 1 Solution February 1, 25 1 Solution to HW1 P2.33 For the system shown in Figure P2.33 on p. 119 of the text, find T(s) = Y 2 (s)/r 1 (s). Determine a relationship that
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationLaboratory 10 Forced Equations and Resonance
MATLAB sessions: Laboratory 9 Laboratory Forced Equations and Resonance ( 3.6 of the Edwards/Penney text) In this laboratory we take a deeper look at second-order nonhomogeneous equations. We will concentrate
More informationLab 3: Poles, Zeros, and Time/Frequency Domain Response
ECEN 33 Linear Systems Spring 2 2-- P. Mathys Lab 3: Poles, Zeros, and Time/Frequency Domain Response of CT Systems Introduction Systems that are used for signal processing are very often characterized
More informationSymbolic Solution of higher order equations
Math 216 - Assignment 4 - Higher Order Equations and Systems of Equations Due: Monday, April 16. Nothing accepted after Tuesday, April 17. This is worth 15 points. 10% points off for being late. You may
More informationMatlab Examples. Dayalbagh Educational Institute. From the SelectedWorks of D. K. Chaturvedi Dr.
Dayalbagh Educational Institute From the SelectedWorks of D. K. Chaturvedi Dr. Winter July 10, 2012 Matlab Examples D. K. Chaturvedi, Dr., Dayalbagh Educational Institute Available at: https://works.bepress.com/dk_chaturvedi/54/
More informationGEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM. COURSE: ECE 3084A (Prof. Michaels)
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-14 COURSE: ECE 3084A (Prof. Michaels) NAME: STUDENT #: LAST, FIRST Write your name on the front page
More informationM A : Ordinary Differential Equations
M A 2 0 5 1: Ordinary Differential Equations Essential Class Notes & Graphics D 19 * 2018-2019 Sections D07 D11 & D14 1 1. INTRODUCTION CLASS 1 ODE: Course s Overarching Functions An introduction to the
More informationBfh Ti Control F Ws 2008/2009 Lab Matlab-1
Bfh Ti Control F Ws 2008/2009 Lab Matlab-1 Theme: The very first steps with Matlab. Goals: After this laboratory you should be able to solve simple numerical engineering problems with Matlab. Furthermore,
More informationProblem Set #7 Solutions Due: Friday June 1st, 2018 at 5 PM.
EE102B Spring 2018 Signal Processing and Linear Systems II Goldsmith Problem Set #7 Solutions Due: Friday June 1st, 2018 at 5 PM. 1. Laplace Transform Convergence (10 pts) Determine whether each of the
More informationUniversiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 3: DISCRETE TIME SYSTEM IN TIME DOMAIN
Universiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 3: DISCRETE TIME SYSTEM IN TIME DOMAIN Pusat Pengajian Kejuruteraan Komputer Dan Perhubungan Universiti Malaysia Perlis Discrete-Time
More information24, B = 59 24, A = 55
Math 128a - Homework 8 - Due May 2 1) Problem 8.4.4 (Page 555) Solution: As discussed in the text, the fourth-order Adams-Bashforth formula is a formula of the type x n+1 = x n + h[af n + Bf n 1 + Cf n
More informationLesson 11: Mass-Spring, Resonance and ode45
Lesson 11: Mass-Spring, Resonance and ode45 11.1 Applied Problem. Trucks and cars have springs and shock absorbers to make a comfortable and safe ride. Without good shock absorbers, the truck or car will
More informationNumerical Integration of Ordinary Differential Equations for Initial Value Problems
Numerical Integration of Ordinary Differential Equations for Initial Value Problems Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu These slides are a
More information9. Introduction and Chapter Objectives
Real Analog - Circuits 1 Chapter 9: Introduction to State Variable Models 9. Introduction and Chapter Objectives In our analysis approach of dynamic systems so far, we have defined variables which describe
More informationEXPERIMENTALLY DETERMINING THE TRANSFER FUNCTION OF A SPRING- MASS SYSTEM
EXPERIMENTALLY DETERMINING THE TRANSFER FUNCTION OF A SPRING- MASS SYSTEM Lab 8 OBJECTIVES At the conclusion of this experiment, students should be able to: Experimentally determine the best fourth order
More informationAssignment 6, Math 575A
Assignment 6, Math 575A Part I Matlab Section: MATLAB has special functions to deal with polynomials. Using these commands is usually recommended, since they make the code easier to write and understand
More informationVector Fields and Solutions to Ordinary Differential Equations using MATLAB/Octave
Vector Fields and Solutions to Ordinary Differential Equations using MATLAB/Octave Andreas Stahel 5th December 27 Contents Vector field for the logistic equation 2 Solutions of ordinary differential equations
More informationINC 341 Feedback Control Systems: Lecture 2 Transfer Function of Dynamic Systems I Asst. Prof. Dr.-Ing. Sudchai Boonto
INC 341 Feedback Control Systems: Lecture 2 Transfer Function of Dynamic Systems I Asst. Prof. Dr.-Ing. Sudchai Boonto Department of Control Systems and Instrumentation Engineering King Mongkut s University
More information2 Background: Fourier Series Analysis and Synthesis
Signal Processing First Lab 15: Fourier Series Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in the Pre-Lab section before
More informationMAT 275 Laboratory 6 Forced Equations and Resonance
MATLAB sessions: Laboratory 6 MAT 275 Laboratory 6 Forced Equations and Resonance In this laboratory we take a deeper look at second-order nonhomogeneous equations. We will concentrate on equations with
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationFall 2003 Math 308/ Numerical Methods 3.6, 3.7, 5.3 Runge-Kutta Methods Mon, 27/Oct c 2003, Art Belmonte
Fall Math 8/ Numerical Methods.6,.7,. Runge-Kutta Methods Mon, 7/Oct c, Art Belmonte Summary Geometrical idea Runge-Kutta methods numerically approximate the solution of y = f (t, y), y(a) = y by using
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52
1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n
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 informationAMS 27L LAB #6 Winter 2009
AMS 27L LAB #6 Winter 2009 Symbolically Solving Differential Equations Objectives: 1. To learn about the MATLAB Symbolic Solver 2. To expand knowledge of solutions to Diff-EQs 1 Symbolically Solving Differential
More informationHOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, Maximum score 7.0 pts.
HOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, 2010. Maximum score 7.0 pts. Three problems are to be solved in this homework assignment. The
More informationUNIVERSIDADE FEDERAL DO AMAZONAS (UFAM) MATERIAL DE AULA PRÁTICA 01 MODELAGEM E SIMULAÇÃO MATEMÁTICA DE SISTEMAS FÍSICOS
UNIVERSIDADE FEDERAL DO AMAZONAS (UFAM) IDENTIFICAÇÃO DA DISCIPLINA Professor Pierre Vilar Dantas Disciplina Laboratório de Sistemas de Controle Horário Sextas, 16h00min às 18h00min Data 17/08/2018 Turma(s)
More informationMAE143 A - Signals and Systems - Winter 11 Midterm, February 2nd
MAE43 A - Signals and Systems - Winter Midterm, February 2nd Instructions (i) This exam is open book. You may use whatever written materials you choose, including your class notes and textbook. You may
More information1. Type your solutions. This homework is mainly a programming assignment.
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K S # 6 + 7 Ahmad F. Taha October 22, 2015 READ Homework Instructions: 1. Type your solutions. This homework
More informationM A : Ordinary Differential Equations
M A 2 0 5 1: Ordinary Differential Equations Essential Class Notes & Graphics C 17 * Sections C11-C18, C20 2016-2017 1 Required Background 1. INTRODUCTION CLASS 1 The definition of the derivative, Derivative
More information1 Overview of Simulink. 2 State-space equations
Modelling and simulation of engineering systems Simulink Exercise 1 - translational mechanical systems Dr. M. Turner (mct6@sun.engg.le.ac.uk 1 Overview of Simulink Simulink is a package which runs in the
More informationEEE161 Applied Electromagnetics Laboratory 1
Dr. Milica Marković Applied Electromagnetics Laboratory page 1 EEE161 Applied Electromagnetics Laboratory 1 Instructor: Dr. Milica Marković Office: Riverside Hall 3028 Email: milica@csus.edu Web:http://gaia.ecs.csus.edu/
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 informationDIGITAL SIGNAL PROCESSING LABORATORY
L AB 5 : DISCRETE T IME SYSTEM IN TIM E DOMAIN NAME: DATE OF EXPERIMENT: DATE REPORT SUBMITTED: 1 1 THEORY Mathematically, a discrete-time system is described as an operator T[.] that takes a sequence
More informationMODELLING A MASS / SPRING SYSTEM Free oscillations, Damping, Force oscillations (impulsive and sinusoidal)
DOING PHYSICS WITH MATLAB MODELLING A MASS / SPRING SYSTEM Free oscillations, Damping, Force oscillations (impulsive and sinusoidal) Download Directory: Matlab mscripts osc_harmonic01.m The script uses
More information