TM

Tutorial di Matlab

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM—

 Avvio di Matlab

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 2

 Avvio di Matlab

Si noti che Matlab esegue i calcoli e

memorizza le variabili (scalari, vettoriali o
matriciali) in doppia precisione

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 3

 Matrici

Per lanciare l’Array Editor è sufficiente cliccare due volte sull’oggetto

matrice (A nel nostro caso) della finestra Workspace

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 4

 Editor di Matlab
Cliccando sull’icona “New M-file” della barra strumenti in alto a sinistra viene lanciato
l’editor di Matlab con cui si possono scrivere e salvare dei file funzione o degli script
(sequenze di comandi Matlab) richiamabili in un secondo momento.

Il direttorio di default dove Matlab salva i file dell’utente è: \Work.

Aggiungendo un nuovo direttorio alla lista dei direttori Matlab ricerca i file funzione o di
script anche in quel direttorio. È buona norma evitare di salvare i propri file nel direttorio di
default \Work al fine di mantenere distinti i singoli progetti ed argomenti di ricerca.

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 5

 Vettori

Il comando clear cancella tutte le variabili presenti nel Workspace.

N.B.: quando si inizia a risolvere un nuovo problema è buona norma cancellare tutte le
variabili presenti nel Workspace per evitare inutili errori di riutilizzo o sovrapposizione di
variabili già definite precedentemente.

Per effettuare un prodotto elemento per elemento tra vettori

di uguali dimensioni si utilizza l’operatore: .*

Se alla fine di un’istruzione si aggiunge il segno ; allora

Matlab non presenta a video il risultato dell’istruzione anche se
la esegue e ne memorizza il valore.

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 6

 Altri comandi di Matlab
Il comando clc cancella invece soltanto il contenuto della finestra di Comando. Le variabili
fino a quel punto introdotte NON vengono cancellate.

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 7

 Diagrammi

7 7

6.8 6.8

6.6 6.6

6.4 6.4

6.2 6.2

6 6

5.8 5.8

5.6 5.6

5.4 5.4

5.2 5.2

5 5
1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 8

 Diagrammi

fplot('tanh',[-2 2])

8 1

7.5 0.8

7 0.6

6.5 0.4

6 0.2

5.5 0

5 -0.2

4.5 -0.4

4 -0.6

3.5 -0.8

3 -1
5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 9

 Diagrammi
% anno e cap sono i vettori con i valori delle ascisse ed ordinate
% 'r-o' stile della curva: "r"=red, "-"=linea continua, "o"=pallini per ogni dato
% grid = griglia per entrambi gli assi
% xlabel, ylabel = descrizione per gli assi delle ascisse e ordinate
% title = titolo del grafico
plot(anno,cap,'r-o'),grid,xlabel('anni'),ylabel('€'),title('deposito bancario');
4
x 10 depos ito banca rio
3.5

2.5
€

1.5

1
0 5 10 15 20 25 30 35
anni

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 10

 Help di Matlab

Help Æ MATLAB Help

Vedi anche:
Demos

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 11

 Diagrammi di funzioni

0.8

0.6

0.4

0.2

-0.2

-0.4

-0.6

-0.8

-1
-3 -2 -1 0 1 2 3

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 12

Variabili e funzioni predefinite (intrinseche)
i, j unità immaginaria

pi pigreco

Inf infinito (ad es.: 3/0)

NaN not a number (as es.: 0/0, Inf/Inf)

sin, cos, tan, asin, acos funzioni trigonometriche (in radianti)

sinh, cosh, asinh, acosh funzioni iperboliche

sqrt, log, log10, exp funzioni varie

rem resto della divisione intera (divisione–modulo)

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 13

 Costrutti if for while
if x >= 0 s = 0.; x = 3.;

y = sqrt(x); for k = 1:100 while x < 25.

elseif x == 0 s = s + 1./k; x = x + 2;

y = 0; end end

else disp(x)

y = NaN;

end for k = 1:3:100

...

end

Operatori di confronto: < <= > >= == ~=

Operatori logici: ~ [NOT], & [AND] , | [OR]

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 14

 Costrutto switch
d = floor(3*rand) + 1 z = floor(10*rand)

switch d switch z

case 1 case {1, 3, 5, 7, 9}

disp(‘È un 1’) disp(‘Dispari’)

case 2 case {2, 4, 6, 8}

disp(‘È un 2’) disp(‘Pari’)

otherwise otherwise

disp(‘È un 3’) disp(‘Zero’)

end end

N.B.: per produrre il simbolo: { premere ALT + 123 (tastierino numerico)

per produrre il simbolo: } premere ALT + 125 (tastierino numerico)

