Articulo de referencia

Algoritmo de Zassenhaus

En matemáticas, el algoritmo de Zassenhaus [ 1 ] es un método para calcular una base para la intersección y la suma de dos subespacios de un espacio vectorial . Recibe su nombre...

En matemáticas, el algoritmo de Zassenhaus [ 1 ] es un método para calcular una base para la intersección y la suma de dos subespacios de un espacio vectorial . Recibe su nombre de Hans Zassenhaus , pero no se conoce ninguna publicación suya de este algoritmo. [ 2 ] Se utiliza en sistemas de álgebra computacional . [ 3 ]

Algoritmo

Aporte

Sea V un espacio vectorial y U , W dos subespacios de dimensión finita de V con los siguientes conjuntos generadores :

U=1,,norte{\displaystyle U=\langle u_ {1},\ldots,u_ {n}\rangle}

y

W=w1,,wk.{\displaystyle W=\langle w_{1},\ldots,w_{k}\rangle.}

Finalmente, dejemosB1,,Bmetro{\displaystyle B_{1},\ldots ,B_{m}}sean vectores linealmente independientes de modo quei{\displaystyle u_{i}}ywi{\displaystyle w_{i}}se puede escribir como

i=j=1metroai,jBj{\displaystyle u_{i}=\sum _{j=1}^{m}a_{i,j}B_{j}}

y

wi=j=1metrobi,jBj.{\displaystyle w_{i}=\sum _{j=1}^{m}b_{i,j}B_{j}.}

Producción

El algoritmo calcula la base de la sumaU+W{\displaystyle U+W}y una base de la intersecciónUW{\displaystyle U\cap W}.

Algoritmo

El algoritmo crea la siguiente matriz de bloques de tamaño((norte+k)×(2metro)){\displaystyle ((n+k)\times (2m))}:

(a1,1a1,2a1,metroa1,1a1,2a1,metroanorte,1anorte,2anorte,metroanorte,1anorte,2anorte,metrob1,1b1,2b1,metro000bk,1bk,2bk,metro000){\displaystyle {\begin{pmatrix}a_{1,1}&a_{1,2}&\cdots &a_{1,m}&a_{1,1}&a_{1,2}&\cdots &a_{1,m}\\\vdots &\vdots &&\vdots &\vdots &\vdots &&\vdots \\a_{n,1}&a_{n,2}&\cdots &a_{n,m}&a_{n,1}&a_{n,2}&\cdots &a_{n,m}\\b_{1,1}&b_{1,2}&\cdots &b_{1,m}&0&0&\cdots &0\\\vdots &\vdots &&\vdots &\vdots &\vdots &&\vdots \\b_{k,1}&b_{k,2}&\cdots &b_{k,m}&0&0&\cdots &0\end{pmatrix}}}

Utilizando operaciones elementales de fila , esta matriz se transforma a la forma escalonada de fila . Entonces, tiene la siguiente forma:

(do1,1do1,2do1,metrodoq,1doq,2doq,metro000d1,1d1,2d1,metro000d,1d,2d,metro000000000000){\displaystyle {\begin{pmatrix}c_{1,1}&c_{1,2}&\cdots &c_{1,m}&\bullet &\bullet &\cdots &\bullet \\\vdots &\vdots &&\vdots &\vdots &\vdots &&\vdots \\c_{q,1}&c_{q,2}&\cdots &c_{q,m}&\bullet &\bullet &\cdots &\bullet \\0&0&\cdots &0&d_{1,1}&d_{1,2}&\cdots &d_{1,m}\\\vdots &\vdots &&\vdots &\vdots &\vdots &&\vdots \\0&0&\cdots &0&d_{\ell ,1}&d_{\ell ,2}&\cdots &d_{\ell ,m}\\0&0&\cdots &0&0&0&\cdots &0\\\vdots &\vdots &&\vdots &\vdots &\vdots &&\vdots \\0&0&\cdots &0&0&0&\cdots &0\end{pmatrix}}}

Aquí,{\displaystyle \bullet }representa números arbitrarios y los vectores (dopag,1,dopag,2,,dopag,metro){\displaystyle (c_{p,1},c_{p,2},\ldots ,c_{p,m})}por cadapag{1,,q}{\displaystyle p\in \{1,\ldots ,q\}}y(dpag,1,,dpag,metro){\displaystyle (d_{p,1},\ldots ,d_{p,m})}por cadapag{1,,}{\displaystyle p\in \{1,\ldots ,\ell \}}son distintos de cero.

Entonces(y1,,yq){\displaystyle (y_{1},\ldots ,y_{q})}con

yi:=j=1metrodoi,jBj{\displaystyle y_{i}:=\sum _{j=1}^{m}c_{i,j}B_{j}}

es una base deU+W{\displaystyle U+W} y(z1,,z){\displaystyle (z_{1},\ldots ,z_{\ell })}con

zi:=j=1metrodi,jBj{\displaystyle z_{i}:=\sum _{j=1}^{m}d_{i,j}B_{j}}

