# Curve di correlazione

Le curve di correlazione sono generate dagli algoritmi delle relative funzioni, il loro studio è propedeutico allo sviluppo delle metodologie che consentono la scoperta sonar.
In questa trattazione le curve s'intendono sempre calcolate tra due segnali: sono indicate come curve di correlazione incrociata.
Gli algoritmi che generano le curve, studiati per la scoperta di segnali elettrici coperti dal disturbo nelle apparecchiature sonar, si mettono in pratica tramite dispositivi indicati come correlatori elettronici o software.
Con gli algoritmi di correlazione si riesce ad individuare, ad esempio, in mezzo al rumore di ampiezza $10$, un segnale elettrico di ampiezza $1.$
La caratteristica di discriminazione dei segnali da parte dei citati algoritmi prevede uno studio delle curve degli stessi in base al tipo dei segnali e/o dispositivi d'utilizzazione al fine di ottimizzarne le proprietà in base al tipo d'impiego di necessità.
Lo studio prevede il tracciamento delle curve delle funzioni di correlazione in dipendenza di molte variabili quali:
Banda delle frequenze dei segnali ($0-F\ \ o\ \ F1-F2$ )
Tipo dei segnali (analogici o digitali)
Rapporto tra le ampiezze dei segnali e dei disturbi
Costante di tempo d'integrazione $rc$
Tempo di ritardo tra i segnali da correlare.
Tracciando le curve di correlazione, e lavorando sulle variabili, si può ottenere il profilo desiderato della funzione di correlazione che meglio sa adatta al tipo d'impiego.

## Applicazioni grafico numeriche

Le applicazioni presuppongono che i segnali acustici da correlare, generati da un semovente navale (la sorgente), giungano ad un sistema ricevente di due sensori dai quali prelevare le rispettive tensioni elettriche dei segnali stessi.

### Curva di correlazione analogica C = f(t) banda 0-F

Si tratta di correlazione analogica $C=f(t)$, normalizzata, in banda $0-F$ con $F$ in $Hz$, tra due segnali con ritardo $tc\ in\ \mu s$. tracciata in un reticolo cartesiano con fondo scala delle ascisse pari a F.scala $\ in\ \mu s$.
L'algoritmo di calcolo della funzione è:
$C(t)=\left[{\frac {\sin \ \{2\cdot \pi \cdot F\cdot (t-tc)\}}{2\cdot \pi \cdot F\cdot (t-tc)}}\right]\ \ \$
Con questa serie di dati, ad esempio:
si ottiene il grafico (asse x = tempo) della funzione di correlazione riportato in figura:
La curva mostra il massimo di correlazione alla 12ª divisione delle ascisse corrispondente a $600\ \mu s$ con $C=+1$ e profilo tondeggiante secondo $\operatorname {sen} x/x$.
La larghezza del lobo a $-3\ dB$ è di $.31.8\ \mu s$.
L'ampiezza massima dei lobi secondari è di $.0.1$.

### Curva di correlazione analogica C = f(t) banda F1-F2

Si tratta di correlazione analogica $C=f(t)$, normalizzata, in banda $F1-F2$ con $F$ in $Hz$, tra due segnali con ritardo $tc\ in\ \mu s$. tracciata in un reticolo cartesiano con fondo scala delle ascisse pari a F.scala $\ in\ \mu s$.
L'algoritmo di calcolo della funzione è:
$C(t)=\left[{\frac {\sin \ (2\cdot \pi \cdot DF\cdot \{t-tc\})\cdot \cos \ (2\cdot \pi \cdot Fo\cdot \{t-tc\})}{(2\cdot \pi \cdot DF\cdot \{t-tc\})}}\right]\ \ \$
Dove:
$DF=(F2-F1)/2$
$Fo=(F1+F2)/2$
Con questa serie di dati ad esempio:
si ottiene il grafico della funzione di correlazione di figura:
La curva mostra il massimo di correlazione alla 4ª divisione delle ascisse corrispondente a $200\ \mu s$ con $C=+1$ e profilo tondeggiante secondo $\operatorname {sen} x/x$.
La larghezza del lobo a $-3\ dB$ è di $100\ \mu s$
L'ampiezza massima dei lobi secondari è di $0.01$

