259 Lecture 15 Introduction to MATLAB What is

259 Lecture 15 Introduction to MATLAB

What is MATLAB? o o o MATLAB, which stands for “MATrix LABoratory” is a highperformance language for technical computing. It integrates computation, visualization, and programming in an easy-to-use environment where problems and solutions are expressed in familiar mathematical notation. Typical uses include: n n n n Math and computation Algorithm development Data acquisition Modeling, simulation, and prototyping Data analysis, exploration, and visualization Scientific and engineering graphics Application development, including graphical user interface building 2

What is MATLAB? o MATLAB is an interactive system whose basic data element is an array that does not require dimensioning. o In many university environments, it is the standard instructional tool for introductory and advanced courses in mathematics, engineering, and science. o In industry, MATLAB is the tool of choice for high-productivity research, development, and analysis. 3

What is MATLAB? o MATLAB features a family of add-on application-specific solutions called toolboxes, which allow you to learn and apply specialized technology. n Toolboxes are comprehensive collections of MATLAB functions (M-files) that extend the MATLAB environment to solve particular classes of problems. n Areas in which toolboxes are available include signal processing, control systems, neural networks, fuzzy logic, wavelets, simulation, and many others. 4

The MATLAB System o The MATLAB system consists of these main parts: n Desktop Tools and Development Environment n The MATLAB Mathematical Function Library n The MATLAB Language (syntax) n Graphics n The MATLAB External Interfaces/Application Programming Interface (used to design software) 5

