18.09.2013 Views

Caracterización de Arenas y Gravas con Ondas Elásticas

Caracterización de Arenas y Gravas con Ondas Elásticas

Caracterización de Arenas y Gravas con Ondas Elásticas

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Doctorado en Ingeniería - Facultad <strong>de</strong> Ingeniería - Universidad Nacional <strong>de</strong> Cuyo - Año 2007<br />

Lic. Armando Luis Imhof<br />

TECNICA DE RECONSTRUCCION ALGEBRAICA (ART)<br />

errdist<br />

i, j<br />

:=<br />

∑<br />

j<br />

S<br />

i, j<br />

( S<br />

i, j)<br />

2<br />

Range3Vec( s, e,<br />

i)<br />

:= count ← 1<br />

for i∈s, s + i..<br />

e<br />

V<br />

V ← 0.005<br />

count<br />

count ← count + 1<br />

errdis errdist T<br />

:=<br />

se forma el vector <strong>de</strong> inicio <strong>con</strong> valores unitarios (hay<br />

varias formas <strong>de</strong> hacerlo, ver SIRT). Se arranca <strong>con</strong> valores<br />

iniciales <strong>de</strong> lentitud (slowness) = 0.005 (= 200 m/s)<br />

xesti := Range3Vec( 1, N,<br />

1)<br />

primera entrada estimada <strong>de</strong> x (para 1era iteración)<br />

error entre valores reales y los calculados<br />

er := ymeas − S⋅xesti er T =<br />

1 2 3 4 5 6<br />

1 -3.282·10-3<br />

-3.318·10-3<br />

-3.425·10-3<br />

-3.597·10-3<br />

-3.823·10-3<br />

-4.096·10-3<br />

Rutina para calcular xest para iteración 's'<br />

xestfART( p)<br />

:= s ← 2<br />

i ← 1<br />

x ← xesti<br />

ejemplo <strong>de</strong> salida para<br />

iteración 150.-<br />

while<br />

x<br />

s ≤ p<br />

er ← ymeas − ( Sx ⋅ )<br />

i i i<br />

x x er errdis<br />

i i 〈〉<br />

← + ⋅<br />

s ← s + 1<br />

i ← i + 1 if i < M<br />

i ← 1 otherwise<br />

xestfART( 150)<br />

T =<br />

1 2 3 4 5 6<br />

1 3.564·10-3<br />

2.691·10-3<br />

2.835·10-3<br />

2.649·10-3<br />

2.714·10-3<br />

p := 300<br />

se fija el nro <strong>de</strong> iteraciones en 300.<br />

se transforma los resultados <strong>de</strong> la función en vector<br />

234<br />

2.627·10-3

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!