### Curva di correlazione analogica C(b) in banda 0-F

Si tratta di correlazione analogica $C=f(b)$ normalizzata, in banda $0-F$ con $F\ in\ Hz$, tra due segnali che colpiscono una base con una inclinazione Brq = (b°) in gradi, tracciata in un reticolo cartesiano con scala delle ascisse pari a Fondo scala (a°) in gradi; la lunghezza della Base d è espressa in metri $(m)$.
In questa sezione di calcolo la geometria del sistema ricevente invece del tempo prevede l'angolo $b$ di puntamento con il max atteso per l'angolo $b$°; in questo esercizio le ascisse non sono dimensionate in tempo ma in gradi sessagesimali.
L'algoritmo di calcolo della funzione è:
$C(t)=\left[{\frac {\sin \ \{2\cdot \pi \cdot F\cdot (tv-tb)\}}{2\cdot \pi \cdot F\cdot (tv-tb)}}\right]\ \ \$
dove:
$tv=m\operatorname {sen} b/1530$ (b = variabile indipendente)
$tb=m\operatorname {sen} b^{o}/1530$ (b° = direzione della sorgente)
Con questa serie di dati ad esempio:
otteniamo il grafico della funzione di correlazione:
La curva mostra il massimo di correlazione alla 3ª divisione delle ascisse corrispondente a $6$° con $C=+1$ e profilo tondeggiante secondo $\operatorname {sen} x/x.$
La larghezza del lobo a $-3\ dB$ è di $4$°
L'ampiezza massima dei lobi secondari è di $0.13$

### Curva di correlazione digitale C=f(t) in banda 0-F

Si tratta di correlazione digitale $C=f(t)$ normalizzata, in banda $0-F$ con $F$ in $Hz$, tra due segnali con ritardo $tc\ in\ \mu s$. tracciata in un reticolo cartesiano con fondo scala delle ascisse pari a "F.scala" $\ in\ \mu s$
L'algoritmo di calcolo è:
$C(t)=(2/\pi )\ \arcsin \left[{\frac {\sin \ \{2\cdot \pi \cdot F\cdot (t-tc)\}}{2\cdot \pi \cdot F\cdot (t-tc)}}\right]\ \ \$
Per queste variabili ad esempio:
si ottiene il grafico della funzione di correlazione:
La curva mostra il massimo di correlazione alla 8ª divisione delle ascisse corrispondente a 200 $\mu s$ con $C=+1$ e profilo a cuspide secondo $\arcsin \ x/x.$
La larghezza del lobo a $-3\ dB$ è di circa 5 $\mu s$.
L'ampiezza massima dei lobi secondari è di $0.09$

### Curva di correlazione digitale in banda F1-F2

Si tratta di correlazione digitale $C=f(t)$ normalizzata, in banda $F1-F2$ con $F$ in $Hz$, tra due segnali con ritardo $tc\ in\ \mu s$. tracciata in un reticolo cartesiano con fondo scala delle ascisse pari a "F.scala" $\ in\ \mu s$
L'algoritmo di calcolo è:
$C(t)=C(t)=(2/\pi )\ \arcsin \left[{\frac {\sin \ (2\cdot \pi \cdot DF\cdot \{t-tc\})\cdot \cos \ (2\cdot \pi \cdot Fo\cdot \{t-tc\})}{(2\cdot \pi \cdot DF\cdot \{t-tc\})}}\right]\ \ \$
Dove:
$DF=(F2-F1)/2$
$Fo=(F1+F2)/2$
Con i dati ad esempio:
si ottiene il grafico della funzione di correlazione:
La curva mostra il massimo di correlazione alla 15ª divisione delle ascisse corrispondente a 1500 $\mu s$. con $C=+1$ e profilo a cuspide secondo $\arcsin \ x/x.$
La larghezza del lobo a $-3\ dB$ è di $120\ \mu s$.
L'ampiezza massima dei lobi secondari è di $0.09$

### Curva di correlazione digitale C = f(t) in presenza dei due segnali e del rumore del mare, in banda F1-F2

