Bolle
| Inviato il: 24/07/2017 05:24:08
|
Questo è lo pseudo codice dell'algoritmo che dovrebbe far capire il metodo del calcolo dell'epicentro...con tanto di spiega nei commenti. Di fatto viene utilizzato solo il calcolo dei punti dell'intersezione di due circonferenze...ma per rendere la cosa comprensibile sono dovuto andare di grafica 
'DeltaT è il tempo di ritardo dell'onda che colpisce il secondo sismografo
Sito = "Sabina"
GeoToCartesiane 42.xyxyx, 12.xyxyx, xp, yp
Sito = "Viterbo"
GeoToCartesiane 42.xyxyx, 12.xyxyx, xp2, yp2
'Distanza tra i due sismografi
dis = Sqr((xp - xp2) * (xp - xp2) + (yp - yp2) * (yp - yp2))
V = 9 'Velocità in km/s dell'onda
TMAx = (dis / 1000) / V ' Calcolo il Tempo massimo che può impiegare l'onda a percorrere i due sismografi...dividendo per 1000 porto in km la distanza
'Il primo Raggio del cerchio è sempre la distanza che separa i due sismografi
Raggio1=dis
'Calcolo il raggio del secondo cerchio che dipende dal tempo di ritardo
Raggio2 = ((TMAx - TDeltaT) / TMAx) * dis
'I punti d'intersezione dei due cerchi sono i punti dell'epicentro
DrawCircle( xp2, yp2 , Raggio1)
DrawCircle( xp, yp , Raggio2)
CoppiaPunti=CalcolaIntersezioneCerchi( xp, yp , Raggio2,xp2, yp2 , Raggio1)
Questo è l'esempio con 1.5 Sec di ritardo, cioè Viterbo riceve l'onda dopo un secondo e mezzo:
Ora aumentando o diminuendo contemporaneamente i due raggi dello stesso passo ( che dipende dalla precisione che vogliamo avere... 1km???) otteniamo ,con l'intersezione delle due circonferenze, tutti gli ipotetici punti dell'epicentro, cioè
For i=-100 to 100 step 10
'Il primo Raggio del cerchio è sempre la distanza che sapare i due sismografi
Raggio1=dis+i
'Calcolo il raggio del secondo cerchio che dipende dal tempo di ritardo
Raggio2 = ((TMAx - TDeltaT) / TMAx) * dis +i
DrawCircle( xp2, yp2 , Raggio1)
DrawCircle( xp, yp , Raggio2)
CoppiaPunti(SismografoSabVit,i)=CalcolaIntersezioneCerchi( xp, yp , Raggio2,xp2, yp2 , Raggio1)
next i
Vediamo graficamente la cosa:
DeltaT=0 Sec
DeltaT=1 Sec
DeltaT=2 Sec
Ripetendo l'operazione per gli altri sismografi e confrontando le distanze tra tutti i punti, quella inferiore a tutte sarà il punto dell'epicentro!
Che voto ho preso?
Modificato da Bolle - 24/07/2017, 05:52:56
---------------
Un risultato se non è ripetibile non esiste (by qqcreafis).
|
