PRINCIPIOS DE ÁLGEBRA LINEAL



John Jairo García Mora

Sonia Jaquelliny Moreno Jiménez

Margarita Patiño Jaramillo

Fondo Editorial Pascual Bravo

Título de la obra: Principios de Álgebra lineal
PRINCIPIOS DE ÁLGEBRA LINEAL


Autores:
John Jairo García Mora
Sonia Jaquelliny Moreno Jiménez
Margarita Emilia Patiño Jaramillo



Diseño del contenido: John Jairo García Mora
Revisión del texto: Juan Guillermo Rivera Berrío
Actualización: Joel Espinosa Longi
Código JavaScript para el libro: Joel Espinosa Longi, IMATE, UNAM.
Recursos interactivos: DescartesJS
Fuentes: Lato y UbuntuMono
Fórmulas matemáticas: KaTeX\KaTeX
Núcleo del libro interactivo: julio 2022



Fondo Editorial Pascual Bravo Calle 73 73A-226
PBX: (574) 4480520
Apartado 6564
Medellín, Colombia

ISBN: 978-958-56858-8-8



LICENCIA

Creative Commons Attribution License 4.0.

Prefacio

El álgebra de matrices es una herramienta básica para casi todas las ramas de la matemática, así como para disciplinas afines tales como la física, la ingeniería y la computación, y temas como el encriptamiento de mensajes, el control analógico y digital, solucionar ecuaciones lineales, entre otras.

Las notas no presuponen ningún conocimiento previo de álgebra lineal, pero si requiere de conocimiento de propiedades de los números reales y estructuras que provienen de la aritmética elemental.

Uno de los problemas que soluciona el algebra lineal, es la solución de sistemas de ecuaciones. Pese a que existen métodos iterativos, como el de Gauss, su cálculo manual se vuelve muy laborioso y más expuesto a cometer errores. Una alternativa, es usar la matriz inversa para calcular el vector de incógnitas de las ecuaciones.

Desde este punto de vista, en este texto calcularemos matrices inversas y resolveremos sistemas lineales de ecuaciones ya sea por la denominada matriz inversa o el método de eliminación de Gauss.

xi

Algunas de las escenas interactivas fueron adaptaciones de las escenas que aparecen en el proyecto Un1001, dirigido por José Luis Abreu León de la UNAM.

Los vídeos son autoría propia de los autores realizados con licencia propia del software Camtasia versión 2019 de la empresa TechSmith.

Las imágenes empleadas fueron obtenidas por suscripción a Pngtree2, entre ellas podemos citar la imagen que nos indica la finalización de cada capítulo tal como se puede observar en la siguiente página.

En cada capítulo aparecen las ecuaciones numeradas con un carácter en mayúsculas y un digito.

Por último, aparecen las referencias bibliográficas en las cuales se apoyo el presente libro de nociones de Álgebra Lineal.

xii



Ecuaciones
Sistemas de ecuaciones

¿QUÉ ES UNA ECUACIÓN? ¿QUÉ ES UN SISTEMA DE ECUACIONES?

Una ecuación es una afirmación de que dos expresiones son iguales, cuando se igualan entre sí dos expresiones, y al menos una de ellas contiene una variable, entonces la proposición matemática es una ecuación en una variable, si tiene dos variables diferentes se denomina ecuación en dos variables y así sucesivamente.

Una solución o raíz de una ecuación es cualquier número que, al reemplazar la(s) variable(s), la convierte en una proposición verdadera. Se dice que un número satisface una ecuación si es una solución de la ecuación. Resolver una ecuación significa hallar todas sus soluciones.

Solucionar una ecuación es hallar sus RAÍCES o CEROS.

Un sistema de ecuaciones algebraicas es un conjunto de dos o más ecuaciones con varias incógnitas que conforman un problema matemático que consiste en encontrar los valores de las incógnitas que satisfacen dichas operaciones. En un sistema de ecuaciones algebraicas las incógnitas son valores numéricos menores a la constante (o más generalmente elementos de un CUERPO sobre el que se plantean las ecuaciones.

17

Mientras que en una ecuación diferencial las incógnitas son funciones o distribuciones de un cierto conjunto definido de antemano. Una solución de dicho sistema es por tanto, un valor o una función que substituida en las ecuaciones del sistema hace que éstas se cumplan automáticamente sin que se llegue a una contradicción.

En álgebra abstracta, un CUERPO es una estructura algebraica en la cual las operaciones llamadas adición y multiplicación se pueden realizar y cumplen las propiedades: asociativa, conmutativa y distributiva de la multiplicación respecto de la adición, además de la existencia de inverso aditivo, de inverso multiplicativo y de un elemento neutro para la adición y otro para la multiplicación, los cuales permiten efectuar las operaciones de sustracción y división.

1.1 Ecuaciones lineales

El concepto de ecuación lineal se encuentra asociado a las rectas en el plano que tienen como expresión algebraica:

y=mx±b(E1)\tag{E1}y=mx\pm{b}

En donde se tiene:

  • mm es el valor de la pendiente, significa el valor de la inclinación con respecto al eje xx positivo.
  • bb es la ordenada en el origen (punto de intercepto con el eje yy).

Veamos algunos ejemplos de este tipo de ecuaciones lineales:

    y=x+7y=-x+7
  • Tiene un mm negativa de valor 11 e intercepta el eje de las ordenadas en (0,7)(0, 7).
18
    y=4x3y=4x-3
  • Tiene un mm positiva de valor 44 e intercepta el eje de las ordenadas en (0,3)(0,-3).
  • y=8y=8
  • Tiene m=0m = 0 es paralela al eje de las abscisas y pasa por el punto (0,8)(0, 8).
  • x=5x=5
  • No tiene mm (inclinación con respecto al eje xx) es paralela al eje de las ordenadas y pasa por el punto (5,0)(5, 0).

En este texto trabajaremos con este tipo de ecuaciones lineales: con dos variables, ecuaciones lineales con una sola variable implica que la otra variable vale CERO, como en los ejemplos 3 y 4.

1.2 Sistemas de Ecuaciones lineales

En general, un sistema con mm ecuaciones lineales y nn incógnitas puede ser escrito en forma normal como:

a11x1+a12x2+...+a1nxn=b1a21x1+a22x2+...+a2nxn=b2...+...+...+...=...am1x1+am2x2+...+a3nxn=bm(E2)\tag{E2}\begin{matrix}{a_{11}x_1}&{+}&{a_{12}x_2}&{+}&{...}&{+}&{a_{1n}x_n}&{=}&{b_1}\\ {a_{21}x_1}&{+}&{a_{22}x_2}&{+}&{...}&{+}&{a_{2n}x_n}&{=}&{b_2}\\ {...}&{+}&{...}&{+}&{...}&{+}&{...}&{=}&{...}\\ {a_{m1}x_1}&{+}&{a_{m2}x_2}&{+}&{...}&{+}&{a_{3n}x_n}&{=}&{b_m}\end{matrix}

Un sistema lineal de dos ecuaciones con dos incógnitas tiene la forma que podemos apreciar en la ecuación (3) de la siguiente página:

19
a11x+a12y=b1a21x+a22y=b2(E3)\tag{E3}\begin{matrix}{a_{11}x}&{+}&{a_{12}y}&{=}&{b_1}\\ {a_{21}x}&{+}&{a_{22}y}&{=}&{b_2}\end{matrix}

En la solución de este sistema se deben tener en cuenta algunas propiedades del álgebra:

  • Si a=ba = b y c=dc = d, entonces a+c=b+da + c = b + d. Ello se traduce en que al sumar dos ecuaciones se obtiene una tercera ecuación que cumple las características de las ecuaciones empleadas para sumar.
  • Si a=ba = b y cc es cualquier número real entonces ac=bcac = bc. Ello nos permite hallar ecuaciones equivalentes.

Al momento de hallar las RAÍCES de un sistema lineal 2×22×2 nos podemos encontrar que:

  • Posee una solución única. Las líneas rectas se interceptan en un punto de coordenadas (x,y)(x,y). En este caso el sistema de ecuaciones se denomina consistente y las ecuaciones son independientes
  • Posee infinitas soluciones. Cuando los coeficientes y el término independiente de una de las ecuaciones es múltiplo o submúltiplo de los coeficientes de la otra ecuación. Al igual que el sistema anterior es consistente pero las ecuaciones son dependientes, expresado de otra forma los pares de números reales corresponde a los puntos de una recta.
  • No tiene solución. Las ecuaciones representan rectas paralelas. El sistema se denomina inconsistente

En la escena interactiva de la siguiente página podemos visualizar estos resultados.Los anteriores resultados también ocurrir en el espacio tridimensional con las ecuaciones de tres planos de ese espacio.

20

Escena interactiva

Haz clic en el botón que aparece en la esquina superior derecha para ampliar e interactuar con el recurso.

1.3 Métodos de solución de un sistema 2×22×2

1.3.1. Solución por el método de igualación

Para ello seguimos los siguientes pasos empleando las propiedades de las igualdades:

21
  1. Despejamos la misma variable en ambas ecuaciones, de allí obtenemos las ecuaciones 3 y 4.
  2. Igualamos los resultados obtenidos.
  3. Solucionamos para la única variable de esa ecuación que resulta de esa igualación.
  4. El valor de la variable obtenida se reemplaza en cualquiera de las ecuaciones expresadas hasta el momento.
  5. Solucionamos para esa variable.
  6. Comprobamos los resultados en las ecuaciones originales.

Dadas las ecuaciones aa y bb siguientes hallar los valores de xx e yy para que se conviertan en una identidad:

(a)5x+3y=113(b)x7y=5\begin{matrix} {(a)}&{5x}&{+}&{3y}&{=}\frac{11}{3} \\ {(b)}&{-x} & {-}&{7y}&{=}{-5} \end{matrix}

Paso 1: despejar la misma variable en ambas ecuaciones:

De (a)(a) despejamos xx:

5x+3y3y=1133y5x=1133y5x+3y-3y=\frac{11}{3}-3y\Longrightarrow5x=\frac{11}{3}-3y 5x5=1133y5x=1133y5x=119y15\frac{5x}{5}=\frac{\frac{11}{3}-3y}{5}\Longrightarrow x=\frac{\frac{11}{3}-3y}{5}\Longrightarrow x=\frac{11-9y}{15}

De (b)(b) despejamos xx:

x+x7y=x57y=x5-x+x-7y=x-5\Longrightarrow-7y=x-5
22
7y+5=x5+5x=57y-7y+5=x-5+5\Longrightarrow x=5-7y

Paso 2: igualamos esos resultados obtenidos:

119y15=57y119y15(15)=(57y)(15)\frac{11-9y}{15}=5-7y\Longrightarrow\frac{11-9y}{15}(15)=(5-7y)(15) 119y=75105y119y+105y=75105y+105y11-9y=75-105y\Longrightarrow 11-9y+105y=75-105y+105y 11+96y=751111+96y=751196y=6411+96y=75\Longrightarrow 11-11+96y=75-11\Longrightarrow96y=64

Paso 3: despejar para esa variable.

96y96=6496y=6496y=23\frac{96y}{96}=\frac{64}{96}\Longrightarrow y=\frac{64}{96}\Longrightarrow y=\frac{2}{3}

Paso 4: sustituimos el valor de esa variable en cualquiera de las ecuaciones del paso 1 y solucionamos para esa variable (paso 5).

x=5723x=5143x=13x=5-7*\frac{2}{3}\Longrightarrow x=5-\frac{14}{3}\Longrightarrow x=\frac{1}{3}

Paso 6: sustituimos el valor de las variables en cualquiera de las ecuaciones dadas inicialmente.

513+323=?53+63=?113=1135*\frac{1}{3}+3*\frac{2}{3}=?\Longrightarrow \frac{5}{3}+\frac{6}{3}=?\Longrightarrow \frac{11}{3}=\frac{11}{3}

1.3.2. Solución por el método de sustitución

Para ello seguimos los siguientes pasos empleando las propiedades de las igualdades:

  1. Despejamos una variable en una de las ecuaciones.
  2. Reemplazamos esa dato en la OTRA ecuación.
23
  1. Solucionamos para la única variable de esa ecuación.
  2. El valor de la variable obtenida se reemplaza en cualquiera de las ecuaciones expresadas hasta el momento.
  3. Solucionamos para esa variable.

Sean las ecuaciones:

(a)11x+14y=20(b)3x10y=36\begin{matrix} {(a)}&{11x}&{+}&{14y}&{=}-20 \\ \\ {(b)}&{3x} & {-}&{10y}&{=}\quad{36} \end{matrix}

Paso 1: despejamos una de las variables en cualquiera de las ecuaciones, de (b) despejaremos la variable xx.

3x10y+10y=36+10y3x=36+10y3x-10y+10y=36+10y\Longrightarrow 3x=36+10y 3x13=(36+10y)13x=36+10y3x=12+10y33x*\frac{1}{3}=(36+10y)*\frac{1}{3}\Longrightarrow x=\frac{36+10y}{3}\Longrightarrow x=12+\frac{10y}{3}

Paso 2: el valor obtenido lo reemplazamos en la ecuación 1 (La OTRA ecuación):

11(12+10y3)+14y=2011\left(12+\frac{10y}{3}\right)+14y=-20 132+1103y+14y=20132132+1523y=20132132+\frac{110}{3}y+14y=-20\Longrightarrow 132-132+\frac{152}{3}y=-20-132 1523y=152\frac{152}{3}y=-152
24

Paso 3. Despejamos la variable de la ecuación obtenida.

1523y=1521523y3152=1523152y=3\frac{152}{3}y=-152\Longrightarrow\frac{152}{3}y*\frac{3}{152}=-152*\frac{3}{152} \Longrightarrow y=-3

Paso 4. El valor de la variable y se reemplaza en la ecuación obtenida en el paso 1 y resolvemos (Paso 5).

x=12+103(3)=1210=2x=12+\frac{10}{3}*{(-3)}\Longrightarrow=12-10=2

Paso 6. Comprobamos en la ecuación (a) los resultados obtenidos:

11(2)+14(3)=?2242=?20=2011(2)+14(-3)=?\Longrightarrow 22-42=?\Longrightarrow -20=-20

Podemos entonces afirmar que el par ordenado (2,3)(2, -3) es el punto común a las dos rectas citadas como se observa en la figura construida con GeoGebra:

25

Interactividad: solución por sustitución

En la siguiente escena aparecen dos ecuaciones, haz clic en el botón paso a paso para ver el proceso. También encontraras dos botones de color rojo para que selecciones la variable a despejar y un botón con el rótulo "Cambiar": con él cambiarás la ecuación en la cual se ha de despejar la variable.

1.3.3. Solución por reducción o eliminación

Se trata de reducir el sistema a una sola ecuación:

26
  1. Sumamos las ecuaciones para determinar si una de las variables se cancela.
  2. De no suceder lo descrito en el paso anterior, procedemos a multiplicar las ecuaciones por números enteros de tal manera que una de esas variables se cancele.
  3. Solucionamos para la única variable de esa ecuación obtenida.
  4. El valor de la variable obtenida se reemplaza en cualquiera de las ecuaciones expresadas hasta el momento.
  5. Solucionamos para esa variable.

Dadas las ecuaciones aa y bb siguientes hallar los valores de pp y de qq para que se conviertan en una identidad:

(a)4x+3y=10(b)5x6y=7\begin{matrix} {(a)}&{4x}&{+}&{3y}&{=}\quad{10} \\ \\ {(b)}&{5x} & {-}&{6y}&{=}{-7} \end{matrix}

Paso 1: sumamos las ecuaciones para saber si alguna de las variables se cancela:

4x+3y=105x6y=7\begin{matrix}&{4x}&{+}&{3y}&{=}&{10}\\ &{5x}&{-}&{6y}&{=}&{-7}\\ \hline \end{matrix} 9x3y=3\begin{matrix}{9x}&{-}&{3y}&{=}&{3}\end{matrix}

Paso 2. Como la suma continúa con las dos variables procedemos a multiplicar ambas ecuaciones por números enteros de tal manera que al sumar se consiga lo propuesto.

27

Lo ideal es aprovechar la variable de ambas ecuaciones que posean signo diferente, para nuestro caso observamos que la variable yy nos facilita el propósito. Procedemos a multiplicar por números enteros

(a)(2)8x+6y=20(a)*(2)\quad{8x}\quad{+}\quad{6y}\quad{=}\quad{20} (b)(1)5x6y=7(b)*(1)\quad{5x}\quad{-}\quad{6y}\quad{=}\quad{-7} 8x+6y=205x6y=7\begin{matrix}&{8x}&{+}&{6y}&{=}&{20}\\ &{5x}&{-}&{6y}&{=}&{-7}\\ \hline \end{matrix} 13x+0y=13\begin{matrix}{13x}&{+}&{0y}&{=}&{13}\end{matrix}

Ahora solucionamos para la única variable en la ecuación:

13x=1313x(113)=13(113)x=113x=13\Longrightarrow 13x\left(\frac{1}{13}\right)=13\left(\frac{1}{13}\right)\Longrightarrow x=1

El próximo paso es reemplazar la variable xx por su valor en cualquiera de las ecuaciones originales para determinar el valor de la otra variable:

5(1)6y=756y=76y=755\left(1\right)-6y=-7\Longrightarrow 5-6y=-7\Longrightarrow -6y=-7-5 6y=126y=126y(16)=12(112)-6y=-12\Longrightarrow -6y=-12 \Longrightarrow -6y\left(\frac{1}{6}\right)=-12\left(\frac{1}{12}\right) y=2y=2

Ahora comprobemos:

4(1)+3(2)=?4+6=1010=104(1)+3(2)=?\Longrightarrow 4+6=10 \Longrightarrow 10=10
28

Vídeo

Este método a veces requiere que se complemente con herramientas algebraicas: en el vídeo observaremos la solución de un sistema 2×22×2 que contiene las variables en el denominador de una expresión racional.

Práctica con los ejercicios de la siguiente página.

Utiliza el método propuesto en cada numeral.

29

1.4 Sistema de ecuaciones lineales 3×33×3

Interpretación geométrica de este sistema

Un sistema en el espacio tridimensional R3ℝ^3, es una ecuación de la forma:

Ax+By+Cz+D=0(E4)\tag{E4} Ax + By + Cz + D = 0

En esa expresión los coeficientes AA, BB, CC y DRD \in ℝ y al menos uno de ellos es de CERO y representa un plano. En un sistema de dos ecuaciones en tres variables se pueden presentar las siguientes situaciones:

30
  1. Dos planos que se cortan formando una línea recta, en ese caso y dado que una recta tiene infinitos puntos, el sistema tiene infinitas soluciones.
  2. Dos planos paralelos formando un sistema sin solución.
  3. Un único plano, es también un sistema con infinitas soluciones dado que el plano tiene infinitos puntos.

Un sistema de tres ecuaciones o 3×33×3 representa las siguientes posibilidades:

  1. Cuando el sistema tiene solución única, son tres planos que se cortan en un solo punto.
  2. Tres planos paralelos cuando el sistema no tiene solución.
  3. Dos planos coincidentes en uno solo y que es paralelo al tercer plano y el sistema no tiene solución.
  4. Tres planos coincidentes, el sistema tiene infinitas soluciones, tantas como tiene uno solo de los planos.
  5. Dos planos paralelos y un tercer plano que los corta en dos rectas paralelas, entonces el sistema no tiene solución.
  6. Tres planos que se cortan en una recta común, podemos expresar que el sistema tiene infinitas soluciones como lo son los puntos de la recta.

En la escena interactiva de la siguiente página se puede observar el comportamiento básico de esos planos.

31

1.5 El método de Gauss

Una generalización del método de reducción, utilizado para resolver sistemas de dos ecuaciones con dos incógnitas, es el método de Gauss. El método consiste en triangular el sistema, y para conseguirlo, utilizamos los criterios de equivalencia que detallaremos a continuación, de manera que consigamos hacer 0 los valores que están debajo de los elementos de la diagonal principal.

SISTEMA EQUIVALENTE. Se dice que dos sistemas de ecuaciones son equivalentes cuando tienen las mismas soluciones.

32

Expresado de otra manera, son aquellos sistemas donde toda solución del primero es solución del segundo y viceversa. Si dos ecuaciones lineales son equivalentes entonces tienen el mismo número de incógnitas aunque no necesariamente igual número de ecuaciones. Es evidente que si se cambia el orden de las ecuaciones, el sistema resultante no sólo es equivalente al inicial, sino que es el mismo.

Lo anterior ya lo habíamos practicado en el método de reducción, allí nos dimos cuenta que si en un sistema lineal de ecuaciones sustituimos una ecuación cualquiera por la ecuación que resulta al multiplicarla por un número cualquiera no nulo y sumarla, miembro a miembro, a otras ecuaciones del sistema (después de multiplicar estas últimas por un número cualquiera) el nuevo sistema lineal que obtenemos es equivalente al primero (tiene las mismas soluciones que el primero).

CRITERIOS DE EQUIVALENCIA, cuando determinamos sistemas de ecuaciones equivalentes debemos tener en cuenta:

Si se multiplican los dos miembros de una ecuación de un sistema por un número distinto de cero, resulta otro sistema equivalente al dado.

Si a una ecuación de un sistema se le suma o resta otra ecuación del mismo, resulta otro sistema equivalente.

OBSERVACIONES

  • Multiplicar una ecuación por cero equivale a suprimirla.
  • El 1º criterio se utiliza para conseguir que los coeficientes de una incógnita, en dos ecuaciones, sean iguales en valor absoluto (salvo el signo) y así poder eliminarla sumando ecuaciones (2º criterio).
33
  • Si los números son enteros conviene elegir como coeficiente común de una incógnita el mcm (mínimo común múltiplo) de los coeficientes para que los cálculos sean más sencillos.
  • Si al sumar ecuaciones de un sistema resulta una ecuación incompatible, el sistema también lo es.

Es de anotar que si al aplicar el método de Gauss que veremos a continuación llegamos a una ecuación del tipo 0x+0y+0z=00x + 0y + 0z = 0, se suprime. Si quedan menos ecuaciones que incógnitas, el sistema tiene infinitas soluciones y se conoce como sistema indeterminado. Veamos el esquema triangular que debemos obtener por el método de Gauss:

Esquema de reducción de Gauss.

Solucionemos el sistema aplicando el método de Gauss.

{1)2x3y+5z=62)4x+2y+12z=163)3x7y4z=4\begin{cases} 1)\qquad 2x - 3y + 5z & = & -6 \\ 2)\qquad 4x + 2y + 12z & = & 16\\ 3)\qquad 3x - 7y - 4z & = & 4\\ \end{cases}

Ahora hallaremos ecuaciones equivalentes entre la 1ª la 3ª de tal manera que al sumarlas el valor de 3x3x se convierta en 0x0x.

Ello implica que debemos determinar el mcm entre los coeficientes de la variable xx de ambas ecuaciones:

34

{1)(3)3)(2){1)6x+9y15z=183)6x14y8z=8R)0x5y23z=26\begin{cases} 1)*(-3)\\ 3)*(2)\\ \\ \end{cases}\begin{cases} 1)\qquad -6x + 9y - 15z & = & 18 \\ 3)\qquad \quad 6x - 14y - 8z & = & 8\\ R)\qquad \quad {\color{red}0x} - 5y - 23z & = & 26\\ \end{cases}

Reescribimos nuestro sistema:

{1)2x3y+5z=62)4x+2y+12z=163)0x5y23z=26\begin{cases} 1)\qquad 2x - 3y + 5z & = & -6 \\ 2)\qquad 4x + 2y + 12z & = & 16\\ 3)\qquad {\color{red}0x} - 5y - 23z & = & 26\\ \end{cases}

Nuestro siguiente paso es hallar ecuaciones equivalentes entre la 1ª la 2ª de tal manera que al sumarlas el valor de 4x4x se convierta en 0x0x. Sumaremos la ecuación 1 multiplicada por (2-2) y la ecuación 2 multiplicada por 11.

{1)(2)2)(1){1)4x+6y10z=123)4x+2y+12z=16R)0x+8y+2z=28\begin{cases} 1)*(-2)\\ 2)*(1)\\ \\ \end{cases}\begin{cases} 1)\qquad -4x + 6y - 10z & = & 12 \\ 3)\qquad \quad 4x + 2y + 12z & = & 16\\ R)\qquad \quad {\color{red}0x} + 8y + 2z & = & 28\\ \end{cases}

Reescribimos nuestro sistema:

{1)2x3y+5z=62)0x+8y+2z=283)0x5y23z=26\begin{cases} 1)\qquad 2x - 3y + 5z & = & -6 \\ 2)\qquad {\color{red}0x} + 8y + 2z & = & 28\\ 3)\qquad {\color{red}0x} - 5y - 23z & = & 26\\ \end{cases}

Nuestro siguiente paso es hallar ecuaciones equivalentes entre la 2ª la 3ª de tal manera que al sumarlas el valor de 5y-5y se convierta en 0y0y. Sumaremos la ecuación 2 multiplicada por 55 y la ecuación 3 multiplicada por 88.

35
{2)(5)3)(8){1)0x+40y+10z=1403)0x40y184z=208R)0x+0y174z=348\begin{cases} 2)*(5)\\ 3)*(8)\\ \\ \end{cases}\begin{cases} 1)\qquad 0x + 40y + 10z & = & 140 \\ 3)\qquad 0x - 40y - 184z & = & 208\\ R)\qquad {\color{red}0x} + {\color{red}0y} - 174z & = & 348\\ \end{cases}

Reescribimos nuestro sistema:

{1)2x3y+5z=62)0x+8y+2z=283)0x+0y174z=348\begin{cases} 1)\qquad 2x - 3y + 5z & = & -6 \\ 2)\qquad {\color{red}0x} + 8y + 2z & = & 28\\ 3)\qquad {\color{red}0x} + {\color{red}0y} - 174z & = & 348\\ \end{cases}

Tomemos la tercera ecuación y solucionemos:

174z(1174)=348(1174)z=2 -174z\left(\frac{1}{-174}\right)=348\left(\frac{1}{-174}\right)\Longrightarrow {\Large{z}={\color{red}-2}}

Ahora reemplazamos z=2z=-2 en la ecuación 2:

