Scribd
Upload
83 views37 pages

An Introduction To Matlab: Kocaeli University Assist. Prof. Dr. Osman Büyük Digital Signal Processing Applications

This document provides an introduction to Matlab. It discusses what Matlab is, its main features and interfaces. It covers topics like variables, arrays, matrices, basic math operations, plotting functions, and control structures. It also demonstrates how to write user-defined functions in Matlab using M-files. The document aims to give readers an overview of Matlab's capabilities and help them get started with the fundamentals of using Matlab for tasks like signal processing.

Uploaded by

Md.Arifur Rahman
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
83 views37 pages

An Introduction To Matlab: Kocaeli University Assist. Prof. Dr. Osman Büyük Digital Signal Processing Applications

This document provides an introduction to Matlab. It discusses what Matlab is, its main features and interfaces. It covers topics like variables, arrays, matrices, basic math operations, plotting functions, and control structures. It also demonstrates how to write user-defined functions in Matlab using M-files. The document aims to give readers an overview of Matlab's capabilities and help them get started with the fundamentals of using Matlab for tasks like signal processing.

Uploaded by

Md.Arifur Rahman
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 37

An Introduction to Matlab

Kocaeli University
Assist. Prof. Dr. Osman Byk
Digital Signal Processing Applications

What is Matlab?





Matlab = Matrix Laboratory


Easy to learn
Robust
Powerful functions and toolboxes






fft() : fast fourier transform


plot()
wavread() : For reading speech files
lpc() : Linear predictive coding
hist() and many others ...

What is Matlab?


Overview
Matlab
Series of
Matlab
commands
m-files

functions
Input
Output
capability

Command
Line
Command execution
like DOS command
window

mat-files

Data
storage/
loading

Features


Command Window

type commands

Current Directory

View folders and m-files

Workspace

View program variables

Double click on a variable
to see it in the Array Editor

Command History

view past commands

save a whole session
using diary

Variables


No need for types. i.e.,


int a;
double b;
float c;

All variables are created with double precision unless


specified and they are matrices.
Example:
>>x=5;
>>x1=2;

After these statements, the variables are 1x1 matrices


with double precision

Array, Matrix


A vector

x = [1 2 5 1]
x =
1

Transpose

y = x

y =
1
2
5
1

A matrix : Rows are seperated by semicolumn


x = [1 2 3; 5 1 4; 3 2 -1]
x =
1
5
3

2
1
2

3
4
-1

Long Array, Matrix


t =1:10

t =

k =2:-0.5:-1

k =

1.5

0.5

= [1:4; 5:8]

x =
1
5

2
6

3
7

4
8

-0.5

-1

10

Generating Vectors from functions




zeros(M,N) MxN matrix of 0s

ones(M,N)

MxN matrix of 1s

rand(M,N)

MxN matrix of uniformly


distributed random
numbers on (0,1)

x = zeros(1,3)
x =
0
0
0
x = ones(1,3)
x =
1
1
1
x = rand(1,3)
x =
0.9501 0.2311 0.6068

Matrix Index



The matrix indices begin from 1 (not 0 as in C)


The matrix indices must be positive integer

Given:

A(-2), A(0)
Error: ??? Subscript indices must either be real positive integers or logicals.
A(4,2)
Error: ??? Index exceeds matrix dimensions.

Exercise


Print the element in row 3 and column 2

Print all elements in row 2

Print first and second elements in column 2

Concatenation of Matrices


x = [1 2], y = [4 5], z=[0 0]


A = [x y]
1

B = [x ; y]
1 2
4 5
C = [x y ;z]
Error:
??? Error using ==> vertcat CAT arguments dimensions are not consistent.

Operators (arithmetic)
+
*
/
^

addition
subtraction
multiplication
division
power
complex conjugate transpose

Matrices Operations

Given A and B:

Addition

Subtraction

Product

Transpose

Operators (Element by Element)

.* element-by-element multiplication
./ element-by-element division
.^ element-by-element power

The use of . Element Operation


A = [1 2 3; 5 1 4; 3 4 -1]
A=
1 2 3
5 1 4
3 4 -1

x = A(1,:)
x=

c=x./y

d = x .^2

b=

c=
0.33 0.5 -3

d=

y = A(3 ,:)
y=

1 2 3

b = x .* y

3 8 -3
3 4 -1

K= x^2
Erorr:
??? Error using ==> mpower Matrix must be square.
B=x*y
Erorr:
??? Error using ==> mtimes Inner matrix dimensions must agree.

Exercise



x=[1 2;3 4];


z=[1 2 3];
1)
2)
3)
4)
5)

y=[5 6;7 8];


t=[4 5 6];

Sum of x and y
Product of x and y
Element by element product of x and y
Element by element division of z over t
Element by element division of z over t, but for
the first two elements only

Basic Task: Plot the function sin(x)


between 0x4


Create an x-array of 100 samples between 0


and 4. May use linspace function.
>>x=linspace(0,4*pi,100);

Calculate sin(.) of the x-array


1
0.8

>>y=sin(x);

0.6
0.4
0.2
0

Plot the y-array

-0.2
-0.4
-0.6

>>plot(y)

-0.8
-1

10

20

30

40

50

60

70

80

90

100

Plot the function e-x/3sin(x) between


0x4


Create an x-array of 100 samples between 0


