MATLAB

Tutorial Course

1
 Contents
1. Session One
 What is Matlab?
 MATLAB Parts
 MATLAB Desktop
 Matrices
 Numerical Arrays
 String Arrays
 Elementary Math
 Logical Operators
 Math Functions
 Polynomials and Interpolation
2  Importing and Exporting Data
 Contents Continued
 Graphics Fundamentals
 2D plotting
 Subplots
 3D plotting
 Specialized Plotting
 Editing and Debugging M-files
2. Session Two
 Script and Function Files
 Basic Parts of an M-file
 Flow Control Statements
 M-file Programming
3
 Contents Continued
 Data types
 Multidimensional Arrays
 Structures
 Cell Arrays
 Nonlinear Numerical Functions
 Ordinary Differential Equations (ODE)
 Handle Graphics
 Graphic Objects
 Graphical User Interface (GUI)

4
 What is MATLAB?
high-performance software
 Computation
 Visualization
 Easy-to-use environment.
high-level language
 Data types
 Functions
 Control flow statements
 Input/output
 Graphics
 Object-oriented programming capabilities
5
 MATLAB Parts

Developed Environment
Programming Language
Graphics
Toolboxes
Application Program Interface

6
 Toolboxes
Collections of functions to solve problems of
several applications.
 DSP Toolbox
 Image Toolbox
 Wavelet Toolbox
 Neural Network Toolbox
 Fuzzy Logic Toolbox
 Control Toolbox
 Communication Toolbox

7
 MATLAB Desktop Tools
Command Window
Command History
Help Browser
Workspace Browser
Editor/Debugger
Launch Pad
8
 Calculations at the Command
Line
MATLAB as a calculator Assigning Variables
»» -5/(4.8+5.32)^2
-5/(4.8+5.32)^2 »» aa == 2; Semicolon
2;
ans
ans == suppresses
»» bb == 5;
5;
-0.0488
-0.0488 screen output
»» (3+4i)*(3-4i) »» a^b
a^b
(3+4i)*(3-4i)
ans ans
ans ==
ans == Results
25
25 32
32 assigned to
»» cos(pi/2)
cos(pi/2) »» xx == 5/2*pi; “ans” if name
5/2*pi; not specified
ans
ans == »» yy == sin(x)
6.1230e-017 sin(x)
6.1230e-017 yy ==
»» exp(acos(0.3))
exp(acos(0.3))
ans 11
ans ==
3.5470
3.5470 »» zz == asin(y)
asin(y) () parentheses for
zz == function inputs
1.5708
1.5708
A Note about Workspace:
Numbers stored in double-precision floating point format
9
 General Functions
whos: List current variables
clear: Clear variables and functions from
memory
Close: Closes last figures
cd: Change current working directory
dir: List files in directory
echo: Echo commands in M-files
format: Set output format
10
 Getting help
help command (>>help)
lookfor command (>>lookfor)
Help Browser (>>doc)
helpwin command (>>helpwin)
Search Engine
Printable Documents
 “Matlabroot\help\pdf_doc\”
Link to The MathWorks
11
 Matrices
Entering and Generating Matrices
Subscripts
Scalar Expansion
Concatenation
Deleting Rows and Columns
Array Extraction
Matrix and Array Multiplication

12
 Entering Numeric Arrays
»» a=[1
a=[1 2;3
2;3 4]
4]
aa == Use square
11 22 brackets [ ]
33 44
Row separator »» b=[-2.8,
b=[-2.8, sqrt(-7),
sqrt(-7), (3+5+6)*3/4]
(3+5+6)*3/4]
semicolon (;) bb ==
-2.8000
-2.8000 00 ++ 2.6458i
2.6458i 10.5000
10.5000
Column separator » b(2,5) = 23
» b(2,5) = 23
space / comma (,)
bb ==
-2.8000
-2.8000 00 ++ 2.6458i
2.6458i 10.5000
10.5000 00 00
00 00 00 00 23.0000
23.0000

•Any MATLAB expression can be entered as a matrix

element
•Matrices must be rectangular. (Set undefined elements to
13 zero)
 The Matrix in MATLAB
