Geologia e Sismologia - Studio di Sismografi e Precursori
come eliminare il 50Hz

VISUALIZZA L'ALBUM

Pag: (3)   [1]   2   3    (Ultimo)


FinePagina

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 20/06/2020 17:25:23

per iniziare un segnale che si suppone della stessa ampiezza del segnale a 50Hz

senza mascheratura e con mascheratura dati

poi un segnale con ampiezza 1/1000 rispetto a quello di disturbo a 50 Hz
senza mascheratura e con mascheratura dati

poi un segnale 1/100.000 rispetto a quello di disturbo a 50 Hz
senza mascheratura e con mascheratura dati



Modificato da qqcreafis - 20/06/2020, 23:23:25


---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 20/06/2020 21:55:04

il segnale e il 50 hz hanno la stessa ampiezza, il filtro passa basso se la cava alla grande con la mascheratura



Immagine Allegata: unomascherato.png
 
ForumEA/U/unomascherato.png



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 20/06/2020 22:09:42

se si elimina la maschera si nota solo un piccolo peggioramento (N.B sono tagliati i primi e gli ultimi dati )



Immagine Allegata: unononmascherato.png
 
ForumEA/U/unononmascherato.png



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 20/06/2020 22:14:35

paasiamo ora ad un segnale di 1/1000 dell'ampiezza dei 50 Hz



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 20/06/2020 22:33:08

con mascheratura l'uscita del filtro è ancora molto buona

N.B. l'ampiezza del segnale è 1/1000 del 50Hz



Immagine Allegata: unosumillemascherato.png
 
ForumEA/U/unosumillemascherato.png



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 20/06/2020 22:50:35

se non si maschera ci sono casini



Immagine Allegata: unosumillenonmascherato.png
 
ForumEA/U/unosumillenonmascherato.png



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 20/06/2020 23:21:15

passiamo ora ad un 1/100.000 con maschera il risultato è ancora buono



Immagine Allegata: unosucentomilamascherato.png
 
ForumEA/U/unosucentomilamascherato.png



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

Bolle

Avatar
GigaWatt


AMMINISTRATORE
Msg:9362

Stato:



Inviato il: 21/06/2020 12:57:17

Prof, sono curioso di capire dove approdi.

Cmq con il SW dell'EArthQQuake già facciamo tutto questo...questo è un pezzo estratto dal simulatore che ci fa vedere il dettaglio di cosa si combina con il segnale...la parte centrale (colore nero racchiuso tra le linee verticali viola e blu) è quello che si prende in considerazione...processamento dopo processamento.

http://www.energialternativa.info/public/newforum/ForumEA/U/RicostruzioneSegnaleCentrato.png



E' lo spike che ci crea problemi...cioè quando arriva un picco non periodico...quello ci fa male!



---------------
Un risultato se non è ripetibile non esiste (by qqcreafis).




 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 22/06/2020 21:00:34


CITAZIONE

Prof, sono curioso di capire dove approdi.
mi sto riferendo a segnali veramente annegati in profondità nel 50Hz

comunque spiego in dettaglio il metodo
1) si maschera(un set di 1024 (o quello che si crede)dati viene moltiplicato per una funzione che è zero ai bordi ed uno al centro , ho usato ad es sen(pi.greco()/1023*(i-1))

2) si trasforma
3) si f iltra (si annullano i canali al di sopra dello spettro del segnale)
4) si antitrasforma
5) si smaschera (dividendo per il valore della maschera, ovviamente tagliando i bordi)

e poi ovviamente si tiene solo la parte centrale e si riaggiorna il set di dati ecc ecc



Modificato da qqcreafis - 22/06/2020, 21:25:36


---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 22/06/2020 23:04:09

E' lo spike che ci crea problemi...cioè quando arriva un picco non periodico...quello ci fa male![/quote]

quanto è largo lo spike?



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

Bolle

Avatar
GigaWatt


AMMINISTRATORE
Msg:9362

Stato:



Inviato il: 23/06/2020 08:10:32

Ce ne è di tutti i tipi...tendenzialmente l'ampiezza varia molto e la durata è legata all'ampiezza...questo è un esempio:



