Analytical Mechanics: MATLAB
|
|
- Lesley Flynn
- 5 years ago
- Views:
Transcription
1 Analytical Mechanics: MATLAB Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 1 / 47
2 Agenda Agenda 1 Vector and Matrix 2 Graph 3 Ordinary Differential Equations 4 Optimization 5 Parameter Passing 6 Random Numbers 7 Summary Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 2 / 47
3 Agenda Problem We drive a 2-DOF open loop menipulator based on joint PID control. Let us simulate the motion of the manipulator. l2 y lc2 link 2 joint 1 lc1 l1 θ1 θ2 joint 2 link 1 x Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 3 / 47
4 Agenda Problem Step 1. Step 2. Step 3. Step 4. Derive equations of motion (kinematics / dynamics) Numerically solve the derived equations of motion Describe the derived numerical solution by graphs or movies (visualization) Analyze the simulated motion Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 4 / 47
5 Agenda Problem Soloving a set of simultaneous linear equations [ ] [ ] H11 H 12 ω1 = H 12 H 22 ω 2 [ f1 f 2 ] Solving a set of ordinary differential equations [ ] [ ] ω1 α1 (θ = 1, θ 2, ω 1, ω 2 ) ω 2 α 2 (θ 1, θ 2, ω 1, ω 2 ) Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 5 / 47
6 Agenda What is MATLAB? 1 Software for numerical calculation 2 can handle vectors or matrices directly 3 Functions such as ODE solvers and optimization 4 Toolboxes for various applications 5 both programming and interactive calculation Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 6 / 47
7 Agenda What is MATLAB? MATLAB environment MATLAB Total Academic Headcount (TAH) MATLAB with all toolboxes is available April March Information matlab.html rainbow/service/software_matlab_student.html Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 7 / 47
8 Agenda What is MATLAB? Install MATLAB into your own PC or mobile Sample programs are on the web of the class Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 8 / 47
9 Vector and Matrix Vector and Matrix Column vector x = [ 2; 3; -1 ]; Row vector y = [ 2, 3, -1 ]; Matrix A = [ 4, -2, 1;... -2, 5, 2;... -2, 3, 2 ]; Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 9 / 47
10 Vector and Matrix Vector and Matrix Symbol... implies that the sentense continues. Column vector x = [ 2;... 3; ]; Column vector x = [ 2; 3; -1 ]; Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 10 / 47
11 Vector and Matrix Vector and Matrix Multiplication p = A*x; q = y*a; >> p p = >> Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 11 / 47
12 Vector and Matrix Vector and Matrix Multiplication p = A*x; q = y*a; >> q q = >> Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 12 / 47
13 Vector and Matrix Matrix operations >> A A = >> A(3,2) ans = 3 Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 13 / 47
14 Vector and Matrix Matrix operations >> A A = >> A(3,2) = 6; >> A A = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 14 / 47
15 Vector and Matrix Matrix operations >> A(3,:) ans = >> A(:,2) ans = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 15 / 47
16 Vector and Matrix Matrix operations >> A A = >> A(:,2) = [ 0; 2; 1 ]; >> A A = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 16 / 47
17 Vector and Matrix Matrix operations >> A A = >> A(3,:) = [ 3, -5, -1 ]; >> A A = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 17 / 47
18 Vector and Matrix Matrix operations >> A A = >> B = A([1,3],:); >> B B = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 18 / 47
19 Vector and Matrix Matrix operations >> A A = >> C = A(:,[2,1]); >> C C = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 19 / 47
20 Vector and Matrix Basic row operations A(3,:) = 5*A(3,:); multiply the 3rd row by 5 A(1,:) = A(1,:) + 4*A(2,:); add 4-times of the 2nd row to the 1st row A([3,1],:) = A([1,3],:); exchange the 1st and the 3rd rows Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 20 / 47
21 Vector and Matrix Solving simultaneous linear equation A = [ 4, -2, 1;... -2, 5, 2;... -2, 3, 2 ]; p = [ 1; 9; 3 ]; Solve a simultaneous linear equation Ax = p >> x = A\p; >> x x = >> A*x Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 21 / 47
22 Vector and Matrix Solving simultaneous linear equation operator \ is general but less effective when coefficient matrix is positive-definite and symmetric, apply Cholesky decomposition inertia matrices are positive-definite and symmetric Cholesky decomposition positive-definite and symmetric matrix M can be decomposed as M = U T U where U is an upper trianglular matrix. Mx = b = U T Ux = b = { U T y = b Ux = y Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 22 / 47
23 Vector and Matrix Cholesky decomposition program Cholesky.m fprintf( Cholesky decomposition\n ); M = [ 4, -2, -2;... -2, 2, 0;... -2, 0, 3 ]; U = chol(m); U U *U Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 23 / 47
24 Vector and Matrix Cholesky decomposition >> Cholesky Cholesky decomposition U = ans = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 24 / 47
25 Vector and Matrix Cholesky decomposition program b = [ 4; 2; -7 ]; y = U \b; x = U\y; x result x = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 25 / 47
26 Graph Graph >> x = [0:10] x = >> f = x.*x f = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 26 / 47
27 Graph Graph >> plot(x,f) Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 27 / 47
28 Graph Graph >> t = [0:0.1:10] t = >> x = sin(t) x = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 28 / 47
29 Graph Graph >> plot(t,x) Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 29 / 47
30 Ordinary Differential Equations Solving Ordinary Differential Equations van der Pol equation ẍ 2(1 x 2 )ẋ + x = 0 q = [ x v { ẋ = v v = 2(1 x 2 )v x ] [, q = f (t, q) = v 2(1 x 2 )v x ] Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 30 / 47
31 Ordinary Differential Equations Solving Ordinary Differential Equations File van_der_pol.m describes function f (t, q) function dotq = van_der_pol (t, q) x = q(1); v = q(2); dotx = v; dotv = 2*(1-x^2)*v - x; dotq = [dotx; dotv]; end File name van der Pol should be consistent to function name van der Pol. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 31 / 47
32 Ordinary Differential Equations Solving Ordinary Differential Equations Program solve_van_der_pol.m interval = 0.00:0.10:10.00; qinit = [ 2.00; 0.00 ]; [time, q] = ode45(@van_der_pol, interval, qinit); Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 32 / 47
33 Ordinary Differential Equations Solving Ordinary Differential Equations Draw a graph of time t and variable x plot(time, q(:,1), - ); Draw a graph of time t and variable v plot(time, q(:,2), - ); - solid line -- broken line -. dot-dash line : dotted line Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 33 / 47
34 Ordinary Differential Equations Solving Ordinary Differential Equations graph of time t and variable x Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 34 / 47
35 Ordinary Differential Equations Solving Ordinary Differential Equations graph of time t and variable v Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 35 / 47
36 Optimization Optimization Minimizing Rosenbrock function File Rosenbrock.m minimize f (x 1, x 2 ) = 100(x 2 x 2 1 ) 2 + (1 x 1 ) 2 function f = Rosenbrock( x ) x1 = x(1); x2 = x(2); f = 100*(x2 - x1^2)^2 + (1 - x1)^2; end Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 36 / 47
37 Optimization Optimization File Rosenbrock_minimize.m Result xinit = [ -1.2; 1.0 ]; [xmin, fmin] = fminsearch(@rosenbrock, xinit); xmin fmin >> Rosenbrock_minimize xmin = fmin = e-10 Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 37 / 47
38 Parameter Passing ODE with Parameter ordinary differential equation ẍ + bẋ + 9x = 0 where b is a parameter ẋ = v v = bv 9x Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 38 / 47
39 Parameter Passing Global Variable Function Program function dotq = damped_vibration (t, q) global b; x = q(1); v = q(2); dotx = v; dotv = -b*v - 9*x; dotq = [dotx; dotv]; end global b; interval = [0,10]; qinit = [2.00;0.00]; b = 1.00; [time,q] = ode45(@damped_vibration,interval,qinit); Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 39 / 47
40 Parameter Passing Nested Function Function with arguments of time, state variable vector, and parameter Program function dotq = damped_vibration_param (t, q, b) x = q(1); v = q(2); dotx = v; dotv = -b*v - 9*x; dotq = [dotx; dotv]; end interval = [0,10]; qinit = [2.00;0.00]; b = 1.00; damped_vibration damped_vibration_param (t,q,b); [time,q] = ode45(damped_vibration,interval,qinit); Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 40 / 47
41 Parameter Passing Global Variable vs Nested Function Global Variable Simple program Global variables may conflict against local variables Nested Function Somewhat complicated Must perform function defininion whenever paramerer values change Never conflict with other variables Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 41 / 47
42 Random Numbers Uniform Random Numbers Uniform Random Numbers in interval ( 0, 1 ) rng( shuffle, twister ); for k=1:10 x = rand; s = num2str(x); disp(s); end Symbol shuffle generates different random numbers whenever the program runs. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 42 / 47
43 Random Numbers Uniform Random Numbers Uniform Random Numbers in interval ( 0, 1 ) rng(0, twister ); for k=1:10 x = rand; s = num2str(x); disp(s); end specifying seed 0 generates unique random numbers whenever the program runs. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 43 / 47
44 Random Numbers dice.m function k = dice() % simulating a dice x = rand; if x < 1/6.00 k = 1; elseif x < 2/6.00 k = 2; elseif x < 3/6.00 k = 3; elseif x < 4/6.00 k = 4; elseif x < 5/6.00 k = 5; else k = 6; end end Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 44 / 47
45 Random Numbers dice run.m for i=1:10 s = []; for j=1:10 k = dice(); s = [s,, num2str(k)]; end disp(s); end Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 45 / 47
46 Random Numbers >> dice_run >> dice_run Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 46 / 47 dice run.m
47 Summary Summary Numerical calculation using MATLAB linear calculation (vectors and matrices) solving simultaneous linear equations solving ordinary differential equations numerically optimization parameter passing random numbers Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: MATLAB 47 / 47
Analytical Mechanics: Variational Principles
Analytical Mechanics: Variational Principles Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Variational Principles 1 / 71 Agenda
More informationKinetic energy. Analytical Mechanics: Link Mechanisms. Kinetic energy. Agenda. Kinematics of two link open mechanism.
Analtical Mechanics: Link Mechanisms Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Kinetic energ velocit of the center of mass of link : angular velocit of link : kinetic energ of link : ẋ c l c θ S
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationJanuary 18, 2008 Steve Gu. Reference: Eta Kappa Nu, UCLA Iota Gamma Chapter, Introduction to MATLAB,
Introduction to MATLAB January 18, 2008 Steve Gu Reference: Eta Kappa Nu, UCLA Iota Gamma Chapter, Introduction to MATLAB, Part I: Basics MATLAB Environment Getting Help Variables Vectors, Matrices, and
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 informationNumerical Methods Lecture 2 Simultaneous Equations
CGN 42 - Computer Methods Numerical Methods Lecture 2 Simultaneous Equations Topics: matrix operations solving systems of equations Matrix operations: Adding / subtracting Transpose Multiplication Adding
More informationET3-7: Modelling I(V) Introduction and Objectives. Electrical, Mechanical and Thermal Systems
ET3-7: Modelling I(V) Introduction and Objectives Electrical, Mechanical and Thermal Systems Objectives analyse and model basic linear dynamic systems -Electrical -Mechanical -Thermal Recognise the analogies
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 537 Homewors Friedland Text Updated: Wednesday November 8 Some homewor assignments refer to Friedland s text For full credit show all wor. Some problems require hand calculations. In those cases do
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2013 PROBLEM SET 2
STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2013 PROBLEM SET 2 1. You are not allowed to use the svd for this problem, i.e. no arguments should depend on the svd of A or A. Let W be a subspace of C n. The
More informationLecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material.
Lecture 1 1 1 Department of Mathematics Dartmouth College 3/28/07 Outline Course Overview http://www.math.dartmouth.edu/~m23s07 Matlab Ordinary differential equations Definition An ordinary differential
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
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 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 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 informationSection Matrices and Systems of Linear Eqns.
QUIZ: strings Section 14.3 Matrices and Systems of Linear Eqns. Remembering matrices from Ch.2 How to test if 2 matrices are equal Assume equal until proved wrong! else? myflag = logical(1) How to test
More informationMATH 3511 Lecture 1. Solving Linear Systems 1
MATH 3511 Lecture 1 Solving Linear Systems 1 Dmitriy Leykekhman Spring 2012 Goals Review of basic linear algebra Solution of simple linear systems Gaussian elimination D Leykekhman - MATH 3511 Introduction
More informationMath 552 Scientific Computing II Spring SOLUTIONS: Homework Set 1
Math 552 Scientific Computing II Spring 21 SOLUTIONS: Homework Set 1 ( ) a b 1 Let A be the 2 2 matrix A = By hand, use Gaussian elimination with back c d substitution to obtain A 1 by solving the two
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation / 59 Agenda Agenda
More informationScientific Computing
Scientific Computing Direct solution methods Martin van Gijzen Delft University of Technology October 3, 2018 1 Program October 3 Matrix norms LU decomposition Basic algorithm Cost Stability Pivoting Pivoting
More informationGet acquainted with the computer program, The Quadratic Transformer. When you're satisfied that you understand how it works, try the tasks below.
Weaving a Parabola Web with the Quadratic Transformer In this activity, you explore how the graph of a quadratic function and its symbolic expression relate to each other. You start with a set of four
More informationLaw of Trichotomy and Boundary Equations
Law of Trichotomy and Boundary Equations Law of Trichotomy: For any two real numbers a and b, exactly one of the following is true. i. a < b ii. a = b iii. a > b The Law of Trichotomy is a formal statement
More informationET3-7: Modelling II(V) Electrical, Mechanical and Thermal Systems
ET3-7: Modelling II(V) Electrical, Mechanical and Thermal Systems Agenda of the Day 1. Resume of lesson I 2. Basic system models. 3. Models of basic electrical system elements 4. Application of Matlab/Simulink
More information5.1 Banded Storage. u = temperature. The five-point difference operator. uh (x, y + h) 2u h (x, y)+u h (x, y h) uh (x + h, y) 2u h (x, y)+u h (x h, y)
5.1 Banded Storage u = temperature u= u h temperature at gridpoints u h = 1 u= Laplace s equation u= h u = u h = grid size u=1 The five-point difference operator 1 u h =1 uh (x + h, y) 2u h (x, y)+u h
More informationWe wish to solve a system of N simultaneous linear algebraic equations for the N unknowns x 1, x 2,...,x N, that are expressed in the general form
Linear algebra This chapter discusses the solution of sets of linear algebraic equations and defines basic vector/matrix operations The focus is upon elimination methods such as Gaussian elimination, and
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More informationCS 219: Sparse matrix algorithms: Homework 3
CS 219: Sparse matrix algorithms: Homework 3 Assigned April 24, 2013 Due by class time Wednesday, May 1 The Appendix contains definitions and pointers to references for terminology and notation. Problem
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 informationGenerating Functions (Revised Edition)
Math 700 Fall 06 Notes Generating Functions (Revised Edition What is a generating function? An ordinary generating function for a sequence (a n n 0 is the power series A(x = a nx n. The exponential generating
More informationENGIN 211, Engineering Math. Complex Numbers
ENGIN 211, Engineering Math Complex Numbers 1 Imaginary Number and the Symbol J Consider the solutions for this quadratic equation: x 2 + 1 = 0 x = ± 1 1 is called the imaginary number, and we use the
More informationLinear Algebra Section 2.6 : LU Decomposition Section 2.7 : Permutations and transposes Wednesday, February 13th Math 301 Week #4
Linear Algebra Section. : LU Decomposition Section. : Permutations and transposes Wednesday, February 1th Math 01 Week # 1 The LU Decomposition We learned last time that we can factor a invertible matrix
More information1.Chapter Objectives
LU Factorization INDEX 1.Chapter objectives 2.Overview of LU factorization 2.1GAUSS ELIMINATION AS LU FACTORIZATION 2.2LU Factorization with Pivoting 2.3 MATLAB Function: lu 3. CHOLESKY FACTORIZATION 3.1
More informationChapter 3 + some notes on counting the number of degrees of freedom
Chapter 3 + some notes on counting the number of degrees of freedom Minimum number of independent parameters = Some number of dependent parameters minus the number of relationships (equations) you can
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationCEE 361 Methods of Structural Analysis Fall Semester, 2000
Lab 7 Lecture Through the course of today s lecture and this handout, the four functions required for the completion of laboratory #7 will be gone over. Ideally, the information contained in this document,
More informationResearch Reports on Mathematical and Computing Sciences
ISSN 1342-284 Research Reports on Mathematical and Computing Sciences Exploiting Sparsity in Linear and Nonlinear Matrix Inequalities via Positive Semidefinite Matrix Completion Sunyoung Kim, Masakazu
More informationDetermining a span. λ + µ + ν = x 2λ + 2µ 10ν = y λ + 3µ 9ν = z.
Determining a span Set V = R 3 and v 1 = (1, 2, 1), v 2 := (1, 2, 3), v 3 := (1 10, 9). We want to determine the span of these vectors. In other words, given (x, y, z) R 3, when is (x, y, z) span(v 1,
More informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
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 informationA primer on matrices
A primer on matrices Stephen Boyd August 4, 2007 These notes describe the notation of matrices, the mechanics of matrix manipulation, and how to use matrices to formulate and solve sets of simultaneous
More informationLesson Rigid Body Dynamics
Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body
More informationLinear Systems of n equations for n unknowns
Linear Systems of n equations for n unknowns In many application problems we want to find n unknowns, and we have n linear equations Example: Find x,x,x such that the following three equations hold: x
More information5. Surface Temperatures
5. Surface Temperatures For this case we are given a thin rectangular sheet of metal, whose temperature we can conveniently measure at any points along its four edges. If the edge temperatures are held
More information6. 3D Kinematics DE2-EA 2.1: M4DE. Dr Connor Myant
DE2-EA 2.1: M4DE Dr Connor Myant 6. 3D Kinematics Comments and corrections to connor.myant@imperial.ac.uk Lecture resources may be found on Blackboard and at http://connormyant.com Contents Three-Dimensional
More informationSystem Simulation using Matlab
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
More informationCOMS 6100 Class Notes
COMS 6100 Class Notes Daniel Solus September 20, 2016 1 General Remarks The Lecture notes submitted by the class have been very good. Integer division seemed to be a common oversight when working the Fortran
More informationMath 3191 Applied Linear Algebra
Math 191 Applied Linear Algebra Lecture 1: Inner Products, Length, Orthogonality Stephen Billups University of Colorado at Denver Math 191Applied Linear Algebra p.1/ Motivation Not all linear systems have
More informationGaussian Elimination -(3.1) b 1. b 2., b. b n
Gaussian Elimination -() Consider solving a given system of n linear equations in n unknowns: (*) a x a x a n x n b where a ij and b i are constants and x i are unknowns Let a n x a n x a nn x n a a a
More information(x + 1)(x 2) = 4. x
dvanced Integration Techniques: Partial Fractions The method of partial fractions can occasionally make it possible to find the integral of a quotient of rational functions. Partial fractions gives us
More information1 Positive definiteness and semidefiniteness
Positive definiteness and semidefiniteness Zdeněk Dvořák May 9, 205 For integers a, b, and c, let D(a, b, c) be the diagonal matrix with + for i =,..., a, D i,i = for i = a +,..., a + b,. 0 for i = a +
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationNumerical Methods for Inverse Kinematics
Numerical Methods for Inverse Kinematics Niels Joubert, UC Berkeley, CS184 2008-11-25 Inverse Kinematics is used to pose models by specifying endpoints of segments rather than individual joint angles.
More informationMechanics of Structures (CE130N) Lab 3
UNIVERSITY OF CALIFORNIA AT BERKELEY CE 130N, Spring 2009 Department of Civil and Environmental Engineering Prof. S. Govindjee and Dr. T. Koyama Structural Engineering, Mechanics and Materials Lab 3 1
More informationINTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS
INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS % Many processes or systems involving the rates of change (derivatives) of variables can be described by one or more equations which describe the behavior
More informationLab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations
Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in
More informationExercise 1b: Differential Kinematics of the ABB IRB 120
Exercise 1b: Differential Kinematics of the ABB IRB 120 Marco Hutter, Michael Blösch, Dario Bellicoso, Samuel Bachmann October 5, 2016 Abstract The aim of this exercise is to calculate the differential
More informationDYNAMICS OF PARALLEL MANIPULATOR
DYNAMICS OF PARALLEL MANIPULATOR The 6nx6n matrices of manipulator mass M and manipulator angular velocity W are introduced below: M = diag M 1, M 2,, M n W = diag (W 1, W 2,, W n ) From this definitions
More informationExample: Inverted pendulum on cart
Chapter 11 Eample: Inverted pendulum on cart The figure to the right shows a rigid body attached by an frictionless pin (revolute) joint to a cart (modeled as a particle). Thecart slides on a horizontal
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 information+ MATRIX VARIABLES AND TWO DIMENSIONAL ARRAYS
+ MATRIX VARIABLES AND TWO DIMENSIONAL ARRAYS Matrices are organized rows and columns of numbers that mathematical operations can be performed on. MATLAB is organized around the rules of matrix operations.
More informationRobotics I. February 6, 2014
Robotics I February 6, 214 Exercise 1 A pan-tilt 1 camera sensor, such as the commercial webcams in Fig. 1, is mounted on the fixed base of a robot manipulator and is used for pointing at a (point-wise)
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 informationCSCI 239 Discrete Structures of Computer Science Lab 6 Vectors and Matrices
CSCI 239 Discrete Structures of Computer Science Lab 6 Vectors and Matrices This lab consists of exercises on real-valued vectors and matrices. Most of the exercises will required pencil and paper. Put
More informationDynamics of Machines and Mechanisms
Dynamics of Machines and Mechanisms Introduction SPACAR Mechanical Automation (WA) Room: HR Z 2.29 Phone: (053) 489 2557 Email: R.G.K.M.Aarts@utwente.nl Info: http://www.wa.ctw.utwente.nl/lectures/113173/
More informationAPPPHYS 217 Tuesday 6 April 2010
APPPHYS 7 Tuesday 6 April Stability and input-output performance: second-order systems Here we present a detailed example to draw connections between today s topics and our prior review of linear algebra
More informationPrincipal Components Analysis (PCA)
Principal Components Analysis (PCA) Principal Components Analysis (PCA) a technique for finding patterns in data of high dimension Outline:. Eigenvectors and eigenvalues. PCA: a) Getting the data b) Centering
More information7x 5 x 2 x + 2. = 7x 5. (x + 1)(x 2). 4 x
Advanced Integration Techniques: Partial Fractions The method of partial fractions can occasionally make it possible to find the integral of a quotient of rational functions. Partial fractions gives us
More informationAM205: Assignment 2. i=1
AM05: Assignment Question 1 [10 points] (a) [4 points] For p 1, the p-norm for a vector x R n is defined as: ( n ) 1/p x p x i p ( ) i=1 This definition is in fact meaningful for p < 1 as well, although
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 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 informationNUMERICAL METHODS. lor CHEMICAL ENGINEERS. Using Excel', VBA, and MATLAB* VICTOR J. LAW. CRC Press. Taylor & Francis Group
NUMERICAL METHODS lor CHEMICAL ENGINEERS Using Excel', VBA, and MATLAB* VICTOR J. LAW CRC Press Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup,
More informationSolving Equations Quick Reference
Solving Equations Quick Reference Integer Rules Addition: If the signs are the same, add the numbers and keep the sign. If the signs are different, subtract the numbers and keep the sign of the number
More informationCS 221 Lecture 9. Tuesday, 1 November 2011
CS 221 Lecture 9 Tuesday, 1 November 2011 Some slides in this lecture are from the publisher s slides for Engineering Computation: An Introduction Using MATLAB and Excel 2009 McGraw-Hill Today s Agenda
More informationDynamics Algorithms for Multibody Systems
International Conference on Multi Body Dynamics 2011 Vijayawada, India. pp. 351 365 Dynamics Algorithms for Multibody Systems S. V. Shah, S. K. Saha and J. K. Dutt Department of Mechanical Engineering,
More informationPhysics 110. Electricity and Magnetism. Professor Dine. Spring, Handout: Vectors and Tensors: Everything You Need to Know
Physics 110. Electricity and Magnetism. Professor Dine Spring, 2008. Handout: Vectors and Tensors: Everything You Need to Know What makes E&M hard, more than anything else, is the problem that the electric
More informationNew Mexico Tech Hyd 510
Vectors vector - has magnitude and direction (e.g. velocity, specific discharge, hydraulic gradient) scalar - has magnitude only (e.g. porosity, specific yield, storage coefficient) unit vector - a unit
More informationNumerical Linear Algebra
Chapter 3 Numerical Linear Algebra We review some techniques used to solve Ax = b where A is an n n matrix, and x and b are n 1 vectors (column vectors). We then review eigenvalues and eigenvectors and
More information6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski. Solutions to Problem Set 1 1. Massachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Solutions to Problem Set 1 1 Problem 1.1T Consider the
More informationPLC Papers Created For:
PLC Papers Created For: Quadratics intervention Deduce quadratic roots algebraically 1 Grade 6 Objective: Deduce roots algebraically. Question 1. Factorise and solve the equation x 2 8x + 15 = 0 Question
More informationSection 8.3 Partial Fraction Decomposition
Section 8.6 Lecture Notes Page 1 of 10 Section 8.3 Partial Fraction Decomposition Partial fraction decomposition involves decomposing a rational function, or reversing the process of combining two or more
More informationNumerical Linear Algebra
Numerical Linear Algebra Direct Methods Philippe B. Laval KSU Fall 2017 Philippe B. Laval (KSU) Linear Systems: Direct Solution Methods Fall 2017 1 / 14 Introduction The solution of linear systems is one
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 informationAutomatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year
Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
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 informationChapter Practice Test Name: Period: Date:
Name: Period: Date: 1. Draw the graph of the following system: 3 x+ 5 y+ 13 = 0 29 x 11 y 7 = 0 3 13 y = x 3x+ 5y+ 13= 0 5 5 29x 11y 7 = 0 29 7 y = x 11 11 Practice Test Page 1 2. Determine the ordered
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 533 Homeworks Spring 07 Updated: Saturday, April 08, 07 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. Some homework assignments refer to the textbooks: Slotine
More informationTMA4125 Matematikk 4N Spring 2017
Norwegian University of Science and Technology Institutt for matematiske fag TMA15 Matematikk N Spring 17 Solutions to exercise set 1 1 We begin by writing the system as the augmented matrix.139.38.3 6.
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 informationMATLAB Project: LU Factorization
Name Purpose: To practice Lay's LU Factorization Algorithm and see how it is related to MATLAB's lu function. Prerequisite: Section 2.5 MATLAB functions used: *, lu; and ludat and gauss from Laydata4 Toolbox
More informationNumerical Linear Algebra
Numerical Linear Algebra Decompositions, numerical aspects Gerard Sleijpen and Martin van Gijzen September 27, 2017 1 Delft University of Technology Program Lecture 2 LU-decomposition Basic algorithm Cost
More informationProgram Lecture 2. Numerical Linear Algebra. Gaussian elimination (2) Gaussian elimination. Decompositions, numerical aspects
Numerical Linear Algebra Decompositions, numerical aspects Program Lecture 2 LU-decomposition Basic algorithm Cost Stability Pivoting Cholesky decomposition Sparse matrices and reorderings Gerard Sleijpen
More informationLinear System Theory
Linear System Theory - MATLAB Exercise Prof. Robert X. Gao Electromechanical Systems Laboratory Department of Mechanical Engineering Outline Review System Modeling Procedure Example Simple Problem Example
More information[POLS 8500] Review of Linear Algebra, Probability and Information Theory
[POLS 8500] Review of Linear Algebra, Probability and Information Theory Professor Jason Anastasopoulos ljanastas@uga.edu January 12, 2017 For today... Basic linear algebra. Basic probability. Programming
More informationIn this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:
Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this
More informationModel Order Reduction via Matlab Parallel Computing Toolbox. Istanbul Technical University
Model Order Reduction via Matlab Parallel Computing Toolbox E. Fatih Yetkin & Hasan Dağ Istanbul Technical University Computational Science & Engineering Department September 21, 2009 E. Fatih Yetkin (Istanbul
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 informationECEn 483 / ME 431 Case Studies. Randal W. Beard Brigham Young University
ECEn 483 / ME 431 Case Studies Randal W. Beard Brigham Young University Updated: December 2, 2014 ii Contents 1 Single Link Robot Arm 1 2 Pendulum on a Cart 9 3 Satellite Attitude Control 17 4 UUV Roll
More informationSolution Set 3, Fall '12
Solution Set 3, 86 Fall '2 Do Problem 5 from 32 [ 3 5 Solution (a) A = Only one elimination step is needed to produce the 2 6 echelon form The pivot is the in row, column, and the entry to eliminate is
More information