Columns
(n)
1 2 3 4 5
A= 4
1
10
6
1
11
6
16
2
21
A (2,4)
1
2 7 12 17 22
2 8 1.2 9 4 25

Rows (m) 3 7.2 3 5 8

7 13
1 18
11 23 A (17)
4 0 4
0.5 9 4 14
5 19
56 24
5 10 15 20 25
5 23 83 13 0 10 Rectangular Matrix:
Scalar: 1-by-1 array
Vector: m-by-1 array
1-by-n array
Matrix: m-by-n array
14
 Entering Numeric Arrays
Scalar expansion »» w=[1
w=[1 2;3
2;3 4]
4] ++ 55
ww ==
66 77
88 99
Creating »» xx == 1:5
1:5
sequences: xx ==
colon operator (:) 11 22 33 44 55
»» yy == 2:-0.5:0
2:-0.5:0
yy ==
2.0000
2.0000 1.5000
1.5000 1.0000
1.0000 0.5000
0.5000 00
»» zz == rand(2,4)
rand(2,4)
Utility functions for
zz ==
creating matrices.
0.9501
0.9501 0.6068
0.6068 0.8913
0.8913 0.4565
0.4565
0.2311
0.2311 0.4860
0.4860 0.7621
0.7621 0.0185
0.0185

15
 Numerical Array
Concatenation
Use [ ] to combine »» a=[1
a=[1 2;3
2;3 4]
4]
existing arrays as aa == Use square
matrix “elements” 11 22
brackets [ ]
33 44
»» cat_a=[a,
cat_a=[a, 2*a;
2*a; 3*a,
3*a, 4*a;
4*a; 5*a,
5*a, 6*a]
6*a]
Row separator:
cat_a
cat_a ==
semicolon (;) 11 22 22 44
33 44 66 88
Column separator: 33 66 44 88 4*a
space / comma (,) 99 12
12 12
12 16
16
55 10
10 66 12
12
15
15 20
20 18
18 24
24

Note:
16 The resulting matrix must be rectangular
 Deleting Rows and Columns
»» A=[1
A=[1 55 9;4
9;4 33 2.5;
2.5; 0.1
0.1 10
10 3i+1]
3i+1]
AA ==
1.0000
1.0000 5.0000
5.0000 9.0000
9.0000
4.0000
4.0000 3.0000
3.0000 2.5000
2.5000
0.1000
0.1000 10.0000
10.0000 1.0000+3.0000i
1.0000+3.0000i
»» A(:,2)=[]
A(:,2)=[]
AA ==
1.0000
1.0000 9.0000
9.0000
4.0000
4.0000 2.5000
2.5000
0.1000
0.1000 1.0000
1.0000 ++ 3.0000i
3.0000i
»» A(2,2)=[]
A(2,2)=[]
???
??? Indexed
Indexed empty
empty matrix
matrix assignment
assignment is
is not
not allowed.
allowed.

17
 Array Subscripting / Indexing

1 2 3 4 5
A= 4
1
10
6
1
11
6
16
2
21

1
2
2 8 1.2 7 9 12 4 17
25 22
A(1:5,5) A(1:end,end)
3 7.2 3 5 8
7 13 1 18
11 23 A(:,5) A(:,end)
A(21:25) A(21:end)’
A(3,1) 4 0 4
0.5 9 4 14 5 19
56 24
A(3) 5 10 15 20 25
5 23 83 13 0 10
A(4:5,2:3)
A([9 14;10 15])

18
 Matrix Multiplication
»» aa == [1
[1 22 33 4;
4; 55 66 77 8];
8]; [2x4]
»» bb == ones(4,3);
ones(4,3); [4x3]
»» cc == a*b
a*b [2x4]*[4x3] [2x3]
cc ==
10
10 10
10 10
10
26
26 26
26 26
26 a(2nd row).b(3rd column)

Array Multiplication
»» aa == [1
[1 22 33 4;
4; 55 66 77 8];
8];
»» bb == [1:4;
[1:4; 1:4];
1:4];
»» cc == a.*b
a.*b
cc ==
11 44 99 16
16
55 12
12 21
21 32
32 c(2,4) = a(2,4)*b(2,4)

