Tengo este sistema de matrices de ecuaciones lineales que es inconsistente y no tiene ninguna solución, pero estoy tratando de resolverlo utilizando una octava en forma de matriz.
A = [1, -1, 2; 1, 0, 1; 2, -3, 5; 3, 2, -1] B = [4; 6; 4; 1]
Cuando lo procese así AB
me da esto
ans = [- 5.3182; 14.2273; 11.5000]
pero este ejemplo es bien conocido por ser inconsistente por qué me da respuesta
Respuestas
1 para la respuesta № 1Interesante pregunta. Sin duda, como usted sabe, en primer lugar está tratando de resolver 4 ecuaciones en 3 variables independientes. Esto solo es posible cuando la cuarta ecuación se puede deducir de la primera 3. En segundo lugar, la parte cuadrada del sistema: las primeras 3 ecuaciones (filas) están degeneradas, con cero determinante: pruebe
det(A(1:3,:))
Entonces, ¿qué está haciendo el solucionador numérico en este caso? Estoy bastante seguro de que la respuesta está en el manual:
4. Si la matriz no es cuadrada, o cualquiera de los solucionadores anteriores marca una matriz singular o casi singular, encuentre una solución de mínimos cuadrados utilizando la función LAPACK xGELSD.