MATLAB
MATLAB
MATLAB
What is MATLAB?
MATLAB, short for Matrix Laboratory, is a powerful programming language and interactive environment
primarily designed for numerical computations, data visualization, and algorithm development. It's widely
used in engineering, science, and mathematics due to its efficiency in handling matrices and arrays.
Key Features
Matrix-based computations: MATLAB excels at performing operations on matrices and arrays,
making it ideal for linear algebra, statistics, and data analysis.
Interactive environment: Allows for rapid prototyping, testing, and exploration of ideas.
Extensive function library: Offers a vast collection of built-in functions for various mathematical
operations, signal processing, image processing, and more.
Visualization capabilities: Provides tools for creating high-quality plots, graphs, and animations to
visualize data and results.
Programming language: Supports both scripting and procedural programming, enabling the
development of complex algorithms and applications.
Integration with other languages: Interacts with C, C++, Java, and other languages for extending
functionality.
Common Applications
Signal processing: Analyzing and manipulating signals from various sources.
Image and video processing: Processing, analyzing, and manipulating images and videos.
Control systems: Designing and simulating control systems for various applications.
Financial modeling: Building models for financial analysis and risk assessment.
Machine learning: Developing and training machine learning models.
Scientific computing: Performing complex numerical computations in various scientific fields.
Getting Started
To use MATLAB, you can:
Purchase a MATLAB license: Provides access to the full suite of features and support.
Use a MATLAB student version: Offers a discounted version for students.
Try MATLAB Online: Access a web-based version for basic tasks.
Use Octave: A free, open-source alternative with similar syntax and capabilities.
MATLAB Environment
The MATLAB environment consists of several key components:
Command Window: Where you enter commands and view results.
Workspace: Displays variables and their values.
Command History: Stores previously entered commands.
Current Folder: Shows the current working directory.
Editor: For creating and editing MATLAB code.
MATLAB as a Calculator
MATLAB operates like a sophisticated calculator. You input commands at the >> prompt and get immediate
results.
Example:
5+5 Output: ans = 10
Key Points:
Interpreted language: Commands are executed line by line.
Basic arithmetic: Use +, -, *, /, and ^ for addition, subtraction, multiplication, division, and
exponentiation, respectively.
Special values: pi for π, Inf for infinity, i or j for imaginary unit, NaN for 'Not a Number'.
Suppressing Output
Use a semicolon (;) at the end of a line to prevent MATLAB from displaying the result.
Example:
x = 3;
y = x + 5;
Adding Comments
Use the percentage symbol (%) to indicate a comment.
Example:
x = 9 % Assign the value 9 to x
Remember:
MATLAB is case-sensitive.
Indentation is not required but improves readability.
Common Errors
Syntax errors: Incorrect use of operators or symbols.
Runtime errors: Errors that occur during program execution, like division by zero.
Logical errors: Errors in the program's logic leading to incorrect results.
\ Left-division operator.
/ Right-division operator.
. Decimal point.
= Assignment operator.
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:
Understanding Variables
Fundamental data structure: In MATLAB, every variable is essentially an array or matrix.
Creation: Assign values using the equal sign (=).
Referencing: Use variable names in calculations and operations.
Default variable: If an expression's result isn't assigned, MATLAB stores it in ans.
Multiple assignments: Assign values to multiple variables in a single line.
Managing Variables
Listing variables:
o who: Displays variable names.
o whos: Provides detailed information about variables (size, type, memory usage).
Clearing variables:
o clear: Removes all variables from the workspace.
o clear var1 var2: Removes specific variables.
Controlling Number Display
Format options:
o format short: Default, four decimal places.
o format long: 16 digits.
o format bank: Two decimal places.
o format short e: Exponential notation with four decimal places.
o format long e: Exponential notation with 15 digits.
o format rat: Rational approximation.
Key Points:
MATLAB is dynamically typed, meaning you don't need to declare variable types beforehand.
Variable names are case-sensitive.
Efficient memory management is crucial for large datasets.
Creating Vectors and Matrices in MATLAB
Vectors
A vector is a one-dimensional array of numbers.
Row vectors:
Elements separated by spaces or commas.
Example:
row_vector = [1, 2, 3, 4];
Column vectors:
Elements separated by semicolons.
Example:
column_vector = [1; 2; 3; 4];
Matrices
A matrix is a two-dimensional array of numbers.
Rows are separated by semicolons.
Elements within a row are separated by spaces or commas.
Example:
matrix = [1 2 3; 4 5 6; 7 8 9];
m=
1 2 3
4 5 6
7 8 9
Key points:
MATLAB is a matrix-based language, so vectors are essentially special cases of matrices.
You can access individual elements using indexing (e.g., matrix(2,3) accesses the element in the
second row, third column).
There are other ways to create vectors and matrices, such as using the colon operator (:), built-in
functions like ones, zeros, and rand, and more.
MATLAB - Commands
MATLAB is an interactive program for numerical computation and data visualization. You can enter a
command by typing it at the MATLAB prompt '>>' on the Command Window.
In this section, we will provide lists of commonly used general MATLAB commands.
Commands for Managing a Session
MATLAB provides various commands for managing a session. The following table provides all such commands
−
Command Purpose
%d Format as an integer.
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 −
Clicking the Run button on the editor window or
Just typing the filename (without extension) in the command prompt: >> prog1
The command window prompt displays the result −
6170
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
uint8
2
8-bit unsigned integer
int16
3
16-bit signed integer
uint16
4
16-bit unsigned integer
int32
5
32-bit signed integer
uint32
6
32-bit unsigned integer
int64
7
64-bit signed integer
uint64
8
64-bit unsigned integer
single
9
single precision numerical data
double
10
double precision numerical data
logical
11
logical values of 1 or 0, represent true and false respectively
char
12
character data (strings are stored as vector of characters)
cell array
13
array of indexed cells, each capable of storing an array of a different dimension and data type
structure
14 C-like structures, each structure having named fields capable of storing an array of a different
dimension and data type
function handle
15
pointer to a function
16 user classes
objects constructed from a user-defined class
java classes
17
objects constructed from a Java class
Example
Create a script file with the following code −
mat2cell Convert array to cell array with potentially different sized cells
is Detect state
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)
When you run the file, it produces the following result −
x=3
ans = 0
ans = 1
ans = 1
ans = 1
ans = 1
x = 23.540
ans = 0
ans = 1
ans = 1
ans = 1
ans = 1
x=
1 2 3
ans = 0
ans = 1
ans = 1
ans = 0
x = Hello
ans = 0
ans = 0
ans = 1
ans = 0
ans = 0
Following is the general form of a typical decision making structure found in most of
the programming languages −
MATLAB provides following types of decision making statements. Click the following
links to check their detail −
Sr.N
Statement & Description
o.
if...else...end statement
2 An if statement can be followed by an optional else statement, which executes
when the boolean expression is false.
If... elseif...elseif...else...end statements
3 An if statement can be followed by one (or more) optional elseif... and
an else statement, which is very useful to test various conditions.
nested if statements
4
You can use one if or elseif statement inside another if or elseif statement(s).
switch statement
5
A switch statement allows a variable to be tested for equality against a list of values.
Programming languages provide various control structures that allow for more
complicated execution paths.
MATLAB provides following types of loops to handle looping requirements. Click the
following links to check their detail −
Sr.N
Loop Type & Description
o.
while loop
1 Repeats a statement or group of statements while a given condition is true. It tests the
condition before executing the loop body.
for loop
2 Executes a sequence of statements multiple times and abbreviates the code that manages
the loop variable.
nested loops
3
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.
Sr.N
Control Statement & Description
o.
break statement
1 Terminates the loop statement and transfers execution to the statement immediately
following the loop.
continue statement
2 Causes the loop to skip the remainder of its body and immediately retest its condition prior
to reiterating.
Vectors in MATLAB: A Comprehensive Overview
Creating Vectors
Row vectors: Elements separated by spaces or commas.
Matlab
row_vector = [1, 2, 3, 4];
Use code with caution.
Column vectors: Elements separated by semicolons.
Matlab
column_vector = [1; 2; 3; 4];
Use code with caution.
Accessing Vector Elements
Indexing: Use parentheses to access specific elements.
Matlab
v = [1, 2, 3, 4];
element = v(2); % Access the second element
Use code with caution.
Colon operator: Access a range of elements.
Matlab
sub_vector = v(2:4); % Access elements from index 2 to 4
Use code with caution.
Vector Operations
Addition and subtraction: Add or subtract corresponding elements of vectors with the same size.
Matlab
result = vector1 + vector2;
Use code with caution.
Scalar multiplication: Multiply each element by a scalar.
Matlab
result = scalar * vector;
Use code with caution.
Transpose: Convert a row vector to a column vector, or vice versa.
Matlab
column_vector = row_vector';
Use code with caution.
Appending vectors: Combine vectors using square brackets.
Matlab
new_vector = [vector1, vector2];
Use code with caution.
Magnitude: Calculate the length of a vector using the norm function.
Matlab
magnitude = norm(vector);
Use code with caution.
Dot product: Calculate the dot product of two vectors using the dot function.
Matlab
dot_product = dot(vector1, vector2);
Use code with caution.
Vectors with uniformly spaced elements: Use the colon operator to create vectors with evenly
spaced elements.
Matlab
vector = 1:2:10; % Create a vector from 1 to 10 with a step of 2
Use code with caution.
Visualizing Vectors
Plot: Use the plot function to visualize vectors graphically.
Matlab
plot(vector);
Use code with caution.
Example:
Matlab
% Create vectors
v1 = [1, 2, 3];
v2 = [4, 5, 6]';
% Vector operations
sum_vector = v1 + v2;
scalar_product = 2 * v1;
transpose_v2 = v2';
appended_vector = [v1, v2'];
magnitude_v1 = norm(v1);
dot_product_v1v2 = dot(v1, v2);
% Visualization
plot(v1);
MATLAB Matrices:
Creating and Manipulating Matrices
Creating a matrix: Enclose elements in square brackets, separating rows with semicolons and
elements within a row with spaces or commas.
Matlab
matrix = [1 2 3; 4 5 6; 7 8 9];
Accessing elements: Use indices to access individual elements.
Matlab
element = matrix(2, 3); % Access element at row 2, column 3
Accessing rows and columns: Use the colon operator to access entire rows or columns.
Matlab
row = matrix(2, :); % Access the second row
column = matrix(:, 3); % Access the third column
Creating submatrices: Use colon notation to specify ranges of rows and columns.
Matlab
submatrix = matrix(2:3, 1:2); % Create a 2x2 submatrix
Deleting rows or columns: Assign an empty matrix ([]) to the desired row or column.
Matlab
matrix(:, 3) = []; % Delete the third column
Matrix Operations
Addition and subtraction: Add or subtract corresponding elements of matrices with the same
dimensions.
Matlab
result = matrix1 + matrix2;
Scalar multiplication: Multiply each element of a matrix by a scalar.
Matlab
result = scalar * matrix;
Matrix multiplication: Multiply matrices according to linear algebra rules.
Matlab
result = matrix1 * matrix2;
Transpose: Create a new matrix by interchanging rows and columns.
Matlab
transpose_matrix = matrix';
Concatenation: Combine matrices horizontally or vertically using square brackets.
Matlab
horizontal_concat = [matrix1, matrix2];
vertical_concat = [matrix1; matrix2];
Determinant: Calculate the determinant of a square matrix using the det function.
Matlab
determinant = det(matrix);
Inverse: Calculate the inverse of a square matrix using the inv function.
Matlab
inverse_matrix = inv(matrix);
Special Matrices
Identity matrix: Create a square matrix with ones on the main diagonal using eye.
Matlab
identity_matrix = eye(3);
Zeros matrix: Create a matrix filled with zeros using zeros.
Matlab
zeros_matrix = zeros(2, 3);
Ones matrix: Create a matrix filled with ones using ones.
Matlab
ones_matrix = ones(3, 3);
Random matrix: Create a matrix with random elements using rand.
Matlab
random_matrix = rand(4, 4);
MATLAB Arrays: A Comprehensive Overview
Understanding Arrays
In MATLAB, all variables are fundamentally arrays, even scalars.
Vectors: One-dimensional arrays (row or column).
Matrices: Two-dimensional arrays.
Multidimensional arrays: Arrays with more than two dimensions.
Special Arrays
Zeros: Create an array filled with zeros using zeros(m,n).
Ones: Create an array filled with ones using ones(m,n).
Identity matrix: Create a square identity matrix using eye(n).
Random numbers: Create an array of random numbers between 0 and 1 using rand(m,n).
Magic square: Create a square matrix with equal row, column, and diagonal sums using magic(n).
Multidimensional Arrays
Creation: Extend a matrix by adding dimensions.
Matlab
A(:,:,2) = [1, 2, 3; 4, 5, 6; 7, 8, 9];
Use code with caution.
Accessing elements: Use multiple indices.
Matlab
element = A(2, 3, 1);
Use code with caution.
Function Purpose
transpose Transpose
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
MATLAB executes the statement and returns the following result −
ans =
100 95 90 85 80 75 70 65 60 55 50
Let us take another example −
0:pi/8:pi
MATLAB executes the statement and returns the following result −
ans =
Columns 1 through 7
0 0.3927 0.7854 1.1781 1.5708 1.9635 2.3562
Columns 8 through 9
2.7489 3.1416
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.
is a vector in four-dimensional array A. The vector includes A(i,j,k,1), A(i,j,k,2), A(i,j,k,3), and
A(i,j,k,:)
so on.
is all the elements of A, regarded as a single column. On the left side of an assignment
A(:) 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) % second column of A
A(:,2:3) % second and third column of A
A(2:3,2:3) % second and third rows and second and third columns
When you run the file, it displays the following result −
A=
1 2 3 4
4 5 6 7
7 8 9 10
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.
You can choose to store any number or array of numbers as integers or as single-precision numbers.
All numeric types support basic array operations and mathematical operations.
Conversion to Various Numeric Data Types
MATLAB provides the following functions to convert to various numeric data types −
Function Purpose
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 −
38 23 45
x=
38 23 45
x=
{
[1,1] = 38
[1,2] = 23
[1,3] = 45
}
Smallest and Largest Integers
The functions intmax() and intmin() return the maximum and minimum values that can be represented with
all types of integer numbers.
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 −
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
Functions for storing text in character arrays, combine character arrays, etc.
blanks Create string of blank characters
Functions for changing string to upper- or lowercase, creating or removing white space
A = pi*1000*ones(1,5);
sprintf(' %f \n %.2f \n %+.2f \n %12.2f \n %012.2f \n', A)
When you run the file, it displays the following result −
ans = 3141.592654
3141.59
+3141.59
3141.59
000003141.59
Joining Strings
Create a script file and type the following code into it −
disc;
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);
end % end of function quadratic2
You can call the above function from command prompt as −
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.
Private functions reside in subfolders with the special name private.
They are visible only to functions in the parent folder.
Example
Let us rewrite the quadratic function. This time, however, the disc function calculating the discriminant, will
be a private function.
Create a subfolder named private in working directory. Store the following function file disc.m in it −
function dis = disc(a,b,c)
%function calculates the discriminant
dis = sqrt(b^2 - 4*a*c);
end % end of sub-function
Create a function quadratic3.m in your working directory and type the following code in it −
function [x1,x2] = quadratic3(a,b,c)
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);
end % end of quadratic3
You can call the above function from command prompt as −
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 −
function avg = average(nums)
global TOTAL
avg = sum(nums)/TOTAL;
end
Create a script file and type the following code in it −
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
MATLAB - Data Import
Importing data in MATLAB means loading data from an external file. The importdata function allows loading
various data files of different formats. It has the following five forms −
Sr.No. Function & Description
A = importdata(filename)
1
Loads data into array A from the file denoted by filename.
A = importdata('-pastespecial')
2
Loads data from the system clipboard rather than from a file.
A = importdata(___, delimiterIn)
3 Interprets delimiterIn as the column separator in ASCII file, filename, or the clipboard data. You
can use delimiterIn with any of the input arguments in the above syntaxes.
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.
Our text file weeklydata.txt looks like this −
SunDay MonDay TuesDay WednesDay ThursDay FriDay SaturDay
95.01 76.21 61.54 40.57 55.79 70.28 81.53
73.11 45.65 79.19 93.55 75.29 69.87 74.68
60.68 41.85 92.18 91.69 81.32 90.38 74.51
48.60 82.14 73.82 41.03 0.99 67.22 93.18
89.13 44.47 57.63 89.36 13.89 19.88 46.60
Create a script file and type the following code in it −
filename = 'weeklydata.txt';
delimiterIn = ' ';
headerlinesIn = 1;
A = importdata(filename,delimiterIn,headerlinesIn);
% View data
for k = [1:7]
disp(A.colheaders{1, k})
disp(A.data(:, k))
disp(' ')
end
When you run the file, it displays the following result −
SunDay
95.0100
73.1100
60.6800
48.6000
89.1300
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.
Copy the following lines to the clipboard −
Mathematics is simple
Create a script file and type the following code −
A = importdata('-pastespecial')
When you run the file, it displays the following result −
A=
'Mathematics is simple'
Low-Level File I/O
The importdata function is a high-level function. The low-level file I/O functions in MATLAB allow the most
control over reading or writing data to a file. However, these functions need more detailed information about
your file to work efficiently.
MATLAB provides the following functions for read and write operations at the byte or character level −
Function Description
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 −
M = fscanf(fid, '%*s %*s\n%*s %*s %*s %*s\nM=%d\n\n', 1);
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.
Create a script file and type the following code in it −
filename = '/data/myfile.txt';
rows = 7;
cols = 5;
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
MATLAB - Data Output
Data export (or output) in MATLAB means to write into files. MATLAB allows you to use your data in another
application that reads ASCII files. For this, MATLAB provides several data export options.
You can create the following type of files −
Rectangular, delimited ASCII data file from an array.
Diary (or log) file of keystrokes and the resulting text output.
Specialized ASCII file using low-level functions such as fprintf.
MEX-file to access your C/C++ or Fortran routine that writes to a particular text file format.
Apart from this, you can also export data to spreadsheets.
There are two ways to export a numeric array as a delimited ASCII data file −
Using the save function and specifying the -ascii qualifier
Using the dlmwrite function
Syntax for using the save function is −
save my_data.out num_array -ascii
where, my_data.out is the delimited ASCII data file created, num_array is a numeric array and −ascii is the
specifier.
Syntax for using the dlmwrite function is −
dlmwrite('my_data.out', num_array, 'dlm_char')
where, my_data.out is the delimited ASCII data file created, num_array is a numeric array and dlm_char is
the delimiter character.
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];
save array_data1.out num_array -ascii;
type array_data1.out
dlmwrite('array_data2.out', num_array, ' ');
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
1234
4567
7890
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
Or export the cell array using low-level file I/O functions.
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.
For example, let us write the word 'hello' to a file −
h = 'hello';
save textdata.out h -ascii
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
Writing to Diary Files
Diary files are activity logs of your MATLAB session. The diary function creates an exact copy of your session
in a disk file, excluding graphics.
To turn on the diary function, type −
diary
Optionally, you can give the name of the log file, say −
diary logdata.out
To turn off the diary function −
diary off
You can open the diary file in a text editor.
Exporting Data to Text Data Files with Low-Level I/O
So far, we have exported numeric arrays. However, you may need to create other text files, including
combinations of numeric and character data, nonrectangular output files, or files with non-ASCII encoding
schemes. For these purposes, MATLAB provides the low-level fprintf function.
As in low-level I/O file activities, before exporting, you need to open or create a file with the fopen function
and get the file identifier. By default, fopen opens a file for read-only access. You should specify the
permission to write or append, such as 'w' or 'a'.
After processing the file, you need to close it with fclose(fid) function.
The following example demonstrates the concept −
Example
Create a script file and type the following code in it −
% Table Header
fprintf(fid, 'Log Function\n\n');
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
Define the function, y = f(x)
Call the plot command, as plot(x, y)
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.
Create a script file and type the following code −
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.
Create a script file and type the following code −
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 −
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
f(x) = 3x4 + 2x3+ 7x2 + 2x + 9 and
g(x) = 5x3 + 9x + 2
Create a script file and type the following code −
x = [-10 : 0.01: 10];
y = 3*x.^4 + 2 * x.^3 + 7 * x.^2 + 2 * x + 9;
g = 5 * x.^3 + 9 * x + 2;
plot(x, y, 'r', x, g, 'g')
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] )
The following example shows this −
Example
Create a script file and type the following code −
x = [0 : 0.01: 10];
y = exp(-x).* sin(2*x + 3);
plot(x, y), axis([0 10 -1 1])
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.
Syntax for the command is −
subplot(m, n, p)
where, m and n are the number of rows and columns of the plot array and p specifies 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)
Create a script file and type the following code −
x = [0:0.01:5];
y = exp(-1.5*x).*sin(10*x);
subplot(1,2,1)
plot(x,y), xlabel('x'),ylabel('exp(–1.5x)*sin(10x)'),axis([0 5 -1 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
Three dimensional plots
Drawing Bar Charts
The bar command draws a two dimensional bar chart. Let us take up an example to demonstrate the idea.
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.
Create a script file and type the following code −
x = [1:10];
y = [75, 58, 90, 87, 50, 85, 92, 75, 60, 95];
bar(x,y), xlabel('Student'),ylabel('Score'),
title('First Sem:')
print -deps graph.eps
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.
MATLAB provides a contour function for drawing contour maps.
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.
The meshgrid command is used for generating a matrix of elements that give the range over x and y along
with the specification of increment in each case.
Let us plot our function g = f(x, y), where −5 ≤ x ≤ 5, −3 ≤ y ≤ 3. Let us take an increment of 0.1 for both the
values. The variables are set as −
[x,y] = meshgrid(–5:0.1:5, –3:0.1:3);
Lastly, we need to assign the function. Let our function be: x2 + y2
Create a script file and type the following code −
[x,y] = meshgrid(-5:0.1:5,-3:0.1:3); %independent variables
g = x.^2 + y.^2; % our function
contour(x,y,g) % call the contour function
print -deps graph.eps
When you run the file, MATLAB displays the following contour map −
You can also use the mesh command to generate a three-dimensional surface. However, the surf command
displays both the connecting lines and the faces of the surface in color, whereas, the mesh command creates
a wireframe surface with colored lines connecting the defining points.
MATLAB - Algebra
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.
We will also discuss factorizing and simplification of algebraic expressions.
Solving Basic Algebraic Equations in MATLAB
The solve function is used for solving algebraic equations. In its simplest form, the solve function takes the
equation enclosed in quotes as an argument.
For example, let us solve for x in the equation x-5 = 0
solve('x-5=0')
MATLAB will execute the above statement and return the following result −
ans =
5
You can also call the solve function as −
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
If the equation involves multiple symbols, then MATLAB by default assumes that you are solving for x,
however, the solve function has another form −
solve(equation, variable)
where, you can also mention the 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
Solving Basic Algebraic Equations in Octave
The roots function is used for solving algebraic equations in Octave and you can write above examples as
follows −
For example, let us solve for x in the equation x-5 = 0
roots([1, -5])
Octave will execute the above statement and return the following result −
ans = 5
You can also call the solve function as −
y = roots([1, -5])
Octave will execute the above statement and return the following result −
y=5
Solving Quadratic Equations in MATLAB
The solve function can also solve higher order equations. It is often used to solve quadratic equations. The
function returns the roots of the equation in an array.
The following example solves the quadratic equation x2 -7x +12 = 0. Create a script file and type the following
code −
eq = 'x^2 -7*x + 12 = 0';
s = solve(eq);
disp('The first root is: '), disp(s(1));
disp('The second root is: '), disp(s(2));
When you run the file, it displays the following result −
The first root is:
3
The second root is:
4
Solving Quadratic Equations in Octave
The following example solves the quadratic equation x2 -7x +12 = 0 in Octave. Create a script file and type the
following code −
A = [5, 9; 3, -6];
b = [5;4];
A\b
When you run the file, it displays the following result −
ans =
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.
Create a script file and type the following code −
syms x %symbolic variable x
syms y %symbolic variable x
% 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 run the file, it displays the following result −
ans =
x^2 + 4*x - 45
ans =
x^4 + x^3 - 43*x^2 + 23*x + 210
ans =
2*cos(x)*sin(x)
ans =
cos(x)*cos(y) - sin(x)*sin(y)
ans =
x^4 - 7*x^3
ans =
x^6 - 8*x^5 + 15*x^4
Expanding and Collecting Equations in Octave
You need to have symbolic package, which provides expand and the collect function to expand and collect an
equation, respectively. The following example demonstrates the concepts −
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.
Create a script file and type the following code −
% first of all load the package, make sure its installed.
pkg load 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), 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)
Factorization and Simplification of Algebraic Expressions
The factor function factorizes an expression and the simplify function simplifies an expression. The following
example demonstrates the concept −
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))
When you run the file, it displays the following result −
ans =
(x - y)*(x^2 + x*y + y^2)
ans =
[ (x - y)*(x + y), (x + y)*(x^2 - x*y + y^2)]
ans =
x^2 + 4
MATLAB - Calculus
MATLAB provides various ways for solving problems of differential and integral calculus, solving differential
equations of any degree and calculation of limits. Best of all, you can easily plot the graphs of complex
functions and check maxima, minima and other stationery points on a graph by solving the original function,
as well as its derivative.
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.
In the next chapter Differential, we will compute derivative of an expression and find the local maxima and
minima on a graph. We will also discuss solving differential equations.
Finally, in the Integration chapter, we will discuss integral calculus.
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
limit((x^3 + 5)/(x^4 + 7))
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
Let's take another example,
limit(x^2 + 5, 3)
MATLAB will execute the above statement and return the following result −
ans =
14
Calculating Limits using Octave
Following is Octave version of the above example using symbolic package, try to execute and compare the
result −
pkg load symbolic
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
Verification of Basic Properties of Limits
Algebraic Limit Theorem provides some basic properties of limits. These are as follows −
l2 =
17
lAdd =
34
lSub =
0
lMult =
289
lDiv =
1
Verification of Basic Properties of Limits using Octave
Following is Octave version of the above example using symbolic package, try to execute and compare the
result −
pkg load symbolic
symbols
x = sym("x");
f = (3*x + 5)/(x-3);
g = x^2 + 1;
l1 = subs(f, x, 4)
l2 = subs (g, x, 4)
lAdd = subs (f+g, x, 4)
lSub = subs (f-g, x, 4)
lMult = subs (f*g, x, 4)
lDiv = subs (f/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
Left and Right Sided Limits
When a function has a discontinuity for some particular value of the variable, the limit does not exist at that
point. In other words, limits of a function f(x) has discontinuity at x = a, when the value of limit, as x
approaches x from left side, does not equal the value of the limit as x approaches from right side.
This leads to the concept of left-handed and right-handed limits. A left-handed limit is defined as the limit as
x -> a, from the left, i.e., x approaches a, for values of x < a. A right-handed limit is defined as the limit as x ->
a, from the right, i.e., x approaches a, for values of x > a. When the left-handed limit and right-handed limit
are not equal, the limit does not exist.
Let us consider a function −
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 −
By plotting the graph of the function and showing the discontinuity.
By computing the limits and showing that both are different.
The left-handed and right-handed limits are computed by passing the character strings 'left' and 'right' to the
limit command as the last argument.
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
MATLAB provides the diff command for computing symbolic derivatives. In its simplest form, you pass the
function you want to differentiate to diff command as an argument.
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
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
t = sym("t");
f = 3*t^2 + 2*t^(-2);
differentiate(f,t)
Octave executes the code and returns the following result −
ans =
-(4.0)*t^(-3.0)+(6.0)*t
Verification of Elementary Rules of Differentiation
Let us briefly state various equations or rules for differentiation of functions and verify these rules. For this
purpose, we will write f'(x) for a first order derivative and f"(x) for a second order derivative.
Following are the rules for differentiation −
Rule 1
For any functions f and g and any real numbers a and b are the derivative of the function −
h(x) = af(x) + bg(x) with respect to x is given by −
h'(x) = af'(x) + bg'(x)
Rule 2
The sum and subtraction rules state that if f and g are two functions, f' and g' are their derivatives
respectively, then,
(f + g)' = f' + g'
(f - g)' = f' - g'
Rule 3
The product rule states that if f and g are two functions, f' and g' are their derivatives respectively, then,
(f.g)' = f'.g + g'.f
Rule 4
The quotient rule states that if f and g are two functions, f' and g' are their derivatives respectively, then,
(f/g)' = (f'.g - g'.f)/g2
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)
f = (t^2 + 3)*(sqrt(t) + t^3)
der2 = diff(f)
f = (x^2 - 2*x + 1)*(3*x^3 - 5*x^2 + 2)
der3 = diff(f)
f = (2*x^2 + 3*x)/(x^3 + 1)
der4 = diff(f)
f = (x^2 + 1)^17
der5 = diff(f)
f = (t^3 + 3* t^2 + 5*t -9)^(-6)
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
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
x = sym("x");
t = sym("t");
f = (x + 2)*(x^2 + 3)
der1 = differentiate(f,x)
f = (t^2 + 3)*(t^(1/2) + t^3)
der2 = differentiate(f,t)
f = (x^2 - 2*x + 1)*(3*x^3 - 5*x^2 + 2)
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)
f = (t^3 + 3* t^2 + 5*t -9)^(-6)
der6 = differentiate(f,t)
Octave executes the code and returns the following result −
f=
(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
ca.x ca.x.ln c.a (ln is natural logarithm)
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
Following is Octave equivalent of the above calculation −
pkg load symbolic
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 = 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)
Octave executes the code and returns the following result −
y=
exp(x)
ans =
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)
Computing Higher Order Derivatives
To compute higher derivatives of a function f, we use the syntax diff(f,n).
Let us compute the second derivative of the function y = f(x) = x .e-3x
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)
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
x = sym("x");
f = x*Exp(-3*x);
differentiate(f, x, 2)
Octave executes the code and returns the following result −
ans =
(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.
Create a script file and type the following code into it −
syms x
y = 3*sin(x)+7*cos(5*x); % defining the function
lhs = diff(y,2)+y; %evaluting the lhs of the equation
rhs = -5*cos(2*x); %rhs of the equation
if(isequal(lhs,rhs))
disp('Yes, the equation holds true');
else
disp('No, the equation does not hold true');
end
disp('Value of LHS is: '), disp(lhs);
When you run the file, it displays the following result −
No, the equation does not hold true
Value of LHS is:
-168*cos(5*x)
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
x = sym("x");
y = 3*Sin(x)+7*Cos(5*x); % defining the function
lhs = differentiate(y, x, 2) + y; %evaluting the lhs of the equation
rhs = -5*Cos(2*x); %rhs of the equation
if(lhs == rhs)
disp('Yes, the equation holds true');
else
disp('No, the equation does not hold true');
end
disp('Value of LHS is: '), disp(lhs);
Octave executes the code and returns the following result −
No, the equation does not hold true
Value of LHS is:
-(168.0)*cos((5.0)*x)
Finding the Maxima and Minima of a Curve
If we are searching for the local maxima and minima for a graph, we are basically looking for the highest or
lowest points on the graph of the function at a particular locality, or for a particular range of values of the
symbolic variable.
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.
To find the stationary points of a function we differentiate, we need to set the derivative equal to zero and
solve the equation.
Example
Let us find the stationary points of the function f(x) = 2x3 + 3x2 − 12x + 17
Take the following steps −
First let us enter the function and plot its graph.
syms x
y = 2*x^3 + 3*x^2 - 12*x + 17; % defining the function
ezplot(y)
MATLAB executes the code and returns the following plot −
ezplot(y)
print -deps graph.eps
Our aim is to find some local maxima and minima on the graph, so let us find the local maxima and minima
for the interval [-2, 2] on the graph.
syms x
y = 2*x^3 + 3*x^2 - 12*x + 17; % defining the function
ezplot(y, [-2, 2])
MATLAB executes the code and returns the following plot −
x = sym('x');
y = inline("2*x^3 + 3*x^2 - 12*x + 17");
x = sym("x");
y = 2*x^3 + 3*x^2 - 12*x + 17;
g = differentiate(y,x)
Octave executes the code and returns the following result −
g=
-12.0+(6.0)*x+(6.0)*x^(2.0)
Let us solve the derivative function, g, to get the values where it becomes zero.
s = solve(g)
MATLAB executes the code and returns the following result −
s=
1
-2
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
x = sym("x");
y = 2*x^3 + 3*x^2 - 12*x + 17;
g = differentiate(y,x)
roots([6, 6, -12])
Octave executes the code and returns the following result −
g=
-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.
subs(y, 1), subs(y, -2)
MATLAB executes the code and returns the following result −
ans =
10
ans =
37
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
x = sym("x");
y = 2*x^3 + 3*x^2 - 12*x + 17;
g = differentiate(y,x)
roots([6, 6, -12])
subs(y, x, 1), subs(y, x, -2)
ans =
10.0
ans =
37.0-4.6734207789940138748E-18*I
Therefore, The minimum and maximum values on the function f(x) = 2x3 + 3x2 − 12x + 17, in the interval [-2,2]
are 10 and 37.
Solving Differential Equations
MATLAB provides the dsolve command for solving differential equations symbolically.
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.
It returns a symbolic solution with a set of arbitrary constants that MATLAB labels C1, C2, and so on.
You can also specify initial and boundary conditions for the problem, as comma-delimited list following the
equation as −
dsolve('eqn','cond1', 'cond2',…)
For the purpose of using dsolve command, derivatives are indicated with a D. For example, an equation like
f'(t) = -2*f + cost(t) is entered as −
'Df = -2*f + cos(t)'
Higher derivatives are indicated by following D by the order of the derivative.
For example the equation f"(x) + 2f'(x) = 5sin3x should be entered as −
'D2y + 2Dy = 5*sin(3*x)'
Let us take up a simple example of a first order differential equation: y' = 5y.
s = dsolve('Dy = 5*y')
MATLAB executes the code and returns the following result −
s=
C2*exp(5*t)
Let us take up another example of a second order differential equation as: y" - y = 0, y(0) = -1, y'(0) = 2.
dsolve('D2y - y = 0','y(0) = -1','Dy(0) = 2')
MATLAB executes the code and returns the following result −
ans =
exp(t)/2 - (3*exp(-t))/2
MATLAB - Integration
Integration deals with two essentially different types of problems.
In the first type, derivative of a function is given and we want to find the function. Therefore, we
basically reverse the process of differentiation. This reverse process is known as anti-differentiation,
or finding the primitive function, or finding an indefinite integral.
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.
Finding Indefinite Integral Using MATLAB
By definition, if the derivative of a function f(x) is f'(x), then we say that an indefinite integral of f'(x) with
respect to x is f(x). For example, since the derivative (with respect to x) of x2 is 2x, we can say that an
indefinite integral of 2x is x2.
In symbols −
f'(x2) = 2x, therefore,
∫ 2xdx = x2.
Indefinite integral is not unique, because derivative of x2 + c, for any value of a constant c, will also be 2x.
This is expressed in symbols as −
∫ 2xdx = x2 + c.
Where, c is called an 'arbitrary constant'.
MATLAB provides an int command for calculating integral of an expression. To derive an expression for the
indefinite integral of a function, we write −
int(f);
For example, from our previous example −
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)
When you run the file, it displays the following result −
ans =
piecewise([n == -1, log(x)], [n ~= -1, x^(n + 1)/(n + 1)])
f=
sin(n*t)
ans =
-cos(n*t)/n
ans =
(a*sin(pi*t))/pi
ans =
a^x/log(a)
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)
pretty(int(1 - 10*x + 9 * x^2))
int((3 + 5*x -6*x^2 - 7*x^3)/2*x^2)
pretty(int((3 + 5*x -6*x^2 - 7*x^3)/2*x^2))
Note that the pretty function returns an expression in a more readable format.
When you run the file, it displays the following result −
ans =
sin(x)
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
7x 3x 5x x
- ---- - ---- + ---- + --
12 5 8 2
we write,
int(x, a, b)
For example, to calculate the value of we write −
int(x, 4, 9)
MATLAB executes the above statement and returns the following result −
ans =
65/2
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
x = sym("x");
f = x;
c = [1, 0];
integral = polyint(c);
32.500
An alternative solution can be given using quad() function provided by Octave as follows −
pkg load symbolic
symbols
f = inline("x");
[a, ierror, nfneval] = quad(f, 4, 9);
x = sym("x");
f = x^3 - 2*x +5;
c = [1, 0, -2, 5];
integral = polyint(c);
5.7500
An alternative solution can be given using quad() function provided by Octave as follows −
pkg load symbolic
symbols
x = sym("x");
f = inline("x^3 - 2*x +5");
Area:
0.3326
Following is Octave equivalent of the above calculation −
pkg load symbolic
symbols
x = sym("x");
f = inline("x^2*cos(x)");
ezplot(f, [-4,9])
print -deps graph.eps
p = [1 7 0 -5 9];
polyval(p,4)
MATLAB executes the above statements and returns the following result −
ans = 693
MATLAB also provides the polyvalm function for evaluating a matrix polynomial. A matrix polynomial is
a polynomial with matrices as variables.
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
Finding the Roots of Polynomials
The roots function calculates the roots of a polynomial. For example, to calculate the roots of our polynomial
p, type −
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:
MATLAB - Transforms
MATLAB provides command for working with transforms, such as the Laplace and Fourier transforms.
Transforms are used in science and engineering as a tool for simplifying analysis and look at data from
another angle.
For example, the Fourier transform allows us to convert a signal represented as a function of time to a
function of frequency. Laplace transform allows us to convert a differential equation to an algebraic equation.
MATLAB provides the laplace, fourier and fft commands to work with Laplace, Fourier and Fast Fourier
transforms.
The Laplace Transform
The Laplace transform of a function of time f(t) is given by the following integral −
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 transform turns differential equations into algebraic ones. To compute a Laplace transform of a
function f(t), write −
laplace(f(t))
Example
In this example, we will compute the Laplace transform of some commonly used functions.
Create a script file and type the following code −
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))
When you run the file, it displays the following result −
ans =
1/s^2
ans =
2/s^3
ans =
362880/s^10
ans =
1/(b + s)
ans =
w/(s^2 + w^2)
ans =
s/(s^2 + w^2)
The Inverse Laplace Transform
MATLAB allows us to compute the inverse Laplace transform using the command ilaplace.
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))
When you run the file, it displays the following result −
ans =
t^6/720
ans =
2*exp(-t*w)
ans =
cos(2*t)
ans =
ilaplace(exp(-b*t), t, x)
ans =
sin(t*w)
ans =
cos(t*w)
The Fourier Transforms
Fourier transforms commonly transforms a mathematical function of time, f(t), into a new function,
sometimes denoted by or F, whose argument is frequency with units of cycles/s (hertz) or radians per
second. The new function is then known as the Fourier transform and/or the frequency spectrum of the
function f.
Example
Create a script file and type the following code in it −
syms x
f = exp(-2*x^2); %our function
ezplot(f,[-2,2]) % plot of our function
FT = fourier(f) % Fourier transform
When you run the file, MATLAB plots the following graph −
Area:
0.3326
But to give area of the same curve in Octave, you will have to make use of symbolic package as follows −
pkg load symbolic
symbols
x = sym("x");
f = inline("x^2*cos(x)");
ezplot(f, [-4,9])
print -deps graph.eps
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.
Let us drag a 'Sine Wave' block into the model.
Select 'Sinks' from the library and drag a 'Scope' block into the model.
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)
You should get the below graph from the scope.