3.2: Espacio nulo
- Page ID
- 113108
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)El espacio nulo de una\(n\) matriz\(m\) -by-\(A\) es la colección de aquellos vectores en\(\mathbb{R}^{n}\) que se\(A\) mapea al vector cero en\(\mathbb{R}^m\). Más precisamente,
\[\mathcal{N}(A) = \{x \in \mathbb{R}^n | Ax = 0\} \nonumber\]
Ejemplo de espacio nulo
Como ejemplo, examinamos la matriz\(A\)
\[A = \begin{pmatrix} {0}&{1}&{0}&{0}\\ {-1}&{0}&{1}&{0}\\ {0}&{0}&{0}&{1} \end{pmatrix} \nonumber\]
Es bastante fácil ver que el espacio nulo de esta matriz es:
\[\mathcal{N}(A) = \{t \begin{pmatrix} {1}\\ {0}\\ {1}\\ {0} \end{pmatrix} | t \in \mathbb{R}\} \nonumber\]
Esta es una línea en\(\mathbb{R}^{4}\)
El espacio nulo responde a la pregunta de la singularidad de las soluciones a\(S \textbf{x} = \textbf{f}\). Para, si\(S \textbf{x} = \textbf{f}\) y\(S \textbf{y} = \textbf{f}\) entonces\(S(\textbf{x}-\textbf{y}) = S\textbf{x}-S\textbf{y} = \textbf{f}-\textbf{f} = 0\) y así\((\textbf{x}-\textbf{y}) \in \mathcal{N}(S)\). De ahí que una solución a\(S \textbf{x} = \textbf{f}\) sea única si, y sólo si,\(\mathcal{N} S = \{0\}\)
Método para encontrar la base
Expongamos ahora una base para el espacio nulo de una matriz arbitraria A. Observamos que resolver\(A \textbf{x} = 0\) es resolver\(A_{red} \textbf{x} = 0\). Con respecto a este último, suponemos que
\[\{c_{j} | j = \{1, \cdots, r\}\} \nonumber\]
son los índices de las columnas pivotantes y que
\[\{c_{j} | j = \{r+1, \cdots, n\}\} \nonumber\]
son los índices de las columnas no pivotantes. En consecuencia, definimos las variables de pivote r
\[\{x_{c_{j}} | j = \{1, \cdots, r\}\} \nonumber\]
y las variables\(n-r\) libres
\[\{x_{c_{j}} | j = \{r+1, \cdots, n\}\} \nonumber\]
Se resuelve\(A_{red} \textbf{x} = 0\) expresando cada una de las variables pivotantes en términos de las variables no pivotantes o libres. En el ejemplo anterior,\(x_{1}, x_{2}\), y\(x_{4}\) son pivotar mientras\(x_{3}\) está libre. Resolviendo para el pivote en términos de lo libre, encontramos\(x_{4} = 0, x_{3} = x_{1}, x_{2} = 0\), o, escrito como vector,
\[\textbf{x} = x_{3} \begin{pmatrix} {1}\\{0}\\{1}\\{0} \end{pmatrix} \nonumber\]
donde\(x_{3}\) es gratis. Como\(x_{3}\) rangos sobre todos los números reales, la x anterior traza una línea en\(\mathbb{R}_{4}\). Esta línea es precisamente el espacio nulo de\(A\). Resumiendo estos cálculos llegamos a:
Supongamos que\(A\) es m-by-n con índices de pivote\(\{c_{j} | j = \{1, \cdots, r\}\}\) e índices libres\(\{c_{j} | j = \{r+1, \cdots, n\}\}\). Una base para se\(\mathcal{N}(A)\) pueden construir\(n-r\) vectores\(\{z^{1}, z^{2}, \cdots,z^{n-r}\}\) donde\(z^{k}\), y solo\(z^k\) posee un distinto de cero en su\(c_{r+k}\) componente.
Una observación de MATLAB
Como es habitual, MATLAB tiene una manera de hacer nuestras vidas más simples. Si ha definido una matriz A y desea encontrar una base para su espacio nulo, simplemente llame a la función null (A)
. Una pequeña nota sobre esta función: si se agrega una bandera extra, 'r'
, como en null (A, 'r')
, entonces la base se muestra “racionalmente” en lugar de puramente matemáticamente. Las páginas de ayuda de MATLAB definen la diferencia entre los dos modos como el modo racional que es útil pedagógicamente y el modo matemático de más valor (¡jadeo!) matemáticamente.
Reflexiones finales sobre los espacios nulos
Hay mucho más para encontrar espacios nulos; suficiente, de hecho, para garantizar otro módulo. Un aspecto importante y el uso de los espacios nulos es su capacidad para informarnos sobre la singularidad de las soluciones. Si utilizamos el espacio de columna para determinar la existencia de una solución\(\textbf{x}\) a la ecuación\(A \textbf{x} = b\). Una vez que sabemos que existe una solución es una cuestión perfectamente razonable querer saber si esta solución es o no la única solución a este problema. La regla dura y rápida es que una solución\(\textbf{x}\) es única si y solo si el espacio nulo de\(A\) está vacío. Una manera de pensar sobre esto es considerar que si\(A \textbf{x} = 0\) no tiene una solución única entonces, por linealidad, tampoco lo hace\(A \textbf{x} = b\). Por el contrario, si\((Az = 0) \wedge (z \ne 0) \wedge (A \textbf{y} = b)\) entonces\(A(z+\textbf{y}) = b\) también.