Informatica per Scienze Geologiche LT a.a.2017-2018 Esercitazione Analisi dati Pozzo Trebiciano Parte B2 Parte B2 Slides preparate da Dott. Pastorutti Docente: Prof. Carla Braitenberg, Dipartimento Matematica e Geoscienze, Via Weiss 1, Università di Trieste E-mail: [email protected]Tel. 040 5582258
12
Embed
Laboratorio Informatica LT 2017 2018 B2 - units.it...(vhufl]lr olyhoor dftxd qhoo¶$elvvr gl 7uhelfldqr $qql gl rvvhuyd]lrqh xq fdpslrqh rjql rud 8qlwd¶ gl plvxud p vrsud lo olyhoor
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Informatica per Scienze Geologiche LTa.a.2017-2018
Esercitazione Analisi dati PozzoTrebiciano
Parte B2Parte B2
Slides preparate da Dott. Pastorutti
Docente: Prof. Carla Braitenberg, Dipartimento Matematica e Geoscienze, Via Weiss 1, Università di Trieste
Continua il programma: qui la parte che crea il grafico
figure;
plot(liv(:,1),liv(:,6))
title('Livello acqua nel pozzo Abisso Trebiciano');
xlabel('tempo in anni');
ylabel('livello (m)');ylabel('livello (m)');
Ricerca dei massimi localinella serie livello_AbissoTrebiciano
%% valori estremiload('livello_AbissoTrebiciano_Tdec_CORRETTO.mat');liv=livello_tdec;n=length(liv(:,1));[xmax, kmax] = max(liv(:,6));tmax=liv(kmax,2:4);disp(['max value at ' ,num2str(tmax),...
%% graficofigura1 = figure;assi1 = axes('Parent',figura1);plot(assi1,liv(:,1),liv(:,6))hold ontitle('Livello acqua nel pozzo Abisso Trebiciano');xlabel('tempo in anni');ylabel('livello (m)');
Ricerca dei massimi localinella serie livello_AbissoTrebiciano
...[xmax, kmax] = max(liv(:,6));...
• Massimo assoluto: possiamo usare xmax, kmax già calcolati prima
scatter(assi1,liv(kmax,1),xmax)
x y
NB: la figura ottenuta prima non va chiusa,
altrimenti:
x y
Ricerca dei massimi localinella serie livello_AbissoTrebiciano
function [M,I] = maxsort(in,n)M = zeros(1,n); I = M;for i=1:n
[M(i),I(i)] = max(in);in(I(i)) = 0;
endend
• cercare i primi n valori più grandi: definiamo una funzione
1234567 end
• Salvate la function in un file con lo stesso nome della chiamata, quindi maxsort.m
«Ottieni il valore M e l’indice I degli n valori più grandi di in.»
Nella riga 2 vengono creati M ed I, vettori lunghi n, composti di soli zeri.
Le righe da 3 a 6 appartengono a un ciclo for: ad ogni iterazione viene calcolato l’n-esimo massimo e il suo indice, quindi quell’elemento di in viene sostituito con zero.
7
Ricerca dei massimi localinella serie livello_AbissoTrebiciano