Il rumore nel sonar e la curva gaussiana


Il rumore che agisce contro la scoperta sonar è presente nel mezzo a causa dello stato del mare e altre cause; il rumore dell'acqua in superficie dipendente dal vento e/o dalla pioggia, il riverbero e il rumore dalle macchine proprie della nave penalizzano il segnale a volte in modo molto pesante.

lezione
lezione
Il rumore nel sonar e la curva gaussiana
Tipo di risorsa Tipo: lezione
Materia di appartenenza Materia: Analisi dei disturbi nell'ambiente marino
Avanzamento Avanzamento: lezione completa al 100%

Il tema in oggetto è stato esaminato nella 2^ lezione della materia di Analisi dei disturbi nell'ambiente marino lezione della quale la presente ne è un complemento.

L'algoritmo di Gauss è proposto di seguito in modo analitico cercando di fare chiarezza per chi si occupa di sonar e non di statistica.

La curva gaussiana

modifica

E' definita con il termine di cui al titolo una particolare funzione matematica che è in grado di rappresentare il comportamento di un notevolissimo numero di variabili naturali.

Tra le variabili citate anche le tensioni di rumore che sono convogliate, tramite i trasduttori, nei circuiti di elaborazione del sonar.

Dato che la caratteristica della tensione di rumore generata da un sensore acustico immerso in mare è tale che il suo livello di picco istantaneo può essere, casualmente, superiore a   volte il proprio valore efficace o inferiore a   volte tale valore, la curva di Gauss indica qual è la probabilità percentuale (  ) che ciò accada.

Con   si avrà il   di probabilità che ciò si verifichi, con   il   ecc.

La conoscenza delle citate percentuali è utile per la regolazione della soglia di rivelazione del sonar dalla quale dipendono, sia le probabilità di scoperta del bersaglio, sia le probabilità di falso allarme.

L'algoritmo che definisce la curva di Gauss è il seguente:

 

le variabili che compaiono sono quattro e le denominazioni loro attribuite, in questo caso specifico, sono

strettamente attinenti alle tensioni di rumore che si trovano nei circuiti del sonar :

  = variabile dipendente che esprime la probabilità percentuale sopra menzionata.

  = variabile indipendente che corre lungo l'asse delle ascisse, da   espressa in valore efficace.

  = valore efficace della tensione di rumore in esame.

  = valore medio della tensione di rumore in esame.

La curva gaussiana semplificata

modifica

Alla semplificazione dell'algoritmo indicato in precedenza giocano a favore le caratteristiche del rumore; questo infatti ha il valor medio uguale a zero annullando la quarta variabile della formula che perciò

si presenta così:

 

Se assumiamo ora che il valore efficace del rumore sia pari ad   ( possono essere, indifferentemente, dai   ai  ), l'algoritmo si semplifica ulteriormente come sotto indicato:

 

Esercizio di calcolo

modifica

Implementando in un calcolatore l'algoritmo semplificato otteniamo la sottostante routine di calcolo che consente di eseguire rapidamente le verifiche della procedura:


y = 1 / Sqr(6.28) * Exp(-x ^ 2 / 2)  

Print y   

Se supponiamo ad esempio che il rumore in esame abbia l'ampiezza   si può calcolare qual è la percentuale di probabilità che un livello istantaneo di picco sia  


  • La probabilità che in una tensione di rumore di   alcuni picchi possano raggiungere livelli pari alla metà del suo valore efficace si computa:

per   si ha:   pari al   del tempo.


  • La probabilità che in una tensione di rumore di   alcuni picchi possano raggiungere livelli pari ad una volta del suo valore efficace si computa:

per   si ha:   pari al   del tempo.


  • La probabilità che in una tensione di rumore di   alcuni picchi possano raggiungere livelli pari a due volte del suo valore efficace si computa:

per   si ha:   pari al   del tempo.


  • La probabilità che in una tensione di rumore di   alcuni picchi possano raggiungere livelli pari a tre volte del suo valore efficace si computa:

per   si ha :   pari al   del tempo.


  • La probabilità che in una tensione di rumore di   alcuni picchi possano raggiungere livelli pari a quattro volte del suo valore efficace si computa:

per   si ha:   pari al   del tempo.

Il grafico della curva gaussiana

modifica

Di notevole interesse è tracciare completamente la curva gaussiana tra limiti di   non molto elevati per poter ancora apprezzare i valori di  .

Questo si può fare con la routine seguente accompagnata da un adatto reticolo grafico ed istruzione pset:


for x = -k1 to k2 step .001  
	
y = 1 / Sqr(6.28) * Exp(-x ^ 2 / 2)  
	
next x 

Un esempio di tale grafica è sotto riportato per gli estremi di   definiti tra  


 
Curva gaussiana; distribuzione normale


Bibliografia

modifica
  • James J. Faran Jr e Robert Hills Jr, Correlators for signal reception, in Office of Naval Research (contract n5 ori-76 project order x technical memorandum no. 27), Cambridge, Massachusetts, Acoustics Research Laboratory Division of Applied Science Harvard University, 1952.
  • M. Sheldon Ross Probabilità e statistica per l'ingegneria e le scienze, Apogeo Trento, 2003
  • M Stephen Stigler, Mathematical Statistics in the Early States,The Annals of Statistics v.6 n.2 pp. 239–265|, 1978