8y+2(2)=288y4=288y4+4=28+48y + 2(-2)=28 \Longrightarrow 8y-4=28 \Longrightarrow 8y-4+{\color{red}4}=28+{\color{red}4} 8y=328y(18)=32(18)y=48y=32 \Longrightarrow 8y\left(\frac{1}{8}\right)=32 \left(\frac{1}{8}\right)\Longrightarrow {\Large{y}={\color{red}4}}

Por último, reemplacemos y=4y=4 y z=2z=-2 en la ecuación 1:

2x3(4)+5(2)=62x1210=62x-3(4)+5(-2)=-6 \Longrightarrow 2x-12-10=-6 2x22=62x22+22=6+222x=16 \Longrightarrow 2x-22=-6 \Longrightarrow 2x-22+{\color{red}22}=-6+{\color{red}22} \Longrightarrow 2x=16 2x=162x(12)=16(12)x=82x=16 \Longrightarrow 2x\left(\frac{1}{2}\right)=16 \left(\frac{1}{2}\right)\Longrightarrow {\Large{x}={\color{red}8}}

El sistema planteado tiene solución única: (8,4,2)(8, 4, 2)

36


En los numerales 1 a 10 utilice el método de Gauss para solucionar el sistema dado:

1){x+y+2z=192x+4y3z=13x+6y=5z2){2x+y=zx+3y+2z=5x+4y+9=y 1) \begin{cases} x + y + 2z = 19 \\ 2x + 4y - 3z = 1\\ 3x + 6y = 5z\\ \end{cases}\qquad \qquad 2) \begin{cases} 2x + y = z \\ x + 3y + 2z = 5\\ x + 4y + 9 = y\\ \end{cases}
37
3){8a+bc=24a+3b+2c=114a3b5c=264){5x+2y+2z=36x+5y+4z=457x7y+7z=21 3) \begin{cases} 8a + b - c = 2 \\ 4a + 3b + 2c = 11\\ -4a - 3b - 5c = -26\\ \end{cases}\qquad \qquad 4) \begin{cases} -5x + 2y + 2z = -3 \\ 6x + 5y + 4z = 45\\ 7x - 7y + 7z = 21\\ \end{cases} 5){2x7y+5z=13x4y+2z=135x4y+7z=66){5m4n+7p=383m+3n+3p=32m10n5p=17 5) \begin{cases} 2x - 7y + 5z = 1 \\ 3x - 4y + 2z = \frac{1}{3}\\ 5x - 4y + 7z = 6\\ \end{cases}\qquad \qquad \qquad 6) \begin{cases} 5m - 4n + 7p = -38\\ 3m + 3n + 3p = -3\\ 2m - 10n - 5p = -17\\ \end{cases} 7){2a+3b+4c=44b+7c=a+235a+b+5c=38){5m5n+5p=1543m+n2p=144m+3n3p=234 7) \begin{cases} 2a + 3b + 4c = 4 \\ 4b + 7c = a + 23\\ 5a + b + 5c = 3\\ \end{cases}\qquad \qquad \qquad 8) \begin{cases} 5m - 5n + 5p = - \frac{15}{4}\\ -3m + n - 2p = \frac{1}{4}\\ 4m + 3n - 3p = \frac{23}{4}\\ \end{cases} 9){4a+3b3c=47ab12c=1202a+2b5c=610){2h+2j5k=104h+3j2k=78h+jk=15 9) \begin{cases} 4a + 3b - 3c = 47\\ a - b - 12c = -120\\ 2a + 2b - 5c = -6\\ \end{cases}\qquad \qquad \qquad 10) \begin{cases} 2h + 2j - 5k = 10\\ -4h + 3j - 2k = 7\\ 8h + j - k = -15\\ \end{cases}
38



Vectores
Cantidades vectoriales

Vectores - Definición

Los vectores y la historia.

El Álgebra Lineal estudia los conjuntos denominados espacios vectoriales o espacios euclidianos. Un espacio euclidiano es el conjunto de nn-adas ordenadas, también conocido por espacio nn-dimensional y se denota por Rnℝ^n este es una sucesión de n números reales ejemplo (a1,a2,...,an)(a_1,a_2,...,a_n) donde los vectores Rnℝ^n se clasifican así:

  • R1ℝ^1, es el espacio en una dimensión: línea recta.
  • R2ℝ^2, es el espacio bidimensional representado por pares ordenados.
  • R3ℝ^3, es el espacio tridimensional representado por ternas ordenadas.
  • Rnℝ^n, es el espacio en n dimensiones y representado por nn-adas ordenadas.

Un vector es una forma de representar determinados objetos en el plano y en el espacio, es un conjunto ordenado de números escritos en forma de renglones o en forma de columnas.

43

Un Vector renglón de nn componentes es un conjunto ordenado de nn números escritos así:

(x1,x2,...,xn)(V1)\large {\tag{V1}(x_1, x_2, ... , x_n)}

2.1 Interpretación geométrica y algebraica de vector en R2ℝ^2

Sean AA y BB dos puntos cualesquiera en el plano, el segmento de recta dirigido de AA hasta BB y el conjunto de todos los segmentos equivales a él se denominan vectores. Desde el punto de vista algebraico, un vector en el plano es el par ordenado de números reales (a,b)(a, b). Los números aa y bb se denominan componentes del vector:

AB=v=(a,b)(V2)\large {\tag{V2}\overline{AB}=\overrightarrow{v}=(a, b)}

En la gráfica, podemos observar que no obstante los segmentos de recta posean igual punto de inicio e igual punto final, la dirección del vector es diferente.

Vectores de igual magnitud, dirección y sentidos opuestos.
44

2.2 Componentes rectangulares de un vector

Dado que un vector es la representación del conjunto de segmentos de rectas equivalentes, se define como módulo, magnitud o medida del vector a la longitud de cualquiera de esos segmentos con igual dirección y sentido.

Componentes rectangulares de un vector.

Ese módulo es la resultante de aplicar el denominado Teorema de Pitágoras, puesto que las componentes del vector son las proyecciones del mismo en los ejes coordenados. Esas componentes se obtienen de la siguiente manera:

  • La componente en xx es la diferencia entre la coordenada xx del punto de aplicación (punto final del vector) y la coordenada xx del punto de inicio.
  • La componente en yy es la diferencia entre la coordenada yy del punto de aplicación (punto final del vector) y la coordenada yy del punto de inicio.
45

El módulo de un vector vv es un escalar y se representa por: v(V3)\large{\tag{V3}|v|}

La dirección de esos vectores equivalentes es la dirección del vector equipolente (vector que tiene su punto de inicio en el origen de coordenadas del plano cartesiano), esta dado por el ángulo cuya tangente es la razón entre la componente en y dividida por la componente en xx. Es el ángulo que forma con el eje xx positivo.

v=[x2x1,y2y1]=[](V4)\tag{V4}{\large{\overrightarrow{v}=[ x_2-x_1,y_2-y_1 ]=[ \, ]}} v=vx2+vy2(V5)\tag{V5}{\large|\overrightarrow{v}|=\sqrt{{\overrightarrow{v_x}}^2+{\overrightarrow{v_y}}^2}} θ=Tan1(vyvx)(V6)\tag{V6}{\theta=Tan^{-1}\left(\frac{\overrightarrow{v_y}}{\overrightarrow{v_x}}\right)}

Un vector uu tiene su punto final en el cuarto cuadrante del plano cartesiano y sus coordenadas están representadas por el par ordenado (4,6)(-4,-6) y su punto de origen está ubicado en las coordenadas (8,1).(8,-1).

Determinar el valor de las componentes rectangulares, su magnitud y su dirección.

  1. Determinar y graficar el vector equipolente.
  2. u=[(48),(61)]=[12,5]u = [(-4-8), (-6--1)] = [-12, -5]
46
Vector equipolente.
  1. Hallemos el valor de la magnitud de ese vector:
  2. v=(12)2+(5)2=169=13U|\overrightarrow{v}|=\sqrt{{(-12)}^2+{(-5)}^2}=\sqrt{169}=13U
  3. Determinemos el valor del ángulo, teniendo en cuenta que el vector está en el IV cuadrante y es necesario adicionarle π\pi (180°180°).
  4. θ=Tan1(512)0.394+π3.534202.62°\theta = Tan^{-1}\left(\frac{-5}{-12}\right)\cong 0.394+\pi \cong 3.534 \cong 202.62°
  5. Hallemos las componentes rectangulares a partir del siguiente grafico:
Ángulo resultante.
47

Podemos observar que es necesario recurrir a otras relaciones trigonométricas del triángulo rectángulo: el seno y el coseno, veamos:

θ=Sen1(vy13)vy=13Senθ4.97U\theta = Sen^{-1}\left(\frac{\overrightarrow{v_y}}{13}\right) ∴ \overrightarrow{v_y}= 13*Sen \theta \cong -4.97 U θ=Cos1(vx13)vx=13Cosθ12.01U\theta = Cos^{-1}\left(\frac{\overrightarrow{v_x}}{13}\right) ∴ \overrightarrow{v_x}=13*Cos \theta \cong -12.01 U

Para un vector w\overrightarrow{w} dado determinar:

  1. El valor de la coordenada xx del punto de aplicación (x,3)(x,3) para que el valor de su módulo sea igual a 157\sqrt{157} si las coordenadas del punto de partida son (1,3)(1,-3)
  2. La dirección de ese vector w.\overrightarrow{w}.

Coordenada en xx:

w=(x1)2+(3(3))2=157|\overrightarrow{w}|=\sqrt {(x-1)^2+(3-(-3))^2}=\sqrt{157} (x1)2+(6)2=157x22x+1+36=157(x-1)^2+(6)^2=157 \Longrightarrow x^2-2x+1+36=157 x22x+37157=157157x22x120=0x^2-2x+37-{\color{red}{157}}=157-{\color{red}{157}} \Longrightarrow x^2-2x-120=0

Por fórmula general:

x=(2)±(2)24(1)(120)2(1)x=12ox=10x=\frac{-(-2)\pm \sqrt{(-2)^2-4(1)(-120)}}{2(1)}\Longrightarrow x=12 \quad o \quad x=-10
48

La dirección del vector w\overrightarrow{w}

θ=Tan1(611)=28.61°\theta= Tan^{-1}\left(\frac{6}{11}\right)= 28.61°

Interactivo

Luego de recordar la información haz clic en el botón Cerrar y podrás visualizar un nuevo vector desplazando los puntos azules.

49

  • En los numerales 1 - 6 determine las coordenadas del punto final o del punto de inicio según el caso y grafique el vector.
    1. u=[3,2]u = [3, 2]; punto inicial en (5,2)(-5, 2)
    2. u=[5,2]u = [-5, -2]; punto inicial en (2,1)(2, -1)
    3. u=[5/2,5/2]u = [5/2, 5/2]; punto inicial en (5/2,7/2)(5/2, -7/2)
    4. u=[0,14]u = [0, 14]; punto final en (6,6)(-6, 6)
    5. u=[13,7]u = [13, -7]; punto final en (3/4,3)(3/4, -3)
    6. u=[3/3,2/5]u = [√3/3, 2/√5]; punto final en (1,1)(1, 1)
  • En los numerales 7 - 14 determine el valor del módulo de los vectores dados cuyos puntos final e inicial se dan:
    1. punto inicial en (4,2)(4, -2); punto final en (5,2)(-5, 2)
    2. punto inicial en (8,5)(-8, 5); punto final en (5,8)(5, -8)
    3. punto inicial en (3/4,1/2)(3/4, -1/2); punto final en (7,12)(7, 12)
    4. punto inicial en (0.75,1)(0.75, 1); punto final en (2.5,2.25)(-2.5, 2.25)
    5. punto inicial en (2,7)(-2, -7); punto final en (0,22/5)(0, 22/√5)
    6. punto inicial en (2/6,3/2)(2/√6, 3/2); punto final en (5/3,2)(-5/3, -2)
    7. punto inicial en (1/2,1/7)(1/2, -1/7); punto final en (3,22/8)(3, 22/√8)
    8. punto inicial en (4/5,3/2)(4/√5, 3/2); punto final en (5/3,20)(-5/3, -20)
  • En los numerales 15 - 21 determine el valor del ángulo de los vectores dados cuyos puntos final e inicial se dan:
    1. punto inicial en (4,2)(4, -2); punto final en (5,2)(-5, 2)
    2. punto inicial en (8,5)(-8, 5); punto final en (5,8)(5, -8)
    3. punto inicial en (3/4,1/2)(3/4, -1/2); punto final en (7,12)(7, 12)
    4. punto inicial en (0.75,1)(0.75, 1); punto final en (2.5,2.25)(-2.5, 2.25)
50
  1. punto inicial en (2,7)(-2, -7); punto final en (0,22/5)(0, 22/√5)
  2. punto inicial en (2/6,3/2)(2/√6, 3/2); punto final en (5/3,2)(-5/3, -2)
  3. punto inicial en (1/2,1/7)(1/2, -1/7); punto final en (3,22/8)(3, 22/√8)

2.3 Multiplicación de un vector por un escalar

Simplemente es un vector cuyas componentes son proporcionales al escalar dado y se presentan las siguientes opciones:

  • Un vector de menor magnitud, igual dirección e igual sentido si el escalar es mayor de CERO y menor que 11.
  • Un vector de mayor magnitud, igual dirección e igual sentido si el escalar es mayor que 11.
  • Un vector de menor magnitud, su dirección aumenta 180°180° y sentido contrario si el escalar dado es menor que CERO y mayor de 1-1.
  • Un vector de mayor magnitud, su dirección se incrementa 180°180° y sentido contrario si el escalar dado es menor de 1-1.

Una fuerza equivalente a 200200 Newton se aplica a un objeto formando un ángulo de 30°30° con la horizontal, determinar el incremento o disminución de sus componentes si la esa fuerza se duplica. ¿Será que la magnitud de esas componentes también se duplica?

  1. Determinemos el valor de las componentes originales:
  2. Sen30°=Fy200NFy=200NSen30°Fy=100NSen 30°=\frac{\vec{F_y}}{200N} ∴ {\vec{F_y}}=200N*Sen 30° ∴{\vec{F_y}}=100N
51
Cos30°=Fx200NFx=200NCos30°Fx=173.20NCos 30°=\frac{\vec{F_x}}{200N} ∴ {\vec{F_x}}=200N*Cos 30° ∴{\vec{F_x}}=173.20N
  1. Determinemos las componentes si multiplicamos la magnitud del vector por el escalar 22: 2200Newton=400Newton2*200 Newton = 400 Newton.
  2. Sen30°=Fy400NFy=400NSen30°Fy=200NSen 30°=\frac{\vec{F_y}}{400N} ∴ {\vec{F_y}}=400N*Sen 30° ∴{\vec{F_y}}=200N Cos30°=Fx400NFx=400NCos30°Fx=346.40NCos 30°=\frac{\vec{F_x}}{400N} ∴ {\vec{F_x}}=400N*Cos 30° ∴{\vec{F_x}}=346.40N

2.4 El producto escalar

El producto escalar de dos vectores es un número real que resulta al multiplicar el producto de sus módulos por el coseno del ángulo que forman, es también conocido como producto interno o producto punto.

Geométricamente, es el producto de dos magnitudes euclidianas de dos vectores u=(a1,b1)u = (a_1, b_1) y v=(a2,b2)v = (a_2, b_2) por el coseno del ángulo θθ formado por esos dos vectores:

uv=uvCosθ(V8)\tag {V8} \large{uv = |u||v|Cos θ}

Analíticamente el producto punto se expresa como:

uv=a1a2+b1b2(V2)\tag {V2}\large {uv = a_1a_2 + b_1b_2}

θθ es el ángulo no negativo más pequeño de esos dos vectores equipolentes o sea que tienen el origen de coordenadas como su punto inicial y estará comprendido en el intervalo [0,π][0, π].

52

Según esta definición, podemos calcular el ángulo

θ=Cos1(uvuv)(V9)\tag{V9} θ=Cos^{-1}\left(\frac{uv}{|u||v|}\right)

En la gráfica vemos el vector u=(4,6)u = (4,-6) y el vector v=(3,6)v = (3,6), determinar el ángulo θθ formado por ellos:

Ángulo entre dos vectores
  1. Calculemos inicialmente el producto punto entre los vectores:
  2. uv=(4)(3)+(6)(6)=24uv=(4)(3)+(6)(-6)=-24
  3. La magnitud de uu
  4. u=(4)2+(6)2=52|u|=\sqrt{(4)^2+(-6)^2}=\sqrt{52}
53
  1. El módulo de vv
  2. u=(3)2+(6)2=45|u|=\sqrt{(3)^2+(6)^2}=\sqrt{45}
  3. Aplicamos la fórmula para el ángulo
  4. θ=Cos1(245245)=Cos1(242340) θ=Cos^{-1}\left(\frac{-24}{|\sqrt{52}||\sqrt{45}|}\right)=Cos^{-1}\left(\frac{-24}{|\sqrt{2340}|}\right) θ=Cos1(0.496831) θ=Cos^{-1}(-0.496831) θ=119.74°=2.089Rad θ=119.74°=2.089 Rad

2.5 Vectores paralelos

Dado que dos vectores uu y vv son paralelos si el ángulo entre ellos es 00 o π\pi, además de ello tienen la misma dirección o direcciones opuestas.

Vamos a determinar las coordenadas (5,b1)(-5, b1) del vector uu de u=29|u|= \sqrt{29} para que sea paralelo al vector v=(10,4)v = (10,4).

  1. El punto de partida es la fórmula:
  2. Cos(θ)=uvuvCos(\theta)=\frac{uv}{|u||v|}
  3. Ese coseno debe ser igual a 11 (00) o 1-1 (π\pi), por lo tanto:
  4. 1=uvuv1=uvuv\large {1=\frac{uv}{|u||v|}\quad \vee \quad 1=\frac{uv}{|u||v|}}
54
  1. En el denominador nos hace falta el módulo de vv.
  2. v=(10)2+(4)2=116|v|=\sqrt{(10)^2+(4)^2}=\sqrt{116}
  3. Reemplazamos v|v| en el denominador:
  4. 1=uv291161=uv29116 {1=\frac{uv}{|\sqrt{29}||\sqrt{116}|}\quad \vee \quad -1=\frac{uv}{|\sqrt{29}||\sqrt{116}|}} 1=uv58C1=uv58 {1=\frac{uv}{58}C -1=\frac{uv}{58}}
  5. Eliminemos denominadores en ambas igualdades:
  6. 1(58)=uv58(58)1(58)=uv58(58) 1(58)=\frac{uv}{58}(58)\quad \vee \quad -1(58)=\frac{uv}{58}(58) 58=uv58=uv 58=uv\quad \vee \quad -58=uv
  7. Con los datos conocidos reemplacemos el producto punto
  8. 58=(10)(5)+4b158=50+4b158=(10)(-5)+4b_1\quad ∴ \quad 58=-50+4b_1 \vee 58=(10)(5)+4b158=50+4b1-58=(10)(-5)+4b_1\quad ∴ \quad -58=-50+4b_1
  9. Despejamos b1b_1:
  10. (58+50)(14)=(50+50+4b1)(14) (58+50)\left(\frac{1}{4}\right)=(-50+50+4b_1)\left(\frac{1}{4}\right) 27=b1 27=b_1
55
(58+50)(14)=(50+50+4b1)(14) (-58+50)\left(\frac{1}{4}\right)=(-50+50+4b_1)\left(\frac{1}{4}\right) 2=b1 -2=b_1
  1. Verificamos con cuál de esos valores se obtiene la misma dirección del vector vv
  2. Tan(v)=410=25Tan(v)=275Tan(v)=\frac{4}{10}=\frac{2}{5}\vee Tan(v)=\frac{27}{-5}

Observamos que la respuesta correcta es b1=2\large {b_1=-2}

2.6 Vectores ortogonales o perpendiculares

Dos vectores uu y vv diferentes de CERO, son ortogonales si el ángulo entre ellos es π/2π/2, ello equivale a decir que el coseno del ángulo vale CERO o que el producto punto vale CERO.

2.7 Proyección de un vector sobre otro vector

Proyección de un vector sobre otro vector
56

Sean uu y vv dos vectores diferentes de CERO, entonces la proyección de uu sobre vv es un vector denominado como:

proyvu=uvv2v(V10)\tag{V10}\large {proy_vu=\frac{u \cdot v}{|v|^2}\cdot v}

Interactivo

Mueve los puntos AA y BB para que puedas visualizar la proyección de uu sobre vv.

57

Vamos a determinar las coordenadas del vector que proyecta el vector u=(6,5)u = (-6,5) sobre el vector v=(2,2)v = (2,-2)

  1. Paso 1. Determinar el módulo del vector vv:
  2. v=(2)2+(2)2=8|v|=\sqrt{(2)^2+(-2)^2}=\sqrt{8}
  3. Paso 2. Hallar el valor del producto escalar de los vectores
  4. uv=(6)(2)+(5)(2)=22u \cdot v=(-6)(2)+(5)(-2)=-22
  5. Aplicar la fórmula dada en la definición:
  6. proyvu=2282[2,2]=114[2,2] {proy_vu=\frac{-22}{|\sqrt{8}|^2}\cdot [2,2]=-\frac{11}{4} \cdot [2,2]} proyvu=[112,112]=[5.5,5.5] proy_vu=[-\frac{11}{2},\frac{11}{2}]=[-5.5, 5.5]

2.8 Vectores en R3ℝ^3

Un vector en el espacio tridimensional analizado desde el punto de vista algebraico es una terna o triada de números reales (x,y,z)(x, y, z) donde los números xx, yy y zz se denominan componentes del vector.

2.9 Módulo de un vector en R3ℝ^3

Dados los puntos P=(x1,y1,z1)P = (x_1, y_1, z_1) y Q=(x2,y2,z2)Q = (x_2, y_2, z_2), podemos definir las coordenadas del vector uu que va desde PP hasta QQ como:

58
u=[(x2x1),(y2y1),(z2z1)]u=[(x_2-x_1),(y_2-y_1),(z_2-z_1)]

Al igual que los vectores en R2ℝ^2 el módulo será la distancia entre en punto inicial y el punto final del vector:

u=(x2x1)2+(y2y1)2+(z2z1)2|u|=\sqrt{(x_2-x_1)^2+(y_2-y_1)^2+(z_2-z_1)^2}

Interactivo

Ahora realizaremos algunos cálculos con respecto a los vectores en el espacio tridimensional.

59

Calcular el módulo del vector uu que tiene su origen en P=(3,2,4)P = (-3, 2, 4) y va hasta el punto Q=(4,0,4)Q = (4, 0, -4):

(4(3))2+(02)2+((4)4)2\sqrt{(4-(3))^2+(0-2)^2+((-4)-4)^2} (7)2+(2)2+(8)2=161\sqrt{(7)^2+(-2)^2+(-8)^2}=\sqrt{161}

Hemos calculado la norma del vector u=[7,2,8]u=[7, -2,-8]

2.10 Dirección de un vector en R3ℝ^3

La disposición más común para representar un vector en el espacio sigue la denominada regla de la mano derecha, donde los dedos extendidos señalan la dirección positiva de los ejes coordenados.

Disposición de la regla de la mano derecha
60

Los ejes coordenados así dispuestos determinan tres planos coordenados y se denominan plano xyxy (nuestro plano habitual en dos dimensiones), el plano xzxz y el plano yzyz.

Planos coordenados en el espacio

2.11 Cosenos directores

Los ángulos del vector tridimensional
61

Los vectores en el espacio tridimensional forman un ángulo con cada uno de los ejes coordenados.

En una base ortonormal (Una base es ortogonal si los vectores de la base son perpendiculares entre sí), se llaman cosenos directores del vector u=[ux,uy,uz]u = [u_x, u_y, u_z], a los cosenos de los ángulos que forma el vector u con los vectores de la base.

Dado el vector u=[ux,uy,uz]u = [u_x, u_y, u_z], cada una de esas componentes forma un ángulo θ\theta que se determina por:

θx=Cos1(uxu);θy=Cos1(uyu)(V11)\tag{V11} \theta_x=Cos^{-1}\left(\frac{u_x}{|u|}\right);\theta_y=Cos^{-1}\left(\frac{u_y}{|u|}\right)θz=Cos1(uzu)\theta_z=Cos^{-1}\left(\frac{u_z}{|u|}\right)

Con esos cosenos directores se cumple que:

Cos2(θx)+Cos2(θy)+Cos2(θz)=1(V12)\tag{V12} Cos^2(\theta_x)+Cos^2(\theta_y)+Cos^2(\theta_z)=1

Determinar los ángulos que el vector u=[1,4,3]u = [-1, 4, 3] forma con los ejes coordenados.

  1. El primer paso es reconocer o calcular (cuando se conocen las ternas de los puntos de inicio y final del vector) las componentes ux=1u_x = -1, uy=4u_y = 4, uz=3u_z = 3
  2. Paso 2: calcular la magnitud del vector:
  3. u=sqrt(1)2+(4)2+(3)2=26|u|=sqrt{(-1)^2+(4)^2+(3)^2}=\sqrt{26}
62
  1. Paso 3. Calcular los ángulos:
  2. θx=Cos1(126)=101.30°\theta_x=Cos^{-1}\left(\frac{-1}{\sqrt{26}}\right)=101.30° θy=Cos1(426)=38.32°\theta_y=Cos^{-1}\left(\frac{4}{\sqrt{26}}\right)=38.32° θz=Cos1(326)=53.96°\theta_z=Cos^{-1}\left(\frac{3}{\sqrt{26}}\right)=53.96°
  3. Verifiquemos la propiedad que expresa que la suma de los cuadrados de los cosenos directores es igual a UNO.
  4. (126)2+(326)2+(426)2=\left(\frac{-1}{\sqrt{26}}\right)^2+\left(\frac{3}{\sqrt{26}}\right)^2+\left(\frac{4}{\sqrt{26}}\right)^2= (126)+(926)+(1626)=2626=1\left(\frac{1}{26}\right)+\left(\frac{9}{26}\right)+\left(\frac{16}{26}\right)=\frac{26}{26}=1

2.12 Vectores unitarios

