Metodi numerici con elementi

di Programmazione
A.A. 2013-2014

Introduzione al MatLab
VI parte
1
Docente: Vittoria Bruni

Email: [email protected]

Ufficio: Via A. Scarpa,

Pal. B, I piano, Stanza n. 16
Tel. 06 49766648

Ricevimento: Giovedì 14.00-15.00

2
Il materiale didattico è disponibile sul sito http://ingaero.uniroma1.it/
nella pagina dedicata al corso Metodi Numerici con elementi di
Programmazione

Per consultazione: Getting Started with MatLab – The mathworks

www.mathworks.com

3
 Files .m
Matlab permette di definire nuove funzioni

I comandi (lista di istruzioni) che definiscono la funzione devono

essere inseriti in un file con estensione .m

• Function files: nome_file.m

• è possibile passare variabili in input
• al termine dell’esecuzione della funzione tutte le variabili locali
utilizzate vengono cancellate
• vengono restituite una o più variabili di output
 Esempio
Esercizio: Definire la funzione stat che calcoli la media e la deviazione
standard di un vettore x

Si crea il file di testo stat.m (nome del file in cui si memorizza la funzione)
con il seguente contenuto
-----------------------------------------------------------------------------------
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli
% elementi del vettore x
n = length(x);
mean = sum(x)/n;
stdev = sqrt(sum((x-mean).^2)/n);
------------------------------------------------------------------------------------------------
E si chiama dal Command Window nel modo seguente
>>x=rand(50,1);
>>[mean,stdev]=stat(x);
 Creazione di m-files

• Un file .m può essere generato con qualsiasi editor di testo

ASCII (ad es. notepad di windows o l’editor di Matlab)

• Per creare un nuovo m-file (come per lo script file):

– Menù File->New->M-File (oppure function file)

 Function files
• La prima riga del file definisce la sintassi della funzione
function [output]=myfunc(input)
output = elenco delle variabili richieste in output separate da virgole
Se non si richiede alcuna variabile in output si scrive solo []
subito dopo il comando function
input = elenco delle variabili di input separate da virgole
myfunc = nome della funzione
Nota: Il nome con cui viene richiamata la funzione dal
Command Window è quello del file in cui è memorizzata
Si consiglia chiamare il file .m in cui si memorizza la
funzione con il nome della funzione myfunc

• Nelle righe successive è consigliabile scrivere lo help della funzione

• Dopo lo help si scrive la sequenza di istruzioni che, a partire dagli input,

calcola gli output richiesti
 Function files
- Le variabili di input e di output possono essere variabili semplici,
vettori, matrici o altre strutture dati

- Il nome della funzione deve sempre iniziare con una lettera dello
alfabeto, può contenere numeri e underscore ma non può contenere
spazi

- La funzione termina con l’ultima istruzione (nelle versioni più recenti

di Matlab si consiglia l’uso del comando end) oppure al primo eventuale
comando return

- Dopo la prima riga di comando, si possono scrivere dei commenti alla

funzione
Per esempio, si può descrivere cosa fa la funzione, le variabili di
input richieste e quelle di output restituite. Tali commenti vanno
preceduti ad ogni riga dal simbolo % e vengono visualizzati come
help della funzione
 Esempio

-----------------------------------------------------------------------------------
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli elementi
% del vettore x (dato di input)

n = length(x);
mean = sum(x)/n;
stdev = sqrt(sum((x-mean).^2)/n);
------------------------------------------------------------------------------------------------
 Esempio

Nome della funzione (coincide al nome

dato al file .m)

-----------------------------------------------------------------------------------
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli elementi
% del vettore x (dato di input)

n = length(x);
mean = sum(x)/n;
stdev = sqrt(sum((x-mean).^2)/n);
------------------------------------------------------------------------------------------------
 Esempio
Elenco delle variabili da restituire in
output

Elenco delle variabili di input

-----------------------------------------------------------------------------------
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli elementi
% del vettore x (dato di input)

n = length(x);
mean = sum(x)/n;
stdev = sqrt(sum((x-mean).^2)/n);
------------------------------------------------------------------------------------------------
 Esempio

-----------------------------------------------------------------------------------
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli elementi
% del vettore x (dato di input)

Help della funzione

n = length(x);
mean = sum(x)/n;
Sequenza di istruzioni
stdev = sqrt(sum((x-mean).^2)/n);
------------------------------------------------------------------------------------------------
 Esempio

-----------------------------------------------------------------------------------
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli elementi
% del vettore x (dato di input)