Octave – A Free MATLAB Clone o One of the “drawbacks” to MATLAB is that it is very expensive (~$500 for an academic license). o One possible alternative to MATLAB is OCTAVE, which is a free clone of MATLAB (i. e. Octave mimics the functionality of MATLAB, often in the exact same way). o For more information on Octave (technically GNU Octave), see: n Octave downloads from Source Forge (both Windows and Mac OS X installers can be found here: http: //octave. sourceforge. net/ n Octave Homepage: http: //www. gnu. org/software/octave/ 6

MATLAB Desktop Current Directory Window Enter Commands here in the Command Window Command History Window 7

Octave (CLI) – “Old Interface” 8

Octave (GUI) – “New Interface” 9

Calling MATLAB Help o To call Help in MATLAB, either click on the MATLAB Help link, or type “helpdesk” in the Command Window. n Octave does NOT have a built-in help file, so “helpdesk” doesn’t apply! o Notice that the Command History window keeps track of all of the commands you type! 10

Ending a MATLAB Session o To end a MATLAB session, type “exit” or “quit”. 11

MATLAB Basics o If you forget how to use a command, but can remember the command name, such as “plot”, type “help plot”. o A description of the command, along with syntax, options, related commands, and examples will appear! 12

MATLAB Basics o In MATLAB, another way to get information about a command is via “doc”, which gives more detailed information than “help”. o Try “doc plot”. 13

MATLAB Basics o Many MATLAB functions are created from M-files, which are programs that contain a set of commands. o These commands are run in order to implement the function. o To see the contents of an m-file, use the “type” command. o This command does not work as well in Octave! o Try “type linspace”. 14

MATLAB Basics o To re-enter a command, use the up or down arrow, or type the first few letters of the command, and hit return! o To clear out your Command Window, use “clc”. o Commands in the Command History can be re-entered by double-clicking on them! n The Octave GUI must be on for this to work! 15

Entering Matrices o o o Matrices are the basic “building blocks” of MATLAB - essentially everything MATLAB does is based on calculations involving a matrix. The most basic way to enter a matrix is by typing in the entries, inside square brackets, with row entries separated by spaces and new rows designated by semicolons. Try entering these! 16

The Colon Operator o Another way to enter a matrix is via the colon, which is used to define increments within a vector. o Try each of the following: n x = [1: 5] n y = [1: 0. 1: 5] n z = [0: pi/4: 2*pi] 17

The Colon Operator o We can “pull out” parts of a matrix with the colon! o Try each of the following: o A(2, : ) is row 2 of matrix A. o A(: , 3) is column 3 of matrix A. o B(1: 2: 3, : ) is rows 1 and 3 of matrix B. 18

The Colon Operator o A(2, : ) = [] removes the second row from A. o To insert the row back in, we can concatenate matrices via A = [A(1, : ); 0 -1 6; A(2, : )]. o Try w = [x x]. o For more information on the colon operator, type “help colon”. 19

Matrix Operations o o o To re-display matrix A, type “A”. Matrix Operations: Addition: + Subtraction: Matrix Multiplication: * Matrix Power: ^ Transpose – i. e. switch columns and rows (conjugate transpose in the complex case): ’ To transpose a matrix with complex entries, use. ’ Inverse of a square matrix A: inv(A). Matrix Left Division: (usage AB = inv(A)*B) Matrix Right Division: / (usage A/B = A*inv(B)) o o o o o Example 1: Try each of the following! A B A+A A-A A+B A*B B*A inv(A) A*inv(A) A/A AB inv(A)*B A^3 A*A*A A’ A(1, 1) = 2+i A’ A. ’ 20

Matrix Operations - Example 1 (cont. ) 21

Array Operations o o o Scalar Multiplication: k*A with k a scalar. “Normal” arithmetic works as expected with the usual order of operations – each number is a 1 x 1 matrix. Array Multiplication, i. e. entry -wise multiplication of matrices of the same size: A. *B Array Division, i. e. entry-wise division of matrices of the same size: A. /B Array Exponentiation, i. e. entry-wise exponentiation of matrices of the same size: A. ^B o o o Example 2: Try each of the following! A(1, 1)=1 C = 5*A C(: , 1) = [1 1 1] A. *C A. /C A. *B A. ^C A. *inv(A) A. ^3 1+2/3*4^7 22

Array Operations – Example 2 (cont. ) 23

Special Matrices and Matrix Functions “eye(n)” gives an n x n identity matrix “zeros(m, n)” gives an m x n matrix of 0’s. “ones(m, n)” gives an m x n matrix of 1’s. “det(A)” computes the determinant of an n x n matrix A. (Recall that for A to be invertible, det(A) 0 must hold!) o To clear matrix A out of the Workspace, use “clear A”. To clear ALL variables out of the Workspace, use “clear all”. o o 24

Linear Systems of Equations o MATLAB can be used to solve systems of linear equations. o Example 3: Solve the system x + 2 y = 1 3 x + 4 y = -1 o Solution: This system can be written in matrix form AX=b with: 25

Linear Systems of Equations o Example 3 (cont. ) o Check if det(A) 0. o If so, then X = inv(A)*b will follow! 26

Making a Table of Values for a Function o o o o MATLAB can be used to evaluate functions at input values entered as a vector! For multiplication, division, and powers, use the array versions of these operators. Example 4: Try the following commands to make a table of values of the function y = (x*sin(x 2)-4)/(x+1). x = [1: 0. 5: 5]; % The ; suppresses output to the screen! y = (x. *sin(x. ^2)-4). /(x+1); % Array operators used here. [x; y] % Concatenate the two vectors. [x; y]’ % Transpose the matrix with x in column 1 and y in column 2 to get our table! 27

Creating an M-file o o o o One of the most powerful features of MATLAB is the Mfile, which is a file containing a set of commands that can be executed by MATLAB. M-files syntax: *. m These can be created as a text file and saved as an M-file. To call the editor to make an Mfile, type “edit”. M-files can also be made from within MATLAB via pull-down menus File->New->Script or by highlighting the commands in the Command History window. n In the Octave GUI, choose File->New Script. Create an M-file that contains the commands used to make the table in Example 4 above! Save it on the Desktop as example 4. m 28

MATLAB’s Path and Running an Mfile o o To run the M-file we created, we need to put the directory containing the file in MATLAB’s path. To do this, type “addpath ‘DIRECTORY’ ”, where the ‘DIRECTORY’ is where the file is located, for example, a command like this should be used: “addpath ‘C: UsersMAKDesktop’” In MATLAB, this can also be done via the pull-down menus File->Set Path. To run the set of commands, type the name of the M-file in the Command Window! 29

The Plot Command o MATLAB’s “plot” command is used to draw graphs of functions. o This is done by plotting (x, y) ordered pairs on a coordinate plane. o To plot the table created above, use the command “plot(x, y)”. o Notice that the graph we get is very jagged – to smooth it out, add more x-values and corresponding y-values to the vectors x and y! 30

The Plot Command o o o To the right is the same graph as above, with the commands: x = [1: 0. 05: 5]; y = (x. *sin(x. ^2)-4). /(x+1); plot(x, y) To add a title, axes labels, and a grid to the graph, use the “title”, “xlabel”, “ylabel”, and “grid” commands! title(‘y = f(x)’) xlabel(‘x’) ylabel(‘y’) grid on Try these, then put them into your example 4 M-file. Try plotting y = x*cos(x) on the interval [- , 3 ]. 31

Plotting More Than One Function o o o To plot more than one function on the same graph, use the plot command with the sets of input and output variables listed in order: plot(x, y, x 1, y 1). Colors and plot styles can also be specified – to see a complete list of the available options, type “help plot” or “doc plot”. Here is how to plot the functions f(x) = x^2 and f’(x) = 2 x on the same graph! x = [-1: 0. 1: 2]; y = x. *x; y 1 = 2*x; plot(x, y, ’ro’, x, y 1, ’b-’) title(‘y = f(x) and y = f’’(x)’) xlabel(‘x’) ylabel(‘y’) legend(’y=x^2’, ’y’’=2 x’) 32

Piecewise Defined Functions o o o In order to plot a piecewise defined function, we concatenate the separate pieces. x 1 = [-2: 0. 1: 0]; y 1 = x 1. ^2 + 1; x 2 = [0: 0. 1: 2]; y 2 = sin(x 2) +1; x = [x 1 x 2]; y = [y 1 y 2]; plot(x, y, ’g’) xlabel(‘x’) ylabel(‘y’) title(’f(x)=x^2, if x is in [-2, 0]; sin(x) + 1, if x is in (0, 2]’) 33

References o Using MATLAB in Calculus by Gary Jenson o MATLAB Help Files – Getting Started o MATLAB Tutorial from the University of Utah at http: //www. math. utah. edu/lab/ms/matl ab/matlab. html o Octave n http: //octave. sourceforge. net/ n http: //www. gnu. org/software/octave/ 34