En R2ℝ^2 existen dos vectores que nos permiten especificar cualquier vector en función de ellos, son vectores de magnitud 11, el primero de ellos es el vector (1,0)(1,0) que se denota por ii. El otro vector es el vector (0,1)(0,1) denotado por jj. Así, el vector v=(3,7)v = (3,7) puede expresarse como v=3i+7jv = 3i + 7j.

De igual manera en el espacio tridimensional existen tres vectores unitarios que nos permiten representar cualquier vector:

63
  1. El vector (1,0,0)(1, 0 ,0) denotado por ii
  2. El vector (0,1,0)(0, 1 ,0) denotado por jj
  3. El vector (0,0,1)(0, 0, 1) denotado por kk

En la gráfica se puede ver más detenidamente las descripción y posición de los vectores unitarios.

Vectores unitarios en el espacio tridimensional

Sea uu un vector unitario y vv otro vector diferente de CERO, entonces podemos determinar un vector unitario que tiene la misma dirección de vv mediante:

u=vv(V13)\tag{V13} u=\frac{v}{|v|}

Determinar los vectores unitarios que tengan la misma dirección de:

a) mm = (5,  6,  4)(-5,\; 6,\; 4),\quad b) nn = (0,  3,  7)(0,\; -3,\; 7) y,\quad c) pp = (18,38,58)\large {(\frac{1}{8},-\frac{3}{8},\frac{5}{8})}.

64

a)m=(5)2+(6)2+(4)2=77\qquad |m|=\sqrt{(-5)^2+(6)^2+(4)^2}=\sqrt{77}

u=[5,6,4]77=[577,677,477]\large {u=\frac{[-5,6,4]}{\sqrt{77}}}=\left[-\frac{5}{\sqrt{77}},\frac{6}{\sqrt{77}},\frac{4}{\sqrt{77}}\right]

b)n=(0)2+(3)2+(7)2=58\qquad |n|=\sqrt{(0)^2+(-3)^2+(7)^2}=\sqrt{58}

u=[0,3,7]58=[0,358,758]\large {u=\frac{[0,-3,7]}{\sqrt{58}}}=\left[0,-\frac{3}{\sqrt{58}},\frac{7}{\sqrt{58}}\right]

c)p=(18)2+(38)2+(58)2=3564=358\qquad |p|=\sqrt{(\frac{1}{8})^2+(-\frac{3}{8})^2+(\frac{5}{8})^2}=\sqrt{\frac{35}{64}}=\frac{\sqrt{35}}{8}

u=835[18,38,58]=[135,335,535]\large {u=\frac{8}{\sqrt{35}}} \cdot \left[\frac{1}{8},-\frac{3}{8},\frac{5}{8}\right]=\left[\frac{1}{\sqrt{35}},-\frac{3}{\sqrt{35}},\frac{5}{\sqrt{35}}\right]

2.13 Álgebra de vectores

Suma de vectores

Sean los vectores u=a1i+b1j+c1zu = a_1i + b_1j + c_1z y v=a2i+b2j+c2zv = a_2i + b_2j + c_2z, el vector suma al que llamaremos r=u+vr = u + v es igual a la suma de las componentes respectivas de los vectores uu y vv, se expresa como: d=(a1+a2)i+(b1+b2)j+(c1+c2)zd = (a_1 + a_2)i + (b_1+b_2)j + (c_1+c_2)z. Otra forma de expresar esa suma es =(a1+a2,b1+b2,c1+c2) = (a_1 + a_2, b_1 + b_2 , c_1 + c_2).

Diferencia de vectores

Sean los vectores u=a1i+b1j+c1zu = a_1i + b_1j + c_1z y v=a2i+b2j+c2zv = a_2i + b_2j + c_2z, el vector diferencia al que llamaremos d=uvd = u - v (será la suma de uu con el opuesto de vv) será igual a la diferencia de las componentes

65

respectivas de los vectores uu y vv, y se expresa como: d=(a1a2)i+(b1b2)j+(c1c2)zd = (a_1 - a_2)i + (b_1 - b_2)j + (c_1 - c_2)z. Otra forma de expresar esa diferencia es d=(a1a2,b1b2,c1c2)d = (a_1 - a_2, b_1 - b_2 , c_1 - c_2).

Dado que el método gráfico no posee la precisión ideal para realizar estas operaciones, ni con el denominado cabeza con cola y el método del paralelogramo, lo ideal es realizar la suma algebraica de los componentes rectangulares de todos los vectores involucrados en la operación.

Dados el vector u=3i4j+2ku = 3i -4j + 2k y el vector v=i+j2kv = -i + j -2k, determinar el módulo del vector resultante de la siguiente operación: 3u+2v3u + 2v.

  1. Multiplicar cada vector por el escalar indicado
  2. 3u=3[3i4j+2k]=[9i12j+6k]3u=3[3i -4j + 2k]=[9i-12j+6k] 2v=2[i+j2k]=[2i+2j4k]2v=2[ -i + j -2k]=[-2i+2j-4k]
  3. Realizar la operación algebraica entre los componentes del paso anterior
  4. 3u+2v=[92]i+[12+2]j+[64]k=7i10j+2k3u+2v=[9-2]i+[-12+2]j+[6-4]k=7i-10j+2k
  5. Determinar el valor exacto del módulo del vector resultante
  6. 3u+2v=(7)2+(10)2+(2)2=153|3u+2v|=\sqrt{(7)^2+(-10)^2+(2)^2}=\sqrt{153}
66

En la escena se encuentran tres vectores que puedes mover a placer con el desplazamiento de sus puntos de inicio y de aplicación. A la derecha de los vectores podrás observar unas diferentes operaciones de suma y resta que podemos realizar con ellos.

2.14 Producto Vectorial o producto cruz

También conocido como producto interno es una operación que se define solo en R3ℝ^3 puesto que a diferencia del producto punto cuyo resultado es un escalar, el resultado aquí es otro vector. Para el vector u=a1i+b1j+c1ku = a_1i + b_1j + c_1k y el vector u=a2i+b2j+c2ku = a_2i + b_2j + c_2k, el producto vectorial está definido como:

67

uu x vv = (b1c2c1b2)i+(c1a2a1c2)j+(a1b2b1a2)k(V11)(b_1c_2 - c_1b_2)i + (c_1a_2 - a_1c_2)j + (a_1b_2 - b_1a_2)k \quad {(V11)}

El resultado es un vector perpendicular a los vectores que se multiplican, y por lo tanto normal al plano que los contiene. Debido a su capacidad de obtener un vector perpendicular a otros dos vectores, cuyo sentido varía de acuerdo al ángulo formado entre estos dos vectores, esta operación es aplicada con frecuencia para resolver problemas matemáticos, físicos o de ingeniería como el torque de una fuerza, la velocidad angular y la dirección del campo magnético entre algunas de ellas.

Producto vectorial

La dirección de esa resultante está dada por la regla de la mano derecha, si denominamos el producto vectorial como el vector ww, entonces uu, vv y ww forman (en ese orden) la tríada de la mano derecha. Si hacemos rotar el vector uu en sentido contrario a las manecillas del reloj hasta que sea colineal con el vector vv. la resultante es el vector ww, en caso contrario de hacer rotar el vector vv en el mismo sentido de las manecillas del reloj hasta que sea colineal con uu, la resultante sería w.-w.

u×v=(u×v)(V14)\tag{V14} u \times v=-(u\times v)
68

Si θ\theta es el ángulo entre los vectores uu y vv (siempre menor o igual a 180°180°), entonces la magnitud de su producto cruz esta dado por:

u×v=uvSen(θ)(V15)\tag{V15} |u \times v| = |u| |v| Sen(\theta)

Los productos vectoriales de los vectores unitarios ii, jj, kk están dados por:

i×i=0i×j=ki×k=ji \times i = 0 \qquad i \times j = k \qquad i \times k = - j j×j=0j×i=kj×k=ij \times j = 0 \qquad j \times i = - k \qquad j \times k = i k×k=0k×i=jk×j=ik \times k = 0 \qquad k \times i = j \qquad k \times j = - i

Considerando que los vectores uu y vv tienen el mismo origen, existe la interpretación desde la geometría elemental que son los lados adyacentes de un paralelogramo.

Determinar el ángulo que forman y el área del paralelogramo cuyos lados adyacentes son los vectores u=[7,2,0]u = [7,2,0] y v=[6,7,0]v = [6,7,0]. Dado que los vectores se encuentran en el plano xyxy, el producto vectorial será en la dirección kk. Recordemos que:

θ=Cos1(uvuv)\theta =Cos^{-1}\left(\frac{u \cdot v}{|u||v|}\right)
  1. Calcular el producto punto entre ambos vectores
  2. uv=(7)(6)+(2)(7)=56u \cdot v=(7)(6)+(2)(7)=56
69
  1. Calcular los módulos de uu y de vv
  2. u=(7)2+(2)2=56|u|=\sqrt{(7)^2+(2)^2}=\sqrt{56} v=(6)2+(7)2=85|v|=\sqrt{(6)^2+(7)^2}=\sqrt{85}
  3. Calcular el ángulo
  4. θ=Cos1(565685)=Cos1(564505)\theta =Cos^{-1}\left(\frac{56}{|\sqrt{56}||\sqrt{85}|}\right)=Cos^{-1}\left(\frac{56}{\sqrt{4505}}\right) θ=Cos1(0.8343...)=33.45°\theta=Cos^{-1}(0.8343...)=33.45°
  5. Calcular el área del paralelogramo
  6. =[(2)(0)(0)(7)]i+[(0)(6)(7)(0)]j+[(7)(7)(2)(6)]k=[(2)(0)-(0)(7)]i+[(0)(6)-(7)(0)]j+[(7)(7)-(2)(6)]k u×v=[4912]k=37ku \times v=[49-12]k=37k Aˊrea=u×v=sqrt(37)2=37U2Área=|u \times v|=sqrt{(37)^2}=37 U^2

Consideremos ahora el producto cruz de uu y vv y el vector resultante ww, si esos vectores no se encuentran en el mismo plano, se convierten en los lados de un paralelepípedo en el espacio (volumen).

El cálculo de ese volumen se conoce como el triple producto escalar de uu y vv. Existe la condición de que No sean nulos, No paralelos de dos en dos dos y NO coplanares. Ese volumen es el resultado del producto punto de uno de los vectores por el producto escalar de los otros dos.

70
Volumen del paralelepípedo
Volumen=w(u×)v)U3(V16)\tag{V16} Volumen=|w \cdot (u \times) v)|U^3

Dados los vectores u=i,5j,3ku = i, -5j, -3k; v=i,j,2kv = -i, -j, 2k y w=iw = i, determinar si el volumen existe y su valor.

  1. Hallemos el producto vectorial del vector uu y el vector vv:
  2. u×v=[(5)(2)(3)(1)]iu \times v = [ (-5)(2) - (-3)(-1) ] i +[(3)(1)(1)(2)]j+[(1)(1)(5)(1)]k+ [ (-3)(-1) - (1)(2) ] j + [ (1)(-1) - (-5)(-1) ] k u×v=[103]i+[32]j+[15]ku \times v = [ -10 - 3 ] i + [ 3 - 2 ] j + [ -1 - -5 ] k u×v=13i+j6ku \times v = -13 i + j - 6 k
  3. Calculemos el valor absoluto del producto punto entre w=iw = i, 0j0j, 0k0k con el producto vectorial del paso anterior v=13i+j6kv = -13 i + j - 6 k
  4. w(u×v)=(1)(13)+(0)(1)+(0)(6)=13=13U3| w \cdot (u \times v ) | = | (1)(-13) + ( 0)(1) + ( 0)(-6) | = | -13 | = 13 U^3

71

Lo anterior es diferente al triple producto cruz de 3 vectores uu, vv y ww en R3ℝ^3 que está dado por:

u×(v×w)=(uw)v(uv)w(V17)\tag{V17} u \times ( v \times w) = (u \cdot w)v - (u \cdot v)w

Otros conceptos a considerar en los espacios vectoriales son 1) Dependencia lineal y, 2) Independencia lineal. Ello lo podemos expresar como: dos vectores en un espacio vectorial son linealmente dependientes y solo si uno de ellos es múltiplo escalar del otro.

Sean v1,v2,...,vnv_1, v_2, ..., v_n vectores en un espacio vectorial, se dice que esos vectores son linealmente dependientes si existen nn escalares c1,c2,...,cnc_1, c_2, ... , c_n no todos iguales a CERO y que cumplan:

c1v1+c2v2++cnvn=0(V18)\tag{V18} c_1v_1 + c_2v_2 + \dotso + c_nv_n = 0

  • Realizar las operaciones indicadas en los numerales 1-6 empleando los vectores:
  • u=[4,4,0],w=[4,3,2],a=[5,1,3]u=[4,4,0], w=[4,-3,2], a=[5,-1,3] b=[0,3,0],c=[4,3,0],e=[4,4,3]b=[0,-3,0], c=[4,3,0], e=[-4,4,3]
    1. 3w2c+b+a\quad 3w-2c+b+a
    2. ab+cd+2w\quad a-b+c-d+2w
    3. 3e+2w+2ce\quad 3e+2w+2c-e
    4. 5a+d2e+b\quad 5a+d-2e+b
    5. ed2wb\quad e-d-2w-b
    6. 2a+2d2e+2b\quad 2a+2d-2e+2b
72
  • En los numerales 7-14 determine si los vectores dados son ortogonales, paralelos o no poseen ninguna de esas condiciones.
    1. u=9i+4j;m=4i+9j\quad u=9i+4j;m=-4i+9j
    2. a=10j7k;d=6j3k\quad a=10j-7k;d=6j-3k
    3. e=i+5k;b=5i+k\quad e=i+5k;b=-5i+k
    4. h=4i3j;a=10i7j\quad h=4i-3j;a=10i-7j
    5. d=10j+3k;e=7j+4k\quad d=10j+3k;e=7j+4k
    6. a=6i2j;b=6i6j\quad a=6i-2j;b=6i-6j
    7. o=23j;p=7i\quad o=-23j;p=7i
    8. u=10i+3j;v=10i3j\quad u=10i+3j;v=-10i-3j
  • En los numerales 15-20 determine los vectores unitarios que poseen igual dirección que la resultante de los vectores indicados en cada operación.
    1. (3i3j+k)+(j+5k)\quad (3i-3j+k)+(j+5k)
    2. (10i7k)(i+6j3k)\quad (10i-7k)-(i+6j-3k)
    3. (ij2k)+(i+4k)\quad (-i-j-2k)+(i+4k)
    4. (2i2j+k)(ij+k)\quad (2i-2j+k)-(i-j+k)
    5. (3j)+(4i5j+2k)\quad (3j)+(4i-5j+2k)
    6. (2i+5k)+(2i5j)\quad (2i+5k)+(2i-5j)
  • En los numerales 21-26 determine el área del paralelogramo que forman los vértices adyacentes dados.
    1. (1,2,3);(2,1,0);(0,4,0)\quad (1,-2,3);(2,1,0);(0,4,0)
    2. (7,5,9);(3,6,5);(2,1,3)\quad (7,-5,9);(-3,-6,-5);(2,-1,-3)
73
    1. (4,8,10);(1,8,7);(5,7,5)\quad (4,8,10);(1,-8,-7);(-5,7,-5)
    2. (1,3,2);(1,2,1);(0,3,1)\quad (1,-3,2);(1,2,-1);(0,3,1)
    3. (1,2,2);(2,1,0);(0,0,0)\quad (1,2,-2);(-2,1,0);(0,0,0)
    4. (2,2,1);(2,1,0);(3,4,4)\quad (-2,-2,-1);(2,-1,0);(3,4,4)
  • En los numerales 27-32 determine el producto vectorial de los vectores dados.
    1. u=3i3j;v=j+5k\quad u=3i-3j;v=j+5k
    2. u=10i7k;v=i+6j\quad u=10i-7k; v=i+6j-
    3. u=i2k;v=i+4k\quad u=-i-2k;v=i+4k
    4. u=2i2j;v=ij\quad u=2i-2j; v=-i-j
    5. u=3j;v=5j+2k\quad u=3j;v=-5j+2k
    6. u=2i+5k;v=2i5j\quad u=2i+5k;v=2i-5j
  • En los numerales 33-38 determine la proyección del vector uu sobre el vv.
    1. u=3i3j;v=i+5j\quad u=3i-3j;v=-i+5j
    2. u=10i7j;v=i+6j\quad u=10i-7j; v=i+6j-
    3. u=5i2j;v=i+4j\quad u=-5i-2j;v=i+4j
    4. u=2i2j;v=i+3j\quad u=-2i-2j; v=-i-+3j
    5. u=3j;v=5j+2j\quad u=-3j;v=-5j+2j
    6. u=2i+5j;v=2i5j\quad u=2i+5j;v=2i-5j
74



Matrices
Tipos de matriz
Operaciones
Determinantes
Matrices: algo de historia.

Un concepto clave en álgebra lineal es el de matriz, definida como un arreglo bidimensional de números. Dado que puede definirse tanto la suma como el producto de matrices, en mayor generalidad se dice que son elementos de un anillo. Una matriz se representa por medio de una letra mayúscula (A,B,...A,B,...) y sus elementos con la misma letra en minúscula (a,b,...a,b,...), con un doble subíndice donde el primero indica la fila y el segundo la columna a la que pertenece.

Los elementos individuales de una matriz m×nm × n, a menudo denotados por aiai, ajaj, donde el máximo valor de sus elementos (i,j)(i,j) en ii es mm, y el máximo valor de jj es nn. Siempre que la matriz tenga el mismo número de filas y de columnas que otra matriz, estas se pueden sumar o restar elemento por elemento.

Las matrices se utilizan para múltiples aplicaciones y sirven, en particular, para representar los coeficientes de los sistemas de ecuaciones lineales que vimos en el capítulo anterior o para representar transformaciones lineales dada una base. En este último caso, las matrices desempeñan el mismo papel que los datos de un vector para las aplicaciones lineales.

Las matrices y los determinantes son herramientas del álgebra que facilitan el ordenamiento de datos, así como su manejo. Los conceptos de matriz y todos los relacionados fueron desarrollados básicamente en el siglo XIX por matemáticos como los ingleses J.J. Sylvester y

79

Arthur Cayley y el Irlandés William Hamilton. Las matrices se encuentran en aquellos ámbitos en los que se trabaja con datos regularmente ordenados y aparecen en situaciones propias de las Ciencias Sociales , Económicas y Biológicas.

Definición

Se denomina Matriz a una forma rectangular compuesta de filas (mm) y columnas (nn), donde cada elemento lleva asociado un doble índice, el primero (i=1,2,...,mi = 1, 2, ..., m) para indicar la fila y el segundo (j=1,2,...,nj = 1, 2, ..., n) para indicar la columna. El orden de la matriz (m×nm \times n) expresa las filas y columnas que componen la matriz. Se designa como m×nm \times n o como un conjunto de sus elementos (aijaij). Si I={1,2,...,m}I=\{1, 2, ..., m\} y J={1,2,...,n}J=\{1, 2, ..., n\}, una Matriz de mm filas y nn columnas y coeficientes reales se representa como M(m×n;R)M(m\times n; R).

3.1 Matriz sobre un cuerpo RR

Un cuerpo RR es un conjunto no vacío, junto con dos operaciones binarias internas que normalmente denominamos: 1°) Operación suma representada por ++ y, 2°) Operación producto representada por ·, operaciones que verifican las siguientes propiedades:

  1. Conmutatividad de la suma y el producto.
  2. Asociatividad de la suma y del producto.
  3. Existe un elemento neutro para la suma y un elemento neutro para el producto.
  4. Todos los elementos tienen un opuesto y todo elemento diferente del neutro de la suma tiene un inverso.
  5. El producto es distributivo por ambos lados con respecto a la suma.
80

3.2 Representación Matricial

A=(aij)=[a11a12a1na21a22a2nam1am2amn](M1) \tag{M1}A=(a_{ij})=\left[ \begin{array}{rrr} a_{11} & a_{12} & \dotsb & a_{1n}\\ a_{21} & a_{22} & \dotsb & a_{2n} \\ \vdots & \vdots & \dotsb & \vdots\\ a_{m1} & a_{m2} & \dotsb & a_{mn} \\ \end{array} \right]

3.3 Tipos de matriz

Vamos a construir una matriz AA de orden 3×33 \times 3 que cumpla las siguientes condiciones:

81
  1. aij=5aij = 5. Con la condición de que i=ji = j
  2. aij=0aij = 0. Con la condición de que iji ≠ j

Paso 1. Construimos una matriz genérica de 3×33 \times 3

A=[a11a12a13a21a22a23a31a32a33] A=\left[ \begin{array}{rrr} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \\ \end{array} \right]

Paso 2. Aplicamos la primera condición: aij=5a_{ij} = 5. Con la condición de que i=ji = j. Dicho de otra forma cambiaremos los elementos a11a_{11}, a22a_{22} y a33a_{33} por 55.

A=[5a12a13a215a23a31a325] A=\left[ \begin{array}{ccc} 5 && a_{12} && a_{13} \\ a_{21} && 5 && a_{23} \\ a_{31} && a_{32} && 5 \\ \end{array} \right]

Paso 3. Aplicamos la segunda condición de que iji ≠ j. Dicho de otra forma cambiaremos los elementos a12a_{12}, a13a_{13}, a21a_{21}, a23a_{23}, a31a_{31}, a32a_{32} por 00.

A=[500050005] A=\left[ \begin{array}{ccc} 5 & 0 & 0 \\ 0 & 5 & 0 \\ 0 & 0 & 5 \\ \end{array} \right]

Acabamos de construir una matriz escalar donde todos los elementos de la diagonal principal son iguales a 55.

  1. Construir una matriz 3×53 \times 5 donde todos los componentes sean CERO menos los que se detallan en la imagen:
82
  • Construir una matriz 4×44 \times 4 con las siguientes condiciones:
  • Los elementos de la diagonal principal sean UNO
  • Los demás componentes aija_{ij} sean las suma de i+ji + j

Operaciones con matrices

3.4 Suma de matrices.

La suma de dos matrices es otra matriz que resulta de sumar los elementos que corresponden a la misma fila y a la misma columna. Ambas matrices deben ser del mismo orden.

A=(aij),B=(bij)Mm,n;A = (a_{ij}), B = (b_{ij}) ∈ Mm, n ; A+B=C=(cij)Mm,nA + B = C = (c_{ij}) ∈ Mm, n tal que cij=aij+bij(M2)\tag{M2}c_{ij}=a_{ij}+b_{ij}

Sumar las matrices AA y BB que se presentan en la siguiente página.

83

3.5 Multiplicación una matriz por un escalar.

El producto de una matriz por un escalar es otra matriz que resulta de multiplicar el escalar por cada uno de los elementos de la matriz.

84

Dadas A=(aij)Mm,nyφKφA=(φaij)A = (a_{ij}) ∈ Mm, n y φ ∈ K → φA=(φa_{ij})

Sea φ=¾\large {φ = ¾}, determinar la matriz resultante de φA\large{φA}.

85

3.6 Multiplicación de matrices

Únicamente podemos multiplicar dos matrices AA y BB si el número de columnas de AA es igual al número de filas de BB, podemos expresar que AA y BB son compatibles bajo la multiplicación. De lo contrario se denominan incompatibles.

Sea A=aijA=a_{ij} un matriz m×nm \times n, y B=bijB=b_{ij} una matriz n×pn \times p, entonces el producto de AA y BB es una matriz m×pm \times p, C=cijC=c_{ij}, caracterizada por:

cij\large c_{ij} = (renglón ii de AA)*(columna jj de BB)

La siguiente escena es una adaptación del proyecto Un100.

86

Aunque esta aplicación puede realizarse sin emplear matrices, veremos su utilidad en la situación propuesta. Para una temporada escolar que se avecina, un vendedor al menudeo solicita cotización (en miles de unidades) a tres distribuidores mayoristas de los artículos que más le solicitan sus clientes y realiza la siguiente tabla expresada en dólares:

Vamos a determinar con cual de los mayoristas es más económico para el vendedor si compra 12 millares de cuadernos, 7 millares de colores, 5 millares de carpetas, 25 millares de lápices y 6 millares de borradores.

Sea AA la matriz de las cotizaciones y BB la matriz de la posible compra por el vendedor minorista:

AB=[305423012703154028107530738301385][1252567] AB=\left[ \begin{array}{ccccc} 305 & 42 & 30 & 12 & 70 \\ 315 & 40 & 28 & 10 & 75 \\ 307 & 38 & 30 & 13 & 85 \\ \end{array} \right] \left[ \begin{array}{c} 12 \\ 5 \\ 25 \\ 6 \\ 7 \\ \end{array} \right]

El costo de la compra con el mayorista AA:

(305)(12)+(42)(5)+(30)(25)+(12)(6)+(70)(7)=5182(305)(12)+(42)(5)+(30)(25)+(12)(6)+(70)(7)=5182

El costo de la compra con el mayorista BB:

(315)(12)+(40)(5)+(28)(25)+(10)(6)+(75)(7)=52265(315)(12)+(40)(5)+(28)(25)+(10)(6)+(75)(7)=52265
87

El costo de la compra con el mayorista CC:

(307)(12)+(38)(5)+(30)(25)+(13)(6)+(85)(7)=5297(307)(12)+(38)(5)+(30)(25)+(13)(6)+(85)(7)=5297

De acuerdo a lo anterior, para esa compra al vendedor minorista le es más económico comprarle al mayorista AA, un estudio completo implica los niveles de ganancia, el impuesto al valor agregado, entre muchas variables que no son competencia de este texto.

3.7 Traspuesta de una matriz

La siguiente escena es otra adaptación del proyecto Un100 descrito al inicio de este texto:

88

Si A=aijA = a_{ij} con m×nm \times n, entonces At=ajiA^t = a_{ji} con n×mn \times m. Es decir, la transpuesta de una matriz se obtiene intercambiando las filas y las columnas de esta.

Propiedades de la traspuesta de una matriz

3.8 Rango de una matriz

El propósito con las transformaciones elementales es convertir una matriz concreta en otra matriz mas fácil de estudiar.

89

En concreto, siempre será posible conseguir una matriz escalonada: "Sea AA una matriz y FF una fila de AA. Diremos que FF es nula si todos los números de FF coinciden con el cero.