n = length(x);
mean = sum(x)/n;
Sequenza di istruzioni
stdev = sqrt(sum((x-mean).^2)/n);
------------------------------------------------------------------------------------------------
Definizione delle variabili
da restitutire in output
 Esempio
Elenco delle variabili da restituire in
output
Nome della funzione (coincide al nome
dato al file .m)
Elenco delle variabili di input
-----------------------------------------------------------------------------------
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli elementi
% del vettore x (dato di input)

Help della funzione

n = length(x);
mean = sum(x)/n;
Sequenza di istruzioni
stdev = sqrt(sum((x-mean).^2)/n);
------------------------------------------------------------------------------------------------
Definizione delle variabili
da restitutire in output
 Esempio
------------------------------------------------------------------------------------------------
Scritta e salvata nel file stat.m, si può chiamare dal Command Window nel
modo seguente
Prima di chiamare la funzione è
>>x=rand(50,1); necessario definire la variabile da
dare in input
>>[mean,stdev]=stat(x);
La funzione stat è chiamata sulla
variabile x e assegna gli output
alle variabili mean e stdev che
sono create solo quando è invocata
la funzione stat

Le istruzioni precedenti sono equivalenti a

E’ importante l’ordine e il tipo
>>y=rand(50,1); delle variabili di input e output e
>>[my,sdy]=stat(y); non il nome (che può essere
diverso da quello usato nel file .m)
 Function files
- Una funzione può essere chiamata dal Command window oppure
all’interno di uno script o di un’altra funzione
Oss: Per usare una funzione all’interno di uno script, la funzione
deve essere nella stessa directory dello script o deve essere
una funzione predefinita (contenuta nel path)
- Le variabili definite all’interno della funzione sono locali e si
perdono al termine dell’esecuzione della funzione, tranne le
variabili richieste in output

Per esempio la variabile n della funzione stat non sarà contenuta nella
lista delle variabili dello workspace perché non è uno degli output
function [mean,stdev]=stat(x)
% Calcola la media (mean) e la deviazione standard (stdev) degli elementi
% del vettore x (dato di input)
n = length(x);
mean = sum(x)/n;
stdev = sqrt(sum((x-mean).^2)/n);
 Esercizi
Scrivere una funzione che calcoli il fattoriale di un numero
intero n. Il nome della funzione è fatt

Riprendiamo lo script myfactorial.m Come si deve modificare

in modo tale che diventi una funzione?
% my_factorial.m
% script per il calcolo del fattoriale di un intero n introdotto dall’utente.
% La variabile fatt contiene il valore n!

n = input(‘introduci il numero intero positivo n = ’);

fact = 1;
for i = 2:n
fact = fact * i;
end
 Esercizi
Scrivere una funzione che calcoli il fattoriale di un numero
intero n. Il nome della funzione è fatt

Riprendiamo lo script myfactorial.m Come si deve modificare

in modo tale che diventi una funzione?
% my_factorial.m
% script per il calcolo del fattoriale di un intero n introdotto dall’utente.
% La variabile fatt contiene il valore n!

n = input(‘introduci il numero intero positivo n = ’);

fact = 1;
for i = 2:n
da restituire come output
fact = fact * i;
end
function [fact] = fatt(n)
Esercizi
% fatt(n)
% funzione per il calcolo del fattoriale di un intero n
% INPUT
% n = numero intero
% OUTPUT
% fact = fattoriale di n (n!)
fact = 1;
for i = 2:n
fact = fact * i;
end

Dal Command window

>> f = fatt(6);
>> disp(f)
 Esercizio
Scrivere una function matlab (fun_eqlin.m) per il calcolo della soluzione
dell’ equazione lineare di primo grado ax+b=0
Soluzione: I parametri di input della funzione sono i coefficienti a e b
dell’equazione lineare. Lo output della funzione è la soluzione x

function [x] = fun_eqlin(a,b)

% [x] = fun_eqlin(a,b)
% Funzione per il calcolo della soluzione dell’ equazione di primo grado ax+b=0
% La funzione stampa un messaggio nei casi di equazione impossibile o indeterminata
% INPUT
% a = primo coefficiente dell’equazione
% b = termine noto dell’equazione
% OUTPUT
% x = soluzione dell’equazione lineare
if a~=0
x=-b/a;
elseif b~=0
disp('equazione impossibile')
else
disp('equazione indeterminata') 20
end
 Esercizio

Chiamando fun_eqlin.m dal Command Window

>> x = fun_eqlin(10,2);
>> disp(x)
0.2000

