# Calcolo parametri acustici del sonar

Il calcolo dei parametri acustici del sonar è stato propedeutico alle fasi di progettazione dei sonar per i sottomarini classi Toti e Sauro; per le basi elettroacustiche (sistemi di trasduttori) e per le componenti elettroniche (sistemi di elaborazione segnali e gruppi accessori).
Il calcolo dei parametri acustici veniva fatto, nel lontano 1960, utilizzando strumenti quali la "Pascalina" (calcolatore meccanico a rotazione manuale), il regolo calcolatore a stecca, tecnigrafo (per l'esecuzione delle sommatorie di vettori necessarie per il calcolo della direttività delle basi idrofoniche) e l'indispensabile volume "Handbook of Mathematical Functions", con un pesante impegno di tempo, sono oggi rapidamente fattibili con il calcolatore software senza la necessità dello sviluppo di algoritmi non sempre facilmente manipolabili.
È da tempo che i metodi accennati sono stati sostituiti da sistemi di calcolo su PC, veloci ed affidabili, secondo quanto di seguito illustrato.

## Sistema di calcolo su P.C.

Il sistema di calcolo è sviluppato in un programma per PC che, una volta lanciato, presenta la pagina di selezione delle funzioni di calcolo mostrata in figura con la quale scegliere il tema voluto.
Le funzioni sono nell'ordine:
* Livello spettrale del rumore del mare in funzione della frequenza e dello stato dell'ambiente
* Livello spettrale del rumore emesso da un C.T. in funzione della frequenza e della velocità del mezzo
* Livello dell'attenuazione per propagazione in funzione della frequenza e della distanza
* Caratteristica di direttività e guadagno di un trasduttore circolare in funzione della frequenza e del diametro
* Caratteristica di direttività e guadagno di una base lineare in funzione della frequenza e della lunghezza
* Distanza di un bersaglio mediante riflessione dal fondo in funzione dell'angolo di depressione e della profondità
* Guadagno di direttività di una sorgente rettangolare/quadrata in funzione della frequenza e delle dimensioni
* Portata di scoperta di un sonar passivo in funzione di: (Banda-SL-NL-DI-RC-d-Propag.)
* Portata di scoperta di un sonar attivo in funzione di: (Ftx-BWrx-SL-TS-NL-DI-ti-d-Propag.)

## Il file di calcolo

Il programma di calcolo può essere scaricato liberamente, su uno dei due link:
* Calcolatore parametri Sonar, sonar-info.info.
* Calcolatore parametri Sonar (ZIP), github.com.

## Procedure di calcolo

Un insieme di esercizi per il calcolo di nove funzioni relative ai sistemi sonar, attivi e passivi è illustrato confrontando gli algoritmi classici con gli analoghi processi automatici

### Rumore del mare
#### Algoritmi classici

Con riferimento alla voce rumore del mare si riporta l'algoritmo di calcolo:
Il rumore del mare, indicato con la sigla $NL$, è variabile in ampiezza in funzione della frequenza e dello stato del mare secondo l'espressione approssimata:
$NL=k-20\cdot \log _{10}{(f^{5/6})}$.
Determinazione di k :
| SS | 0 | 1/2 | 1 | 2 | 4 | 6 |
| k | 44 | 50 | 55 | 62 | 66 | 71 |

#### Calcolo automatico

Con il calcolatore il valore di $(NL)$ è di acquisizione immediata; la schermata di lavoro per la valutazione del rumore del mare è mostrata in figura:
Calcolo del livello spettrale $(NL)$ del rumore del mare nelle seguenti condizioni:
* Frequenza operativa $fo=3\ kHz$
* Stato del mare $SS=2$
Il valore di $fo$ si digita, in $kHz$ , nell'apposita casella.
Per inserire lo stato del mare $(SS)$ si deve cliccare sul cursore della finestra "Stato del mare" e selezionare il valore; nel nostro caso $SS=2$
Cliccando su "Calcolo" nella finestra a fianco compare il livello calcolato di rumore spettrale $NL=53.9\ dB/\mu Pa/{\sqrt {Hz}}$
Si deve osservare che i valori di $SS$ non possono essere scelti al di fuori dei dati della lista.
Per valutare con approssimazione un valore non in lista si deve cliccare "Grafici" ottenendo una nuova schermata con i grafici a colori dei $6$ stati del mare predefiniti, tra i quali poter valutare il valore di $NL$ più vicino al dato di $SS$ fuori tabella.
I valori dei grafici sono estrapolati dal regolo Sonar Performance Calculator della Raytheon.
Osservazione:
valori di $SS;NL$ sono orientativi e pertanto non richiedono calcoli di precisione.

### Livello spettrale del rumore di un semovente
#### Algoritmi classici

Con riferimento alla voce rumori dei semoventi navali l'algoritmo di calcolo di $SL$:
$SL=k-20\cdot \log _{10}{(f^{1.108})}$
$f$: frequenza in $kHz$
$k$: coefficiente dipendente dalla velocità della nave i nodi ($kn$) secondo la tabella:
| Kn | 10 | 15 | 20 | 25 |
| k | 127.5 | 138 | 145.4 | 151.4 |

#### Calcolo automatico

Con il calcolatore il valore di $SL$ è di acquisizione immediata:
Calcolo del livello spettrale $(SL)$ del rumore irradiato da un cacciatorpediniere (CT) nelle seguenti condizioni:
* Frequenza operativa $fo=10\ kHz$
* Velocità CT: $V=10\ nodi$
Il valore di $fo$ si digita, in $kHz$, nell'apposita casella.
Per inserire la velocità del CT si deve cliccare sul cursore della finestra "Velocità del CT" e selezionare il valore; nel nostro caso $V=10\ nodi$
Cliccando su "Calcolo" nella finestra a fianco compare il livello calcolato di rumore spettrale $SL=105.3\ dB/\mu Pa/{\sqrt {Hz}}/1m$
Si deve osservare che i valori di SS non possono essere scelti al di fuori dei dati della lista.
Per valutare con approssimazione un valore non in lista si deve cliccare "Grafici" ottenendo una nuova schermata con i grafici a colori del $SL$ per $4$ valori di velocità predefiniti, tra i quali poter valutare il valore di $SL$ più vicino al dato di $V$ fuori tabella.
I valori dei grafici sono estrapolati dal regolo "Sonar Performance Calculator" della Raytheon.
Oervazioni:
* I valori di $V\ e\ SL$ sono orientativi e non richiedono pertanto calcoli di precisione.
* I livelli di rumore emessi dal CT son relativi ad unità operative negli anni 40-45.

### Attenuazione del suono in mare
#### Algoritmi classici

Con riferimento alla voce Trasmissione del suono in mare si riportano gli algoritmi di calcolo:
$TL=60+20\cdot \log _{10}{R}$ (attenuazione per divergenza sferica)
$TL=60+10\cdot \log _{10}{R}$ (attenuazione per divergenza sferico - cilindrica)
$\alpha =\left[{\frac {0.1\cdot f^{2}}{1+f^{2}}}\right]+\left[{\frac {40\cdot f^{2}}{4100+f^{2}}}\right]+\left[{\frac {2.75\cdot f^{2}}{10^{4}}}\right]$ (attenuazione per assorbimento)
dove:
$\alpha \ in\ dB/km$
$f\ in\ kHz$

#### Calcolo automatico

Calcolo dell'attenuazione del suono, $(TL)$ , per propagazione nelle seguenti condizioni:
* Frequenza operativa $fo=7\ kHz$
* Distanza dalla sorgente $Dist.=15\ km$
* Tipo di propagazione: Sferica
Il valore di $fo=7\ kHz$ si digita nell'apposita casella, così il valore della distanza $D=15\ km$
Per inserire il tipo di propagazione si deve cliccare sul cursore della finestra e selezionare "Sferica"
Cliccando su "Calcolo" nelle finestre a destra compaiono sia il valore di $TL$ dovuto alla divergenza: $TL=\ 83.5\ dB$, sia l'attenuazione dovuta all'assorbimento $TL=8.8\ dB$
L'attenuazione che dovrà essere considerata per le fasi di progettazione sarà la somma dei due valori $83.5+8.8=92.3\ dB$; questi valori sono per unico percorso.
Ricordare che in questo tipo di computazioni le frazioni di decibel non vengono considerate ai fini progettuali. Cliccando successivamente su "grafici" si possono osservare le curve generali di attenuazione teoriche:
* $TL$ per propagazione sferica ,in rosso
* $TL$ per propagazione sferico/cilindrica, in blu
* $TL$ per assorbimento,in nero

### Direttività trasduttore circolare
#### Algoritmi classici

Non sono disponibili voci su questo tema.

#### Calcolo automatico

Direttività di un trasduttore circolare con le seguenti caratteristiche:
* Frequenza operativa $fo=10\ kHz$
* Diametro $d=0.5\ m$
Il valore di $fo=10\ kHz$ si digita nell'apposita casella, così il valore del diametro $d=0.5\ m.$
Cliccando su "Calcolo" viene presentata, in doppia scala lineare, la curva di direttività in funzione dell'angolo (in figura è tracciata soltanto una metà della curva) ; nelle finestre in basso a destra compaiono:
* Larghezza del lobo a $-3\ dB:\alpha =19$°
* Larghezza del lobo a $-6\ dB:\alpha =25.4$°
* Guadagno di direttività $DI=20.2\ dB$

### Direttività base acustica lineare
#### Algoritmi classici

Non sono disponibili voci su questo tema.

#### Calcolo automatico

Calcolo della direttività di una base lineare con le seguenti caratteristiche:
* Frequenza operativa $fo=1\ kHz$
* Lunghezza $L=10\ m$
Il valore di $fo=1\ kHz$ si digita nell'apposita casella, così il valore della lunghezza $L=10\ m$
Cliccando su "Calcolo" viene presentata, in doppia scala lineare, la curva di direttività in funzione dell'angolo (in figura è tracciata soltanto una metà della curva.) ; nelle finestre in basso a destra compaiono:
* Larghezza del lobo a $-3\ dB:\alpha =8$°
* Larghezza del lobo a $-6\ dB:\alpha =10.8$°
* Guadagno di direttività $DI=11.2\ dB$

### Misura della distanza del bersaglio
#### Algoritmi classici

Con riferimento alla voce Misura della distanza con il sonar si riportano gli algoritmi di calcolo
$R={\frac {2\cdot P}{\tan \beta }}$
dove:
$R=$ distanza del bersaglio
$P=$ quota del sottomarino
$\beta =$ angolo di depressione
Se i due mezzi si trovano a quote diverse; uno alla quota $P$ e l'altro alla quota $P-dq$ il calcolo di $R$ si esegue con la formula:
$R={\sqrt[]{[({\frac {2\cdot P-dq}{\tan \beta }})^{2}+dq^{2}]}}$

#### Calcolo automatico

Calcolo della distanza di un bersaglio attivo nelle seguenti condizioni operative :
* Profondità del fondo $h=180\ m$
* Angolo di depressione misurato $\beta =6$°
Il valore di $h=180\ m$ si digita nell'apposita casella, così il valore dell'angolo di depressione $\beta =6$°
Cliccando su "Calcolo" viene presentata, in doppia scala lineare, la curva generale della variazione della distanza $d$ in funzione dell' angolo di depressione per la profondità di $h=180\ m.$ impostati; nella finestra in centro compare la distanza in metri relativa alla situazione ipotizzata:
$d=3425\ m.$

### Guadagno base idrofonica rettangolare
#### Algoritmi classici

Non sono disponibili voci su questo tema.

#### Calcolo automatico

Calcolo del guadagno di una sorgente rettangolare (trasduttore) con le seguenti caratteristiche :
* Frequenza di lavoro $fo=12\ kHz$
* Dimensioni $50\ cm$ x $25\ cm$
Inserite le variabili cliccando su "Calcolo", nella finestra apposita, viene visualizzato il guadagno del trasduttore:
$DI=19.8\ dB$

### Portata di scoperta del sonar passivo
#### Algoritmi classici

Con riferimento alla voce Portata sonar passivo si riporta l'algoritmo di calcolo: ${\begin{cases}TL=60+20\cdot \log _{10}{R}+\alpha \cdot R\\TL=SL+DI-NL-DT+10\cdot \log _{10}{BW}\end{cases}}$

#### Calcolo automatico

Calcolo della distanza di scoperta teorica di un sonar passivo, in condizioni di propagazione normale, le cui variabili acustiche sono di seguito indicate:
* Banda di ricezione $BW=2000\ Hz-34000\ Hz$
* Livello spettrale emesso dal bersaglio $SL=100\ dB/\mu Pa/{\sqrt {Hz}}/1m$
* Rumore spettrale del mare $NL=23\ dB/\mu Pa/{\sqrt {Hz}}$
* Guadagno della base ricevente $DI=21\ dB$
* Costante di tempo dell'integratore $RC=1\ S.$
* Parametro $d=10$ dipende dalla probabilità di scoperta Priv. e dalla probabilità di falso allarme Pfa.
* Propagazione ipotizzata $Sferica$
Le variabili sopra indicate devono essere digitate nelle caselle di "IMPOSTAZIONE DATI".
Cliccando il pulsante "Calcolo" si risolve il problema voluto che indica in $35\ km$ la distanza di scoperta del bersaglio.

### Portata di scoperta di un sonar attivo
#### Algoritmi classici

Con riferimento alla voce Portata sonar attivo si riporta l'algoritmo di calcolo della portata per propagazione sferica per la componente attiva del sonar si ottiene dalla soluzione del sistema trascendente in $R$ dove tutte le variabili sono espresse in decibel:
${\begin{cases}TL=120+40\cdot \log _{10}{R}+2\cdot \alpha \cdot R\\TL=LI+TS-NL-DT+DI\end{cases}}$

#### Calcolo automatico

Calcolo della portata di scoperta teorica di un sonar attivo, in condizioni di propagazione normale, le cui variabili acustiche sono di seguito indicate:
* Frequenza di emissione $Ftx=9000\ Hz$
* Banda di ricezione $Bw=500\ Hz$
* Livello di emissione $SL=240\ dB/\mu Pa$
* Rumore spettrale del mare $NL=45\ dB/\mu Pa/{\sqrt {Hz}}$
* Guadagno della base ricevente $DI=29dB$
* Costante di tempo dell'integratore $RC=0.01\ S.$
* Parametro $d=7$ dipende dalla probabilità di scoperta Priv. e dalla probabilità di falso allarme Pfa.
* Forza del bersaglio $TS=10\ dB$
* Propagazione ipotizzata = $Sferico/Cilindrica$
Le variabili devono essere digitate nelle caselle di "IMPOSTAZIONE DATI"; cliccando il pulsante "Calcolo" si risolve il problema voluto che indica in $31.8\$ la distanza di scoperta del bersaglio..