19
 Matrix Manipulation Functions
• zeros: Create an array of all zeros
• ones: Create an array of all ones
• eye: Identity Matrix
• rand: Uniformly distributed random numbers
• diag: Diagonal matrices and diagonal of a matrix
• size: Return array dimensions
• fliplr: Flip matrices left-right
• flipud: Flip matrices up and down
• repmat: Replicate and tile a matrix
20
 Matrix Manipulation Functions
• transpose (’): Transpose matrix
• rot90: rotate matrix 90
• tril: Lower triangular part of a matrix
• triu: Upper triangular part of a matrix
• cross: Vector cross product
• dot: Vector dot product
• det: Matrix determinant
• inv: Matrix inverse
• eig: Evaluate eigenvalues and eigenvectors
• rank: Rank of matrix
21
 Character Arrays (Strings)
»» str
str == 'Hi
'Hi there,'
there,'
str
str ==
Hi
Hi there,
there,
»» str2
str2 == 'Isn''t
'Isn''t MATLAB
MATLAB great?'
great?'
str2
str2 ==
Isn't
Isn't MATLAB
MATLAB great?
great?

str = H i t h e r e , 1x9 vector

22
 String Array
Concatenation
Using [ ] operator: »» str
str ='Hi
='Hi there,';
there,'; 1x9 vectors
Each row must be »» str1='Everyone!';
str1='Everyone!';
same length »» new_str=[str,
new_str=[str, '' ',', str1]
str1]
new_str
new_str ==
Row separator:
Hi
Hi there,
there, Everyone!
Everyone! vectors
1x19 vector
semicolon (;)
»» str2
str2 == 'Isn''t
'Isn''t MATLAB
MATLAB great?';
great?';
Column separator: »» new_str2=[new_str;
new_str2=[new_str; str2]
str2]
space / comma (,) new_str2
new_str2 ==
Hi
Hi there,
there, Everyone!
Everyone!
Isn't MATLAB great?
2x19 matrix
Isn't MATLAB great?

For strings of different length:

• STRVCAT »» new_str3
new_str3 == strvcat(str,
strvcat(str, str2)
str2)
• char new_str3
new_str3 ==
Hi
Hi there,
there, 2x19 matrix
Isn't
Isn't MATLAB
MATLAB great?
great? (zero padded)
23
 Working with String Arrays
String Comparisons
 strcmp: compare whole strings
 strncmp: compare first ‘N’ characters
 findstr: finds substring within a larger string
Converting between numeric & string arrays:
 num2str: convert from numeric to string array
 str2num: convert from string to numeric array

24
 Elementary Math
Logical Operators

Math Functions

Polynomial and Interpolation

25
 Logical Operations

= = equal to »» Mass
Mass == [-2
[-2 10
10 NaN
NaN 30
30 -11
-11 Inf
Inf 31];
31];
»» each_pos
each_pos == Mass>=0
Mass>=0
> greater than
each_pos
each_pos ==
< less than 00 11 00 11 00 11 11
»» all_pos
all_pos == all(Mass>=0)
all(Mass>=0)
>= Greater or equal all_pos =
all_pos =
<= less or equal 00
»» all_pos
all_pos == any(Mass>=0)
any(Mass>=0)
~ not all_pos
all_pos ==
11
& and
»» pos_fin
pos_fin == (Mass>=0)&(isfinite(Mass))
(Mass>=0)&(isfinite(Mass))
| or pos_fin
pos_fin ==
00 11 00 11 00 00 11
isfinite(), etc. . . .
all(), any() Note:
• 1 = TRUE
find
• 0 = FALSE
26
 Elementary Math Function
• abs, sign: Absolute value and Signum
Function
• sin, cos, asin, acos…: Triangular functions
• exp, log, log10: Exponential, Natural and
Common (base 10) logarithm
• ceil, floor: Round toward infinities
• fix: Round toward zero

27
Elementary Math Function
round: Round to the nearest integer
gcd: Greatest common devisor
lcm: Least common multiple
sqrt: Square root function
real, imag: Real and Image part of
complex
rem: Remainder after division
 Elementary Math Function