and 4.
>>x=linspace(0,4*pi,100);

Calculate sin(.) of the x-array


>>y=sin(x);

Calculate e-x/3 of the x-array


>>y1=exp(-x/3);

Multiply the arrays y and y1


>>y2=y*y1;

Plot the function e-x/3sin(x) between


0x4


Multiply the arrays y and y1 correctly


>>y2=y.*y1;

Plot the y2-array


0.7

>>plot(y2)

0.6
0.5
0.4
0.3
0.2
0.1
0
-0.1
-0.2
-0.3

10

20

30

40

50

60

70

80

90

100

Display Facilities
0.7
0.6

plot(.)

0.5

x=linspace(0,4*pi,100);
y=sin(x);
y1=exp(-x/3);
y2=y.*y1;
plot(x,y2)

0.4
0.3
0.2
0.1
0
-0.1
-0.2

stem(.)

-0.3

stem(y2)
stem(x,y2)

10

20

30

40

50

60

70

80

90

100

0.7
0.6
0.5

subplot(.)

subplot(3,1,1)
plot(x,y);
subplot(3,1,2)
plot(x,y1)
subplot(3,1,3)
plot(x,y2)

0.4
0.3
0.2
0.1
0
-0.1
-0.2
-0.3

10

20

30

40

50

60

70

80

90

100

Display Facilities


title(.)
This is the sinus function

>>title(This is the sinus function)

1
0.8

xlabel(.)

0.6
0.4

>>xlabel(x (secs))
sin(x)

0.2

ylabel(.)

0
-0.2
-0.4
-0.6

>>ylabel(sin(x))

legend(.)
>>legend(Sinus function)

-0.8
-1

10

20

30

40

50
60
x (secs)

70

80

90

100

Exercise


Write a program which plots 2


waves in the same graph.














Two cycles of sin wave


Two cycles of cos wave.

t=0:0.1:4*pi;
x=sin(t);
y=cos(t);
plot(t,x,'r+-',t,y,'bo-.');
figure;
plot(t,x,'r+-');
hold on;
plot(t,y,'bo-.)
legend('sin','cos');
hold off;

1
0.8
0.6
0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
-1

10

12

14

Operators (relational, logical)











== Equal to
~= Not equal to
< Strictly smaller
> Strictly greater
<= Smaller than or equal to
>= Greater than equal to
& And operator
| Or operator

Flow Control








if
for
while
break
continue
return
.

Control Structures


If Statement Syntax

if (Condition_1)
Matlab Commands
elseif (Condition_2)
Matlab Commands
elseif (Condition_3)
Matlab Commands
else
Matlab Commands
end

Some Dummy Examples


if ((a>3) & (b==5))
Some Matlab Commands;
end
if (a<3)
Some Matlab Commands;
elseif (b~=5)
Some Matlab Commands;
end
if (a<3)
Some Matlab Commands;
else
Some Matlab Commands;
end

Control Structures
Some Dummy Examples


For loop syntax

for i=Index_Array
Matlab Commands
end

for i=1:100
Some Matlab Commands;
end
for j=1:3:200
Some Matlab Commands;
end
for m=13:-0.2:-21
Some Matlab Commands;
end
for k=[0.1 0.3 -13 12 7 -9.3]
Some Matlab Commands;
end

Control Structures


While Loop Syntax

while (condition)
Matlab Commands
end

Dummy Example
while ((a>3) & (b==5))
Some Matlab Commands;
end

Exercise
Write a for program to compute 0.5+1+1.5+2++50


1.
2.
3.

4.

With a for loop


With a while loop
Without loop : Might use sum(.). If possible avoid loops in
Matlab.
Without loop, without sum(.)

a=[0.5:0.5:50];
x=sum(a)
a=[0.5:0.5:50];
b=ones(length(a),1);
x=a*b

Use of M-File
Click to create
a new M-File

Extension .m
A text file containing script or function or program to run

Use of M-File

Save file as Denem430.m

If you include ; at the


end of each statement,
result will not be shown
immediately

Writing User Defined Functions




Functions are m-files which can be executed by


specifying some inputs and supply some desired outputs.
The code telling the Matlab that an m-file is actually a
function is
function out1=functionname(in1)
function out1=functionname(in1,in2,in3)
function [out1,out2]=functionname(in1,in2)

You should write this command at the beginning of the


m-file and you should save the m-file with a file name
same as the function name

Writing User Defined Functions




Example
 Write a function : out=squarer (A, ind)
 Takes the square of the input matrix, A, if the input
indicator, ind, is equal to 1
 Takes the element by element square of A if ind is equal
to 2
Same Name

Useful Commands


% is the neglect sign for Matlab

It is equivalent to \\ in C

Anything after it on the same line is neglected by Matlab


compiler

Add your comments to your functions with %

It is a good programming practice to add comments to


your codes

Usefull commands


Sometimes slowing down the execution is


done deliberately for observation purposes
You can use the command pause for this
purpose
pause %wait until any key
pause(3) %wait 3 seconds

Useful Commands


The following two commands are very


frequently used
When you need some help, ask to Matlab
>>help functionname

>>lookfor keyword

Exercise


Letss implement dft(.) in Matlab





With loops
Without loops


In Matlab, try to avoid loops when possible since they


are slow

Questions?

Thank you!

You might also like