Si FF es no nula, llamamos PIVOTE de FF al primer número distinto de cero de FF contando de izquierda a derecha en ese renglón. Una MATRIZ ESCALONADA es aquella que verifica las siguientes propiedades:

  1. El primer coeficiente no nulo de cada fila es 11, y se llama el pivote de la fila.
  2. Todas las filas nulas (caso de existir) se encuentran en la parte inferior de la matriz.
  3. El pivote de cada fila no nula se encuentra estrictamente más a la derecha que el pivote de la fila inmediatamente superior.

Una matriz escalonada se dice reducida si la columna arriba de cada pivote también es de ceros.

Veamos las siguientes matrices: AA no es escalonada, BB es escalonada por (2) y CC es escalonada por (3).

A=[100010013]  B=[102016000]  C=[132015] A=\left[ \begin{array}{rrr} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 1 & 3 \\ \end{array} \right] \; B=\left[ \begin{array}{rrr} 1 & 0 & 2 \\ 0 & 1 & 6 \\ 0 & 0 & 0 \\ \end{array} \right] \; C=\left[ \begin{array}{rrr} 1 & 3 & 2 \\ 0 & 1 & 5 \\ \end{array} \right]

En una matriz escalonada EE se define su RANGO como rg(E)rg(E) como el número de filas de esa matriz.

Podemos transformar una matriz no escalonada en una matriz escalonada mediante las siguientes transformaciones por filas:

90
  1. Intercambiar la posición de dos filas.
  2. Multiplicar una fila por un número real diferente de cero.
  3. Sustituir una fila por el resultado de sumarle a dicha fila otra fila que ha sido previamente multiplicada por un número cualquiera.

Partiendo de cualquier matriz AA se puede llegar, mediante una cantidad finita de transformaciones elementales, a una matriz escalonada EE. En el siguiente video se aplican algunas transformaciones:

91

Rango de una matriz

3.9 Determinante de una matriz

Para cada matriz cuadrada AA, podemos asociar un número llamado determinante de AA. Al hallar el determinante de una matriz se sustituyen los corchetes por barras verticales Se dice que un determinante de una matriz de n×nn \times n es un determinante de orden nn o un determinante de nn-ésimo orden. Se expresa como detAdet A ó A|A|.

Para una matriz cuadrada 2×22 \times 2 el determinante se calcula:

92
A=a11a12a21a22=(a11)(a22)(a12)(a21)(M3)\large {|A|=\tag {M3}\begin{vmatrix} a_{11} & a_{12}\\ a_{21} & a_{22} \end{vmatrix}=(a_{11})(a_{22})-(a_{12})(a_{21})}

Geométricamente el determinante se interpreta como el área orientada del paralelogramo que determinan los vectores fila linealmente independientes:

a=[a11,a12]b=[b11,b12]\large {\vec{a}=[a_{11},a_{12}]\qquad \vec{b}=[b_{11},b_{12}]}
Interpretación geométrica del determinante 2×22\times 2

Si A=12|A| = - 12, determinar el valor de a21a_{21} sabiendo que a11a_{11} tiene un valor igual a 5-5, a12a_{12} es igual 88 y a22a_{22} es igual a 1212

12=(5)(12)(8)(a12)\large {-12= (-5)(12)-(8)(a_{12})} 12=608(a12)12+60=60+608(a12)\large {-12= -60-8(a_{12}) ∴ -12+60=-60+60-8(a_{12})}
93
48=8(a12)48(18)=8(a12)(18)\large {48=-8(a_{12})\quad ∴ \quad 48\left(\frac{1}{-8}\right)=-8(a_{12})\left(\frac{1}{-8}\right)} 6=a12\large {-6=a_{12}}

Como ya conocemos la manera de hallar el determinante de una matriz cuadrada 2×22 \times 2 entenderemos con mayor propiedad el procedimiento que se describe a continuación y que explicaremos en el vídeo:

3.10 Determinante de una matriz 3×33 \times 3

94

Sea AA una matriz cuadrada de orden 33, se llama determinante de AA al número que se obtiene calculando la suma de los productos de los elementos de la diagonal principal menos los productos la diagonal secundaria.

Otra forma práctica de recordar la definición es la siguiente: Se escriben a la derecha (o debajo) de la matriz las dos primeras líneas. La diagonal principal y sus dos paralelas llevan el signo ++ (positivo), la diagonal secundaria y sus dos paralelas llevan el signo - (negativo).

3.11 Propiedades de los determinantes

  1. Si los elementos de una fila de una matriz de un determinante se pueden descomponer en dos sumandos, entonces su determinante es igual a la suma de dos determinantes que tienen en esa fila los primeros y segundos sumandos respectivamente, y en las demás los mismos elementos que el determinante inicial.
  2. Si se multiplican todos los elementos de una fila o de una columna por un número en una matriz cuadrada, el determinante queda multiplicado por ese número:
  3. det(KF1,F2,F3)=Kdet(F1,F2,F3)det (K · F_1, F_2, F_3) = K · det (F_1, F_2, F_3)
  4. Si tenemos dos matrices cuadradas A y B del mismo orden se cumple que:
  5. det(AB)=det(A)det(B)det (A · B) = det (A) · det (B)
  6. Si intercambiamos dos filas o dos columnas de una matriz cuadrada, su determinante cambia de signo con respecto al inicial:
95
    det(F1,F2,F3)=det(F1,F2,F3)det (F_1, F_2, F_3) = - det (F_1, F_2, F_3)
  1. Si una matriz cuadrada tiene una fila o una columna con todos sus elementos nulos, su determinante vale 00.
  2. det(0,C2,C3)=0det (0, C_2, C_3) = 0
  3. Si una matriz cuadrada tiene dos filas o dos columnas iguales, su determinante vale 00.
  4. det(F1,F1,F2)=0det (F_1, F_1, F_2)= 0

    Comprobemos calculado el determinante de la matriz AA:

    A=[245245020] A=\left[ \begin{array}{ccc} 2 & 4 & 5 \\ 2 & 4 & 5 \\ 0 & -2 & 0 \\ \end{array} \right]
  5. Si dos filas o dos columnas de una matriz cuadrada son proporcionales, su determinante será nulo:
  6. det(F1,KF1,F2)=0det (F_1, K · F_1, F_2)= 0
  7. Si una fila (columna) de una matriz cuadrada es una combinación lineal de las restantes filas (columnas), su determinante vale 00. Veamos el ejemplo:
  8. det(F1,F2,K1F1)+K2F2=0det (F_1, F_2, K_1 · F_1) + K_2 · F_2 = 0

Calcular el determinante de la matriz AA.

96
A=[24512181613] A=\left[ \begin{array}{rrr} 2 & 4 & 5 \\ 1 & 2 & -1 \\ 8 & 16 & 13 \\ \end{array} \right] 2[(2)(13)(1)(16)]4[(1)(13)(1)(8)]+5[(1)(16)(2)(8)]2[(2)(13)-(-1)(16)]-4[(1)(13)-(-1)(8)]+5[(1)(16)-(2)(8)] 2[26+16]4[13+8]+5[1616]2[26+16]-4[13+8]+5[16-16] 2[42]4[21]+5[0]8484=02[42]-4[21]+5[0]\Longrightarrow 84-84=0 A=0\large{|A|=0}
  1. Si a una fila o columna de una matriz cuadrada se le suma otra paralela , su determinante no varía.
  2. det(C1+C2,C2,C3)=det(C1,C2,C3)+det(C2,C2,C3)det (C_1+ C_2, C_2, C_3) = det (C_1, C_2, C_3) + det ( C_2, C_2, C_3) =det(C1,C2,C3) = det (C_1, C_2, C_3)
  3. Si a una fila o columna de una matriz cuadrada se le suma otra paralela multiplicada por un número, su determinante no varía.
  4. det(F1+KF2,F2,F3)=det(F1,F2,F3)+det(K1F2,F2,F3)det (F_1+K · F_2, F_2, F_3) = det (F_1, F_2, F_3) + det ( K_1 · F_2, F_2, F_3) =det(F1,F2,F3) = det (F_1, F_2,F_3)
  5. El determinante de una matriz y el de su traspuesta son iguales:
  6. A=At|A| = |A^t|
  7. El determinante de una suma de matrices no es la suma de sus determinantes:
97
A+B=/A+B|A + B| \quad{{=}\llap{/\,}}\quad |A| + |B|

Desarrollo de un determinante por los elementos de una línea.

Definición 1. Si A=(aij)A = (a_{ij}) es una matriz cuadrada de orden nn, se llama menor complementario del elemento aija_{ij}, y se representa por MijM_{i j}, al determinante de la submatriz que se obtiene al suprimir de AA la fila ii y la columna jj.

Nota. El adjunto de un elemento es igual a su menor complementario si la suma de subíndices es par, y a su opuesto si es impar. Ejemplo: El adjunto del elemento a21=1a_{21}=-1, de la matriz que vemos a continuación es A21=25A_{2 1} = -25.

A=[134101417]M21=[3417]=25 A=\left[ \begin{array}{rrr} 1 & 3 & 4 \\ -1 & 0 & 1 \\ 4 & 1 & -7 \\ \end{array} \right] ∴ M_{21}=\left[ \begin{array}{rrr} 3 & 4 \\ 1 & -7 \\ \end{array} \right]=-25

Definición 2. Se llama adjunto del elemento aija_{i j} a:

Aij=(1)i+jMij\large {A_{i j} = (-1)i + j M_{i j}}

Proposición 1. Un determinante es igual a la suma de los productos de los elementos de una línea cualquiera por sus respectivos adjuntos.

Desarrollo de un determinante en una matriz 4×44 \times 4 por adjuntos de una fila o de una columna.

98

Calcular el determinante de la matriz AA. Seleccionamos la fila o columna que incluya la mayor cantidad de elementos nulos:

A=[3311211110120201] A=\left[ \begin{array}{cccc} 3 & 3 & 1 & 1 \\ -2 & 1 & -1 & 1 \\ 1 & 0 & -1 & 2 \\ 0 & 2 & 0 & 1 \\ \end{array} \right]

Dada la matriz AA podemos visualizar que la fila 44 posee la mayor cantidad de elementos nulos, por lo tanto el cálculo de A|A| debe considerar el signo de sus adjuntos:

Adj=[++++++++] Adj=\left[ \begin{array}{cccc} + & - & + & - \\ - & + & - & + \\ + & - & + & - \\ - & + & - & + \\ \end{array} \right]

Calculemos el adjunto de a41a_{41} que debe ir precedido de signo negativo:

Adj(a41)=(-)(0)[311111012]=(-)(0)(6)=0 Adj(a_{41})=({\color{red}{\textbf{-}}})(0)\left[ \begin{array}{rrr} 3 & 1 & 1 \\ 1 & -1 & 1 \\ 0 & -1 & 2 \\ \end{array} \right]=({\color{red}{\textbf{-}}})(0)(-6)=0

Calculemos el adjunto de a42a_{42} que debe ir precedido de signo positivo:

Adj(a42)=(+)(0)[311211112]=(+)(2)(5)=10 Adj(a_{42})=({\color{red}{\textbf{+}}})(0)\left[ \begin{array}{rrr} 3 & 1 & 1 \\ -2 & -1 & 1 \\ 1 & -1 & 2 \\ \end{array} \right]=({\color{red}{\textbf{+}}})(2)(5)=10
99

Calculemos el adjunto de a43a_{43} que debe ir precedido de signo negativo:

Adj(a43)=(-)(0)[311211102]=(-)(0)(20)=0 Adj(a_{43})=({\color{red}{\textbf{-}}})(0)\left[ \begin{array}{rrr} 3 & 1 & 1 \\ -2 & 1 & 1 \\ 1 & 0 & 2 \\ \end{array} \right]=({\color{red}{\textbf{-}}})(0)(20)=0

Calculemos el adjunto de a44a_{44} que debe ir precedido de signo positivo:

Adj(a44)=(+)(0)[331211101]=(+)()(13)=13 Adj(a_{44})=({\color{red}{\textbf{+}}})(0)\left[ \begin{array}{rrr} 3 & 3 & 1 \\ -2 & 1 & -1 \\ 1 & 0 & -1 \\ \end{array} \right]=({\color{red}{\textbf{+}}})()(-13)=-13

El determinante de esa matriz es la suma de los valores calculados, he aquí la conveniencia de utilizar la fila o columna que mayor cantidad de valores nulos posea.

A=Adj(a41)+Adj(a42)+Adj(a43)+Adj(a44)|A|=Adj(a_{41})+Adj(a_{42})+Adj(a_{43})+Adj(a_{44}) A=0+10+013=3\large {|A|=0+10+0-13=-3}

3.12 Inversa de una matriz

Una matriz cuadrada AA, se dice invertible o que tiene matriz inversa, si existe una matriz cuadrada del mismo orden, que se denota por A1A^{-1}, tal que:

A1A=aA1=In(M4)\tag{M4} \large {A^{-1}* A = a * A^{-1} = I_{n}}

No toda matriz cuadrada es invertible. Cuando AA tiene inversa se dice que AA es regular; en caso contrario diremos que AA es singular.

100

Propiedades de la matriz inversa

  1. La inversa de A1A^{-1} es AA: (propiedad involutiva)
  2. La inversa de un producto de dos matrices es el producto de las inversas de esas matrices en orden inverso.
  3. (AB)1=B1A1A,BMm×n(A*B)^{-1} = B^{-1} * A^{-1} ∀ A,B ∈ M_{m \times n}; AA, BB Regulares.

  4. La matriz inversa solo existe si el determinante es diferente de CERO.

La inversa de una matriz 2×22 \times 2 por su definición y su equivalente práctica:

A1=1A(adjA)t=1A[a11a12a21a22](M5) \tag {M5} A^{-1}=\frac{1}{|A|}(adj A)^t=\frac{1}{|A|}\left[ \begin{array}{rrr} a_{11} & -a_{12} \\ -a_{21} & a_{22} \\ \end{array} \right]

A=[abcd]A1=1adbc[dbca] A=\left[ \begin{array}{cc} a & b \\ c & d \\ \end{array} \right]\Longrightarrow A^{-1}=\frac{1}{ad-bc}\left[ \begin{array}{rr} d & -b \\ -c & a \\ \end{array} \right]

Dada la matriz AA, calcular la matriz inversa a partir de la definición:

Si multiplicamos ambas matrices para obtener la matriz identidad 2×22 \times 2 tenemos dos sistemas de ecuaciones lineales que permiten esclarecer el porqué de la parte bb de la ecuación (M7)(M7):

2a+3b=13a+2b=0}2c+3d=03c+2d=1} \left. \begin{array}{rcr} 2a+3b & = & 1 \\ -3a+2b & = & 0 \end{array} \right\}\qquad \left. \begin{array}{rcr} 2c+3d & = & 0 \\ -3c+2d & = & 1 \end{array} \right\}
101

Solucionando ambos sistemas se obtiene que:

a=213b=313c=313d=213a=\frac{2}{13}\qquad b=\frac{3}{13} \qquad c=-\frac{3}{13} \qquad d=-\frac{2}{13} [2332][213313313213]=[1001] \left[ \begin{array}{rrr} 2 & 3 \\ -3 & 2\\ \end{array} \right]{\large \left[ \begin{array}{rrr} \frac{2}{13} & -\frac{3}{13} \\ \\ -\frac{3}{13} & \frac{2}{13} \\ \end{array} \right]}=\left[ \begin{array}{rrr} 1 & 0 \\ 0 & 1 \\ \end{array} \right]

En nuestro siguiente ejercicio vamos a comprobar que dada la matriz AA y multiplicada por A1A^{-1} se obtiene la matriz identidad.

Dadas las matrices AA y A1A^{-1} comprobar si su producto es la matriz identidad:

A=[234111201]A1=[1311311123253] A=\left[ \begin{array}{rrr} 2 & -3 & 4 \\ 1 & 1 & -1 \\ 2 & 0 & 1 \\ \end{array} \right]\qquad A^{-1}=\left[ \begin{array}{rrr} \frac{1}{3} & 1 & -\frac{1}{3} \\ \\ 1 & 1 & -1 \\ \\ -\frac{2}{3} & -2 & \frac{5}{3}\\ \end{array} \right] a11=(2)(13)+(3)(1)+(4)(23)=1\qquad{\color{red}{a_{11}}}=(2)\left(\frac{1}{3}\right)+(-3)(-1)+(4)\left(-\frac{2}{3}\right)=1

a12=(2)(1)+(3)(2)+(4)(2)=0\qquad{\color{red}{a_{12}}}=(2)(1)+(-3)(-2)+(4)(-2)=0

a13=(2)(13)+(3)(2)+(4)(53)=0\qquad{\color{red}{a_{13}}}=(2)\left(-\frac{1}{3}\right)+(-3)(2)+(4)\left(\frac{5}{3}\right)=0

a21=(1)(13)+(1)(1)+(1)(23)=0\qquad{\color{red}{a_{21}}}=(1)\left(\frac{1}{3}\right)+(1)(-1)+(-1)\left(-\frac{2}{3}\right)=0

102

a22=(1)(1)+(1)(2)+(1)(2)=1\qquad{\color{red}{a_{22}}}=(1)(1)+(1)(-2)+(-1)(-2)=1

a23=(1)(13)+(1)(2)+(1)(53)=0\qquad{\color{red}{a_{23}}}=(1)\left(-\frac{1}{3}\right)+(1)(2)+(-1)\left(\frac{5}{3}\right)=0

a31=(2)(13)+(0)(1)+(1)(23)=0\qquad{\color{red}{a_{31}}}=(2)\left(\frac{1}{3}\right)+(0)(-1)+(1)\left(-\frac{2}{3}\right)=0

a32=(2)(1)+(0)(2)+(1)(2)=0\qquad{\color{red}{a_{32}}}=(2)(1)+(0)(-2)+(1)(-2)=0

a33=(2)(13)+(0)(1)+(1)(53)=1\qquad{\color{red}{a_{33}}}=(2)\left(-\frac{1}{3}\right)+(0)(-1)+(1)\left(\frac{5}{3}\right)=1

Hemos comprobado que:

[234111201][1311311123253]=[100010001] \left[ \begin{array}{rrr} 2 & -3 & 4 \\ 1 & 1 & -1 \\ 2 & 0 & 1 \\ \end{array} \right]\left[ \begin{array}{rrr} \frac{1}{3} & 1 & -\frac{1}{3} \\ \\ 1 & 1 & -1 \\ \\ -\frac{2}{3} & -2 & \frac{5}{3}\\ \end{array} \right]=\left[ \begin{array}{rrr} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \\ \end{array} \right]

3.13 Matriz inversa por el método de Gauss

El cálculo de la matriz inversa por el método de Gauss supone transformar una matriz en otra equivalente por filas para lograr:

  1. Formar una matriz de orden n×2nn \times 2n tal que las primeras columnas sean las de la matriz AA y las otras nn las de la matriz identidad de orden nn.
  2. Mediante las transformaciones elementales de las filas de una matriz, convertir la matriz anterior en otra que tenga en las n primeras columnas la matriz identidad y en las nn últimas otra matriz que prescitamente será A1A^{-1}.
103

Una matriz aumentada, o matriz ampliada, se obtiene al combinar dos matrices.

Para hallar la matriz inversa combinaremos la matriz cuadrada dada con la respectiva matriz identidad, se representa como: (AI)(A|I) tal y como se muestra a continuación:

A|II|A1(M6)\tag{M6} A\textbf{|}I\rightarrow I\textbf{|}A^{-1}

Ejemplo de matriz 2×22 \times 2 ampliada con la matriz identidad 2×22 \times 2.

A|I=[32|1041|01] A\textbf{|}I=\left[ \begin{array}{rrr} 3 & 2 & \textbf{|} & \color{red}{1} & \color{red}{0}\\ 4 & 1 & \textbf{|} & \color{red}{0} & \color{red}{1}\\ \end{array} \right]

A continuación en el vídeo ilustraremos el procedimiento para hallar la matriz inversa de un sistema 2×22 \times 2:

104

Ahora veamos mediante este ejercicio las transformaciones necesarias para hallar la matriz inversa de una matriz 3×33 \times 3

A=[021202110] A=\left[ \begin{array}{rrr} 0 & 2 & -1 \\ 2 & 0 & -2 \\ -1 & 1 & 0 \\ \end{array} \right] A|I=[021|100202|010110|001] A\textbf{|}I=\left[ \begin{array}{ccccc} 0 & 2 & -1 & \textbf{|} & \color{red}{1} & \color{red}{0} & \color{red}{0}\\ 2 & 0 & -2 & \textbf{|} & \color{red}{0} & \color{red}{1} & \color{red}{0}\\ -1 & 1 & 0 & \textbf{|} & \color{red}{0} & \color{red}{0} & \color{red}{1}\\ \end{array} \right]

Observamos que todos los elementos de la diagonal principal son CERO y por ello debemos cambiar F1F_1 como punto de partida:

F1=F1+F2A|I=[223|110202|010110|001]F_1=F_1+F_2 \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} \color{red}{2} & \color{red}{2} & \color{red}{-3} & \textbf{|} & \color{red}{1} & \color{red}{1} & \color{red}{0}\\ 2 & 0 & -2 & \textbf{|} & 0 & 1 & 0\\ -1 & 1 & 0 & \textbf{|} & 0 & 0 & 1\\ \end{array} \right] F2=F2F1A|I=[223|110021|100110|001]F_2=F_2-F_1 \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} 2 & 2 & -3 & \textbf{|} & 1 & 1 & 0\\ \color{red}{0} & \color{red}{-2} & \color{red}{1} & \textbf{|} & \color{red}{-1} & \color{red}{0} & \color{red}{0}\\ -1 & 1 & 0 & \textbf{|} & 0 & 0 & 1\\ \end{array} \right] F3=2F3+F1A|I=[223|110021|100043|112]F_3=2F_3+F_1 \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} 2 & 2 & -3 & \textbf{|} & 1 & 1 & 0\\ 0 & -2 & 1 & \textbf{|} & -1 & 0 & 0\\ \color{red}{0} & \color{red}{4} & \color{red}{-3} & \textbf{|} & \color{red}{1} & \color{red}{1} & \color{red}{2}\\ \end{array} \right] F3=F3+2F2A|I=[223|110021|100001|112]F_3=F_3+2F_2 \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} 2 & 2 & -3 & \textbf{|} & 1 & 1 & 0\\ 0 & -2 & 1 & \textbf{|} & -1 & 0 & 0\\ \color{red}{0} & \color{red}{0} & \color{red}{-1} & \textbf{|} & \color{red}{-1} & \color{red}{1} & \color{red}{2}\\ \end{array} \right]
105
F1=F1+3F2A|I=[240|210021|100001|112]F_1=F_1+3F_2 \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} \color{red}{2} & \color{red}{-4} & \color{red}{0} & \textbf{|} & \color{red}{-2} & \color{red}{1} & \color{red}{0}\\ 0 & -2 & 1 & \textbf{|} & -1 & 0 & 0\\ 0 & 0 & -1 & \textbf{|} & -1 & 1 & 2\\ \end{array} \right] F2=F2+F3A|I=[240|210020|212001|112]F_2=F_2+F_3 \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} 2 & -4 & 0 & \textbf{|} & -2 & 1 & 0\\ \color{red}{0} & \color{red}{-2} & \color{red}{0} & \textbf{|} & \color{red}{-2} & \color{red}{1} & \color{red}{2}\\ 0 & 0 & 1 & \textbf{|} & -1 & 1 & 2\\ \end{array} \right] F1=F12F2A|I=[200|214020|212001|112]F_1=F_1-2F_2 \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} \color{red}{2} & \color{red}{0} & \color{red}{0} & \textbf{|} & \color{red}{2} & \color{red}{-1} & \color{red}{-4}\\ 0 & -2 & 0 & \textbf{|} & -2 & 1 & 2\\ 0 & 0 & -1 & \textbf{|} & -1 & 1 & 2\\ \end{array} \right] F1=F12A|I=[100|1122020|212001|112]F_1=\frac{F_1}{2} \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} \color{red}{1} & \color{red}{0} & \color{red}{0} & \textbf{|} & \color{red}{1} & \color{red}{-\frac{1}{2}} & \color{red}{-2}\\ 0 & -2 & 0 & \textbf{|} & -2 & 1 & 2\\ 0 & 0 & -1 & \textbf{|} & -1 & 1 & 2\\ \end{array} \right] F2=F22A|I=[100|1122010|0121001|112]F_2=\frac{F_2}{-2} \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} 1 & 0 & 0 & \textbf{|} & 1 & -\frac{1}{2} & -2\\ \color{red}{0} & \color{red}{1} & \color{red}{0} & \textbf{|} & \color{red}{0} & \color{red}{-\frac{1}{2}} & \color{red}{-1}\\ 0 & 0 & -1 & \textbf{|} & -1 & 1 & 2\\ \end{array} \right] F3=F31A|I=[100|1122010|0121001|112]F_3=\frac{F_3}{-1} \longrightarrow A\textbf{|}I=\left[ \begin{array}{ccccc} 1 & 0 & 0 & \textbf{|} & 1 & -\frac{1}{2} & -2\\ 0 & 1 & 0 & \textbf{|} & 0 & -\frac{1}{2} & -1\\ \color{red}{0} & \color{red}{0} & \color{red}{1} & \textbf{|} & \color{red}{-1} & \color{red}{-1} & \color{red}{2}\\ \end{array} \right]

Lo realizado en a través del ejercicio desarrollado implica que la matriz AA se convierte en la matriz identidad y la matriz BB (Identidad) se convierte en la matriz inversa de la matriz AA.

106
A1=[11220121112] A^{-1}=\left[ \begin{array}{ccc} 1 & -\frac{1}{2} & -2\\ \\ 0 & -\frac{1}{2} & -1\\ \\ -1 & -1 & 2 \end{array} \right]

3.14 Matriz inversa por cofactores

La matriz adjunta, de adjuntos o de cofactores de la matriz AA, que denotamos por Adj(A)Adj(A) o AA*, es la matriz cuyo elemento (i,j)(i , j) (fila ii y columna jj) es el adjunto adiad i, j=(1)i+jdet(Ai,j)j = ( - 1 ) i + j · det( A i , j ) donde la matriz Ai,jA i , j es la matriz que resulta al quitar a la matriz AA la fila ii y la columna jj.

