Introduzione a

MatLab

Prof. Sebastiano Battiato

MatLab Environment
 MATLAB

Matlab = Matrix Laboratory

Originally a user interface for numerical linear algebra
routines (Lapak/Linpak)
Commercialized 1984 by The Mathworks
Since then heavily extended (defacto-standard)

Alternatives Complements
Matrix-X Maple (symbolic)
Octave(free; GNU) Mathematica (symbolic)
Lyme (free; Palm)
 MATLAB

z The MATLAB environment is command oriented somewhat

like UNIX. A prompt appears on the screen and a
MATLAB statement can be entered. When the <ENTER>
key is pressed, the statement is executed, and another
prompt appears.

z If a statement is terminated with a semicolon ( ; ), no

results will be displayed. Otherwise results will appear
before the next prompt.

z The following slide is the text from a MATLAB screen.

 MATLAB

To get started, type one of these commands: helpwin,

helpdesk, or demo

» a=5;
» b=a/2

b=

2.5000

»
 MATLAB Variable Names

z Variable names ARE case sensitive

z Variable names can contain up to 63 characters (as of

MATLAB 6.5 and newer)

z Variable names must start with a letter followed by letters,

digits, and underscores.
MATLAB Special Variables

ans Default variable name for results

pi Value of π
inf ∞
NaN Not a number e.g. 0/0
i (and j) i=j= −1
eps Smallest incremental number
realmin The smallest usable positive real number
realmax The largest usable positive real number
MATLAB Math & Assignment
Operators

Power ^ or .^ a^b or a.^b

Multiplication * or .* a*b or a.*b
Division / or ./ a/b or a./b
or \ or .\ b\a or b.\a
NOTE: 56/8 = 8\56

- (unary) + (unary)
Addition + a + b
Subtraction - a - b
Assignment = a = b (assign b to a)
 Other MATLAB symbols

>> prompt
... continue statement on next line
, separate statements and data
% start comment which ends at end of line
; (1) suppress output
(2) used as a row separator in a matrix
: specify range
 MATLAB Help System

Search for appropriate function

>> lookfor keyword

Rapid help with syntax and function definition

>> help function

An advanced hyperlinked help system is launched by

>> helpdesk

Complete manuals as PDF files

 MATLAB Matrices

z MATLAB treats all variables as matrices. For our purposes

a matrix can be thought of as an array, in fact, that is how it
is stored.

z Vectors are special forms of matrices and contain only one

row OR one column.

z Scalars are matrices with only one row AND one column
 MATLAB Matrices

z A matrix with only one row AND one column is a scalar. A

scalar can be created in MATLAB as follows:

» a_value=23

a_value =

23
 MATLAB Matrices

z A matrix with only one row is called a row vector. A row

vector can be created in MATLAB as follows (note the
commas):

» rowvec = [12 , 14 , 63]

rowvec =

12 14 63
 MATLAB Matrices

z A matrix with only one column is called a column vector. A

column vector can be created in MATLAB as follows (note
the semicolons):

» colvec = [13 ; 45 ; -2]

colvec =

13
45
-2
 MATLAB Matrices

z A matrix can be created in MATLAB as follows (note the

commas AND semicolons):

» matrix = [1 , 2 , 3 ; 4 , 5 ,6 ; 7 , 8 , 9]

matrix =

1 2 3
4 5 6
7 8 9
 Extracting a Sub-Matrix

z A portion of a matrix can be extracted and stored in a

smaller matrix by specifying the names of both matrices
and the rows and columns to extract. The syntax is:

sub_matrix = matrix ( r1 : r2 , c1 : c2 ) ;

where r1 and r2 specify the beginning and ending rows and

c1 and c2 specify the beginning and ending columns to be
extracted to make the new matrix.
 MATLAB Matrices

A column vector can be Here we extract column 2 of

extracted from a matrix. As the matrix and make a
an example we create a column vector:
matrix below:

» matrix=[1,2,3;4,5,6;7,8,9] » col_two=matrix( : , 2)

matrix = col_two =

1 2 3 2
4 5 6 5
7 8 9 8
 MATLAB Matrices

A row vector can be Here we extract row 2 of the

extracted from a matrix. As matrix and make a row
an example we create a vector. Note that the 2:2
matrix below: specifies the second row and
the 1:3 specifies which
» matrix=[1,2,3;4,5,6;7,8,9] columns of the row.