es una base deUW{\displaystyle U\cap W}.

Prueba de corrección

Primero, definimosπ1:V×VV,(a,b)a{\displaystyle \pi _{1}:V\times V\to V,(a,b)\mapsto a}ser la proyección al primer componente.

Dejar H:={(,)U}+{(w,0)wW}V×V.{\displaystyle H:=\{(u,u)\mid u\in U\}+\{(w,0)\mid w\in W\}\subseteq V\times V.} Entoncesπ1(H)=U+W{\displaystyle \pi _{1}(H)=U+W}y H(0×V)=0×(UW){\displaystyle H\cap (0\times V)=0\times (U\cap W)}.

También,H(0×V){\displaystyle H\cap (0\times V)}es el núcleo deπ1|H{\displaystyle {\pi _{1}|}_{H}}, la proyección restringida a H. Por lo tanto,oscuro(H)=oscuro(U+W)+oscuro(UW){\displaystyle \dim(H)=\dim(U+W)+\dim(U\cap W)}.

El algoritmo de Zassenhaus calcula una base de H. En las primeras m columnas de esta matriz, hay una base.yi{\displaystyle y_{i}}deU+W{\displaystyle U+W}.

Las filas del formulario(0,zi){\displaystyle (0,z_{i})}(conzi0{\displaystyle z_{i}\neq 0}) están obviamente enH(0×V){\displaystyle H\cap (0\times V)}. Debido a que la matriz está en forma escalonada por filas , también son linealmente independientes. Todas las filas que son diferentes de cero ((yi,){\displaystyle (y_{i},\bullet )}y(0,zi){\displaystyle (0,z_{i})}) son una base de H , por lo que hayoscuro(UW){\displaystyle \dim(U\cap W)}semejantezi{\displaystyle z_{i}}s. Por lo tanto, elzi{\displaystyle z_{i}}s forman una base deUW{\displaystyle U\cap W}.

Ejemplo

Consideremos los dos subespaciosU=(1101),(0011){\displaystyle U=\left\langle \left({\begin{array}{r}1\\-1\\0\\1\end{array}}\right),\left({\begin{array}{r}0\\0\\1\\-1\end{array}}\right)\right\rangle }yW=(5033),(0532){\displaystyle W=\left\langle \left({\begin{array}{r}5\\0\\-3\\3\end{array}}\right),\left({\begin{array}{r}0\\5\\-3\\-2\end{array}}\right)\right\rangle }del espacio vectorialR4{\displaystyle \mathbb {R} ^{4}}.

Utilizando la base estándar , creamos la siguiente matriz de dimensión(2+2)×(24){\displaystyle (2+2)\times (2\cdot 4)}:

(11011101001100115033000005320000).{\displaystyle \left({\begin{array}{rrrrrrrr}1&-1&0&1&&1&-1&0&1\\0&0&1&-1&&0&0&1&-1\\\\5&0&-3&3&&0&0&0&0\\0&5&-3&-2&&0&0&0&0\end{array}}\right).}

Utilizando operaciones elementales de fila , transformamos esta matriz en la siguiente matriz:

(10000101001100001101){\displaystyle \left({\begin{array}{rrrrrrrrr}1&0&0&0&&\bullet &\bullet &\bullet &\bullet \\0&1&0&-1&&\bullet &\bullet &\bullet &\bullet \\0&0&1&-1&&\bullet &\bullet &\bullet &\bullet \\\\0&0&0&0&&1&-1&0&1\end{array}}\right)}(Algunas entradas han sido reemplazadas por "{\displaystyle \bullet }"porque son irrelevantes para el resultado."

Por lo tanto ((1000),(0101),(0011)){\displaystyle \left(\left({\begin{array}{r}1\\0\\0\\0\end{array}}\right),\left({\begin{array}{r}0\\1\\0\\-1\end{array}}\right),\left({\begin{array}{r}0\\0\\1\\-1\end{array}}\right)\right)}es una base deU+W{\displaystyle U+W}, y ((1101)){\displaystyle \left(\left({\begin{array}{r}1\\-1\\0\\1\end{array}}\right)\right)}es una base deUW{\displaystyle U\cap W}.

Véase también

Referencias

  1. Luks, Eugene M.; Rákóczi, Ferenc; Wright, Charles RB (abril de 1997), "Algunos algoritmos para grupos de permutación nilpotentes", Journal of Symbolic Computation , 23 (4): 335–354 , doi : 10.1006/jsco.1996.0092.
  2. ^ Fischer, Gerd (2012), Lernbuch Lineare Algebra und Analytische Geometrie (en alemán), Vieweg+Teubner , págs. 207–210 , doi : 10.1007/978-3-8348-2379-3 , ISBN  978-3-8348-2378-6
  3. The GAP Group (13 de febrero de 2015), "24 Matrices" , Manual de Referencia GAP, Versión 4.7 , consultado el 11 de junio de 2015.
  • «Mathematik-Online-Lexikon: Zassenhaus-Algorithmus» (en alemán) . Consultado el 15 de septiembre de 2012 .