In questo esempio si vede la presenza della 50hz che viene spazzata via...si vede anche il segnale ricostruito (di colore nero) a sinistra della finestrella centrale...il tratto nero della parte centrale è il segnale che viene ricostruito ad ogni battuta e si aggiunge al precedente...per questo motivo è anche chiaro sul motivo del ritardo di 5 secondi che abbiamo (1024*10ms/2).



---------------
Un risultato se non è ripetibile non esiste (by qqcreafis).




 

Bolle

Avatar
GigaWatt


AMMINISTRATORE
Msg:9362

Stato:



Inviato il: 18/04/2021 18:08:14

Prof,con quasi un anno di ritardo ....ho implementato la tua maschera per vedere cosa accade sui dati reali.

Allora...se la facciamo girare da sola senza i vari algoritmi e filtri di stabilizzazione...succede il solito casino e non si nota nulla...mentre se si applica insieme a tutto l'ambaradam (filtri ed algoritmi esotici) si nota un leggero miglioramento...ma considerando che si pensava di essere arrivati al top...


-TOP CODE-

-BOTTOM CODE-





http://www.energialternativa.info/public/newforum/ForumEA/U/MaskNoMask1.png



Vediamo un dettaglio zummolato:

http://www.energialternativa.info/public/newforum/ForumEA/U/MaskNoMask2.png




Poca roba? Sempre meglio che niente!

Ma cosa accade se arriva uno spike? Bè ...qua la ricostruzione migliora nettamente...anche se...non è proprio il top ma cmq c'è un bel salto in avanti!!

http://www.energialternativa.info/public/newforum/ForumEA/U/SegnaleMaskNoMask.png



questo con i dati in presa diretta del geofono e il pennino ...

http://www.energialternativa.info/public/newforum/ForumEA/U/GeofonoMaskNoMask.png




Il segnale sopra è lo stesso di questo video...cosa dire? La maschera diminuisce sensibilmente gli artefatti...c'è un bel miglioramento...grande prof!!!Faccine/Smileys_0006.gif



Modificato da Bolle - 18/04/2021, 18:55:49


---------------
Un risultato se non è ripetibile non esiste (by qqcreafis).




 

Bolle

Avatar
GigaWatt


AMMINISTRATORE
Msg:9362

Stato:



Inviato il: 25/11/2021 15:33:29

Facciamo un'altro passetto avanti ... mentre si provava il sismografo di NE mi sono detto...ma se invece della maschera del sin(x) utilizzassimo una maschera sin(x)^y???
Tatààààààà:

http://www.energialternativa.info/public/newforum/ForumEA/U/Sin_X_AllaY.jpg



Le rappresentazioni sono sin^1 (blu) sin^2 sin^3 sin^4 sin^10(blu) e sin^22(arancione) .


Bè...con sin(x)^5...otteniamo un eccellente risultato!

Forse andrebbe la pena provare anche qualche gaussiana...ba bo bu!



Modificato da Bolle - 30/11/2021, 13:08:09


---------------
Un risultato se non è ripetibile non esiste (by qqcreafis).




 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 27/11/2021 21:24:30


CITAZIONE

Forse andrebbe la pena provare anche qualche gaussiana...ba bo bu!


giusto!



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 

qqcreafis

Avatar
PetaWatt


Utente
Msg:10028

Stato:



Inviato il: 27/11/2021 23:04:56

P.T.

incredibile si riesce ad eliminare il 50 Hz anche con un rapporto
ampiezza segnale /ampiezza 50Hz =10^-6

(1 volt di 50 Hz e 1microvolt di segnale)

il risultato è quasi perfetto a con mask=3



Immagine Allegata: PT.png
 
ForumEA/U/PT.png



---------------
ODE AD UNO STUDENTE MERITEVOLE (Sermone)
Allora Xyz è stato bravo! noi dobbiamo riconoscere ...

 
 InizioPagina
 

Pag: (3)   [1]   2   3    (Ultimo)

Versione Completa!

Home page     TOP100-SOLAR    Home page forum