>> x = fun_eqlin(0,2);
equazione impossibile

>> x = fun_eqlin(0,0);
equazione indeterminata

21
 Esercizio
Scrivere una function Matlab (fun_eq2grado) per il calcolo delle radici della
equazione di secondo grado ax2+bx+c=0
La funzione deve ricevere come input i tre coefficienti dell’equazione e
restituire come output le soluzioni

function [x] = fun_eq2grado(a,b,c)

% [x] = fun_eq2grado(a,b,c)
% Calcola le radici dell’ equazione ax^2+bx+c=0 e le memorizza nel vettore x
% INPUT
% a = primo coefficiente dell’equazione
% b = secondo coefficiente dell’equazione
% c = terzo coefficiente dell’equazione (termine noto)
% OUTPUT
% x = vettore delle soluzioni
if a==0,
disp(‘l’’equazione e’’ di primo grado’),
x(1) = fun_eqlin(b,c);
else, D = b^2-4*a*c;
if D < 0, disp(‘le radici non sono reali’),
elseif D==0, x(1) = -b/(2*a);
else, x(1) = (-b-sqrt(D))/(2*a);, x(2) = (-b+sqrt(D))/(2*a);, end 22
end
 Esercizio
Scrivere una function Matlab (fun_eq2grado) per il calcolo delle radici della
equazione di secondo grado ax2+bx+c=0
La funzione deve ricevere come input i tre coefficienti dell’equazione e
restituire come output le soluzioni

function [x] = fun_eq2grado(a,b,c)

% [x] = fun_eq2grado(a,b,c)
% Calcola le radici dell’ equazione ax^2+bx+c=0 e le memorizza nel vettore x
% INPUT
% a = primo coefficiente dell’equazione
% b = secondo coefficiente dell’equazione
% c = terzo coefficiente dell’equazione (termine noto)
% OUTPUT
% x = vettore delle soluzioni
if a==0,
disp(‘l’’equazione e’’ di primo grado’), è una funzione
x(1) = fun_eqlin(b,c); esterna
else, D = b^2-4*a*c;
if D < 0, disp(‘le radici non sono reali’),
elseif D==0, x(1) = -b/(2*a);
else, x(1) = (-b-sqrt(D))/(2*a);, x(2) = (-b+sqrt(D))/(2*a);, end 23
end
 Esercizio
Chiamando fun_eq2grado.m dal Command Window
>> x=fun_eq2grado(1,2,3);
>> disp(x)
-2 -1

>> x=fun_eq2grado(0,3,2);
l'equazione e' di primo grado
>> disp(x)
-0.6667

>> x=fun_eq2grado(1,2,1);
>> disp(x)
-1.00000

>> x=fun_eq2grado(1,1,1);
le radici non sono reali

24
 Function files
- Una funzione può contenere un’altra funzione (sotto-funzioni) mentre
uno script può solo richiamare funzioni memorizzate in altri file .m
Una sotto-funzione può essere richiamata solo dalla funzione
all’interno della quale è definita
function [mean,stdev]=stat(x) function [m] = calcola_media(x,n)
% Calcola la media (mean) e % calcola_media calcola il valor
% la deviazione standard (stdev) % medio degli elementi del vettore x
% di lunghezza n
% degli elementi del vettore x (input)
% INPUT
n = length(x); % x = vettore di numeri
mean = calcola_media(x,n); % n = lunghezza di x
stdev = calcola_stddev(x,n); %
% OUTPUT
% m = media degli elementi di x
function [st] = calcola_stddev(x,n) if length(x)=n,
mean = calcola_media(x,n); error(‘controllare la lunghezza di x’),
y = (x-mean).^2; end
m=sum(x)/n;
st = sqrt(calcola_media(y,n));
 Function files
File stat.m File calcola_media.m

function [mean,stdev]=stat(x) function [m] = calcola_media(x,n)

% Calcola la media (mean) e % calcola_media calcola il valor
% la deviazione standard (stdev) % medio degli elementi del vettore x
% di lunghezza n
% degli elementi del vettore x (input)
% INPUT
n = length(x); % x = vettore di numeri
mean = calcola_media(x,n); % n = lunghezza di x
%
stdev = calcola_stddev(x,n);
% OUTPUT
% m = media degli elementi di x
function [st] = calcola_stddev(x,n) if length(x)=n,
mean = calcola_media(x,n); error(‘controllare la lunghezza di x’),
end
y = (x-mean).^2;
m=sum(x)/n;
st = sqrt(calcola_media(y,n));
 Function files