per produrre il simbolo: | premere ALT + 124 (tastierino numerico)

per produrre il simbolo: ~ premere ALT + 126 (tastierino numerico)

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 15
 Output formattato su video
>> fprintf('Pigreco = %g ',pi)
Pigreco = 3.14159
>> fprintf('Pigreco = %15.2f ',pi)
Pigreco = 3.14
>> fprintf('Pigreco = %15.4f ',pi)
Pigreco = 3.1416
>> fprintf('Pigreco = %15.8f ',pi)
Pigreco = 3.14159265
>> fprintf('Pigreco = %15.10f ',pi)
Pigreco = 3.1415926536
>> fprintf('Pigreco = %20.16f ',pi)
Pigreco = 3.1415926535897931
>> fprintf('Pigreco = %15.2e ',pi)
Pigreco = 3.14e+000
>> fprintf('Pigreco = %15.16e ',pi)
Pigreco = 3.1415926535897931e+000

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 16

 Output formattato su file

unit = fopen(‘NomeFile.ext’, ‘w’);

fprintf(unit, ‘\nPigreco = %15.10f ',pi);

... ... ...

fclose(unit);

Utilizzando i comandi fwrite e fread è anche possibile scrivere, ‘w’, o leggere, ‘r’,
su o da file binari.

Si veda anche (utilizzando l’help) i comandi: save e load.

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 17

 Buone norme di programmazione
Qualora si debba scrivere un numero reale è opportuno nei vari linguaggi di
programmazione aggiungere il punto decimale.
Esempio: (1./4.) anziché (1/4)

Infatti, anche se Matlab perdona l’utente sbadato, non altrettanto avviene in altri linguaggi
di programmazione quali il Fortran o il C++.
Programmando in Fortran ad esempio: 1/4 = 0 in quanto la divisione viene eseguita tra
due cifre intere ed il risultato è automaticamente posto uguale ad un intero.
Attenzione quindi: Pr^(1/3) = 1. qualsiasi sia il valore di Pr.

La scrittura corretta è viceversa: Pr^(1./3.).

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 18

 Buone norme di programmazione
Altro rischio programmativo è rappresentato dagli operatori:
* oppure .*

/ oppure ./

^ oppure .^

Gli operatori: *, /, ^ operano a livello vettoriale–matriciale in linea con l’analisi classica.

Viceversa gli operatori: .*, ./, .^ operano a livello vettoriale–matriciale sui singoli elementi
dei vettori o matrici.

Per effettuare un prodotto tra matrici in senso classico (prodotto righe per colonne) si
utilizza l’istruzione: C = A * B;

Viceversa per effettuare il prodotto elemento per elemento: C = A .* B;

N.B.: Per evitare ambiguità tra l’utilizzo del punto come separatore decimale oppure come
operatore congiunto ai simboli: *, /, ^ è opportuno mantenere gli operatori: .*, ./, .^ ben
separati dalle variabili sulle quali essi operano. Esempio: 1. / 4. * pi * a ^ 2

N.B.: elevando a potenza intera non si deve aggiungere il punto. Esempio: a^2, b^5

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 19

 HOW TO
HT1

Domanda: come produrre dei punti equispaziati secondo n intervalli comprendendo gli
estremi?

Risposta: innanzitutto evitare assolutamente il costrutto Matlab™:

x = [xLow:delta:xUp].

Ad esempio la seguente istruzione è scorretta: x = [0.1: 0.01: 1.2] a causa della

precisione limitata della CPU, anche in doppia precisione.

Viceversa è corretto il seguente algoritmo:

Input: xLow, xUp, nIntervalli

Algoritmo corretto:
delta = (xUp – xLow) / nIntervalli;

for i = 1: nIntervalli + 1

x(i) = xLow + delta * (i-1);

next i

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 20

 HOW TO
HT1 continua

Il seguente algoritmo è meno corretto in quanto accumula gli errori di troncamento legati
all’operazione di somma. I termini x(i) sono cioè meno precisi.

Algoritmo poco preciso:

delta = (xUp – xLow) / nIntervalli;

x(1) = xLow;

for i = 1, nIntervalli

x(i+1) = x(i) + delta;

next i

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 21

 HOW TO
HT2

Domanda: come produrre n punti equispaziati che comprendano anche gli estremi di un
intervallo?

Risposta: innanzitutto evitare assolutamente il costrutto Matlab™:

x = [xLow:delta:xUp].

È corretto il seguente algoritmo:

Input: xLow, xUp, nPunti

Algoritmo:
delta = (xUp – xLow) / (nPunti – 1);

for i = 1: nPunti

x(i) = xLow + delta * (i-1);

next i

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 22

 HOW TO
HT3
Domanda: come funziona esattamente il costrutto while ?

Risposta: contrariamente a quanto qualcuno può pensare, il costrutto while continua ad

eseguire le istruzioni in esso contenute MENTRE la sua condizione è VERA.

Esempio: data la serie armonica (somma degli inversi dei numeri naturali) per sapere
quanti termini sono necessari prima di raggiungere almeno il valore 5 si scriverà:
n = 0;

somma = 0.;

while somma < 5.

n = n + 1;

somma = somma + 1. / n;

end

disp([‘ n = ’, num2str(n)]);

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 23

 HOW TO
HT3 continua

Sarebbe un grave errore utilizzare l’espressione:

while somma = 5.

Parimenti è un’imprecisione utilizzare l’espressione:

while somma <= 5.

in quanto il problema chiede chiaramente: “prima di raggiungere almeno il valore 5”. Ciò
sta a significare che il ciclo while deve terminare non appena è stato raggiunto o superato
il valore 5. Supponiamo, cosa non vera, che la somma dei primi 83 reciproci dei numeri
naturali conducesse esattamente al valore 5. Se si utilizzasse l’istruzione:
while somma <= 5.

la procedura while sarebbe autorizzata a iterare ancora una volta per sommare il termine
84-esimo che quindi andrebbe oltre il valore 5 richiesto dal problema.

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 24

 HOW TO
HT4
Domanda: come determinare il numero di righe e colonne di una matrice?
Risposta: occorre utilizzare l’istruzione size.

Esempio:
A = zeros(3,5)
nRows = size(A,1);
nCols = size(A,2);
disp(['nRows = ',num2str(nRows),' nCols = ',num2str(nCols)]);

Output:
A =
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
nRows = 3 nCols = 5

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 25

 HOW TO
HT5

Domanda: come determinare il numero di dimensioni di una matrice?

Risposta: occorre utilizzare le istruzioni size e length.

Esempio:
A = zeros(3,5);

nRowsCols = size(A)

nDimensions = length(nRowsCols)

Output:
nRowsCols =

3 5

nDimensions =

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 26

 HOW TO
HT6

Domanda: come salvare uno script (sequenza di istruzioni) in Matlab™ ?

Risposta: innanzitutto occorre comprendere su quale direttorio Matlab sta puntando

tramite il comando: pwd.

Salvando il file nel direttorio riportato da Matlab™ tramite il comando pwd sarà poi possibile
eseguirlo dalla linea di comando digitando il solo nome del file stesso.

Esempio:
>> pwd

ans =

X:\MATLAB6.5\work

Se si crea un nuovo file e lo si salva con il nome: prova.m, sarà sufficiente lanciare il
comando: prova per eseguirlo.

Se si vuole cambiare direttorio dove salvare e far eseguire i file di Matlab™ è possibile
utilizzare il comando chdir nomedir. Esempio: chdir c:\myfiles\ese1\

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 27

 HOW TO
HT7
Domanda: come disegnare con Matlab™ quattro diagrammi in un’unica finestra?
Risposta: si utilizza l’istruzione subplot(n,m,i) che permette di disegnare n righe ed m
colonne di diagrammi puntando nella fattispecie sul diagramma i-esimo contando da
sinistra a destra e dall’alto in basso. Il diagramma è poi ottenuto tramite la solita istruzione
plot.

Esempio:
subplot(2,2,1); % 4 diagrammi (2 x 2), punto al #1
plot(x1,y1),title ‘diagramma 1';
subplot(2,2,2); % punto al #2
plot(x2,y2),title ‘diagramma 2';
subplot(2,2,3); % punto al #3
plot(x3,y3),title ‘diagramma 3';
subplot(2,2,4); % punto al #4
plot(x4,y4),title ‘diagramma 4';

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 28

 HOW TO
HT8
Domanda: come si determina la parte intera di un numero in virgola mobile?
Risposta: provare i seguenti comandi di Matlab™: fix, round, floor, ceil.

Esempio:
x = [-5.8, -3.2, -2.5, 0., 0.5, 2.1, 3.9]
fixX = fix(x)
roundX = round(x)
floorX = floor(x)
ceilX = ceil(x)

Output:
x = -5.8000 -3.2000 -2.5000 0 0.5000 2.1000 3.9000

fixX = -5 -3 -2 0 0 2 3

roundX = -6 -3 -3 0 1 2 4

floorX = -6 -4 -3 0 0 2 3

ceilX = -5 -3 -2 0 1 3 4
Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 29
 HOW TO
HT9
Domanda: come si fa a sapere in quale direttorio Matlab™ sta operando?
Risposta: si utilizza il comando: pwd
Esempio:
>> pwd
ans =
E:\MyDocuments\Corsi\CDPDIC\Sources\Matlab\Ese3

HT10
Domanda: come si fa a sapere il nome dei file nel direttorio in cui si sta operando?
Risposta: si utilizza il comando: dir
Esempio:
>> dir
DvdSolveL.m Ese34.m Ese35.m Ese38.m UseDvdSolveA.m

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 30

 HOW TO
HT11
Domanda: come si fa a sapere a quale direttorio appartiene una certa funzione?
Risposta: si utilizza il comando: which FUN (oppure anche which FUN –ALL)
Esempio:
>> which fsolve
C:\MATLAB6p5\toolbox\optim\fsolve.m

HT12
Domanda: come si fa a cambiare direttorio tramite linea di comando?
Risposta: si utilizza il comando: cd
Esempio:
>> cd c:\windows\system32

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 31

 HOW TO
HT13
Domanda: come si fa a calcolare l’integrale analitico di una funzione?
Risposta: si utilizzano i comandi: syms e int
Esempio:
>> syms x
>> int(x^2)
ans =
1/3*x^3

HT14
Domanda: come si fa a calcolare la derivata analitica di una funzione
Risposta: si utilizzano i comandi: syms e diff
Esempio:
>> syms x
>> diff(x^3)
ans =
3*x^2
Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 32
 HOW TO
HT15

Domanda: come si fa a specificare i limiti inferiore e superiore di un asse in un diagramma

prodotto con il comando plot ?

Risposta: si utilizzano le istruzioni: XLim e/o YLim e il comando set

Esempio:

plot(xF,yF,‘r-'),grid,title('Funzione di Lagrange’);

set(gca,'YLim',[-0.5 1.5]);

Spiegazione: gca è il puntatore alla finestra attuale dove vengono presentati i grafici.
L’istruzione YLim indica che si vuole specificare un intervallo per l’asse delle ordinate. Il
vettore successivo [-0.5 1.5] definisce i valori minimo e massimo per tale asse.

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 33

 HOW TO
HT16

Domanda: come si fa a mostrare i caratteri speciali nelle scritte (ad esempio titolo ed assi di
un diagramma) ?

Risposta: si utilizzano i caratteri speciali TeX che premettono un carattere di backslash

Esempio:

plot(xF,yF,‘r-'), title('Funzione \alpha’); % lettera greca alfa

Altri esempi:

\beta \gamma \sim \infty \circ \Gamma \Delta \pm \forall

\exists \approx \neq

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 34

 HOW TO
HT17

Domanda: se durante l’esecuzione del codice si incontra un problema insormontabile e si

desidera bloccare il programma come si fa?

Risposta: si utilizza il comando exit

Esempio:

if portata < 0.

disp(‘La portata è negativa. Impossibile proseguire...’);

exit

end

Spiegazione: l’istruzione exit blocca l’esecuzione del programma, esce dallo script o dalla
funzione in cui Matlab stava operando e torna al prompt dei comandi.

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 35

 HOW TO
HT18

Domanda: come eseguire una serie di istruzioni iterativamente fintantoché una delle due
condizioni non è più vera?

Risposta: si utilizza il comando while con l’operatore di confronto &

Esempio:

while (b-a) < EPSI & iConto < MAX_ITER

iConto = iConto + 1;

... % serie di istruzioni iterative

end

Spiegazione: il ciclo while viene eseguito mentre le condizioni poste sono vere. Se ad
esempio si vuole proseguire con le iterazioni, fintantoché l’intervallo (b-a) è maggiore di
EPSI o fintantoché il numero massimo di iterazioni MAX_ITER non è stato raggiunto il ciclo
while riportato nell’esempio è ciò che fa al caso nostro.
Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 36
 HOW TO
HT19

Domanda: esiste un modo per memorizzare una informazione in una variabile di tipo
stringa?

Risposta: Sì. È molto semplice. Basta assegnare la stringa alfanumerica tra virgolette
semplici ad una nuova variabile.

Esempio:

unaStringa = ‘Soluzione non accettabile’;

messaggio = ‘Il programma abortisce per una grave errore.’;

disp(unaStringa);

disp(messaggio);

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 37

 Bibliografia
• http://www.eece.maine.edu/mm/matweb.html

• http://spicerack.sr.unh.edu/~mathadm/tutorial/software/matlab/

• http://www.engin.umich.edu/group/ctm/basic/basic.html

• http://www.mines.utah.edu/gg_computer_seminar/matlab/matlab.html

• http://www.math.ufl.edu/help/matlab-tutorial/

• http://www.indiana.edu/~statmath/math/matlab/

Davide Manca – Calcoli di Processo dell’Ingegneria Chimica – Politecnico di Milano TM— 38