Matlab Tutorial
Matlab Tutorial
Audience
This tutorial has been prepared for the beginners to help them understand
basic to advanced functionality of MATLAB. After completing this tutorial
you will find yourself at a moderate level of expertise in using MATLAB
from where you can take yourself to next levels.
Prerequisites
We assume you have a little knowledge of any computer programming and
understand concepts like variables, constants, expression, statements,
etc. If you have done programming in any other high-level programming
language like C, C++ or Java, then it will be very much beneficial and
learning MATLAB will be like a fun for you.
MATLAB - Overview
It has numerous built-in commands and math functions that help you in
mathematical calculations, generating plots, and performing numerical
methods.
Linear Algebra
Algebraic Equations
Non-linear Functions
Statistics
Data Analysis
Integration
Transforms
Curve Fitting
Features of MATLAB
Following are the basic features of MATLAB −
It provides built-in graphics for visualizing data and tools for creating custom
plots.
Uses of MATLAB
MATLAB is widely used as a computational tool in science and engineering
encompassing the fields of physics, chemistry, math and all engineering
streams. It is used in a range of applications including −
Control Systems
Computational Finance
Computational Biology
MATLAB - Environment Setup
After downloading the installer the software can be installed through few
clicks.
Understanding the MATLAB Environment
MATLAB development IDE can be launched from the icon created on the
desktop. The main working window in MATLAB is called the desktop. When
MATLAB is started, the desktop appears in its default layout −
Current Folder − This panel allows you to access the project folders and files.
Command Window − This is the main area where commands can be entered
at the command line. It is indicated by the command prompt (>>).
Workspace − The workspace shows all the variables created and/or imported
from files.
Command History − This panel shows or return commands that are entered
at the command line.
MATLAB - Basic Syntax
Hands on Practice
Type a valid expression, for example,
5 + 5
When you click the Execute button, or type Ctrl+E, MATLAB executes it
immediately and the result returned is −
ans = 10
When you click the Execute button, or type Ctrl+E, MATLAB executes it
immediately and the result returned is −
ans = 9
Another example,
When you click the Execute button, or type Ctrl+E, MATLAB executes it
immediately and the result returned is −
ans = 1
Another example,
When you click the Execute button, or type Ctrl+E, MATLAB executes it
immediately and the result returned is −
ans = Inf
warning: division by zero
Another example,
732 * 20.3
When you click the Execute button, or type Ctrl+E, MATLAB executes it
immediately and the result returned is −
ans = 1.4860e+04
For example,
x = 3;
y = x + 5
When you click the Execute button, or type Ctrl+E, MATLAB executes it
immediately and the result returned is −
y = 8
Adding Comments
The percent symbol (%) is used for indicating a comment line. For
example,
You can also write a block of comments using the block comment operators
% { and % }.
The MATLAB editor includes tools and context menu items to help you add,
remove, or change the format of comments.
\ Left-division operator.
/ Right-division operator.
. Decimal point.
= Assignment operator.
Name Meaning
Inf Infinity.
pi The number π
Naming Variables
Variable names consist of a letter followed by any number of letters, digits
or underscore.
MATLAB is case-sensitive.
Variable names can be of any length, however, MATLAB uses only first N
characters, where N is given by the function namelengthmax.
For example,
save myfile
You can reload the file anytime later using the load command.
load myfile
MATLAB - Variables
MATLAB will execute the above statement and return the following result
−
x = 3
It creates a 1-by-1 matrix named x and stores the value 3 in its element.
Let us check another example,
MATLAB will execute the above statement and return the following result
−
x = 4
Once a variable is entered into the system, you can refer to it later.
When an expression returns a result that is not assigned to any variable, the
system assigns it to a variable named ans, which can be used later.
For example,
sqrt(78)
MATLAB will execute the above statement and return the following result
−
ans = 8.8318
sqrt(78);
9876/ans
MATLAB will execute the above statement and return the following result
−
ans = 1118.2
x = 7 * 8;
y = x * 7.89
MATLAB will execute the above statement and return the following result
−
y = 441.84
Multiple Assignments
You can have multiple assignments on the same line. For example,
a = 2; b = 7; c = a * b
MATLAB will execute the above statement and return the following result
−
c = 14
MATLAB will execute the above statement and return the following result
−
Your variables are:
a ans b c
MATLAB will execute the above statement and return the following result
−
Attr Name Size Bytes Class
==== ==== ==== ==== =====
a 1x1 8 double
ans 1x70 757 cell
b 1x1 8 double
c 1x1 8 double
The clear command deletes all (or the specified) variable(s) from the
memory.
clear x % it will delete x, won't display anything
clear % it will delete all variables in the workspace
% peacefully and unobtrusively
Long Assignments
Long assignments can be extended to another line by using an ellipses
(...). For example,
initial_velocity = 0;
acceleration = 9.8;
time = 20;
MATLAB will execute the above statement and return the following result
−
final_velocity = 196
For example −
format long
x = 7 + 10/3 + 5 ^ 1.2
MATLAB will execute the above statement and return the following result−
x = 17.2319816406394
Another example,
format short
x = 7 + 10/3 + 5 ^ 1.2
MATLAB will execute the above statement and return the following result
−
x = 17.232
The format bank command rounds numbers to two decimal places. For
example,
format bank
daily_wage = 177.45;
weekly_wage = daily_wage * 6
MATLAB will execute the above statement and return the following result
−
weekly_wage = 1064.70
For example,
format short e
4.678 * 4.9
MATLAB will execute the above statement and return the following result
−
ans = 2.2922e+01
format long e
x = pi
MATLAB will execute the above statement and return the following result
−
x = 3.141592653589793e+00
The format rat command gives the closest rational expression resulting
from a calculation. For example,
format rat
4.678 * 4.9
MATLAB will execute the above statement and return the following result
−
ans = 34177/1491
Creating Vectors
Row vectors
Column vectors
For example,
r = [7 8 9 10 11]
MATLAB will execute the above statement and return the following result
−
r =
7 8 9 10 11
Another example,
r = [7 8 9 10 11];
t = [2, 3, 4, 5, 6];
res = r + t
MATLAB will execute the above statement and return the following result
−
res =
9 11 13 15 17
MATLAB will execute the above statement and return the following result
−
c =
7
8
9
10
11
Creating Matrices
A matrix is a two-dimensional array of numbers.
m = [1 2 3; 4 5 6; 7 8 9]
MATLAB will execute the above statement and return the following result
−
m =
1 2 3
4 5 6
7 8 9
MATLAB - Commands
Command Purpose
Command Purpose
Command Purpose
The fscanf and fprintf commands behave like C scanf and printf
functions. They support the following format codes −
%s Format as a string.
%d Format as an integer.
The format function has the following forms used for numeric display −
Command Purpose
Plotting Commands
MATLAB provides numerous commands for plotting graphs. The following
table shows some of the commonly used commands for plotting −
Command Purpose
In previous chapters, you have learned how to enter commands from the
MATLAB command prompt. MATLAB also allows you to write series of
commands into a file and execute the file as complete unit, like writing a
function and calling it.
The M Files
MATLAB allows writing two kinds of program files −
Scripts − script files are program files with .m extension. In these files, you
write series of commands, which you want to execute together. Scripts do not
accept inputs and do not return any outputs. They operate on data in the
workspace.
You can use the MATLAB editor or any other text editor to create
your .mfiles. In this section, we will discuss the script files. A script file
contains multiple sequential lines of MATLAB commands and function calls.
You can run a script by typing its name at the command line.
If you are using the command prompt, type edit in the command prompt.
This will open the editor. You can directly type edit and then the filename
(with .m extension)
edit
Or
edit <filename>
The above command will create the file in default MATLAB directory. If you
want to store all program files in a specific folder, then you will have to
provide the entire path.
Let us create a folder named progs. Type the following commands at the
command prompt (>>) −
If you are creating the file for first time, MATLAB prompts you to confirm
it. Click Yes.
Alternatively, if you are using the IDE, choose NEW -> Script. This also
opens the editor and creates a file named Untitled. You can name and save
the file after typing the code.
NoOfStudents = 6000;
TeachingStaff = 150;
NonTeachingStaff = 20;
Total = NoOfStudents + TeachingStaff ...
+ NonTeachingStaff;
disp(Total);
After creating and saving the file, you can run it in two ways −
Just typing the filename (without extension) in the command prompt: >>
prog1
Example
Create a script file, and type the following code −
a = 5; b = 7;
c = a + b
d = c + sin(b)
e = 5 * d
f = exp(-d)
When the above code is compiled and executed, it produces the following
result −
c = 12
d = 12.657
e = 63.285
f = 3.1852e-06
MATLAB - Data Types
If the variable already exists, then MATLAB replaces the original content
with new content and allocates new storage space, where necessary.
For example,
Total = 42
The above statement creates a 1-by-1 matrix named 'Total' and stores the
value 42 in it.
The following table shows the most commonly used data types in MATLAB
−
1 int8
2 uint8
3 int16
4 uint16
6 uint32
7 int64
8 uint64
9 single
10 double
11 logical
12 char
13 cell array
14 structure
pointer to a function
16 user classes
17 java classes
Example
Create a script file with the following code −
n = 2345
d = double(n)
un = uint32(789.50)
rn = 5678.92347
c = int32(rn)
When the above code is compiled and executed, it produces the following
result −
str = Hello World!
n = 2345
d = 2345
un = 790
rn = 5678.9
c = 5679
Function Purpose
Following table provides the functions for determining the data type of a
variable −
Function Purpose
is Detect state
Example
Create a script file with the following code −
x = 3
isinteger(x)
isfloat(x)
isvector(x)
isscalar(x)
isnumeric(x)
x = 23.54
isinteger(x)
isfloat(x)
isvector(x)
isscalar(x)
isnumeric(x)
x = [1 2 3]
isinteger(x)
isfloat(x)
isvector(x)
isscalar(x)
x = 'Hello'
isinteger(x)
isfloat(x)
isvector(x)
isscalar(x)
isnumeric(x)
ans = 0
ans = 1
ans = 1
ans = 0
x = Hello
ans = 0
ans = 0
ans = 1
ans = 0
ans = 0
MATLAB - Operators
Arithmetic Operators
Relational Operators
Logical Operators
Bitwise Operations
Set Operations
Arithmetic Operators
MATLAB allows two different types of arithmetic operations −
The matrix operators and array operators are differentiated by the period
(.) symbol. However, as the addition and subtraction operation is same for
matrices and arrays, the operator is same for both cases. The following
table gives brief description of the operators −
Show Examples
1 +
Addition or unary plus. A+B adds the values stored in variables A and B.
A and B must have the same size, unless one is a scalar. A scalar can be
added to a matrix of any size.
2 -
3 *
4 .*
5 /
Slash or matrix right division. B/A is roughly the same as B*inv(A). More
precisely, B/A = (A'\B')'.
6 ./
Array right division. A./B is the matrix with elements A(i,j)/B(i,j). A and B
must have the same size, unless one of them is a scalar.
7 \
8 .\
Array left division. A.\B is the matrix with elements B(i,j)/A(i,j). A and B
must have the same size, unless one of them is a scalar.
9 ^
10 .^
Array power. A.^B is the matrix with elements A(i,j) to the B(i,j) power.
A and B must have the same size, unless one of them is a scalar.
11 '
12 .'
Relational Operators
Relational operators can also work on both scalar and non-scalar data.
Relational operators for arrays perform element-by-element comparisons
between two arrays and return a logical array of the same size, with
elements set to logical 1 (true) where the relation is true and elements set
to logical 0 (false) where it is not.
1 <
Less than
2 <=
3 >
Greater than
4 >=
5 ==
Equal to
6 ~=
Not equal to
Logical Operators
MATLAB offers two types of logical operators and functions −
Bitwise Operations
Bitwise operators work on bits and perform bit-by-bit operation. The truth
tables for &, |, and ^ are as follows −
0 1 0 1 1
1 1 1 1 0
1 0 0 1 1
Assume if A = 60; and B = 13; Now in binary format they will be as follows
−
A = 0011 1100
B = 0000 1101
-----------------
~A = 1100 0011
Function Purpose
Set Operations
MATLAB provides various functions for set operations, like union,
intersection and testing for set membership, etc.
Show Examples
1 intersect(A,B)
Set intersection of two arrays; returns the values common to both A and
B. The values returned are in sorted order.
2 intersect(A,B,'rows')
Treats each row of A and each row of B as single entities and returns the
rows common to both A and B. The rows of the returned matrix are in
sorted order.
3 ismember(A,B)
4 ismember(A,B,'rows')
Treats each row of A and each row of B as single entities and returns a
vector containing 1 (true) where the rows of matrix A are also rows of B.
Elsewhere, it returns 0 (false).
5 issorted(A)
Returns logical 1 (true) if the elements of A are in sorted order and logical
0 (false) otherwise. Input A can be a vector or an N-by-1 or 1-by-N cell
array of strings. A is considered to be sorted if Aand the output of
sort(A) are equal.
6 issorted(A, 'rows')
7 setdiff(A,B)
Sets difference of two arrays; returns the values in A that are not in B.
The values in the returned array are in sorted order.
8 setdiff(A,B,'rows')
Treats each row of A and each row of B as single entities and returns the
rows from A that are not in B. The rows of the returned matrix are in
sorted order.
9 setxor
10 union
11 unique
4 nested if statements
5 switch statement
There may be a situation when you need to execute a block of code several
number of times. In general, statements are executed sequentially. The
first statement in a function is executed first, followed by the second, and
so on.
1 while loop
3 nested loops
You can use one or more loops inside any another loop.
MATLAB supports the following control statements. Click the following links
to check their detail.
1 break statement
2 continue statement
Causes the loop to skip the remainder of its body and immediately retest
its condition prior to reiterating.
MATLAB - Vectors
Row vectors
Column vectors
Row Vectors
Row vectors are created by enclosing the set of elements in square
brackets, using space or comma to delimit the elements.
r = [7 8 9 10 11]
MATLAB will execute the above statement and return the following result
−
r =
7 8 9 10 11
Column Vectors
Column vectors are created by enclosing the set of elements in square
brackets, using semicolon to delimit the elements.
MATLAB will execute the above statement and return the following result
−
c =
7
8
9
10
11
v(3)
MATLAB will execute the above statement and return the following result
−
ans = 3
When you reference a vector with a colon, such as v(:), all the components
of the vector are listed.
v(:)
MATLAB will execute the above statement and return the following result
−
ans =
1
2
3
4
5
6
rv = [1 2 3 4 5 6 7 8 9];
sub_rv = rv(3:7)
MATLAB will execute the above statement and return the following result
−
sub_rv =
3 4 5 6 7
Vector Operations
In this section, let us discuss the following vector operations −
Transpose of a Vector
Appending Vectors
Magnitude of a Vector
Vector Dot Product
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8]
MATLAB will execute the above statement and return the following result
−
a =
1 2 3 4 5
2 3 4 5 6
3 4 5 6 7
4 5 6 7 8
For example, to refer to the element in the 2nd row and 5th column, of the
matrix a, as created in the last section, we type −
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8];
a(2,5)
MATLAB will execute the above statement and return the following result
−
ans = 6
Let us create a column vector v, from the elements of the 4 th row of the
matrix a −
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8];
v = a(:,4)
MATLAB will execute the above statement and return the following result
−
v =
4
5
6
7
You can also select the elements in the mth through nth columns, for this
we write −
a(:,m:n)
Let us create a smaller matrix taking the elements from the second and
third columns −
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8];
a(:, 2:3)
MATLAB will execute the above statement and return the following result
−
ans =
2 3
3 4
4 5
5 6
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8];
a(:, 2:3)
MATLAB will execute the above statement and return the following result
−
ans =
2 3
3 4
4 5
5 6
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8];
sa = a(2:3,2:4)
MATLAB will execute the above statement and return the following result
−
sa =
3 4 5
4 5 6
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8];
a( 4 , : ) = []
MATLAB will execute the above statement and return the following result
−
a =
1 2 3 4 5
2 3 4 5 6
3 4 5 6 7
a = [ 1 2 3 4 5; 2 3 4 5 6; 3 4 5 6 7; 4 5 6 7 8];
a(: , 5)=[]
MATLAB will execute the above statement and return the following result
−
a =
1 2 3 4
2 3 4 5
3 4 5 6
4 5 6 7
Example
In this example, let us create a 3-by-3 matrix m, then we will copy the
second and third rows of this matrix twice to create a 4-by-3 matrix.
Create a script file with the following code −
a = [ 1 2 3 ; 4 5 6; 7 8 9];
new_mat = a([2,3,2,3],:)
Matrix Operations
In this section, let us discuss the following basic and commonly used
matrix operations −
Division of Matrices
Transpose of a Matrix
Concatenating Matrices
Matrix Multiplication
Determinant of a Matrix
Inverse of a Matrix
MATLAB - Arrays
For example −
zeros(5)
MATLAB will execute the above statement and return the following result
−
ans =
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
For example −
ones(4,3)
MATLAB will execute the above statement and return the following result
−
ans =
1 1 1
1 1 1
1 1 1
1 1 1
For example −
eye(4)
MATLAB will execute the above statement and return the following result
−
ans =
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
For example −
rand(3, 5)
MATLAB will execute the above statement and return the following result
−
ans =
0.8147 0.9134 0.2785 0.9649 0.9572
0.9058 0.6324 0.5469 0.1576 0.4854
0.1270 0.0975 0.9575 0.9706 0.8003
A Magic Square
A magic square is a square that produces the same sum, when its
elements are added row-wise, column-wise or diagonally.
magic(4)
MATLAB will execute the above statement and return the following result
−
ans =
16 2 3 13
5 11 10 8
9 7 6 12
4 14 15 1
Multidimensional Arrays
An array having more than two dimensions is called a multidimensional
array in MATLAB. Multidimensional arrays in MATLAB are an extension of
the normal two-dimensional matrix.
Generally to generate a multidimensional array, we first create a two-
dimensional array and extend it.
a = [7 9 5; 6 1 9; 4 3 2]
MATLAB will execute the above statement and return the following result
−
a =
7 9 5
6 1 9
4 3 2
a(:, :, 2)= [ 1 2 3; 4 5 6; 7 8 9]
MATLAB will execute the above statement and return the following result
−
a =
ans(:,:,1) =
0 0 0
0 0 0
0 0 0
ans(:,:,2) =
1 2 3
4 5 6
7 8 9
For example,
b = rand(4,3,2)
MATLAB will execute the above statement and return the following result
−
b(:,:,1) =
0.0344 0.7952 0.6463
0.4387 0.1869 0.7094
0.3816 0.4898 0.7547
0.7655 0.4456 0.2760
b(:,:,2) =
0.6797 0.4984 0.2238
0.6551 0.9597 0.7513
0.1626 0.3404 0.2551
0.1190 0.5853 0.5060
Where,
Example
Create a script file and type the following code into it −
a = [9 8 7; 6 5 4; 3 2 1];
b = [1 2 3; 4 5 6; 7 8 9];
c = cat(3, a, b, [ 2 3 1; 4 7 8; 3 9 0])
Array Functions
MATLAB provides the following functions to sort, rotate, permute, reshape,
or shift array contents.
Function Purpose
transpose Transpose
Examples
The following examples illustrate some of the functions mentioned above.
y = rand(3, 4, 5, 2);
numel(s) % no of elements in s
b =
7 8 9
1 2 3
4 5 6
c =
8 9 7
2 3 1
5 6 4
Sorting Arrays
Create a script file and type the following code into it −
sort(v) % sorting v
Cell Array
Cell arrays are arrays of indexed cells where each cell can store an array
of a different dimensions and data types.
The cell function is used for creating a cell array. Syntax for the cell
function is −
C = cell(dim)
C = cell(dim1,...,dimN)
D = cell(obj)
Where,
C is the cell array;
dim1, ... , dimN are scalar integers that specify the dimensions of C;
Example
Create a script file and type the following code into it −
c = cell(2, 5);
Enclosing the indices in braces {}, to refer to the data within individual cells
When you enclose the indices in first bracket, it refers to the set of cells.
For example −
c(1:2,1:2)
MATLAB will execute the above statement and return the following result
−
ans =
{
[1,1] = Red
[2,1] = 1
[1,2] = Blue
[2,2] = 2
}
You can also access the contents of cells by indexing with curly braces.
For example −
c{1, 2:4}
MATLAB will execute the above statement and return the following result
−
ans = Blue
ans = Green
ans = Yellow
MATLAB - Colon Notation
If you want to create a row vector, containing integers from 1 to 10, you
write −
1:10
MATLAB executes the statement and returns a row vector containing the
integers from 1 to 10 −
ans =
1 2 3 4 5 6 7 8 9 10
If you want to specify an increment value other than one, for example −
100: -5: 50
0:pi/8:pi
You can use the colon operator to create a vector of indices to select rows,
columns or elements of arrays.
The following table describes its use for this purpose (let us have a matrix
A) −
Format Purpose
A(:,:) is the equivalent two-dimensional array. For matrices this is the same
as A.
A(:) is all the elements of A, regarded as a single column. On the left side
of an assignment statement, A(:) fills A, preserving its shape from
before. In this case, the right side must contain the same number of
elements as A.
Example
Create a script file and type the following code in it −
A = [1 2 3 4; 4 5 6 7; 7 8 9 10]
A(2:3,2:3) % second and third rows and second and third columns
ans =
2
5
8
ans =
2 3
5 6
8 9
ans =
5 6
8 9
MATLAB - Numbers
MATLAB supports various numeric classes that include signed and
unsigned integers and single-precision and double-precision floating-point
numbers. By default, MATLAB stores all numeric values as double-
precision floating point numbers.
Function Purpose
Example
Create a script file and type the following code −
x =
x =
38 23 45
x =
38 23 45
x =
38 23 45
x =
38 23 45
Example
Let us extend the previous example a little more. Create a script file and
type the following code −
x = num2cell(x)
38 23 45
x =
38 23 45
x =
{
[1,1] = 38
[1,2] = 23
[1,3] = 45
}
Both the functions take the integer data type as the argument, for
example, intmax(int8) or intmin(int64) and return the maximum and
minimum values that you can represent with the integer data type.
Example
The following example illustrates how to obtain the smallest and largest
values of integers. Create a script file and write the following code in it −
Both the functions when called with the argument 'single', return the
maximum and minimum values that you can represent with the single-
precision data type and when called with the argument 'double', return the
maximum and minimum values that you can represent with the double-
precision data type.
Example
The following example illustrates how to obtain the smallest and largest
floating point numbers. Create a script file and write the following code in
it −
realmin('single'), realmax('single'))
realmin('double'), realmax('double'))
MATLAB - Strings
Creating a character string is quite simple in MATLAB. In fact, we have
used it many times. For example, you type the following in the command
prompt −
MATLAB will execute the above statement and return the following result
−
my_string = Tutorials Point
whos
MATLAB will execute the above statement and return the following result
−
Name Size Bytes Class Attributes
my_string 1x16 32 char
Example
Create a script file and type the following code into it −
str_back_to_char= char(str_ascii)
str_back_to_char = char(str_16bit)
84 117 116 111 114 105 97 108 39 115 32 80 111 105 110 116
84 117 116 111 114 105 97 108 39 115 32 80 111 105 110 116
str_back_to_char = Tutorial's Point
Using the MATLAB concatenation operator [] and separating each row with a
semicolon (;). Please note that in this method each row must contain the same
number of characters. For strings with different lengths, you should pad with
space characters as needed.
Using the char function. If the strings are of different lengths, char pads the
shorter strings with trailing blanks so that each row has the same number of
characters.
Example
Create a script file and type the following code into it −
Using the string concatenation function, strcat. This method removes trailing
spaces in the inputs.
Example
Create a script file and type the following code into it −
MATLAB cell array can hold different sizes and types of data in an array.
Cell arrays provide a more flexible way to store strings of varying length.
The cellstr function converts a character array into a cell array of strings.
Example
Create a script file and type the following code into it −
disp(profile)
Function Purpose
Examples
The following examples illustrate some of the above-mentioned string
functions −
FORMATTING STRINGS
Create a script file and type the following code into it −
A = pi*1000*ones(1,5);
JOINING STRINGS
Create a script file and type the following code into it −
first_names = strtok(students)
COMPARING STRINGS
Create a script file and type the following code into it −
if (strcmp(str1, str2))
else
end
MATLAB - Functions
A function is a group of statements that together perform a task. In
MATLAB, functions are defined in separate files. The name of the file and
of the function should be the same.
Functions can accept more than one input arguments and may return more
than one output arguments.
Example
The following function named mymax should be written in a file
named mymax.m. It takes five numbers as argument and returns the
maximum of the numbers.
Create a function file, named mymax.m and type the following code in it
−
max = n1;
max = n2;
end
max = n3;
end
max = n4;
end
end
The first line of a function starts with the keyword function. It gives the
name of the function and order of arguments. In our example,
the mymax function has five input arguments and one output argument.
The comment lines that come right after the function statement provide
the help text. These lines are printed when you type −
help mymax
MATLAB will execute the above statement and return the following result
−
This function calculates the maximum of the
five numbers given as input
MATLAB will execute the above statement and return the following result
−
ans = 89
Anonymous Functions
An anonymous function is like an inline function in traditional programming
languages, defined within a single MATLAB statement. It consists of a
single MATLAB expression and any number of input and output arguments.
You can define an anonymous function right at the MATLAB command line
or within a function or script.
This way you can create simple functions without having to create a file
for them.
Example
In this example, we will write an anonymous function named power, which
will take two numbers as input and return first number raised to the power
of the second number.
result1 = power(7, 3)
Primary functions can be called from outside of the file that defines them,
either from command line or from other functions, but sub-functions
cannot be called from command line or other functions, outside the
function file.
Sub-functions are visible only to the primary function and other sub-
functions within the function file that defines them.
Example
Let us write a function named quadratic that would calculate the roots of
a quadratic equation. The function would take three inputs, the quadratic
co-efficient, the linear co-efficient and the constant term. It would return
the roots.
%constant term
d = disc(a,b,c);
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);
quadratic(2,4,-4)
MATLAB will execute the above statement and return the following result
−
ans = 0.7321
Nested Functions
You can define functions within the body of another function. These are
called nested functions. A nested function contains any or all of the
components of any other function.
Nested functions are defined within the scope of another function and they
share access to the containing function's workspace.
Example
Let us rewrite the function quadratic, from previous example, however,
this time the disc function will be a nested function.
d = sqrt(b^2 - 4*a*c);
disc;
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);
quadratic2(2,4,-4)
MATLAB will execute the above statement and return the following result
−
ans = 0.73205
Private Functions
A private function is a primary function that is visible only to a limited
group of other functions. If you do not want to expose the implementation
of a function(s), you can create them as private functions.
Example
Let us rewrite the quadratic function. This time, however, the disc function
calculating the discriminant, will be a private function.
% a quadratic equation.
%constant term
d = disc(a,b,c);
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);
quadratic3(2,4,-4)
MATLAB will execute the above statement and return the following result
−
ans = 0.73205
Global Variables
Global variables can be shared by more than one function. For this, you
need to declare the variable as global in all the functions.
If you want to access that variable from the base workspace, then declare
the variable at the command line.
The global declaration must occur before the variable is actually used in a
function. It is a good practice to use capital letters for the names of global
variables to distinguish them from other variables.
Example
Let us create a function file named average.m and type the following code
in it −
global TOTAL
avg = sum(nums)/TOTAL;
end
global TOTAL;
TOTAL = 10;
n = [34, 45, 25, 45, 33, 19, 40, 34, 38, 42];
av = average(n)
When you run the file, it will display the following result −
av = 35.500
1 A = importdata(filename)
2 A = importdata('-pastespecial')
Loads data from the system clipboard rather than from a file.
3 A = importdata(___, delimiterIn)
Loads data from ASCII file, filename, or the clipboard, reading numeric
data starting from line headerlinesIn+1.
Returns the detected delimiter character for the input ASCII file
in delimiterOut and the detected number of header lines
in headerlinesOut, using any of the input arguments in the previous
syntaxes.
By default, Octave does not have support for importdata() function, so you will have
to search and install this package to make following examples work with your
Octave installation.
Example 1
Let us load and display an image file. Create a script file and type the
following code in it −
filename = 'smile.jpg';
A = importdata(filename);
image(A);
When you run the file, MATLAB displays the image file. However, you must
store it in the current directory.
Example 2
In this example, we import a text file and specify Delimiter and Column
Header. Let us create a space-delimited ASCII file with column headers,
named weeklydata.txt.
filename = 'weeklydata.txt';
headerlinesIn = 1;
A = importdata(filename,delimiterIn,headerlinesIn);
% View data
for k = [1:7]
disp(A.colheaders{1, k})
disp(A.data(:, k))
disp(' ')
end
MonDay
76.2100
45.6500
41.8500
82.1400
44.4700
TuesDay
61.5400
79.1900
92.1800
73.8200
57.6300
WednesDay
40.5700
93.5500
91.6900
41.0300
89.3600
ThursDay
55.7900
75.2900
81.3200
0.9900
13.8900
FriDay
70.2800
69.8700
90.3800
67.2200
19.8800
SaturDay
81.5300
74.6800
74.5100
93.1800
46.6000
Example 3
In this example, let us import data from clipboard.
Mathematics is simple
MATLAB provides the following functions for read and write operations at
the byte or character level −
Function Description
The fgetl and fgets functions read one line of a file at a time, where a newline
character separates each line.
The fread function reads a stream of data at the byte or bit level.
Example
We have a text data file 'myfile.txt' saved in our working directory. The file
stores rainfall data for three months; June, July and August for the year
2012.
The data in myfile.txt contains repeated sets of time, month and rainfall
measurements at five places. The header data stores the number of
months M; so we have M sets of measurements.
M = 3
12:00:00
June-2012
17.21 28.52 39.78 16.55 23.67
19.15 0.35 17.57 NaN 12.01
17.92 28.49 17.40 17.06 11.09
9.59 9.33 NaN 0.31 0.23
10.46 13.17 NaN 14.89 19.33
20.97 19.50 17.65 14.45 14.00
18.23 10.34 17.95 16.46 19.34
09:10:02
July-2012
12.76 16.94 14.38 11.86 16.89
20.46 23.17 NaN 24.89 19.33
30.97 49.50 47.65 24.45 34.00
18.23 30.34 27.95 16.46 19.34
30.46 33.17 NaN 34.89 29.33
30.97 49.50 47.65 24.45 34.00
28.67 30.34 27.95 36.46 29.34
15:03:40
August-2012
17.09 16.55 19.59 17.25 19.22
17.54 11.45 13.48 22.55 24.01
NaN 21.19 25.85 25.05 27.21
26.79 24.98 12.23 16.99 18.67
17.54 11.45 13.48 22.55 24.01
NaN 21.19 25.85 25.05 27.21
26.79 24.98 12.23 16.99 18.67
We will import data from this file and display this data. Take the following
steps −
Open the file with fopen function and get the file identifier.
Describe the data in the file with format specifiers, such as '%s' for a string,
'%d' for an integer, or '%f' for a floating-point number.
To skip literal characters in the file, include them in the format description. To
skip a data field, use an asterisk ('*') in the specifier.
For example, to read the headers and return the single value for M, we write
−
By default, fscanf reads data according to our format description until it does
not find any match for the data, or it reaches the end of the file. Here we will
use for loop for reading 3 sets of data and each time, it will read 7 rows and
5 columns.
We will create a structure named mydata in the workspace to store data read
from the file. This structure has three fields - time, month, and raindata array.
filename = '/data/myfile.txt';
rows = 7;
cols = 5;
fid = fopen(filename);
% read the file headers, find M (number of months)
for n = 1:M
mydata(n).raindata = ...
end
for n = 1:M
disp(mydata(n).time), disp(mydata(n).month)
disp(mydata(n).raindata)
end
fclose(fid);
09:10:02
July-2012
12.7600 NaN 34.0000 33.1700 24.4500
16.9400 24.8900 18.2300 NaN 34.0000
14.3800 19.3300 30.3400 34.8900 28.6700
11.8600 30.9700 27.9500 29.3300 30.3400
16.8900 49.5000 16.4600 30.9700 27.9500
20.4600 47.6500 19.3400 49.5000 36.4600
23.1700 24.4500 30.4600 47.6500 29.3400
15:03:40
August-2012
17.0900 13.4800 27.2100 11.4500 25.0500
16.5500 22.5500 26.7900 13.4800 27.2100
19.5900 24.0100 24.9800 22.5500 26.7900
17.2500 NaN 12.2300 24.0100 24.9800
19.2200 21.1900 16.9900 NaN 12.2300
17.5400 25.8500 18.6700 21.1900 16.9900
11.4500 25.0500 17.5400 25.8500 18.6700
Diary (or log) file of keystrokes and the resulting text output.
There are two ways to export a numeric array as a delimited ASCII data
file −
Example
The following example demonstrates the concept. Create a script file and
type the following code −
num_array = [ 1 2 3 4 ; 4 5 6 7; 7 8 9 0];
type array_data1.out
type array_data2.out
When you run the file, it displays the following result −
1.0000000e+00 2.0000000e+00 3.0000000e+00 4.0000000e+00
4.0000000e+00 5.0000000e+00 6.0000000e+00 7.0000000e+00
7.0000000e+00 8.0000000e+00 9.0000000e+00 0.0000000e+00
1 2 3 4
4 5 6 7
7 8 9 0
Please note that the save -ascii command and the dlmwrite function does
not work with cell arrays as input. To create a delimited ASCII file from
the contents of a cell array, you can
Either, convert the cell array to a matrix using the cell2mat function
If you use the save function to write a character array to an ASCII file, it
writes the ASCII equivalent of the characters to the file.
h = 'hello';
type textdata.out
MATLAB executes the above statements and displays the following result.
which is the characters of the string 'hello' in 8-digit ASCII format.
1.0400000e+02 1.0100000e+02 1.0800000e+02 1.0800000e+02 1.1100000e+02
diary
Optionally, you can give the name of the log file, say −
diary logdata.out
After processing the file, you need to close it with fclose(fid) function.
Example
Create a script file and type the following code in it −
x = 0:10:100;
y = [x; log(x)];
% Table Header
fclose(fid);
0.000000 -Inf
10.000000 2.302585
20.000000 2.995732
30.000000 3.401197
40.000000 3.688879
50.000000 3.912023
60.000000 4.094345
70.000000 4.248495
80.000000 4.382027
90.000000 4.499810
100.000000 4.605170
MATLAB - Plotting
To plot the graph of a function, you need to take the following steps −
Define x, by specifying the range of values for the variable x, for which the
function is to be plotted
Following example would demonstrate the concept. Let us plot the simple
function y = x for the range of values for x from 0 to 100, with an
increment of 5.
x = [0:5:100];
y = x;
plot(x, y)
When you run the file, MATLAB displays the following plot −
Let us take one more example to plot the function y = x2. In this example,
we will draw two graphs with the same function, but in second time, we
will reduce the value of increment. Please note that as we decrease the
increment, the graph becomes smoother.
x = [1 2 3 4 5 6 7 8 9 10];
x = [-100:20:100];
y = x.^2;
plot(x, y)
When you run the file, MATLAB displays the following plot −
x = [-100:5:100];
y = x.^2;
plot(x, y)
The xlabel and ylabel commands generate labels along x-axis and y-axis.
The grid on command allows you to put the grid lines on the graph.
The axis equal command allows generating the plot with the same scale
factors and the spaces on both axes.
Example
Create a script file and type the following code −
x = [0:0.01:10];
y = sin(x);
Example
Create a script file and type the following code −
x = [0 : 0.01: 10];
y = sin(x);
g = cos(x);
Code Color
w White
k Black
b Blue
r Red
c Cyan
g Green
m Magenta
y Yellow
Example
Let us draw the graph of two polynomials
g(x) = 5x3 + 9x + 2
g = 5 * x.^3 + 9 * x + 2;
When you run the file, MATLAB generates the following graph −
Setting Axis Scales
The axis command allows you to set the axis scales. You can provide
minimum and maximum values for x and y axes using the axis command
in the following way −
axis ( [xmin xmax ymin ymax] )
Example
Create a script file and type the following code −
x = [0 : 0.01: 10];
When you run the file, MATLAB generates the following graph −
Generating Sub-Plots
When you create an array of plots in the same figure, each of these plots
is called a subplot. The subplot command is used for creating subplots.
where, m and n are the number of rows and columns of the plot array
and pspecifies where to put a particular plot.
Each plot created with the subplot command can have its own
characteristics. Following example demonstrates the concept −
Example
Let us generate two plots −
y = e−1.5xsin(10x)
y = e−2xsin(10x)
x = [0:0.01:5];
y = exp(-1.5*x).*sin(10*x);
subplot(1,2,1)
y = exp(-2*x).*sin(10*x);
subplot(1,2,2)
plot(x,y),xlabel('x'),ylabel('exp(–2x)*sin(10x)'),axis([0 5 -1 1])
When you run the file, MATLAB generates the following graph −
MATLAB - Graphics
This chapter will continue exploring the plotting and graphics capabilities
of MATLAB. We will discuss −
Drawing bar charts
Drawing contours
Example
Let us have an imaginary classroom with 10 students. We know the
percent of marks obtained by these students are 75, 58, 90, 87, 50, 85,
92, 75, 60 and 95. We will draw the bar chart for this data.
x = [1:10];
y = [75, 58, 90, 87, 50, 85, 92, 75, 60, 95];
bar(x,y), xlabel('Student'),ylabel('Score'),
title('First Sem:')
When you run the file, MATLAB displays the following bar chart −
Drawing Contours
A contour line of a function of two variables is a curve along which the
function has a constant value. Contour lines are used for creating contour
maps by joining points of equal elevation above a given level, such as
mean sea level.
Example
Let us generate a contour map that shows the contour lines for a given
function g = f(x, y). This function has two variables. So, we will have to
generate two independent variables, i.e., two data sets x and y. This is
done by calling the meshgrid command.
When you run the file, MATLAB displays the following contour map −
Let us modify the code a little to spruce up the map
set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2)
When you run the file, MATLAB displays the following contour map −
As before, to define g, we first create a set of (x,y) points over the domain
of the function using the meshgrid command. Next, we assign the
function itself. Finally, we use the surf command to create a surface plot.
Example
Let us create a 3D surface map for the function g = xe-(x2 + y2)
[x,y] = meshgrid(-2:.2:2);
g = x .* exp(-x.^2 - y.^2);
surf(x, y, g)
When you run the file, MATLAB displays the following 3-D map −
So far, we have seen that all the examples work in MATLAB as well as its
GNU, alternatively called Octave. But for solving basic algebraic equations,
both MATLAB and Octave are little different, so we will try to cover MATLAB
and Octave in separate sections.
MATLAB will execute the above statement and return the following result
−
ans =
5
y = solve('x-5 = 0')
MATLAB will execute the above statement and return the following result
−
y =
5
You may even not include the right hand side of the equation −
solve('x-5')
MATLAB will execute the above statement and return the following result
−
ans =
5
solve(equation, variable)
For example, let us solve the equation v – u – 3t2 = 0, for v. In this case,
we should write −
solve('v-u-3*t^2=0', 'v')
MATLAB will execute the above statement and return the following result
−
ans =
3*t^2 + u
roots([1, -5])
Octave will execute the above statement and return the following result −
ans = 5
y = roots([1, -5])
Octave will execute the above statement and return the following result −
y = 5
s = solve(eq);
solve('(x-3)^2*(x-7)=0')
MATLAB will execute the above statement and return the following result
−
ans =
3
3
7
In case of higher order equations, roots are long containing many terms.
You can get the numerical value of such roots by converting them to
double. The following example solves the fourth order equation x4 − 7x3 +
3x2 − 5x + 9 = 0.
s = solve(eq);
Please note that the last two roots are complex numbers.
Solving Higher Order Equations in Octave
The following example solves the fourth order equation x4 − 7x3 + 3x2 −
5x + 9 = 0.
s = roots(v);
5x + 9y = 5
3x – 6y = 4
s.x
s.y
When you run the file, it displays the following result −
ans =
22/19
ans =
-5/57
In same way, you can solve larger linear systems. Consider the following
set of equations −
x + 3y -2z = 5
3x + 5y + 6z = 7
2x + 4y + 3z = 8
5x + 9y = 5
3x – 6y = 4
A = [5, 9; 3, -6];
b = [5;4];
A \ b
1.157895
-0.087719
In same way, you can solve larger linear systems as given below −
x + 3y -2z = 5
3x + 5y + 6z = 7
2x + 4y + 3z = 8
Expanding and Collecting Equations in
MATLAB
The expand and the collect function expands and collects an equation
respectively. The following example demonstrates the concepts −
When you work with many symbolic functions, you should declare that
your variables are symbolic.
% expanding equations
expand((x-5)*(x+9))
expand((x+2)*(x-3)*(x-5)*(x+7))
expand(sin(2*x))
expand(cos(x+y))
% collecting equations
collect(x^3 *(x-7))
collect(x^4*(x-3)*(x-5))
When you work with many symbolic functions, you should declare that
your variables are symbolic but Octave has different approach to define
symbolic variables. Notice the use of Sin and Cos, which are also defined
in symbolic package.
symbols
x = sym ('x');
y = sym ('y');
z = sym ('z');
% expanding equations
expand((x-5)*(x+9))
expand((x+2)*(x-3)*(x-5)*(x+7))
expand(Sin(2*x))
expand(Cos(x+y))
% collecting equations
collect(x^3 *(x-7), z)
collect(x^4*(x-3)*(x-5), z)
When you run the file, it displays the following result −
ans =
-45.0+x^2+(4.0)*x
ans =
210.0+x^4-(43.0)*x^2+x^3+(23.0)*x
ans =
sin((2.0)*x)
ans =
cos(y+x)
ans =
x^(3.0)*(-7.0+x)
ans =
(-3.0+x)*x^(4.0)*(-5.0+x)
Example
Create a script file and type the following code −
syms x
syms y
factor(x^3 - y^3)
factor([x^2-y^2,x^3+y^3])
simplify((x^4-16)/(x^2-4))
This chapter will deal with problems of calculus. In this chapter, we will
discuss pre-calculus concepts i.e., calculating limits of functions and
verifying the properties of limits.
Calculating Limits
MATLAB provides the limit function for calculating limits. In its most basic
form, the limit function takes expression as an argument and finds the
limit of the expression as the independent variable goes to zero.
For example, let us calculate the limit of a function f(x) = (x3 + 5)/(x4 +
7), as x tends to zero.
syms x
MATLAB will execute the above statement and return the following result
−
ans =
5/7
The limit function falls in the realm of symbolic computing; you need to
use the syms function to tell MATLAB which symbolic variables you are
using. You can also compute limit of a function, as the variable tends to
some number other than zero. To calculate lim x->a(f(x)), we use the limit
command with arguments. The first being the expression and the second
is the number, that x approaches, here it is a.
For example, let us calculate limit of a function f(x) = (x-3)/(x-1), as x
tends to 1.
limit((x - 3)/(x-1),1)
MATLAB will execute the above statement and return the following result
−
ans =
NaN
limit(x^2 + 5, 3)
MATLAB will execute the above statement and return the following result
−
ans =
14
symbols
x = sym("x");
subs((x^3+5)/(x^4+7),x,0)
Octave will execute the above statement and return the following result −
ans =
0.7142857142857142857
g(x) = x2 + 1.
Example
Create a script file and type the following code into it −
syms x
f = (3*x + 5)/(x-3);
g = x^2 + 1;
l1 = limit(f, 4)
l2 = limit (g, 4)
lAdd = limit(f + g, 4)
lSub = limit(f - g, 4)
lMult = limit(f*g, 4)
l2 =
17
lAdd =
34
lSub =
0
lMult =
289
lDiv =
1
symbols
x = sym("x");
f = (3*x + 5)/(x-3);
g = x^2 + 1;
l1 = subs(f, x, 4)
l2 = subs (g, x, 4)
Octave will execute the above statement and return the following result −
l1 =
17.0
l2 =
17.0
lAdd =
34.0
lSub =
0.0
lMult =
289.0
lDiv =
1.0
f(x) = (x - 3)/|x - 3|
We will show that limx->3 f(x) does not exist. MATLAB helps us to establish
this fact in two ways −
Example
Create a script file and type the following code into it −
f = (x - 3)/abs(x-3);
ezplot(f,[-1,5])
l = limit(f,x,3,'left')
r = limit(f,x,3,'right')
When you run the file, MATLAB draws the following plot
After this following output is displayed −
l =
-1
r =
1
MATLAB - Differential
For example, let us compute the derivative of the function f(t) = 3t2 + 2t-
2
Example
Create a script file and type the following code into it −
syms t
f = 3*t^2 + 2*t^(-2);
diff(f)
When the above code is compiled and executed, it produces the following
result −
ans =
6*t - 4/t^3
symbols
t = sym("t");
f = 3*t^2 + 2*t^(-2);
differentiate(f,t)
Rule 1
For any functions f and g and any real numbers a and b are the derivative
of the function −
Rule 2
The sum and subtraction rules state that if f and g are two functions, f'
and g' are their derivatives respectively, then,
Rule 3
The product rule states that if f and g are two functions, f' and g' are their
derivatives respectively, then,
Rule 4
The quotient rule states that if f and g are two functions, f' and g' are
their derivatives respectively, then,
Rule 5
The polynomial or elementary power rule states that, if y = f(x) = xn,
then f' = n. x(n-1)
A direct outcome of this rule is that the derivative of any constant is zero,
i.e., if y = k, any constant, then
f' = 0
Rule 6
The chain rule states that, derivative of the function of a function h(x) =
f(g(x)) with respect to x is,
h'(x)= f'(g(x)).g'(x)
Example
Create a script file and type the following code into it −
syms x
syms t
f = (x + 2)*(x^2 + 3)
der1 = diff(f)
der2 = diff(f)
der3 = diff(f)
f = (2*x^2 + 3*x)/(x^3 + 1)
der4 = diff(f)
f = (x^2 + 1)^17
der5 = diff(f)
der6 = diff(f)
When you run the file, MATLAB displays the following result −
f =
(x^2 + 3)*(x + 2)
der1 =
2*x*(x + 2) + x^2 + 3
f =
(t^(1/2) + t^3)*(t^2 + 3)
der2 =
(t^2 + 3)*(3*t^2 + 1/(2*t^(1/2))) + 2*t*(t^(1/2) + t^3)
f =
(x^2 - 2*x + 1)*(3*x^3 - 5*x^2 + 2)
der3 =
(2*x - 2)*(3*x^3 - 5*x^2 + 2) - (- 9*x^2 + 10*x)*(x^2 - 2*x + 1)
f =
(2*x^2 + 3*x)/(x^3 + 1)
der4 =
(4*x + 3)/(x^3 + 1) - (3*x^2*(2*x^2 + 3*x))/(x^3 + 1)^2
f =
(x^2 + 1)^17
der5 =
34*x*(x^2 + 1)^16
f =
1/(t^3 + 3*t^2 + 5*t - 9)^6
der6 =
-(6*(3*t^2 + 6*t + 5))/(t^3 + 3*t^2 + 5*t - 9)^7
symbols
x = sym("x");
t = sym("t");
f = (x + 2)*(x^2 + 3)
der1 = differentiate(f,x)
der2 = differentiate(f,t)
der3 = differentiate(f,x)
f = (2*x^2 + 3*x)/(x^3 + 1)
der4 = differentiate(f,x)
f = (x^2 + 1)^17
der5 = differentiate(f,x)
der6 = differentiate(f,t)
(2.0+x)*(3.0+x^(2.0))
der1 =
3.0+x^(2.0)+(2.0)*(2.0+x)*x
f =
(t^(3.0)+sqrt(t))*(3.0+t^(2.0))
der2 =
(2.0)*(t^(3.0)+sqrt(t))*t+((3.0)*t^(2.0)+(0.5)*t^(-0.5))*(3.0+t^(2.0))
f =
(1.0+x^(2.0)-(2.0)*x)*(2.0-(5.0)*x^(2.0)+(3.0)*x^(3.0))
der3 =
(-2.0+(2.0)*x)*(2.0-(5.0)*x^(2.0)+(3.0)*x^(3.0))+((9.0)*x^(2.0)-(10.0)*x)*(1.0+x^(2.0)-
(2.0)*x)
f =
(1.0+x^(3.0))^(-1)*((2.0)*x^(2.0)+(3.0)*x)
der4 =
(1.0+x^(3.0))^(-1)*(3.0+(4.0)*x)-(3.0)*(1.0+x^(3.0))^(-2)*x^(2.0)*((2.0)*x^(2.0)+(3.0)*x)
f =
(1.0+x^(2.0))^(17.0)
der5 =
(34.0)*(1.0+x^(2.0))^(16.0)*x
f =
(-9.0+(3.0)*t^(2.0)+t^(3.0)+(5.0)*t)^(-6.0)
der6 =
-(6.0)*(-9.0+(3.0)*t^(2.0)+t^(3.0)+(5.0)*t)^(-7.0)*(5.0+(3.0)*t^(2.0)+(6.0)*t)
Derivatives of Exponential, Logarithmic
and Trigonometric Functions
The following table provides the derivatives of commonly used
exponential, logarithmic and trigonometric functions −
Function Derivative
ex ex
ln x 1/x
lncx 1/x.ln c
xx xx.(1 + ln x)
sin(x) cos(x)
cos(x) -sin(x)
sec(x) sec(x).tan(x)
csc(x) -csc(x).cot(x)
Example
Create a script file and type the following code into it −
syms x
y = exp(x)
diff(y)
y = x^9
diff(y)
y = sin(x)
diff(y)
y = tan(x)
diff(y)
y = cos(x)
diff(y)
y = log(x)
diff(y)
y = log10(x)
diff(y)
y = sin(x)^2
diff(y)
y = cos(3*x^2 + 2*x + 1)
diff(y)
y = exp(x)/sin(x)
diff(y)
When you run the file, MATLAB displays the following result −
y =
exp(x)
ans =
exp(x)
y =
x^9
ans =
9*x^8
y =
sin(x)
ans =
cos(x)
y =
tan(x)
ans =
tan(x)^2 + 1
y =
cos(x)
ans =
-sin(x)
y =
log(x)
ans =
1/x
y =
log(x)/log(10)
ans =
1/(x*log(10))
y =
sin(x)^2
ans =
2*cos(x)*sin(x)
y =
cos(3*x^2 + 2*x + 1)
ans =
-sin(3*x^2 + 2*x + 1)*(6*x + 2)
y =
exp(x)/sin(x)
ans =
exp(x)/sin(x) - (exp(x)*cos(x))/sin(x)^2
symbols
x = sym("x");
y = Exp(x)
differentiate(y,x)
y = x^9
differentiate(y,x)
y = Sin(x)
differentiate(y,x)
y = Tan(x)
differentiate(y,x)
y = Cos(x)
differentiate(y,x)
y = Log(x)
differentiate(y,x)
%y = Log10(x)
%differentiate(y,x)
y = Sin(x)^2
differentiate(y,x)
y = Cos(3*x^2 + 2*x + 1)
differentiate(y,x)
y = Exp(x)/Sin(x)
differentiate(y,x)
exp(x)
y =
x^(9.0)
ans =
(9.0)*x^(8.0)
y =
sin(x)
ans =
cos(x)
y =
tan(x)
ans =
1+tan(x)^2
y =
cos(x)
ans =
-sin(x)
y =
log(x)
ans =
x^(-1)
y =
sin(x)^(2.0)
ans =
(2.0)*sin(x)*cos(x)
y =
cos(1.0+(2.0)*x+(3.0)*x^(2.0))
ans =
-(2.0+(6.0)*x)*sin(1.0+(2.0)*x+(3.0)*x^(2.0))
y =
sin(x)^(-1)*exp(x)
ans =
sin(x)^(-1)*exp(x)-sin(x)^(-2)*cos(x)*exp(x)
f = x*exp(-3*x);
diff(f, 2)
MATLAB executes the code and returns the following result −
ans =
9*x*exp(-3*x) - 6*exp(-3*x)
symbols
x = sym("x");
f = x*Exp(-3*x);
differentiate(f, x, 2)
(9.0)*exp(-(3.0)*x)*x-(6.0)*exp(-(3.0)*x)
Example
In this example, let us solve a problem. Given that a function y = f(x) =
3 sin(x) + 7 cos(5x). We will have to find out whether the equation f"
+ f = -5cos(2x) holds true.
syms x
if(isequal(lhs,rhs))
else
end
symbols
x = sym("x");
if(lhs == rhs)
else
end
For a function y = f(x) the points on the graph where the graph has zero
slope are called stationary points. In other words stationary points are
where f'(x) = 0.
syms x
ezplot(y)
symbols
x = sym('x');
ezplot(y)
syms x
symbols
x = sym('x');
g = diff(y)
MATLAB executes the code and returns the following result −
g =
6*x^2 + 6*x - 12
symbols
x = sym("x");
g = differentiate(y,x)
s = solve(g)
symbols
x = sym("x");
g = differentiate(y,x)
roots([6, 6, -12])
-12.0+(6.0)*x^(2.0)+(6.0)*x
ans =
-2
1
This agrees with our plot. So let us evaluate the function f at the
critical points x = 1, -2. We can substitute a value in a symbolic function
by using the subs command.
symbols
x = sym("x");
g = differentiate(y,x)
roots([6, 6, -12])
Therefore, The minimum and maximum values on the function f(x) = 2x3 +
3x2 − 12x + 17, in the interval [-2,2] are 10 and 37.
The most basic form of the dsolve command for finding the solution to a
single equation is
dsolve('eqn')
where eqn is a text string used to enter the equation.
You can also specify initial and boundary conditions for the problem, as
comma-delimited list following the equation as −
dsolve('eqn','cond1', 'cond2',…)
MATLAB - Integration
The second type of problems involve adding up a very large number of very
small quantities and then taking a limit as the size of the quantities approaches
zero, while the number of terms tend to infinity. This process leads to the
definition of the definite integral.
Definite integrals are used for finding area, volume, center of gravity,
moment of inertia, work done by a force, and in numerous other
applications.
In symbols −
∫ 2xdx = x2.
∫ 2xdx = x2 + c.
syms x
int(2*x)
MATLAB executes the above statement and returns the following result −
ans =
x^2
Example 1
In this example, let us find the integral of some commonly used
expressions. Create a script file and type the following code in it −
syms x n
int(sym(x^n))
f = 'sin(n*t)'
int(sym(f))
syms a t
int(a*cos(pi*t))
int(a^x)
Example 2
Create a script file and type the following code in it −
syms x n
int(cos(x))
int(exp(x))
int(log(x))
int(x^-1)
int(x^5*cos(5*x))
pretty(int(x^5*cos(5*x)))
int(x^-5)
int(sec(x)^2)
ans =
exp(x)
ans =
x*(log(x) - 1)
ans =
log(x)
ans =
(24*cos(5*x))/3125 + (24*x*sin(5*x))/625 - (12*x^2*cos(5*x))/125 + (x^4*cos(5*x))/5 -
(4*x^3*sin(5*x))/25 + (x^5*sin(5*x))/5
2 4
24 cos(5 x) 24 x sin(5 x) 12 x cos(5 x) x cos(5 x)
----------- + ------------- - -------------- + ------------
3125 625 125 5
3 5
4 x sin(5 x) x sin(5 x)
------------- + -----------
25 5
ans =
-1/(4*x^4)
ans =
tan(x)
2
x (3 x - 5 x + 1)
ans =
- (7*x^6)/12 - (3*x^5)/5 + (5*x^4)/8 + x^3/2
6 5 4 3
7 x 3 x 5 x x
- ---- - ---- + ---- + --
12 5 8 2
Finding Definite Integral Using MATLAB
By definition, definite integral is basically the limit of a sum. We use
definite integrals to find areas such as the area between a curve and the
x-axis and the area between two curves. Definite integrals can also be
used in other situations, where the quantity required can be expressed as
the limit of a sum.
The int function can be used for definite integration by passing the limits
over which you want to calculate the integral.
To calculate
we write,
int(x, a, b)
int(x, 4, 9)
MATLAB executes the above statement and returns the following result −
ans =
65/2
symbols
x = sym("x");
f = x;
c = [1, 0];
integral = polyint(c);
32.500
symbols
f = inline("x");
Example 1
Let us calculate the area enclosed between the x-axis, and the curve y =
x3−2x+5 and the ordinates x = 1 and x = 2.
a = int(f, 1, 2)
x = sym("x");
integral = polyint(c);
5.7500
symbols
x = sym("x");
Example 2
Find the area under the curve: f(x) = x2 cos(x) for −4 ≤ x ≤ 9.
f = x^2*cos(x);
ezplot(f, [-4,9])
a = int(f, -4, 9)
Area:
0.3326
symbols
x = sym("x");
f = inline("x^2*cos(x)");
ezplot(f, [-4,9])
MATLAB - Polynomials
p = [1 7 0 -5 9];
Evaluating Polynomials
The polyval function is used for evaluating a polynomial at a specified
value. For example, to evaluate our previous polynomial p, at x = 4, type
−
p = [1 7 0 -5 9];
polyval(p,4)
MATLAB executes the above statements and returns the following result −
ans = 693
For example, let us create a square matrix X and evaluate the polynomial
p, at X −
p = [1 7 0 -5 9];
X = [1 2 -3 4; 2 -5 6 3; 3 1 0 2; 5 -7 3 8];
polyvalm(p, X)
MATLAB executes the above statements and returns the following result −
ans =
2307 -1769 -939 4499
2314 -2376 -249 4695
2256 -1892 -549 4310
4570 -4532 -1062 9269
p = [1 7 0 -5 9];
r = roots(p)
MATLAB executes the above statements and returns the following result −
r =
-6.8661 + 0.0000i
-1.4247 + 0.0000i
0.6454 + 0.7095i
0.6454 - 0.7095i
The function poly is an inverse of the roots function and returns to the
polynomial coefficients. For example −
p2 = poly(r)
MATLAB executes the above statements and returns the following result −
p2 =
Columns 1 through 3:
Columns 4 and 5:
p = polyfit(x,y,n)
Example
Create a script file and type the following code −
% and plot the estimate over the real data values for comparison:
x2 = 1:.1:6;
y2 = polyval(p,x2);
plot(x,y,'o',x2,y2)
grid on
When you run the file, MATLAB displays the following result −
p =
4.1056 -47.9607 222.2598 -362.7453 191.1250
MATLAB provides the laplace, fourier and fft commands to work with
Laplace, Fourier and Fast Fourier transforms.
Laplace transform is also denoted as transform of f(t) to F(s). You can see
this transform or integration process converts f(t), a function of the
symbolic variable t, into another function F(s), with another variable s.
laplace(f(t))
Example
In this example, we will compute the Laplace transform of some commonly
used functions.
syms s t a b w
laplace(a)
laplace(t^2)
laplace(t^9)
laplace(exp(-b*t))
laplace(sin(w*t))
laplace(cos(w*t))
ans =
2/s^3
ans =
362880/s^10
ans =
1/(b + s)
ans =
w/(s^2 + w^2)
ans =
s/(s^2 + w^2)
For example,
ilaplace(1/s^3)
MATLAB will execute the above statement and display the result −
ans =
t^2/2
Example
Create a script file and type the following code −
syms s t a b w
ilaplace(1/s^7)
ilaplace(2/(w+s))
ilaplace(s/(s^2+4))
ilaplace(exp(-b*t))
ilaplace(w/(s^2 + w^2))
ilaplace(s/(s^2 + w^2))
ans =
2*exp(-t*w)
ans =
cos(2*t)
ans =
ilaplace(exp(-b*t), t, x)
ans =
sin(t*w)
ans =
cos(t*w)
Example
Create a script file and type the following code in it −
syms x
When you run the file, MATLAB plots the following graph −
The following result is displayed −
FT =
(2^(1/2)*pi^(1/2)*exp(-w^2/8))/2
ezplot(FT)
f = ifourier(-2*exp(-abs(w)))
MATLAB will execute the above statement and display the result −
f =
-2/(pi*(x^2 + 1))
MATLAB - GNU Octave Tutorial
MATLAB vs Octave
Most MATLAB programs run in Octave, but some of the Octave programs
may not run in MATLAB because, Octave allows some syntax that MATLAB
does not.
For example, MATLAB supports single quotes only, but Octave supports
both single and double quotes for defining strings. If you are looking for a
tutorial on Octave, then kindly go through this tutorial from beginning
which covers both MATLAB as well as Octave.
Compatible Examples
Almost all the examples covered in this tutorial are compatible with
MATLAB as well as Octave. Let's try following example in MATLAB and
Octave which produces same result without any syntax changes −
This example creates a 3D surface map for the function g = xe -(x2 + y2).
[x,y] = meshgrid(-2:.2:2);
g = x .* exp(-x.^2 - y.^2);
surf(x, y, g)
Non-compatible Examples
Though all the core functionality of MATLAB is available in Octave, there
are some functionality for example, Differential & Integration Calculus,
which does not match exactly in both the languages. This tutorial has tried
to give both type of examples where they differed in their syntax.
f = x^2*cos(x);
ezplot(f, [-4,9])
a = int(f, -4, 9)
Area:
0.3326
But to give area of the same curve in Octave, you will have to make use
of symbolic package as follows −
symbols
x = sym("x");
f = inline("x^2*cos(x)");
ezplot(f, [-4,9])
Simulink supports −
system-level design
simulation
There are several other add-on products provided by MathWorks and third-
party hardware and software products that are available for use with
Simulink.
Using Simulink
To open Simulink, type in the MATLAB work space −
simulink
Simulink opens with the Library Browser. The Library Browser is used
for building simulation models.
On the left side window pane, you will find several libraries categorized on
the basis of various systems, clicking on each one will display the design
blocks on the right window pane.
Building Models
To create a new model, click the New button on the Library Browser's
toolbar. This opens a new untitled model window.
Model elements are added by selecting the appropriate elements from the
Library Browser and dragging them into the Model window.
Alternately, you can copy the model elements and paste them into the
model window.
Examples
Drag and drop items from the Simulink library to make your project.
For the purpose of this example, two blocks will be used for the simulation
- A Source (a signal) and a Sink (a scope). A signal generator (the
source) generates an analog signal, which will then be graphically
visualized by the scope(the sink).
Begin by dragging the required blocks from the library to the project
window. Then, connect the blocks together which can be done by dragging
connectors from connection points on one block to those of another.
Drag a signal line from the output of the Sine Wave block to the input of
the Scope block.
Run the simulation by pressing the 'Run' button, keeping all parameters
default (you can change them from the Simulation menu)