Geoestadistica trabajo

...

-

CASOS PARA LAS FUNCIONES AUXULIARES

EN UNA DIMENSION

CASO 1

[pic 17][pic 18]

Donde: [pic 19]

Donde siguiendo evaluando para cada h tenemos:[pic 20][pic 21]

[pic 22]

[pic 23]

[pic 24]

[pic 25]

[pic 26]

Entonces obtenemos el valor la variabilidad para cada hk.[pic 27]

Finalmente sumamos y obtenemos el promedio de todas las variabilidades encontradas:

[pic 28]

X(L) = (A,AB)= 1/L ()[pic 29][pic 30]

CASO 2

[pic 31][pic 32][pic 33]

Donde: [pic 34]

Analizando igual que en el caso anterior obtenemos:

[pic 35]

Equivale a :

[pic 36]

X(L) = (AB,AB)= 1/L ()[pic 37][pic 38]

EN DOS DIMESIONES

CASO A= H(L,l)=(B,ABCD)[pic 40][pic 39]

CASO B H(L,l)=(AC,BD)[pic 41]

CASO C H(L,l)=(AC,ABCD)[pic 42]

CASO D H(L,l)=(AC,AB)[pic 43]

CASO E H(L,l)=(ABCD,ABCD)[pic 44]

-

PROGRAMACION PARA EL CALCULO DE VARIOGRAMA PROMEDIO Y CURVAS ISOVALORICAS

Codificación en Prog. Matlab R2015a:

CASOS UNA DIMENSION

CASO 1

[pic 45]

[pic 46]

[pic 47]

CAS0 2[pic 48]

[pic 49]

Caso A

clc;

clear;

disp('PROGRAMA PARA HACER UN VARIOGRAMA PROMEDIO')

disp('TIPO A')

disp('para una malla de 5x5 datos')

disp('para un modelo esférico')

a=[0:5];

b=[0:5];

c=input('ingresar la constante C = ');

alcance=input('ingresar el alcance a = ');

%PROCEDIMIENTO

a1=length(a);

b1=length(b);

for a=1:a1

for b=1:b1

rx=zeros(a,b);

ry=zeros(a,b);

r=zeros(a,b);

for i=1:a

for j=1:b

ry(i,j)=ry(i,j)+j;

rx(i,j)=rx(i,j)+i;

end

end

h=ones(a,b);

for p=1:a

for q=1:b

h(p,q)=h(p,q)*(sqrt(((a-rx(p,q))^2)+(ry(p,q)-1)^2));

end

end

v=zeros(a,b);

for l=1:a

for k=1:b

if h(l,k)==0

v(l,k)=(c*((((3/2)*(h(l,k)/alcance))-((1/2)*((h(l,k)/alcance)^3)))));

end

if h(l,k)

v(l,k)=(c*((((3/2)*(h(l,k)/alcance))-((1/2)*((h(l,k)/alcance)^3)))));

else

v(l,k)=c;

end

end

end

%FUNCION VARIOGRAMA PROMEDIO GENERAL

VPG(a,b)=((sum(sum(v)))/(2*a*b));

end

end

%CURVAS ISOVALORICAS

figure(1)

[x,y]=meshgrid(0:5);

z=VPG;

[c1,h1]=contour(x,y,z,20), axis square;

clabel(c1,h1);

title('CURVAS ISOVALORICAS PARA LA FUNCION TIPO A');

grid on;

shg;

Caso B

clc;

clear;

disp('PROGRAMA PARA HACER UN VARIOGRAMA PROMEDIO');

disp('TIPO B');

disp('para una malla de 5x5 datos');

disp('para un modelo esferico')

a=[0:5];

b=[0:5];

y=[0:5];

c=input('ingresar