AJO_1617_Ajustelineal_Ejemplo1_triangulo_sinpesos

hace 185 días por sevillad1617

L = vector([99,53-200,49]).column() A = matrix([[1,0],[-1,-1],[0,1]]) N = transpose(A)*A T = transpose(A)*L X = N^(-1)*T X.n() 
       
[98.6666666666667]
[48.6666666666667]
[98.6666666666667]
[48.6666666666667]
LmasV = A*X+vector([0,200,0]).column() LmasV.n() 
       
[98.6666666666667]
[52.6666666666667]
[48.6666666666667]
[98.6666666666667]
[52.6666666666667]
[48.6666666666667]