System Simulation using Matlab
|
|
- Candace Copeland
- 5 years ago
- Views:
Transcription
1 EE4314 Fall 2008 System Simulation using Matlab The purpose of this laboratory work is to provide experience with the Matlab software for system simulation. The laboratory work contains a guide for solving the following general problem: Given the dynamics of a system in the form of a set of differential equations, use Matlab to find what n p m the response of the system is, i.e. x( t) R or y( t) R, to a given input, u( t) R, considering certain initial values for the states of the system, x(0) = x0. I. Simulation of dynamical systems using ode integration in Matlab The first method used to obtain in simulation the response of a system involves the use of the ode functions in Matlab. The ode -type of functions provide solutions to Ordinary Differential Equations, such as the ones which describe the systems in our consideration. We will restrict our study to the use of the ode23 Matlab function. To find out more about the ode functions type for example help ode23 in the Matlab command window and read the help file. The first paragraph in the help file describes how the ode23 function can be called and what parameters it returns: ODE23 Solve non-stiff differential equations, low order method. [TOUT,YOUT] = ODE23(ODEFUN,TSPAN,Y0) with TSPAN = [T0 TFINAL] integrates the system of differential equations y' = f(t,y) from time T0 to TFINAL with initial conditions Y0. ODEFUN is a function handle. For a scalar T and a vector Y, ODEFUN(T,Y) must return a column vector corresponding to f(t,y). Each row in the solution array YOUT corresponds to a time returned in the column vector TOUT. To obtain solutions at specific times T0,T1,...,TFINAL (all increasing or all decreasing), use TSPAN = [T0 T1... TFINAL]. The first parameter that the ode23 function requires, i.e. ODEFUN, is in fact a handle to a Matlab function which describes the dynamics of the system to be simulated, i.e. the differential equations y' = f(t,y). Say for example that we are given the system described by the differential equation: m a q + m ga cos( q ) = τ This is the differential equation which describes the dynamics of a one-link planar rotational stiff robotic arm, in an environment without friction. The q 1 is the angle of the robotic arm with respect to the horizontal, τ is the torque that is rotating the arm, a 1 is the length of the arm, m 1 is the mass of the arm and g is the acceleration of gravity. This equation has been obtained using Lagrange s equation (as discussed in class). The ode23 function can only integrate first order differential equations, and for this reason we will introduce a second state to be able to describe this system with a system of two first order differential equations. Thus, with the notations, x1 = q1, x2 = q 1, u = τ, we obtain 1
2 x 1 = x2 g 1. x 2 = cos( x1 ) + u a 2 1 m1a 1 Using these two first order differential equations we can now write the Matlab function which describes the dynamics of this system. function xdot=robot(t,x); %the x input variable of the function is in fact a vector containing the values %of the states x1 and x2; x1=x(1) and x2=x(2) g=9.8 ; a1=1 ; m1=10 ; %defines the values for the constant parameters in the system dynamics xdot= [x(2) ; -(g/a1)*cos(x(1))+u/(m1*a1^2)]; %defines a vector which is in fact the derivative of the states x1_dot and x2_dot Note that there is a variable which we have not defined. That is u, which is the input function. Say for example that we what to simulate the behavior of the system in response to a input torque -0.2t τ = 20*e cos(2 πt), which is a weighted attenuated cosine function with a period of one second. Then in our program we have to insert a line which will calculate the value of the input at each given moment in time. We now have function xdot=robot(t,x); g=9.8 ; a1=1 ; m1=10 ; u=20*exp(-t*0.2)*cos(2*pi*t); % notice that the time is the other input variable in out function xdot= [x(2) ; -(g/a1)*cos(x(1))+u/(m1*a1^2)]; At this point we have set up a function which describes the system dynamics. Save this file with the name robot.m in your working folder (you can save it directly in your personal folders). In order to simulate our system, and actually get the system response, we have to use the ode23 function. The second parameter of the ode function is the time interval during which we want to obtain the system response. Say for example that we want to see the system states during 10 seconds. Then we define the simulation time span as [0 10]. The last parameter that the ode function requires is the initial conditions of the system. Say for example that the robotic arm is in position indicated by a pi/2rad angle and 0 rad/sec angular velocity. Then Y0 is [pi/2 ; 0]. The ode23 function will return two matrices, which in the help file are named TOUT and YOUT. The first matrix is in fact the time vector which contains all the time points between 0 and 10 at which the system state, which is YOUT, was measured. In our case since the system has 2 states then YOUT will be a matrix with 2 columns (one for each state) and the same number of lines as the TOUT vector. We can now write a program that will simulate our system. (this is a separate file that the file that you created for the robot dynamics) 2
3 function main clear all; % clears the workspace close all; % closes all open figures clc; % clears the command window % these three commands are not required for your program to work tint= [0 10] ; % defines the time interval [t0 tf], we will simulate % the system for 10 seconds x0= [pi/2 0]' ; % initial conditions [t,x]= ode23( robot, tint, x0); plot(t,x); % plots the system states versus time But this type of plot does not give much information to an untrained eye. To get a nicer visual result one can add the following to the main program a1=1; x1=cos(x(:,1))*a1; y1=sin(x(:,1))*a1; figure(3); for i=1:length(x1), figure(3); plot(0,0,'o',x1(i),y1(i),'o',[0 x1(i)],[0 y1(i)],'-r'); axis([ ]); pause(t(i+1)-t(i)); end The answers to the following questions should be turned in at the end of the class. Use the provided page at the end of this document. A. Explain, in the provided space on the last page, what is the effect of each of the code lines in the above program section. B. Make a simulation of the system considering 0 input (u(t)=0) and initial conditions given as [0 0] (which is: zero angle and zero angular velocity). Explain the behavior of the system. II. Simulation of dynamical systems using block diagrams in Matlab Simulink To simulate the dynamics of the same one link robot system one can also use the Simulink environment and the block diagram description of the system. Start by typing simulink in the Matlab command window. The Simulink library browser opens. Create a new simulink file. Start by inserting as many integrators as the number of states of the system. The simulink file will look like the one in Figure 1 (as the system has two states). 3
4 Figure 1 The output of the first integrator is x 1 and the output of the second integrator is x 2. Knowing that x 1 = x2, one can simply add the connection line between the input of the first integrator and the output of the second one (as presented in Figure 2) Figure 2 g 1 x x u a m a The input of the second integrator is given by the relation 2 = cos( 1) Note that it depends on the input signal u and the output of the first integrator. First we will group these two signals x 1, u in a vector of signals [ x1 u ]. For this you can use the Mux block which can be found in the signal routing library in Simulink (as per Figure 3). This new vector of signals will serve as input to a block which will have to calculate the value of ẋ 2 g 1 x x u a m a according to the equation 2 = cos( 1) You can find such a block in the user defined functions tab of the simulink library. You now obtain the schematic in Figure 4. 4
5 Figure 3 Figure 4 To introduce the function that the Fcn block has to perform you can double-click on the Fcn block and insert the appropriate function in the window that appeared. Figure 5 In this case the relation to be inserted is (-9.8/1)*cos(u(1))+(1/(10*1^2))*u(2). In this relation u(1) carries the information on x1(t) and u(2) carries the information on u(t). We now need to define our input signal u( t ). One can define this input signal in many different ways using a number of combinations of Simulink blocks. One of them is given in Figure 6, where the function block denoted Fcn1 performs the following calculation depending on its input 20*exp(-0.2*u(1))*cos(2*pi*u(1)). 5
6 Figure 6 This simulink schematic for simulating our robotic system is now complete. We only need to visualize the result. We thus insert in the schematic a Scope block which will present the two states. Figure 7 You can now run the simulation by hitting the play button in the menu. It will simulate the system behavior for 10 seconds. If you double-click on the scope block, a window will appear which will present the evolution of the two states of the system; just like the following line of Matlab code did. plot(t,x); % plots the system states versus time You can now save this data from Simulink and use it for plot in Matlab. If you hit the parameters button (second button) in the scope menu the following window will appear (which has 2 tabs). Figure 8 6
7 In the data history tab select the save data to workspace tick. In this case, after you run the simulation one more time, in the Matlab workspace will appear a data structure called ScopeData. To access the time vector you use the name ScopeData.time and to access the values of the states you use the name ScopeData.signals.values. If now you run the following code you will plot the result, which was presented in the Scope, in a Matlab figure. x=scopedata.signals.values; t=scopedata.time; plot(t,x); Make a simulation of the system considering 0 input (u(t)=0) and initial conditions given as [0 0] (which is: zero angle and zero angular velocity). To introduce the initial states in the simulink block diagram, double-click the two integrators and introduce in the appropriate field the values for the initial states. Having understood how the ode functions and simulink block diagrams can be used to simulate a system, given its differential equation dynamics, now you can work on the homework problems. 7
8 Name: A. Explain, in the provided space, what each of the lines of code in this program section is doing. (Comment the code.) function main tint= [0 10] ; % defines the time interval [t0 tf], we will simulate % the system for 10 seconds x0= [pi/2 0]' ; % initial conditions [t,x]= ode23( robot, tint, x0); plot(t,x); % plots the system states versus time a1=3; x1=cos(x(:,1))*a1; y1=sin(x(:,1))*a1; figure(3); for i=1:length(x1), figure(3); plot(0,0,'o',x1(i),y1(i),'o',[0 x1(i)],[0 y1(i)],'-r'); axis([ ]); pause(t(i+1)-t(i)); end B. Explain the behavior of the one link robotic system in response to 0 input (u(t)=0) and initial conditions given as [0 0] (which is: zero angle and zero angular velocity). 8
ROBOTICS Laboratory Problem 02
ROBOTICS 2015-2016 Laboratory Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with or without friction along the horizontal
More informationMSMS Matlab Problem 02
MSMS 2014-2015 Matlab Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with friction along the horizontal rail; the cart supports
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 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 informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationSTATE VARIABLE (SV) SYSTEMS
Copyright F.L. Lewis 999 All rights reserved Updated:Tuesday, August 05, 008 STATE VARIABLE (SV) SYSTEMS A natural description for dynamical systems is the nonlinear state-space or state variable (SV)
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 informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationSimulink Tutorial 1 CPE562
Simulink Tutorial 1 CPE562 Week 1 Introduction to Simulink Familiarization with Simulink blocks Sources: Constants Sinks: Display Operations: Sum, Product, Add, Divide. Mathematical operations involving
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 informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 1 Dianne P. O Leary c 2008 1 The
More informationLaboratory 1. Solving differential equations with nonzero initial conditions
Laboratory 1 Solving differential equations with nonzero initial conditions 1. Purpose of the exercise: - learning symbolic and numerical methods of differential equations solving with MATLAB - using Simulink
More informationLinear System Theory
Linear System Theory - Introduction to Simulink Prof. Robert X. Gao Electromechanical Systems Laboratory Department of Mechanical Engineering Outline Block Diagram Introduction Launching Simulink Modeling
More informationOrdinary Differential Equations (ODEs) Background. Video 17
Ordinary Differential Equations (ODEs) Background Video 17 Daniel J. Bodony Department of Aerospace Engineering University of Illinois at Urbana-Champaign In this video you will learn... 1 What ODEs are
More information(t) ), ( ) where t is time, T is the reaction time, u n is the position of the nth car, and the "sensitivity coefficient" λ may depend on un 1(
Page 1 A Model of Traffic Flow Everyone has had the experience of sitting in a traffic jam, or of seeing cars bunch up on a road for no apparent good reason MATLAB and SIMULINK are good tools for studying
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 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 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 informationMATLAB BASICS. Instructor: Prof. Shahrouk Ahmadi. TA: Kartik Bulusu
MATLAB BASICS Instructor: Prof. Shahrouk Ahmadi 1. What are M-files TA: Kartik Bulusu M-files are files that contain a collection of MATLAB commands or are used to define new MATLAB functions. For the
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 informationNumerical solution of ODEs
Péter Nagy, Csaba Hős 2015. H-1111, Budapest, Műegyetem rkp. 3. D building. 3 rd floor Tel: 00 36 1 463 16 80 Fax: 00 36 1 463 30 91 www.hds.bme.hu Table of contents Homework Introduction to Matlab programming
More informationMatlab Course. Anna Kristine Wåhlin. Department of Geophysics, University of Oslo. January Matlab Course p.1/??
Matlab Course Anna Kristine Wåhlin Department of Geophysics, University of Oslo January 2003 Matlab Course p.1/?? Numerical estimate of the derivative An estimate of the time derivative of dataset at time
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE201: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE21: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationMATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction
MATLAB Ordinary Differential Equation (ODE) solver for a simple example 1. Introduction Differential equations are a convenient way to express mathematically a change of a dependent variable (e.g. concentration
More informationUsing MATLAB with the Convolution Method
ECE 350 Linear Systems I MATLAB Tutorial #5 Using MATLAB with the Convolution Method A linear system with input, x(t), and output, y(t), can be described in terms of its impulse response, h(t). x(t) h(t)
More informationLab 5a: Pole Placement for the Inverted Pendulum
Lab 5a: Pole Placement for the Inverted Pendulum November 1, 2011 1 Purpose The objective of this lab is to achieve simultaneous control of both the angular position of the pendulum and horizontal position
More informationLab 1g: Horizontally Forced Pendulum & Chaotic Motion
58:080 Experimental Engineering OBJECTIVE Lab 1g: Horizontally Forced Pendulum & Chaotic Motion The objective of this lab is to study horizontally forced oscillations of a pendulum. This will be done trough
More informationBMEN 398: MATLAB Module: Higher Order Differential Equations; SIMULINK Fall 2005 Updated 8/21/2005 Hart
BMEN 398: MATLAB Module: Higher Order Differential Equations; SIMULINK Fall 2005 Updated 8/21/2005 Hart Higher Order ODEs: (Rao, 2002) Although we can now write MATLAB code to find numerical solutions
More informationModel of a DC Generator Driving a DC Motor (which propels a car)
Model of a DC Generator Driving a DC Motor (which propels a car) John Hung 5 July 2011 The dc is connected to the dc as illustrated in Fig. 1. Both machines are of permanent magnet type, so their respective
More informationEXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE (V_3)
TA name Lab section Date TA Initials (on completion) Name UW Student ID # Lab Partner(s) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE (V_3) 121 Textbook Reference: Knight, Chapter 13.1-3, 6. SYNOPSIS In
More informationStability and Frequency Response of Linear Systems
ECE 350 Linear Systems I MATLAB Tutorial #4 Stability and Frequency Response of Linear Systems This tutorial describes the MATLAB commands that can be used to determine the stability and/or frequency response
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 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 informationLab 6a: Pole Placement for the Inverted Pendulum
Lab 6a: Pole Placement for the Inverted Pendulum Idiot. Above her head was the only stable place in the cosmos, the only refuge from the damnation of the Panta Rei, and she guessed it was the Pendulum
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 informationLab 3: Quanser Hardware and Proportional Control
Lab 3: Quanser Hardware and Proportional Control The worst wheel of the cart makes the most noise. Benjamin Franklin 1 Objectives The goal of this lab is to: 1. familiarize you with Quanser s QuaRC tools
More informationSolving Differential Equations Using MATLAB
Solving Differential Equations Using MATLAB Abraham Asfaw aasfaw.student@manhattan.edu November 28, 2011 1 Introduction In this lecture, we will follow up on lecture 2 with a discussion of solutions to
More informationSimulink Modeling Tutorial
Simulink Modeling Tutorial Train system Free body diagram and Newton's law Model Construction Running the Model Obtaining MATLAB Model In Simulink, it is very straightforward to represent a physical system
More informationRotary Motion Servo Plant: SRV02. Rotary Experiment #11: 1-DOF Torsion. 1-DOF Torsion Position Control using QuaRC. Student Manual
Rotary Motion Servo Plant: SRV02 Rotary Experiment #11: 1-DOF Torsion 1-DOF Torsion Position Control using QuaRC Student Manual Table of Contents 1. INTRODUCTION...1 2. PREREQUISITES...1 3. OVERVIEW OF
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 information2 Solving Ordinary Differential Equations Using MATLAB
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
More informationENGG 5402 Course Project: Simulation of PUMA 560 Manipulator
ENGG 542 Course Project: Simulation of PUMA 56 Manipulator ZHENG Fan, 115551778 mrzhengfan@gmail.com April 5, 215. Preface This project is to derive programs for simulation of inverse dynamics and control
More informationRotary Motion Servo Plant: SRV02. Rotary Experiment #01: Modeling. SRV02 Modeling using QuaRC. Student Manual
Rotary Motion Servo Plant: SRV02 Rotary Experiment #01: Modeling SRV02 Modeling using QuaRC Student Manual SRV02 Modeling Laboratory Student Manual Table of Contents 1. INTRODUCTION...1 2. PREREQUISITES...1
More informationExperiment A11 Chaotic Double Pendulum Procedure
AME 21216: Lab I Fall 2017 Experiment A11 Chaotic Double Pendulum Procedure Deliverables: Checked lab notebook, plots with captions Background Measuring and controlling the angular position and velocity
More informationSystem Parameters and Frequency Response MAE 433 Spring 2012 Lab 2
System Parameters and Frequency Response MAE 433 Spring 2012 Lab 2 Prof. Rowley, Prof. Littman AIs: Brandt Belson, Jonathan Tu Technical staff: Jonathan Prévost Princeton University Feb. 21-24, 2012 1
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/ME/AE324: Dynamical Systems. Chapter 4: Block Diagrams
EE/ME/AE324: Dynamical Systems Chapter 4: Block Diagrams and Computer Simulation Block Diagrams A block diagram is an interconnection of: Blocks representing math operations Wires representing signals
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 informationTime Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors
Time Response of Dynamic Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Multi-dimensional trajectories Numerical integration Linear and nonlinear systems Linearization
More informationGraphical User Interface (GUI) for Torsional Vibration Analysis of Rotor Systems Using Holzer and MatLab Techniques
Basrah Journal for Engineering Sciences, vol. 14, no. 2, 2014 255 Graphical User Interface (GUI) for Torsional Vibration Analysis of Rotor Systems Using Holzer and MatLab Techniques Dr. Ameen Ahmed Nassar
More informationMotion with Constant Acceleration
Motion with Constant Acceleration INTRODUCTION Newton s second law describes the acceleration of an object due to an applied net force. In this experiment you will use the ultrasonic motion detector to
More informationSRV02-Series Rotary Experiment # 1. Position Control. Student Handout
SRV02-Series Rotary Experiment # 1 Position Control Student Handout SRV02-Series Rotary Experiment # 1 Position Control Student Handout 1. Objectives The objective in this experiment is to introduce the
More information1. Consider the 1-DOF system described by the equation of motion, 4ẍ+20ẋ+25x = f.
Introduction to Robotics (CS3A) Homework #6 Solution (Winter 7/8). Consider the -DOF system described by the equation of motion, ẍ+ẋ+5x = f. (a) Find the natural frequency ω n and the natural damping ratio
More informationDC-motor modelling and parameter identification
DC-motor modelling and parameter identification This version: November 1, 2017 Name: LERTEKNIK REG P-number: Date: AU T O MA RO TI C C O N T L Passed: LINKÖPING Chapter 1 Introduction The purpose of this
More informationMEAM 520. More Velocity Kinematics
MEAM 520 More Velocity Kinematics Katherine J. Kuchenbecker, Ph.D. General Robotics, Automation, Sensing, and Perception Lab (GRASP) MEAM Department, SEAS, University of Pennsylvania Lecture 12: October
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
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 informationMotion II. Goals and Introduction
Motion II Goals and Introduction As you have probably already seen in lecture or homework, and if you ve performed the experiment Motion I, it is important to develop a strong understanding of how to model
More informationFigure 12.1: A simple pendulum
Chapter 12 A Simple Pendulum by Brian Patterson In this module you will use DIYModeling to build a simulation of a simple pendulum. The basic ideas can be extended to other types of pendulums, such as
More information1 M62 M62.1 CONSERVATION OF ANGULAR MOMENTUM FOR AN INELASTIC COLLISION
1 M62 M62.1 CONSERVATION OF ANGULAR MOMENTUM FOR AN INELASTIC COLLISION PRELAB: Before coming to the lab, you must write the Object and Theory sections of your lab report and include the Data Tables. You
More informationLab Partner(s) TA Initials (on completion) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE
TA name Lab section Date TA Initials (on completion) Name UW Student ID # Lab Partner(s) EXPERIMENT 7: ANGULAR KINEMATICS AND TORQUE 117 Textbook Reference: Walker, Chapter 10-1,2, Chapter 11-1,3 SYNOPSIS
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 informationMathematical Modelling Using Simulink
Experiment Two Mathematical Modelling Using Simulink Control Systems Laboratory Dr. Zaer Abo Hammour Dr. Zaer Abo Hammour Control Systems Laboratory 1. Mathematical Model Definition A mathematical model
More informationWILEY. Differential Equations with MATLAB (Third Edition) Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg
Differential Equations with MATLAB (Third Edition) Updated for MATLAB 2011b (7.13), Simulink 7.8, and Symbolic Math Toolbox 5.7 Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg All
More informationGeneral Physics I Lab. M1 The Atwood Machine
Purpose General Physics I Lab In this experiment, you will learn the basic operation of computer interfacing and use it in an experimental study of Newton s second law. Equipment and components Science
More information-ASTR 204 Application of Astroimaging Techniques
Lab 5 - JPL Ephemeris, Binary Maker 3 In Lab 5 we ll look at another ephemeris generating tool; Horizons Web-Interface from JPL, and Binary Maker 3 program for generating radial velocity curves and 3-D
More informationPhysics 8 Monday, October 28, 2013
Physics 8 Monday, October 28, 2013 Turn in HW8 today. I ll make them less difficult in the future! Rotation is a hard topic. And these were hard problems. HW9 (due Friday) is 7 conceptual + 8 calculation
More informationNewton's 2 nd Law. . Your end results should only be interms of m
Newton's nd Law Introduction: In today's lab you will demonstrate the validity of Newton's Laws in predicting the motion of a simple mechanical system. The system that you will investigate consists of
More informationVideo Analysis Inertial and non-inertial reference frames
Video Analysis Inertial and non-inertial reference frames Apparatus Tracker software (free; download from http://www.cabrillo.edu/ dbrown/tracker/) video: two-carts.mov from http://physics.highpoint.edu/
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More information(W: 12:05-1:50, 50-N202)
2016 School of Information Technology and Electrical Engineering at the University of Queensland Schedule of Events Week Date Lecture (W: 12:05-1:50, 50-N202) 1 27-Jul Introduction 2 Representing Position
More informationApplying MapCalc Map Analysis Software
Applying MapCalc Map Analysis Software Generating Surface Maps from Point Data: A farmer wants to generate a set of maps from soil samples he has been collecting for several years. Previously, he would
More informationIntroduction to MATLAB
Introduction to MATLAB Violeta Ivanova, Ph.D. Educational Technology Consultant MIT Academic Computing violeta@mit.edu http://web.mit.edu/violeta/www/iap2006 Topics MATLAB Interface and Basics Linear Algebra
More informationState Feedback Controller for Position Control of a Flexible Link
Laboratory 12 Control Systems Laboratory ECE3557 Laboratory 12 State Feedback Controller for Position Control of a Flexible Link 12.1 Objective The objective of this laboratory is to design a full state
More informationTUTORIAL: STATE VARIABLES and MATLAB
TUTORIAL TUTORIAL: STATE VARIABLES and MATLAB Time-domain analysis of circuits with more than one L and C is difficult because it requires solution of characteristic equations higher than second degree.
More informationThe University of Hong Kong Department of Physics. Physics Laboratory PHYS3350 Classical Mechanics Experiment No The Physical Pendulum Name:
The University of Hong Kong Department of Physics Physics Laboratory PHYS3350 Classical Mechanics Experiment No. 3350-2 The Physical Pendulum Name: University No: Introduction One of the practical uses
More informationEE 474 Lab Part 2: Open-Loop and Closed-Loop Control (Velocity Servo)
Contents EE 474 Lab Part 2: Open-Loop and Closed-Loop Control (Velocity Servo) 1 Introduction 1 1.1 Discovery learning in the Controls Teaching Laboratory.............. 1 1.2 A Laboratory Notebook...............................
More informationUsing Simulink to analyze 2 degrees of freedom system
Using Simulink to analyze 2 degrees of freedom system Nasser M. Abbasi Spring 29 page compiled on June 29, 25 at 4:2pm Abstract A two degrees of freedom system consisting of two masses connected by springs
More informationFinite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals
Finite Element Modules for Enhancing Undergraduate Transport Courses: Application to Fuel Cell Fundamentals Originally published in 007 American Society for Engineering Education Conference Proceedings
More informationTIphysics.com. Physics. Pendulum Explorations ID: By Irina Lyublinskaya
Pendulum Explorations ID: 17 By Irina Lyublinskaya Time required 90 minutes Topic: Circular and Simple Harmonic Motion Explore what factors affect the period of pendulum oscillations. Measure the period
More information1. Write the symbolic representation and one possible unit for angular velocity, angular acceleration, torque and rotational inertia.
ROTATIONAL DYNAMICS Pre-Lab Questions Page Name: Class: Roster Number: Instructor: 1. Write the symbolic representation and one possible unit for angular velocity, angular acceleration, torque and rotational
More informationChapter 2: Numeric, Cell, and Structure Arrays
Chapter 2: Numeric, Cell, and Structure Arrays Topics Covered: Vectors Definition Addition Multiplication Scalar, Dot, Cross Matrices Row, Column, Square Transpose Addition Multiplication Scalar-Matrix,
More informationExperiment P28: Conservation of Linear and Angular Momentum (Smart Pulley)
PASCO scientific Physics Lab Manual: P28-1 Experiment P28: Conservation of Linear and Angular Momentum (Smart Pulley) Concept Time SW Interface Macintosh File Windows File rotational motion 45 m 500 or
More informationDP Physics Torque Simulation
DP Physics Torque Simulation Name Go to Phet Simulation: ( http://phet.colorado.edu/simulations/sims.php?sim=torque ) Part I: Torque 1. Click the tab at the top that says torque 2. Set the force equal
More informationAE2610 Introduction to Experimental Methods in Aerospace DYNAMIC RESPONSE OF A 3-DOF HELICOPTER MODEL
AE2610 Introduction to Experimental Methods in Aerospace DYNAMIC RESPONSE OF A 3-DOF HELICOPTER MODEL Objectives The primary objective of this experiment is to introduce the student to the measurement
More information2.004 Dynamics and Control II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.004 Dynamics and Control II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts Institute
More informationLab 1 Uniform Motion - Graphing and Analyzing Motion
Lab 1 Uniform Motion - Graphing and Analyzing Motion Objectives: < To observe the distance-time relation for motion at constant velocity. < To make a straight line fit to the distance-time data. < To interpret
More informationOctave. Tutorial. Daniel Lamprecht. March 26, Graz University of Technology. Slides based on previous work by Ingo Holzmann
Tutorial Graz University of Technology March 26, 2012 Slides based on previous work by Ingo Holzmann Introduction What is? GNU is a high-level interactive language for numerical computations mostly compatible
More informationTutorial. Getting started. Sample to Insight. March 31, 2016
Getting started March 31, 2016 Sample to Insight CLC bio, a QIAGEN Company Silkeborgvej 2 Prismet 8000 Aarhus C Denmark Telephone: +45 70 22 32 44 www.clcbio.com support-clcbio@qiagen.com Getting started
More informationAppendix A Prototypes Models
Appendix A Prototypes Models This appendix describes the model of the prototypes used in Chap. 3. These mathematical models can also be found in the Student Handout by Quanser. A.1 The QUANSER SRV-02 Setup
More informationMecanum-Wheeled Vehicle Base
Mecanum-Wheeled Vehicle Base Dan Fisher in partial fulfillment of ETLS 789 Simulation and Visualization of Dynamic Systems University of St. Thomas Dr. Michael Hennessey Fall 2014 Due December 18, 2014
More informationLecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202)
J = x θ τ = J T F 2018 School of Information Technology and Electrical Engineering at the University of Queensland Lecture Schedule Week Date Lecture (M: 2:05p-3:50, 50-N202) 1 23-Jul Introduction + Representing
More informationYou w i ll f ol l ow these st eps : Before opening files, the S c e n e panel is active.
You w i ll f ol l ow these st eps : A. O pen a n i m a g e s t a c k. B. Tr a c e t h e d e n d r i t e w i t h t h e user-guided m ode. C. D e t e c t t h e s p i n e s a u t o m a t i c a l l y. D. C
More information17-Nov-2015 PHYS MAXWELL WHEEL. To test the conservation of energy in a system with gravitational, translational and rotational energies.
Objective MAXWELL WHEEL To test the conservation of energy in a system with gravitational, translational and rotational energies. Introduction A wheel is suspended by two cords wrapped on its axis. After
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 informationSome tools and methods for determination of dynamics of hydraulic systems
Some tools and methods for determination of dynamics of hydraulic systems A warm welcome to the course in Hydraulic servo-techniques! The purpose of the exercises given in this material is to make you
More informationA Simple Introduction to EES Version (Handout version 5.1) Copyright C. S. Tritt, Ph.D. September 20, 2005
A Simple Version 7.441 (Handout version 5.1) Copyright C. S. Tritt, Ph.D. September 20, 2005 The BE-381 textbook, Cengel & Turner, 2 ed., comes with a limited version of the EES software package. The academic
More informationSolution for Assignment 4
SE3XO3/CS4X3 11F4.1 Solution for Assignment 4 Due. Dec. 1, Thursday, 13:3. 1. (4 marks) The integral defining the error function erf(x) = 2 x e t2 dt π is fairly easy to evaluate numerically. Write a MATLAB/Octave
More informationComputer Aided Control Design
Computer Aided Control Design Project-Lab 3 Automatic Control Basic Course, EL1000/EL1100/EL1120 Revised August 18, 2008 Modified version of laboration developed by Håkan Fortell and Svante Gunnarsson
More information