• max, min: Maximum and Minimum of arrays
• mean, median: Average and Median of arrays
• std, var: Standard deviation and variance
• sort: Sort elements in ascending order
• sum, prod: Summation & Product of Elements
• trapz: Trapezoidal numerical integration
• cumsum, cumprod: Cumulative sum, product
• diff, gradient: Differences and Numerical
Gradient

28
 Polynomials and
Interpolation
Polynomials
 Representing
 Roots (>> roots)
 Evaluation (>> polyval)
 Derivatives (>> polyder)
 Curve Fitting (>> polyfit)
 Partial Fraction Expansion (residue)
Interpolation
 One-Dimensional (interp1)
 Two-Dimensional (interp2)
29
 Example
polysam=[1
polysam=[1 00 00 8];
8];
roots(polysam)
roots(polysam)
ans
ans ==
-2.0000
-2.0000
1.0000
1.0000 ++ 1.7321i
1.7321i
1.0000 - 1.7321i
1.0000 - 1.7321i
Polyval(polysam,[0
Polyval(polysam,[0 11 2.5
2.5 44 6.5])
6.5])
ans
ans ==
8.0000
8.0000 9.0000
9.0000 23.6250
23.6250 72.0000
72.0000 282.6250
282.6250
polyder(polysam)
polyder(polysam)
ans
ans ==
33 00 00
[r
[r pp k]=residue(polysam,[1
k]=residue(polysam,[1 22 1])
1])
rr == 33 77
pp == -1
-1 -1-1
kk == 11 -2-2

30
 Example
xx == [0:
[0: 0.1:
0.1: 2.5];
2.5];
yy == erf(x);
erf(x);
pp == polyfit(x,y,6)
polyfit(x,y,6)
pp ==
0.0084
0.0084 -0.0983
-0.0983 0.4217
0.4217 -0.7435
-0.7435 0.1471
0.1471 1.1064
1.1064 0.0004
0.0004

interp1(x,y,[0.45
interp1(x,y,[0.45 0.95
0.95 2.2
2.2 3.0])
3.0])
ans
ans ==
0.4744
0.4744 0.8198
0.8198 0.9981
0.9981 NaN
NaN
31
 Importing and Exporting
Data
Using the Import Wizard

Using Save and Load command

save
save fname
fname load
load fname
fname
save
save fname
fname xx yy zz load
load fname
fname xx yy zz
save
save fname
fname -ascii
-ascii load
load fname
fname -ascii
-ascii
save
save fname
fname -mat
-mat load
load fname
fname -mat
-mat

32
 Input/Output for Text File

•Read formatted data, reusing the

format string N times.
»»[A1…An]=textread(filename,format,N)
[A1…An]=textread(filename,format,N)

•Import and Exporting Numeric Data

with General ASCII delimited files
»» MM == dlmread(filename,delimiter,range)
dlmread(filename,delimiter,range)

33
 Input/Output for Binary File
fopen: Open a file for input/output
fclose: Close one or more open files
fread: Read binary data from file
fwrite: Write binary data to a file
fseek: Set file position indicator
»» fid= fopen('mydata.bin',, 'wb');
fid= fopen('mydata.bin' 'wb');
»» fwrite (fid,eye(5) , 'int32');
fwrite (fid,eye(5) , 'int32');
»» fclose
fclose (fid);
(fid);
»» fid=
fid= fopen('mydata.bin',, 'rb');
fopen('mydata.bin' 'rb');
»» M= fread(fid, [5 5], 'int32')
M= fread(fid, [5 5], 'int32')
»» fclose
fclose (fid);
(fid);

34
 Graphics
Fundamentals

35
 Graphics
Basic Plotting
plot, title, xlabel, grid,
legend, hold, axis
Editing Plots
Property Editor
Mesh and Surface Plots
meshgrid, mesh, surf,
colorbar, patch, hidden
Handle Graphics
36
 2-D Plotting
Syntax:
plot(x1,
plot(x1, y1,
y1, 'clm1',
'clm1', x2,
x2, y2,
y2, 'clm2',
'clm2', ...)
...)