Ello implica que para la matriz de adjuntos de un sistema 2×22 \times 2 debemos tener en cuenta los siguientes signos de sus componentes:

Adj(A)=A=[++]\bf {Adj(A)=A^{*}=\left[ \begin{array}{cc} + & -\\ - & +\\ \end{array} \right]}

Se cumple que la matriz inversa de AA, A1A^{-1} se puede escribir en función de su adjunta como

A1=(Adj(A))tA=(At)A(M7)\tag{M7} \bf {A^{-1}=\frac{(Adj(A))^t}{|A|}=\frac{(A^t)^*}{|A|}}

Donde la notación de la potencia tt expresa la traspuesta de AA.

Dada la matriz AA determinar su inversa por cofactores:

107

A=[5113]\bf {A=\left[ \begin{array}{cr} 5 & -1\\ 1 & 3\\ \end{array} \right]}
  1. Paso 1. Verificar si el determinante de la matriz A es =/{=}\llap{/\,} de CERO:
  2. A=[(5)(3)][(1)(1)]=16|A|=[(5)(3)]-[(-1)(1)]=16
  3. Paso 2. Determinar la matriz Adj(A)Adj (A) por cofactores:
  4. Adj11=(1)2det(3)=3\qquad Adj_{11}=(-1)^2det(3)=3

    Adj12=(1)3det(1)=1\qquad Adj_{12}=(-1)^3det(1)=-1

    Adj21=(1)3det(1)=1\qquad Adj_{21}=(-1)^3det(-1)=1

    Adj22=(1)4det(5)=5\qquad Adj_{22}=(-1)^4det(5)=5

    Adj(A)=A=[3115]\bf {Adj(A)=A^*=\left[ \begin{array}{cr} 3 & -1\\ 1 & 5\\ \end{array} \right]}
  5. Paso 3. Establecemos (A)t(A^*)^t
  6. (Adj(A))t=(A)t=[3115]\bf {(Adj(A))^t=(A^*)^t=\left[ \begin{array}{cr} 3 & 1\\ -1 & 5\\ \end{array} \right]}
  7. Paso 4. Determinar la matriz A1A^{-1}:
  8. A1=116[3115]=[316116116516]A^{-1}=\frac{1}{16}\left[ \begin{array}{cr} 3 & 1\\ -1 & 5\\ \end{array} \right]=\Large {\left[ \begin{array}{cr} \frac{3}{16} & \frac{1}{16}\\ -\frac{1}{16} & \frac{5}{16}\\ \end{array} \right]}
108

Al igual que la adjunta de un sistema 2×22 \times 2, la adjunta de un sistema 3×33 \times 3 implica que debemos tener en cuenta los signos de sus componentes como se presenta a continuación:

Adj(A)=[a11a12a13a21a22a23a31a32a33]=[+++++]\bf {Adj(A)=\left[ \begin{array}{ccc} a_{11} & a_{12} & a_{13}\\ a_{21} & a_{22} & a_{23}\\ a_{31} & a_{32} & a_{33} \end{array} \right]=\left[ \begin{array}{ccc} + & - & +\\ - & + & -\\ + & - & + \end{array} \right]}

La matriz de cofactores tendría la siguiente estructura donde se resalta el signo correspondiente a cada elemento:

Adjunta de una matriz 3×33 \times 3

Dada la matriz AA determinar su inversa por cofactores:

A=[222110034]\bf {A=\left[ \begin{array}{ccc} 2 & 2 & 2\\ 1 & -1 & 0\\ 0 & 3 & 4 \end{array} \right]}
  1. Paso 1. Hallar el determinante aplicando la regla de Sarrus.
109

A=222110034221103|A|=\begin{vmatrix} 2 & 2 & 2\\ 1 & -1 & 0\\ 0 & 3 & 4 \end{vmatrix}\begin{matrix} 2 & 2 \\ 1 & -1 \\ 0 & 3 \end{matrix} A=[(2)(1)(4)+(2)(0)(0)+(2)(1)(3)]|A|=[(2)(-1)(4)+(2)(0)(0)+(2)(1)(3)] [(0)(1)(2)+(3)(0)(2)+(4)(1)(2)]=10-[(0)(-1)(2)+(3)(0)(2)+(4)(1)(2)]=-10
  1. Paso 2. Dado que el determinante es de CERO, procedemos con la determinación de la matriz de cofactores:
  2. a11=(1)1+1[(1)(4)(0)(3)]=4\qquad a_{11}=(-1)^{1+1}[(-1)(4)-(0)(3)]=-4

    a12=(1)1+2[(1)(4)(0)(0)]=4\qquad a_{12}=(-1)^{1+2}[(1)(4)-(0)(0)]=-4

    a13=(1)1+3[(1)(3)(1)(0)]=3\qquad a_{13}=(-1)^{1+3}[(1)(3)-(1)(0)]=3

    a21=(1)2+1[(2)(4)(3)(2)]=2\qquad a_{21}=(-1)^{2+1}[(2)(4)-(3)(2)]=-2

    a22=(1)2+2[(2)(4)(0)(2)]=8\qquad a_{22}=(-1)^{2+2}[(2)(4)-(0)(2)]=8

    a23=(1)2+3[(2)(3)(2)(0)]=6\qquad a_{23}=(-1)^{2+3}[(2)(3)-(2)(0)]=-6

    a31=(1)3+1[(2)(0)(1)(2)]=2\qquad a_{31}=(-1)^{3+1}[(2)(0)-(-1)(2)]=2

    a32=(1)3+2[(2)(0)(1)(2)]=2\qquad a_{32}=(-1)^{3+2}[(2)(0)-(1)(2)]=2

    a33=(1)3+3[(2)(1)(2)(1)]=4\qquad a_{33}=(-1)^{3+3}[(2)(-1)-(2)(1)]=-4

  3. Paso 3. Ahora hallamos (A)t(A*)^t = (Adj(A))t(Adj (A))^t:
110

A=[443286224];(Adj(A))t=[422482364]\bf {A=\left[ \begin{array}{ccc} -4 & -4 & 3\\ -2 & 8 & -6\\ 2 & 2 & -4 \end{array} \right];(Adj(A))^t=\left[ \begin{array}{ccc} -4 & -2 & 2\\ -4 & 8 & 2\\ 3 & -6 & -4 \end{array} \right]}
  1. Paso 4. Hallamos la matriz inversa:
A1=110[422482364]=\bf {A^{-1}=\frac{1}{-10}\left[ \begin{array}{ccc} -4 & -2 & 2\\ -4 & 8 & 2\\ 3 & -6 & -4 \end{array} \right]=} A1=[2515152545153103525]\bf {A^{-1}=\left[ \begin{array}{ccc} \large{\frac{2}{5}} & \large{\frac{1}{5}} & \large{-\frac{1}{5}}\\ \\ \large{\frac{2}{5}} & \large{-\frac{4}{5}} & \large{-\frac{1}{5}}\\ \\ \large{-\frac{3}{10}} & \large {\frac{3}{5}} & \large{-\frac{2}{5}} \end{array} \right]}

Al principio del libro se realizó un recuento de cómo resolver muchos problemas matemáticos a través del planteamiento de un sistema de ecuaciones lineales. Existen diversos métodos para resolverlos (recodemos los conocidos método de igualación, substitución e igualación).

Pero, si dichos sistemas están formados por gran cantidad de ecuaciones e incógnitas, el aplicar los métodos anteriores resulta ser una tarea sumamente dificultosa y que, muy probablemente, conducirá a resultados erróneos. Para solventar este problema, los sistemas de ecuaciones lineales se pueden plantear matricialmente y resolverlos haciendo uso de la matriz inversa.

En el próximo capitulo solucionaremos algunas aplicaciones de la vida real empleando la matriz inversa.

111
112



Aplicaciones
Sistemas de ecuaciones
Vectores
Criptografía
116

4.1 Solución de ecuaciones

Las matrices constituyen una herramienta que nos permite saber rápidamente si un sistema de ecuaciones lineales (SEL) tiene soluciones y el tipo de soluciones: un sistema puede no tener solución, tener una única solución o tener infinitas soluciones.

Para estas aplicaciones se requiere comprender los conceptos relacionados con las matrices: matrices regular, inversa y adjunta y el rango de una matriz.

Una vez clasificado el SEL se procede a calcular, en caso de haberlas, sus soluciones. Para ello disponemos de varios métodos: El método de Gauss, la Regla de Cramer y el de la matriz inversa si el SEL es compatible determinado, y la eliminación de Gauss o de Gauss-Jordan si es un sistema indeterminado. Para esa solución trabajamos el SEL como un sistema matricial (matriz de coeficientes y matriz o vector de términos independientes).

A partir de su matriz ampliada del SEL (matriz de coeficientes con el vector de términos independientes), el teorema de Rouché-Frobenius, nos permite clasificar un SEL según el tipo de solución.

En esta sección emplearemos las matrices para solucionar sistemas de ecuaciones de nn ecuaciones lineales y nn variables.

Para solucionar ecuaciones lineales vamos a recurrir a:

  • Verificación del determinante de los coeficientes el cual debe ser diferente de CERO.
  • Multiplicación de matrices.
  • Igualdad de matrices.
117

Un sistema lineal lo podemos escribir como una ecuación de matrices:

AX=BA \cdot X = B

Donde AA es la matriz de coeficientes del sistema, XX es la matriz de las variables (para este ecuación la variable es xx) y BB es la matriz de los términos independientes.

A=[a11a12a1na21a22a2nam1am2amn];X=[x1x2xn];B=[b1b2bn]A=\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22} & \cdots & a_{2n}\\ \vdots & \vdots & \vdots & \vdots\\ a_{m1} & a_{m2} & \cdots & a_{mn}\\ \end{bmatrix};X=\begin{bmatrix} x_1\\ x_2\\ \vdots\\ x_n\\ \end{bmatrix};B=\begin{bmatrix} b_1\\ b_2\\ \vdots\\ b_n\\ \end{bmatrix}

A)

Sistema lineal 2×22 \times 2, ilustraremos el proceso de solución con un ejemplo paso a paso y una aplicación:

Resolver el sistema:

E13x12y=27E25x+3y=11\begin{array}{rrrr} E_1&&&3x&-&12y&=-&27 \\ E_2&&-&5x&+&3y&=&11 \end{array}

Paso 1. Escribir la ecuación de matrices:

[31253][xy]=[2711] \begin{bmatrix} 3 & -12\\ -5 & 3\\ \end{bmatrix} \begin{bmatrix} x\\ y\\ \end{bmatrix}= \begin{bmatrix} -27\\ 11\\ \end{bmatrix}

Paso 2. Verificar si el determinante de la matriz de coeficientes es diferente de CERO.

det(A)=A=(3)(3)(5)(12)=51det(A)=|A|=(3)(3)-(-5)(-12)=-51
118

Paso 3. Determinar la matriz inversa de la matriz de coeficientes, se emplea uno de los métodos descritos en el capitulo anterior y comprobamos que se obtiene la matriz identidad:

[31253][35123]=151[510051]=[1001] \begin{bmatrix} 3 & -12\\ -5 & 3\\ \end{bmatrix} \begin{bmatrix} 3 & 5\\ 12 & 3\\ \end{bmatrix}=-\frac{1}{51} \begin{bmatrix} -51 & 0\\ 0 & -51\\ \end{bmatrix}= \begin{bmatrix} 1 & 0\\ 0 & 1\\ \end{bmatrix}

Ello implica que:

[xy]=151[3511251551351][2711]=[12] \begin{bmatrix} x\\ y\\ \end{bmatrix}= -\frac{1}{51}\begin{bmatrix} -\frac{3}{51} & -\frac{12}{51}\\ &\\ -\frac{5}{51} & -\frac{3}{51}\\ \end{bmatrix}\begin{bmatrix} -27\\ 11\\ \end{bmatrix}=\begin{bmatrix} -1\\ 2\\ \end{bmatrix}

Ahora comprobemos en E1E_1y en E2E_2:

3(1)12(2)=?324=273(-1)-12(2)=?\Longrightarrow -3-24=-27 5(1)+3(2)=?5+6=11-5(-1)+3(2)=?\Longrightarrow 5+6=11

Aplicación

En un remate de aduanas ofrecen un lote compuesto de: 400400 celulares tipo A1A_1 y 1500 celulares del tipo Z3Z_3 por 110500110500 dólares. Otro lote ofrecido consta de 15001500 celulares del tipo A1A_1 y 600600 celulares del tipo Z3Z_3 por 138000138000 dólares. ¿Cuál es el precio de cada tipo de celular?

Establecemos el sistema de ecuaciones y emplearemos como variables la misma denominación de los tipos de celular:

E1400A1+1500Z3=110500E21500A1+600Z3=138000\begin{array}{rrrr} E_1&&&400A_1&+&1500Z_3&=110500 \\ E_2&&&1500A_1&+&600Z_3&=138000 \end{array}
119

Paso 1. Establecemos la ecuación de matrices:

[40015001500600][A1Z3]=[110500138000]\begin{bmatrix} 400 & 1500\\ 1500 & 600\\ \end{bmatrix} \begin{bmatrix} A_1\\ Z_3\\ \end{bmatrix}=\begin{bmatrix} 110500\\ 138000\\ \end{bmatrix}

Paso 2. Verificar si el determinante de la matriz de coeficientes es diferente de CERO.

A=(400)(600)(1500)(1500)A=2400002250000|A|=(400)(600)-(1500)(1500)\Longrightarrow |A|=240000-2250000 A=2010000\Longrightarrow |A|=-2010000

Paso 3. Determinar la matriz inversa de la matriz de coeficientes, se emplea uno de los métodos descritos en el capitulo anterior y comprobamos que se obtiene la matriz identidad:

[40015001500600][60015001500400] \begin{bmatrix} 400 &1500\\ 1500 &600\\ \end{bmatrix} \begin{bmatrix} 600 & -1500\\ -1500 & 400\\ \end{bmatrix}\Longrightarrow =12010000[2010000002010000]=[1001]=-\frac{1}{2010000} \begin{bmatrix} -2010000 & 0\\ 0 & -2010000\\ \end{bmatrix}=\begin{bmatrix} 1 & 0\\ 0 & 1\\ \end{bmatrix}

Este arreglo matricial nos lleva a:

[A1Z3]=12010000[60015001500400][110500138000] \begin{bmatrix} A_1\\ Z_3\\ \end{bmatrix}= -\frac{1}{2010000}\begin{bmatrix} 600 & -1500\\ &\\ -1500 & 400\\ \end{bmatrix}\begin{bmatrix} 110500\\ 138000\\ \end{bmatrix} [A1Z3]=1201000[140700000110550000]=[7055] \begin{bmatrix} A_1\\ Z_3\\ \end{bmatrix}=-\frac{1}{201000} \begin{bmatrix} -140700000\\ -110550000\\ \end{bmatrix}=\begin{bmatrix} 70\\ 55\\ \end{bmatrix}

Le dejamos al lector la comprobación de las ecuaciones mediante la sustitución de A1A_1 con el valor de 7070 dólares y de Z3Z_3 con un costo de 5555 dólares.

120

B)

Sistema lineal 3×33 \times 3, al igual que en el apartado AA, ilustraremos el proceso de solución con un ejemplo paso a paso y una aplicación:

En el siguiente ejercicio hallaremos la solución de un sistema 3×33 \times 3

Hallar los valores de x1x_1, x2x_2 y de x3x_3 que satisfacen el siguiente sistema:

4x1+5x28x3=97x112x2+5x3=23x1+9x2+10x3=61\begin{alignedat}{3} 4&x_1+&5&x_2-8&x_3 = -9 \\ 7&x_1-&12&x_2+5&x_3 = -2\\ -3&x_1+&9&x_2+10&x_3=61 \end{alignedat}

Paso 1. Establecemos la ecuación de matrices:

[45871253910][x1x2x3]=[9261] \begin{bmatrix} 4 & 5 & -8\\ 7 & -12 & 5\\ -3 & 9 & 10\\ \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3\\ \end{bmatrix}=\begin{bmatrix} -9\\ -2\\ 61\\ \end{bmatrix}

Paso 2. Verificar si el determinante de la matriz de coeficientes es diferente de CERO mediante los elementos de la primera línea:

A=4[(12)(10)(9)(5)]|A|=4[(-12)-(10)-(9)(5)] 5[(7)(10)(3)(5)]-5[(7)(10)-(-3)(5)] 8[(7)(9)(3)(12)]=1301-8[(7)(9)-(-3)(-12)]=-1301

Paso 3. Hallar A1A^{-1} ya que A0|A| ≠ 0, iniciamos con la matriz de adjuntos:

121

a11=(1)1+1[(12)(10)(9)(5)]=165a_{11}=(-1)^{1+1}[(-12)(10)-(9)(5)]=-165 a12=(1)1+2[(7)(10)(3)(5)]=85a_{12}=(-1)^{1+2}[(7)(10)-(-3)(5)]=-85 a13=(1)1+3[(7)(9)(3)(12)]=27a_{13}=(-1)^{1+3}[(7)(9)-(-3)(-12)]=27 a21=(1)2+1[(5)(10)(9)(8)]=122a_{21}=(-1)^{2+1}[(5)(10)-(9)(-8)]=-122 a22=(1)2+2[(4)(10)(3)(8)]=16a_{22}=(-1)^{2+2}[(4)(10)-(-3)(-8)]=16 a23=(1)2+3[(4)(9)(3)(5)]=51a_{23}=(-1)^{2+3}[(4)(9)-(-3)(5)]=-51 a31=(1)3+1[(5)(5)(12)(8)]=71a_{31}=(-1)^{3+1}[(5)(5)-(-12)(-8)]=-71 a32=(1)3+2[(4)(5)(7)(8)]=76a_{32}=(-1)^{3+2}[(4)(5)-(7)(-8)]=-76 a33=(1)3+3[(4)(12)(7)(5)]=83a_{33}=(-1)^{3+3}[(4)(-12)-(7)(5)]=-83 Adj(A)=A=[16585271221651717683]Adj(A)=A^{*}=\begin{bmatrix} -165 & -85 & 27\\ -122 & 16 & -51\\ -71 & -76 & -83\\ \end{bmatrix}

Paso 4. Ahora su traspuesta:

(Adj(A))T=(A)T=[16512271851676275183](Adj(A))^{T}=(A^{*})^{T}=\begin{bmatrix} -165 & -122 & -71\\ -85 & 16 & -76\\ 27 & -51 & -83\\ \end{bmatrix}

Paso 5. Determinemos A1A^{-1} mediante:

A1=1A(Adj(A))TA^{-1}=\frac{1}{|A|}(Adj(A))^{T}
122

A1=11301[16512271851676275183]A^{-1}=\frac{1}{-1301}\begin{bmatrix} -165 & -122 & -71\\ -85 & 16 & -76\\ 27 & -51 & -83\\ \end{bmatrix}

A1=[16513011221301711301851301161301761301271301511301831301]A^{-1}=\begin{bmatrix} \frac{165}{1301} & \frac{122}{1301} & \frac{71}{1301}\\ \frac{85}{1301} & \frac{16}{1301} & \frac{76}{1301}\\ -\frac{27}{1301} & \frac{51}{1301} & \frac{83}{1301}\\ \end{bmatrix}

Paso 6. Solucionemos la ecuación matricial:

[x1x2x3]=[16513011221301711301851301161301761301271301511301831301][9261]\begin{bmatrix} x_1\\ x_2\\ x_3\\ \end{bmatrix}=\begin{bmatrix} \frac{165}{1301} & \frac{122}{1301} & \frac{71}{1301}\\ \frac{85}{1301} & \frac{16}{1301} & \frac{76}{1301}\\ -\frac{27}{1301} & \frac{51}{1301} & \frac{83}{1301}\\ \end{bmatrix}\begin{bmatrix} -9\\ -2\\ 61\\ \end{bmatrix} [x1x2x3]=[1651301(9)+1221301(2)+711301(61)851301(9)+161301(2)+761301(61)271301(9)+511301(2)+831301(61)]=[234]\begin{bmatrix} x_1\\ x_2\\ x_3\\ \end{bmatrix}=\begin{bmatrix} \frac{165}{1301}(-9) + \frac{122}{1301}(-2) + \frac{71}{1301}(61)\\ \frac{85}{1301}(-9) + \frac{16}{1301}(-2) + \frac{76}{1301}(61)\\ -\frac{27}{1301}(-9) + \frac{51}{1301}(-2) + \frac{83}{1301}(61)\\ \end{bmatrix}=\begin{bmatrix} 2\\ 3\\ 4\\ \end{bmatrix}

El sistema tiene como solución: x1=2x_1 = 2 ; x2=3x_2 = 3 ; x3=4x_3 = 4, es consistente determinado.

Aplicación

Una ecuación cuadrática tiene como gráfica una parábola representada por el polinomio:

p(x)=a0±a1x±a2x2p(x)=a_0\pm a_1x\pm a_2x^{2}
123

Vamos a determinar los coeficientes del polinomio que representa la parábola que pasa por los puntos dados en la figura:

Cada par ordenado debe satisfacer la ecuación del polinomio dado, por lo tanto se forma el sistema:

a0+2a1+4a2=1a0+5a1+25a2=24a0+6a1+9a2=1\begin{alignedat}{3} &a_0+&2&a_1+4&a_2 = 1 \\ &a_0+&5&a_1+25&a_2 = 2\\ 4&a_0+&6&a_1+9&a_2=1 \end{alignedat}

Vídeo solución del problema planteado:

124

Independiente de su valor posicional en el sistema decimal, la suma de las cifras un número menor de 10001000 y mayor de 100100 es igual a 1717. Si al número original se le restan 297297 unidades el número se invierte posicionalmente. Además, en el número original la suma de las cifras las unidades y la de las centenas excede en 55 unidades a la cifra de las decenas. ¿Cuál es ese número?

Establecemos las variables

A=A102B=B101C=C100A = A*10^{2} \qquad\qquad B = B*10^{1} \qquad \qquad C = C* 10^{0}

Acorde con ello, la traducción del lenguaje Vernáculo genera las siguientes ecuaciones:

125

Ecuación 1:

a+b+c=17a+b+c=17

Ecuación 2:

100A+10B+C297=100C+10B+A100A+10B+C-297=100C+10B+A

Ecuación 3:

A+C=B+5A+C=B+5

El sistema generado luego de simplificar es entonces:

A+B+C=17A+B+C=17 AC=3A-C=3 AB+C=5A-B+C=5

Ahora la matriz ampliada AB\color{red}A|B:

[1111710131115] \left[ \begin{array}{lrrcr} 1 & 1 & 1 & \color{blue}\vert & \color{red}17\\ 1 & 0 & -1 & \color{blue}\vert & \color{red}3\\ 1 & -1 & 1 & \color{blue}\vert & \color{red}5\\ \end{array} \right]

Solución por el método de Gauss escalonando la matriz ampliada:

F2=F2F3F_2=F_2-F_3 \Longrightarrow [1111701221115] \left[ \begin{array}{lrrcr} 1 & 1 & 1 & \color{blue}\vert & \color{red}17\\ 0 & 1 & -2 & \color{blue}\vert & \color{red}-2\\ 1 & -1 & 1 & \color{blue}\vert & \color{red}5\\ \end{array} \right]
126
F3=F3F1F_3=F_3-F_1 \Longrightarrow [11117012202012] \left[ \begin{array}{lrrcr} 1 & 1 & 1 & \color{blue}\vert & \color{red}17\\ 0 & 1 & -2 & \color{blue}\vert & \color{red}-2\\ 0 & -2 & 0 & \color{blue}\vert & \color{red}-12\\ \end{array} \right] F3=F3+2F2F_3=F_3+2F_2 \Longrightarrow [11117012200416] \left[ \begin{array}{lrrcr} 1 & 1 & 1 & \color{blue}\vert & \color{red}17\\ 0 & 1 & -2 & \color{blue}\vert & \color{red}-2\\ 0 & 0 & -4 & \color{blue}\vert & \color{red}-16\\ \end{array} \right]

El nuevo sistema reducido:

A+B+C=17B2C=24C=16}\left. \begin{array}{rcr} A+B+C & = & 17 \\ B-2C & = & -2 \\ -4C & = & -16 \end{array} \right\}

De este sistema reducido solucionamos de atrás hacia adelante:

4C=164C(14)=16(14)C=4-4C=-16 \Longrightarrow -4C \left(\frac{1}{-4}\right)=-16\left(\frac{1}{-4}\right)\Longrightarrow C=\color{red}4

Conocido el valor de CC tenemos:

B2C=2B2(4)=2B8=2B-2C=-2 \Longrightarrow B-2(4)=-2 \Longrightarrow B-8=-2 B8+8=2+8B=6B - 8 + 8 = -2 + 8 \Longrightarrow B = \color{red}6

Conocido los valores de BB y de CC tenemos:

A+6+4=17A+10=17A+1010=1710A + 6 + 4 = 17 \Longrightarrow A + 10 = 17 \Longrightarrow A + 10 - 10 = 17 - 10 A=7A = \color{red}7
127

El número buscado es:

A(102)+b¿B(101)+C(100)A(10^2)+b¿B(10^1)+C(10^0)\Longrightarrow 7(100)+6(10)+4(1)=700+60+4=7647(100)+6(10)+4(1)=700+60+4=764

Otro de los campos dónde existen los sistemas de ecuaciones es la solución de circuitos eléctricos, en el siguiente vídeo veremos una solución por el método de mallas.

4.2 Vectores y Matrices

1. Los espacios vectoriales

En el capítulo 2 se aplicaron ciertas propiedades a los conjuntos de

128

vectores en el plano (R2ℝ^2) y los vectores en el espacio (R3ℝ^3) de tal manera que, si sumamos dos o más vectores en el plano su resultante es otro vector en el mismo espacio, de igual manera sucede lo mismo con vectores en el espacio, expresado de otra forma se aplican las propiedades conmutativa y asociativa para la suma y la propiedad distributiva para la multiplicación por un escalar.