matrix = » rowvec=matrix(2 : 2 , 1 : 3)

1 2 3 rowvec =
4 5 6
7 8 9 4 5 6
 Reading Data from files

z MATLAB supports reading an entire file and creating a

matrix of the data with one statement.
>> load mydata.dat; % loads file into matrix.
% The matrix may be a scalar, a vector, or a
% matrix with multiple rows and columns. The
% matrix will be named mydata.
>> size (mydata) % size will return the number
% of rows and number of
% columns in the matrix
>> length (myvector) %length will return the total
% no. of elements in myvector
 Plotting with MATLAB

z MATLAB will plot one vector vs. another. The first one will
be treated as the abscissa (or x) vector and the second as
the ordinate (or y) vector. The vectors have to be the same
length.
z MATLAB will also plot a vector vs. its own index. The index
will be treated as the abscissa vector. Given a vector “time”
and a vector “dist” we could say:
>> plot (time, dist) %plotting versus time
>> plot (dist) %plotting versus index
 Plotting with MATLAB

z There are commands in MATLAB to "annotate" a plot to put

on axis labels, titles, and legends. For example:
>> % To put a label on the axes we would use:
>> xlabel ('X-axis label')
>> ylabel ('Y-axis label')

>> % To put a title on the plot, we woulduse:

>> title ('Title of my plot')
 Plotting with MATLAB

z Vectors may be extracted from matrices. Normally, we

wish to plot one column vs. another. If we have a matrix
“mydata” with two columns, we can obtain the columns as
a vectors with the assignments as follows:

>> first_vector=mydata(:,1); % First column

>> second_vector=mydata(:,2); % Second one
>> % and we can plot the data
>> plot ( first_vector , second_vector )
 Some Useful MATLAB
commands

who List known variables

whos List known variables plus their size
help Ex: >> help sqrt Help on using sqrt
lookfor Ex: >> lookfor sqrt Search for
keyword sqrt in m-files
what Ex: >> what a: List MATLAB files in a:
clear Clear all variables from work space
clear x y Clear variables x and y from work space
clc Clear the command window
 Some Useful MATLAB
commands

what List all m-files in current directory

dir List all files in current directory
ls Same as dir
type test Display test.m in command window
delete test Delete test.m
cd a: Change directory to a:
chdir a: Same as cd
pwd Show current directory
which test Display current directory path to test.m
MATLAB Relational Operators

z MATLAB supports six relational operators.

Less Than <

Less Than or Equal <=
Greater Than >
Greater Than or Equal >=
Equal To ==
Not Equal To ~=
 MATLAB Logical Operators

z MATLAB supports three logical operators.

not ~ % highest precedence

and & % equal precedence with or
or | % equal precedence with and
 MATLAB Logical
Functions
MATLAB also supports some logical functions.
xor (exclusive or) Ex: xor (a, b)
Where a and b are logical expressions. The xor operator
evaluates to true if and only if one expression is true and the
other is false. True is returned as 1, false as 0.

any (x) returns 1 if any element of x is nonzero

all (x) returns 1 if all elements of x are nonzero
isnan (x) returns 1 at each NaN in x
isinf (x) returns 1 at each infinity in x
finite (x) returns 1 at each finite value in x
 Matlab Selection Structures

z An if - elseif - else structure in MATLAB.

Note that elseif is one word.

if expression1 % is true
% execute these commands
elseif expression2 % is true
% execute these commands
else % the default
% execute these commands
end
MATLAB Repetition Structures

z A for loop in MATLAB for x = array

for x = 1: 0.5 : 10
% execute these commands
end

z A while loop in MATLAB while expression

while x <= 10
% execute these commands
end
 Scalar - Matrix Addition

» a=3;
» b=[1, 2, 3;4, 5, 6]
b=
1 2 3
4 5 6
» c= b+a % Add a to each element of b
c=
4 5 6
7 8 9
Scalar - Matrix Subtraction

» a=3;
» b=[1, 2, 3;4, 5, 6]
b=
1 2 3
4 5 6
» c = b - a %Subtract a from each element of b
c=
-2 -1 0
1 2 3
 Scalar - Matrix Multiplication

» a=3;
» b=[1, 2, 3; 4, 5, 6]
b=
1 2 3
4 5 6
» c = a * b % Multiply each element of b by a
c=
3 6 9
12 15 18
 Scalar - Matrix Division

