Articulo de referencia

Núcleo (álgebra lineal)

Un ejemplo de núcleo: el operador lineal. L : ( incógnita , y ) ⟶ ( incógnita , incógnita ) {\displaystyle L:(x,y)\longrightarrow (x,x)} transforma todos los puntos en el ( incó...

Un ejemplo de núcleo: el operador lineal.L:(incógnita,y)(incógnita,incógnita){\displaystyle L:(x,y)\longrightarrow (x,x)}transforma todos los puntos en el(incógnita=0,y){\displaystyle (x=0,y)}línea hasta el punto cero(0,0){\displaystyle (0,0)}, por lo tanto, forman el núcleo del operador lineal

En matemáticas , el núcleo de una aplicación lineal , también conocido como espacio nulo , es la parte del dominio que se mapea al vector cero del codominio ; el núcleo es siempre un subespacio lineal del dominio. [ 1 ] Es decir, dada una aplicación lineal L : VW entre dos espacios vectoriales V y W , el núcleo de L es el espacio vectorial de todos los elementos v de V tales que L ( v ) = 0 , donde 0 denota el vector cero en W , [ 2 ] o de forma más simbólica:  ker(L)={vVL(v)=0}=L1(0).{\displaystyle \ker(L)=\left\{\mathbf {v} \in V\mid L(\mathbf {v} )=\mathbf {0} \right\}=L^{-1}(\mathbf {0} ).}

Propiedades

Núcleo e imagen de una aplicación lineal L de V a W

El núcleo de L es un subespacio lineal del dominio V. [ 3 ] [ 2 ]

En el mapa linealL:VW,{\displaystyle L:V\to W,}dos elementos de V tienen la misma imagen en W si y solo si su diferencia se encuentra en el núcleo de L , es decir,

L(v1)=L(v2)L(v1v2)=0.{\displaystyle L\left(\mathbf {v} _{1}\right)=L\left(\mathbf {v} _{2}\right)\quad \iff \quad L\left(\mathbf {v} _{1}-\mathbf {v} _{2}\right)=\mathbf {0} .}

De esto se deduce, por el primer teorema de isomorfismo , que la imagen de L es isomorfa al cociente de V por el núcleo: soy(L)V/ker(L).{\displaystyle \operatorname {im} (L)\cong V/\ker(L).}En el caso en que V sea de dimensión finita , esto implica el teorema de rango-nulidad : oscuro(kerL)+oscuro(soyL)=oscuro(V).{\displaystyle \dim(\ker L)+\dim(\operatorname {im} L)=\dim(V).} donde el términoEl rango se refiere a la dimensión de la imagen deL,oscuro(soyL),{\displaystyle \dim(\operatorname {estoy} L),}mientrasLa nulidad se refiere a la dimensión del núcleo deL,oscuro(kerL).{\displaystyle \dim(\ker L).}[ 4 ] Es decir, Rango(L)=oscuro(soyL) y Nulidad(L)=oscuro(kerL),{\displaystyle \operatorname {Rango} (L)=\dim(\operatorname {im} L)\qquad {\text{ y }}\qquad \operatorname {Nulidad} (L)=\dim(\ker L),} de modo que el teorema de rango-nulidad pueda reformularse como Rango(L)+Nulidad(L)=oscuro(dominioL).{\displaystyle \operatorname {Rango} (L)+\operatorname {Nulidad} (L)=\dim \left(\operatorname {dominio} L\right).}

Cuando V es un espacio con producto interno , el cocienteV/ker(L){\displaystyle V/\ker(L)}puede identificarse con el complemento ortogonal en V deker(L){\displaystyle \ker(L)}. Esta es la generalización a operadores lineales del espacio fila , o coimagen , de una matriz.

Generalización a módulos

La noción de núcleo también tiene sentido para los homomorfismos de módulos , que son generalizaciones de espacios vectoriales donde los escalares son elementos de un anillo , en lugar de un cuerpo . El dominio de la aplicación es un módulo, y el núcleo constituye un submódulo . En este caso, los conceptos de rango y nulidad no son necesariamente aplicables.

En el análisis funcional

Si V y W son espacios vectoriales topológicos tales que W es de dimensión finita, entonces un operador lineal L : VW es continuo si y solo si el núcleo de L es un subespacio cerrado de V.

Representación como multiplicación de matrices

Consideremos un mapa lineal representado como una matriz A de m × n con coeficientes en un campo K (típicamenteR{\displaystyle \mathbb {R} }odo{\displaystyle \mathbb {C} }), es decir, operando sobre vectores columna x con n componentes sobre K. El núcleo de esta aplicación lineal es el conjunto de soluciones de la ecuación A x = 0 , donde 0 se entiende como el vector cero . La dimensión del núcleo de A se llama nulidad de A. En notación de constructor de conjuntos , norte(A)=Nulo(A)=ker(A)={incógnitaKnorteAincógnita=0}.{\displaystyle \operatorname {N} (A)=\operatorname {Null} (A)=\operatorname {ker} (A)=\left\{\mathbf {x} \in K^{n}\mid A\mathbf {x} =\mathbf {0} \right\}.} La ecuación matricial es equivalente a un sistema homogéneo de ecuaciones lineales : Aincógnita=0a11incógnita1+a12incógnita2++a1norteincógnitanorte=0a21incógnita1+a22incógnita2++a2norteincógnitanorte=0 ametro1incógnita1+ametro2incógnita2++ametronorteincógnitanorte=0.{\displaystyle A\mathbf {x} =\mathbf {0} \;\;\Leftrightarrow \;\;{\begin{alignedat}{7}a_{11}x_{1}&&\;+\;&&a_{12}x_{2}&&\;+\;\cdots \;+\;&&a_{1n}x_{n}&&\;=\;&&&0\\a_{21}x_{1}&&\;+\;&&a_{22}x_{2}&&\;+\;\cdots \;+\;&&a_{2n}x_{n}&&\;=\;&&&0\\&&&&&&&&&&\vdots \ \;&&&\\a_{m1}x_{1}&&\;+\;&&a_{m2}x_{2}&&\;+\;\cdots \;+\;&&a_{mn}x_{n}&&\;=\;&&&0{\text{.}}\\\end{alignedat}}} Por lo tanto, el núcleo de A es el mismo que el conjunto solución de las ecuaciones homogéneas anteriores.

Propiedades del subespacio

El núcleo de una matriz m × n A sobre un cuerpo K es un subespacio lineal de K n . Es decir, el núcleo de A , el conjunto Null( A ) , tiene las siguientes tres propiedades:

  1. Null( A ) siempre contiene el vector cero , ya que A 0 = 0 .
  2. Si x ∈ Null( A ) e y ∈ Null( A ) , entonces x + y ∈ Null( A ) . Esto se deduce de la propiedad distributiva de la multiplicación de matrices sobre la suma.
  3. If x ∈ Null(A) and c is a scalarcK, then cx ∈ Null(A), since A(cx) = c(Ax) = c0 = 0.

The row space of a matrix

The product Ax can be written in terms of the dot product of vectors as follows: Ax=[a1xa2xamx].{\displaystyle A\mathbf {x} ={\begin{bmatrix}\mathbf {a} _{1}\cdot \mathbf {x} \\\mathbf {a} _{2}\cdot \mathbf {x} \\\vdots \\\mathbf {a} _{m}\cdot \mathbf {x} \end{bmatrix}}.}

Here, a1, ... , am denote the rows of the matrix A. It follows that x is in the kernel of A, if and only if x is orthogonal (or perpendicular) to each of the row vectors of A (since orthogonality is defined as having a dot product of 0).

The row space, or coimage, of a matrix A is the span of the row vectors of A. By the above reasoning, the kernel of A is the orthogonal complement to the row space. That is, a vector x lies in the kernel of A, if and only if it is perpendicular to every vector in the row space of A.

The dimension of the row space of A is called the rank of A, and the dimension of the kernel of A is called the nullity of A. These quantities are related by the rank–nullity theorem[4]rank(A)+nullity(A)=n.{\displaystyle \operatorname {rank} (A)+\operatorname {nulity} (A)=n.}

Left null space

The left null space, or cokernel, of a matrix A consists of all column vectors x such that xTA = 0T, where T denotes the transpose of a matrix. The left null space of A is the same as the kernel of AT. The left null space of A is the orthogonal complement to the column space of A, and is dual to the cokernel of the associated linear transformation. The kernel, the row space, the column space, and the left null space of A are the four fundamental subspaces associated with the matrix A.

Nonhomogeneous systems of linear equations

The kernel also plays a role in the solution to a nonhomogeneous system of linear equations: Ax=bora11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2 am1x1+am2x2++amnxn=bm{\displaystyle A\mathbf {x} =\mathbf {b} \quad {\text{or}}\quad {\begin{alignedat}{7}a_{11}x_{1}&&\;+\;&&a_{12}x_{2}&&\;+\;\cdots \;+\;&&a_{1n}x_{n}&&\;=\;&&&b_{1}\\a_{21}x_{1}&&\;+\;&&a_{22}x_{2}&&\;+\;\cdots \;+\;&&a_{2n}x_{n}&&\;=\;&&&b_{2}\\&&&&&&&&&&\vdots \ \;&&&\\a_{m1}x_{1}&&\;+\;&&a_{m2}x_{2}&&\;+\;\cdots \;+\;&&a_{mn}x_{n}&&\;=\;&&&b_{m}\\\end{alignedat}}} If u and v are two possible solutions to the above equation, then A(uv)=AuAv=bb=0{\displaystyle A(\mathbf {u} -\mathbf {v} )=A\mathbf {u} -A\mathbf {v} =\mathbf {b} -\mathbf {b} =\mathbf {0} } Thus, the difference of any two solutions to the equation Ax = b lies in the kernel of A.

It follows that any solution to the equation Ax = b can be expressed as the sum of a fixed solution v and an arbitrary element of the kernel. That is, the solution set to the equation Ax = b is {v+xAv=bxNull(A)},{\displaystyle \left\{\mathbf {v} +\mathbf {x} \mid A\mathbf {v} =\mathbf {b} \land \mathbf {x} \in \operatorname {Null} (A)\right\},} Geometrically, this says that the solution set to Ax = b is the translation of the kernel of A by the vector v. See also Fredholm alternative and flat (geometry).

Illustration

The following is a simple illustration of the computation of the kernel of a matrix (see § Computation by Gaussian elimination, below for methods better suited to more complex calculations). The illustration also touches on the row space and its relation to the kernel.

Consider the matrix A=[235423].{\displaystyle A={\begin{bmatrix}2&3&5\\-4&2&3\end{bmatrix}}.} The kernel of this matrix consists of all vectors (x, y, z) ∈ R3 for which [235423][xyz]=[00],{\displaystyle {\begin{bmatrix}2&3&5\\-4&2&3\end{bmatrix}}{\begin{bmatrix}x\\y\\z\end{bmatrix}}={\begin{bmatrix}0\\0\end{bmatrix}},} which can be expressed as a homogeneous system of linear equations involving x, y, and z: 2x+3y+5z=0,4x+2y+3z=0.{\displaystyle {\begin{aligned}2x+3y+5z&=0,\\-4x+2y+3z&=0.\end{aligned}}}

The same linear equations can also be written in matrix form as: [23504230].{\displaystyle \left[{\begin{array}{ccc|c}2&3&5&0\\-4&2&3&0\end{array}}\right].}

Through Gauss–Jordan elimination, the matrix can be reduced to: [101/1600113/80].{\displaystyle \left[{\begin{array}{ccc|c}1&0&1/16&0\\0&1&13/8&0\end{array}}\right].}

Rewriting the matrix in equation form yields: x=116zy=138z.{\displaystyle {\begin{aligned}x&=-{\frac {1}{16}}z\\y&=-{\frac {13}{8}}z.\end{aligned}}}

The elements of the kernel can be further expressed in parametric vector form, as follows: [xyz]=c[1/1613/81](where cR){\displaystyle {\begin{bmatrix}x\\y\\z\end{bmatrix}}=c{\begin{bmatrix}-1/16\\-13/8\\1\end{bmatrix}}\quad ({\text{where }}c\in \mathbb {R} )}

Since c is a free variable ranging over all real numbers, this can be expressed equally well as: [xyz]=c[12616].{\displaystyle {\begin{bmatrix}x\\y\\z\end{bmatrix}}=c{\begin{bmatrix}-1\\-26\\16\end{bmatrix}}.} The kernel of A is precisely the solution set to these equations (in this case, a line through the origin in R3). Here, the vector (−1,−26,16)T constitutes a basis of the kernel of A. The nullity of A is therefore 1, as it is spanned by a single vector.

The following dot products are zero: [235][12616]=0and[423][12616]=0,{\displaystyle {\begin{bmatrix}2&3&5\end{bmatrix}}{\begin{bmatrix}-1\\-26\\16\end{bmatrix}}=0\quad \mathrm {and} \quad {\begin{bmatrix}-4&2&3\end{bmatrix}}{\begin{bmatrix}-1\\-26\\16\end{bmatrix}}=0,} which illustrates that vectors in the kernel of A are orthogonal to each of the row vectors of A.

These two (linearly independent) row vectors span the row space of A—a plane orthogonal to the vector (−1,−26,16)T.

With the rank 2 of A, the nullity 1 of A, and the dimension 3 of A, we have an illustration of the rank-nullity theorem.

Examples

  • If L: RmRn, then the kernel of L is the solution set to a homogeneous system of linear equations. As in the above illustration, if L is the operator: L(x1,x2,x3)=(2x1+3x2+5x3,4x1+2x2+3x3){\displaystyle L(x_{1},x_{2},x_{3})=(2x_{1}+3x_{2}+5x_{3},\;-4x_{1}+2x_{2}+3x_{3})} then the kernel of L is the set of solutions to the equations 2x1+3x2+5x3=04x1+2x2+3x3=0{\displaystyle {\begin{alignedat}{7}2x_{1}&\;+\;&3x_{2}&\;+\;&5x_{3}&\;=\;&0\\-4x_{1}&\;+\;&2x_{2}&\;+\;&3x_{3}&\;=\;&0\end{alignedat}}}
  • Let C[0,1] denote the vector space of all continuous real-valued functions on the interval [0,1], and define L: C[0,1] → R by the rule L(f)=f(0.3).{\displaystyle L(f)=f(0.3).} Then the kernel of L consists of all functions fC[0,1] for which f(0.3) = 0.
  • Let C(R) be the vector space of all infinitely differentiable functions RR, and let D: C(R) → C(R) be the differentiation operator: D(f)=dfdx.{\displaystyle D(f)={\frac {df}{dx}}.} Then the kernel of D consists of all functions in C(R) whose derivatives are zero, i.e. the set of all constant functions.
  • Let R be the direct product of infinitely many copies of R, and let s: RR be the shift operators(x1,x2,x3,x4,)=(x2,x3,x4,).{\displaystyle s(x_{1},x_{2},x_{3},x_{4},\ldots )=(x_{2},x_{3},x_{4},\ldots ).} Then the kernel of s is the one-dimensional subspace consisting of all vectors (x1, 0, 0, 0, ...).
  • If V is an inner product space and W is a subspace, the kernel of the orthogonal projectionVW is the orthogonal complement to W in V.

Computation by Gaussian elimination

La base del núcleo de una matriz se puede calcular mediante la eliminación gaussiana .

Para ello, dada una matriz A de m × n , primero construimos la matriz aumentada por filas.[AI],{\displaystyle {\begin{bmatrix}A\\\hline I\end{bmatrix}},}donde I es la matriz identidad n × n .

Calculando su forma escalonada de columna mediante eliminación gaussiana (o cualquier otro método adecuado), obtenemos una matriz[Bdo].{\displaystyle {\begin{bmatrix}B\\\hline C\end{bmatrix}}.}Una base del núcleo de A consiste en las columnas no nulas de C tales que la columna correspondiente de B es una columna nula .

De hecho, el cálculo puede detenerse tan pronto como la matriz superior esté en forma escalonada por columnas: el resto del cálculo consiste en cambiar la base del espacio vectorial generado por las columnas cuya parte superior es cero.

Por ejemplo, supongamos que A=[103028015014000179000000].{\displaystyle A={\begin{bmatrix}1&0&-3&0&2&-8\\0&1&5&0&-1&4\\0&0&0&1&7&-9\\0&0&0&0&0&0\end{bmatrix}}.} Entonces [AI]=[103028015014000179000000100000010000001000000100000010000001].{\displaystyle {\begin{bmatrix}A\\\hline I\end{bmatrix}}={\begin{bmatrix}1&0&-3&0&2&-8\\0&1&5&0&-1&4\\0&0&0&1&7&-9\\0&0&0&0&0&0\\\hline 1&0&0&0&0&0\\0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&1&0&0\\0&0&0&0&1&0\\0&0&0&0&0&1\end{bmatrix}}.}

Al colocar la parte superior en forma escalonada de columnas mediante operaciones de columna en toda la matriz se obtiene [Bdo]=[100000010000001000000000100328010514000100001079000010000001].{\displaystyle {\begin{bmatrix}B\\\hline C\end{bmatrix}}={\begin{bmatrix}1&0&0&0&0&0\\0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&0&0&0\\\hline 1&0&0&3&-2&8\\0&1&0&-5&1&-4\\0&0&0&1&0&0\\0&0&1&0&-7&9\\0&0&0&0&1&0\\0&0&0&0&0&1\end{bmatrix}}.}

Las últimas tres columnas de B son columnas cero. Por lo tanto, los tres últimos vectores de C , [351000],[210710],[840901]{\displaystyle \left[\!\!{\begin{array}{r}3\\-5\\1\\0\\0\\0\end{array}}\right],\;\left[\!\!{\begin{array}{r}-2\\1\\0\\-7\\1\\0\end{array}}\right],\;\left[\!\!{\begin{array}{r}8\\-4\\0\\9\\0\\1\end{array}}\right]} son una base del núcleo de A.

Prueba de que el método calcula el núcleo: Dado que las operaciones de columna corresponden a la postmultiplicación por matrices invertibles, el hecho de que[AI]{\displaystyle {\begin{bmatrix}A\\\hline I\end{bmatrix}}}se reduce a[Bdo]{\displaystyle {\begin{bmatrix}B\\\hline C\end{bmatrix}}}significa que existe una matriz invertiblePAG{\displaystyle P}de tal manera que[AI]PAG=[Bdo],{\displaystyle {\begin{bmatrix}A\\\hline I\end{bmatrix}}P={\begin{bmatrix}B\\\hline C\end{bmatrix}},}conB{\displaystyle B}en forma escalonada de columna. Por lo tantoAPAG=B{\displaystyle AP=B},IPAG=do{\displaystyle IP=C}, yAdo=B{\displaystyle AC=B}. Un vector columnav{\displaystyle \mathbf {v} }pertenece al núcleo deA{\displaystyle A}(eso esAv=0{\displaystyle A\mathbf {v} =\mathbf {0} }) si y solo siBw=0,{\displaystyle B\mathbf {w} =\mathbf {0} ,}dóndew=PAG1v=do1v{\displaystyle \mathbf {w} =P^{-1}\mathbf {v} =C^{-1}\mathbf {v} }. ComoB{\displaystyle B}está en forma escalonada de columna,Bw=0{\displaystyle B\mathbf {w} =\mathbf {0} }, si y solo si las entradas no nulas dew{\displaystyle \mathbf {w} }corresponden a las columnas cero deB{\displaystyle B}. Multiplicando pordo{\displaystyle C}, se puede deducir que este es el caso si y solo siv=dow{\displaystyle \mathbf {v} =C\mathbf {w} }es una combinación lineal de las columnas correspondientes dedo{\displaystyle C}.

Cálculo numérico

El problema de calcular el núcleo en un ordenador depende de la naturaleza de los coeficientes.

Coeficientes exactos

Si los coeficientes de la matriz son números exactos, la forma escalonada por columnas de la matriz se puede calcular con el algoritmo de Bareiss de forma más eficiente que con la eliminación gaussiana. Resulta aún más eficiente utilizar la aritmética modular y el teorema chino del resto , lo que reduce el problema a varios problemas similares sobre campos finitos (esto evita la sobrecarga derivada de la no linealidad de la complejidad computacional de la multiplicación de enteros).

Para coeficientes en un campo finito, la eliminación gaussiana funciona bien, pero para las matrices grandes que aparecen en criptografía y en el cálculo de bases de Gröbner , se conocen mejores algoritmos, que tienen aproximadamente la misma complejidad computacional , pero son más rápidos y se comportan mejor con el hardware informático moderno .

Cálculo de punto flotante

Para matrices cuyas entradas son números de punto flotante , el problema de calcular el núcleo solo tiene sentido para matrices cuyo número de filas es igual a su rango: debido a los errores de redondeo , una matriz de punto flotante casi siempre tiene un rango completo , incluso cuando es una aproximación de una matriz de rango mucho menor. Incluso para una matriz de rango completo, es posible calcular su núcleo solo si está bien condicionada , es decir, si tiene un número de condición bajo . [ 5 ]

Incluso para una matriz de rango completo bien condicionada, la eliminación gaussiana no se comporta correctamente: introduce errores de redondeo demasiado grandes para obtener un resultado significativo. Dado que el cálculo del núcleo de una matriz es un caso particular de resolución de un sistema homogéneo de ecuaciones lineales, el núcleo puede calcularse con cualquiera de los diversos algoritmos diseñados para resolver sistemas homogéneos. Un software de vanguardia para este propósito es la biblioteca Lapack .

Véase también

Notas y referencias

  1. Weisstein, Eric W. "Núcleo" . mathworld.wolfram.com . Consultado el 9 de diciembre de 2019 .
  2. 1 2 "Núcleo (espacio nulo) | Brilliant Math & Science Wiki" . brilliant.org . Consultado el 09/12/2019 .
  3. El álgebra lineal, tal como se analiza en este artículo, es una disciplina matemática muy consolidada con numerosas fuentes. Casi todo el material de este artículo se puede encontrar en Lay (2005) , Meyer (2001) y las conferencias de Strang.
  4. 1 2 Weisstein, Eric W. "Teorema de rango-nulidad" . mathworld.wolfram.com . Consultado el 9 de diciembre de 2019 .
  5. "Copia archivada" (PDF) . Archivado del original (PDF) el 29-08-2017 . Recuperado el 14-04-2015 .{{cite web}}: CS1 mantenimiento: copia archivada como título ( enlace )

Bibliografía

  • Axler, Sheldon Jay (1997), Álgebra lineal bien hecha (2.ª  ed.), Springer-Verlag, ISBN 0-387-98259-0.
  • Lay, David C. (2005), Álgebra lineal y sus aplicaciones (3.ª  ed.), Addison Wesley, ISBN 978-0-321-28713-7.
  • Meyer, Carl D. (2001), Análisis matricial y álgebra lineal aplicada , Sociedad de Matemáticas Industriales y Aplicadas (SIAM), ISBN 978-0-89871-454-8Archivado del original el 31 de octubre de 2009.
  • Poole, David (2006), Álgebra lineal: una introducción moderna (2.ª  ed.), Brooks/Cole, ISBN 0-534-99845-3.
  • Anton, Howard (2005), Álgebra lineal elemental (versión de aplicaciones) (9.ª  ed.), Wiley International.
  • Leon, Steven J. (2006), Álgebra lineal con aplicaciones (7.ª  ed.), Pearson Prentice Hall.
  • Lang, Serge (1987). Álgebra lineal . Springer. ISBN 9780387964126.
  • Trefethen, Lloyd N.; Bau, David III (1997), Álgebra lineal numérica , SIAM, ISBN 978-0-89871-361-9.