Example:
x=[0:0.1:2*pi];
x=[0:0.1:2*pi];
y=sin(x);
y=sin(x);
z=cos(x);
z=cos(x);
plot(x,y,x,z,'linewidth',2)
plot(x,y,x,z,'linewidth',2)
title('Sample
title('Sample Plot','fontsize',14);
Plot','fontsize',14);
xlabel('X
xlabel('X values','fontsize',14);
values','fontsize',14);
ylabel('Y
ylabel('Y values','fontsize',14);
values','fontsize',14);
legend('Y
legend('Y data','Z
data','Z data')
data')
grid
grid on
on

37
 Sample Plot

Title

Ylabel
Grid

Legend

Xlabel
38
 Subplots
Syntax: subplot(rows,cols,index)
subplot(rows,cols,index)

»»subplot(2,2,1);
subplot(2,2,1);
»» ……

»»subplot(2,2,2)
subplot(2,2,2)
»» ...
...

»»subplot(2,2,3)
subplot(2,2,3)
»» ...
...

»»subplot(2,2,4)
subplot(2,2,4)
»» ...
...

39
 Surface Plot Example
xx == 0:0.1:2;
0:0.1:2;
yy == 0:0.1:2;
0:0.1:2;
[xx,
[xx, yy]
yy] == meshgrid(x,y);
meshgrid(x,y);
zz=sin(xx.^2+yy.^2);
zz=sin(xx.^2+yy.^2);
surf(xx,yy,zz)
surf(xx,yy,zz)
xlabel('X
xlabel('X axes')
axes')
ylabel('Y axes')
ylabel('Y axes')

40
 3-D Surface Plotting
contourf-colorbar-plot3-waterfall-contour3-mesh-surf

41
 Specialized Plotting Routines
bar-bar3h-hist-area-pie3-rose

42
 Editing and Debugging M-
Files
What is an M-File?
The Editor/Debugger
Search Path
Debugging M-Files
 Types of Errors (Syntax Error and Runtime
Error)
 Using keyboard and “ ; ” statement
 Setting Breakpoints
 Stepping Through
 Continue, Go Until Cursor, Step, Step In, Step Out
 Examining Values
 Selecting the Workspace
 Viewing Datatips in the Editor/Debugger
 Evaluating a Selection

43
 Debugging
Select
Workspace
Set Auto-
Breakpoints

tips

44
 Programming and
Application
Development

45
 Script and Function Files
• Script Files
• Work as though you typed commands into
MATLAB prompt
• Variable are stored in MATLAB workspace

• Function Files
• Let you make your own MATLAB Functions
• All variables within a function are local
• All information must be passed to functions as
parameters
46 • Subfunctions are supported
 Basic Parts of a Function M-File
Output Arguments Function Name Input Arguments

Online Help function y = mean (x)

% MEAN Average or mean value.
% For vectors, MEAN(x) returns the mean value.
% For matrices, MEAN(x) is a row vector
% containing the mean value of each column.
[m,n] = size(x);
if m == 1
Function Code
m = n;
end
y = sum(x)/m;

47
 Flow Control Statements

if Statement
if
if ((attendance
((attendance >=
>= 0.90)
0.90) && (grade_average
(grade_average >=
>= 60))
60))
pass
pass == 1;
1;
end;
end;

while Loops
eps
eps == 1;
1;
while
while (1+eps)
(1+eps) >> 11
eps
eps == eps/2;
eps/2;
end
end
eps
eps == eps*2
eps*2

48
 Flow Control Statements
for Loop
aa == zeros(k,k)
zeros(k,k) %% Preallocate
Preallocate matrix
matrix
for
for mm == 1:k
1:k
for
for nn == 1:k
1:k
a(m,n)
a(m,n) == 1/(m+n
1/(m+n -1);
-1);
end
end
end
end

switch Statement
method
method == 'Bilinear';
'Bilinear';
switch
switch lower(method)
lower(method)
case
case {'linear','bilinear'}
{'linear','bilinear'}
disp('Method
disp('Method is
is linear')
linear')
case
case 'cubic'
'cubic'
disp('Method
disp('Method is
is cubic')
cubic')
otherwise
otherwise
disp('Unknown
disp('Unknown method.')
method.')
end
end
Method
Method is
is linear
49 linear
 M-file Programming Features
SubFunctions
Varying number of input/output arguments
Local and Global Variables
Obtaining User Input
 Prompting for Keyboard Input
 Pausing During Execution
Errors and Warnings
 Displaying error and warning Messages
Shell Escape Functions (! Operator)
Optimizing MATLAB Code
 Vectorizing loops
 Preallocating Arrays
50
 Function M-file
function
function rr == ourrank(X,tol)
ourrank(X,tol)
%% rank
rank ofof aa matrix
matrix Multiple Input Arguments
ss == svd(X);
svd(X); use ( )
if
if (nargin
(nargin ==== 1)
1) »r=ourrank(rand(5),.1);
tol
tol == max(size(X))
max(size(X)) ** s(1)*
s(1)* eps;
eps;
end
end
rr == sum(s
sum(s >> tol);
tol);

function
function [mean,stdev]
[mean,stdev] == ourstat(x)
ourstat(x)
[m,n]
[m,n] == size(x);
size(x);
Multiple Output
if
if mm ==
== 11
Arguments, use [ ]
m = n;
»[m std]=ourstat(1:9); endm = n;
end
mean
mean == sum(x)/m;
sum(x)/m;
stdev
stdev == sqrt(sum(x.^2)/m
sqrt(sum(x.^2)/m –– mean.^2);
mean.^2);
51
 Data Types

Numeric Arrays
Multidimensional Arrays
Structures and Cell Arrays
52
 Multidimensional Arrays
The first references array dimension
1, the row.
»» AA == pascal(4);
pascal(4);
The second references dimension 2, » A(:,:,2) = magic(4)
» A(:,:,2) = magic(4)
the column. A(:,:,1)
A(:,:,1) ==
The third references dimension 3, 11 11 11 11
The page. 11 22 33 44
1 0 0 0 11 33 66 10
10
0 1 0 0
11 44 10
10 20
20
0 0 1 0
0 0 0 1 A(:,:,2)
A(:,:,2) ==
16
16 22 33 13
13
0 0 0 0 Page N
55 11
11 10
10 88
16 20 30 130 0
110 1
100 80 0 99 77 66 12
12
1 15 1
9 70 60 120 0
44 14 15 11
1 2 3 4 14 15
1 3
4 14 15
6 10
1
»» A(:,:,9)
A(:,:,9) ==
1 4 10 20 diag(ones(1,4));
diag(ones(1,4));
53 Page 1
 Structures
•Arrays with named data containers called fields.
»» patient.name='John
patient.name='John Doe';
Doe';
»» patient.billing
patient.billing == 127.00;
127.00;
»» patient.test=
patient.test= [79
[79 75
75 73;
73;
180
180 178
178 177.5;
177.5;
220
220 210
210 205];
205];

•Also, Build structure arrays using the struct function.

•Array of structures
»» patient(2).name='Katty
patient(2).name='Katty Thomson';
Thomson';
»» Patient(2).billing
Patient(2).billing == 100.00;
100.00;
»» Patient(2).test=
Patient(2).test= [69
[69 25
25 33;
33; 120
120 128
128 177.5;
177.5; 220
220
210
210 205];
205];
54
 Cell Arrays
•Array for which the elements are cells and can hold
other MATLAB arrays of different types.

»» A(1,1)
A(1,1) == {[1
{[1 44 3;
3;
00 55 8;
8;
77 22 9]};
9]};
»» A(1,2)
A(1,2) == {'Anne
{'Anne Smith'};
Smith'};
»» A(2,1)
A(2,1) == {3+7i};
{3+7i};
»» A(2,2)
A(2,2) == {-pi:pi/10:pi};
{-pi:pi/10:pi};

•Using braces {} to point to elements of cell array

•Using celldisp function to display cell array
55
 Nonlinear Numerical Functions
• inline function »» ff == inline('3*sin(2*x.^2)','x')
inline('3*sin(2*x.^2)','x')
ff ==
Inline
Inline function:
function:
Use char function f(x)
f(x) == 3*sin(2*x.^2)
3*sin(2*x.^2)
to convert inline »» f(2)
f(2)
object to string ans
ans ==
2.9681
2.9681
• Numerical Integration using quad
»» QQ == quad('1./(x.^3-2*x-5)',0,2);
quad('1./(x.^3-2*x-5)',0,2);
»» FF == inline('1./(x.^3-2*x-5)');
inline('1./(x.^3-2*x-5)');
»» QQ == quad(F,0,2);
quad(F,0,2);
»» QQ == quad('myfun',0,2)
quad('myfun',0,2)
Note:
quad function use adaptive function y = myfun(x)

56 Simpson quadrature y = 1./(x.^3-2*x-5);

 Nonlinear Numerical
Functions
fzero finds a zero of a single variable function
[x,fval]= fzero(fun,x0,options)
 fun is inline function or m-function
fminbnd minimize a single variable function
on a fixed interval. x1<x<x2

fminsearch
[x,fval]= minimize
fminbnd(fun,x1,x2,options)
a several variable
function
[x,fval]= fminsearch(fun,x0,options)
Use optimset to determine options parameter.

57 options = optimset('param1',value1,...)
 Ordinary Differential
Equations
(Initial Value Problem)
An explicit ODE with initial value:

Using ode45 for non-stiff functions and

ode23t for stiff functions.
[t,y] = solver(odefun,tspan,y0,options)
function dydt = odefun(t,y) Initialvlue

[initialtime finaltime]

58 •Use odeset to define options parameter

 ODE Example:

function dydt=myfunc(t,y)
dydt=zeros(2,1);
dydt(1)=y(2);
dydt(2)=(1-y(1)^2)*y(2)-y(1);

»» [t,y]=ode45('myfunc',[0
[t,y]=ode45('myfunc',[0 20],[2;0])
20],[2;0])
3

Note:
1

Help on odeset to set options

0

for more accuracy and other

-1

useful utilities like drawing

-2

results during solving.

59 -3
0 2 4 6 8 10 12 14 16 18 20
 Handle Graphics
Graphics in MATLAB consist of objects:
 root, figure, axes, image, line, patch,
rectangle, surface, text, light
Creating Objects
Setting Object Properties Upon
Creation
Obtaining an Object’s Handles
Knowing Object Properties
Modifying Object Properties
 Using Command Line
 Using Property Editor
60
 Graphics Objects
Root
object

Figure
object
UIMenu
objects UIControl
Axes object
objects

Surface
object

Line
objects

Text
objects 61
 Obtaining an Object’s Handle
1. Upon Creation
h_line = plot(x_data, y_data, ...)

2. Utility Functions
What is the current object?
0 - root object handle • Last object created
• OR
gcf - current figure handle
• Last object clicked
gca- current axis handle
gco- current object handle

3. h_obj
FINDOBJ
= findobj(h_parent, 'Property', 'Value', ...)

63 Default = 0 (root object)

 Modifying Object Properties
• Obtaining a list of current properties:
get(h_object)

• Obtaining a list of settable properties:

set(h_object)

• Modifying an object’s properties

 Using Command Line
set(h_object,'PropertyName','New_Value',...)

 Using Property Editor

64
 Graphical User Interface
What is GUI?
What is figure and *.fig file?
Using guide command
GUI controls
GUI menus

65
 Axes static text
Frames

Checkbox Slider

Edit text

66 Radio Buttons Push Buttons

 Guide Editor

Property Inspector

Result Figure
67
 Conclusion
 Matlab is a language of technical computing.
 Matlab, a high performance software, a high-
level language
 Matlab supports GUI, API, and …
 Matlab Toolboxes best fits different applications
 Matlab …

68
 Getting more help
• Contact http://www.mathworks.com/support
• You can find more help and FAQ about
mathworks products on this page.
• Contact comp.soft-sys.matlab Newsgroup
• Using Google Groups Page to Access this page
http://groups.google.com/

69
 ?

70