È vero però che è ben scomodo mettere i nostri propri file .m lì dentro, in una cartella distante dal desktop e già piena di file...si rischia di fare confusione. Il nostro sistemaoperativo stupido ci viene però in aiuto un'altra volta:
Fate click(aaargh) col tasto destro sul suddetto collegamento, scegliete proprietà, nella maschera che si apre digitate(oppure copia incolla, tutto rigorosamente col mouse)nel campo "DA" il percorso completo della cartella in cui bramate tenere i vostri elaborati. Io ad esempio tengo una cartella sul desktop , quindi nel mio caso il percorso è:"c:/documents and settings/giulio/desktop/scriptoctave".
Ora quando riaprirete Octave utilizzando questo collegamento egli riconscerà come home la directory che avete appena settato, quindi con un pratico

dal prompt di octave potrete aprire il file nomefile.m che si trova nella vostra pratica cartella sul desktop, o dove più vi aggrada.Inoltre, una volta terminata l'esecuzione del medesimo, potrete continuare a lavorare sulle variabili che sono state calcolate da Octave, oppure lanciare un altro script, senza bisogno di riavviare octave(che pure, almeno sul mio computer) non è immediato a partire.

Ho inoltre recentemente scoperto che è possibile incollare il testo nel prompt di octave, il che si rivela spesso utile e pratico. Per farlo non funzionano ctrl V , ctrl shift V,shift ins o qualunque altra combinazione vi venga in mente, in Windows XP si può fare(a quanto ne so io)solo nel seguente modo:cliccando col tasto destro sul pulsante della finestra di octave nella barra delle applicazioni, selezionando modifica e poi incolla.
È inoltre molto utile la caratteristica di Octave di ricordarsi i comandi inseriti dal terminale. Questo significa che se avete appena scritto "A=1*4-log(10)/(4-a+b^2)"e vi serve riscriverlo e/o modificarlo(e farlo rieseguire a Octave), potete premere la frecciainsu, e riappariranno così al prompt i comandi che avete inserito recentemente.Ancora più utile può essere scrivere le prime lettere del comando che state cercando e premere la freccia insu, in questo modo appariranno solo i comandi già digitati che iniziano con quelle lettere.

Vediamo ora nel dettaglio le funzioni che abbiamo usato in classe:

Il numero degli elementi contenuti in un array si chiama "lunghezza dell'array". Per conoscerlo è possibile usare la funzione "length(arg)", dove l'argomento "arg" èl'array di cui vogliamo conoscere la lunghezza:

La funzione sin(arg) restituisce il seno dell'argomento "arg", che è l'angolo espresso in radianti:

Esiste in matlab la variabile di sistema "pi" che vale PIgreco (3.141592653589793...), per cui:

Volendo calcolare al tempo 0.04 secondi il valore di una funzione sinusoidale di ampiezza 1 e frequenza 237:

Abbiamo visto più sopra che è possibile far eseguire con un solo comando la stessa operazione su più valori, facendola eseguire non su uno scalare, ma su un array:

Vediamo a cosa ci può servire questo per i nostri scopi:
abbiamo più volte avuto necessità di calcolare il valore di una funzione sinusoidale con determinate frequenza e ampiezza in un certo intervallo di tempo, per poi poterla elaborare in qualche maniera(sovracampionare, sottocampionare, sottoquantizzare, puslwitmodulare etc etc)
La prima cosa da stabilire è la frequenza di campionamento. L'inverso di questa è poi il sample-increment, ovvero l'intervallo di tempo che intercorre tra un campionee il successivo. Conoscendo il sample-increment e la durata che vogliamo che abbia il nostro campione possiamo crearci un vettore contenente tutti gli "istanti" in cuivogliamo calcolare il valore della sinusoide. Una volta creato questo array possiamo agevolmente calcolare i valori come abbiamo appena visto, basta decidere la frequenzache deve avere la nostra sinusoide:

