Utente:Funzioni di correlazione/Sandbox666
E' di notevole interesse l'esame del comportamento dei ricevitori in correlazione nei fasci preformati in presenza di segnale affetto da rumore; un modo realistico di tale esame è fattibile mediante simulazioni dinamiche al P.C. sviluppate per un piccolo numero di fasci generati da processi di correlazione a segnali limitati in ampiezza (correlazione per moltiplicazione di segni).
In questa lezione svilupperemo l'argomento e vedremo come sono articolate le istruzioni per la simulazione dell'intero processo con le seguenti routine:
- generazione dei rumori
- generazione dei segnali
- regolazione del rapporto
- limitazione di rumori e segnali
- correlazione tra i segnali
- formazione di alcuni fasci preformati
- presentazione video dinamica del fenomeno
Una volta appreso il meccanismo di simulazione un adatto file eseguibile consentirà di prendere confidenza con la procedura per studiarne i particolari comportamenti in fase dinamica.
Sulla formula principale
modificaNella sezione precedente abbiamo visto che l'algoritmo:
![{\displaystyle Crt=\left({\frac {2}{\pi }}\right)ArcSen\left[{\frac {1}{1+\left({\frac {Ni}{Si}}\right)^{2}}}\right]\qquad )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b593f9c037df36af073221f71381eeb3221711a8)
mostra come varia l'ampiezza delle funzione di correlazione in dipendenza del rapporto ; questa grandezza viene espressa con un numero pari al livello di una tensione continua.
Detta espressione è valida per correlazione a coincidenza di polarità.
La formula è stata verificata in via sperimentale ed ha dimostrato un'aderenza assoluta al fenomeno fisico così come mostra la figura 1:
nella figura, a tratto rosso, è tracciata la funzione teorica e, con punteggiata, l'andamento dell'ampiezza della sperimentale (simulata tramite software) ricavate in dipendenza del rapporto variabile da
La struttura dei fasci preformati
modificaCon riferimento allo schema riportato nella Lezione 11^ della materia Sistemi riceventi in correlazione possiamo immaginare un insieme di simulazione che generi segnali e rumori e costruisca un fascio preformato atto a discriminare il segnale in mezzo al rumore.
La struttura di simulazione secondo detto schema, opportunamente modificato, è riportata nella figura 2
In figura 2 si evidenziano:
in blu le componenti per la simulazione dei rumori e del segnale:
- generazione dei rumori e scorrelati tra loro
- generazione del segnale e dei segnali inquinati dal rumore
- regolazione del rapporto
in nero le altre componenti per la simulazione del correlatore che genera il fascio preformato:
- limitazione di rumori e dei segnali
- correlazione tra i segnali e integrazione
Per rendere la simulazione dei fasci preformati non eccessivamente pesante per il software ipotizziamo un numero di fasci pari a il centrale è puntato sul segnale simulato inquinato dal rumore come si vede in figura 2, gli altri otto, uguali tra loro, si suppongono colpiti dal solo rumore correlato secondo lo schema di figura 3:
In questa impostazione le strutture di simulazione ipotizzano uno scenario subacqueo come quello riportato in figura 4 dove, a causa di lobi direttivi sensibilmente spaziati tra loro, si può assumere che il fascio centrale (colore rosso) riceva segnale + disturbo mentre gli fasci laterali (colore nero) ricevano, praticamente soltanto disturbi.
I fasci sono puntati rispettivamente:
*primo a ° *secondo a ° *terzo a ° *quarto a ° *centrale a ° *primo a ° *secondo a ° *terzo a ° *quarto a °
Secondo quanto detto il livello di correlazione all'uscita del fascio centrale segue la legge indicata nella 1), mentre il livello di correlazione alle uscite degli altri fasci è nullo mancando il segnale, questi presentano soltanto la varianza dovuta al rumore che segue la legge:
La varianza , all'uscita dei correlatore, altro non è che il rumore in valore efficace espresso con un numero.
Anche il livello del fascio centrale, oltre ad avere ampiezza variabile secondo la prima espressione data, presenta una varianza dovuta al rumore e definita sempre dall'espressione di di 2).
I fasci preformati sono presentati sullo schermo del P.C. con visione tipo A e sono interpolati linearmente mediando tra le ampiezze dei fasci contigui.
Il pannello di presentazione e controllo
modificaIl pannello in oggetto, governato dal file eseguibile dinamico.exe, è mostrato in figura 5 e su di esso si vedono, oltre le tracce dei fasci preformati, i diversi comandi per la gestione del software di simulazione dinamica; li esaminiamo di seguito:
Il file per la presentazione dinamica ed i controlli è scaricabile all'indirizzo: Fasci preformati e il disturbo oppure su distfasci
Comando variazione del rapporto segnale/disturbo
modificaLa variazione di si esegue via software agendo sul blocco 4 di figura 2 (l'attenuatore dell'ampiezza del segnale).
Il comando del livello d'attenuazione comprende una variazione del rapporto che si estende da a passi di
L'impostazione del rapporto è possibile tramite il Vscroll indicato con la sigla (VSN) in figura 5, il rapporto di impostato è leggibile nell'apposito label indicato con la sigla (LSN).
Il software riceve tramite (VSN) il valore di impostato che agisce nel calcolo della funzione 1).
Comando puntamento del bersaglio
modificaIl puntamento del bersaglio, tramite indice verticale verde chiaro, si esegue mediante lo spostamento della slitta di Hscroll indicata con la sigla ( HBRQ ), il comando si estende da ° a ° a passi di °.
L'angolo impostato è leggibile nell'apposito label indicato con la sigla ( LBRQ ).
Il puntamento del bersaglio può essere fatto in automatico premendo il bottone per l'inseguimento ( BBRQ )
Variazione della costante di tempo d'integrazione
modificaLa variazione della costante di tempo si esegue via software agendo sul blocco (int) di figura 2, il software la riceve tramite il combo indicato con la sigla ( CRC ) che agisce nel calcolo della funzione 2).
Il valore di impostato è presentato sul label ( LRC ).
Variazione dell'ampiezza dei tracciato
modificaIl tracciato di figura 5 può essere amplificato o ridotto d'ampiezza a comando dei due bottoni indicati con ( BGUA ).
Regolazione livello soglia
modificaAl tracciato di figura 5 può essere sovrapposta una fascia orizzontale da utilizzare come livello di soglia (fascia e regolatore del livello sono visibili soltanto a comando)
Utilizzo del file eseguibile dinamico.exe
modificaCon il file di cui al titolo è possibile eseguire numerosi rilievi interessanti che mostrano le problematiche che si incontrano quando si opera con sistemi a fasci preformati in presenza dell'inevitabile rumore del mare.
misura con Si/Ni = + 20 dB e RC = 0.1 s <
modificaSecondo la figura 1, con , il livello della correlazione sul fascio a ° è oltre il del valore max raggiungibile e la varianza è trascurabile; questa è una tra le condizioni ottimali per la scoperta del bersaglio. come si vede in figura 6.
In figura l'indice di collimazione è posizionato sul massimo del fascio e su ( LBRQ ) si leggono ° di rilevamento.
Nella schermata dinamica vista su P.C. si può osservare come per i fasci adiacenti a ° la varianza produce sensibili ondulazioni su tutti gli fasci secondo la funzione 2), ciò non è un problema per la scoperta del bersaglio dato che questo è già rivelato dal fascio a °.
Per meglio apprezzare la varianza sugli fasci si può incrementare, momentaneamente, il guadagno con il bottone + ( BGUA ), in tal caso, restando il livello di correazione sul fascio °, potrà portare alla saturazione della presentazione con l'impossibilità di una collimazione precisa della direzione del bersaglio
misura con Si/Ni = + 12 dB e RC = 0.1 s
modificaSecondo la figura 1, con , il livello della correlazione sul fascio a ° e circa il del valore max raggiungibile e la varianza è molto piccola; questa è ancora tra le condizioni ottimali per la scoperta del bersaglio come si vede in figura 7.
In figura come esempio di mobilità dell'indice di collimazione lo si è posizionato su di un fascio a sinistra del max e su ( LBRQ ) si leggono °.
Nella schermata dinamica vista su P.C. si può osservare come per i fasci adiacenti a ° la varianza produca ondulazioni simili al caso precedente su tutti gli fasci; iniziano ad evidenziarsi leggere oscillazioni sul fascio a °.
misura con Si/Ni = + 6 dB e RC = 0.1 s
modificaSecondo la figura 1, con , il livello della correlazione sul fascio a ° e circa il del valore max raggiungibile e la varianza è piccola; questa è ancora tra le condizioni ottimali per la scoperta del bersaglio. come si vede in figura 8.
In figura l'indice di collimazione è posizionato su di un fascio a destra del max e su ( LBRQ ) si leggono °.
Nella schermata dinamica vista su P.C. si può osservare come per i fasci adiacenti a ° la varianza produca ondulazioni simili al caso precedente su tutti gli fasci, sono sensibili le ondulazioni sul fascio a ° ma non ne pregiudicano ancora il corretto rilevamento
misura con Si/Ni = 0 dB e RC = 0.1 s \ ( RC = 0.2 s.)
modificaSecondo la figura 1, con , il livello della correlazione sul fascio a ° è circa il del valore max raggiungibile e la varianza è sensibile; questa è ancora tra le condizioni buone per la scoperta del bersaglio. come si vede in figura 9.
In figura l'indice di collimazione è posizionato sul fascio a ° e su ( LBRQ )
si leggono ° di rilevamento.
misura con Si/Ni = -6 dB e RC = (da modificare osservando il video del P.C.)
modificaSecondo la figura 1, con , il livello della correlazione sul fascio a ° e circa il del valore max raggiungibile e la varianza è elevata; questa è la condizione che obbliga l'operatore a scegliere una adatta costante d'integrazione per consentire il rilevamento del bersaglio come si vede in figura 10.
In figura l'indice di collimazione è posizionato sul fascio a ° e su ( LBRQ ) si leggono ° di rilevamento.
Nella schermata dinamica vista su P.C. si può osservare come per i fasci adiacenti a ° la varianza produce ondulazioni superiori al caso precedente su tutti gli 8 fasci dato che, per veder meglio il tracciato, è stato incrementato di poco il guadagno + ( BGUA )
Per ridurre sensibilmente le ondulazioni sul fascio a °, per non pregiudicare la precisione di rilevamento, si può impostare il valore di con una visibile riduzione della varianza su tutti i fasci.
Livelli del rapporto Si/Ni inferiori a - 6 dB
modificaNegli esercizi precedenti il rapporto ha assunto valori decrescenti, da , tali che il fascio a ° è risultato sempre e nettamente superiore al disturbo degli altri fasci; ne abbiamo però constatato il peggioramento dal punto di vista delle ondulazioni e dell'ampiezza del fascio stesso.
Infatti riducendosi il rapporto l'ampiezza del fascio segue l'andamento decrescente della 1) restando invece inalterata la varianza espressa dalla 2); quest'ultima può ridursi soltanto aumentando la costante d'integrazione così come fatto in precedenza.
La prevalenza costante dell'ampiezza del fascio a ° rispetto agli altri ci consente di affermare che in ogni caso visto in precedenza la probabilità di rivelazione del fascio è stata sempre del e la probabilità di falso allarme è risultata sempre
Proseguendo con la riduzione del rapporto vedremo cambiare le accennate condizioni entrando nelle problematiche relative alla scoperta sonar a grandi distanze, problematiche già affrontate in via teorica nelle lezioni precedenti.
misura con Si/Ni = -12 dB e RC = 0.8 s.
modificaSecondo la figura 1, con , il livello della correlazione sul fascio a ° è soltanto il circa del valore max raggiungibile e la varianza è molto elevata rispetto al livello del fascio; questa condizione obbliga l'operatore a scegliere una adatta costante d'integrazione ed un opportuno valore di guadagno + ( BGUA ) per consentire il rilevamento del bersaglio come si vede in figura 11.
In figura l'indice di collimazione è posizionato sul fascio a ° e su ( LBRQ ) si leggono ° di rilevamento
Nella schermata dinamica vista su P.C. si può osservare come nei fasci adiacenti a ° la varianza produca ondulazioni a volte superiori al fascio a ° con l'insorgenza dell'incertezza della direzione di provenienza del suono: ecco quindi comparire le variabili probabilistiche accennate in precedenza , e l'introduzione del criterio della "soglia".
misura con Si/Ni = -18 dB e RC = 1.8 s.
modificaSecondo la figura 1, con , il livello della correlazione sul fascio a ° è soltanto l' circa del valore max raggiungibile e la varianza è elevatissima rispetto ad esso; questa è la condizione che obbliga l'operatore a scegliere una adatta costante d'integrazione ed un opportuno valore di guadagno + ( BGUA ) per consentire il rilevamento del bersaglio come si vede in figura 12.
In figura l'indice di collimazione è posizionato sul fascio a ° e su ( LBRQ ) si leggono ° di rilevamento.
Nella schermata dinamica vista su P.C. si può osservare come nei fasci adiacenti a ° la varianza produca ondulazioni a volte superiori al fascio a ° ( ecco la comparsa del falso allarmemcon Pfa > 0 ) e che a volte il fascio a ° è completamente assente ( ecco l'insorgere della rivelazione mancata del bersaglio con ).
Il valore di di quest'ultimo esercizio è, per esperienza diretta, il valore minore possibile, ottenibile con correlatori a coincidenza di polarità, accettando una coppia di valori probabilistici quali: e
osservazioni sulla misura con Si/Ni = -18 dB e RC = 1.8 s
modificaPer stabilire più facilmente quante volte il fascio a ° supera i livelli degli altri fasci o, per vedere se altri fasci superano il livello del primo è utile apporre una soglia grafica, variabile a comando con spostamento verticale.
Tale soglia è presentabile cliccando sul lato estremo a sinistra dello schermo, ed eliminabile cliccando sul lato estremo destro.
Una volta chiamata la soglia compare sul lato sinistro dello schermo uno scroll verticale nominato "S", contemporaneamente al centro della presentazione dei fasci compare una fascia chiara che può spostarsi in alto od in basso a comando di "S" così come si vede in figura 12-a.
Nella presentazione l'indice è sulla direzione del fascio a ° ma nell'attimo della copiatura dell'immagine il segnale era assente ( caso di )
Livelli di segnali eccedenti la soglia si hanno per fasci adiacenti non puntati sul bersaglio acustico (caso di ).
Conclusioni
modificaGli esercizi svolti indirizzano ad un impiego aperto dell'eseguibile "dinamico.exe" con il quale è possibile impostare i più diversi valori di accompagnati da opportune scelte del guadagno + ( BGUA ) per ottimizzare la migliore visione del bersaglio a °.
Nei casi in cui il valore di richieda la valutazione delle variabili probabilistiche , l'esercizio si fa più interessante ma complicato dal fatto che stabilire l'entità delle variabili citate diventa difficile; un'idea di come affrontare il problema della soglia di rivelazione è illustrato nella 17^ lezione della materia Principi, sistemi e metodologie per la localizzazione subacquea passiva .
Per chiudere questa lezione è interessante esaminare il riflesso di quanto esposto in precedenza alla luce della previsione di portata di un sonar.
Se esaminiamo la variazione d'ampiezza del fascio a ° e il comportamento degli altri fasci preformati immaginandoli facenti parte di un sonar le cui caratteristiche siano quelle elencate nella zona di sinistra di figura 13 possiamo vedere come, tali grandezze, subordinino, in dipendenza del loro valore, la previsione di portata del sonar come mostra la figura citata:
In figura la curva blu, normalizzata, indica la variazione d'ampiezza del fascio a ° in funzione della distanza , dato che la variabile , oltre che dai valori indicati a sinistra della figura stessa, dipende da .
La retta orizzontale rossa indica invece il livello di , anch'esso normalizzato, espresso in picco-picco.
Il grafico individua tre intervalli caratteristici di :
- Per il fascio ° avendo un livello decisamente superiore a quello degli altri fasci ha sempre una e una
- Per compreso tra la superiorità del fascio a ° rispetto agli altri fasci è incerta ed altrettanto incerta la definizione della coppia delle variabili probabilistiche.
- Per è l'ampiezza del livello degli altri fasci che supera l'ampiezza del fascio a ° implicacando sempre
Note
modifica
Bibliografia
modifica- R. J. Urick, Principles of underwater sound, 3ª ed., Mc Graw – Hill, 1968.
- C. Del Turco, La correlazione , Collana scientifica ed. Moderna La Spezia,1993