AJO_Ajustelineal_Ejemplo1_triangulo_sinpesos

110 days ago by sevillad1819

L = vector([99,53,49]).column() # observaciones C = vector([0,-200,0]).column() # constantes geom├ętricas L = L+C A = matrix([[1,0],[-1,-1],[0,1]]) P = identity_matrix(3) # esto es sin pesos #P = diagonal_matrix([10,1,1]) #P = diagonal_matrix([10,10,10]) N = transpose(A)*P*A T = transpose(A)*P*L X = N^(-1)*T X.n() 
       
[98.6666666666667]
[48.6666666666667]
[98.6666666666667]
[48.6666666666667]
LmasV = A*X-C LmasV.n() 
       
[98.6666666666667]
[52.6666666666667]
[48.6666666666667]
[98.6666666666667]
[52.6666666666667]
[48.6666666666667]
A*X-L 
       
[-1/3]
[-1/3]
[-1/3]
[-1/3]
[-1/3]
[-1/3]