3 a lezione - laboratorio a.a 2004-2005 Corso di Laurea Ingegneria MECCANICA.

Post on 01-May-2015

217 views 0 download

Transcript of 3 a lezione - laboratorio a.a 2004-2005 Corso di Laurea Ingegneria MECCANICA.

33aa lezione - laboratorio lezione - laboratorio

a.a 2004-2005

Corso di Laurea Ingegneria MECCANICA

Esercizio1 - parte aEsercizio1 - parte a

Inserire nel file ciclo_if.m considerato nella lezione 2, un comando di input per gestire in maniera interattiva la dimensione della matrice; usare poi il comando disp per visualizzare il risultato.Costruire le matrici di dimensioni N=5 e N=4 (nell’ordine) ed analizzare i risultati.

Soluzione esercizio1- parte aSoluzione esercizio1- parte a% File ciclo_if.m

N=input('Inserisci N ');

for I=1:N

for J=1:N

if I == J

A(I,J) = 2;

elseif abs(I-J)== 1

A(I,J) = -1;

else

A(I,J) = 0;

end

end

end

disp('A = ');disp(A)

Output file Output file ciclo_if.mciclo_if.m>> clear all

>> ciclo_if

Inserisci N 5

A =

2 -1 0 0 0

-1 2 -1 0 0

0 -1 2 -1 0

0 0 -1 2 -1

0 0 0 -1 2

>> ciclo_if

Inserisci N 4

A =

2 -1 0 0 0

-1 2 -1 0 0

0 -1 2 -1 0

0 0 -1 2 -1

0 0 0 -1 2

?????• Ad ogni chiamata viene aggiornata nei termini, ma non

nella dimensione, la matrice A contenuta nel workspace• Il nome della variabile di output non può essere modificato

senza modificare il file

Esercizio1 – parte bEsercizio1 – parte b

Creare un file function, dal nome fun_ciclo_if.m, che contenga tutte le istruzioni del file ciclo_if.m, tranne la prima (N=input(‘Inserisci N’)).

Utilizzare il file fun_ciclo_if.m per costruire le matrici con dimensioni N=5 e N=4 (nell’ordine) ed analizzare i risultati.

Soluzione esercizio1 – parte bSoluzione esercizio1 – parte b

function A=fun_ciclo_if(N)

for I=1:N

for J=1:N

if I == J

A(I,J) = 2;

elseif abs(I-J)== 1

A(I,J) = -1;

else

A(I,J) = 0;

end

end

end

Output file Output file fun_ciclo_if.mfun_ciclo_if.m>> clear all

>> A=fun_ciclo_if(5);

A =

2 -1 0 0 0

-1 2 -1 0 0

0 -1 2 -1 0

0 0 -1 2 -1

0 0 0 -1 2

>> A=fun_ciclo_if(4);

A =

2 -1 0 0

-1 2 -1 0

0 -1 2 -1

0 0 -1 2

>>

• Ad ogni chiamata la matrice A, contenuta nel workspace, viene sostituita

• Si possono avere tutte e due le matrici dando un nome diverso alla variabile di output nella chiamata della funzione

Esercizio 2 Esercizio 2

Scrivere un file function che assegni a y la seguente espressione:

y=(x/2.4)^3-2*x/2.4+cos(pi*x/2.4);

Si calcoli, quindi, il valore di y per:

• x=2.5,

• x=[2.3 4.5 6.2].

File function e risultati esercizio 2 File function e risultati esercizio 2

>> y=p(2.5)y = -1.9445>> x=[2.3 4.5 6.2];>> y=p(x)??? Error using ==> ^Matrix must be square.Error in ==> C:\MATLAB6p1\work\p.mOn line 2 ==> y=(x/2.4)^3-2*x/2.4+cos(pi*x/2.4);

function y =p(x)

y=(x/2.4)^3-2*x/2.4+ …cos(pi* x/2.4);

Correzione file esercizio 2 Correzione file esercizio 2

>> y=p(x)

y =

-2.0280 3.7657 11.8147

>>

>> y=p(x)

y =

-2.0280 3.7657 11.8147

>>

function y =p(x)y=(x/2.4).^3-2*x/2.4+cos(pi*x/2.4);

Altro modo di risolvere l’esercizio 2 Altro modo di risolvere l’esercizio 2

>> y=p(x)

y =

-2.0280 3.7657 11.8147

>>

>> y=p(x)

y =

-2.0280 3.7657 11.8147

>>

function y =p(x)f='(x/2.4).^3-2*x/2.4+cos(pi*x/2.4)';y=eval(f);

Esercizio 3Esercizio 3

A- Scrivere un file script che consenta di disegnare, in [0,4] le funzioni:

y=3*sin(pi*x) e y= exp(-0.2*x) sulla stessa figura. Si consideri la