Definición: un espacio vectorial real V es un conjunto de objetos en R2ℝ^2 y R3ℝ^3 con dos operaciones que obedecen las reglas de suma de vectores y multiplicación por un escalar.

Es de resaltar que esos objetos se denominan vectores pero no siempre serán objetos que tienen dirección, magnitud y sentido como las fuerzas, velocidades y aceleraciones tal como veremos en algunos ejercicios resueltos.

Ejemplos de espacios vectoriales:

  • Vectores en el plano. Dimensión 2
129
  • Vectores en el espacio. Dimensión 3
  • Matrices2×2_{2 \times 2}. Dimensión 4
  • A=[a11a12a21a22]\Large A=\begin{bmatrix} a_{11} & a_{12}\\ \\ a_{21} & a_{22}\\ \end{bmatrix}
  • Matrices3×3_{3 \times 3}. Dimensión 9
  • A=[a11a12a13a21a22a23a31a32a33]\Large A=\begin{bmatrix} a_{11} & a_{12} & a_{13}\\ \\ a_{21} & a_{22} & a_{23}\\ \\ a_{31} & a_{32} & a_{33}\\ \end{bmatrix}
130
  • Polinomio de 1er grado. Dimensión 2
  • P1(x)=x7\large P_1(x) = x-7
  • Polinomio de 2o grado. Dimensión 3
  • P2(x)=x2+x7\large P_2(x) = x^2 + x - 7
  • Polinomio de 3er grado. Dimensión 4
  • P3(x)=7x3+x2+x7\large P_3(x) = 7x^3 + x^2 + x - 7

Para que un conjunto dado sea un espacio vectorial VV es necesario que se cumplan los siguientes axiomas:

  1. Si tanto xx como yVy ∈ V, (x+y)V∴ (x + y) ∈ V
  2. Para todos los vectores z,y,zV,(x+y)+z=x+(y+z)z, y, z ∈ V, (x + y) + z = x + (y + z)
  3. un vector 0V0 ∈ V tal que para xV,x+0=0+x=x∀ x ∈ V, x + 0 = 0 + x = x
  4. Si xVx ∈ V, un vector (x)V(-x) ∈ V, tal que x+(x)=0x + (-x) = 0
  5. Si tanto xx como yVy ∈ V, x+y=y+x;xy=yxx + y = y + x; xy = yx
  6. Si xVx ∈ V y αα es un escalar, αxV∴ αx ∈ V
  7. Si tanto xx como yVy ∈ V y αα es un escalar, α(x+y)=αx+αy∴ α (x + y) = αx + αy
  8. Si xVx ∈ V, αα y ββ son escalares, (α+β)x=αx+βx∴ (α + β)x = αx + βx
  9. Si xVx ∈ V, αα y ββ son escalares, α(βx)=(αβ)x∴ α (βx) = (αβ)x
  10. vector xVx ∈ V, 1x=x1x = x

Dado un conjunto VV, con solo demostrar que no cumple uno solo de los anteriores axiomas es razón suficiente para declarar que no es un espacio vectorial.

131

Espacio vectorial trivial. Se denomina así al conjunto V=0V = { 0 }, puesto que cumple todos los axiomas descritos ya que el CERO es su único elemento.

Verificar y la recta que pasa por los puntos A(1,0A(1,0) y B(4,0)B(4,0) constituye un espacio vectorial.

De la Geometría Analítica tenemos que su ecuación es:

y1x0=0140y=x4+1\frac{y-1}{x-0}=\frac{0-1}{4-0}\Longrightarrow y=-\frac{x}{4}+1

Ahora construyamos una tabla de valores que nos permitirá realizar la comprobación de los axiomas de los espacios vectoriales:

Dado que V=(x,y)y=0.25x+1V = {(x,y)|y = -0.25x + 1}, ahora vamos a comprobar que:

(x1,y1)+(x2,y2)=(x1+x2,y1+y2)(x_1,y_1)+(x_2,y_2)=(x_1+x_2,y_1+y_2) (1,0)+(4,0)=(1+4,0+0)=(5,0)(1,0)+(4,0)=(1+4,0+0)=(5,0)

Podemos observar que el par ordenado (5,0)(5,0) no se encuentra en los pares ordenados de la tabla, es decir no cumple la propiedad de cerradura para la suma. Por lo tanto:

V=(x,y)y=0.25x+1V = {(x,y)| y = -0.25x + 1}

NO es un espacio vectorial

132

Verificar si las matrices invertibles AA y BB forman o no forman un espacio vectorial.

A=[123011102] A=\left[ \begin{array}{rrr} 1 & 2 & 3 \\ 0 & 1 & -1 \\ 1 & 0 & -2 \\ \end{array} \right] B=[210101123]\color{#0000FF} B=\left[ \begin{array}{rrr} 2 & 1 & 0 \\ 1 & 0 & 1 \\ 1 & 2 & 3 \\ \end{array} \right]

Recordemos que una matriz es invertible si su determinante es diferente de CERO.

A=123011102123011A=[4][3]=7|A|=\begin{vmatrix} 1 & 2 & 3 \\ 0 & 1 & -1\\ 1 & 0 & -2\\ 1 & 2 & 3 \\ 0 & 1 & -1 \end{vmatrix}\Longrightarrow|A|=[-4]-[3]=-7 B=210101123210101B=[1][7]=6\color{#0000FF}|B|=\begin{vmatrix} 2 & 1 & 0 \\ 1 & 0 & 1\\ 1 & 2 & 3\\ 2 & 1 & 0 \\ 1 & 0 & 1 \end{vmatrix}\Longrightarrow|B|=[1]-[7]=-6
  1. El primer axioma se cumple ya que la suma de las matrices A3×3A_{3 \times 3} y B3×3B_{3 \times 3} dan como resultado una matriz 3×33 \times 3.
133
  1. En caso de sumar una matriz C3×3C_{3 \times 3} a las matrices AA y BB podemos aplicar la propiedad asociativa y el resultado será otra matriz 3×33 \times 3.
  2. Si a la matriz A3×3A_{3 \times 3} le sumamos su matriz opuesta A3×3-A_{3 \times 3} el resultado será una matriz nula 3×33 \times 3, de igual manera sucede con la matriz BB y su opuesta.
  3. Si a la suma de las matrices A3×3A_{3 \times 3} y B3×3B_{3 \times 3} le sumamos una matriz nula 3×33 \times 3, el resultado de ello será una matriz 3×33 \times 3.
  4. Dado que la multiplicación de matrices no es conmutativa, el resultado es diferente lo que implica que el axioma de conmutatividad no se cumple, veamos:
A×B=7711022036A \times B= \begin{vmatrix} 7 & 7 & 11 \\ 0 & -2 & -2\\ 0 & -3 & -6 \end{vmatrix}

Es diferente de:

B×A=255221445B \times A= \begin{vmatrix} 2 & 5 & 5 \\ 2 & 2 & 1\\ 4 & 4 & -5 \end{vmatrix}

Podemos afirmar que si V=MmnV = M_{mn} denota el conjunto de matrices de orden m×nm \times n con componentes reales, entonces con la suma de matrices y multiplicación por un escalar usuales se puede verificar que MmnM_{mn} es un espacio vectorial cuyo neutro aditivo es la matriz de ceros de dimensiones m×nm \times n es un espacio vectorial pero no la multiplicación de matrices puesto que no es conmutativa.

Podemos concluir que las matrices dadas no forman un espacio vectorial.

134

Ejercicio resuelto. Vídeo

Sub-espacios vectoriales

Se dice que WW es un subespacio vectorial de VV si WW es un subconjunto no vacío de VV, y WW es un espacio vectorial, junto con las operaciones de suma \oplus entre vectores y multiplicación por un escalar \otimes definidas para VV.

Regla de los sub-espacios vectoriales

Un subconjunto no vacío WW de un espacio vectorial VV es un subespacio de VV si se cumplen las dos reglas de cerradura.

135

Reglas que permitan determinar si un subconjunto NO vacío es un sub-espacio vectorial, esas dos reglas las podemos discriminar como:

  1. Condición 1. 0W0 ∈ W
  2. Condición 2. Si xx y yWy ∈ W, (x+y)W∴ (x + y) ∈ W
  3. Condición 3. Si xWx ∈ W, α(x)Wα∴ α(x) ∈ W ∀ α

Determinar si W={(x1,x2,x3)R3    5x13x3=0}W = \{(x_1, x_2, x_3) \in ℝ^3 \; | \; 5x_1 − 3x_3 = 0 \} es un sub-espacio vectorial de R3ℝ^3.

  • Condición 1: ¿el vector (0,0,0)W(0,0,0) ∈ W?
  • (5)(0)(3)(0)=00=0(5)(0) − (3)(0) = 0 − 0 = 0
  • Condición 2: Dados W1=(x1,y1,z1)W_1 = (x_1, y_1, z_1); W2=(x2,y2,z2)W_2 = (x_2, y_2, z_2) tal que W∈ W.
W1+W2=(x1+x2,y1+y2,z1+z2)W_1+W_2 = (x_1 + x_2, y_1 + y_2, z_1 + z_2)

Este vector lo sustituimos en la ecuación implícita dada:

5(x1+x2)3(z1+z2)=05(x_1 + x_2) − 3(z_1 + z_2) = 0

Aplicamos la propiedad distributiva:

5x1+5x23z13z2=05 x_1 + 5 x_2 − 3 z_1 − 3 z_2 = 0

Como W1W_1, W2WW_2 ∈ W, 5x13z1=0∴ 5 x_1− 3 z_1 = 0 (1); 5x23z2>=05 x_2 − 3 z_2> = 0 (2) y reescribiendo la ecuación de la línea anterior tenemos:

5x13z2+5x23z2=0;0+0=0\color{red}{5 x_1 − 3 z_2} + \color{blue}{5 x_2 − 3 z_2} = 0 ; 0 + 0 = 0
136
  • Condición 3.
α(x1,y1,z1)=(αx1,αy1,αz1)\alpha (x_1, y_1, z_1) = (\alpha x_1, \alpha y_1, \alpha z_1)

Reemplazando en la ecuación implícita:

5αx13αx3=05 \alpha x_1 − 3 \alpha x_3 = 0

Obteniendo un factor común:

α(5x13x3)=0;0=0\alpha (5x_1 − 3 x_3) =0;0=0

Dado que hemos demostrado que se cumplen las tres condiciones, podemos afirmar que WW es un sub-espacio vectorial de VV.

Operaciones con Sub-espacios vectoriales

Si SS y WW son sub-espacios de VV que cumplen las operaciones de y se realizan las siguientes operaciones entre ellos:

  1. Unión
SW=uVuSuWS ∪ W = {u ∈ V | u ∈ S ∨ u ∈ W}

En general la unión de sub-espacios NO es un sub-espacio de VV.

  1. Intersección

La intersección de sub-espacios es un sub-espacio de VV.

Si SW=0S ∩ W = {0}

SW=uVuSuWS ∩ W = {u ∈ V | u ∈ S ∧ u ∈ W}
137

A esa intersección se le denomina suma directa y se expresa como SWS ⊕ W. A este se le denomina sub-espacio trivial formado por el vector nulo.

  1. Suma
S+W=uVu=(u1+u1)u1Su2WS + W = {u ∈ V | u = (u1 + u1) ∧ u1 ∈ S ∧ u2 ∈ W}

Vídeo

Ahora podemos hablar de la combinación lineal en esos subespacios vectoriales.

138

Combinación lineal de vectores

Una combinación lineal de vectores v1v_1, v2v_2, v3v_3, ......, vnv_n es en un espacio vectorial VV, es la suma que tiene la forma:

α1v1+α2>v2+α3v3+...+αnvn\alpha_1v_1 + \alpha_2>v_2 + \alpha_3v_3 + ... + \alpha_nv_n

Donde α1,α2,α3,...,αnα_1, α_2, α_3, ..., α_n son escalares.

En un espacio vectorial VV, los vectores v1,v2,v3,...,vnv_1, v_2, v_3, ..., v_n generan a VV si todo vector en VV se puede expresar como una combinación lineal de v1,v2,v3,...,vnv_1, v_2, v_3, ..., v_n

Expresado de otra manera, dados dos vectores: v1v_1 y v2v_2, y dos números: aa y bb, el vector resultante de av1+bv2av_1 + bv_2 se dice que es una combinación lineal de v1v_1 y v2v_2.

Una combinación lineal de dos o más vectores es el vector que se obtiene al sumar esos vectores multiplicados por sendos escalares. Cualquier vector se puede poner como combinación lineal de otros dos que tengan diferente dirección. Esta combinación lineal es única.

Dados v1=[3,2]v_1 = [-3,2] y v2=[3,4]v_2 = [3,4], hallar la combinación lineal m=3v1+2v2m = 3v_1 + 2v_2

m=3[3,2]+2[3,4]=[9,6]+[6,8]m = 3[-3,2] + 2[3,4] = [-9,6] + [6,8]\Longrightarrow m=[9+6,6+8]m=[3,14]m = [-9+6, 6+8]\Longrightarrow m = [-3,14]

En la gráfica de la página siguiente se aprecia la combinación lineal de los vectores v1v_1 y v2v_2.

139

¿Existirán escalares C1C_1 y C2C_2 de tal manera que el vector

v=[2,6]v = [2, 6]

sea la combinación lineal de los vectores

v1=[4,1]yv2=[1,2]v_1 = [4, 1] \quad y \quad v_2 = [-1, 2] [26]=C1[41]+C2[12]\begin{bmatrix} -2 \\ 6 \end{bmatrix}=C_1\begin{bmatrix} 4\\ 1 \end{bmatrix}+C_2\begin{bmatrix} -1 \\ 2 \end{bmatrix}

Establecemos el sistema de ecuaciones:

{4C1C2=2C1+2C2=6\begin{cases} 4C_1 - C_2 = -2 \\ C_1 + 2C_2 = 6 \end{cases}

La matriz ampliada y la reducción por Gauss es:

[412126] \left [ \begin{array}{rrcr} 4 & -1 & \color{blue}{\vert} &-2\\ 1 & 2 & \color{blue}{\vert} &6\\ \end{array} \right]
140

Reduciendo F2=4F2F1F_2 = 4F_2 - F_1 se obtiene:

[4120926] \left [ \begin{array}{rrcr} 4 & -1 & \color{blue}{\vert} &-2\\ 0 & 9 & \color{blue}{\vert} &26\\ \end{array} \right]

Escribiendo el nuevo sistema equivalente:

{1)4C1C2=22)9C2=26\begin{cases} 1)&4C_1 - C_2 = -2 \\ 2)&9C_2 = 26 \end{cases}

Despejando en C2C_2 en la ecuación (2) y reemplazando el valor obtenido en la ecuación (1) se obtiene que C1=2/9C_1 = 2/9 y que C2=26/9C_2 = 26/9.

Podemos expresar que los vectores serían combinación lineal cuando v1v_1 es multiplicado por C1C_1 y v2v_2 es multiplicado por C2C_2.

Vamos a determinar si los vectores v1=[1,2,1]v_1 = [1,2,-1], el vector v2=[0,1,2]v_2 = [0,-1,2] y el vector v3=[1,1,1]v_3 = [1,1,-1] pueden ser la combinación lineal del vector descrito como v=[3,7,9]v = [-3,-7,9].

a[121]+b[012]+c[111]=[379] {\color{red}a}\begin{bmatrix} \quad1 \\ \quad2 \\ -1 \end{bmatrix}+{\color{red}b}\begin{bmatrix} \quad0 \\ -1 \\ \quad2 \end{bmatrix}+{\color{red}c}\begin{bmatrix} \quad1 \\ \quad1 \\ -1 \end{bmatrix}=\begin{bmatrix} -3 \\ -7\\ \quad9 \end{bmatrix}

En el vídeo veremos la solución a la inquietud planteada, trataremos de hallar los valores de aa, bb y cc que permitan verificar esa combinación lineal.

Plantearemos el sistema de ecuaciones que describe la situación, además, la matriz ampliada generada:

141
a+c=32ab+c=7a+2bc=9}\left. \begin{array}{rcr} a+c & = & -3 \\ 2a -b+c & = & -7 \\ -a+2b-c & = & 9 \end{array} \right\}\Longrightarrow [101321171219] \left [ \begin{array}{rrrcr} 1 & 0 & 1 &\color{red}{\vert} &-3\\ 2 & -1 & 1 &\color{red}{\vert} &-7\\ -1 & 2 & -1 &\color{red}{\vert} &9\\ \end{array} \right]

Dependencia e independencia lineal

Dos vectores en un espacio vectorial son linealmente dependientes si y sólo si uno de ellos es un múltiplo escalar del otro.

142

Vectores linealmente dependientes. Varios vectores libres del plano se dice que son linealmente dependientes si hay una combinación lineal de ellos que es igual al vector cero, sin que sean cero todos los coeficientes de la combinación lineal. También se cumple el reciproco: si un vector es combinación lineal de otros, entonces todos los vectores son linealmente dependientes.

Propiedades

  1. Si varios vectores son linealmente dependientes, entonces al menos uno de ellos se puede expresar como combinación lineal de los demás.
  2. a1v1+a2v2+...+anvn=0(A1)\tag{A1} a_1v_1+a_2v_2+ ... + a_nv_n=0 v1=a2a1v2a3a1v3v_1=-\frac{a_2}{a_1}v_2-\frac{a_3}{a_1}v_3
  3. Dos vectores del plano son linealmente dependientes si, y sólo si, son paralelos.
  4. Dos vectores libres del plano u=[u1,u2]u = [u_1,u_2] y v=[v1,v2]v = [v_1,v_2] son linealmente dependientes si sus componentes son proporcionales.
143

Si un determinado conjunto de vectores son linealmente dependientes se dice que son vectores LDLD.

Si un determinado conjunto de vectores son linealmente independientes se dice que son vectores LILI.

Seanv1,v2,v3...vkRSean \quad v_1, v_2, v_3 ... v_k \in ℝ a1v1+a2v2+...+akvk=0 a_1v_1+a_2v_2+ ... + a_kv_k = 0 a1[v1v2]+a1[u1u2]+...+ak[w1w2]=[00](A2)\tag{A2} {\color{red}a_1} \begin{bmatrix} v_1 \\ v_2 \\ \end{bmatrix}+{\color{red}a_1} \begin{bmatrix} u_1 \\ u_2 \\ \end{bmatrix}+...+{\color{red}a_k} \begin{bmatrix} w_1 \\ w_2 \\ \end{bmatrix}= \begin{bmatrix} 0\\ 0\\ \end{bmatrix} a1[v1v2v3]+a2[u1u2u3]+...+ak[w1w2w3]=[000]{\color{red}a_1} \begin{bmatrix} v_1 \\ v_2 \\ v_3 \end{bmatrix}+{\color{red}a_2} \begin{bmatrix} u_1 \\ u_2 \\ u_3 \end{bmatrix}+...+{\color{red}a_k} \begin{bmatrix} w_1 \\ w_2 \\ w_3 \end{bmatrix}= \begin{bmatrix} 0\\ 0\\ 0 \end{bmatrix}

Si el sistema de ecuaciones así formado indicará que los vectores son LDLD si y solo sí presenta soluciones no triviales.

Si el sistema de ecuaciones establecido NO tiene solución diremos que los vectores son LILI.

Un conjunto de nn vectores en Rmℝ^m es linealmente independiente si n>mn > m.

Vamos a determinar si v1=[2,2]v_1 = [2,2] y v2=[3,4]v_2 = [-3,4] son LILI o son LDLD.

a1[22]+a2[34]=[00]{\color{red}a_1} \begin{bmatrix} 2 \\ 2 \\ \end{bmatrix}+{\color{red}a_2} \begin{bmatrix} -3\\ \quad 4\\ \end{bmatrix}= \begin{bmatrix} 0\\ 0\\ \end{bmatrix}
144

Ahora emplearemos el método de reducción de Gauss para la matriz ampliada:

[230240]F2=F2F1[230070] \left [ \begin{array}{rrrcr} 2 & -3 &\color{red}{\vert} &0\\ 2 & 4 &\color{red}{\vert} &0\\ \end{array} \right] F_2=F_2-F_1\left [ \begin{array}{rrrcr} 2 & -3 &\color{red}{\vert} &0\\ 0 & 7 &\color{red}{\vert} &0\\ \end{array} \right]

Esta reducción nos expresa en F2F_2 que 7c2=07c_2 = 0 y por lo tanto c2=0c_2 = 0. Reemplazando en F1F_1 nos dice que c2=0c_2 = 0. El resultado así obtenido es la solución trivial que determina que los vectores dados son LILI.

Determinar la dependencia o independencia lineal de los siguientes vectores:

[2412];[102];[416]\begin{bmatrix} -2 \\ \quad 4\\ -12 \end{bmatrix}; \begin{bmatrix} 1 \\ 0\\ 2 \end{bmatrix}; \begin{bmatrix} 4\\ 1\\ 6 \end{bmatrix}

Ahora vamos a formalizar la ecuación teórica con los escalares aa, bb y cc. A continuación nuestra sistema de ecuaciones que resolveremos por la regla de Cramer:

a[2412]+b[102]+c[416]=[000]{\color{red}a}\begin{bmatrix} -2 \\ \quad 4\\ -12 \end{bmatrix}+{\color{red}b} \begin{bmatrix} 1 \\ 0\\ 2 \end{bmatrix}+{\color{red}c} \begin{bmatrix} 4\\ 1\\ 6 \end{bmatrix}= \begin{bmatrix} 0 \\ 0 \\ 0 \end{bmatrix} 2a+b+4c=04a+c=012a+2b+6c=0} \left. \begin{array}{rcr} -2a+b+4c & = & 0 \\ 4a+c & = & 0 \\ -12a+2b+6c & = & 0 \end{array} \right\}

Aplicando Sarrus calculemos el determinante del sistema, necesario para aplicar Cramer:

145
21440112262140122\begin{vmatrix} -2&1&4\\ \quad4&0&1\\ -12&2&6 \end{vmatrix}\begin{matrix} -2 & 1 \\ \quad4 & 0\\ -12 & 2 \end{matrix} det=(2)(0)(6)+(1)(1)(12)+(4)(4)(2)det={(-2)(0)(6)+(1)(1)(-12)+(4)(4)(2)}- (12)(0)(4)+(2)(1)(2)+(6)(4)(1)=0{(-12)(0)(4)+(2)(1)(-2)+(6)(4)(1)}=0

Podemos concluir que si el determinante de la matriz de coeficientes es igual a 00, el sistema no tiene solución. Por lo tanto los vectores dados son LILI.

Determinar si los polinomios P1(x)=x25x+6P_1(x) = x^2 − 5x + 6, P2(x)=x2xP_2(x) = x^2 − x y P3(x)=x2P_3(x) = x − 2 son LILI o son LDLD.

a(x25x+6)+b(2x)+c(x2)=0\color{red}{a} (x^2 − 5x + 6)+\color{red}{b} (^2 − x)+\color{red}{c} (x − 2)=0

Aplicando la propiedad distributiva, agrupando términos semejantes y factorizando se tiene:

ax25ax+6a+bx2bx+cx2c=0ax^2-5ax+6a+bx^2-bx+cx-2c=0 (ax2+bx2)+(5axbx+cx)+(6a2c)=0(ax^2+bx^2)+(-5ax-bx+cx)+(6a-2c)=0 x2(a+b)+x(5ab+c)+x0(6a2c)=0x2+0x1+0x0x^2(a+b)+x(-5a-b+c)+x^0(6a-2c)=0x^2+0x^1+0x^0

Nuestro sistema de ecuaciones:

a+b=05ab+c=06a2c=0}\left. \begin{array}{rrcr} a+b& = & 0 \\ -5a-b+c & = & 0 \\ 6a -2c & = & 0 \end{array} \right\}
146

Hallemos el determinante de la matriz de coeficientes por Sarrus, necesario para solucionar el sistema por la regla de Cramer:

110511602115160\begin{vmatrix} \quad1&\quad1&\quad0\\ - 5&-1&\quad1\\ \quad6&\quad0&- 2 \end{vmatrix}\begin{matrix} \quad1 & \quad1 \\ - 5 & -1\\ \quad 6 & \quad 0 \end{matrix} (1)(1)(2)+(1)(1)(6)+(0)(5)(0){(1)(-1)(2)+(1)(1)(6)+(0)(-5)(0)}- (6)(1)(0)+(0)(1)(1)+(2)(5)(1)=2{(6)(-1)(0)+(0)(1)(1)+(-2)(-5)(1)}=-2 a=0100110020101002=02=0a=\frac{\begin{vmatrix} 0&\quad1&\quad0\\ 0&-1&\quad1\\ 0&\quad0&- 2 \end{vmatrix}\begin{matrix} 0 & \quad1 \\ 0 & -1\\ 0 & \quad 0 \end{matrix}}{-2}=\frac{0}{-2}=0 b=1005016021050602=02=0b=\frac{\begin{vmatrix} \quad1&0&\quad0\\ -5&0&\quad1\\ \quad6&0&- 2 \end{vmatrix}\begin{matrix} \quad 1 & 0 \\ -5 & 0\\ \quad6 & 0 \end{matrix}}{-2}=\frac{0}{-2}=0 c=1105106001151602=02=0c=\frac{\begin{vmatrix} \quad1&\quad1&0\\ -5&-1&0\\ \quad6&\quad0&0 \end{vmatrix}\begin{matrix} \quad 1 & \quad1 \\ -5 & -1\\ \quad6 & \quad0 \end{matrix}}{-2}=\frac{0}{-2}=0

Siendo la solución del sistema trivial podemos decir que los polinomios son LILI.

Determinar la LILI o LDLD de las matrices A1A_1; A2A_2 y A3A_3.

147
A1=[121011];A2=[011111]A3=[212102] A_1=\left[ \begin{array}{rrr} 1 & 2 & 1 \\ 0 & 1 & 1 \\ \end{array} \right];A_2=\left[ \begin{array}{rrr} 0 & 1 & 1 \\ 1 & 1 & 1 \\ \end{array} \right]A_3=\left[ \begin{array}{rrr} 2 & 1 & 2 \\ 1 & 0 & 2 \\ \end{array} \right]

Establecemos la ecuación matricial y reemplazamos:

a[121011]+b[011111]+c[212102]=[000000] {\color{red}{a}}\left[ \begin{array}{rrr} 1 & 2 & 1 \\ 0 & 1 & 1 \\ \end{array} \right]+{\color{red}{b}}\left[ \begin{array}{rrr} 0 & 1 & 1 \\ 1 & 1 & 1 \\ \end{array} \right]+{\color{red}{c}}\left[ \begin{array}{rrr} 2 & 1 & 2 \\ 1 & 0 & 2 \\ \end{array} \right]=\left[ \begin{array}{rrr} 0 & 0 & 0 \\ 0 & 0 & 0 \\ \end{array} \right] [a2aa0aa]+[0bbbbb]+[2cc2cc02c]=[000000] \left[ \begin{array}{rrr} a & 2a & a \\ 0 & a & a \\ \end{array} \right]+\left[ \begin{array}{rrr} 0 & b & b \\ b & b & b \\ \end{array} \right]+\left[ \begin{array}{rrr} 2c & c & 2c \\ c & 0 & 2c \\ \end{array} \right]=\left[ \begin{array}{rrr} 0 & 0 & 0 \\ 0 & 0 & 0 \\ \end{array} \right]

Sumando A1A_1; A2A_2 y A3A_3 y reescribiendo le ecuación matricial:

[a+2c2a+b+ca+b+2cb+ca+ba+b+2c]=[000000] \left[ \begin{array}{rrr} a+2c & 2a+b+c & a+b+2c \\ b+c & a+b & a+b+2c \\ \end{array} \right]=\left[ \begin{array}{rrr} 0 & 0 & 0 \\ 0 & 0 & 0 \\ \end{array} \right]

Hemos obtenido un sistema homogéneo de seis (6) ecuaciones con las variables aa, bb y cc.

E1a+2c=0E22a+b+c=0{\color{red}{E_1}} \quad a+2c=0 \qquad {\color{red}{E_2}} \quad 2a+b+c=0 E3a+b+2c=0E4b+c=0{\color{red}{E_3}} \quad a+b+2c=0 \qquad {\color{red}{E_4}}\quad b+c=0 E5a+b=0E6a+b+2c=0{\color{red}{E_5}} \quad a+b=0 \qquad {\color{red}{E_6}}\quad a+b+2c=0

Solución:

Si de (4) despejamos bb y reemplazamos en (2) obtenemos que 2a=02a = 0, entonces a=0a = 0. Si de (1) despejamos aa y reemplazamos en (3) que es igual a (6) obtenemos que b=0b = 0, y si reemplazamos b=0b = 0 en (4) obtenemos que c=0c = 0. El sistema tiene una solución trivial lo que indica que A1A_1; A2A_2 y A3A_3 son LILI.

148

Generación de espacios vectoriales

Conjunto generador. Los vectores v1v_1, v2v_2, ......, vnv_n de un espacio vectorial VV generan a VV si todo vector en VV se puede escribir como una combinación lineal de los mismos. Para todo vector perteneciente al espacio vectorial VV existen escalares α1α_1, α2α_2, ......, α_n$ tales que:

v=α1v1+α2v2+...+αnvn(A2)\tag{A2}v = α_1v_1 + α_2v_2 + ... + α_nv_n

Conjuntos de vectores generadores en R2ℝ^2 y R3ℝ^3

Conjunto generador en R2ℝ^2
i=[10]i=[10]j=[01]j=[01](A3)\tag{A3} i=\begin{bmatrix} 1\\ 0\\ \end{bmatrix}\quad -i=\begin{bmatrix} -1\\ \quad 0\\ \end{bmatrix}\quad j=\begin{bmatrix} 0\\ 1\\ \end{bmatrix}\quad -j=\begin{bmatrix} \quad0\\ -1\\ \end{bmatrix}
149
Conjunto generador en R3ℝ^3
i=[100]i=[100]j=[010]j=[010] i=\begin{bmatrix} 1\\ 0\\ 0\\ \end{bmatrix}\quad -i=\begin{bmatrix} -1\\ \quad 0\\ \quad 0\\ \end{bmatrix} \quad j=\begin{bmatrix} 0\\ 1\\ 0\\ \end{bmatrix}\quad -j=\begin{bmatrix} \quad 0\\ -1\\ \quad 0\\ \end{bmatrix} z=[001]z=[001](A4) \tag{A4} z=\begin{bmatrix} 0\\ 0\\ 1\\ \end{bmatrix}\quad -z=\begin{bmatrix} \quad 0\\ \quad 0\\ -1\\ \end{bmatrix}

Espacio generado por un conjunto de vectores.

Sean v1v_1, v2v_2, ......, vkv_k, kk vectores de un espacio vectorial VV. El espacio generado por {v1,v2,...,vk}\{v_1, v_2, ... , v_k\} es el conjunto de combinaciones lineales v1v_1, v2v_2, ......, vkv_k. Teniendo presente que α1α_1, α2α_2, ... , αkα_k son escalares arbitrarios tenemos:

150
gen[v1,v2,...,vk]=gen[α1v1+α2v2+...+αkvk](A5) \tag {A5}gen[v_1, v_2,...,v_k]=gen[α_1v_1+α_2v_2+...+α_kv_k]

Para solucionar problemas acerca del conjunto generador de un espacio vectorial se requiere conocer su rango.

El rango de una matriz coincide con el número de vectores independientes que conforman la matriz.

El rgrg de la matriz es igual al número de filas distintas de cero en la matriz después de reducirla a la forma escalonada de la fila usando transformaciones elementales sobre las filas de la matriz.

Veamos en este ejemplo de rango de una matriz M2x3M_{2x3}:

A=[120132]F2=F2F1[120012]rg=2 A=\left[ \begin{array}{rrr} -1 & 2 & 0 \\ -1 & 3 & 2 \\ \end{array} \right] F_2=F_2-F_1 \left[ \begin{array}{rrr} -1 & 2 & 0 \\ \quad 0 & 1 & 2 \\ \end{array} \right] ∴ {\color{red}{rg=2}}

Debido a que AA es una matriz 2×32 \times 3 y luego de la única transformación elemental realizada ninguna fila está formada por CEROS, el rango es igual al número de filas inicial.

Determinemos el rango de la matriz AA

A=[333124001033312] A=\left[ \begin{array}{rrrrr} 3 & -3 & -3 & -1 & 2\\ 4 & 0 & 0 & -1 & 0\\ -3 & 3 & 3 & 1 & 2 \\ \end{array} \right]

Realizar las transformaciones elementales donde la primera a realizar es para lograr que a21a_{21} sea igual a CERO:

F2=F243F1F_2=F_2-\frac{4}{3}F_1
151
A=[33312044138333312]F3=F3F1 A=\left[ \begin{array}{rrrrr} 3 & -3 & -3 & -1 & 2\\ 0 & 4 & 4 & \frac{1}{3} & -\frac{8}{3} \\ -3 & 3 & 3 & 1 & 2 \\ \end{array} \right] \rightarrow F_3=F_3-F_1 A=[33312044138300000]rg=2 A=\left[ \begin{array}{rrrrr} 3 & -3 & -3 & -1 & 2\\ 0 & 4 & 4 & \frac{1}{3} & -\frac{8}{3} \\ 0 & 0 & 0 & 0 & 0 \\ \end{array} \right] ∴ {\color{red}{rg=2}}

Al ser la última fila compuesta de CEROS, el rango es el número de filas diferentes a ese valor o sea que el rango es igual a 22.

Ya definido que es el rango de una matriz retomemos el la idea de conjunto generador de un espacio vectorial.

Normalmente el conjunto generador se denota por la expresión S=[v1,v2,...,vn]S = [v_1, v_2, ... , v_n] o también como gen[v1,v2,...,vn]gen[v_1, v_2, ... , v_n].

Si gen(S)=Vgen(S) = V, entonces se dice que VV es generado por [v1,v2,...,vn][v_1, v_2, ..., v_n] o que SS genera a VV.

Si S=[v1,v2,...,vn]S = [v_1, v_2, ... , v_n] es un conjunto de vectores en un espacio vectorial VV, entonces gen(S)gen(S) es el menor sub-espacio de VV que contiene a SS; ello implica que cualquier otro sub-espacio de VV que contenga a SS debe contener a gen(S)gen(S).

Determinación de gen(S)gen(S)

Determinar un conjunto generador requiere que:

(1)[v1,v2,...,vk]Rn {\color{red}{(1)}}\qquad [v_1, v_2,...,v_k]\subset R^n (2)rg=[v1,v2,...,vk]=n {\color{red}{(2)}}\qquad rg=[v_1, v_2,...,v_k]=n

152

Si solamente tenemos un vector vv, entonces el conjunto generado por vv es el conjunto de todos los vectores paralelos a vv o sea que tienen la forma αvαv:

gen(v)=αvαR(A6) \tag{A6} gen(v) = {α v | α ∈ ℝ}

De manera general podemos expresar que si vv es un es un vector diferente de CERO en Rnℝ^n, el conjunto generado es una recta que pasa por el origen.

Cuando se tienen dos vectores v1v_1 y v2v_2, entonces el conjunto generado por esos dos vectores tiene la forma:

gen(v1,v2)=αv1+βv2α,βR(A7)\tag{A7} gen(v_1, v_2) = {α v1 + β v2| α, β ∈ ℝ}

En general si se tienen dos vectores v1v_1 y v2v_2 diferentes de CERO, el conjunto generado es un plano que pasa por el origen.

Encontremos por medio de restricciones en las entradas de los vectores el conjunto generado por v1v_1 y v2v_2. Buscaremos las condiciones de un vector v=[a,b,c]v = [a,b,c] tal que haya solución al sistema.

gen{v1=[121];v2=[203]} gen\begin{Bmatrix} v_1=\begin{bmatrix} 1 \\ 2 \\ 1 \\ \end{bmatrix}; & v_2=\begin{bmatrix} 2 \\ 0 \\ 3 \end{bmatrix} \\ \end{Bmatrix} [abc]=c1[121]+c2[203]{1)C1+2C2=a2)2C1=b3)C1+3C2=c \begin{bmatrix} a \\ b \\ c\\ \end{bmatrix}=c_1\begin{bmatrix} 1 \\ 2 \\ 1\\ \end{bmatrix}+c_2\begin{bmatrix} 2 \\ 0 \\ 3\\ \end{bmatrix}\Longrightarrow \begin{cases} 1)&C_1 + 2C_2 = a\\ 2)& 2C_1= b\\ 3)&C_1+3C_2=c \end{cases}
153

Veamos la matriz de coeficientes y la respectiva matriz ampliada BB:

A=[122013] A=\left[ \begin{array}{rrr} 1 & 2 \\ 2 & 0 \\ 1 & 3 \\ \end{array} \right] AB=[12a20b13c]A|B= \left [\begin{array}{rrrcr} 1 & 2 &\color{red}{\vert} &a\\ 2 & 0 &\color{red}{\vert} &b\\ 1 & 3 &\color{red}{\vert} &c\\ \end{array} \right]

Ahora realizamos las transformaciones elementales de la matriz ampliada y establecemos su rango, si es diferente al rango de la matriz de coeficientes que es igual a 33 el sistema no tendrá solución:

AB=[12a20b13c]F2=F22F1A|B= \left [\begin{array}{rrrcr} 1 & 2 &\color{red}{\vert} &a\\ 2 & 0 &\color{red}{\vert} &b\\ 1 & 3 &\color{red}{\vert} &c\\ \end{array} \right] \rightarrow F_2 = F_2-2F_1 AB=[12a042a+b13c]F3=F3F1A|B= \left [\begin{array}{rrrcr} 1 & 2 &\color{red}{\vert} & a\\ 0 & -4 &\color{red}{\vert} & -2a+b\\ 1 & 3 &\color{red}{\vert} &c\\ \end{array} \right] \rightarrow F_3 = F_3-F_1 AB=[12a042a+b006a+b+4c]A|B= \left [\begin{array}{rrrcr} 1 & 2 &\color{red}{\vert} &a\\ 0 & -4 &\color{red}{\vert} &-2a+b\\ 0 & 0 &\color{red}{\vert} &-6a+b+4c\\ \end{array} \right]

Observamos que NO existe ninguna fila formado por ceros únicamente en la matriz ampliada, por lo tanto su rango es 33, y el rango de la matriz escalonada es 22:

rg(A)>rg(AB)rg(A)<rg(AB)rg(A)>rg(A|B) \lor rg(A)<{rg(A|B)}, es un sistema incompatible.

154

En la tercera fila encontramos una contradicción y la forma de evitarla es igualando a CERO, es una restricción. Es equivalente a decir que 0=6a+b+4c0 = − 6a + b + 4c es la ecuación de un plano que pasa por el origen.

En la gráfica vemos ese plano que satisface el conjunto generado con restricciones:

Conjunto generado con restricciones.

Vamos a determinar si los vectores v1=[1,1,0]v_1 = [1,1,0]; v2=[1,0,1]v_2 = [1,0,1] y v3=[0,1,1]v_3 = [0,1,1] constituyen un sistema generador en R3ℝ^3.

El primer paso es establecer un vector R3ℝ^3 genérico =[a,b,c]= [a,b,c].

Luego establecemos nuestra ecuación matricial, el sistema de ecuaciones generado y la respectiva matriz ampliada:

155
C1[110]+C2[011]+C3[101]=[abc] C_1\begin{bmatrix} 1 \\ 1 \\ 0\\ \end{bmatrix}+C_2\begin{bmatrix} 0 \\ 1 \\ 1\\ \end{bmatrix}+C_3\begin{bmatrix} 1 \\ 0 \\ 1\\ \end{bmatrix}=\begin{bmatrix} a \\ b\\ c\\ \end{bmatrix} {1)C1+C3=a2)C1+C2=b3)C2+C3=c \begin{cases} 1)&C_1 + C_3 = a\\ 2)&C_1 + C_2 = b\\ 3)&C_2+C_3 =c \end {cases} \Longrightarrow [101a110b011c] \left [ \begin{array}{rrrcr} 1 & 0 & 1 &\color{red}{\vert} &a\\ 1 & 1 & 0 &\color{red}{\vert} &b\\ 0 & 1 & 1 &\color{red}{\vert} &c\\ \end{array} \right]

Nuestro próximo paso es obtener una matriz superior mediante la reducción de Gauss: F2=F2F1[101a011a+b011c] F_2=F_2-F_1 \to \left [ \begin{array}{rrrcr} 1 & 0 & 1 &\color{red}{\vert} &a\\ 0 & 1 & -1 &\color{red}{\vert} &-a+b\\ 0 & 1 & 1 &\color{red}{\vert} &c\\ \end{array} \right] F3=F3F2[101a011a+b002ab+c] F_3=F_3-F_2 \to \left [ \begin{array}{rrrcr} 1 & 0 & 1 &\color{red}{\vert} &a\\ 0 & 1 & -1 &\color{red}{\vert} &-a+b\\ 0 & 0 & 2 &\color{red}{\vert} &a-b+c\\ \end{array} \right]

Lo anterior indica que el rango de la matriz ampliada es 33, el mismo rango de la matriz de coeficientes. Nuestro nuevo sistema de ecuaciones es:

{1)C1+C3=a2)C2C3=a+b3)2C3=ab+c\begin{cases} 1)&C_1 + C_3 = a\\ 2)&C_2 - C_3 = -a+b\\ 3)&2C_3 =a-b+c \end {cases} \Longrightarrow
156

Obtenemos las restricciones de las variables, las que estarán expresadas en términos de las componentes aa, bb y cc:

De la ecuación 3:

C3=ab+c2C_3=\frac{a-b+c}{2}

El valor de C3C_3 en la ecuación 2:

C2ab+c2=a+bC_2-\frac{a-b+c}{2}=-a+b C2=a+b+ab+c2=2a+2b+ab+c2C_2=-a+b+\frac{a-b+c}{2}=\frac{-2a+2b+a-b+c}{2} C2=a+b+c2C_2=\frac{-a+b+c}{2}

El valor de C3C_3 en la ecuación 1:

C1+ab+c2=aC_1+\frac{a-b+c}{2}=a C1=aab+c2=2aa+bc2C_1=a-\frac{a-b+c}{2}=\frac{2a-a+b-c}{2} C1=a+bc2C_1=\frac{a+b-c}{2}

Entonces cualquier vector [a,b,c][a,b,c] es un espacio generado por la combinación lineal expresada como:

[abc]=a+bc2[110]+a+b+c2[101]+ab+c2[011]\begin{bmatrix} a \\ b \\ c\\ \end{bmatrix}=\frac{a+b-c}{2}\begin{bmatrix} 1 \\ 1 \\ 0\\ \end{bmatrix}+\frac{-a+b+c}{2}\begin{bmatrix} 1 \\ 0 \\ 1\\ \end{bmatrix}+\frac{a-b+c}{2}\begin{bmatrix} 0 \\ 1 \\ 1\\ \end{bmatrix}
157

Determinar las componentes del vector v1=[3,2,3]v_1 = [3,-2,-3] en la base BB generada por los vectores v2=[0,1,1]v_2 = [0, 1, 1]; v3=[1,1,1]v_3 = [-1,-1, 1] y v4=[0,1,1]v_4 = [0, -1, 1], es decir si son una combinación lineal que genera un espacio vectorial.

Paso 1. Verificar si es un conjunto generador de un vector genérico [a,b,c][a,b,c], ¿existen valores para las constantes xx, yy y zz?

a[011]+b[111]+c[011]=[xyz]a\begin {bmatrix} 0 \\ 1 \\ 1\\ \end{bmatrix}+b\begin{bmatrix} -1 \\ -1 \\ 1\\ \end{bmatrix}+c\begin{bmatrix} 0 \\ -1 \\ 1\\ \end{bmatrix}=\begin{bmatrix} x \\ y \\ z\\ \end{bmatrix} [0aa]+[bbb]+[0cc]=[xyz]\begin {bmatrix} 0 \\ a \\ a\\ \end{bmatrix}+\begin{bmatrix} -b \\ -b \\ b\\ \end{bmatrix}+\begin{bmatrix} 0 \\ -c \\ c\\ \end{bmatrix}=\begin{bmatrix} x \\ y \\ z\\ \end{bmatrix}

El sistema de ecuaciones generado y la matriz ampliada de ese sistema:

{b=xabc=ya+b+c=z\begin{cases} & - b = x\\ & a - b - c = y\\ & a + b + c =z \end {cases} [010x111y111z] \left [ \begin{array}{rrrcr} 0 & -1 & 0 &\color{red}{\vert} &x\\ 1 & -1 & -1 &\color{red}{\vert} &y\\ 1 & 1 & 1 &\color{red}{\vert} &z\\ \end{array} \right]

Aplicamos Gauss para determinar el valor de xx, yy y zz:

[010x111y111z]F3F1 \left [ \begin{array}{rrrcr} 0 & -1 & 0 &\color{red}{\vert} &x\\ 1 & -1 & -1 &\color{red}{\vert} &y\\ 1 & 1 & 1 &\color{red}{\vert} &z\\ \end{array} \right]\updownarrow \frac{F_3}{F_1}
158
[111z111y010x]F2=F2F1 \left [ \begin{array}{rrrcr} 1 & 1 & 1 &\color{red}{\vert} &z\\ 1 & -1 & -1 &\color{red}{\vert} &y\\ 0 & -1 & 0 &\color{red}{\vert} &x\\ \end{array} \right] \to F_2=F_2-F_1 [111z022yz010x]F3=2F3F2 \left [ \begin{array}{rrrcr} 1 & 1 & 1 &\color{red}{\vert} &z\\ 0 & -2 & -2 &\color{red}{\vert} &y-z\\ 0 & -1 & 0 &\color{red}{\vert} &x\\ \end{array} \right] \to F_3=2F_3-F_2 [111z022yz0022xy+z]F2=F22 \left [ \begin{array}{rrrcr} 1 & 1 & 1 &\color{red}{\vert} & z\\ 0 & -2 & -2 &\color{red}{\vert} & y-z\\ 0 & 0 & 2 &\color{red}{\vert} & 2x-y+z\\ \end{array} \right] \to F_2= \frac{F_2}{-2} [111z011yz20022xy+z]F3=F32 \left [ \begin{array}{rrrcr} 1 & 1 & 1 &\color{red}{\vert} & z\\ 0 & 1 & 1 &\color{red}{\vert} & \Large \frac{y-z}{-2}\\ 0 & 0 & 2 &\color{red}{\vert} & 2x-y+z\\ \end{array} \right] \to F_3= \frac{F_3}{2} [111z011yz20012xy+z2] \left [ \begin{array}{rrrcr} 1 & 1 & 1 &\color{red}{\vert} & z\\ 0 & 1 & 1 &\color{red}{\vert} & \Large \frac{y-z}{-2}\\ 0 & 0 & 1 &\color{red}{\vert} & \Large \frac{2x-y+z}{2}\\ \end{array} \right] {a+b+c=zb+c=yz2c=2xy+z2 \begin{cases} & a + b + c =\Large z\\ \\ & b + c = \Large \frac {y-z}{-2}\\ \\ & c = \Large \frac{2x-y+z}{2} \end {cases}
159
{a=zbcb=yz2cc=2xy+z2\begin{cases} & \Large a = z - b - c\\ \\ & \Large b = \frac {y-z}{-2}-c\\ \\ & \Large c = \frac{2x-y+z}{2} \end {cases}

Ahora reemplazando cc en bb y luego cc y bb en aa tenemos:

{a=zbcb=yz22xy+z2=x+yzc=2xy+z2\begin{cases} & \Large a = z - b - c\\ \\ & \Large b = \frac {y-z}{-2}- \color{blue}{\frac{2x-y+z}{2}}= \color{blue}{-x+y-z}\\ \\ & \Large c = \frac{2x-y+z}{2} \end {cases} {a= z(x+yz)2xy+z2=y+3z2b=yz22xy+z2=x+yzc=2xy+z2\begin{cases} & \Large a =\ z -(-x+y-z) - \frac{2x-y+z}{2}= \color{blue}{\frac{-y+3z}{2}}\\ \\ & \Large b = \frac {y-z}{-2}- {\frac{2x-y+z}{2}}= \large \color{blue}{-x+y-z}\\ \\ & \Large c = \color{blue}{\frac{2x-y+z}{2}} \end {cases}

Podemos expresar entonces que:

[xyy]=y+3z2[011]+(x+yz)[111]+2xy+z2[011]\begin{bmatrix} x \\ y \\ y\\ \end{bmatrix}= \color{red}{\frac{-y+3z}{2}} \color{black} \begin{bmatrix} 0 \\ 1 \\ 1\\ \end{bmatrix}+\color{red}{(-x+y-z)} \color{black} \begin{bmatrix} -1 \\ -1 \\ \quad 1\\ \end{bmatrix}+\color{red}{\frac{2x-y+z}{2}} \color{black} \begin{bmatrix} \quad 0 \\ -1 \\ \quad 1\\ \end{bmatrix}
160

Paso 2. Hallar las coordenadas del vector v1=[3,2,3]v_1 = [3,-2,-3] en la base generada por los vectores dados.

[323]=α[011]+β[111]+γ[011]\begin{bmatrix} \quad 3 \\ -2 \\ -3\\ \end{bmatrix}= \color{red}{\alpha} \color{black} \begin{bmatrix} 0 \\ 1 \\ 1\\ \end{bmatrix}+\color{red}{\beta} \color{black} \begin{bmatrix} -1 \\ -1 \\ \quad 1\\ \end{bmatrix}+\color{red}{\gamma} \color{black} \begin{bmatrix} \quad 0 \\ -1 \\ \quad 1\\ \end{bmatrix} {β=3αβγ=2α+β+γ=3{β=3α=52γ=52\begin{cases} & \large -\beta = 3 \\ & \large \alpha - \beta - \gamma = -2\\ & \large \alpha + \beta + \gamma = -3\\ \end {cases} \Longrightarrow \begin{cases} & \large \beta = -3 \\ & \large \alpha = - \frac{5}{2}\\ & \large \gamma = \frac{5}{2}\\ \end {cases} [323]=[35252]\begin{bmatrix} \quad 3\\ -2 \\ -3\\ \end{bmatrix}= \begin{bmatrix} \color{red}{-3}\\ \\ \color{red}{\Large -\frac{5}{2}}\\ \\ \color{red}{\Large \quad \frac{5}{2}}\\ \end{bmatrix}

Son las coordenadas del vector especificado.

4.3 Criptografía

Introducción. Según lo expresa el diccionario de la Real Academia Española (RAE) criptografía es una combinación de dos vocablos griegos: crypto, que significa “oculto” o “secreto”, y grapho, que significa “escritura”.

Así, la criptografía es el arte de escribir en códigos o de manera enigmática. Un sistema criptográfico de comunicación implica que exista una codificación y una decodificación que tanto el remitente como el destinatario del mensaje conozcan:

161

Una regla de correspondencia específica entre un conjunto de símbolos (como las letras del alfabeto y los signos de puntuación que componen los mensajes) y un conjunto de números enteros.

Una matriz AA no singular especificada, recordando que una matriz es singular si y sólo si su determinante es nulo.

La criptografía es la ciencia que se encarga de diseñar métodos para mantener confidencial a la información que es enviada por un medio inseguro. Casi todos los medios de comunicación son inseguros, es decir, un espía siempre puede intervenir una comunicación, y en tal caso conocer su contenido, alterar el contenido, borrar el contenido, etc.

LA MATEMÁTICA Y LOS MENSAJES SECRETOS

Lester Hill autor del cifrado de sustitución poligráfica.

Lester Hill en el año 1929 propuso un método que consiste en cifrar en bloque mediante el empleo de matrices y otras operaciones algebraicas, el cual surgió como respuesta a la gran vulnerabilidad que presentan los cifrados monográficos, como el cifrado afín, estos eran desencriptados con facilidad empleando técnica estadística.

162

Lo anterior es sólo una posibilidad, también podrían asignarse las letras en orden decreciente o comenzando por el número 3, etc.

Método de encriptación.

A las letras del alfabeto se le asignan los números del 1 al 27. Al espacio en blanco se le asigna el número 0, para poder separar palabras como se ve en la siguiente tabla, la cual emplearemos para nuestro ejemplo:

Tabla 1. Asignación para encriptar nuestro ejemplo.

Fuente propia de los autores

Cuando se habla de cifrado usando la matemática, las matrices inversas se pueden usar para proporcionar un procedimiento simple y efectivo para codificar y decodificar mensajes. Para codificar un mensaje los elementos que se requieren son: un emisor, un receptor, un mensaje y un código.

La referencia a un código, expresa un método de codificación, se refiere a un algoritmo biunívoco (una función biyectiva) que asigna a cada carácter del mensaje otro carácter. Este método hace que el mensaje enviado por el emisor se transforme en una cadena de símbolos ilegibles para el resto de los receptores que no sean legales.

163

Dependiendo de la calidad del método de codificación, el mensaje será más o menos difícil de descifrar si es capturado por receptores no autorizados.

Cualquier matriz cuyos elementos sean enteros positivos y sea invertible se puede usar como matriz de codificación. Si la matriz de código es de n×nn \times n se construye con el mensaje una matriz de nn filas y tantas columnas como sean necesarias, escribiendo los números por columna y rellenando al final con espacios en blanco si fuera necesario. Luego se multiplica a izquierda por la matriz de código y el resultado es el mensaje codificado.

Para recuperar el mensaje se multiplica la matriz anterior a izquierda por la inversa de la matriz de código. Después, el mensaje es convertido a números dividido en vectores fila sin codificar, cada uno con "nn" elementos, teniendo en cuenta que para codificar el mensaje se elige una matriz invertible AA, de orden n×nn \times n y las matrices renglón no codificada se multiplican por AA, para así obtener las matrices renglón codificadas.

Ilustraremos el proceso con el mensaje:

“En el año 2007 el Instituto Tecnológico Pascual Bravo se transformó en Institución Universitaria, y fue uno de los primeros colegios públicos, así mismo, primer colegio técnico del país, con más de 80 años de existencia”

El mensaje se ha dividido en matrices de 1×31 \times 3, ellas son:

[5145];[12115];[163028];[28512]\left[ \begin{array}{rrr} 5 & 14 & 5 \\ \end{array}\right];\left[ \begin{array}{rrr} 12 & 1 & 15 \\\end{array}\right];\left[ \begin{array}{rrr} 16 & 30 & 28 \\\end{array}\right];\left[ \begin{array}{rrr} 28 & 5 & 12 \\\end{array}\right] [91420];[21921];[222116];[2153]\left[ \begin{array}{rrr} 9 & 14 & 20 \\ \end{array}\right];\left[ \begin{array}{rrr} 21 & 9 & 21 \\\end{array}\right];\left[ \begin{array}{rrr} 22 & 21 & 16 \\\end{array}\right];\left[ \begin{array}{rrr} 21 & 5 & 3 \\\end{array}\right] [141612];[1679];[31617];[1203]\left[ \begin{array}{rrr} 14 & 16 & 12 \\ \end{array}\right];\left[ \begin{array}{rrr} 16 & 7 & 9 \\\end{array}\right];\left[ \begin{array}{rrr} 3 & 16 & 17 \\\end{array}\right];\left[ \begin{array}{rrr} 1 & 20 & 3 \\\end{array}\right]
164
[22112];[2191];[231620];[52119]\left[ \begin{array}{rrr} 22 & 1 & 12 \\ \end{array}\right];\left[ \begin{array}{rrr} 2 & 19 & 1 \\\end{array}\right];\left[ \begin{array}{rrr} 23 & 16 & 20 \\\end{array}\right];\left[ \begin{array}{rrr} 5 & 21 & 19 \\\end{array}\right] [11420];[6169];[15165];[14914]\left[ \begin{array}{rrr} 1 & 14 & 20 \\ \end{array}\right];\left[ \begin{array}{rrr} 6 & 16 & 9 \\\end{array}\right];\left[ \begin{array}{rrr} 15 & 16 & 5 \\\end{array}\right];\left[ \begin{array}{rrr} 14 & 9 & 14 \\\end{array}\right] [21921];[2239];[161422];[14923]\left[ \begin{array}{rrr} 21 & 9 & 21 \\ \end{array}\right];\left[ \begin{array}{rrr} 22 & 3 & 9 \\\end{array}\right];\left[ \begin{array}{rrr} 16 & 14 & 22 \\\end{array}\right];\left[ \begin{array}{rrr} 14 & 9 & 23 \\\end{array}\right] [51920];[12115];[163028];[28512]\left[ \begin{array}{rrr} 5 & 19 & 20 \\ \end{array}\right];\left[ \begin{array}{rrr} 12 & 1 & 15 \\\end{array}\right];\left[ \begin{array}{rrr} 16 & 30 & 28 \\\end{array}\right];\left[ \begin{array}{rrr} 28 & 5 & 12 \\\end{array}\right] [51920];[9211];[1991];[37266]\left[ \begin{array}{rrr} 5 & 19 & 20 \\ \end{array}\right];\left[ \begin{array}{rrr} 9 & 21 & 1 \\\end{array}\right];\left[ \begin{array}{rrr} 19 & 9 & 1 \\\end{array}\right];\left[ \begin{array}{rrr} 37 & 26 & 6 \\\end{array}\right] [22522];[14164];[51216];[201719]\left[ \begin{array}{rrr} 22 & 5 & 22 \\ \end{array}\right];\left[ \begin{array}{rrr} 14 & 16 & 4 \\\end{array}\right];\left[ \begin{array}{rrr} 5 & 12 & 16 \\\end{array}\right];\left[ \begin{array}{rrr} 20 & 17 & 19 \\\end{array}\right] [9135];[191620];[31612];[51316]\left[ \begin{array}{rrr} 9 & 13 & 5 \\ \end{array}\right];\left[ \begin{array}{rrr} 19 & 16 & 20 \\\end{array}\right];\left[ \begin{array}{rrr} 3 & 16 & 12 \\\end{array}\right];\left[ \begin{array}{rrr} 5 & 13 & 16 \\\end{array}\right] [1257];[91621];[5314];[9316]\left[ \begin{array}{rrr} 12 & 5 & 7 \\ \end{array}\right];\left[ \begin{array}{rrr} 9 & 16 & 21 \\\end{array}\right];\left[ \begin{array}{rrr} 5 & 3 & 14 \\\end{array}\right];\left[ \begin{array}{rrr} 9 & 3 & 16 \\\end{array}\right] [9316];[4512];[1719];[20316]\left[ \begin{array}{rrr} 9 & 3 & 16 \\ \end{array}\right];\left[ \begin{array}{rrr} 4 & 5 & 12 \\\end{array}\right];\left[ \begin{array}{rrr} 17 & 1 & 9 \\\end{array}\right];\left[ \begin{array}{rrr} 20 & 3 & 16 \\\end{array}\right] [14131];[2045];[36281];[151620]\left[ \begin{array}{rrr} 14 & 13 & 1 \\ \end{array}\right];\left[ \begin{array}{rrr} 20 & 4 & 5 \\\end{array}\right];\left[ \begin{array}{rrr} 36 & 28 & 1 \\\end{array}\right];\left[ \begin{array}{rrr} 15 & 16 & 20 \\\end{array}\right] [455];[25920];[21514];[391]\left[ \begin{array}{rrr} 4 & 5 & 5 \\ \end{array}\right];\left[ \begin{array}{rrr} 25 & 9 & 20 \\\end{array}\right];\left[ \begin{array}{rrr} 21 & 5 & 14 \\\end{array}\right];\left[ \begin{array}{rrr} 3 & 9 & 1 \\\end{array}\right]

La matriz de codificación es la matriz inversa de la matriz AA que detallamos como:

A=[133143434]A=\left[ \begin{array}{rrr} 1 & 3 & 3 \\ 1 & 4 & 3 \\ 4 & 3 & 4 \\ \end{array} \right]

El mensaje que se pretende leer, ya está dividido en matrices de 1×31 \times 3 y la matriz de codificación es la de la derecha, la matriz AA, así entonces, si X=[x1,x2,...,xn]X = [x_1, x_2, ..., x_n] es una matriz no codificada, al ser multiplicada por la matriz AA invertible, entonces Y=XAY = XA, donde YY es la matriz codificada correspondiente, y la matriz para descifrar el mensaje, quien

165

lo recibe, es la inversa de AA, veamos las cuatro primeras:

[5145]×[133143434]=[248677]\left[ \begin{array}{rrr} 5 & 14 & 5 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 1 & 3 & 3 \\ 1 & 4 & 3 \\ 4 & 3 & 4 \\ \end{array} \right]= \left[ \begin{array}{rrr} 24 & 86 & 77 \\ \end{array} \right] [12115]×[133143434]=[185559]\left[ \begin{array}{rrr} 12 & 1 & 15 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 1 & 3 & 3 \\ 1 & 4 & 3 \\ 4 & 3 & 4 \\ \end{array} \right]= \left[ \begin{array}{rrr} 18 & 55 & 59 \\ \end{array} \right] [163028]×[133143434]=[74252250]\left[ \begin{array}{rrr} 16 & 30 & 28 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 1 & 3 & 3 \\ 1 & 4 & 3 \\ 4 & 3 & 4 \\ \end{array} \right]= \left[ \begin{array}{rrr} 74 & 252 & 250 \\ \end{array} \right] [28512]×[133143434]=[45140147]\left[ \begin{array}{rrr} 28 & 5 & 12 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 1 & 3 & 3 \\ 1 & 4 & 3 \\ 4 & 3 & 4 \\ \end{array} \right]= \left[ \begin{array}{rrr} 45 & 140 & 147 \\ \end{array} \right]

El lector puede verificar las demás multiplicaciones con la matriz AA. Por lo tanto, la serie de matrices renglón cifradas es:

[248677];[185559];[74252250]\left[ \begin{array}{rrr} 24 & 86 & 77 \\ \end{array}\right];\left[ \begin{array}{rrr} 18 & 55& 59 \\\end{array}\right];\left[ \begin{array}{rrr} 74 & 252 & 250 \\\end{array}\right] [45140147];[43143149];[51162174]\left[ \begin{array}{rrr} 45 & 140 & 147 \\ \end{array}\right];\left[ \begin{array}{rrr} 43 & 143 & 149 \\\end{array}\right];\left[ \begin{array}{rrr} 51 & 162 & 174 \\\end{array}\right] [59198193];[299290];[42142138]\left[ \begin{array}{rrr} 59 & 198 & 193 \\ \end{array}\right];\left[ \begin{array}{rrr} 29 & 92 & 90 \\\end{array}\right];\left[ \begin{array}{rrr} 42 & 142 & 138 \\\end{array}\right] [32103105];[36124125];[249275]\left[ \begin{array}{rrr} 32 & 103 & 105 \\ \end{array}\right];\left[ \begin{array}{rrr} 36 & 124 & 125 \\\end{array}\right];\left[ \begin{array}{rrr} 24 & 92 & 75\\\end{array}\right] [35106117];[228567];[59193197]\left[ \begin{array}{rrr} 35 & 106 & 117 \\ \end{array}\right];\left[ \begin{array}{rrr} 22 & 85 & 67 \\\end{array}\right];\left[ \begin{array}{rrr} 59 & 193 & 197\\\end{array}\right]
166
[45156154];[35119125];[41139142]\left[ \begin{array}{rrr} 45 & 156 & 154 \\ \end{array}\right];\left[ \begin{array}{rrr} 35 & 119 & 125 \\\end{array}\right];\left[ \begin{array}{rrr} 41 & 139 & 142\\\end{array}\right] [36124113];[37120125];[51162174]\left[ \begin{array}{rrr} 36 & 124 & 113 \\ \end{array}\right];\left[ \begin{array}{rrr} 37 & 120 & 125 \\\end{array}\right];\left[ \begin{array}{rrr} 51 & 162 & 174\\\end{array}\right] [34105111];[46147161];[44151152]\left[ \begin{array}{rrr} 34 & 105 & 111\\ \end{array}\right];\left[ \begin{array}{rrr} 46 & 147 & 161\\\end{array}\right];\left[ \begin{array}{rrr} 44 & 151 & 152\\\end{array}\right] [3111494];[299688];[69233213]\left[ \begin{array}{rrr} 31 & 114 & 94\\ \end{array}\right];\left[ \begin{array}{rrr} 29 & 96 & 88\\\end{array}\right];\left[ \begin{array}{rrr} 69 & 233 & 213\\\end{array}\right] [49152169];[34118106];[33111115]\left[ \begin{array}{rrr} 49 & 152 & 169\\ \end{array}\right];\left[ \begin{array}{rrr} 34 & 118 & 106\\\end{array}\right];\left[ \begin{array}{rrr} 33 & 111 & 115\\\end{array}\right] [56185187];[279486];[355181185]\left[ \begin{array}{rrr} 56 & 185 & 187\\ \end{array}\right];\left[ \begin{array}{rrr} 27 & 94 & 86\\\end{array}\right];\left[ \begin{array}{rrr} 355 & 181 & 185\\\end{array}\right] [31109105];[247588];[247779]\left[ \begin{array}{rrr} 31 & 109 & 105\\ \end{array}\right];\left[ \begin{array}{rrr} 24 & 75 & 88\\\end{array}\right];\left[ \begin{array}{rrr} 24 & 77 & 79\\\end{array}\right] [46154159];[226980];[2887100]\left[ \begin{array}{rrr} 46 & 154 & 159\\ \end{array}\right];\left[ \begin{array}{rrr} 22 & 69 & 80\\\end{array}\right];\left[ \begin{array}{rrr} 28 & 87 & 100\\\end{array}\right] [216875];[278290];[39120133]\left[ \begin{array}{rrr} 21 & 68 & 75\\ \end{array}\right];\left[ \begin{array}{rrr} 27 & 82 & 90\\\end{array}\right];\left[ \begin{array}{rrr} 39 & 120 & 133\\\end{array}\right] [269785];[299192];[65223196]\left[ \begin{array}{rrr} 26 & 97 & 85\\ \end{array}\right];\left[ \begin{array}{rrr} 29 & 91 & 92\\\end{array}\right];\left[ \begin{array}{rrr} 65 & 223 & 196\\\end{array}\right] [51169173];[144747];[54171182]\left[ \begin{array}{rrr} 51 & 169 & 173\\ \end{array}\right];\left[ \begin{array}{rrr} 14 & 47 & 47\\\end{array}\right];\left[ \begin{array}{rrr} 54 & 171 & 182\\\end{array}\right] [40125134];[134840]\left[ \begin{array}{rrr} 40 & 125 & 134\\ \end{array}\right];\left[ \begin{array}{rrr} 13 & 48 & 40\\\end{array}\right]

Para quienes desconozcan la matriz AA es difícil descifrar el mensaje, sin embargo, para un receptor autorizado que conozca la matriz AA descifrarlo es sencillo, simplemente se debe multiplicar por A1A^{-1} las matrices renglón cifradas para recuperar las matrices no cifradas, por lo tanto, YY (matriz codificada) al ser multiplicada por A1A^{-1} para obtener:

YA1=(XA)A1=X(A8)\tag{A8} YA^{-1} = (XA)A^{-1} = X
167

Así entonces, hay que obtener la matriz inversa, la que ya calculada mediante el método de Gauss Jordan es:

A1=[733110101] A^{-1}=\left[ \begin{array}{rrr} 7 & -3 & -3 \\ -1 & 1 & 0 \\ -1 & 0 & 1 \\ \end{array} \right]

Realizado este proceso, se ha obtenido la serie de matrices ya descifradas, con las que se puede leer el mensaje enviado, veamos las primeras cinco (dejamos la verificación del resto al lector):

[248677]×[733110101]=[5145]\left[ \begin{array}{rrr} 24 & 86 & 77 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 7 & -3 & -3 \\ -1 & 1 & 0 \\ -1 & 0 & 1 \\ \end{array} \right]= \left[ \begin{array}{rrr} 5 & 14 & 5 \\ \end{array} \right] [185559]×[733110101]=[12115]\left[ \begin{array}{rrr} 18 & 55 & 59 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 7 & -3 & -3 \\ -1 & 1 & 0 \\ -1 & 0 & 1 \\ \end{array} \right]= \left[ \begin{array}{rrr} 12 & 1 & 15 \\ \end{array} \right] [74252250]×[733110101]=[163028]\left[ \begin{array}{rrr} 74 & 252 & 250 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 7 & -3 & -3 \\ -1 & 1 & 0 \\ -1 & 0 & 1 \\ \end{array} \right]= \left[ \begin{array}{rrr} 16 & 30 & 28 \\ \end{array} \right] [45140147]×[733110101]=[28512]\left[ \begin{array}{rrr} 45 & 140 & 147 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 7 & -3 & -3 \\ -1 & 1 & 0 \\ -1 & 0 & 1 \\ \end{array} \right]= \left[ \begin{array}{rrr} 28 & 5 & 12 \\ \end{array} \right] [43143149]×[733110101]=[91420]\left[ \begin{array}{rrr} 43 & 143 & 149 \\ \end{array} \right] \times \left[ \begin{array}{rrr} 7 & -3 & -3 \\ -1 & 1 & 0 \\ -1 & 0 & 1 \\ \end{array} \right]= \left[ \begin{array}{rrr} 9 & 14 & 20\\ \end{array} \right]

Las matrices que se han descifrado corresponden a las que se definieron inicialmente al obtener las matrices renglón 1×31 \times 3, y así se

168

logra descifrar el mensaje, al remplazar cada uno de los elementos de la matriz por la letra correspondiente.

Resumen

El remitente del mensaje lo codificará por medio de una matriz AA no singular y el destinatario del mensaje codificado lo decodificará por medio de la matriz A1A^{-1} (única). La matriz AA se llama matriz de codificación y A1A^{-1} se llama matriz de decodificación. El mensaje numérico se escribirá ahora como una matriz MM.

La matriz de codificación AA y la matriz de decodificación A1A^{-1} se elaboran de tal manera que:

  1. AA sea NO singular
  2. AA solo posea entradas de números enteros.
  3. A1A^{-1} solo posea entradas de números enteros.

Codificar el mensaje "EL JUEGO DE TRONOS"

  1. Usando la correspondencia natural de la tabla 1, obtenemos:

51201022571604502119161416205 \quad 12 \quad 0 \quad 10 \quad 22 \quad 5 \quad 7 \quad 16 \quad 0 \quad 4 \quad 5 \quad 0 \quad 21 \quad 19 \quad 16 \quad 14 \quad 16 \quad 20

  1. Ahora escribiremos la matriz MM con 18 entradas (corresponden a los 18 caracteres del mensaje a cifrar), recordemos que una matriz m×nm \times n tiene mnmn entradas.

Por lo tanto podemos emplear cualquiera de las matrices 3×63 \times 6, 6×36 \times 3, 9×29 \times 2 o una matriz 2×92 \times 9.

169

Emplearemos dos matrices M1M_1 y M2M_2 de 3×33 \times 3, ello nos permitirá codificar el mensaje empleando una matriz 3×33 \times 3.

M1=[5120102257160]; M_1= \left[ \begin{array}{rrr} 5 & 12 & 0 \\ 10 & 22 & 5 \\ 7 & 16 & 0 \\ \end{array} \right]; M2=[45021916141620] M_2= \left[ \begin{array}{rrr} 4 & 5 & 0 \\ 2 & 19 & 16 \\ 14 & 16 & 20 \\ \end{array} \right]
  1. Elegimos una matriz No singular (el valor de su determinante es de CERO), analicemos si la siguiente matriz AA cumple el criterio.
A=[123012112] A=\left[ \begin{array}{rrr} 1 & 2 & 3 \\ 0 & -1 & -2 \\ -1 & 1 & 2 \\ \end{array} \right]

Por la regla de Sarrus tenemos:

A=[(1)(1)(2)+(2)(2)(1)+(3)(0)(1)]|A|=[(1)(-1)(2)+(2)(-2)(-1)+(3)(0)(1)][(1)(1)(3)+(1)(2)(1)+(2)(0)(2)]=1-[(-1)(-1)(3)+(1)(-2)(1)+(2)(0)(2)]=1
  1. El mensaje codificado es el resultado de multiplicar M1AM_1A y M2AM_2A:
M1A=[5295347211] M_1A=\left[ \begin{array}{rrr} 5 & -2 & -9 \\ 5 & 3 & -4 \\ 7 & -2 & -11 \\ \end{array} \right]
170
M2A=[432141063250] M_2A=\left[ \begin{array}{rrr} 4 & 3 & 2 \\ -14 & 1 & 0 \\ -6 & 32 & 50 \\ \end{array} \right]

El mensaje codificado contiene las mismas 18 entradas pero algunas de ellas son números enteros negativos lo que dificulta el proceso de decodificación si no se conoce la matriz inversa.

  1. Determinemos la matriz A1A^{-1}. En color rojo las entradas de valor negativo.
a11=123012112=1[(1)(2)(1)(2)]=0a_{11}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=1[(-1)(2)-(1)(-2)]=0 a12=123012112=1[(0)(2)(1)(2)]=2a_{12}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=-1[(0)(2)-(-1)(-2)]=2 a13=123012112=1[(0)(1)(1)(1)]=1a_{13}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=1[(0)(1)-(-1)(-1)]=-1 a21=123012112=1[(2)(2)(1)(3)]=1a_{21}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=-1[(2)(2)-(1)(3)]=-1 a22=123012112=1[(1)(2)(1)(3)]=5a_{22}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=1[(1)(2)-(-1)(3)]=5
171
a23=123012112=1[(1)(1)(1)(2)]=3a_{23}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=-1[(1)(1)-(-1)(2)]=-3 a31=123012112=1[(2)(2)(1)(3)]=1a_{31}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=1[(2)(-2)-(-1)(3)]=-1 a32=123012112=1[(1)(2)(0)(3)]=2a_{32}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=-1[(1)(-2)-(0)(3)]=2 a33=123012112=1[(1)(1)(0)(2)]=1a_{33}=\begin{vmatrix} 1 & \color{red}{2} & 3\\ \color{red}{0} & -1 & \color{red}{-2}\\ -1 & \color{red}{1} & 2\\ \end{vmatrix}=1[(1)(-1)-(0)(2)]=-1 A1=1ACt=11[011252131]=[011252131]A^{-1}=\frac{1}{|A|}C^{t}=\frac{1}{1} \left[ \begin{array}{rrr} 0 & -1 & -1 \\ 2 & 5 & 2 \\ -1 & -3 & -1 \\ \end{array} \right]=\left[ \begin{array}{rrr} 0 & -1 & -1 \\ 2 & 5 & 2 \\ -1 & -3 & -1 \\ \end{array} \right]
  1. Ahora multiplicamos las matrices codificadas por la matriz inversa de AA.

Con la matriz M1M_1

M1AA1=[5295347211][011252131]M_1AA^{-1}= \left[ \begin{array}{rrr} 5 & -2 & -9 \\ 5 & 3 & -4 \\ 7 & -2 & -11 \\ \end{array} \right] \left[ \begin{array}{rrr} 0 & -1 & -1 \\ 2 & 5 & 2 \\ -1 & -3 & -1 \\ \end{array} \right] M1AA1=[5120102257160]M_1AA^{-1}=\left[ \begin{array}{rrr} 5 & 12 & 0 \\ 10 & 22 & 5 \\ 7 & 16 & 0 \\ \end{array} \right]
172

Con la matriz M2M_2

M2AA1=[432141063250][011252131]M_2AA^{-1}= \left[ \begin{array}{rrr} 4 & 3 & 2 \\ -14 & 1 & 0 \\ -6 & 32 & 50 \\ \end{array} \right] \left[ \begin{array}{rrr} 0 & -1 & -1 \\ 2 & 5 & 2 \\ -1 & -3 & -1 \\ \end{array} \right] M1AA1=[45021916141620]M_1AA^{-1}=\left[ \begin{array}{rrr} 4 & 5 & 0 \\ 2 & 19 & 16 \\ 14 & 16 & 20 \\ \end{array} \right]

Observamos que mediante la matriz de decodificación se obtiene el mensaje enviado:

51201022571604502119161416205 \quad 12 \quad 0 \quad 10 \quad 22 \quad 5 \quad 7 \quad 16 \quad 0 \quad 4 \quad 5 \quad 0 \quad 21 \quad 19 \quad 16 \quad 14 \quad 16 \quad 20

E L _ J U E G O _ D E _ T R O N O S

En los dos ejemplos desarrollados se nota que no existe una regla que nos indique como se ha de dividir el mensaje de acuerdo a su correspondencia natural.

En el primero se dividieron en matrices de 3×13 \times 1, en el segundo en matrices de 3×33 \times 3, solo se hace necesario que cumpla las propiedades de la multiplicación de matrices.

Empleando el mismo procedimiento realizado, decodifique el mensaje dada la matriz AA de codificación y BB la matriz codificada. Son sentencias atribuidas al llamado "Hombre Universal": Leonardo Da Vinci.

Para efectos de la realización de los siguientes ejercicios propuestos se ha utilizado la tabla N° 1.

173

Haz clic en la imagen para obtener una vista ampliada del ejercicio propuesto.

Ejercicio 1.

Ejercicio 2.

Ejercicio 3.

Ejercicio 4.

174



176
Bibliografía
Nuestro apoyo teórico

BIBLIOGRAFÍA



[1]Del Valle, J. C.-H. (2012). Álgebra lineal para estudiantes de ingeniería y ciencias. México: Mc Graw-Hill. [2]Grossman, S.I. & FLores, J.J. (2012). Álgebra lineal. México: Mc Graw Hill. [3]Kolman, B. (2013). Álgebra Lineal. México: Pearson Educación. [4]Larson, R. E. (1994). Introducción al Álgebra lineal. México: Limusa. [5]Lay, D. C. (2013). Álgebra lineal para cursos con enfoque por competencias. México: Pearson. [6]Nakos, G., & Joyner, D. (1999). Álgebra lineal con aplicaciones. México: Internacional Thomson Editores. [7]Poole, D. ( 2011). Álgebra lineal una introducción moderna. (3 ed.). México: Cengage Learning. [8]Ramirez, M.M. & Velázquez, M.I. (1996). Introducción al álgebra lineal con geometría analítica y vectorial. Medellín: Universidad Nacional. [9]Williams, G. (2002). Álgebra lineal con aplicaciones. México: McGraw-Hill.
179