File stat.m

function [mean,stdev]=stat(x)
% Calcola la media (mean) e Sottofunzione della funzione stat.m
% la deviazione standard (stdev) • calcola_stdev non è un file .m,
% degli elementi del vettore x (input) non compare nella lista dei file
n = length(x); .m della directory di lavoro,
mean = calcola_media(x,n); • può essere usata solo all’interno
di stat.m
stdev = calcola_stddev(x,n);
non è riconosciuta all’esterno

function [st] = calcola_stddev(x,n)

mean = calcola_media(x,n);
y = (x-mean).^2;
st = sqrt(calcola_media(y,n));
 Esercizio

Riscrivere la funzione matlab fun_eq2grado.m per la soluzione di

un’equazione di secondo grado definendo come sua sottofunzione
la funzione per la soluzione di un’equazione lineare di primo grado

28
 Esercizio

Scrivere le funzioni matlab bisezione_fun, punto_unito_fun,

newton_fun e secanti_fun per la soluzione di un’equazione non
lineare rispettivamente con i metodi di bisezione, del punto unito,
di Newton-Raphson e delle secanti usando opportune variabili di
input e output

29
 Esercizio
function [xn,err1,err2,iter] = bisezione_fun (f,a,b,eps)
% [xn,err1,err2,iter] = bisezione_fun (f,a,b,eps)
% cerca la radice della funzione f nell'intervallo [a,b] con precisione eps
% utilizzando un doppio criterio di arresto
%
% INPUT
% f = espressione della funzione di cui si vuole cercare la radice
% a = estremo inferiore dell'intervallo in cui è stata isolata la radice
% b = estremo superiore dell'intervallo in cui è stata isolata la radice
% eps = limite superiore dell'errore da usare come criterio di arresto
%
% OUTPUT
% xn = approssimazione della radice
% err1 = |xn-x(n-1)| differenza in valore assoluto tra due approssimazioni
% successive
% err2 = f(xn) valore della funzione nell'approssimazione xn
% iter = numero di iterazioni eseguite
30
 Esercizio
format long;
%inizializzazione dei parametri
iter = 0; err1 = b-a; err2 = 10^10;, x0=a;

while (err1>eps) & (err2>eps) % usare | se si vuole che la approssimazione

% prodotta verifichi le due condizioni
% contemporaneamente
xn = (a+b)/2;
if ( f(a)*f(xn) < 0 )
b = xn;
elseif ( f(xn)*f(b) < 0 )
a = xn;
end
iter = iter+1;
err1 = abs(xn-x0);
err2 = abs(f(xn));
x0 = xn;
end 31
 Esercizio
function [xn,err,n_iter] = punto_unito_fun(f,x0,eps,maxiter)
% [xn,err,n_iter] = punto_unito_fun(f,x0,eps,maxiter)
% cerca il punto unito della funzione f con precisione eps
% scegliendo x0 come punto iniziale.
% Il procedimento iterativo si interrompe se si raggiunge il
% numero massimo di iterazioni maxiter
%
% INPUT
% f = espressione della funzione di cui si vuole cercare il punto unito
% x0 = punto iniziale
% eps = limite superiore dell'errore da usare come criterio di arresto
% maxiter = numero massimo di iterazioni consentite
% OUTPUT
% xn = approssimazione prodotta dal metodo del punto unito
% err = |xn-x(n-1)|
% n_iter = numero di iterazioni eseguite
32
 Esercizio
format long;

%inizializzazione dei parametri

n_iter = 0; err = 10;

while (err>eps) & (n_iter<maxiter)

n_iter = n_iter+1;
xn = f(x0);
err = abs(xn-x0);
x0 = xn;
end

if n_iter >= maxiter

fprintf('La soluzione prodotta non ha raggiunto la precisione richiesta in
%5d iterazioni',maxiter)
end

33
 Esercizio

Approssimazione Errore Numero di iterazioni eseguite

34
Esercizio

35
36
37
Esercizio

38
Esercizio 1

39
40
41
42
43
Esercizio 2

44
45
46
47
48
49
50
Esercizi

51
Esercizi

52
Esercizio 1

53
54
55
56
57
58
Esercizio2

59
60
61
Dal Command Window

62
63
64
65
66
67
68
69
Esercizio

70
 Esercizio
Per risolvere l’esercizio precedente è necessario prima risolvere il
seguente

71
Esercizio

72
73
74
75
76
77
78
79
80
81
Esercizio

82
83