Visualizza wikitesto di Analisi del comportamento dei rivelatori d'energia sul rumore

{{Risorsa
| tipo = lezione
| materia1 = Il riconoscimento dei bersagli idrofonici in mezzo al disturbo
| avanzamento = 100%
}}

Nell'impiego del rivelatore d'energia per la scoperta di segnali idrofonici si deve considerare che, sia questi, sia il rumore ambiente, provocano all'uscita del dispositivo gli stessi tipi contributi di potenza non differenziabili facilmente; scopo di questo studio è l'analisi degli effetti del rumore ambiente. 

==Analisi di un rivelatore d'energia con disturbo all'ingresso==

===Il disturbo e la curva di Gauss===

Nella rivelazione d'energia quando all'ingresso del rivelatore a guadagno  <math> K </math>   il rumore d'ingresso 
in banda <math> \Delta  f  \ \ \ \ </math> è <math> Vni = 1 \ Veff.</math>  all'uscita '''<math>U</math>''', dopo l'integratore, si ha: 

* una componente continua di ampiezza <math>N </math> (potenza del rumore) <math> N = K \cdot Vni / 0.4 </math>
* una componente alternata, sovrapposta ad <math>N </math>, di valore efficace (deviazione standard) <math>  \delta n = K \cdot Vni  / { \sqrt{( 4 \ RC \ \Delta  f)}} </math>

dove <math> RC </math> è la costante di tempo dell'integratore.

Se ad esempio assumiamo:

:<math> K = 1 \ </math>
:<math> Vni = 1 \ Veff.</math>
:<math> RC = 1 \ s </math>
:<math> \Delta  f = 4000 \ Hz \ \ \ </math>  abbiamo:
:<math> N = 2.5 \ Vcc </math>
:<math>  \delta n =  1 / { \sqrt{( 4 \cdot 1 \cdot \ 4000) }} =  8 \ mVeff</math> 

Il valore efficace della varianza <math> \delta n</math> presenta dei picchi che oscillano in modo casuale. 

Supponiamo ora che, conoscendo a priori il livello della componente del  disturbo <math>N = 2.5 \ Vcc, \ \ \ </math>  si tenti di compensarla <ref>Più è precisa sarà la compensazione minore è l'effetto dell'ondulazione dei picchi della varianza</ref> e che tale compensazione possa essere fatta, all'uscita dell'integratore, con un errore dell'ordine del <math>10 </math>% il livello  della componente del disturbo compensata sarà : <math> 0.1 \ N  = 0.25 \ Vcc </math> 

Il livello della componente continua potrà variare attorno a  <math> 0.1 \ N </math> secondo i picchi della varianza <math> \delta n = 8 \ mV eff </math>.

Da un punto di vista statistico il comportamento del rumore dopo compensazione è illustrato in figura 1:

{{clear}}[[File:1dtcdtc.jpg|thumb|left|400px|]]
{{clear}}

la curva mostra che valori superiori od inferiori a <math>  0.1 \ N, </math> dovuti alla variabilità di <math>  \delta n </math> hanno una probabilità di verificarsi  decrescente secondo l'andamento della curva di Gauss.

===Il criterio della soglia nella scoperta sonar===
Se consideriamo il rivelatore d'energia sopra menzionato come parte integrante di un sonar passivo dal quale ci si aspetta la scoperta di  segnali emessi da semoventi il livello d'uscita '''<math>U</math>''', in presenza del solo disturbo d'ingresso <math>Vni</math>, non deve   allertare l'operatore al sonar al di sopra di una certa probabilità di falso allarme <math> Pfa</math> stabilita a priori.

Fissare una  percentuale di probabilita di falso allarme, ad esempio <math>Pfa = 5 </math>%, significa interdire l'uscita <math>U</math> del rivelatore con una soglia <math>T</math> quando la varianza del disturbo supera il <math>5 </math>% del tempo.

=== I picchi della varianza e la soglia===

Per la ricerca della percentuala dei valori di <math>  0.1 \ N + \delta \ n </math> che possono creare falsi allarmi si deve fare ancora riferimento  alla figura 1 ponendovi il limite di soglia di rivelazione T e la variabile <math>Pfa</math> relativa alla probabilità di falso allarme accettato; la soglia và fissata, generalmente, in modo che la probabilità di falso allarme nel caso peggiore <ref>Nei calcoli di portata del sonar si assume, generalmente <math> Pfa = 10 \% </math></ref> sia <math>  Pfa \le 10</math> %.

Le condizioni statistiche relative alla nuova impostazione delle variabili sono illustrate in figura 2: 

{{clear}}[[File:gausdtc2.jpg|thumb|left|400px|]]
{{clear}}

Nella figura si vede che il rumore dal lato più sfavorevole ha una densità di probabilità centrata attorno a <math> 0.1 \ N </math> e la soglia <math> T</math> va posta nella posizione ( la più spostata a destra ) tale che l'area a destra della soglia stessa sia il <math> 10 </math>%

Per la determinazione del valore della soglia <math>T</math> si procede alla ricerca dei legami
tra le diverse variabili:

indicando con <math>q</math> la differenza tra il livello di <math> T \ e \ ( 0.1 \ N + \delta \ n )</math>   normalizzata su <math> \delta \ n </math> si ha:  

:<math> q = [  T - ( 0.1 \ N +  \delta \ n ) ] / \delta \ n </math>

Si ricorre ora alla formula classica che lega la <math> P(fa) \ a  \ q </math> secondo la funzione riportata su Lezione 2^ della materia '''Sistemi di calcolo automatico per il sonar''' :
:<math>p(FA)=\left(\frac{1}{2}\right)\left\{1-erf\left(\frac{q}{\sqrt{2} }\right)\right\}</math>

risolvibile velocemente in <math>q </math> per <math> P(fa) = 10 \%</math> mediante la sezione 1) del calcolatore di che indica per <math> P(fa) = 10 \% \ ;  q = 1.29. </math>

con questo valore di <math>q </math> e l'equazione:
:<math> 1.29 = [  T - ( 0.1 \ N + \delta \ n ) ] / \delta \ n </math> ; si calcola il valore di <math>T</math>:
:<math> T =   1.29 \ \delta n + ( 0.1 \ N + \delta \ n )  </math>

ponendo per <math>\delta \ n</math> il valore calcolato all'inizio <math> \delta \ n = 8 \ mV  </math> <ref> nel calcolo ricordare :<math> 8 \ mV = 8/1000 \ V </math> </ref> si ha:
:<math>T \approx \ 0.27 \ V .</math>

Con questo valore di soglia e con una compensazione entro <math> 0.1 \ N </math> la probabilità di falso allarme si attesta a <math> Pfa = 10 \% </math> così come era nell'assunto.

==Il disturbo nei rivelatori quadratici asserviti ai fasci preformati==

Se consideriamo le caratteristiche del rivelatore e dell'impostazione della soglia <math>T</math>, valutate in precedenza, possiamo pensarle facenti parte di un insieme di rivelatori asservito ad un sistema di <math>M </math> fasci preformati; ciascun fascio è collegato al proprio rivelatore e all'uscita dell'integratore di ciascuno degli <math>M </math> canali viene eseguita la compensazione a <math>0.1 \ N </math> e l'apposizione della soglia <math> T </math> al livello di <math> T = 0.27.</math>

A questo punto è interessante valutare la probabilità di falso allarme <math>Pfa </math> su tutti gli <math>M </math> canali pensando, ragionevolmente, che la compensazione entro <math>0.1 \ N</math> possa essere sensibilmente diversa da canale a canale.

Ripetendo la procedura di calcolo, per soglia considerata fissa a <math> T = 0.27,</math> ipotizzando una variazione, se pur minima, degli <math>M</math>delle tolleranze i di compensazione vediamo come la <math> Pfa </math> varia:

* per tolleranza di compensazione del <math>10 \% \ la \ Pfa = 10 \% </math>
* per tolleranza di compensazione del <math>9.5 \% \ la \ Pfa = 0.1 \% </math>

Si conclude quindi che non tutti i fasci preformati, ferma la soglia a <math>T =  0.27 </math>, possono avere lo stesso valore di <math>Pfa</math>; dal più elevato  <math>10 \% </math> fissato come impostazione dello studio, a valori nettamente meno penalizzanti per la scoperta sonar.

==note==
<references/>

==Bibliografia==
* 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.

* R.  J. Urick,  ''Principles of underwater sound'', Mc Graw – hill,  3^ ed. 1968
 
* Milton Abramowitz ..''Handbook of mathematical functios'', USA 1970