In questo esercizio la funzione dipende,oltre che dal tempo, anche dal rapporto $s/n$ (rapporto tra segnale e disturbo espresso in decibel) e dalla costante di tempo $rc$ dell'integratore.
Il max è atteso al tempo $tc,$ l'ampiezza di questo dipende da $s/n$ , la varianza da $rc.$
L'algoritmo di calcolo della funzione è:
$C(t)=(2/\pi )\arcsin K\cdot \left\{{\frac {\sin \ [2\cdot \pi \cdot DF\cdot (t-tc)]\cdot \cos \ [2\cdot \pi \cdot Fo\cdot (t-tc)]}{[2\cdot \pi \cdot DF\cdot (t-tc)]}}\right\}\ \ \$
dove $K$ è una variabile dipendente dal rapporto tra l'ampiezza del segnale $si$ e l'ampiezza del disturbo $ni$:
$K=\left[{\frac {1}{1+(ni/si)^{2}}}\right]\ \ \$
$DF=(F2-F1)/2$
$Fo=(F1+F2)/2$
con i dati d'esempio:
otteniamo il grafico della funzione di correlazione:
Si osservi che l'ampiezza della funzione C, a seguito del rapporto $s/n=+4\ dB$ inserito a calcolo, si è ridotta da $1$ a circa $0.5$ e il suo profilo si è modificato da una cuspide ad un andamento tondeggiante, lo spessore della traccia è indicativo della varianza d'uscita dal correlatore.
La curva mostra il massimo di correlazione alla 10ª divisione delle ascisse corrispondente a $400\ \mu s$. con $C=+0.5$ e profilo secondo $\operatorname {sen} x/x$.
La larghezza del lobo a $-3\ dB$ è di $40\ \mu s$.
L'ampiezza massima dei lobi secondari è di $0.06$

#### Rilievo sperimentale

L'effetto dell'alterazione della funzione di correlazione a causa del rumore sul segnale è mostrato nella fotografia rilevata in laboratorio su di un correlatore digitale per le condizioni: $S/N=-14\ dB$

### Curva di correlazione digitale con trasformata di Hilbert HC=f(t), in banda F1-F2

Si tratta di correlazione digitale normalizzata con trasformata di Hilbert $HC=f(t)$, in banda $F1-F2$ in $Hz$, tra due segnali con ritardo $tc\ in\ \mu s$. tracciata in un reticolo cartesiano con fondo scala delle ascisse pari a F.scala $\ in\ \mu s$.
Questa funzione dipende dal tempo e presenta uno zero dove le altre funzioni presentano il max (trasf. di Hilbert). Lo zero è atteso al tempo $tc.$
$C(t)=(2/\pi )\arcsin K\cdot \left\{{\frac {\sin \ [2\cdot \pi \cdot DF\cdot (t-tc)]\cdot \sin \ [2\cdot \pi \cdot Fo\cdot (t-tc)]}{[2\cdot \pi \cdot DF\cdot (t-tc)]}}\right\}\ \ \$
$K=\left[{\frac {1}{1+(ni/si)^{2}}}\right]\ \ \$
$DF=(F2-F1)/2$
$Fo=(F1+F2)/2$
Con i dati d'esempio:
otteniamo il grafico della funzione di anticorrelazione:
La curva mostra il passaggio per lo zero di correlazione alla 10ª divisione delle ascisse corrispondente a 200 $\mu s$. con $C=0$.
La pendenza attorno all'ascissa $x=200\ \mu s$. è di $0.035/1\ \mu s$.

## Tracciabilità delle curve

Il computo delle funzioni di correlazione e il relativo tracciamento delle curve è stato fatto per via automatica, per un numero molto elevato di punti di calcolo, in modo da consentire una buona grafica di presentazione.
Se ciascuna delle curve venisse tracciata, punto dopo punto, applicando manualmente gli algoritmi delle funzioni di correlazione il compito sarebbe improbo data la complessità di detti algoritmi; il risultato vedrebbe, come in figura, curve indicate a valori discreti che non fornirebbero le informazioni necessarie per il loro impiego.
Il problema del calcolo e della tracciabilità delle curve, con migliaia di punti di calcolo, è risolto su Wikiversità con l'impiego di adatto programma di calcolo all'indirizzo Processo di calcolo automatico delle funzioni di correlazione - Wikiversità