partizione x = 0:0.02:4B- Inserire le label per gli assi x, y ed il titolo.

C- Usare gtext per indicare i vari punti di intersezione dei grafici.

D- Memorizzare il file col nome grafico.

Come creare un grafico 2-DCome creare un grafico 2-D

Sintassi per disegnare una curva con: • specifica dei dati nel vettore x e/o y• specifica del colore, e dello stile della

linea

plot(x, y, ‘colore_stilelinea’)plot(x, y, ‘colore_stilelinea’)

plot(x1, y1,’r:’,x2,y2,‘b-’,...)plot(x1, y1,’r:’,x2,y2,‘b-’,...)

• Sintassi per disegnare più curve:

Risultato esercizio3 Risultato esercizio3 titolo

» title('Grafico ed …')

etichetta asse y

» ylabel('y')

etichetta asse x

» xlabel('x')

gtext

»gtext('x3')

File File grafico.mgrafico.m

clear allx=0:0.02:4;y=3*sin(pi*x);plot(x,y,'r'),xlabel('x');ylabel('y');gridhold on y=exp(-.2*x);plot(x,y,'g') %osservare il numero delle intersezionigtext('x1');gtext('x2');gtext('x3');gtext('x4');

title('Grafico ed intersezioni di 3*sin(pi*x) e exp(-0.2*x)')

hold off

clear allx=0:0.02:4;y=3*sin(pi*x);plot(x,y,'r'),xlabel('x');ylabel('y');gridhold on y=exp(-.2*x);plot(x,y,'g') %osservare il numero delle intersezionigtext('x1');gtext('x2');gtext('x3');gtext('x4');

title('Grafico ed intersezioni di 3*sin(pi*x) e exp(-0.2*x)')

hold off

Come operare con l’editor graficoCome operare con l’editor grafico

Accedere alla finestra “property editor” per migliorare il grafico, selezionando, con il mouse, la voce Edit sulla figura e, dal menu che compare, la voce Figure Properties…

Le icone indicate dalle linee, consentono di inserire testo, frecce e linee su una figura

Come inserire un testo sulla figuraCome inserire un testo sulla figura1. Selezionare

l’icona2. Posizionarsi nel

punto desiderato e cliccare

3. Scrivere x_3 per ottenere x3 (opp. x^3 per x3)

4. cliccare fuori dal riquadro per rendere attivo lo scritto

Come Come inserire una freccia sulla figura

1. Selezionare l’icona

2. Posizionarsi nel punto di inizio della freccia, cliccare il tasto sinistro del mouse che si deve trascinare fino al punto di arrivo della freccia.

Esercizio 4Esercizio 4Scrivere un file script che consenta di disegnare, in [-2,2] e su due “figure” distinte, il grafico della funzione: f(x)=exp(-x2)cos(20x) che viene definita nella function fun, utilizzando i comandi MATLAB:• plot per la figura 1 • fplot per la figura 2.

N.B. Si può utilizzare la partizione x=[-2:0.1:2]

Soluzione esercizio 4Soluzione esercizio 4function y=fun(x)y=exp(-x.^2).*cos(20*x);

figure(1)x=[-2:0.1:2];y=fun(x); plot(x,y),title('Comando plot')xlabel('x');ylabel('y');gridfigure(2)I=[-2,2];fplot('fun',I), gridtitle('Comando fplot')xlabel('x');ylabel('y')

Figura 1Figura 1

Figura 2Figura 2

Comando subplotComando subplot%file figurex=0.1:.1:5;subplot(2,3,1);plot(x,x);title('y=x');xlabel('x'); ylabel('y');subplot(2,3,2);plot(x,x.^2);title('y= x^2');xlabel('x'); ylabel('y');subplot(2,3,3),plot(x,x.^3);title('y= x^3');xlabel('x'); ylabel('y');subplot(2,3,4),plot(x,cos(x));title('y=cos(x)');xlabel('x');ylabel('y');subplot(2,3,5),plot(x,cos(2*x));title('y=cos(2x)');xlabel('x');ylabel('y')subplot(2,3,6),plot(x,cos(3*x));title('y=cos(3x)');xlabel('x');ylabel('y')% I colori e lo spessore sono stati aggiunti utilizzando “figure properties”

Risultati file figureRisultati file figure

Per utilizzare i file functionPer utilizzare i file function

1. Salvare le directory: Sis_lin, eq_non_lin,Approx, ecc. … nella directory

C:\MATLAB6p1\work2. Salvare il file che si sta creando in editing e che magari contiene i dati per l’input della function che si vuole utilizzare, nella directory: C:\MATLAB6p1\work3. Mettere in Current Directory il nome della directory che contiene le function Matlab che si vogliono utilizzare ad esempio C:\MATLAB6p1\work\eq_non_lin Buon lavoro!!!!