Si $\mathrm v \in \mbox{span} (\mathrm u_1, \mathrm u_2)$ , entonces hay un $\mathrm x \in \mathbb R^2$ tal que $\mathrm U \mathrm x = \mathrm v$ donde las columnas de $\mathrm U$ son $\mathrm u_1$ y $\mathrm u_2$ . Sin embargo, $\mathrm w \notin \mbox{span} (\mathrm u_1, \mathrm u_2)$ . Así, tenemos un ejemplo del famoso Mínimos cuadrados problema
$$\begin{array}{ll} \text{minimize} & \| \mathrm U \mathrm x - \mathrm w \|_2^2 \end{array}$$
cuya solución es
$$\hat{\mathrm x} := (\mathrm U^T \mathrm U)^{-1} \mathrm U^T \mathrm w$$
El vector de error es
$$\mathrm w - \mathrm U \hat{\mathrm x} = \mathrm w - \mathrm U (\mathrm U^T \mathrm U)^{-1} \mathrm U^T \mathrm w = (\mathrm I_4 - \mathrm U (\mathrm U^T \mathrm U)^{-1} \mathrm U^T) \, \mathrm w$$
donde $\mathrm P := \mathrm U (\mathrm U^T \mathrm U)^{-1} \mathrm U^T$ es la matriz de proyección que se proyecta sobre el espacio de columnas de $\mathrm U$ .
Utilizando SymPy ,
>>> U = Matrix([[1, 1], [1, -1], [1, 1], [1, -1]])
>>> U
⎡1 1 ⎤
⎢ ⎥
⎢1 -1⎥
⎢ ⎥
⎢1 1 ⎥
⎢ ⎥
⎣1 -1⎦
>>> w = Matrix([1, 2, -1, 2])
>>> w
⎡1 ⎤
⎢ ⎥
⎢2 ⎥
⎢ ⎥
⎢-1⎥
⎢ ⎥
⎣2 ⎦
>>> P = U * (U.T * U)**-1 * U.T
>>> P
⎡1/2 0 1/2 0 ⎤
⎢ ⎥
⎢ 0 1/2 0 1/2⎥
⎢ ⎥
⎢1/2 0 1/2 0 ⎥
⎢ ⎥
⎣ 0 1/2 0 1/2⎦
>>> (eye(4) - P) * w
⎡1 ⎤
⎢ ⎥
⎢0 ⎥
⎢ ⎥
⎢-1⎥
⎢ ⎥
⎣0 ⎦
>>> P * w
⎡0⎤
⎢ ⎥
⎢2⎥
⎢ ⎥
⎢0⎥
⎢ ⎥
⎣2⎦
El $2$ -La norma del vector de error es
$$\| \mathrm w - \mathrm U \hat{\mathrm x} \|_2 = \| (\mathrm I_4 - \mathrm U (\mathrm U^T \mathrm U)^{-1} \mathrm U^T) \, \mathrm w \|_2 = \| (\mathrm I_4 - \mathrm P) \, \mathrm w \|_2 = \sqrt 2$$