Edita: He reescrito toda la pregunta porque no estaba clara. No estoy seguro de que esto sea mejor. Este es un problema de varios pasos utilizando MATLAB, así que he tratado de reducirlo un poco sin perder nada.
Sea $A$ sea una matriz aleatoria de 6x6, y $b$ sea una matriz aleatoria de 6x1. Dado que $A$ se generara aleatoriamente, esperaríamos que fuera no singular. El sistema $Ax = b$ debe tener una solución única.
Ahora cambiemos $A$ para que sea singular. Sustituya la tercera columna de $A$ con una combinación lineal de sus dos primeras columnas, a saber $a_3 = 4*a_1 + 3*a_2$ donde $a_1, a_2,$ y $a_3$ son la primera, segunda y tercera columnas de $A$ .
Ahora dejemos que $y$ sea una matriz aleatoria de 6x1 y que $c = Ay$ .
Calcular la forma escalonada reducida $U$ de $[A$ $c]$ .
La variable libre determinada por la forma escalonada debe ser $x_3$ . Examinando el sistema correspondiente a la matriz $U$ deberías ser capaz de determinar la solución correspondiente a $x_3 = 0$ . Sea el vector columna $w$ sea esta solución.
Ahora establece la última columna de $U$ a ceros. Así que $U$ debería corresponder ahora a la forma escalonada reducida de $(A | 0)$ . Utilice $U$ para determinar la solución del sistema homogéneo cuando la variable libre $x_3 = 1$ e introduzca su resultado como vector $z$ .
Establecer $v = w + 3 * z$ . El vector $v$ debe ser una solución al sistema $Ax = c$ . ¿Por qué? Explica. . . . ¿Cuál es el valor de la variable libre $x_3$ para esta solución?
Esta es la parte en la que estoy atascado:
¿Cómo podríamos determinar todas las posibles soluciones del sistema en términos de los vectores $w$ y $z$ ? Explícalo.