Per verificare che la nostra sinusoide è effettivamente una sinusoide possiamo rappresentare graficamente i dati ottenuti, utilizzando la funzione"plot(x,y)"."x" è il vettore che contiene i valori delle ascisse a cui corrispondono sulle ordinate i valori di "y".
In pratica verra plottato un grafico i cui punti avranno le coordinate:
P1 ( x(1), y(1) )
P2 ( x(2), y(2) )
P3 ( x(3), y(3) )
ecc
ecc
ecc
Ovviamente la lunghezza dei due array dovrà essere la medesima, sennò si verificherà un errore.
Ad esempio:

produrrà il seguente grafico:

A questo punto uno si spara in bocca oppure riflette sul fatto che 220500 campioni non possono venire rappresentati molto bene sui 1024 pixel del vs schermo.
A questo punto quindi zoomiamo nella presunta sinusoide con la funzione axis(arg). L'argomento "arg" dev'essere un vettore.I primi due elementi di questo array saranno rispettivamente l'inizio e la fine del range che vogliamo visualizzare sull'asse x. Se aggiungiamo altri due valori,questi saranno l'inizio e la fine del range che vogliamo visualizzare sull'asse y.

Qui possiamo vedere che effettivamente c'è una sinusoide, e che siccome vengono visualizzati 5 periodi in 50 millisecondi, essa avrà un periodo di 10ms, ovvero unafrequenza di 100Hz.

È possibile decidere il colore e lo stile di un grafico aggiungendo un altro argomento dopo i due array. Questo argomento dovrà essere una stringa(quindi racchiuso tra' ' oppure " ") e contenere alcuni dei seguenti caratteri(ne elenco pochi, una lista completa è nell'help):

Esempio:

Avrete notato che ogni volta che chiamate la funzione plot, il grafico viene stampato in Octave in una finestra intitolata "figure 1.0" . Potrebbe essere comodo visualizzare contemporaneamente più grafici. Per fare ciò è necessario utilizzare la funzione figure(n) , dove l'argomento "n" è un numero intero. L'esecuzione di questo comando provoca l'apertura di una nuova finestra di grafico, intitolata "Figure n.0". (di default è figure(1))

Abbiamo visto che ad ogni chiamata della funzione plot, il grafico precedentemente presente nella figura attiva viene cancellato e sostituito con quello nuovo.
Per plottare più grafici su una stessa figura invece, è necessario utilizzare il comando "hold on".
Esempio:

Ricordatevi di digitare "hold off" quando non volete più che il grafico "persista". Consiglio di aggiungere "hold off" alla fine di ogni script in cui si utilizza "hold on"per evitare strani comportamenti quando viene rieseguito.
In alternativa, è possibile plottare più grafici nella stessa figura utilizzando plot con un maggior numero di argomenti.

Se non vengono specificate opzioni di formato il grafico della prima coppia di vettori verrà rapresentato in "-b", mentre quello della seconda coppia in "-g"

Come detto, la trattazione di plot qui è incompleta, e rimando all'help della funzione per approfondimenti.

L'espressione wavread richiede l'argomento "filename" che è il nome del file riff/wave che si vuole caricare all'interno di Octave. Valgono le considerazioni fatte sopra a riguardo della home directory di octave ecc. Se avete seguito le mie indicazioni, sarà sufficiente avere il file wave nella vostra comoda directory sul desktop(o dove piùvia aggrada) e digitare

Così facendo nel vettore y vengono caricati i valori dei campioni contenuti nel file "ooops.wav".Vi consiglio di lavorare con file monofonici, è molto più semplice e non è necessario introdurre le matrici.
Spesso è comodo il seguente utilizzo di wavread, che permette di assegnare alle variabili y, fc, bps rispettivamente i valori dei campioni, la frequenza di campionamentoe il numero di bit per sample. In questa maniera potremo ad esempio usare fc e bps nei nostri calcoli, qualora ci sia necessario.

Potremo fare poi su questo vettore di valori tutti i conti che ci interessano, alla fine produrremo probabilmente un vettore, magari chiamato out, di dati elaborati,che vorremo in qualche maniera salvare come un file audio per poterlo ascoltare.
Con la funzione wavwrite(y,fc,bps,filename) ciò è possibile. "y" dev'essere un vettore contenente i campioni, "fc" è la frequenza di campionamento, "bps" il numero di bit persample con cui si vuole salvare il file, filename è una stringa contenene il nome del file di destinazione.