» a=3;
» b=[1, 2, 3; 4, 5, 6]
b=
1 2 3
4 5 6
»c=b/a % Divide each element of b by a
c=
0.3333 0.6667 1.0000
1.3333 1.6667 2.0000
 MATLAB Toolbox

z MATLAB has a number of add-on software modules,

called toolbox, that perform more specialized
computations.

More than 60 toolboxes!

Image Processing & Matlab
 >> clear
>> imread('citta.jpg');
>> subplot(2,2,1);
>> f=imread('citta.jpg');
>> imshow(f);
>> flipped=f(end:-1:1, :, :);
>> subplot(2,2,2);
300
>> imshow(flipped);
250

200
>> cropped=f(60:300, 60:420,:);
150
>> subplot(2,2,3);
100 >> imshow(cropped);
50 >> subplot(2,2,4);
0
0 100 200 300 400 >> plot(f(:,10,1));
 Tips and Tricks

z Code Optimization
– Preallocating Arrays: zeros(M,N)
– Vectorizing Loops (1D):

The matlab code:

x=0:M-1
f=A*sin(x/2*pi));
Is better/faster than
for x=1:M
f(x)=A*sin((x-1/(2*pi))
end;
 Tips and Tricks

z Code Optimization
– Vectorizing Loops (2D):

It is possible to use the function meshgrid,which has the

syntax:
[C,R]=meshgrid(c,r)

The function transforms the domain specified by row

vectors c and r into arrays C and R that can be used for
the evaluation of functions of two variables. The rows of
the output array C are copies of the vector c and the
columns of the output array R are copies of the vector r.
 Tips and Tricks
>> clear h =
>> c=[0,1]; % y = 0,1
0 1
>> r=[0,1,2]; % x = 0,1,2
1 2
>> [C,R]=meshgrid(c,r);
>> C 4 5
C =
0 1 >>
0 1
0 1
>> R
R =
0 0
1 1
2 2
>> h=R.^2+C.^2
 Example 1

Compute: f(x,y)=A sin(u0x + v0y)

function [ratio_time,f,g]=twodsinnew(A,u0,v0,M,N)
tic %Start timing
for r=1:M
u0_x=u0*(r-1);
for c=1:N
v0_y=v0*(c-1);
f(r,c)=A*sin(u0_x + v0_y);
end;
end;
t1=toc; %End Timing
 Example 1

tic
r=0:M-1;
c=0:N-1;
[C,R]=meshgrid(c,r);
g=A*sin(u0*R+v0*C);
t2=toc; %End timing

ratio_time= t1/(t2+eps) % Eps in case of t2=0

 Example 2: Color Transform
The human eye is more sensitive to luminance than to chrominance. Typically JPEG throw out 3/4 of
the chrominance information before any other compression takes place. This reduces the amount of
information to be stored about the image by 1/2. With all three components fully stored, 4 pixels needs
3 x 4 = 12 component values. If 3/4 of two components are discarded we need 1 x 4 + 2 x 1 = 6 values.

B Y Cr

G
Cb
R

Example
RGB vs YCbCr
 RGB conversion &
subsampling
RGB → CMY
⎡ C ⎤ ⎡1⎤ ⎡ R ⎤
⎢ M ⎥ = ⎢1⎥ − ⎢ G ⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢⎣ Y ⎥⎦ ⎢⎣1⎥⎦ ⎢⎣ B ⎥⎦
RGB → YCbCr
⎡ Y ⎤ ⎡0.299 0.587 0.114 ⎤ ⎡ R ⎤
⎢C ⎥ = ⎢0.596 − 0.275 − 0.321⎥ ⎢G⎥
⎢ b⎥ ⎢ ⎥⎢ ⎥
⎢⎣Cr ⎥⎦ ⎢⎣0.212 − 0.523 0.311 ⎥⎦ ⎢⎣ B ⎥⎦

Subsampling
4:4:4 (no subsampling)
4:2:2 (Cb, Cr horizontal subsampling)
4:2:0 (Cb, Cr horizontal + vertical subsampling)
 Tutorials on Web

www.imageprocessingplace.com/DIPUM/dipum_tutorials/tutorials.htm

See also: Chapters 1-2 Digital Image Processing using Matlab