Álgebra Lineal
Interactiva

Ejercicios y problemas

Augusto Spela

Álgebra Lineal Interactiva

Ejercicios y problemas



Universidad Tecnológica Nacional - FRBA

Fondo Editorial RED Descartes

Córdoba (España)
2024

Título de la obra:
Álgebra Lineal Interactiva. Ejercicios y problemas.


Autor:
Augusto Spela
Universidad tecnológica Nacional - FRBA
Buenos Aires (Argentina)


Código JavaScript para el libro: Joel Espinosa Longi, IMATE, UNAM.
Recursos interactivos: DescartesJS
Imagen de la portada obtenida de la IA Lexica
Fuentes: Lato y UbuntuMono



Red Educativa Digital Descartes
Córdoba (España)
descartes@proyectodescartes.org
https://proyectodescartes.org

Proyecto iCartesiLibri
https://proyectodescartes.org/iCartesiLibri/index.htm

ISBN: 978-84-18834-86-8


Esta obra está bajo una licencia Creative Commons 4.0 internacional: Reconocimiento-No Comercial-Compartir Igual.

Prefacio

Este libro interactivo no busca ser un libro de cabecera, apenas pretende ser un complemento, enfocado en el "hacer", de los buenos textos que existen en la actualidad sobre Álgebra Lineal. Una lista sugerida de los mismos se ofrece en la bibliografía recomendada

Debo hacer un especial agradecimiento a la UDB Matemática de la Universidad Tecnológica Nacional FRBA, que me permitió utilizar parte del excelente material disponible en su sitio web de Álgebra y Geometría Analítica 1

Algunas de las escenas interactivas fueron tomadas y/o adaptadas de escenas realizadas en DescartesJS por la autoras Elena Álvarez Sáiz y María José García Cebrian [1][2], otras escenas son de autoría de la profesora Melissa Méndez Servin, disponibles en la web [3]. También se han adaptado recursos y applets disponibles en Geogebra.

Profundamente agradecido a los académicos Juan Guillermo Rivera Berrío, Ramiro Antonio Lopera Sánchez, Joel Espinosa Longi y demás desarrolladores del Proyecto iCartesiLibri2, por su invaluable ayuda y aporte de conocimientos, brindados siempre de manera generosa.

Este libro puede leerse descargándolo en un dispositivo digital, pero para visualizar algunas imágenes o videos es necesario estar conectado a Internet.
5
Capítulo I

Vectores, rectas y planos

1.1 Introducción a vectores en R3{\mathbb{R}^3}

Puntos en R3{\mathbb{R}^3}

Para ubicar un punto en R3{\mathbb{R}^3} usaremos como sistema de referencia una terna de ejes perpendiculares entre sí:

  • eje x\color{red}{x} (eje de abscisas, en rojo)
  • eje y\color{green}{y} (eje de ordenadas, en verde)
  • eje z\color{blue}{z} (eje de cotas, en azul)

los cuales se cortan en el punto O (origen de coordenadas).

11

En el siguiente esquema se ven los tres planos que quedan determinados:

  • el plano xy\color{blue}{xy} (en azul)
  • el plano xz\color{green}{xz} ( en verde)
  • el plano yz\color{red}{yz} (en rojo)


Estos planos se conocen como planos coordenados. El nombre del plano xyxy viene de que este plano contiene al eje xx y al eje yy. En forma análoga se derivan los nombres de los otros dos planos.

Se puede demostrar que hay dos formas diferentes de armar un sistema de referencia con tres ejes perpendiculares. Una de esas formas se conoce con el nombre de terna derecha (que es la que usaremos en esta materia y la que hemos presentado recién) y la otra como terna izquierda:

12

1.1.1   Vectores en R3{\mathbb{R}^3}

Queda establecido un sistema de coordenadas donde todo punto de R3{\mathbb{R}^3} se define mediante una terna ordenada de números reales: P(x,y,z)P\left( {x,y,z} \right), y tiene asociado un vector posición: p  =  OP  =(x,  y,  z)\vec p\; = \;\overrightarrow {OP} \; = \left( {x,\;y,\;z} \right)

Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian

1.1.2   Operaciones y nociones básicas de vectores en R3{\mathbb{R}^3}

Sean v=(vx,vy,vz)\vec v = \left( {{v_x},{v_y},{v_z}} \right) y w=(wx,wy,wz)\vec w = \left( {{w_x},{w_y},{w_z}} \right) vectores de R3{\mathbb{R}^3}.

A continuación definimos algunas operaciones y nociones básicas:

  • Igualdad: v  =  w            vx=wx  ,      vy=wy  ,      vz=  wz\vec v\; = \;\vec w\;\;\; \Leftrightarrow \;\;\;{v_x} = {w_x}\;,\;\;\;{v_y} = {w_y}\;,\;\;\;{v_z} = \;{w_z}
  • Suma: v+w=(vx+wx,    vy+wy,    vz+wz)\vec v + \vec w = \left( {{v_x} + {w_x},\;\;{v_y} + {w_y},\;\;{v_z} + {w_z}} \right)
13
  • Vector nulo:   0  =  (0,0,0)\; \vec 0\; = \;\left( {0,0,0} \right)
  • Opuesto de v\vec v: v=(vx,vy,vz) - \vec v = \left( { - {v_x}, - {v_y}, - {v_z}} \right)
  • Resta: vw=v+(w)=(vxwx,vywy,vzwz)\vec v - \vec w = \vec v + \left( { - \vec w} \right) = \left( {{v_x} - {w_x},{v_y} - {w_y},{v_z} - {w_z}} \right)
  • El producto de un escalar por un vector se define:

    v=(vx,vy,vz)    ,    kR      ,k.v=(k.vx,k.vy,k.vz)  \vec v = \left( {{v_x},{v_y},{v_z}} \right)\;\;,\;\;k \in \mathbb{R}\;\;\;,k.\vec v = \left( {k.{v_x},k.{v_y},k.{v_z}} \right)\;

    k.vk.\vec v es un vector tal que:

    • Tiene igual dirección que el vector v\vec v
    • Sentido: Si k>0k > 0 entonces v\vec v y k.vk.\vec v tienen el mismo sentido, si k<0  k < 0\;entonces v\vec v y k.vk.\vec v tienen sentido opuesto. Si   k=0\;k = 0, entonces   0.v=0\;0.\vec v = \vec 0.
    • k.v=kvk.\vec v = \left| k \right|\vec v . El módulo del vector k.vk.\vec v es k\left| k \right| veces el módulo del vector v\vec v.

    ¿Cómo es la longitud del vector k.vk.\vec v respecto de la de v  \vec v\;?

    Si k>1\left| k \right| > 1 entonces k.v>v\lVert k.\vec v\rVert > \lVert \vec v\rVert
    Si k<1\left| k \right| \lt 1 entonces k.v<v\lVert k.\vec v\rVert \lt \lVert \vec v\rVert
    Si k=1\left| k \right| = 1 entonces k.v=v\lVert k.\vec v\rVert = \lVert \vec v\rVert

    Notación

    v:  \lVert \vec v\rVert:\; módulo o norma de un vector
    k:  \left| k \right|:\;módulo o valor absoluto de un número real
    14

    La definición de producto de un escalar por un vector permite enunciar una condición para que dos vectores (no nulos) sean paralelos:

    vw          v  =  k.w      con  kR\vec v\parallel \vec w\;\; \Leftrightarrow \;\;\;\vec v\; = \;k.\vec w\;\;\;con\;k \in \mathbb{R}

    Ejemplos resueltos

    1.1.3   Propiedades de la suma de vectores y del producto por un escalar

    Sean u,v,wR3    y    α,  βR\vec u,\vec v,\vec w \in {\mathbb{R}^3}\;\;y\;\;\alpha ,\;\beta \in \mathbb{R}.

    Vimos que: u+vR3\vec u + \vec v \in {\mathbb{R}^3} y αuR3\alpha\, \vec u \in {\mathbb{R}^3}. Estas operaciones verifican las siguientes propiedades:

    1. u+v=v+u        \vec u + \vec v = \vec v + \vec u\;\;\;\;
    2. (u+v)+w=u+(v+w)        \left( {\vec u + \vec v} \right) + \vec w = \vec u + \left( {\vec v + \vec w} \right)\;\;\;\;
    3. u+0=0+u=u          \vec u + \vec 0 = \vec 0 + \vec u = \vec u\;\;\;\;\;
    4. u+(u)=(u)+u=0\vec u + \left( { - \vec u} \right) = \left( { - \vec u} \right) + \vec u = \vec 0
    5. α(u+v)=αu+αv\alpha \left( {\vec u + \vec v} \right) = \alpha \vec u + \alpha \vec v
    6. (α+β)u=αu+βu\left( {\alpha + \beta } \right)\vec u = \alpha \vec u + \beta \vec u
    7. α(βu)=(αβ)u\alpha \left( {\beta \vec u} \right) = \left( {\alpha \beta } \right)\vec u
    8. 1u=u1\vec u = \vec u
    15

    1.1.4   Módulo o norma de un vector en R3{\mathbb{R}^3}

    Nos interesa hallar una fórmula para calcular el módulo o norma de un vector. En R3{\mathbb{R}^3} el módulo es la longitud del vector. Para deducirla usaremos los triángulos rectángulos que quedan determinados tal como se muestra en el siguiente interactivo:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian

    En general escribiremos: v2=vx2+vy2+vz2{\lVert \vec v\rVert^2} = v_x^2 + v_y^2 + v_z^2

    Aplicando raíz cuadrada a ambos miembros:

    v=vx2+vy2+vz2\lVert \vec v\rVert = \sqrt {v_x^2 + v_y^2 + v_z^2}
    16

    1.1.5 Propiedades del módulo o norma de un vector

    1.  v0            v=0    v=0\ \lVert \vec v\rVert \ge 0\;\;\; \wedge \;\;\;\lVert \vec v\rVert = 0 \; \Leftrightarrow \;\vec v = \vec 0
    2.   k.v=k  v,    kR  \;\lVert k.\vec v \rVert = \left| k \right|\;\lVert \vec v\rVert,\;\;k \in \mathbb{R}\;
    3. Desigualdad triangular: v+wv+w\lVert \vec v + \vec w\rVert \le \lVert \vec v \rVert + \lVert \vec w \rVert

    ¿Qué condiciones tienen que cumplir los vectores v\vec v y w\vec w para que se verifique la igualdad: v+w=v+w  \lVert \vec v + \vec w \rVert = \lVert \vec v \rVert + \lVert \vec w\rVert \;?

      Ejemplo

    Sean v=(1,1,2)\vec v = \left( { - 1,1,2} \right) y w=(3,0,4)  \vec w = \left( {3,0, - 4} \right)\; calcular:

    a) v\lVert \vec v\rVert

    b) 2v\lVert - 2\vec v\rVert

    c) w\lVert \vec w\rVert

    d) v+w\lVert \vec v + \vec w\rVert

    Resolución

    v=(1)2+12+22=6\lVert \vec v\rVert = \sqrt {{{\left( { - 1} \right)}^2} + {1^2} + {2^2}} = \sqrt 6

    2v =(2)2+22+42=24=26 \lVert - 2\vec v\rVert\ = \sqrt {{{\left( { - 2} \right)}^2} + {2^2} + {4^2}} = \sqrt {24} = 2\sqrt 6

    w =32+02+(4)2=5\lVert \vec w\rVert\ = \sqrt {{3^2} + {0^2} + {{\left( { - 4} \right)}^2}} = 5

    17

    v+w=(2,1,2)=22+12+(2)2=3\lVert \vec v + \vec w\rVert = \lVert \left( {2,1, - 2} \right)\rVert = \sqrt {{2^2} + {1^2} + {{\left( { - 2} \right)}^2}} = 3

    Observemos que v+wv+w \lVert \vec v + \vec w\rVert \ne \lVert \vec v\rVert + \lVert \vec w\rVert


    1.1.6   Vector determinado por dos puntos

    Dados los puntos A(XA,YA,ZA)A\left( {{X_A},{Y_A},{Z_A}} \right) y B(XB,YB,ZB)B\left( {{X_B},{Y_B},{Z_B}} \right), el vector AB\overrightarrow {AB} , con origen en AA y extremo en BB, puede obtenerse cómo sigue:

    OA+AB=OB\overrightarrow {OA} + \overrightarrow {AB} = \overrightarrow {OB}

    AB=OBOA=(XB,YB,Z,B);(XA,YA,ZA) \Rightarrow \overrightarrow {AB} = \overrightarrow {OB} - \overrightarrow {OA} = \left( {{X_B},{Y_B},{Z_{,B}}} \right) -; \left( {{X_A},{Y_A},{Z_A}} \right)

    AB=(XBXA,YBYA,ZBZA) \Rightarrow \overrightarrow {AB} = \left( {{X_B} - {X_A},{Y_B} - {Y_A},{Z_B} - {Z_A}} \right)

    En este punto, te proponemos que nuevamente trabajes con el interactivo de la página 16.

    18

    1.1.7   Distancia entre dos puntos

    Problema

    ¿Cómo podríamos calcular la distancia entre R(1,1,4)R\left( {1,1,4} \right) y S(3,0,2)S\left( {3,0,2} \right)?

    Para hallar esta distancia armamos el vector RS\overrightarrow {RS} (o el SR)\overrightarrow {SR}) y calculamos su módulo:

    RS=(2,1,2)\overrightarrow {RS} = \left( {2, - 1, - 2} \right)

    RS=22+(1)2+(2)2=4+1+4=9=3\lVert \overrightarrow {RS}\rVert = \sqrt {{2^2} + {{\left( { - 1} \right)}^2} + {{\left( { - 2} \right)}^2}} = \sqrt {4 + 1 + 4} = \sqrt 9 = 3

    d(R,S)=3 \Rightarrow d\left( {R,S} \right) = 3

    En general

    Dados dos puntos A(xA,yA,zA)A\left( {{x_A},{y_A},{z_A}} \right) y B(xB,yB,zB)B\left( {{x_B},{y_B},{z_B}} \right) la distancia entre los mismos se calcula:

    d(A,B)=AB=(xBxA)2+(yByA)2+(zBzA)2d\left( {A,B} \right) = \lVert \overrightarrow {AB}\rVert = \sqrt {{{\left( {{x_B} - {x_A}} \right)}^2} + {{\left( {{y_B} - {y_A}} \right)}^2} + {{\left( {{z_B} - {z_A}} \right)}^2}}
    Problema

    Encontrar, si es posible, todos los puntos del eje zz cuya distancia al punto A(3,2,1)A\left( {3,2,1} \right) es 55.

    19

    Es recomendable hacer una figura de análisis del problema:

    Un punto del eje zz tiene la forma P(0,0,z)P (0,0,z).        Construyamos el vector desde un punto genérico cualquiera del eje zz hasta AA.

    PA=(3,  2,  1z)\overrightarrow {PA} = \left( {3,\;2,\;1 - z} \right)

    Se pide que el módulo (o norma) de PA\overrightarrow {PA} sea 5, entonces:

    PA=32+22+(1z)2=13+12z+z2= \overrightarrow {PA} = \sqrt {{3^2} + {2^2} + {{\left( {1 - z} \right)}^2}} = \sqrt {13 + 1 - 2z + {z^2}} = =142z+z2=5 = \sqrt {14 - 2z + {z^2}} = 5

    25=142z+z2    z22z11=0  25 = 14 - 2z + {z^2}\;\; \Rightarrow {z^2} - 2z - 11 = 0\;

    z=2+44.1.(11)2                z=244.1.(11)2      z = \frac{{2 + \sqrt {4 - 4.1.\left( { - 11} \right)} }}{2}\;\;\;\; \vee \;\;\;\;z = \frac{{2 - \sqrt {4 - 4.1.\left( { - 11} \right)} }}{2}\;\;\;

    z=2+482                z=2482      z = \frac{{2 + \sqrt {48} }}{2}\;\;\;\; \vee \;\;\;\;z = \frac{{2 - \sqrt {48} }}{2}\;\;\;

    z=1+124,46                z=1122,46    z = 1 + \sqrt {12} \cong 4,46\;\;\;\; \vee \;\;\;\;z = 1 - \sqrt {12} \cong - 2,46\;\;

    20

    Hemos llegado a que zz puede tomar dos valores distintos. Entonces existen dos puntos del eje zz cuya distancia al punto A(3,2,1)A\left( {3,2,1} \right) es 5. Son:

    P1(0,  0,  1+12)            P2(0,  0,  112){P_1}\left( {0,\;0,\;1 + \sqrt {12} } \right)\;\;\; \wedge \;\;\;{P_2}\left( {0,\;0,\;1 - \sqrt {12} } \right)

    1.1.8   Expresión canónica de un vector

    Recordemos que todo vector de R2{\mathbb{R}^2} puede expresarse como combinación lineal de los versores canónicos iˇ=(1,0) \check i = \left( {1,0} \right) y jˇ=(0,1) \check j = \left( {0,1} \right).

    v=(x,y) \vec v = \left( {x,y} \right) v=x(1,0)+y(0,1) \vec v = x\left( {1,0} \right) + y\left( {0,1} \right) v=x.iˇ+y.jˇ \vec v = x.\check i + y.\check j

    En forma análoga, todo vector de R3{\mathbb{R}^3} puede expresarse como combinación lineal de los versores canónicos:

    iˇ=(1,0,0)\check i = \left( {1,0,0} \right)

    jˇ=(0,1,0)\check j = \left( {0,1,0} \right)

    kˇ=(0,0,1)\check k = \left( {0,0,1} \right)

    21

    v=(x,y,z)=x(1,0,0)+y(0,1,0)+z(0,0,1)\vec v = \left( {x,y,z} \right) = x\left( {1,0,0} \right) + y\left( {0,1,0} \right) + z\left( {0,0,1} \right)

    v=x.iˇ+y.jˇ+z.kˇ \large \vec v = x.\check {i} + y.\check {j} + z.\check {k}\: (expresión canónica de v\vec v \,)

    1.1.9   Ángulos directores y cosenos directores de un vector

    Se denominan ángulos directores de un vector a los ángulos determinados por el vector y cada uno de los semiejes positivos, como se muestra en la siguiente figura:

    Los cosenos de dichos ángulos se llaman cosenos directores del vector. Aplicando relaciones trigonométricas, podemos obtener los cosenos directores:

    cos(α)=vxv    ,      cos(β)=vyv    ,    cos(γ)=vzv\cos \left( \alpha \right) = \frac{{{v_x}}}{{\lVert \vec v\rVert}}\;\;,\;\;\;cos\left( \beta \right) = \frac{{{v_y}}}{{\lVert \vec v\rVert}}\;\;,\;\;\cos \left( \gamma \right) = \frac{{{v_z}}}{{\lVert \vec v\rVert}}

    22

    Por lo tanto, los ángulos directores son:

    α=arcos(vxv),    β=arcos(vyv),    γ=arcos(vzv)\alpha = arcos\left( {\frac{{{v_x}}}{{\lVert \vec v\rVert}}} \right),\;\;\beta = arcos\left( {\frac{{{v_y}}}{{\lVert \vec v\rVert}}} \right),\;\;\gamma = arcos\left( {\frac{{{v_z}}}{{\lVert \vec v\rVert}}} \right)

    Donde α,β,γ\alpha ,\beta ,\gamma están comprendidos entre 00 y π\pi .

    Propiedad

    cos2(α)+cos2(β)+cos2(γ)=1{\cos ^2}\left( \alpha \right) + {\cos ^2}\left( \beta \right) + {\cos ^2}\left( \gamma \right) = 1

      Ejemplo

    Hallar los ángulos directores de v=(2,0,2) \vec v = \left( {2,0, -2} \right)

    Resolución

    Hallemos el módulo del vector:

    v=22+02+(2)2=8=22\lVert \vec v\rVert = \sqrt {{2^2} + {0^2} + {{\left( { - 2} \right)}^2}} = \sqrt 8 = 2\sqrt 2

    Ahora calculamos los ángulos con el arco coseno de los cocientes:

    α=arccos(222)=45\alpha = \arccos \left( {\frac{2}{{2\sqrt 2 }}} \right) = 45^\circ

    β=arccos(022)=90\beta = \arccos \left( {\frac{0}{{2\sqrt 2 }}} \right) = 90^\circ

    23

    γ=arccos(222)=135\gamma = \arccos \left( {\frac{{ - 2}}{{2\sqrt 2 }}} \right) = 135^\circ

    Veamos una gráfica del vector y sus ángulos directores:

    24

    1.1.10   Versor asociado a un vector

    Dado un vector no nulo v\vec v, se denomina versor asociado al vector unitario (de módulo 1) que tiene igual dirección y sentido que v\vec v.

    Dado v\vec v distinto de 0\vec 0, su versor asociado se obtiene dividiendo a este por su norma:

    v˘=vv \breve {v} = \frac {\vec v} {\lVert \vec v \rVert}

    Tomando en cuenta los cosenos directores,

    v=(vcosα,    vcosβ,    vcosγ)\vec v = \left( {\lVert \vec v\rVert\cos \alpha,\;\;\lVert \vec v\rVert\cos \beta,\;\;\lVert \vec v\rVert\cos \gamma} \right)

    Entonces

    vv=v˘=(cosα,    cosβ,    cosγ) \frac {\vec v} {\lVert \vec v \rVert} = \breve {v} = (\cos \alpha,\;\; \cos \beta,\;\; \cos \gamma)

    Es decir, las componentes del versor v˘\breve v son los cosenos directores de v \vec v .

    25

    1.2   Producto escalar en R3{\mathbb{R}^3}

    Sean u  ,  vR3\vec u\;,\;\vec v \in {\mathbb{R}^3}, y θ\theta el ángulo entre u\vec u y v\vec v, entonces el producto escalar entre u\vec u y v\vec v se define como sigue:

    u.v={u  vcos(θ)    si    u0  v00        si    u=0      v=0                  [2]\large \vec u.\vec v = \begin{cases} {\lVert \vec u\rVert\;\lVert \vec v\rVert\cos \left( \theta \right)\;\;si\;\;\vec u \ne \vec 0 \wedge \;\vec v \ne \vec 0}\\ {0\;\;\;\;si\;\;\vec u = \vec 0\;\; \vee \;\vec v = \vec 0} \end{cases}\;\;\;\;\;\;\;\;\;\left[ 2 \right]

    1.2.1   Propiedades del producto escalar

    1)  u.v  =  v.u\vec u.\vec v\; = \;\vec v.\vec u

    2)  u.(v  +  w)  =  u.v  +  u.w\vec u.\left( {\vec v\; + \;\vec w} \right)\; = \;\vec u.\vec v\; + \;\vec u.\vec w

    3)  k(u.v)=  (ku.v)=  u.(kv),  kRk\left( {\vec u.\vec v} \right) = \;\left( {k\vec u.\vec v} \right) = \;\vec u.\left( {k\vec v} \right),\;k \in \mathbb{R}

    26

    4)   v.v=(vx,vy,vz).  (vx,vy,vz)=vx2+vy2+vz2=v2\:\: \vec v.\vec v = \left( {{v_x},{v_y},{v_z}} \right).\;\left( {{v_x},{v_y},{v_z}} \right) = v_x^2 + v_y^2 + v_z^2 = {\lVert \vec v\rVert^2}

    v2>0  v0 {\lVert \vec v\rVert^2} > 0\;\forall \vec v \ne \vec 0

    De (4) se deduce que: v=v.v\lVert \vec v\rVert = \sqrt {\vec v.\vec v}

      Ejemplo

    Hallar u.v\vec u.\vec v para u=(0,0,1),  v=(0,3,3)\vec u = \left( {0,0,1} \right),\;\vec v = \left( {0,3, - 3} \right)

    Resolución

    Hagamos una gráfica para visualizar el ángulo entre los dos vectores:

    27

    Calculemos los módulos de u\vec u. y v\vec v:

    u=1\lVert \vec u\rVert = 1

    v=32+(3)2=18=32    \lVert \vec v\rVert = \sqrt {{3^2} + {{\left( { - 3} \right)}^2}} = \sqrt {18} = 3\sqrt 2 \;\;

    A partir del gráfico podemos determinar que el ángulo entre los vectores es   θ=135\;\theta = 135^\circ , o en radianes:   θ=34π\;\theta = \frac{3}{4}\pi .

    Calculemos el producto escalar:

    u.v=1.18.cos(135)=32.(22)=3\vec u.\vec v = 1.\sqrt {18} .\cos \left( {135^\circ } \right) = 3\sqrt 2 .\left( { - \frac{{\sqrt 2 }}{2}} \right) = - 3


    Pero no siempre es tan sencillo. Consideremos los vectores:

    u=(3,5,8)  ,            v  =(1,1,1)\vec u = \left( { - 3,5,8} \right)\;,\;\;\;\;\;\;\vec v\; = \left( {1,1,1} \right)

    Si quisiéramos calcular el producto escalar entre u\vec u y v\vec v, deberíamos conocer el ángulo comprendido entre dichos vectores.

    Usando el teorema del coseno se puede deducir otra fórmula para calcular el producto escalar en función de las componentes de los vectores.

    Sean u=(ux,uy,uz),    v=(vx,vy,vz)R3,  \vec u = \left( {{u_x},{u_y},{u_z}} \right),\;\;\vec v = \left( {{v_x},{v_y},{v_z}} \right) \in {\mathbb{R}^3},\; entonces:

    u.v=uxvx+uyvy+uzvz                [3] \boxed{\vec u.\vec v = {u_x}{v_x} + {u_y}{v_y} + {u_z}{v_z}\;}\;\;\;\;\;\;\;\left[ 3 \right]
    28

    Para los vectores dados, resulta:

    u  .v  =  (3).1  +  5  .1  +  8.1  =  10\vec u\;.\vec v\; = \;\left( { -3} \right).1\; + \;5\;.1\; + \;8.1\; = \;10

    1.2.2   Ángulo entre vectores

    Dados u,v\vec u,\vec v vectores no nulos de R3{\mathbb{R}^3}, queremos hallar el ángulo entre ellos.

    Si θ\theta es el ángulo entre u    y    v\vec u\;\;y\;\;\vec v, de las definiciones [2] y [3] de producto escalar resulta:

    cos(θ)=(u.vuv)\cos \left( \theta \right) = \left( {\frac{{\vec u.\vec v}}{{\lVert \vec u\rVert \lVert \vec v\rVert}}} \right)

    θ=  arccos(uxvx+uyvy+uzvzuv),      0θπ{\rm{\theta }} = {\rm{\;arccos}}\left( {\frac{{{u_x}{v_x} + {u_y}{v_y} + {u_z}{v_z}}}{{\lVert \vec u\rVert \lVert \vec v\rVert}}} \right),\;\;\;0 \le \theta \le \pi

    Por ejemplo, si u=(1,1,3)\vec u = \left( {1,1,3} \right)

    v=(1,0,4)  \vec v = \left( { - 1,0,4} \right)\;

    θ=arccos(1.(1)+1.0+3.412+12+32.(1)2+02+42    )\theta = \arccos \left( {\frac{{1.\left( { - 1} \right) + 1.0 + 3.4}}{{\sqrt {{1^2} + {1^2} + {3^2}} .\sqrt {{{\left( { - 1} \right)}^2} + {0^2} + {4^2}} \;\;}}} \right)

    θ=arccos(1111.17    )36,44\theta = \arccos \left( {\frac{{11}}{{\sqrt {11} .\sqrt {17} \;\;}}} \right) \cong 36,44^\circ

    29
    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian

    1.2.3   Condición de perpendicularidad

    Sean u,  v\vec u,\;\vec v no nulos,

    u.v=0  cos(θ)=0θ=π2  \vec u.\vec v = 0 \Leftrightarrow \;\cos \left( \theta \right) = 0 \Leftrightarrow \theta = \frac{\pi }{2}\;

    Esto permite enunciar una condición de perpendicularidad:

    u    vu.v=0\large \boxed{\vec u \; \bot \; \vec v \: \Leftrightarrow \: \vec u.\vec v = 0}

    Ejercicio para el lector

    Dados u=(1,2,3)      \vec u = \left( {1,2,3} \right)\;\;\;y     v=(0,2,5)    \;\;\vec v = \left( {0,2,5} \right)\;\;encontrar todos los vectores perpendiculares a u\vec u y a v\vec v de módulo 3.

    30

    1.2.4 Proyección de un vector en la dirección de otro

    El producto escalar es útil en problemas en los que se tiene interés en descomponer un vector como suma de vectores perpendiculares.

    Dados dos vectores no nulos u\vec u y v\vec v , nos proponemos descomponer u\vec u como suma de un vector paralelo a v    \vec v\;\;y otro perpendicular a v.  \vec v.\; O sea:

    u=u1+u2,u1v  \vec u = \vec u_1 + \vec u_2,\: \vec u_1 \parallel \vec v\; y   u2    v\;\vec u_2 \; \bot \; \vec v


    u1v      u1=kv,  kR          u=kv+u2\vec u_1 \parallel \vec v\; \Leftrightarrow \;\;\vec u_1 = k\vec v,\;k \in \mathbb{R}\;\;\; \Rightarrow \;\;\vec u = k\vec v + \vec u_2

    Podemos aplicar a ambos miembros producto escalar por v\vec v. Teniendo en cuenta que u2  .  v=0\vec u_2\; .\; \vec v = 0 por ser perpendiculares, resulta:

    31
    uv=(kv+u2)v    uv=kvvv2+u2v0  \vec u \cdotp \vec v =(k\, \vec v + \vec u_2)\, \cdotp \vec v \; \Rightarrow \;\vec u \,\cdotp \vec v = k \, \underbrace {\vec v \,\cdotp \vec v}_{{\lVert \vec v \rVert}^2} + \underbrace {\vec u_2 \,\cdotp \vec v}_{0} \; \Rightarrow   uv=kv2  0  k=uvv2 \Rightarrow \;\vec u \,\cdotp \vec v = k \, \underbrace {{\lVert \vec v \rVert}^2}_{\neq \; 0} \; \Rightarrow k =\frac {\vec u \cdotp \vec v} {{\lVert \vec v \rVert}^2}

    Entonces:

    u1=(uvv2)v\vec u_1 = \left( {\frac{{\vec u \, \cdotp \vec v}}{{{\lVert \vec v\rVert^2}}}} \right)\vec v

    Este vector es la proyección de  u  \;\vec u\; en la dirección de   v\;\vec v :

      proyv(u)=(uvv2)v\;{\overrightarrow {proy} _{\vec v}}\left( {\vec u} \right) = \left( {\frac{{\vec u \, \cdotp \vec v}}{{{\lVert \vec v\rVert^2}}}} \right)\vec v

    El vector u2  {\vec u_2}\; puede obtenerse por diferencia:

    u=u1+u2      u2=uu1\vec u = \vec u_1 + \vec u_2 \;\; \Rightarrow \;\vec u_2 = \vec u - \vec u_1

    Recordemos que u2\vec u_2 debe ser perpendicular a   v\;\vec v.

    Para resolver algunos problemas geométricos, es útil calcular el módulo del vector de proyección:

    proyv(u)=uvv2 v=uvv2  v  =  uvv\lVert {\overrightarrow {proy} _{\vec v}}\left( {\vec u} \right)\rVert = \left| {\frac{{\vec u \, \cdotp \vec v}}{{{\lVert \vec v \rVert^2}}}} \right|\ \lVert \vec v \rVert = \frac{{\left| {\vec u \, \cdotp \vec v} \right|}}{{{\lVert \vec v \rVert^2}}}\;\lVert \vec v \rVert\; = \;\frac{{\left| {\vec u \, \cdotp \vec v} \right|}}{{\lVert \vec v \rVert}}

    32

    Ejemplo

    Descomponer u=(1,2,1)\vec u = \left( {1,2,1} \right) como suma de un vector paralelo a v=(0,1,1)\,\vec v = \left( {0,1, - 1} \right) más otro perpendicular a v\vec v.

    Primero buscamos u1\vec u_1 :

    u1=proyv(u)=(u.vv2)v= \vec u_1 = \overrightarrow {proy} _{\vec v} ( \vec u ) = \Bigg( \frac{\vec u.\vec v}{{\lVert \vec v \rVert^2}} \Bigg) \vec v = u1=((1,2,1).(0,1,1)02+12+(1)2)(0,1,1) \vec u_1 = \Bigg( \frac{(1,2,1).(0,1,-1)}{0^2+1^2+(-1)^2} \Bigg) (0,1,-1)

    =0+212  (0,1,1)=(0,  1  2,12) = \frac{{0 + 2 - 1}}{2}\;\left( {0,1, - 1} \right) = \left( {0,\;\frac{1}{{\;2}}, - \frac{1}{2}} \right)

    u2=uu1=(1,2,1)(0,  1  2,12)=(1,  32,  32)\vec u_2 = \vec u - \vec u_1 = \left( {1,2,1} \right) - \left( {0,\;\frac{1}{{\;2}}, - \frac{1}{2}} \right) = \left( {1,\;\frac{3}{2},\; \frac{3}{2}} \right)

    33

    Comprobación:     u2  v\;\;\vec u_2 \, \, \bot \, \vec v

    (1,  32,  32).(0,1,1)=0\left( {1,\;\frac{3}{2},\;\frac{3}{2}} \right).\left( {0,1, -1} \right) = 0


    La siguiente escena de Geogebra nos permite visualizar la proyección de un vector sobre otro en R3\mathbb{R}^3:

    Objeto interactivo de construcción propia
    34

    A continuación dejamos un ejercicio para que pongas a prueba tus conocimientos sobre producto escalar

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    35

    1.3   Producto vectorial y mixto

    1.3.1   Producto vectorial

    Definición

    Para resolver numerosos problemas de Geometría, Física e Ingeniería, interesa construir un vector en R3{\mathbb{R}^3} que sea perpendicular a dos vectores dados.

    O sea: dados u,vR3\vec u,\vec v \in {\mathbb{R}^3}, nos proponemos hallar un vector w\vec w tal que wu  \vec w \bot \vec u\; y     wv\;\overrightarrow {\;w} \bot \vec v.

    El producto vectorial es una operación entre vectores que facilita la obtención de w\vec w.

    Definición: El producto vectorial de u  y  v\vec u\;y\;\vec v, que indicaremos u×v\vec u \times \vec v es un vector

    u×v=w\vec u \times \vec v = \vec w

    Que tiene:

    • Dirección perpendicular a u\vec u y a v\vec v:     wu        wv\;\;\vec w \bot \vec u\;\; \wedge \;\;\vec w \bot \vec v
    • Sentido: regla de la mano derecha: si con la mano derecha se recorre el menor ángulo posible entre u\vec u y v\vec v, el pulgar indica el sentido de w\vec w.
    36
  • Módulo:
  • w=u  v  sen(θ)\lVert \vec w\rVert = \lVert \vec u\rVert\;\lVert \vec v\rVert\;sen\left( \theta \right)

    siendo θ  \theta \;el ángulo comprendido entre u    y    v\vec u\;\;y\;\;\vec v

    Se puede ver que no es una operación conmutativa porque si cambiamos el orden de los vectores, se conservan la dirección y el módulo del producto vectorial pero se invierte su sentido:

    u×v=(v×u)\vec u \times \vec v = - \left( {\vec v \times \vec u} \right) Ejemplo

    a=(0,1,1)\vec a = \left( {0,1, - 1} \right)

    b=(0,0,3)\vec b = \left( {0,0,3} \right)

    Hallar el vector a×b\vec a \times \vec b.

    Dirección: eje xx

    Sentido: semieje positivo de xx

    Módulo: 2.3.sen(135)=2.3.22=3    \sqrt 2 .3.sen\left( {135^\circ } \right) = \sqrt 2 .3.\frac{{\sqrt 2 }}{2} = 3\;\;

    Entonces:

    a×b=(3,0,0)\vec a \times \vec b = \left( {3,0,0} \right)

    37

    1.3.2   Propiedades del producto vectorial

    1. u×v=(v×u)\vec u \times \vec v = - \left( {\vec v \times \vec u} \right)
    2. (u+v)×w=u×w+v×w\left( {u + v} \right) \times w = \vec u \times \vec w + \vec v \times \vec w
    3. (ku)×v=k(u×v)=u  ×  (kv)\left( {k\vec u} \right) \times \vec v = k\left( {\vec u \times \vec v} \right) = \vec u\; \times \;\left( {k\vec v} \right)
    4. v×v=0  ,    pues      v×v=v  v  sen(0)=0\vec v \times \vec v = \vec 0\;,\;\;pues\;\;\;\lVert \vec v \times \vec v\rVert = \lVert \vec v\rVert\;\lVert \vec v\rVert\;sen\left( 0 \right) = 0
    5. uv    u×v=0  ,      pues      sen(0)=sen(180)=0\vec u\parallel \vec v\; \Rightarrow \;\vec u \times \vec v = \vec 0\;,\;\;\;pues\;\;\;sen\left( {0^\circ } \right) = sen\left( {180^\circ } \right) = 0
    6. 0×v  =  v×0  =  0\vec 0 \times \vec v\; = \;\vec v \times \vec 0\; = \;\vec 0

    Si u\vec u y   v\overrightarrow {\;v} son vectores no nulos, podemos enunciar una condición necesaria y suficiente de paralelismo:

    uv      u×v=0      u=kv      con    k  R\vec u\parallel \vec v\;\; \Leftrightarrow \;\vec u \times \vec v = \vec 0\;\; \Leftrightarrow \;\vec u = k\vec v\;\;\;con\;\;k\; \in {\mathbb{R}}

    1.3.3   Fórmula para calcular el producto vectorial

    Dados u=(ux,uy,uz)\vec u = \left( {{u_x},{u_y},{u_z}} \right) y v=(vx,vy,vz)\vec v = \left( {{v_x},{v_y},{v_z}} \right) podemos hallar una fórmula para el producto vectorial expresando los vectores en forma canónica.

    Aplicando propiedades del producto vectorial y considerando los productos entre los versores canónicos, se obtiene la siguiente fórmula:

    u×v=(uy  vzuz  vy,  (ux  vzuz  vx),    ux  vyuy  vx)   \vec u \times \vec v = ( {{u_{y\;}}{v_z} - {u_{z\;}}{v_y},\; - ( {{u_{x\;}}{v_z} - {u_{z\;}}{v_x}}),\;\;{u_{x\;}}{v_y} - {u_{y\;}}{v_x}})\; [1]
    38

    Esta fórmula puede expresarse en forma más sencilla utilizando determinantes, tema que presentaremos brevemente y luego desarrollaremos en la próxima unidad.

    Una matriz es un ordenamiento rectangular de números, como caso particular nos interesan las matrices cuadradas (igual número de filas y de columnas).

    A=(a11a12a21a22)matriz de 2x2 (2 filas y 2 columnas) A =\begin{pmatrix} a_{11} & a_{12} \\ a_{21} & a_{22} \end{pmatrix} \textsf {matriz de 2x2 (2 filas y 2 columnas)} A=(a11a12a13a21a22a23a31a32a33)matriz de 3x3 (3 filas y 3 columnas) A =\begin{pmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{pmatrix} \textsf {matriz de 3x3 (3 filas y 3 columnas)}

    A cada matriz cuadrada puede asignársele un número real que llamaremos su determinante y designaremos como det(A)det\left( A \right) o A\left| A \right| . Para matrices de 2×22 \times 2 y de 3×33 \times 3 el determinante se calcula como sigue:

    A = ( a 11 a 12 a 21 a 22 ) | A | = a 11 . a 22 a 12 . a 21

    A=(a11a12a13a21a22a23a31a32a33) A =\begin{pmatrix} \color{blue} {a_{11}} & \color{fuchsia} {a_{12}} & \color{red} {a_{13}} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{pmatrix} \: \Longrightarrow A=a11a22a23a32a33a12a21a23a31a33+a13a21a22a31a32 \Longrightarrow \: \begin{vmatrix} A \end{vmatrix} = \color{blue} {a_{11}} \color{black} \begin{vmatrix} a_{22} & a_{23} \\ a_{32} & a_{33} \end{vmatrix} - \color{fuchsia} {a_{12}} \color{black} \begin{vmatrix} a_{21} & a_{23} \\ a_{31} & a_{33} \end{vmatrix} + \color{red} {a_{13}} \color{black} \begin{vmatrix} a_{21} & a_{22} \\ a_{31} & a_{32} \end{vmatrix}
    39

    Podemos expresar la fórmula [1][1] utilizando determinantes como sigue:

    u×v=i˘j˘k˘uxuyuzvxvyvz= \vec u \times \vec v = \begin{vmatrix} \breve{i} & \breve{j} & \breve{k} \\ u_x & u_y & u_z \\ v_x & v_y & v_z \end{vmatrix} = u×v=uyuzvyvzi˘uxuzvxvzj˘+uxuyvxvyk˘ \vec u \times \vec v = \begin{vmatrix} u_y & u_z \\ v_y & v_z \end{vmatrix}\breve{i} - \begin{vmatrix} u_x & u_z \\ v_x & v_z \end{vmatrix}\breve{j} + \begin{vmatrix} u_x & u_y \\ v_x & v_y \end{vmatrix}\breve{k}

    Con la notación habitual de ternas, resulta:

    u×v=(uyuzvyvz,    uxuzvxvz,    uxuyvxvy) \vec u \times \vec v = \Bigg( \begin{vmatrix} u_y & u_z \\ v_y & v_z \end{vmatrix},\;\; - \begin{vmatrix} u_x & u_z \\ v_x & v_z \end{vmatrix},\;\; \begin{vmatrix} u_x & u_y \\ v_x & v_y \end{vmatrix} \Bigg)

    Veamos cómo utilizar esta regla práctica para calcular un producto vectorial:

    u=(1,2,3)\vec u = \left( {1,2,3} \right)

    v=(0,2,5)\vec v = \left( {0,2,5} \right)

    u×v=i˘j˘k˘123025=(2325,1305,  1202)=(4,5,2) \vec u \times \vec v = \begin{vmatrix} \breve{i} & \breve{j} & \breve{k} \\ 1 & 2 & 3 \\ 0 & 2 & 5 \end{vmatrix} = \Bigg( \begin{vmatrix} 2 & 3 \\ 2 & 5 \end{vmatrix}, - \begin{vmatrix} 1 & 3 \\ 0 & 5 \end{vmatrix},\; \begin{vmatrix} 1 & 2 \\ 0 & 2 \end{vmatrix} \Bigg)=(4,-5,2) Comprobemos que el vector obtenido es ortogonal a u\vec u y a v\vec v:

    (4,5,2).(1,2,3)=0\left( {4, - 5,2} \right).\left( {1,2,3} \right) = 0

    (4,5,2).(0,2,5)=0\left( {4, - 5,2} \right).\left( {0,2,5} \right) = 0

    40

    Te invitamos a que escribas los vectores del ejercicio anterior en el siguiente applet:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian

    1.3.4 Interpretación geométrica del módulo del producto vectorial

    Consideremos los siguientes vectores y calculemos el módulo del producto vectorial.

    u=(1,2,3),      v=(0,2,5)\vec u = \left( {1,2,3} \right),\;\;\;\vec v = \left( {0,2,5} \right)

    u×v=(4,5,2)\vec u \times \vec v = \left( {4, - 5,2} \right)

    u×v=45=3.5\lVert \vec u \times \vec v \rVert = \sqrt {45} = 3.\sqrt 5

    ¿Qué representa este número 353\sqrt 5 ?

    41

    Dibujemos cualquier par de vectores u\vec u y v\vec v y hagamos el paralelogramo determinado entre ellos:

    Aˊrea  del  paralelogramo=u.hÁrea\;del\;paralelogramo = \lVert \vec u \rVert.h

    sen(θ)=hv  h=vsen(θ)sen\left( \theta \right) = \frac{h}{{\lVert \vec v \rVert}}\; \Rightarrow h = \lVert \vec v \rVert sen\left( \theta \right)

    Aˊrea  del  paralelogramo=u  v  sen(θ)=u×v Área\;del\;paralelogramo = \lVert \vec u \rVert\;\lVert \vec v \rVert\;sen\left( \theta \right) = \lVert \vec u \times \vec v\ \rVert

    Conclusión: Dados dos vectores no paralelos, el módulo de su producto vectorial representa el área del paralelogramo determinado por dichos vectores.

    42
    Ejemplo

    Dados A(1,3,1)A\left( {1,3,1} \right), B(2,3,5)B\left( {2, - 3,5} \right) y C(0,2,1)C\left( {0,2,1} \right) calcular el área del ΔABC{\rm{\Delta }}ABC.

    Podemos pensar al área del triángulo como la mitad del área del paralelogramo:

    Aˊrea  del  triaˊngulo=12.AB×ACÁrea\;del\;triá ngulo = \frac{1}{2}. \lVert \overrightarrow {AB} \times \overrightarrow {AC} \rVert

    AB=(1,6,4)\overrightarrow {AB} = \left( {1, - 6,4} \right)

    AC=(1,1,0)\overrightarrow {AC} = \left( { - 1, - 1,0} \right)

    AB×AC=(4,4,7)\overrightarrow {AB} \times \overrightarrow {AC} = \left( {4, - 4, - 7} \right)

    43

    AB×AC=16+16+49=81=9\lVert \overrightarrow {AB} \times \overrightarrow {AC} \rVert = \sqrt {16 + 16 + 49} = \sqrt {81} = 9

    Aˊrea  del  triaˊngulo=12.9=4,5Área\;del\;triá ngulo = \frac{1}{2}.9 = 4,5

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    44

    Proponemos la siguiente ejercitación a manera de repaso del tema visto.

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    45

    1.3.5   Producto mixto

    Definición

    Dados tres vectores u,v,wR3\vec u,\vec v,\vec w \in {\mathbb{R}^3}, se denomina producto mixto al número real que se obtiene multiplicando u.(v×w)\vec u.\left( {\vec v \times \vec w} \right).

    Para obtener el producto mixto, se calcula primero el producto vectorial y luego el escalar.

    Dejamos como ejercicio para el lector, demostrar que:

    u.(v×w)=(u×v).w\vec u.\left( {\vec v \times \vec w} \right) = \left( {\vec u \times \vec v} \right).\vec w

    Veamos un ejemplo:

    u=(1,2,3)\vec u = \left( {1,2,3} \right)

    v=(0,2,5)\vec v = \left( {0,2,5} \right)

    w=(0,0,2)  \vec w = \left( {0,0,2} \right)\;

    v×w=i˘j˘k˘025002=(4,0,0) \vec v \times \vec w = \begin{vmatrix} \breve{i} & \breve{j} & \breve{k} \\ 0 & 2 & 5 \\ 0 & 0 & 2 \end{vmatrix} = (4,\,0,\,0)

    u.(v×w)=(1,2,3)(4,0,0)=4\vec u.\left( {\vec v \times \vec w} \right) = (1,\,2\,,3)\, \cdotp (4,\,0,\,0)= 4

    Ahora hagamos al revés:

    (u×v).w\left( {\vec u \times \vec v} \right).\vec w

    46
    u×v=i˘j˘k˘123025=(4,5,2) \vec u \times \vec v = \begin{vmatrix} \breve{i} & \breve{j} & \breve{k} \\ 1 & 2 & 3 \\ 0 & 2 & 5 \end{vmatrix} = (4,\,-5,\,2)

    (u×v).w=(4,5,2).(0,0,2)=4\left( {\vec u \times \vec v} \right).\vec w = \left( {4, - 5,2} \right).\left( {0,0,2} \right) = 4

    Existe otro procedimiento para calcular el producto mixto a través de un determinante de 3×3.

    Dados u=(u1,  u2,  u3),  v=(v1,  v2,v3),  w=(w1,  w2,  w3),\vec u = (u_1,\;u_2,\;u_3),\;\vec v = (v_1,\;v_2,v_3),\;\vec w = (w_1,\;w_2,\;w_3), el producto mixto es:

    u(v×w)=(u×v)w=u1u2u3v1v2v3w1w2w3= \vec u \, \cdotp (\vec v \times \vec w) = (\vec u \times \vec v) \, \cdotp \vec w = \begin{vmatrix} u_1 & u_2 & u_3 \\ v_1 & v_2 & v_3 \\ w_1 & w_2 & w_3 \end{vmatrix} = u(v×w)=u1v2v3w2w3u2v1v3w1w3+u3v1v2w1w2 \vec u \, \cdotp (\vec v \times \vec w) = u_1 \begin{vmatrix} v_2 & v_3 \\ w_2 & w_3 \end{vmatrix} - u_2 \begin{vmatrix} v_1 & v_3 \\ w_1 & w_3 \end{vmatrix} + u_3 \begin{vmatrix} v_1 & v_2 \\ w_1 & w_2 \end{vmatrix}

    Calculemos con este método el producto mixto entre:

    u=(1,2,3),  v=(0,2,5),  w=(0,0,2)\vec u = (1,2,3), \; \vec v = (0,2,5), \; \vec w = (0,0,2)

    u(v×w)=123025002=4 \vec u \, \cdotp (\vec v \times \vec w) = \begin{vmatrix} 1 & 2 & 3 \\ 0 & 2 & 5 \\ 0 & 0 & 2 \end{vmatrix} = 4

    ¿Qué significado tiene para la Geometría este número que hemos obtenido con el producto mixto?

    47

    1.3.6   Interpretación geométrica del producto mixto

    Consideremos 3 vectores de R3{\mathbb{R}^3} y construyamos un paralelepípedo (cuerpo cuyas seis caras son paralelogramos):

    Volumen  del  paralelepıˊpedo=aˊrea  de  la  base  .  AlturaVolumen\;del\;paralelepípedo = área\;de\;la\;base\;.\;Altura

    A=u×vA = \lVert \vec u \times \vec v \rVert

    h=cos(θ)  wh = \left| {\cos \left( \theta \right)} \right|\;\lVert \vec w\rVert

    siendo θ\theta el ángulo entre (u×v)    y    w(\vec u \times \vec v) \;\;y\;\;\vec w.

    Observación: cos(θ)\cos \left( \theta \right) podría ser negativo, por eso tomamos su valor absoluto para el cálculo de la altura.

    Entonces:

    V=u×v   w  cos(θ)            [1]V = \lVert \vec u \times \vec v\ \rVert \; \lVert \vec w \rVert \;\left| {\cos \left( \theta \right)} \right|\;\;\;\;\;\;\left[ 1 \right]

    48
    Por otra parte:

    (u×v).w=u×v  w  cos(aˊng(u×v  ,  w))=\left( {\vec u \times \vec v} \right).\vec w = \lVert \vec u \times \vec v \rVert \; \lVert \vec w \rVert \;\cos \left( {á ng\left( {\vec u \times \vec v\;,\;\vec w} \right)} \right) = =u×v  w  cos(θ)        [2] = \lVert \vec u \times \vec v \rVert\;\lVert \vec w \rVert \;\cos \left( \theta \right)\;\;\;\;\left[ 2 \right]

    De [1] y [2] resulta:

    Volumen  del  paralelepipedo=  (u×v).w   \boxed { \text {Volumen\;del\;paralelepipedo} = \left| {\;( {\vec u \times \vec v}).\vec w\;} \right|}

    Retomando el ejemplo con los vectores u=(1,2,3),\vec u = \left( {1,2,3} \right),
    v=(0,2,5),\vec v = \left( {0,2,5} \right), w=(0,0,2)\vec w = \left( {0,0,2} \right); de acuerdo al valor del producto mixto obtenido, el volumen del paralelepípedo determinado por los tres vectores es igual a 4.

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    49

    1.3.7   Coplanaridad

    Consideremos los vectores u=(1,0,3),v=(0,0,2),\vec u = (1,0,3), \vec v = (0,0,2),
    w=(3,0,4) \vec w = (3,0,4). Te proponemos que verifiques que el producto mixto da cero.

    Si el producto mixto es cero, el volumen es 0, o sea que no se forma el paralelepípedo.

    Veamos una gráfica de estos tres vectores:

    Observamos que los tres vectores están en el plano y=0y = 0, es decir que son coplanares.

    50

    Tres vectores u\vec u, v\vec v, w\vec w de R3{\mathbb{R}^3} se denominan coplanares si considerados con un origen común, sus direcciones quedan incluidas en un mismo plano.

    El producto mixto nos permite enunciar una condición de coplanaridad:

    u,v,w    \vec u,\vec v,\vec w\;\;son coplanares     \; \Leftrightarrow \; (u×v).w=0\left( {\vec u \times \vec v} \right).\vec w = 0

    Lo que hemos desarrollado hasta aquí sobre vectores resulta una herramienta potente para el estudio de la geometría de rectas y planos en R3{\mathbb{R}^3}, como veremos a continuación.

    51

    1.4   Ecuaciones del plano

    1.4.1   Deducción de la ecuación general del plano

    Dada una dirección en R3{\mathbb{R}^3}, existen infinitos planos perpendiculares a la misma. Si conocemos además un punto del plano, éste queda determinado de forma única.

    Nos proponemos hallar la ecuación del plano π\pi que pasa por P0(x0,y0,z0){P_0}\left( {{x_0},{y_0},{z_0}} \right) y es perpendicular al vector n=(a,b,c).\vec n = \left( {a,b,c} \right). El vector n\vec n se denomina vector normal del plano.

    ¿Qué condición debe cumplir un punto P(x,y,z)P\left( {x,y,z} \right) para estar en el plano π\pi ? Si armamos el vector   P0P  \overrightarrow {\;{P_0}P} \;, éste debe ser paralelo al plano, o sea perpendicular al vector normal del plano:

    P(x,y,z)π    P0P    n        P0P.n=0P\left( {x,y,z} \right) \in \pi \; \Leftrightarrow \overrightarrow {\;{P_0}P} \; \bot \; \vec n\;\; \Leftrightarrow \;\;\overrightarrow {{P_0}P} .\, \vec n = 0

    (xx0  ,  yy0,  zz0).(a,b,c)=0\left( {x - {x_{0\;}},\;y - {y_0},\;z - {z_0}} \right).\left( {a,b,c} \right) = 0

    a(xx0)+b(yy0)+c(zz0)=0a\left( {x - {x_0}} \right) + b\left( {y - {y_0}} \right) + c\left( {z - {z_0}} \right) = 0

    ax+by+cz+(ax0by0cz0)d=0 ax + by + cz + \underbrace {(-ax_0-by_0-cz_0)}_{d} = 0
    52

    Ecuación general o implícita del plano

    ax+by+cz+d=0 \boxed {ax + by + cz + d = 0}

    1.4.2   Ecuación segmentaria del plano

    Dada la ecuación general de un plano:

    π:    ax+by+cz+d=0\pi :\;\;ax + by + cz + d = 0

    Si a,b,c,d  a,b,c,d\; son distintos de cero, es posible obtener otra ecuación del plano como sigue:

    ax+by+cz=dax + by + cz = - d

    adx+bdy+cdz=1\frac{a}{{ - d}}\,x + \frac{b}{{ - d}}\,y + \frac{c}{{ - d}}\,z = 1

    53

    x(da)+y(db)+z(dc)=1 \frac{x}{\bigg(- \dfrac{d}{a}\bigg)} + \frac{y}{\bigg(- \dfrac{d}{b}\bigg)} + \frac{z}{\bigg(-\dfrac{d}{c}\bigg)} = 1

    Si llamamos p=da,  q=db,  r=dc p = - \dfrac{d}{a},\; q = - \dfrac{d}{b},\; r = - \dfrac{d}{c}

    Resulta:

    Ecuación segmentaria del plano

    xp+yq+zr=1 \boxed {\frac{x}{p} + \frac{y}{q} + \frac{z}{r} = 1}

    Veamos qué indican p,  \,p,\; q,  q,\;y r\,r\,:

    ¿Cuál es la intersección del plano con el eje x?

    y=0,  z=0  x=p     y=0,\;z=0 \; \rArr x = p \; \rArr \; El plano corta al eje xx en (p,0,0).(p,0,0).

    ¿Cuál es la intersección con el eje y?

    (0,q,0)\left( {0,q,0} \right)

    ¿Y con el eje z?

    (0,0,r)\left( {0,0,r} \right)

    Podemos observar que p, q y r indican las intersecciones con los ejes.

    54

    Ejemplo

    2x3y+z6=02x - 3y + z - 6 = 0

    2x3y+z=62x - 3y + z = 6

    2x63y6+z6=1\frac{{2x}}{6} - \frac{{3y}}{6} + \frac{z}{6} = 1

    x3y2+z6=1\frac{x}{3} - \frac{y}{2} + \frac{z}{6} = 1

    Esta ecuación parece segmentaria pero no lo es por el signo negativo.  La reescribimos así:

    x3+y2+z6=1              Ecuacioˊn  segmentaria\frac{x}{3} + \frac{y}{{ - 2}} + \frac{z}{6} = 1\;\;\;\;\;\;\;Ecuación\;segmentaria

    La ecuación segmentaria es práctica para graficar un plano porque muestra los tres puntos de corte con los ejes:

    55

    1.4.3   Ecuación vectorial paramétrica del plano

    Dados dos vectores u=(u1,u2,u3)\vec u = \left( {{u_1},{u_2},{u_3}} \right) y v  =(v1,  v2,  v3)\vec {v\;} = \left( {{v_1},\;{v_2},\;{v_3}} \right) no paralelos y un punto P0  (x0,  y0,  z0){P_0}\;\left( {{x_0},\;{y_0},\;{z_0}} \right), nos proponemos hallar la ecuación del plano π\pi que pasa por P0{P_0} y es paralelo a u\vec u y v\vec v .

    ¿Cómo podemos obtener un vector perpendicular al plano conociendo dos vectores paralelos a dicho plano?

    n=u×v\vec n = \vec u \times \vec v

    Teniendo n  \vec n\; y el punto P0{P_0}, podemos hallar la ecuación implícita o general del plano π\pi como habíamos visto previamente.

    Obtendremos a continuación otro tipo de ecuación del plano, cuya deducción se basa en el concepto de combinación lineal de vectores.

    Si P(x,y,z)P\left( {x,y,z} \right) es un punto cualquiera del plano π\pi , los vectores P0P  ,  u    y      v  \overrightarrow {{P_0}P} \;,\;\vec u\;\;y\;\;\;\vec v\; son coplanares

    56

    Entonces
    α,βR      P0P=α  u+β  v\exists \alpha ,\beta \in \mathbb{R}\;|\;\;\overrightarrow {{P_0}P} = \alpha \;\vec u + \beta \;\vec v

    Esto significa que el vector P0P\overrightarrow {{P_0}P} puede expresarse como combinación lineal de u\vec u y v\vec v, como se muestra en la figura:


    (xx0  ,  yy0,zz0)=α.(u1,u2,u3)  +β  (v1,v2,v3)\left( {x - {x_0}\;,\;y - {y_0},z - {z_0}} \right) = \alpha .\left( {{u_1},{u_2},{u_3}} \right)\; + \beta \;\left( {{v_1},{v_2},{v_3}} \right)

    Por lo tanto:

    (x,y,z)=(x0,y0,z0)+α  (u1,u2,u3)+β  (v1,v2,v3)\left( {x,y,z} \right) = \left( {{x_0},{y_0},{z_0}} \right) + \alpha \;\left( {{u_1},{u_2},{u_3}} \right) + \beta \;\left( {{v_1},{v_2},{v_3}} \right),    con  α,βR\;\;con\;\alpha ,\beta \in \mathbb{R}
    O en notación vectorial:

    Ecuación vectorial paramétrica del plano

    (x,y,z)=OP0+α.u+β.v   \boxed { \left( {x,y,z} \right) = \overrightarrow {O{P_0}} + \alpha .\vec u + \beta .\vec v\; }
    57

    Ejemplo

    Armar la ecuación vectorial paramétrica del plano paralelo a u=(3,1,5)\vec u = \left( {3, - 1,5} \right) y v=(7,3,2)\vec v = \left( {7,3,2} \right) que pasa por el punto P0(0,1,8).{P_0}\left( {0, - 1,8} \right).

    De acuerdo con lo que hemos visto, tenemos toda la información para escribir la ecuación vectorial paramétrica:

    (x,y,z)=(0,1,8)+α(3,1,5)+β(7,3,2);  α,βR\left( {x,y,z} \right) = \left( {0, - 1,8} \right) + \alpha \left( {3, - 1,5} \right) + \beta \left( {7,3,2} \right);\;\alpha ,\beta \in \mathbb{R}

    Nota: Para cada α  y  βR\alpha \;y\;\beta \in \mathbb{R} se obtiene un punto del plano. Por ejemplo si   α=1    y      β=1\;\alpha = 1\;\;y\;\;\;\beta = - 1 se obtiene el punto (x,y,z)=(4,  5,  11)\left( {x,y,z} \right) = \left( { - 4,\; - 5,\;11} \right).

    Busquemos ahora la ecuación general de este plano.

    n=u×v=(3,1,5)×(7,3,2)=(17,  29,16)\vec n = \vec u \times \vec v = \left( {3, - 1,5} \right) \times \left( {7,3,2} \right) = \left( { - 17,\;29,16} \right)

    Luego:         17x+29y+16z+d=0 - 17x + 29y + 16z + d = 0

    Reemplazamos P0{P_0} para obtener dd:

    17.0+29.(1)+16.8+d=0  d=99 - 17.0 + 29.\left( { - 1} \right) + 16.8 + d = 0 \Rightarrow \;d = - 99

    Luego:        17x+29y+16z99=0 - 17x + 29y + 16z - 99 = 0

    que es la ecuación general o implícita del plano.

    58

    1.4.4   De la ecuación general a la ecuación vectorial paramétrica

    Dada la ecuación general de un plano, ¿cómo puede obtenerse una ecuación vectorial paramétrica de dicho plano?

    Consideremos el siguiente ejemplo:

    ω:      2xy+3z+9=0\omega :\;\;\;2x - y + 3z + 9 = 0

    Podemos despejar cualquiera de las variables, por ejemplo yy:

    y=2x+3z+9y = 2x + 3z + 9

    Entonces: ω:      (x,y,z)=(x,  2x+3z+9,  z)\omega :\;\;\;\left( {x,y,z} \right) = \left( {x,\;2x + 3z + 9,\;z} \right)

    Reescribimos como suma de tres vectores, de forma tal que uno de ellos tenga los términos con xx, otro los términos con zz y otro los términos independientes:

    (x,y,z)=(x,2x,0)+(0,3z,z)+(0,9,0)  \left( {x,y,z} \right) = \left( {x,2x,0} \right) + \left( {0,3z,z} \right) + \left( {0,9,0} \right)\;

    (x,y,z)=x(1,2,0)+z(0,3,1)+(0,9,0)  ,  con    x,zR  \left( {x,y,z} \right) = x\left( {1,2,0} \right) + z\left( {0,3,1} \right) + \left( {0,9,0} \right)\;,\;con\;\;x,z \in \mathbb{R}\;

    Si llamamos x=αx = \alpha , z=β,  z = \beta,\; resulta:

    ω:  (x,y,z)=(0,9,0)+α(1,2,0)+β(0,3,1),  con  α,βR\omega:\;\left( {x,y,z} \right) = \left( {0,9,0} \right) + \alpha \left( {1,2,0} \right) + \beta \left( {0,3,1} \right),\;con\;\alpha ,\beta \in \mathbb{R}

    59

    Obtuvimos así una ecuación vectorial paramétrica del plano ω.\omega.

    El lector puede comprobar que:

    i)   los vectores u=(1,2,0)\vec u = (1,2,0)\, y v=(0,3,1)  \,\vec v = (0,3,1)\; son perpendiculares a n=(2,1,3),  \,\vec n = (2,-1,3),\; o sea que son paralelos al plano;

    ii)  P0(0,9,0)ω.P_0 \: (0,9,0) \in \omega.

    Ecuación vectorial del plano. (Applet de Geogebra de elaboración propia.)
    60

    Interactivo que permite visualizar la deducción de la página 52:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian

    Para practicar:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    61

    Ahora que ya hemos visto distintas ecuaciones de planos, podemos detenernos en el caso particular de los que contienen a los ejes coordenados, estudiemos estos casos:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian

    Para practicar:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    62

    1.5   Ángulos y distancias

    1.5.1   Ángulo entre dos planos

    Sean los planos   π1:  a1x+b1y+c1z+d1=0\; \pi_1:\;{a_1}x + {b_1}y + {c_1}z + {d_1} = 0\: y
    π2:  a2x+b2y+c2z+d2=0.\pi_2:\;{a_2}x + {b_2}y + {c_2}z + {d_2} = 0.

    Dichos planos forman dos ángulos suplementarios, como muestra la figura:


    El ángulo entre dos planos es el ángulo entre sus respectivos vectores normales:

    63

    aˊng(π1,π2)=aˊng(n1,n2)á ng\left( {{\pi _1},{\pi _2}} \right) = á ng\left( {\overrightarrow {{n_1}} ,\overrightarrow {{n_2}} } \right)

    Si llamamos θ  \theta \; a dicho ángulo, resulta:

    cos(θ)=n1.n2n1n2\cos \left( \theta \right) = \frac{{\overrightarrow {{n_1}} .\overrightarrow {{n_2}} }}{{\lVert\overrightarrow {{n_1}}\rVert \lVert\overrightarrow {{n_2}}\lVert }}

    Según el sentido de n1    y    n2\overrightarrow {{n_1}} \;\;y\;\;\overrightarrow {{n_2}} , se obtendrá alguno de los dos ángulos suplementarios. Convenimos en tomar el menor de estos ángulos, por lo cual agregamos módulo en la fórmula anterior:

    cos(θ)=n1.n2n1n2\cos \left( \theta \right) = \frac{{\left| {\overrightarrow {{n_1}} .\overrightarrow {{n_2}} } \right|}}{{\lVert\overrightarrow {{n_1}}\rVert \lVert\overrightarrow {{n_2}}\rVert }}

    Ángulo entre dos planos

      θ=arccos(n1.n2n1n2),      0θπ2  \boxed {\:\: \theta = \arccos \Bigg( {\frac{{\left| {\overrightarrow {{n_1}} .\overrightarrow {{n_2}} } \right|}}{{\lVert\overrightarrow {{n_1}} \rVert \lVert\overrightarrow {{n_2}}\rVert }}} \Bigg),\;\;\;0 \le \theta \le \frac{\pi }{2}\:\:}
    Ejemplo

    Dados los planos: π1:      xy+2=0{\pi _1}:\;\;\;x - y + 2 = 0 π2:(x,y,z)=α(1,2,3)+β(0,1,1)  {\pi _2}:\left( {x,y,z} \right) = \alpha \left( {1,2,3} \right) + \beta \left( {0,1,1} \right)\; Calcular ang(π1,π2)ang\left( {{\pi _1},{\pi _2}} \right).

    64

    El plano π2{\pi _2} está dado en forma vectorial paramétrica, para hallar el ángulo pedido necesitamos n2  \overrightarrow {{n_2}} \;:

    n2=(1,2,3)×(0,1,1)=(1,1,1)\overrightarrow {{n_2}} = \left( {1,2,3} \right) \times \left( {0,1,1} \right) = \left( { - 1, - 1,1} \right)

    n1=(1,1,0)  \overrightarrow {{n_1}} = \left( {1, - 1,0} \right)\;

    n1.n2=0\overrightarrow {{n_1}} .\overrightarrow {{n_2}} = 0

    Esto quiere decir que n1    n2\overrightarrow {{n_1}} \; \bot \; \overrightarrow {{n_2}} , entonces el ángulo es θ=90\theta = 90^\circ .

    La definición de ángulo entre planos nos permite enunciar condiciones de perpendicularidad y de paralelismo entre planos.

    1.5.2   Planos perpendiculares y planos paralelos

    Sean   π1  \;{\pi _1}\; y   π2  \;{\pi _2}\; planos de vectores normales   n1  \;\overrightarrow {{n_1}}\; y   n2\;\overrightarrow {{n_2}}
    respectivamente:

    Planos perpendiculares:     π1    π2        n1.n2=0\;\;{\pi _1} \; \bot \; {\pi _2}\;\; \Leftrightarrow \;\;\overrightarrow {{n_1}} .\overrightarrow {{n_2}} = 0

    Planos paralelos:
    π1π2        n1n2        n1=kn2,    kR{\pi _1} \parallel {\pi _2}\;\; \Leftrightarrow \;\;\overrightarrow {{n_1}} \parallel \overrightarrow {{n_2}} \;\; \Leftrightarrow \;\;\overrightarrow {{n_1}} = k\overrightarrow {{n_2}},\;\;k \in \mathbb{R}

    Consideremos por ejemplo:

    π1:    2x3y+z+1=0                      n1=(2,3,1){\pi _1}:\;\;2x - 3y + z + 1 = 0\;\;\;\;\;\;\;\;\;\;\overrightarrow {\;{n_1}} = \left( {2, - 3,1} \right) π2:    4x6y+2z+5=0                      n2=(4,6,2){\pi _2}:\;\;4x - 6y + 2z + 5 = 0\;\;\;\;\;\;\;\;\;\;\;\overrightarrow {{n_2}} = \left( {4, - 6,2} \right) π3:    4x6y+2z+2=0                      n3=(4,6,2){\pi _3}:\;\;4x - 6y + 2z + 2 = 0\;\;\;\;\;\;\;\;\;\;\;\overrightarrow {{n_3}} = \left( {4, - 6,2} \right)
    65

    Como n2=2  n1,  \overrightarrow {n_2} = 2\;\overrightarrow {n_1},\; podemos afirmar que π1{\pi _1} y π2{\pi _2} son paralelos.

    Análogamente, como n3  =  2  n1,  \overrightarrow {n_3}\; = \;2\;\overrightarrow {n_1},\; los planos π1{\pi _1} y π3{\pi _3} también son paralelos. Pero además se verifica que d3=2  d1{d_3} = 2\;{d_1}, por lo cual π1{\pi _1} y π3{\pi _3} son coincidentes, o sea π1=π3{\pi _1} = {\pi _3}.

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    66

    1.5.3   Distancia de un punto a un plano

    Dados un plano π:    ax+by+cz+d=0\pi :\;\;ax + by + cz + d = 0 y un punto A(xA,yA,zA)A\left( {{x_A},{y_A},{z_A}} \right) , nos proponemos calcular la distancia de AA a π\pi .

    La distancia de AA a π\pi es la longitud del segmento, siendo la proyección ortogonal (perpendicular) de AA sobre π\pi .
    Consideremos un punto cualquiera P(x,y,z)P\left( {x,y,z} \right) perteneciente a π{\rm{\pi }}.

    AA=proyn(PA)\overrightarrow {A'A} = {\overrightarrow {proy} _{\vec n}}\left( {\overrightarrow {PA} } \right)

    Entonces dist(A,π)=proyn(PA)dist\left( {A,\pi } \right) = \lVert{\overrightarrow {proy} _{\vec n}}\left( {\overrightarrow {PA} } \right)\rVert siendo PP un punto cualquiera del plano.

    67

    Veamos un ejemplo, dados:

    π:    x+2y+3z+1=0\pi :\;\;x + 2y + 3z + 1 = 0

    A(0,2,1)A\left( {0,2,1} \right)

    Calcular dist(A,π)dist\left( {A,\pi } \right)

    Recordemos que la norma de la proyección de un vector en la dirección de otro se calcula así:

    proyn(PA)=PA.nn\lVert{\overrightarrow {proy} _{\vec n}}\left( {\overrightarrow {PA} } \right) \rVert = \frac{{\left| {\overrightarrow {PA} .\vec n} \right|}}{{\lVert\vec n\rVert}}

    Tomemos un punto cualquiera del plano, por ejemplo P(1,0,0).  P\left( { - 1,0,0} \right).\;Entonces

    PA=(1,2,1)\overrightarrow {PA} = \left( {1,2,1} \right)

    dist(A,π)=proyn(PA)=PA.nn=dist\left( {A,\pi } \right) = \lVert {\overrightarrow {proy} _{\vec n}}\left( {\overrightarrow {PA} } \right)\rVert = \frac{{\left| {\overrightarrow {PA} .\vec n} \right|}}{{\lVert\vec n\rVert}} = PA.nn=(1,2,1).(1,2,3)12+22+32=814 \frac{{\left| {\overrightarrow {PA} .\vec n} \right|}}{{\lVert\vec n\rVert}} = \frac{{\left| {\left( {1,2,1} \right).\left( {1,2,3} \right)} \right|}}{{\sqrt {{1^2} + {2^2} + {3^2}} }} = \frac{8}{{\sqrt {14} }}

    A continuación deduciremos una fórmula que permite calcular en forma muy sencilla la distancia de un punto a un plano.

    68

    Sean:
    π:    ax+by+cz+d=0\pi :\;\;ax + by + cz + d = 0

    A(xA,yA,zA)A\left( {{x_A},{y_A},{z_A}} \right)

    Habíamos visto que:

    d(A,π)=proyn(PA)            siendo  P(x,y,z)πd\left( {A,\pi } \right) = \lVert{\overrightarrow {proy} _{\vec n}}\left( {\overrightarrow {PA} } \right)\rVert\;\;\;\;\;\;siendo\;P\left( {x,y,z} \right) \in \pi

    PA=(xAx,yAy,zAz)\overrightarrow {PA} = \left( {{x_A} - x,{y_A} - y,{z_A} - z} \right)

    n=(a,b,c)\vec n = \left( {a,b,c} \right)

    Entonces:

    proyn(PA)=PA.nn=a(xAx)+b(yAy)+c(zAz)n\lVert{\overrightarrow {proy} _{\vec n}} ( {\overrightarrow {PA} } )\rVert = \frac{{| {\overrightarrow {PA} .\vec n} |}}{{\lVert\vec n\rVert}} = \frac{{| {a( {{x_A} - x} ) + b( {{y_A} - y}) + c( {{z_A} - z})} |}}{{\lVert\vec n\rVert}}

    proyn(PA)=axA+byA+czA+da2+b2+c2\lVert{\overrightarrow {proy} _{\vec n}}\left( {\overrightarrow {PA} } \right) \lVert= \frac{{\left| {a{x_A} + b{y_A} + c{z_A} + d} \right|}}{{\sqrt {{a^2} + {b^2} + {c^2}} }}

    Concluimos que:

    Distancia de punto a plano

      dist(A,π)=a  xA+b  yA+c  zA+da2+b2+c2     \boxed {\;dist\left( {A,\pi } \right) = \frac{{\left| {a\;{x_A} + b\;{y_A} + c\;{z_A} + d} \right|}}{{\sqrt {{a^2} + {b^2} + {c^2}} }}\;\;}
    69

    Retomemos el ejemplo que habíamos desarrollado:

    π:    x+2y+3z+1=0\pi :\;\;x + 2y + 3z + 1 = 0

    A(0,2,1)A\left( {0,2,1} \right)

    De acuerdo con la fórmula demostrada, la distancia es:

    dist(A,π)=  0+4+3+112+22+32=  814dist\left( {A,\pi } \right) = \frac{{\left| {\;0 + 4 + 3 + 1} \right|}}{{\sqrt {{1^2} + {2^2} + {3^2}} }} = \;\frac{8}{{\sqrt {14} }}

    Tal como habíamos calculado antes pero… ¡más fácil!

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    70

    1.5.4   Distancia entre planos paralelos

    Dados dos planos π1{\pi _1} y π2      {\pi _2}\;\;\;paralelos, ¿cómo podemos hallar la distancia entre ambos?

    Consideremos los siguientes planos paralelos:

    π1:  2x3y+z+1=0,        π2:  4x6y+2z+5=0{\pi _1}:\;2x -3y + z + 1 = 0,\;\;\;\; {\pi _2}:\;4x - 6y + 2z + 5 = 0

    Todos los puntos de π1{\pi _1} están a la misma distancia de π2,  {\pi _2},\; por lo tanto podemos elegir un punto cualquiera de π1{\pi _1} y calcular su distancia a π2.{\pi _2}.\, Por ejemplo:   P1(0,0,1)\;{P_1}\left( {0,0, - 1} \right)

    Aplicamos la fórmula de distancia de un punto a un plano:

    d(π1,π2)=d(P1,π2)=4.06.0+2.(1)+542+(6)2+22=356d\left( {{\pi _1},{\pi _2}} \right) = d\left( {{P_1},{\pi _2}} \right) = \frac{{\left| {4.0 - 6.0 + 2.\left( { - 1} \right) + 5} \right|}}{{\sqrt {{4^2} + {{\left( { - 6} \right)}^2} + {2^2}} }} = \frac{3}{{\sqrt {56} }}

    Observación: Si los planos no son paralelos, la distancia entre ambos es 0.

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    71

    1.6   Haz de planos

    Sean π1{\pi _1} y π2{\pi _2} dos planos no paralelos:

    π1:    a1  x+b1  y+c1  z+d1=0{\pi _1}:\;\;{a_1}\;x + {b_1}\;y + {c_1}\;z + {d_1} = 0

    π2:    a2  x+b2  y+c2  z+d2=0{\pi _2}:\;\;{a_2}\;x + {b_{2\;}}y + {c_2}\;z + {d_2} = 0

    La intersección de dos planos no paralelos es una recta. Se denomina haz de planos al conjunto de planos que pasan por dicha recta. Uno podría imaginarse al haz de planos como si fueran las hojas de un libro abierto:

    72

    Haz de planos que pasan por la recta     r=π1π2\;\;r = {\pi _1} \cap {\pi _2}

    Puede demostrarse que la ecuación del haz de planos que pasan por la recta de intersección entre   π1  y    π2  \;{\pi _1}\;y\;\;{\pi _2}\; es la siguiente:

    k1(a1x+b1y+c1z+d1)+k2(a2x+b2y+c2z+d2)=0{k_1}\left( {{a_1}x + {b_1}y + {c_1}z + {d_1}} \right) + {k_2}\left( {{a_2}x + {b_2}y + {c_2}z + {d_2}} \right) = 0 con   k1,k2R\;{k_1},{k_2} \in \mathbb{R}

    Para cada par de valores de k1{k_1} y k2{k_2} se obtiene un plano que pasa por la recta rr.

    Si   k1=0  \;{k_1} = 0\; y   k20,  \;{k_2} \ne 0,\; se obtiene la ecuación del plano   π2.\;{\pi _2}.

    Si   k2=0  \;{k_2} = 0\; y   k10,  \;{k_1} \ne 0,\; se obtiene la ecuación del plano π1.\:{\pi _1}.

    Si suponemos que alguna de las constantes es diferente de cero, por ejemplo k10,  {k_1} \ne 0,\; podemos dividir la ecuación del haz por k1{k_1}:

    k1k1(a1x+b1y+c1z+d1)+k2k1(a2x+b2y+c2z+d2)=0\frac{{{k_1}}}{{{k_1}}}\left( {{a_1}x + {b_1}y + {c_1}z + {d_1}} \right) + \frac{{{k_2}}}{{{k_1}}}\left( {{a_2}x + {b_2}y + {c_2}z + {d_2}} \right) = 0

    Y renombrando k2k1=k{\Large \frac{{{k_2}}}{{{k_1}}}} = k queda:

    (a1x+b1y+c1z+d1)+k(a2x+b2y+c2z+d2)=0\left( {{a_1}x + {b_1}y + {c_1}z + {d_1}} \right) + k\left( {{a_2}x + {b_2}y + {c_2}z + {d_2}} \right) = 0

    Esta expresión se llama "ecuación del haz reducido". ¿De dónde proviene el nombre de «reducido»?

    73

    Falta π2  {\pi _2}\; porque π2  {\pi _2}\; se corresponde con k1=0{k_1} = 0. Por lo tanto, en el haz reducido están todos los planos que pasan por r=π1π2r = {\pi _1} \cap {\pi _2} excepto π2.{\pi _2}.

      Ejemplo

    Dados los planos:

    π1:  x+2y+3z+1=0{\pi _1}:\;x + 2y + 3z + 1 = 0

    π2:    3x5y+z+10=0{\pi _2}:\;\;3x - 5y + z + 10 = 0

    Encontrar la ecuación de un plano que pase por la recta de intersección entre π1  {\pi _1}\; y π2  {\pi _2}\; y que:

    a) Sea paralelo al eje xx

    b) Sea perpendicular al plano x+y+z=0x + y + z = 0

    Se pide encontrar un plano que pase por la recta de intersección entre π1{\pi _1} y π2,  {\pi _2},\; entonces podemos armar el haz de planos que pasa por dicha recta:

    α(x+2y+3z+1)+β(3x5y+z+10)=0\alpha \left( {x + 2y + 3z + 1} \right) + \beta \left( {3x - 5y + z + 10} \right) = 0

    a) Para que el plano sea paralelo al eje xx su vector normal debe ser n=(0,  b,  c)\vec n = \left( {0,\;b,\;c} \right). O sea, el plano debe ser de la forma: by+cz+d=0by + cz + d = 0.

    74

    Reescribimos el haz como sigue:

    (α+3β)ax+(2α5β)by+(3α+β)cz+(α+10β)d=0 \underbrace{(\alpha + 3\beta)}_{\text{a}}x + \underbrace{(2\alpha - 5\beta)}_{\text{b}}y + \underbrace{(3\alpha + \beta)}_{\text{c}}z + \underbrace{(\alpha + 10\beta)}_{\text{d}} = 0

    Para que sea paralelo al eje xx,    α+3β=0    α=3β\;\;\alpha + 3\beta = 0\; \Rightarrow \; \alpha = - 3\beta

    Reemplazando:

    11βy8βz+7β=0 -11\beta y - 8\beta z + 7\beta = 0

    β(11y8z+7)=0\beta \left( { - 11y - 8z + 7} \right) = 0

    π3:  11y8z+7=0{\pi _3}:\; - 11y - 8z + 7 = 0

    b) Ahora queremos un plano del haz que sea perpendicular al plano x+y+z=0x + y + z = 0. Dos planos son perpendiculares si sus vectores normales son perpendiculares.

    Luego:

    (1,1,1)(α+3β,  2α5β,  3α+β)=0  \left( {1,1,1} \right)\left( {\alpha + 3\beta, \; 2\alpha - 5\beta, \; 3\alpha + \beta } \right) = 0\;

    6αβ=06\alpha - \beta = 0

    β=6α\beta = 6\alpha

    75

    Reemplazando:

    19αx28αy+9αz+61α=019\alpha x - 28\alpha y + 9\alpha z + 61\alpha = 0

    α(19x28y+9z+61)=0\alpha \left( {19x - 28y + 9z + 61} \right) = 0

    19x28y+9z+61=019x - 28y + 9z + 61 = 0


    Para practicar:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    76

    1.7   Recta en R3{\mathbb{R}^3}

    1.7.1   Ecuaciones de la recta en R3{\mathbb{R}^3}

    Sabemos que una recta en R2{\mathbb{R}^2} puede expresarse por la ecuación:

    y=ax+by = ax + b

    Pero ¿qué representa esta ecuación en R3{\mathbb{R}^3}? En R3{\mathbb{R}^3} es un plano paralelo al eje zz, y en R2{\mathbb{R}^2} es una recta:

    Para definir un plano es suficiente conocer un vector perpendicular al plano y un punto del mismo. ¿Qué datos permiten definir una recta en R3{\mathbb{R}^3}?

    77

    Para definir en forma vectorial una recta en R3{\mathbb{R}^3}, es suficiente conocer un punto de la recta y un vector director que indique la dirección de la misma, o sea un vector paralelo a la recta.

    1.7.2   Ecuación vectorial de la recta

    Dados un vector v=(v1,v2,v3)\vec v = \left( {{v_1},{v_2},{v_3}} \right) y un punto P0(x0,y0,z0){P_0}\left( {{x_0},{y_0},{z_0}} \right), nos proponemos hallar la ecuación de la recta rr que pasa por el punto P0{P_0} y es paralela al vector v\vec v.

    Consideremos un punto P(x,y,z)P\left( {x,y,z} \right) perteneciente a la recta rr. El vector P0P  \overrightarrow {{P_0}P} \; resultará paralelo al vector director v\vec v:

    P0P=αv\overrightarrow {{P_0}P} = \alpha \vec v

    (xx0,yy0,zz0)=α(v1,v2,v3)\left( {x - {x_0},y - {y_0},z - {z_0}} \right) = \alpha \left( {{v_1},{v_2},{v_3}} \right)

    78

    Ecuación vectorial de la recta

      (x,y,z)=(x0,y0,z0)+α(v1,v2,v3),      αR     \boxed {\;( {x,y,z} ) = ( {{x_0},{y_0},{z_0}} ) + \alpha ( {{v_1},{v_2},{v_3}}),\;\;\;\alpha \in \mathbb{R}\;\;}
    Ejemplo

    Hallar la ecuación vectorial de la recta que pasa por los puntos M(3,2,1)M\left( {3,2,1} \right) y S(1,1,0)S\left( { - 1,1,0} \right).

    Tenemos como datos dos puntos de la recta, entonces los vectores MS\overrightarrow {MS} y SM\overrightarrow {SM} son paralelos a dicha recta. Elegimos uno de ellos como vector director:

    v=MS=(4,1,1)\vec v = \overrightarrow {MS} = \left( { - 4, - 1, - 1} \right)

    Podemos tomar cualquiera de los dos puntos dados cómo punto de paso, por ejemplo M.  M.\; Entonces la ecuación es:

    (x,y,z)=(3,2,1)+α(4,1,1),      αR\left( {x,y,z} \right) = \left( {3,2,1} \right) + \alpha \left( { -4, - 1, - 1} \right),\;\;\;\alpha \in \mathbb{R}

    (ecuación vectorial de la recta MS)

    Para cada valor de αR,  \alpha \in \mathbb{R},\;se obtiene un punto de la recta. Por ejemplo, si α=1\alpha = - 1 se obtiene el punto P1(7,3,2)r.{P_1}\left( {7,3,2} \right) \in r.

    ¿(5,3,1)r  \left( {5, - 3,1} \right) \in r\;?

    79

    Veamos si existe algún valor de α\alpha que verifique esta ecuación vectorial:

    (5,3,1)=(3,2,1)+α(4,1,1)\left( {5, - 3,1} \right) = \left( {3,2,1} \right) + \alpha \left( { - 4, - 1, - 1} \right)

    {34α=  52α=31α=  1 \begin{cases} \begin{align*} 3 - 4 \alpha =& \; 5 \\ 2 - \alpha =& -3 \\ 1- \alpha =& \; 1 \end{align*} \end{cases}

    Este sistema es incompatible, así que el punto no pertenece a la recta.

    ¿Para qué valor de α\alpha se obtiene el punto SS?


    1.7.3   Ecuaciones paramétricas de la recta

    Hemos visto que la ecuación vectorial de una recta es:

    (x,y,z)=(x0,y0,z0)+α(v1,v2,v3)\left( {x,y,z} \right) = \left( {{x_0},{y_0},{z_0}} \right) + \alpha \left( {{v_1},{v_2},{v_3}} \right)

    Por igualdad de vectores:

    {x=x0+α  v1y=y0+α  v2z=z0+α  v3,   αR \begin{cases} \begin{align*} x = x_0 + \alpha \;v_1 \\ y = y_0 + \alpha \;v_2 \\ z = z_0 + \alpha \;v_3 \end{align*} , \:\:\: \alpha \in \mathbb{R} \end{cases}

    Éstas son las ecuaciones cartesianas paramétricas de la recta.

    80

    1.7.4   Ecuaciones simétricas de la recta

    Si   v1,v2,v3  \;{v_1},{v_2},{v_3}\; son distintos de cero, entonces:

    α=xxov1,        α=yyov2,        α=zzov3\alpha = \frac{{x - {x_o}}}{{{v_1}}},\;\;\;\;\alpha = \frac{{y - {y_o}}}{{{v_2}}},\;\;\;\;\alpha = \frac{{z - {z_o}}}{{{v_3}}}

    Igualando, resulta:

    Ecuaciones simétricas de la recta

    xxov1=yyov2=zzov3     \boxed{ \:\frac{{x - x_o}}{{v_1}} = \frac{{y - y_o}}{{v_2}} = \frac{{z - z_o}}{{v_3}}\;\;}
    Ejemplo

    Consideremos la ecuación vectorial de la recta MSMS:

    (x,y,z)=(3,2,1)+α(4,1,1),      αR  \left( {x,y,z} \right) = \left( {3,2,1} \right) + \alpha \left( { - 4, - 1, - 1} \right),\;\;\;\alpha \in \mathbb{R}\;

    ¿Cómo podemos obtener las ecuaciones paramétricas de la recta? Simplemente por igualdad de vectores escribimos:

    {x=34αy=2αz=1α,   αR \begin{cases} \begin{align*} x &= 3 - 4\alpha \\ y &= 2 - \alpha \\ z &= 1 - \alpha \end{align*}, \:\:\: \alpha \in \mathbb{R} \end{cases}

    Ecuaciones paramétricas de la recta MS

    81

    Para obtener las ecuaciones simétricas, despejamos el parámetro e igualamos:

    α=x34,      α=y21,    α=z11\alpha = \frac{{x - 3}}{{ - 4}},\;\;\;\alpha = \frac{{y - 2}}{{ - 1}},\;\;\alpha = \frac{{z - 1}}{{ - 1}}

    x34=y2=z1,  \dfrac{{x - 3}}{4} = y - 2 = z - 1,\; ecuaciones simétricas de la recta MS.

    Revisión:

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    82
    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian

    1.7.5   Recta definida como intersección de dos planos

    Dos planos no paralelos
    π1:a1x+b1y+c1z+d1=0{\pi _1}:{a_1}x + {b_1}y + {c_1}z + {d_1} = 0 π2:a2x+b2y+c2z+d2=0{\pi _2}:{a_2}x + {b_2}y + {c_2}z + {d_2} = 0 determinan al cortarse una recta en R3\mathbb{R^3} que queda expresada por el sistema de ecuaciones lineales:

    {a1x+b1y+c1z+d1=0a2x+b2y+c2z+d2=0 \begin{cases} \begin{align*} {a_1}x + {b_1}y + {c_1}z + {d_1} = 0 \\ {a_2}x + {b_2}y + {c_2}z + {d_2} = 0 \end{align*} \end{cases}
    83
    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    84

    1.8   Recta y plano: intersecciones y ángulos

    1.8.1   Intersección entre recta y plano

    ¿Qué casos pueden presentarse en la intersección entre una recta y un plano?

    Caso 1

    Una recta puede ser concurrente con un plano:

    rπ={P}r \cap \pi = \left\{ P \right\}



    Caso 2

    Una recta puede ser paralela a un plano:
    rπr\parallel \pi rπ=r \cap \pi = \emptyset

    Caso 3

    Una recta puede estar incluida en un plano:
    rπr \subset \pi rπ=rr \cap \pi = r

    85

    Dados:

    π:    2x3y+z+1=0\pi :\;\;2x - 3y + z + 1 = 0

    r1:(x,y,z)=(0,1,3)+α(1,0,1){r_1}:\left( {x,y,z} \right) = \left( {0,1,3} \right) + \alpha \left( {1,0,1} \right)

    ¿Cómo se busca la intersección entre la recta y el plano?

    Escribimos las ecuaciones paramétricas de la recta y las reemplazamos en la ecuación del plano:

    {x=αy=1z=3+α \begin{cases} \begin{align*} x &= \alpha \\ y &= 1 \\ z &= 3 + \alpha \end{align*} \end{cases}

    2α3.1+(3+α)+1=0      α=132\alpha - 3.1 + \left( {3 + \alpha } \right) + 1 = 0\;\; \Rightarrow \;\alpha = - \frac{1}{3}

    Reemplazando el valor del parámetro α\alpha en las ecuaciones de la recta, obtenemos el punto de intersección:

    r1π={  (13,1,83)  }    {r_1} \cap \pi = \left\{ {\;\left( { - \frac{1}{3},1,\frac{8}{3}} \right)\;} \right\}\;\;

    86

    Busquemos ahora la intersección del mismo plano π\pi con la recta

    r2:  (x,y,z)=(0,0,1)+λ(3,2,0){r_2}:\;\left( {x,y,z} \right) = \left( {0,0, - 1} \right) + \lambda \left( {3,2,0} \right)

    Escribimos las ecuaciones paramétricas:

    {x=3λy=2λz=1 \begin{cases} \begin{align*} x &= 3\lambda \\ y &= 2\lambda \\ z &= -1 \end{align*} \end{cases}

    Reemplazamos en la ecuación del plano

    2(3λ)3(2λ)1+1=0      0=0    λ2\left( {3\lambda } \right) - 3\left( {2\lambda } \right) - 1 + 1 = 0\; \Rightarrow \;\;0 = 0\;\;\forall \lambda

    Queda una expresión que es verdadera para todo λ\lambda . Esto significa que todo punto de la recta verifica la ecuación del plano. En este caso podemos afirmar que la recta está incluida en el plano, por lo tanto: r2π=r2  .{r_2} \cap \pi = {r_2}\;.

    87

    Considerando el mismo plano π\pi , hallemos la intersección con la recta

    r3:(x,y,z)=(5,0,0)+t(0,1,3){r_3}:\left( {x,y,z} \right) = \left( {5,0,0} \right) + t\left( {0,1,3} \right)

    Reiterando el procedimiento, resulta:

    {x=5y=tz=3t \begin{cases} \begin{align*} x &= 5 \\ y &= t \\ z &= 3t \end{align*} \end{cases}

    103t+3t+1=0    11=0        10 - 3t + 3t + 1 = 0\; \Rightarrow \;11 = 0\;\;\;\;absurdo

    Este absurdo nos indica que la recta y el plano no tienen ningún punto en común, o sea que la recta es paralela al plano y por lo tanto: rπ=r \cap \pi = \emptyset

    88

    En resumen:

    Para hallar la intersección entre un plano y una recta, se reemplazan las ecuaciones paramétricas de la recta en la ecuación del plano.

    Pueden presentarse tres casos:

    i) Es posible despejar el valor del parámetro, entonces reemplazando este valor en las ecuaciones de la recta se obtiene el punto de intersección. En este caso:

    rπ={P}r \cap \pi = \left\{ P \right\}

    ii) 0=0rπrπ=r0 = 0 \Rightarrow r \subset \pi \Rightarrow r \cap \pi = r

    iii) 0=k(con  k0)Absurdo0 = k\left( {con\;k \ne 0} \right) \Rightarrow Absurdo rπrπ= \Rightarrow r\parallel \pi \Rightarrow r \cap \pi = \emptyset

    1.8.2   Paralelismo entre recta y plano

    ¿Existe una manera de anticipar si una recta es paralela a un plano sin buscar la intersección?

    Una vez más, los vectores resultarán una herramienta potente para la geometría de rectas y planos.

    89
    Observemos la siguiente figura:

    ¿Cómo deben ser el vector normal del plano y el vector director de la recta para que rπr\parallel \pi?

    rπ    v  n    v  .n=0r\parallel \pi \; \Leftrightarrow \; \vec v \: \bot \; \vec n\; \Leftrightarrow \;\vec v\;.\vec n = 0

    ¿Qué ocurre si la recta está incluida en el plano?

    En este caso también se verifica que el vector director de la recta es perpendicular al normal del plano. Pero a diferencia del caso anterior, todos los puntos de la recta están en el plano. Esto nos permite afirmar que:

    90
    rπ    {v    nPrπ r \subset \pi \; \Leftrightarrow \;\begin{cases} \vec v \; \bot \; \vec n \\ P_r \in \pi \end{cases}

    Ejemplo

    Dados el plano π:  x+yz3=0\pi:\;x + y - z - 3 = 0 y la recta
    r:(x,y,z)=(1,0,0)+t(0,2,2)r:\left( {x,y,z} \right) = \left( {1,0,0} \right) + t\left( {0,2,2} \right), comprobar que la recta es paralela al plano. ¿Está incluida en el plano?

    Si la recta es paralela al plano entonces su vector director v\vec v debe ser perpendicular al vector normal del plano n\vec n. Luego n.v\vec n.\vec v debe ser cero:

    (1,1,1)(0,2,2)=22=0\left( {1,1, - 1} \right)\left( {0,2,2} \right) = 2 - 2 = 0

    Para saber si la recta está incluida en el plano veamos si el punto (1,0,0)\left( {1,0,0} \right) satisface la ecuación del plano π\pi :

    1  +003=0  2=0    Absurdo1\; + 0 - 0 - 3 = 0 \Rightarrow \; - 2 = 0\;\;Absurdo

    Como el punto no satisface la ecuación podemos concluir que rr no está incluida en π\pi .

    Ejercicio para el lector

    Sea   π\;\pi el plano paralelo al eje yy que pasa por (1,1,1)\left( {1,1,1} \right) y (1,2,3)\left( {1,2,3} \right),

    91

    y sea la recta

    {xy=0x+kz=2 \begin{cases} \begin{align*} x - y &= 0 \\ x + kz &= 2 \end{align*} \end{cases}

    Hallar, si es posible, el valor de kk para que la recta rr sea paralela a π\pi .
    Si existe kk, analizar si rπr \subset \pi .

    Respuesta:   k=0  \;k=0\; y la recta no pertenece al plano.


    1.8.3   Perpendicularidad entre recta y plano

    Así como encontramos una condición vectorial para el paralelismo entre una recta y un plano, nos preguntamos si existirá una condición para la perpendicularidad. Observemos la siguiente figura:

    92

    ¿Cómo deben ser el vector normal del plano y el vector director de la recta para que r    πr \; \bot \; \pi?

    r    π        v  n      v=kn      r \; \bot \; \pi \;\; \Leftrightarrow \;\;\vec v\parallel \;\vec n\; \Leftrightarrow \;\;\vec v = k\vec n\;\;\;

    Ejemplo

    Dado el plano π:  x3y+z+1=0\pi:\;x-3y + z + 1 = 0 , hallar la ecuación de la recta perpendicular a π\pi que pasa por A(1,0,3)A\left( {1,0,3} \right).

    Como la recta es perpendicular al plano π\pi entonces su vector director es paralelo al vector normal del plano. Podemos tomar:

    v=(1,3,1)\vec v = \left( {1, - 3,1} \right)

    Ya tenemos el vector director y un punto de paso, luego la ecuación vectorial es:

    r:(x,y,z)=(1,0,3)+λ(1,3,1),      λRr:\left( {x,y,z} \right) = \left( {1,0,3} \right) + \lambda \left( {1, - 3,1} \right),\;\;\;\lambda \in \mathbb{R}

    Ejercicio para el lector

    Dado el haz de planos x+yz+2+k(xy+z)=0,x + y -z + 2 + k\left( {x - y + z} \right) = 0, analizar si existe algún plano del haz que sea perpendicular a la recta r:{x1=0yz=0 r: \begin{cases} \begin{align*} x - 1 &= 0 \\ y - z &= 0 \end{align*} \end{cases}

    93

    Para reforzar y practicar:

    Objetos interactivos diseñados por Elena Esperanza Álvarez Sáiz y M. J. García Cebrian
    94

    1.8.4   Ángulo entre recta y plano

    Consideremos el siguiente esquema:

    Sea rr una recta no paralela ni perpendicular a un plano π\pi . Sea la proyección ortogonal de rr sobre   π\;\pi . Se define como ángulo entre rr y   π\;\pi al ángulo agudo que forman rr y rr'.
    ¿Cómo podemos hallar dicho ángulo? Veamos la siguiente figura:

    95

    Podemos calcular el ángulo β\beta entre v  \vec v\; y   n\;\vec n :

    cos(β)=v.nv  n\cos \left( \beta \right) = \frac{{\vec v.\vec n}}{{\lVert\vec v\rVert\;\lVert\vec n\rVert}}

    Como habíamos visto previamente, de acuerdo con las direcciones de v  \vec v\; y n\vec n el ángulo β\beta podría ser mayor de 90º. Agregamos módulo para obtener un ángulo agudo:

    cos(β)=v.nv  n  \cos \left( \beta \right) = \frac{{\left| {\vec v .\vec n} \right|}}{{\lVert\vec v\rVert\;\lVert\vec n\rVert}}\;

    Como α\alpha y β\beta son complementarios se cumple la igualdad sen(α)=cos(β),sen( \alpha ) = \cos ( \beta), entonces:

    sen(α)=  v.n  v  n      sen\left( \alpha \right) = \frac{{\left| {\;\vec v.\vec n\;} \right|}}{{\lVert\vec v\rVert\;\lVert\vec n\rVert}}\;\;\;

    Por lo tanto:

    Ángulo entre recta y plano

        α=arcsen(v.nv  n),    0απ2    \boxed {\;\; \alpha = arcsen \left( \frac{\left| {\vec v .\vec n} \right|}{\lVert\vec v\rVert\;\lVert\vec n\rVert} \right),\;\;0 \le \alpha \le \frac{\pi }{2}\;\;}
    96

    Casos particulares:

    • Si   α=0,  {\rm{\;}}\alpha = 0,\; entonces la recta es paralela al plano (v.n=0\vec v.\vec n = 0).
    • Si α=π2,\alpha = \dfrac{\pi}{2}, entonces la recta es perpendicular al plano (v=kn\vec v = k\vec n)

    Ejemplo

    Hallar el ángulo entre el plano π:  x+z8=0\pi:\;x + z - 8 = 0 y el eje xx.

    Conocemos el vector normal de π\pi y el vector director de la recta. Luego aplicamos directamente la fórmula:

    α=arcsen(  v.n  v  n)=arcsen((1,0,0).(1,0,1)2)=\alpha = {\rm{arcsen}}\left( {\frac{{\left| {\overrightarrow {\;v} .\overrightarrow {n\;} } \right|}}{{\lVert\vec v\rVert\;\lVert\vec n\rVert}}} \right) = {\rm{arcsen}}\left( {\frac{{\left| {\left( {1,0,0} \right).\left( {1,0,1} \right)} \right|}}{{\sqrt 2 }}} \right) = α=arcsen(12)=π4\alpha = {\rm{arcsen}}\left( {\frac{1}{{\sqrt 2 }}} \right) = \frac{\pi }{4}
    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    97

    1.8.5   Intersección de rectas en R3{\mathbb{R}^3}

    Sabemos que dos rectas en R2{\mathbb{R}^2} o bien se cortan en un único punto o bien son paralelas.

    Pero en R3{\mathbb{R}^3}, además de estos dos casos, existen rectas que ni se cortan ni son paralelas: son las rectas alabeadas.

    Tenemos entonces tres casos en R3{\mathbb{R}^3}:

    Caso 1: Rectas concurrentes o incidentes

    r1r2={P}{r_1} \cap {r_2} = \left\{ P \right\}



    Caso 2: Rectas paralelas
    r1r2    v1=kv2{r_1}\parallel {r_2}\; \Leftrightarrow \;{\vec v_1} = k{\vec v_2}

    Las rectas paralelas podrían ser coincidentes. Para verificar si dos rectas paralelas son coincidentes basta con ver si un punto de una de ellas pertenece o no a la otra recta.

    98
    Caso 3: Rectas alabeadas

    Existe otra posición posible para las rectas en R3{\mathbb{R}^3}. Consideremos el siguiente esquema en el que las rectas r1{r_1}, r2{r_2} y r3{r_3} contienen a las aristas de un cubo.

    Las rectas r1{r_1} y r3{r_3} son paralelas. En cambio r1{r_1} y r2{r_2}, que no son paralelas ni concurrentes, se denominan alabeadas.

    r1r_1 y r2r_2 son alabeadas     r1r2r1r2=\; \Leftrightarrow \; r_1 \nparallel r_2 \land r_1 \cap r_2 = \empty


    99

    1.8.6   Plano que contiene a dos rectas

    Dos rectas en R3{\mathbb{R}^3} se denominan coplanares si existe un plano que contiene a ambas rectas.

    Habíamos visto que dos rectas en R3{\mathbb{R}^3} pueden ser concurrentes, paralelas o alabeadas. Veamos en qué casos es posible encontrar un plano que las contenga:

    Caso 1: Rectas concurrentes

    Veamos el siguiente gráfico que muestra dos rectas concurrentes:

    Dadas dos rectas concurrentes r1{r_1} y r2{r_2}, ¿cómo podemos encontrar el vector normal del plano?

    v1×v2=n{\vec v_1} \times {\vec v_2} = \vec n

    Para completar la ecuación del plano, tomamos un punto P de cualquiera de las dos rectas.

    100

    Ejercicio de aplicación

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    101
    Caso 2: Rectas paralelas

    Veamos el siguiente gráfico que muestra dos rectas paralelas y el plano que las contiene:

    Dadas dos rectas paralelas, existe un plano que las contiene pero no podemos hallar el vector normal como en el caso de las rectas concurrentes. ¿Por qué?

    Caso 3: Rectas alabeadas

    Dos rectas alabeadas no pueden estar contenidas en un mismo plano.

    Ejemplos resueltos
    102

    1.8.7   Ángulo entre dos rectas

    Definición: El ángulo entre dos rectas de R3{\mathbb{R}^3} es el ángulo entre sus vectores directores.

    Sean las rectas r1{r_1} y r2{r_2} con vectores directores v1\vec v_1 y v2.\vec v_2. Entonces:

    cos(α)=v1.v2v1  v2 \cos \left( \alpha \right) = \frac{\vec v_1 .\vec v_2 }{\lVert\vec v_1 \rVert \;\lVert\vec v_2\rVert \ }

    Usamos la misma convención que para ángulo entre planos y para ángulo entre recta y plano, y aplicamos módulo:

    cos(α)=v1.v2v1  v2\cos \left( \alpha \right) = \frac{\left| \vec v_1 .\vec v_2 \right|}{\lVert\vec v_1\rVert \;\lVert\vec v_2\rVert }

    Por lo tanto:

    Ángulo entre dos rectas

        α=arcos(v1.v2v1   v2 ),      0απ2      \boxed {\;\;\alpha = arcos\left( {\frac{{\left| {\overrightarrow {{v_1}} .\overrightarrow {{v_2}} } \right|}}{{\lVert\overrightarrow {{v_1}}\rVert \ \;\lVert\overrightarrow {{v_2}}\rVert \ }}} \right),\;\;\;0 \le \alpha \le \frac{\pi }{2}\;\;\;}

    Casos particulares:

    Si α=0,  \,\alpha = 0,\; entonces las rectas son paralelas (v1=k  v2)(\vec v_1 = k\;\vec v_2) .

    Si α=π2,  \alpha = \dfrac{\pi}{2},\; entonces las rectas son perpendiculares (v1.v2=0)(\vec v_1 .\vec v_2 = 0).

    103

    Observación: La definición de ángulo que hemos adoptado no toma en cuenta si las rectas se cortan o no.

    Ejemplo

    Sean,

    r1:   (x,y,z)=  λ(1,0,0),      λRr2:    (x,y,z)=(0,1,0)+  t(0,0,1),    tR \begin{matrix} \begin{align*} &{r_1}:\;\:\left( {x,y,z} \right) = \;\lambda\left( {1,0,0} \right),\;\;\;\lambda \in \mathbb{R} \\ &{r_2}:\;\;\left( {x,y,z} \right) = \left( {0,1,0} \right) + \; t\left( {0,0,1} \right),\;\;t \in \mathbb{R} \end{align*} \end{matrix}

    Vimos en un ejemplo anterior que estas rectas son alabeadas. ¿Cuál es el ángulo entre ellas?

    Como     v1.v2=0,  \;\;\vec {v_1} .\vec {v_2} = 0,\; las rectas son perpendiculares, o sea que

    α=π2\alpha = \dfrac{\pi}{{2}}
    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    104

    1.9   Distancias y proyecciones

    1.9.1   Proyecciones ortogonales

    Proyección de un punto sobre un plano

    Dados un plano π\pi y un punto AA no perteneciente a dicho plano, la proyección ortogonal de AA sobre π\pi es el punto AπA' \in \pi tal que AA\overrightarrow {AA'} es un vector perpendicular a π\pi.


    A=proyπ(A)    A' = pro{y_\pi }\left( A \right)\;\;si y sólo si     {AπAA    π\;\; \begin{cases} A' \in \pi \\ \overrightarrow{AA'} \; \bot \; \pi \end{cases}

    Ejemplo

    Dados el plano π:  2xy+3z+1=0\pi: \; 2x - y + 3z + 1 = 0 y el punto A  (4,  1,  3)A \;(4,\;1,\;-3), obtener la proyección ortogonal del punto sobre el plano.

    105

    Para obtener buscamos la recta perpendicular a π\pi que pasa por AA. El vector normal del plano es paralelo a esta recta, por lo tanto podemos tomarlo como vector director:

    r:(x,y,z)=(4,1,3)+λ(2,1,3)r:\left( {x,y,z} \right) = \left( {4,1, - 3} \right) + \lambda \left( {2, - 1,3} \right)

    Ahora buscamos la intersección de la recta con el plano, reemplazando las ecuaciones paramétricas de la recta en la ecuación del plano:

    {x=4+2λy=1λz=3+λ \begin{cases} \begin{align*} x &= 4 + 2\lambda \\ y &= 1 - \lambda \\ z &= -3 + \lambda \end{align*} \end{cases}

    2(4+2λ)(1λ)+3.(3+3λ)+1=02\left( {4 + 2\lambda } \right) - \left( {1 - \lambda } \right) + 3.\left( { - 3 + 3\lambda } \right) + 1 = 0

    8+4λ1+λ9+9λ+1=08 + 4\lambda - 1 + \lambda - 9 + 9\lambda + 1 = 0

    14λ1=0λ=11414\lambda - 1 = 0 \Rightarrow \lambda = \frac{1}{{14}}

    Reemplazamos en la ecuación de la recta para obtener las coordenadas del punto A  A'\;:

    (x,y,z)=(4,1,3)+114(2,1,3)=\left( {x,y,z} \right) = \left( {4,1, - 3} \right) + \frac{1}{{14}}\left( {2, - 1,3} \right) = (x,y,z)=(4,1,3)+(17,114,314)=(297,1314,3914)     (x,y,z)= \left( {4,1, - 3} \right) + \left( {\frac{1}{7}, - \frac{1}{{14}},\frac{3}{{14}}} \right) = \left( {\frac{{29}}{7},\frac{{13}}{{14}}, - \frac{{39}}{{14}}} \right)\;\;

    106

    El lector puede verificar que A  (297,1314,3914)π  A'\;\bigg( {\dfrac{29}{7},\dfrac{13}{14}, - \dfrac{39}{14}}\bigg) \in \pi \; y además que  AA    π  .\;\overrightarrow {AA'} \; \bot \; \pi \;.

    Si Aπ  ,A \in \pi \;, ¿cuál es la proyección de AA sobre  π\;\pi ?

    Proyección de una recta sobre un plano

    Dados una recta rr y un plano π\pi , nos interesa obtener la proyección ortogonal de la recta sobre el plano.

    En general, la proyección de una recta sobre un plano es otra recta rr':

    proyπ(r)=rpro{y_\pi }\left( r \right) = r'

    Para hallar rr'\, podemos proyectar dos puntos de la recta sobre el plano. Si llamamos AA y BB a dichos puntos, resulta

    proyπ(A)=A,        proyπ(B)=B        pro{y_\pi }\left( A \right) = A',\;\;\;\;pro{y_\pi }\left( B \right) = B'\;\;\;\;

    107

    La recta rr' que buscamos es la recta determinada por AA' y BB'.

    Pero teniendo en cuenta que el punto de intersección entre la recta y el plano pertenece a rr', es suficiente proyectar un solo punto de r  r\; para que rr' quede definida.

    También existe un caso especial: si la recta es perpendicular al plano, su proyección es un punto.

    rπproyπ(r)=A=rπr \bot \pi \Rightarrow pro{y_\pi }\left( r \right) = A = r \cap \pi

    Ejemplo

    Dado π:x+y+z3=0\pi : x + y + z - 3 = 0
    a) Hallar la proyección de la recta r:(x,y,z)=λ(0,2,1)r:\left( {x,y,z} \right) = \lambda \left( {0,2,1} \right) sobre π\pi
    b) Dada la recta  s:(x,y,z)=(1,0,0)+t.(2,b,c)\;s:\left( {x,y,z} \right) = \left( {1,0,0} \right) + t.\left( { - 2,b,c} \right), hallar los valores de bb y cc para que la proyección de ss sobre π\pi sea un punto. ¿Cuál es dicho punto?

    108

    Resolución

    a) Para hallar la recta r=proyπ(r)r' = pro{y_\pi }\left( r \right) buscamos dos puntos que pertenezcan a rr'.

    Uno de los puntos puede ser el de intersección:

    2λ+λ3=0  λ=1(x,y,z)=(0,2,1)2\lambda + \lambda - 3 = 0 \Rightarrow \;\lambda = 1 \Rightarrow \left( {x,y,z} \right) = \left( {0,2,1} \right)

    Tomemos el punto (0,0,0)r\left( {0,0,0} \right) \in r. Para proyectarlo sobre el plano buscamos la recta perpendicular al plano que pasa por (0,0,0)\left( {0,0,0} \right):

    (x,y,z)=(0,0,0)+α(1,1,1)\left( {x,y,z} \right) = \left( {0,0,0} \right) + \alpha \left( {1,1,1} \right)

    (x,y,z)=α(1,1,1)\left( {x,y,z} \right) = \alpha \left( {1,1,1} \right)

    Reemplazando en la ecuación del plano:

    α+α+α3=0α=1\alpha + \alpha + \alpha - 3 = 0 \Rightarrow \alpha = 1

    Entonces la proyección del punto (0,0,0)\left( {0,0,0} \right) sobre π\pi da (1,1,1)  \left( {1,1,1} \right)\;

    La recta rr' queda definida por los puntos (0,2,1)\left( {0,2,1} \right) y (1,1,1)\left( {1,1,1} \right):

    r:(x,y,z)=(0,2,1)+γ(1,1,0)r':\left( {x,y,z} \right) = \left( {0,2,1} \right) + \gamma \left( {1, - 1,0} \right)

    109

    b) Para que la proyección de ss sobre π\pi sea un punto, el vector director de ss debe ser paralelo al vector normal del plano. Debe existir un kRk \in \mathbb{R} tal que:

    (1,1,1)=k(2,b,c)\left( {1,1,1} \right) = k\left( { - 2,b,c} \right)

    De aquí deducimos que deben ser b=c=2b = c = - 2. La ecuación de la recta queda:
    s:  (x,y,z)=(1,0,0)+t(2,2,2)s:\;\left( {x,y,z} \right) = \left( {1,0,0} \right) + t\left( { - 2, - 2, - 2} \right) Para hallar el punto reemplazamos las ecuaciones paramétricas de la recta en la ecuación del plano:

    s:  {x=12ty=2tz=2t   tR,      π:x+y+z3=0 s:\; \begin{cases} \begin{align*} x &= 1 - 2t \\ y &= -2t \\ z &= -2t \end{align*} \:\:\: t \in \mathbb{R}, \;\;\; \pi:\, x+y+z-3=0 \end{cases}

    12t2t2t3=026t=0t=131 - 2t - 2t - 2t - 3 = 0 \Rightarrow - 2 - 6t = 0 \Rightarrow t = - \frac{1}{3}

    El punto es: P(53,23,23)P\left( {\frac{5}{3},\frac{2}{3},\frac{2}{3}} \right) Por lo tanto: proyπ(s)=  sπ=  P(53,23,23)pro{y_\pi }\left( s \right) = \;s \cap \pi = \;P\left( {\frac{5}{3},\frac{2}{3},\frac{2}{3}} \right)

    110

    En el siguiente interactivo ejercitamos la proyección de una recta sobre un plano.

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    111

    1.9.2   Planos proyectantes de una recta

    Los planos proyectantes de una recta son aquellos planos que incluyen a la recta y son perpendiculares a los planos coordenados.

    Ejemplo 1

    Hallar los planos proyectantes de la siguiente recta:

    r:  {x2y+z+7=0x+y+z1=0 r: \; \begin{cases} \begin{align*} x - 2y + z + 7 &= 0 \\ -x + y + z - 1 &= 0 \end{align*} \end{cases}

    Resolución
    Teniendo en cuenta que la recta está definida como intersección de dos planos, una forma práctica de hallar los planos proyectantes es construir el haz de planos que pasan por dicha recta:

    112

    Construyamos el haz de planos que pasa por rr:

    H:α(x2y+z+7)+β(x+y+z1)=0\mathcal{H}:\alpha \left( {x - 2y + z + 7} \right) + \beta \left( { - x + y + z - 1} \right) = 0

    Ahora distribuimos y reordenamos:

    (αβ)ax+(2α+β)by+(α+β)cz+(7αβ)d \underbrace{(\alpha - \beta)}_{a}x + \underbrace{(-2\alpha + \beta)}_{b}y + \underbrace{(\alpha + \beta)}_{c}z + \underbrace{(7\alpha - \beta)}_{d}

    El plano proyectante perpendicular al plano xyxy (o sea, paralelo al eje zz) es un plano de la forma ax  +  by  +  d  =  0ax\; + \;by\; + \;d\; = \;0. Por lo tanto:

    α+β=0α=  β\alpha + \beta = 0 \Rightarrow \alpha = - \;\beta

        π1:    2x+3y8=0 \Rightarrow \;\;{\pi _1}:\;\; - 2x + 3y - 8 = 0

    El plano proyectante perpendicular al plano xzxz (o sea, paralelo al eje y) es un plano de la forma ax  +  cz  +  d  =  0ax\; + \;cz\; + \;d\; = \;0. Por lo tanto:

    2α+β=0β=2α - 2\alpha + \beta = 0 \Rightarrow \beta = 2\alpha

      π2:    x+3z+5=0 \Rightarrow \;{\pi _2}:\;\; - x + 3z + 5 = 0

    ¿Cuál es el plano proyectante perpendicular al plano yzyz?

    113
    Ejemplo 2

    Hallar los planos proyectantes de la recta que pasa por los puntos A(1,2,3)A\left( {1,2,3} \right) y B(3,1,1)B\left( {3,1, - 1} \right).

    Resolución

    La ecuación vectorial de la recta que contiene a los puntos AA y BB es r:(x,y,z)=(1,2,3)+t(2,1,4)r: (x,y,z) = (1,2,3) +t(2, - 1, - 4)

    Supongamos que llamamos π1\pi _1 al plano proyectante que es perpendicular al plano xy.  xy.\; ¿Qué condiciones debe cumplir n1\overrightarrow {{n_1}} \,?

    Entonces: (a,b,0).(2,1,4)=02ab=0b=2a\left( {a,b,0} \right).\left( {2, - 1, - 4} \right) = 0 \Rightarrow 2a - b = 0 \Rightarrow b = 2a

    Con a=1      y    b=2\,a = 1\;\;\;y\;\;b = 2 , resulta:

    Teniendo en cuenta que rπ1r \subset {\pi _1} , reemplazamos un punto de la recta (por ejemplo AA) para obtener dd.

    El plano buscado es: π1:  x+2y5=0{\pi _1}:\;x + 2y - 5 = 0

    Dejamos a cargo del lector comprobar que los otros planos proyectantes son:

    π2:  2x+z5=0{\pi _2}:\;2x + z - 5 = 0\, y   π3:4y+z+5=0\;{\pi _3}: - 4y + z + 5 = 0

    114
    Ejemplo 3

    Retomemos el ejemplo anterior, siendo
    r:  (x,y,z)=(1,2,3)+  t(2,1,4)r:\;\left( {x,y,z} \right) = \left( {1,2,3} \right) + \;t\left( {2, - 1, - 4} \right)

    Cuando las componentes del vector director de una recta son distintas de cero, podemos expresarla a través de sus ecuaciones simétricas. En este caso:

    r:          x12=y21=z34r:\;\;\;\;\;\frac{{x - 1}}{2} = \frac{{y - 2}}{{ - 1}} = \frac{{z - 3}}{{ - 4}}

    A partir de las ecuaciones simétricas se deducen tres igualdades, cada una de las cuales se corresponde con la ecuación de un plano proyectante de la recta:

    x12=y21    (x1)=2(y2)   \dfrac{x-1}{2}=\dfrac{y-2}{-1}\; \Rightarrow\;-(x-1)=2(y-2)\;\Rightarrow

    x+2y5=0            \quad\quad\quad\quad \Rightarrow x + 2y - 5 = 0\;\;\;\;\;\; Plano proyectante π1\pi_1

    x12=z34    4(x1)=2(z3)   \dfrac{x-1}{2}= \dfrac{z-3}{-4}\; \Rightarrow\;-4(x-1)=2(z-3)\;\Rightarrow

    2x+z5=0            \quad\quad\quad\quad \Rightarrow 2x + z - 5 = 0\;\;\;\;\;\; Plano proyectante π2\pi_2

    y21=z34    4(y2)=(z3)   \dfrac{y-2}{-1} = \dfrac{z-3}{-4}\; \Rightarrow\;-4(y-2)=-(z-3)\;\Rightarrow

    4y+z+5=0            \quad\quad\quad\quad \Rightarrow -4y +z + 5 = 0\;\;\;\;\;\; Plano proyectante π3\pi_3

    115

    El siguiente gráfico muestra la recta y su plano proyectante π1:{\pi _1}:

    Recta y un plano proyectante

    En la figura puede observarse que el plano proyectante π1{\pi _1} es el plano determinado por la recta  r\;r y su proyección rr' sobre el plano coordenado xyxy. Lo mismo puede afirmarse de los otros dos planos proyectantes.

    116
    Actividad (para reflexionar)

    Hallar los planos proyectantes de cada una de las siguientes rectas:

    r:    (x,y,z)  =  (1,2,3)  +  t(1,2,0)r:\;\;\left( {x,y,z} \right)\; = \;\left( {1,2,3} \right)\; + \;t\left( {1, - 2,0} \right)

    s:      (x,y,z)  =  (2,2,1)  +  λ(0,0,1)s:\;\;\;\left( {x,y,z} \right)\; = \;\left( {2,2,1} \right)\; + \;\lambda \left( {0,0,1} \right)

    A partir de los resultados obtenidos, ¿podrías establecer qué condición debe cumplirse para que una recta tenga dos de sus planos proyectantes iguales?

    ¿Existe algún caso en que no esté definido alguno de los planos proyectantes de una recta?

    1.9.3   Distancias

    a) Distancia punto-recta en R3{\mathbb{R}^3}

    Dados un punto AA y una recta rr de vector director v\vec v, queremos hallar la distancia entre AA y rr (con ArA \notin r). Sea ArA' \in r tal que AA\overrightarrow {AA'} es perpendicular a la recta.

    d(A,r)=AAd\left( {A,r} \right) = \lVert\overrightarrow {AA'} \rVert

    Consideremos un punto PrP \in r y el vector PA\overrightarrow {PA} , y construyamos el paralelogramo determinado por v\vec v y PA\overrightarrow {PA} , tal como lo muestra la figura:

    117


    El segmento AAAA' es la altura del paralelogramo. Si llamamos S al área de dicho paralelogramo, resulta:

    S=b.h=v.d                [1]S = b.h = \lVert\vec v\rVert.d\;\;\;\;\;\;\;\;\left[ 1 \right]

    Recordemos que el área del paralelogramo es igual al módulo del producto vectorial:

    S=PA×v              [2]S = \lVert\overrightarrow {PA} \times \vec v \rVert\;\;\;\;\;\;\;\left[ 2 \right]

    Igualando [1]\left[ 1 \right] y [2]\left[ 2 \right] podemos despejar dd:

    v.d=PA×v\lVert \vec v\rVert.d = \lVert\overrightarrow {PA} \times \vec v\rVert

              d(A,r)=PA×vv\;\; \Rightarrow \;\;\;d\left( {A,r} \right) = \frac{{\lVert\overrightarrow {PA} \times \vec v\rVert}}{{\lVert\vec v\rVert}}
    118
    Ejemplo

    Calcular la distancia entre la recta

    r:(x,y,z)=(1,0,1)+λ(2,1,1)r:(x,y,z) = (1,0,- 1) + \lambda (- 2,1,1) y el punto A(3,1,1)A\left( {3, - 1,1} \right).

    Resolución

    Aplicamos la fórmula:

      d(A,r)=PA×vv\;d\left( {A,r} \right) = \frac{{\lVert\overrightarrow {PA} \times \vec v\rVert}}{{\lVert\vec v\rVert}}

    Tomamos un punto cualquiera de la recta, por ejemplo P(1,0,1)P\left( {1,0, - 1} \right) y formamos:

    PA=(2,1,2)\overrightarrow {PA} = \left( {2, - 1,2} \right)

    d(A,r)=(2,1,2)×(2,1,1)  (2,1,1)=(3,6,0)6  =456d\left( {A,r} \right) = \frac{{\lVert\left( {2, - 1,2} \right) \times \left( { - 2,1,1} \right)\rVert}}{{\;\lVert\left( { - 2,1,1} \right)\rVert}} = \frac{{\lVert\left( { - 3, - 6,0} \right)\rVert}}{{\sqrt 6 \;}} = \frac{{\sqrt {45} }}{{\sqrt 6 }}

    d(A,r)=2,74 d(A,r) = \cong 2,74
    Objeto interactivo diseñado por Elena E. Álvarez Sáiz y Mª José García Cebrian
    119
    b) Distancia entre dos rectas paralelas

    La fórmula que hemos visto permite calcular la distancia entre dos rectas paralelas. Veamos un ejemplo:

    Dadas las rectas r1:(x,y,z)=(1,0,0)+t(2,1,1)  {r_1}:\left( {x,y,z} \right) = \left( {1,0,0} \right) + t\left( {2, - 1,1} \right)\; y
    r2:{y+z=2x+ky=0 r_2: \: \begin{cases} y + z = 2 \\ x + ky = 0 \end{cases}

    Hallar kk tal que r1r2{r_1}\parallel {r_2} y calcular d(r1,r2).d\left( {{r_1},{r_2}} \right).

    Recordemos la condición para que dos rectas sean paralelas:

    r1r2  v1=α  v2          [1]{r_1}\parallel {r_2}\; \Leftrightarrow \overrightarrow {{v_1}} = \alpha \;\overrightarrow {{v_2}} \;\;\;\;\;\left[ 1 \right]

    La recta r2{r_2} está definida como intersección de dos planos. Si hacemos el producto vectorial de los vectores normales tendremos un vector director de la recta:

    (0,1,1)×(1,k,0)=(k,1,1)\left( {0,1,1} \right) \times \left( {1,k,0} \right) = \left( { - k,1, - 1} \right)

    Por [1]:

    (k,  1,  1)=α(2,  1;1)    {k=2α1=α1=α(-k, \; 1,\; -1) = \alpha (2,\;-1;1)\; \Rightarrow \; \begin{cases} \begin{align*} -k &= 2\alpha \\ 1 &= - \alpha \\ -1 &= \alpha \end{align*} \end{cases}

    k=2 \Rightarrow k = 2

    120

    Para calcular la distancia, tomemos dos puntos cualesquiera de las rectas y construyamos el vector P1P2\overrightarrow {{P_1}{P_2}} :


    P1(1,0,0)r1,      P2(0,0,2)r2,      P1P2=(1,0,2){P_1}\left( {1,0,0} \right) \in {r_1},\;\;\;{P_2}\left( {0,0,2} \right) \in {r_2},\;\;\;\overrightarrow {{P_1}{P_2}} = \left( { - 1,0,2} \right)

    La distancia entre las rectas será la distancia entre P2{P_2} y r1{r_1}:

    d(r1,r2)=d(P2,r1)=P1P2×v1v1d\left( {{r_1},{r_2}} \right) = d\left( {{P_2},{r_1}} \right) = \frac{{\overrightarrow {{P_1}{P_2}} \times \overrightarrow {{v_1}} }}{{\overrightarrow {{v_1}} }}

    P1P2×v1=(1,0,2)×(2,1,1)=(2,5,1)\overrightarrow {{P_1}{P_2}} \times {\vec v_1} = \left( { - 1,0,2} \right) \times \left( {2, - 1,1} \right) = \left( {2,5,1} \right)

    P1P2×v=30 \lVert \overrightarrow {{P_1}{P_2}} \times \vec v \rVert = \sqrt {30}

    d(r1,r2)=306=5d\left( {{r_1},{r_2}} \right) = \frac{{\sqrt {30} }}{{\sqrt 6 }} = \sqrt 5

    Como puede observarse, con las mismas herramientas resolvimos un problema diferente.

    121
    c) Distancia entre rectas alabeadas

    Dadas dos rectas r1{r_1} y r2{r_2} no paralelas, nos proponemos calcular la distancia entre ambas:

    La mínima distancia entre dos rectas alabeadas r1{r_1} y r2{r_2} se obtiene al proyectar el vector P1P2\overrightarrow {{P_1}{P_2}} sobre la dirección perpendicular a ambas rectas, dada por v1×v2\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} :

    d(r1,r2)=proyv1×v2(P1P2)d\left( {{r_1},{r_2}} \right) = {\lVert\overrightarrow {proy} _{\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} }}\left( {\overrightarrow {{P_1}{P_2}} } \right)\rVert

    Recordemos que:

    proyv(u)=u.vv\lVert{\overrightarrow {proy} _{\vec v}}\left( {\vec u} \right)\rVert = \frac{{\left| {\vec u.\vec v} \right|}}{{\lVert\vec v\rVert}}

    122

    Entonces podemos obtener una fórmula para la distancia entre r1{r_1} y r2{r_2}:

    Distancia entre rectas alabeadas

          d(r1,r2)=P1P2.(v1×v2)v1×v2       \boxed {\;\;\;d\left( {{r_1},{r_2}} \right) = \frac{{\left| {\overrightarrow {{P_1}{P_2}} .\left( {\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} } \right)} \right|}}{{\lVert\overrightarrow {{v_1}} \times \overrightarrow {{v_2}}\rVert }}\;\;\;}

    ¿Qué significa que la distancia dé 00?

    1.9.4   Condición de coplanaridad

    Habíamos visto que tanto las rectas concurrentes como las paralelas son coplanares (existe un plano que las contiene).

    • Si las rectas se cortan, la distancia entre ellas es cero y por lo tanto: (v1×v2)P1P2=0\left( {\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} } \right)\overrightarrow { \cdot {P_1}{P_2}} = 0
    • Si las rectas son paralelas, v1×v2=0\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} = \vec 0 y por lo tanto también se cumple que (v1×v2)P1P2=0\left( {\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} } \right)\overrightarrow { \cdot {P_1}{P_2}} = 0.

    Entonces estamos en condiciones de enunciar una condición de coplanaridad entre dos rectas:

    r1{r_1} y r2{r_2} son coplanares     (v1×v2).P1P2=0 \; \Leftrightarrow \; \left( {\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} } \right).\overrightarrow {{P_1}{P_2}} = 0
    123

      Ejemplo

    Dadas las rectas:

    r1:(x,y,z)=(3,2,5)+λ(0,1,2){r_1}:\left( {x,y,z} \right) = \left( {3,2,5} \right) + \lambda \left( {0, - 1,2} \right) r2:  {2x+y+3z+2=0x+2y4z+3=0{r_2}:\; \begin{cases} \begin{align*} 2x + y + 3z + 2 = 0 \\ - x + 2y - 4z + 3 = 0 \end{align*} \end{cases}

    Calcular d(r1,r2)d\left( {{r_1},{r_2}} \right)

    Resolución

    Debemos verificar que se trata de rectas alabeadas. Busquemos la dirección de la recta r2  {r_2}\;:

    v2=(2,1,3)×(1,2,4)=(10,5,5)\overrightarrow {{v_2}} = \left( {2,1,3} \right) \times \left( { - 1,2, - 4} \right) = \left( { - 10,5,5} \right)

    Vemos que las rectas no son paralelas porque sus vectores directores no son paralelos. Luego podemos utilizar la fórmula para distancia entre rectas alabeadas:

    d(r1,r2)=P1P2.(v1×v2)v1×v2d\left( {{r_1},{r_2}} \right) = \frac{{\left| {\overrightarrow {{P_1}{P_2}} .\left( {\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} } \right)} \right|}}{{\lVert\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} \rVert}}

    Para hallar  P2  r2  \;{P_2}\; \in {r_2}\;fijamos z=0z = 0, por ejemplo, y averiguamos los valores de xx e yy resolviendo el sistema de ecuaciones:

    124
    {2x+y+2=0x+2y+3=0        x=15,      y=85 \begin{cases} \begin{align*} 2x + y + 2 = 0 \\ - x + 2y + 3 = 0 \end{align*} \end{cases} \;\; \Rightarrow \;\;x = - \frac{1}{5},\;\;\;y = - \frac{8}{5}

    Obtenemos así P2(15,85,0)r2{P_2}\left( { - \dfrac{1}{5}, - \dfrac{8}{5},0} \right) \in {r_2}

    Calculamos:

    P1P2=(165,  185,5)\overrightarrow {{P_1}{P_2}} = \left( { - \frac{{16}}{5},\; - \frac{{18}}{5}, - 5} \right) v1×v2=(0,1,2)×(10,5,5)=(15,  20,10)\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} = \left( {0, - 1,2} \right) \times \left( { - 10,5,5} \right) = \left( { - 15,\; - 20, -10} \right) v1×v2=(15)2+(20)2+(10)2=725 \lVert\overrightarrow {{v_1}} \times \overrightarrow {{v_2}}\rVert = \sqrt {{{\left( { - 15} \right)}^2} + {{\left( { - 20} \right)}^2} + {{\left( { - 10} \right)}^2}} = \sqrt {725} P1P2.(v1×v2)=(165,  185,5).(15,20,10)\overrightarrow {{P_1}{P_2}} .\left( {\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} } \right) = \left( { - \frac{{16}}{5},\; - \frac{{18}}{5}, - 5} \right).\left( { - 15, - 20, - 10} \right) P1P2.(v1×v2)=48+72+50=170 \overrightarrow {{P_1}{P_2}} .\left( {\overrightarrow {{v_1}} \times \overrightarrow {{v_2}} } \right) = 48 + 72 + 50 = 170 d(r1,r2)=1707256,31d\left( {{r_1},{r_2}} \right) = \frac{{170}}{{\sqrt {725} }} \cong 6,31
    125

    Applets para repasar y ejercitar distancia entre rectas:

    Objetos interactivos diseñados por E. Álvarez Sáiz y Mª José García Cebrian
    126
    Ejercicio para el lector 1

    Dadas las rectas:

    r1:(x,y,z)=(2,0,1)+λ(3,1,2){r_1}:\left( {x,y,z} \right) = \left( {2, 0, -1} \right) + \lambda \left( {3, 1, 2} \right) r2:  {x+z2=0y=2{r_2}:\; \begin{cases} x + z - 2 = 0 \\ y = 2 \end{cases}

    Hallar la ecuación del plano π\pi que contiene a r1{r_1} y es paralelo a   r2\;{r_2}.

    Respuesta: x5y+z1=0\: x - 5y + z - 1 = 0

    Nota: en el siguiente applet, uno de los ejemplos corresponde a este ejercicio, pero r1r_1 y r2r_2 están escritas en ecuaciones simétricas. No necesariamente es el primer ejemplo que aparece.

    Objeto interactivo diseñado por Elena E. Álvarez Sáiz y Mª José García Cebrian
    127
    Ejercicio para el lector 2

    Dadas las rectas

    r1:   (x,y,z)=  t(1,1,2)r2:    (x,y,z)=(1,k,0)+  λ  (2,1,0) \begin{matrix} \begin{align*} &{r_1}:\;\:\left( {x,y,z} \right) = \;t\left( {1,1,2} \right) \\ &{r_2}:\;\;\left( {x,y,z} \right) = \left( {1,k,0} \right) + \;\lambda \;\left( {2,1,0} \right) \end{align*} \end{matrix}

    Hallar el valor de kk para que las rectas sean coplanares y encontrar el plano que las contiene.

    Respuesta:   k=12,    π:2x4y+z=0\:\: k = \dfrac{1}{2}, \;\; \pi: 2x - 4y + z = 0

    128

    1.10   Ejercicios surtidos de fin de capítulo

    Nota: Todos los objetos interactivos utilizados en este apartado fueron diseñados por E. Álvarez Sáiz y Mª José García Cebrian

    Punto simétrico respecto de un plano

    129

    Distancia recta plano

    Plano por punto medio de un segmento

    130

    Ángulo entre rectas

    Ángulo entre recta y plano

    Distancia de un punto a un eje

    131

    Aplicación de haz de planos

    Determine el plano que contiene a la intersección de los planos

    P1:x+yz=1,    P2:2xy+z=2 P_1: x + y - z = 1,\;\; P_2: 2x - y + z = 2 y pasa por el punto P(2,1,1).P (2,-1,1).

    Planteamos una ecuación simplificada del haz de planos

    P1+kP2=0P_1 + k P_2 = 0 x+yz1+k(2xy+z2)=0 x + y - z - 1 + k(2x - y + z - 2) = 0

    El punto P(2,1,1)   P (2,-1,1)\; pertenece al haz de planos, por lo que los valores del punto se reemplazan para encontrar el valor de kk :
    2111+k(4+1+12)=0 2 - 1 - 1 - 1 + k (4 + 1 + 1 - 2) = 0
    1+4k=0-1 + 4k = 0
    k=14k =\dfrac {1}{4}

    Sustituimos el valor de kk en la ecuación de la intersección de los planos:

    x+yz1+14(2xy+z2)=0 x + y - z - 1 + \dfrac {1}{4} (2x - y + z - 2) = 0

    Simplificando la ecuación del plano obtenida:

    P:2x+yz=2\boxed {P: 2x + y - z = 2}
    132
    Objeto interactivo de construcción propia

    Pasaje de ecuación vectorial del plano a ecuación general

    Objeto interactivo diseñado por Elena Esperanza Álvarez Sáiz y Mª José García Cebrian
    133
    Capítulo II

    Matrices, determinantes y SEL

    2.1   Matrices

    2.1.1   Definición de Matriz

    Una matriz AA de m×nm \times n es un ordenamiento o arreglo rectangular de escalares dispuestos en mm filas y nn columnas. Para designar a cada uno de los m.nm.n elementos de la matriz se utiliza un doble subíndice que indica el número de fila y número de columna que le corresponde en el arreglo:

    A=(a11a12a1na21a22a2nam1am2amn) \mathbf{A} = \begin{pmatrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{m1} & a_{m2} & \cdots & a_{mn} \end{pmatrix} Así,   a34\;{a_{34}} es el elemento ubicado en la fila tres y la columna cuatro y en general aij{a_{ij}} es el elemento de la matriz AA que está en la fila ii y en la columna jj.
    Las matrices suelen designarse con letras mayúsculas: se anota ARmxnA \in {\mathbb{R}^{mxn}} para indicar que es una matriz con mm filas y nn columnas cuyos elementos son números reales. Se indican con paréntesis o con corchetes:

    139

    Por ejemplo una matriz de dos filas y tres columnas se puede escribir así: A=(320241),AR2×3 \mathbf{A} = \begin{pmatrix} 3 & 2 & 0 \\ -2 & 4 & 1 \end{pmatrix}, \quad \mathbf{A} \in {\mathbb{R}^{2 \times 3}}

    En este caso, diremos que el tamaño u orden de AA es 2×32 \times 3.

    2.1.2   Matriz columna

    Podemos pensar los vectores como casos particulares de matrices: C=(201) matriz  o  vector  columnaCR3×1 \mathbf{C} = \begin{pmatrix} 2 \\ 0 \\ 1 \end{pmatrix} \ matriz\;o\;vector\;columna \quad \mathbf{C} \in {\mathbb{R}^{3 \times 1}}

    2.1.3   Matriz fila

    O también: F=(2          0          1)    matriz  o  vector  fila        FR1×3F = \left( {2\;\;\;\;\;0\;\;\;\;\;1} \right)\;\;matriz\;o\;vector\;fila\;\;\;\;F \in {\mathbb{R}^{1 \times 3}}

    2.1.4  Matriz nula

    La matriz nula es aquélla cuyos elementos son todos ceros. La simbolizamos con OO ó NN)

    2.1.5   Igualdad de matrices

    Dos matrices son iguales si son del mismo orden (tamaño) y sus elementos respectivos son iguales. A,BRmxn          A=B    aij=bij    i,jA,B \in {\mathbb{R}^{mxn}}\;\;\;\;\;A = B\; \Leftrightarrow \;{a_{ij}} = {b_{ij}}\;\;\forall i,j

    140

    2.2   Operaciones con matrices

    2.2.1   Suma de matrices

    Sean A,BRm×nA,B \in {\mathbb{R}^{m \times n}} entonces: A+B=CRm×n          cij=aij+bij      i,jA + B = C \in {\mathbb{R}^{m \times n}}\;\;|\;\;\;{c_{ij}} = {a_{ij}} + {b_{ij}}\;\;\;\forall i,j

    Ejemplo

    2.2.2   Producto de un escalar por una matriz

    Sean ARm×n,  kRA \in {\mathbb{R}^{m \times n}},\;k \in \mathbb{R}, entonces: kA=B  Rm×n          bij=kaij      i,jkA = B\; \in {\mathbb{R}^{m \times n}}\;\;|\;\;\;{b_{ij}} = k{a_{ij}}\;\;\;\forall i,j Por ejemplo,

    si A=(011232) \mathbf{A} = \begin{pmatrix} 0 & -1 & 1 \\ -2 & 3 & 2 \end{pmatrix} ,  entonces   3A=(033696)\mathbf{3A} = \begin{pmatrix}0 & -3 & 3\\ -6 & 9 & 6\end{pmatrix}

    Cuando k=1k = -1, obtenemos la matriz opuesta de A: A= (1)A=(011232) \mathbf{-A}= \ (-1) \mathbf{A} = \begin{pmatrix} 0 & 1 & -1 \\ 2 & -3 & -2 \end{pmatrix}

    141
    Podemos así definir la diferencia (resta) entre dos matrices del mismo tamaño: AB=A+(B)A-B = A + \left( {-B} \right) O sea: AB=CRm×n              cij=aijbijA-B = C \in {\mathbb{R}^{m \times n}}\;\;\;\;|\;\;\;{c_{ij}} = {a_{ij}}-{b_{ij}}

    Ejercicio resuelto

    2.2.3   Propiedades de la suma de matrices y del producto por un escalar

    Sean A,B,CRmxn    y    α,  βRA,B,C \in {\mathbb{R}^{mxn}}\;\;y\;\;\alpha ,\;\beta \in \mathbb{R}. Vimos que: A+BRmxnA + B \in {\mathbb{R}^{mxn}} y αARmxn\alpha A \in {\mathbb{R}^{mxn}}. Estas operaciones verifican las siguientes propiedades:

    1. A+B=B+A        A + B = B + A\;\;\;\;
    2. (A+B)+C=A+(B+C)        \left( {A + B} \right) + C = A + \left( {B + C} \right)\;\;\;\;
    3. A+O=O+A=A          A + O = O + A = A\;\;\;\;\;
    4. A+(A)=(A)+A=OA + \left( { -A} \right) = \left( { -A} \right) + A = O
    5. α(A+B)=αA+αB\alpha \left( {A + B} \right) = \alpha A + \alpha B
    6. (α+β)A=αA+βA\left( {\alpha + \beta } \right)A = \alpha A + \beta A
    7. α(βA)=(αβ)A\alpha \left( {\beta A} \right) = \left( {\alpha \beta } \right)A
    8. 1A=A1A = A
    142

    Puede observarse la analogía entre estas propiedades y las que habíamos enunciado en el capítulo anterior para vectores de R3{\mathbb{R}^3}.

    2.2.4   Producto de matrices

    Sean ARm×nA \in {\mathbb{R}^{m \times n}} y BRn×pB \in {\mathbb{R}^{n \times p}}, o sea que se cumple que la cantidad de columnas de la primera matriz es igual a la cantidad de filas de la segunda: Entonces el producto es:

    Una forma alternativa de expresar el producto es:

    A B = C R m ×p | c i j = k = 1 n a i k b k j AB = C \in {\mathbb{R}^{m \times p}}\;\;\;|\;\;\;\;\;{c_{ij}} = \mathop \sum \limits_{k = 1}^n {a_{ik}}\;{b_{kj}}

    Por ejemplo, al multiplicar una matriz A\mathbf{A} de 3×23\times2 por una matriz B\mathbf{B} de 2×32\times3, tendremos como resultado una matriz de dimensión (3×2)(2×3)=3×3(3\times\cancel{2})(\cancel{2}\times3) = 3\times3.

    143
    Objeto interactivo diseñado por Mª José García Cebrian

    2.2.5   Propiedades del producto

    1) El producto de matrices no es conmutativo. ABBAAB\not=BA
    2) (AB)C=A(BC)\left( {AB} \right)C = A\left( {BC} \right) asociatividad
    3) (A+B)  C=AC+BC\left( {A + B} \right)\;C = AC + BC distributividad a derecha
      P(Q+R)=PQ+PR     P\left( {Q + R} \right) = PQ + PR \; distributividad a izquierda
    4) (kA)B=k(AB)=A(kB)    ,    kR\left( {kA} \right)B = k\left( {AB} \right) = A\left( {kB} \right)\;\;,\;\;k \in \mathbb{R}
    5) OA=O      y      AO=O,  OA = O\;\;\;y\;\;\;AO = O,\; siendo OO la matriz nula

    2.2.6   Matriz traspuesta

    La traspuesta de una matriz ARm×nA \in {\mathbb{R}^{m \times n}}, que indicamos como At{A^t}, es la matriz de n×mn \times m que se obtiene a partir de AA cambiando las filas por las columnas.

    Objeto interactivo diseñado por Mª José García Cebrian
    144

    2.2.7   Propiedades de la matriz traspuesta

    1) (A+B)t=At+Bt{\left( {A + B} \right)^t} = {A^t} + {B^t}
    2) (kA)t=kAt,    kR{\left( {kA} \right)^t} = k{A^t},\;\;k \in \mathbb{R}
    3) (At)t=A{\left( {{A^t}} \right)^t} = A
    4) (AB)t=BtAt{\left( {AB} \right)^t} = {B^t}{A^t}

    2.2.8   Matrices cuadradas: identidad e inversa

    Si AA es una matriz con m=nm = n se dice que es una matriz cuadrada. La matriz identidad I\mathbf{I} es una matriz cuadrada cuyos elementos de la diagonal principal son iguales a 11 y todos los demás son 00. Por ejemplo la matriz identidad de 3×33 \times 3 es

    I=(100010001) \mathbf{I} = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} Una matriz cuadrada AA es inversible, o regular , si existe otra matriz cuadrada (su inversa), que escribiremos A1A^{-1}, que cumple: A1.A=A.A1=I. A^{-1}.A = A.A^{-1} = I. Cuando AA no tiene inversa se dice que AA es singular.

    Objeto interactivo diseñado por Mª José García Cebrian
    145

    2.2.9   Propiedades de la matriz inversa

    Sean A,BRnxnA,B \in {\mathbb{R}^{nxn}} inversibles.  Entonces:

    1) ABAB es inversible y su inversa es: (AB)1=B1A1 \left(AB\right)^{-1}=B^{-1}A^{-1}

    Esto significa que la inversa de ABAB es B1A1{B^{-1}}{A^{-1}}, es decir que el producto de matrices inversibles es inversible.
    ¿Ocurre lo mismo con la suma de matrices inversibles?

    2) (kA)1=1k    A1      (k0){\left( {kA} \right)^{-1}} = \frac{1}{{k\;}}\;{A^{-1}}\;\;\;\left( {k \ne 0} \right)

    3) (At)1=(A1)t{\left( {{A^t}} \right)^{-1}} = {\left( {{A^{-1}}} \right)^t}

    2.2.10   Algunos tipos de matrices

    Objeto interactivo diseñado por John Jairo García Mora et al.
    146

      Ejercicios y problemas

    1) En una ciudad hay tres empresas de telefonía celular (A, B y C) que controlan el mercado.

    Inicialmente cada empresa tiene una fracción de la clientela que denominaremos a0{a_0}, b0{b_0} y c0{c_0}.

    Entonces resulta: a0  +  b0  +  c0=1{a_0}\; + \;{b_0}\; + \;{c_0} = 1 (no hay otras empresas)
    La figura resume el porcentaje de clientes que cambian de empresa durante un período de seis meses.

    Este modelo matemático se basa en los siguientes supuestos:
    – El porcentaje de cambio entre las empresas se mantiene constante con el tiempo.
    Los clientes seguirán siendo consumidores de una de estas tres empresas.
    – No se incorporan nuevos clientes al sistema.

    Llamemos X0=(a0b0c0)X_0=\begin{pmatrix} a_0 \\ b_0\\c_0\end{pmatrix} al vector de estado inicial, y X1=(a1b1c1)X_1=\begin{pmatrix} a_1\\ b_1\\c_1\end{pmatrix} al vector que indica la fracción de la clientela que corresponde a cada empresa al cabo de un semestre.

    Veamos cómo puede obtenerse X1{X_1} a partir de X0{X_0} .

    147

    De acuerdo con la figura, podemos deducir que al cabo de un período (semestre) la empresa A conservará 70% de su clientela propia.

    ¿Qué porcentaje de su clientela conservarán las empresas B y C al cabo de un semestre?

    2) Analice si las siguientes proposiciones son verdderas o falsas. Justifique la respuesta.

    148

    2.3   Determinantes

    ¿Cómo podemos saber si una matriz tiene inversa? El determinante de una matriz proporciona información para responder a esta pregunta.
    Cuando vimos producto vectorial y mixto, habíamos definido determinantes de orden 2 y de orden 3. Recordamos aquí las fórmulas presentadas:
    A cada matriz cuadrada puede asignársele un número real que llamaremos su determinante y designaremos como det(A) oA.det(A)\space o \begin{vmatrix} A\end{vmatrix}.

    Para matrices 2×22×2 y 3×33×3 el determinante se calcula como sigue: A = ( a11 a 12 a21 a22 ) | A | = | a 11 a 12 a 21 a 22 | = a 11 a 22 a 12 a 21

    A=(a11a12a13a21a22a23a31a32a33)A =\begin{pmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33}\end{pmatrix} \rArr A=a11a22a23a32a33a12a21a23a31a33+a13a21a22a31a32\begin{vmatrix} A\end{vmatrix} = a_{11} \begin{vmatrix} a_{22} & a_{23} \\ a_{32} & a_{33} \end{vmatrix} - a_{12} \begin{vmatrix} a_{21} & a_{23} \\ a_{31} & a_{33} \end{vmatrix} + a_{13} \begin{vmatrix} a_{21} & a_{22} \\ a_{31} & a_{32} \end{vmatrix}

    Observación: El determinante no está definido para matrices rectangulares.
    La regla de Sarrus es una forma práctica de calcular determinantes, sólo aplicable para matrices de 3×33 \times 3.

    149

    Consideremos el siguiente esquema en el cual agregamos al final de una matriz de 3×33 \times 3 las filas 1 y 2. El determinante se calcula sumando los productos indicados por las flechas que que van de izquierda a derecha y restando los productos indicados por las flechas que van de derecha a izquierda:

    Objeto interactivo diseñado por Mª José García Cebrian

    2.3.1 Desarrollo de un determinante por cofactores

    Dada ARnxnA \in {\mathbb{R}^{nxn}}, se denomina menor Mij{M_{ij}} a la submatriz de (n1)×  (n1)\left( {n - 1} \right) \times \;\left( {n - 1} \right) que se obtiene a partir de AA eliminando la fila ii y la columna jj.

    Se denomina cofactor Cij{C_{ij}} del elemento aij{a_{ij}} al producto de (1)i+j{\left( { - 1} \right)^{i + j}} por el determinante de la matriz menor Mij  {M_{ij}}\;:
    Cij=(1)i+jdet(Mij){C_{ij}} = {\left( { - 1} \right)^{i + j}}\det \left( {{M_{ij}}} \right)

    150

    Por ejemplo para la matriz A=(314121321)A =\begin{pmatrix} 3 & 1 & 4 \\ -1 & 2 & 1 \\ 3 & 2 & 1\end{pmatrix} Las matrices menores y los cofactores de la primera fila son:

    Observación: (1)(i+j)={1si i+j  es  par1si i+j  es  impar (-1)^{(i+j)} = \begin{cases} 1 &\text{si } i + j \; es \; par \\ -1 &\text{si } i + j \; es \; impar\end{cases}Con estas definiciones previas, estamos en condiciones de enunciar el desarrollo de un determinante de orden nn.
    El determinante de una matriz puede calcularse utilizando los cofactores de cualquier fila o cualquier columna.

    Desarrollo por fila ii:
    det(A)=ai1Ci1+ai2Ci2++ainCin=j=1naijCij  \text{det}(A) = a_{i1} C_{i1} + a_{i2} C_{i2} + \dots + a_{in} C_{in} = \sum_{j = 1}^{n} a_{ij} C_{ij} \ Este cálculo se puede hacer tomando cualquiera de las filas de A A , o sea:   i=1,  2,  ,  n\;i = 1,\;2,\; \ldots ,\;n

    151

    Desarrollo por columna jj:
    det(A)=a1jC1j+a2jC2j++anjCnj=i=1naijCij  \text{det}(A) = a_{1j} C_{1j} + a_{2j} C_{2j} + \dots + a_{nj} C_{nj} = \sum_{i = 1}^{n} a_{ij} C_{ij} \ Este cálculo se puede hacer tomando cualquier columna de A A , o sea:   j=1,  2,  ,  n\;j = 1,\;2,\; \ldots ,\;n

    Objeto interactivo diseñado por Mª José García Cebrian

    2.3.2   Determinante de una matriz triangular

    Si AR3×3A \in {\mathbb{R}^{3 \times 3}} es triangular: A=(a11a12a130a22a2300a33)A =\begin{pmatrix} a_{11} & a_{12} & a_{13} \\ 0 & a_{22} & a_{23} \\ 0 & 0 & a_{33}\end{pmatrix}
    152

    Podemos calcular el determinante de AA por la primera columna:

    A=a11a22a230a33=a11a22a33 \begin{vmatrix} A\end{vmatrix} = a_{11} \begin{vmatrix} a_{22} & a_{23} \\ 0 & a_{33} \end{vmatrix} = a_{11} a_{22} a_{33}

    Hemos llegado a la siguiente conclusión: si AR3×3 A \in {\mathbb{R}^{3 \times 3}} es triangular, su determinante es el producto de los elementos de la diagonal principal.
    Esto se puede generalizar para matrices de cualquier orden: si ARn×nA \in {\mathbb{R}^{n \times n}} es triangular, su determinante es el producto de los elementos de la diagonal principal.

    ARn×n  triangulardet(A)=a11.a22ann A \in {\mathbb{R}^{n \times n}} \; triangular \Rightarrow \det \left( A \right) = {a_{11}}.{a_{22}} \ldots {a_{nn}} Teniendo en cuenta que la matriz identidad es triangular, se deduce que: det(I)=1det(I)=1

    2.3.3  Propiedades de los determinantes

    En general (es decir, a menos que la matriz sea triangular o tenga alguna otra cualidad especial), el cálculo de determinantes por medio del desarrollo por cofactores no es eficiente por el número de operaciones que implica cuando se trabaja con matrices grandes.

    Existen propiedades que facilitan dicho cálculo y que analizaremos aquí en particular para matrices de orden 3, pero se extienden a las columnas y a matrices de cualquier orden.

    153

    Observación: Si ARn×nA \in {\mathbb{R}^{n \times n}}, en algunos casos utilizaremos notación:

    A=(A1    A2          An)A = \left( {{A_1}\;\;{A_2}\;\; \ldots \;\;\;{A_n}} \right)

    donde Aj{A_j} representa la columna jj de la matriz AA.

      Ejercicios y problemas

    1) Demostrar que si kNk \in \mathbb{N},   det(Ak)=[det(A)]k\; \text{det}(A^k) = [\text{det}(A)]^k

    2) Demostrar que si ARn×n  A \in {\mathbb{R}^{n \times n}}\; y kR\,k \in \mathbb{R} entonces det(kA)=kndet(A)\quad \text{det}(kA) = k^n\,\text{det}(A)

    3) Sea AR3×3A \in {\mathbb{R}^{3 \times 3}}\, y   det(A)=2\; \text{det}(A) = 2\, calcular   det(32A3)\; \text{det}(\frac 32\,A^3)\,

    154

    4) Sea A=(A1A2A3)R3×3,  det(A)=k0  A = \begin{pmatrix} A_{1} & A_{2} & A_{3}\end{pmatrix} \in {\mathbb{R}^{3 \times 3}},\; \text{det}(A)=k\neq 0 \; y
    B=(A2B1A1+3A3B2A1A2B3) B = \big( \underbrace{A_{2}}_{B_{1}} \quad \underbrace{A_{1}+3A_{3}}_{B_{2}} \quad \underbrace{A_{1}-A_{2}}_{B_{3}}\, \big) Mostrar que B es inversible y calcular det(2AtB1) \text{det}(-2A^tB^{-1})



    2.3.4  Obtención de la inversa a través de la adjunta

    Para calcular la matriz inversa primero necesitamos definir a la matriz adjunta de una matriz cuadrada AA.

    Llamaremos Cof(A)Cof(A) a la matriz de cofactores de AA, es decir, cada elemento (i,j)(i,j) de Cof(A)Cof(A) es el cofactor (ver §2.3.1)

    Cij=(1)i+j  det(Mij).C_{ij} = (-1)^{i+j}\; \text{det}({M}_{ij}). Entonces, podemos definir a la matriz adjunta de AA como a la traspuesta de Cof(A)Cof(A). Adj(A)=(Cof(A))tAdj(A) =(Cof(A))^{t}

    Si ARn×n  A \in {\mathbb{R}^{n \times n}} \; puede demostrarse que:

    A.Adj(A)=det(A).IAdj(A).A=det(A).I A.Adj(A)=\text{det}(A).I \quad \land \quad Adj(A).A=\text{det}(A).I
    155
    Por lo tanto podemos afirmar que si det(A)0,   \text{det}(A)\neq 0,\; la inversa de AA es: A1=1det(A)Adj(A) A^{-1}=\frac{1}{\text{det}(A)}\, Adj(A)

    En particular, calcular la matriz inversa de una matriz de 2×22\times2 es muy sencillo.

    Consideremos la matriz AA de 2×22\times2 con det(A)0\text{det(A)} \not = 0, se tiene

    A1=1a11a22a12a21[a22a12a21a11]. A^{-1} = \frac{1}{a_{11}a_{22} - a_{12}a_{21}} \begin{bmatrix} a_{22} & -a_{12} \\ -a_{21} & a_{11} \end{bmatrix}.

    Ejemplo del cálculo de la inversa de una matriz.

    Objeto interactivo diseñado por Melissa Méndez Servin
    Se observa cada paso usando los botónes Prev, Sig y Reset (para regresar al primer paso). O bien, ver el desarrollo con otra matriz presionando sobre el botón Otro. Esperar a que automáticamente el applet haga los cálculos.
    156

    2.4 Matrices y sistemas de ecuaciones lineales

    A continuación disponemos de una presentación donde se exponen los contenidos vinculados a sistemas de ecuaciones lineales, su forma matricial, el método de eliminación de Gauss y el teorema de Rouché-Frobenius.

    El material original de la presentación es autoría de la profesora María José García Cebrian, (disponible en https://proyectodescartes.org/iCartesiLibri/materiales_didacticos/AlgebraLinealBachillerato-JS/index.html ) y de los profesores Isabel Pustilnik y Federico Gómez, contenido de su página web Matrices y sistemas de ecuaciones lineales


    157

    2.4.1  Sistemas homogéneos

    Un sistema es homogéneo si todos los términos independientes son cero.
    Estos sistemas siempre son compatibles ya que r(A) = r(A'), todos tienen al menos la solución x1 = x2 = ... = xn = 0, llamada solución trivial, y lo que interesa es saber si hay más soluciones.


    Objeto interactivo diseñado por Mª José García Cebrian
    158

    2.4.2  Sistemas con parámetros

    Son sistemas en los que algunos coeficientes y/o términos independientes se desconocen. Se trata de familias de sistemas de ecuaciones y para cada valor del parámetro o los parámetros se tiene un sistema diferente.
    En estos sistemas se analizan los valores de el o los parámetros para distintos valores de los rangos de A y de A'. Se hallan dichos rangos para estos valores del parámetro y se aplica el Teorema de Rouché-Frobenius.

    Objeto interactivo diseñado por Mª José García Cebrian
    159

      Ejercicios y problemas

    3) A continuación se presentan más ejercicios para practicar. Todos tienen sus respectivas soluciones.

    160
    4)

    5)

    Sean el plano π:  x+y2=0 \pi :\;x + y - 2 = 0 \, y la recta r:{2y+3kz2=03x+ky3z4=0r: \begin{cases} 2y+3kz-2=0 \\ 3x+ky-3z-4=0\end{cases}

    Obtener los valores de k para los cuales:

    a) la recta corta al plano en un único punto;
    b) la recta no interseca al plano;
    c) la recta está incluida en el plano.

    Sugerencia: Este ejercicio puede responderse sin hacer cálculos, teniendo en cuenta los resultados del ejemplo anterior.

    161
    Capítulo III

    Espacios vectoriales

    3.1   Espacios y subespacios vectoriales

    Los vectores geométricos y las matrices presentan similitudes. Las propiedades de la suma y del producto por un escalar son idénticas en ambos conjuntos. Podemos generalizar el concepto de vector a partir de estas propiedades en común.

    3.1.1   Definición de espacio vectorial

    Un espacio vectorial es un conjunto no vacío VV de objetos (vectores), en el que se definen dos operaciones: la suma y el producto por un escalar; sujetas a diez axiomas que se dan a continuación. u,v,wV  \forall u, v, w \in V \; y α,βR,\,\forall {\alpha}, {\beta} \in \R \,, los axiomas deben ser válidos. Llamamos u+vu\, + \,v a la suma de vectores en VV, y αv\alpha v al producto de un escalar α\alpha por un vector vVv \in V.
    1. u+v      Vu + v\;\; \in \;V
    2. u+v=v+uu + v = v + u
    3. (u+v)+w=u+(v+w)\left( {u + v} \right) + w = u + \left( {v + w} \right)
    4. Existe un vector nulo 0VV{0_V} \in V tal que v+0V=vv + {0_V} = v
    5. vV,  \forall v \in V, \; existe un opuesto (v)V\left( { -v} \right) \in V tal que v+(v)=0Vv + \left( {- v} \right) = {0_V}
    6. αvV\alpha v \in V
    7. α(u+v)=αu+αv\alpha \left( {u + v} \right) = \alpha u + \alpha v
    8. (α+β)v=αv+βv\left( {\alpha + \beta } \right)v = \alpha v + \beta v
    9. α(βv)=(αβ)v\alpha \left( {\beta v} \right) = \left( {\alpha \beta } \right)v
    10. 1v=v1v = v
    167

    Observación: En la definición anterior, cuando decimos «escalares» nos estamos refiriendo a números reales. En este caso, se dice que VV es un espacio vectorial real.
    Pero es posible que los escalares pertenezcan a otro conjunto numérico.

    Ejemplo 1

    De acuerdo con las propiedades que vimos en el primer capitulo, podemos afirmar que R3{\mathbb{R}^3} es un espacio vectorial.
    Los espacios Rn{\mathbb{R}^n} , con n1n \ge 1 , son los ejemplos principales de espacios vectoriales. La intuición geométrica desarrollada para R3{R^3} nos ayudará a entender y visualizar muchos conceptos de esta unidad.
    Los vectores de R  n{\mathbb{R}^{\;n}} son n-uplas de números reales, o sea:

    Rn={(x1,  x2,  ,  xn)  ,      con    xiR}{\mathbb{R}^n} = \left\{ {\left( {{x_1},\;{x_2},\; \ldots ,\;{x_n}} \right)\;,\;\;\;con\;\;{x_i} \in \mathbb{R}} \right\}

    En Rn{\mathbb{R}^n} , la suma de vectores y el producto por un escalar se definen así:

    Sean   u=(u1,u2,,un)      y      v=(v1,v2,vn)  Rn  \;u = \left( {{u_1},{u_2}, \ldots ,{u_n}} \right)\;\;\;y\;\;\;v = \left( {{v_1},{v_2}, \ldots {v_n}} \right)\; \in {\mathbb{R}^n}\;

    u+v=(u1+v1,  u2+v2,  ,  un+vn)Rnu + v = \left( {{u_1} + {v_1},\;{u_2} + {v_2},\; \ldots ,\;{u_n} + {v_n}} \right) \in {\mathbb{R}^n}

    αv=(αv1,  αv2,  ,  αvn)Rn\alpha v = \left( {\alpha {v_1},\;\alpha {v_2},\; \ldots ,\;\alpha {v_n}} \right) \in {\mathbb{R}^n}

    Puede comprobarse que las operaciones definidas verifican los axiomas de espacio vectorial.

    168

    Ejemplo 2

    De acuerdo con las propiedades enunciadas en la segunda unidad, para cada mm y nn, Rmxn{\mathbb{R}^{mxn}} es un espacio vectorial.
    Tenemos por ejemplo R2×3{\mathbb{R}^{2×3}}, espacio vectorial cuyos vectores son las matrices de 2×32 \times 3.

    Ejemplo 3

    Llamemos P2{P_2} al conjunto de polinomios de grado menor o igual que 2, incluyendo el polinomio nulo.
    Recordemos la suma de polinomios y la multiplicación por un escalar:
    Dados p(x)=  ao+a1x+a2x2P2p\left( x \right) = \;{a_o} + {a_1}x + {a_2}{x^2} \in {P_2}
    y q(x)=bo+b1x+b2x2P2q\left( x \right) = {b_o} + {b_1}x + {b_2}{x^2} \in {P_2}
    Definimos las operaciones:
    (p+q)(x)=p(x)+q(x)=({p + q}) (x) =p(x) + q(x) = =(ao+bo)+(a1+b1)x+(a2+b2)x2  P2= ({a_o} + {b_o}) + ({a_1} + {b_1})x + ({a_2} + {b_2}){x^2}\; \in {P_2}

    (αp)(x)=αp(x)=(αao)+(αa1)x+(αa2)x2  P2\left( {\alpha p} \right)\left( x \right) = \alpha p\left( x \right) = \left( {\alpha {a_o}} \right) + \left( {\alpha {a_1}} \right)x + \left( {\alpha {a_2}} \right){x^2}\; \in {P_2}

    Puede demostrarse que estas operaciones verifican todos los axiomas de espacio vectorial.

    En particular, el vector nulo en este espacio es el polinomio nulo, es decir el polinomio cuyos coeficientes son todos iguales a cero.

    169

    Generalizando, para cualquier n0n \ge 0 , el conjunto Pn{P_n} de todos los polinomios de grado menor o igual que nn (incluyendo el polinomio nulo) es un espacio vectorial.

    Observación:

    ¿Por qué no definimos Pn{P_n} como el conjunto de polinomios de grado exactamente igual a nn? Si lo definiéramos así, no sería un espacio vectorial como se muestra en el siguiente ejemplo:

    p(x)=x2  p\left( x \right) = {x^2}\; y q(x)=x2+1q\left( x \right) = - {x^2} + 1 son polinomios de grado 2, pero la suma es un polinomio de grado cero. Entonces no se verificaría el primer axioma de espacio vectorial (la suma de vectores de un espacio vectorial VV debe estar en VV).

    170

    3.1.2   Propiedades de los espacios vectoriales

    A partir de los axiomas de espacios vectoriales, pueden demostrarse estas propiedades que resultan «naturales»:

    Propiedad 1 0  u=0V0\;u = {0_V}

    Propiedad 2 α  0V=0V\alpha \;{0_V} = {0_V}

    Propiedad 3 (α)u=(αu)\left( { - \alpha } \right)u = - \left( {\alpha u} \right)

    En particular, para   α=1\; \alpha =-1 :      (1)u=u\;\;\; \left( { - 1} \right)u = - u

    Propiedad 4 α  u=0V      α=0      u=0V\alpha \;u = {0_V}\; \Rightarrow \;\;\alpha = 0\; \vee \;\;u = {0_V} Veamos cómo puede demostrarse esta última propiedad:
    Si α=0\alpha = 0 , se cumple la proposición.
    Si α0\alpha \ne 0 , podemos multiplicar por 1α\frac{1}{\alpha } :

    α  u=0V  1αα  u=1α  0Vu=0V  \alpha \;u = {0_V}\; \Rightarrow \frac{1}{\alpha }\alpha \;u = \frac{1}{\alpha }\;{0_V} \Rightarrow u = {0_V}\;

    171

    3.1.3   Subespacios vectoriales

    Definición

    Sea VV un espacio vectorial y WW un subconjunto no vacío de VV.

    WW es un subespacio de VV si WW es en sí mismo un espacio vectorial con las mismas operaciones (suma de vectores y producto por un escalar) definidas en VV.

    Ejemplo

    W={(x1,x2)  R2:  x2=3x1}W = \left\{ {\left( {{x_1},{x_2}} \right)\; \in {R^2}:\;{x_2} = 3{x_1}} \right\} ¿es un subespacio de R2{\mathbb{R}^2}?

    Primero analicemos el conjunto WW. Son todos vectores de R2{\mathbb{R}^2} tales que la segunda componente es el triple de la primera:

    (x1,3x1)=x1(1,3)\left( {{x_1},3{x_1}} \right) = {x_1}\left( {1,3} \right)

    WW es la recta que pasa por el origen y tiene vector director (1,3), o sea la recta de ecuación y = 3x.

    Para decidir si WW es un subespacio de R2,{\mathbb{R}^2}, habría que verificar que se cumplan los axiomas del 1 al 10. El lector puede comprobar que todos se cumplen en este caso.

    Pero en general no es necesario verificar los axiomas porque existen condiciones para determinar si un subconjunto WW de un espacio vectorial VV es un subespacio, son las siguientes:

    172

    Condiciones necesarias y suficientes para caracterizar subespacios

    Sea WW un subconjunto de un espacio vectorial VV (WV)\left( {W \subseteq V} \right).
    WW es subespacio de VV si y sólo si se cumplen las siguientes condiciones:
    a.  0V{0_V} está en WW.
    b.  Si uu y vv están en WW, entonces u+vu + v está en WW.
    c.  Si uu está en WW y kk es un escalar, kuku está en WW.

    Observaciones

    1.  La condición (a) asegura que W no es vacío. La mejor manera de comprobar si W es un subespacio es buscar primero si contiene al vector nulo. Si 0V{0_V} está en W, entonces deben verificarse las propiedades (b) y (c). Si   0V\;{0_V} no está en W, W no puede ser un subespacio y no hace falta verificar las otras propiedades.
    2.  Las propiedades a, b y c corresponden a los axiomas 4, 1 y 6 de espacios vectoriales.
    3.  Los axiomas 2, 3, 7, 8, 9 y 10 de espacio vectorial se cumplen para WW porque éste es un subconjunto de VV. Puede decirse que WW «hereda» esas propiedades de VV.
    4.  Faltaría comprobar que cada vector de WW tiene su opuesto en WW (axioma 5 de espacios vectoriales):
    Teniendo en cuenta la condición (c) de subespacios,
    si tomamos k=1,  k = -1,\; resulta:
    Para cada uW  ,(1)u=u  Wu \in W\;,\left( {-1} \right)u = -u\; \in W.
    Y por lo tanto cada vector de WW tiene su opuesto en WW.

    173

    De las observaciones anteriores se deduce que las condiciones (a), (b) y (c) son suficientes para demostrar que WW es un espacio vectorial, y por lo tanto subespacio de VV.

    3.1.4   Subespacios triviales

    Si VV es un espacio vectorial, entonces VV es un subespacio de sí mismo.
    Lo mismo sucede con {0V}\left\{ {{0_V}} \right\} ya que   0V+0V=0V\;{0_V} + {0_V} = {0_V}
    y k0V=0V  para  cualquier  k  realk{0_V} = {0_V}\;para\;cualquier\;k\;real

    Se dice que {0V}\left\{ {{0_V}}\right\} y VV son los subespacios triviales de VV.

      Ejercitación sobre subespacios

    Ejemplo 1

    Consideremos el conjunto W={(x,y)R2  /  xy=0}W = \left\{ {\left( {x,y} \right) \in {\mathbb{R}^2}\;/\;xy = 0} \right\}, ¿Es un subespacio de R2{\mathbb{R}^2}?

    Se cumple (a) pues (0,0)W\left( {0,0} \right) \in W

    No se cumple (b) porque la suma de dos vectores de WW puede no estar en WW, por ejemplo:

    (1,0)+(0,1)=(1,1)W\left( {1,0} \right) + \left( {0,1} \right) = \left( {1,1} \right) \notin W

    Entonces WW no es un subespacio de R2{\mathbb{R}^2}.

    174

    Ejemplo 2

    Consideremos el conjunto W={(x,y)R2  /  x=0}W = \left\{ {\left( {x,y} \right) \in {\mathbb{R}^2}\;/\;x = 0} \right\}. Es decir, la recta de ecuación x=0.  x = 0.\;¿Es un subespacio de R2{\mathbb{R}^2}?

    Se cumple (a) pues (0,0)W\left( {0,0} \right) \in W

    Se cumple (b) pues la suma de dos vectores de W,  W,\; está en WW:

    (0,y1)+(0,y2)=(0,y1+y2)\left( {0,{y_1}} \right) + \left( {0,{y_2}} \right) = \left( {0,{y_1} + {y_2}} \right)

    Se cumple (c) pues el producto de un vector de WW por un número real está en WW:

    k(0,y)=(0,ky)k\left( {0,y} \right) = \left( {0,ky} \right)

    Luego WW es subespacio de R2{\mathbb{R}^2}.

    Ejemplo 3

    Consideremos el conjunto W={(x,y)R2/x2y2=0}W = \left\{ {\left( {x,y} \right) \in {\mathbb{R}^2}/{x^2} -{y^2} = 0} \right\}. ¿Es un subespacio de R2{\mathbb{R}^2}?

    x2y2=0y=xx=y{x^2} - {y^2} = 0 \Leftrightarrow y = x \vee x = -y

    Se cumple (a) pues (0,0)W\left( {0,0} \right) \in W

    175

    No se cumple (b) porque la suma de dos vectores de WW puede no estar en WW, por ejemplo:

    (1,1)+(1,1)=(2,0)W\left( {1,1} \right) + \left( {1, -1} \right) = \left( {2,0} \right) \notin W

    Entonces WW no es un subespacio de R2{\mathbb{R}^2}.

    Ejemplo 4

    Consideremos el conjunto W={(x,y,z)R3/x+y+2z=0}W = \left\{ {\left( {x,y,z} \right) \in {\mathbb{R}^3}/x + y + 2z = 0} \right\}. Es decir un plano que pasa por el origen. ¿Es un subespacio de R3{\mathbb{R}^3}?

    De la ecuación del plano se deduce que: x=y2zx = - y - 2z
    Por lo tanto los vectores que pertenecen a WW responden a la forma (y2z,y,z)\left( {-y -2z,y,z} \right) con y,zRy,z \in \mathbb{R}.

    Se cumple (a) pues (0,0,0)W\left( {0,0,0} \right) \in W

    Se cumple (b) pues la suma de dos vectores del plano, sigue estando en ese plano:

    (y2z,  y,z)+(y2z,y,z)=\left( { -y -2z,\;y,z} \right) + \left( {-y’ -2z’,y’,z’} \right)= =(  (y+y)2(z+z)  ,  y+y,z+z  )= \left( {\; -\left( {y + y’} \right) -2\left( {z + z’} \right)\;,\;y + y’,z + z’\;} \right)

    Se cumple (c)\left( c \right) pues k(y2z,y,z)=(ky2kz,ky,kz)Wk\left( { -y -2z,y,z} \right) = \left( { -ky -2kz,ky,kz} \right) \in W

    Entonces WW es subespacio de R3{\mathbb{R}^3}.

    176

    Ejemplo 5

    Consideremos el conjunto W={pP2  /  p(0)=0}W = \left\{ {p \in {P_2}\;/\;p\left( 0 \right) = 0} \right\}.
    Es decir, los polinomios de grado menor o igual que dos (incluyendo el polinomio nulo) tales que evaluados en 00 dan por resultado 00. ¿Es un subespacio de P2{P_2}?

    Se cumple (a) pues el polinomio nulo pertenece a WW.

    Recordemos la definición de suma de funciones y de producto de un real por una función:

    (f  +  g)  (x)  =  f(x)  +  g(x)\left( {f\; + \;g} \right)\;\left( x \right)\; = \;f\left( x \right)\; + \;g\left( x \right), para todo xx perteneciente al dominio de ff y de gg

    (kf)(x)=kf(x)\left( {kf} \right)\left( x \right) = kf\left(x \right) para todo xx perteneciente al dominio de ff.

    Los polinomios son funciones, por lo tanto si consideramos p  ,  qWp\;,\;q \in W, resulta:

    (p+q)(0)=p(0)+q(0)=0+0=0      p+qW\left( {p + q} \right)\left( 0 \right) = p\left( 0 \right) + q\left( 0 \right) = 0 + 0 = 0\; \Rightarrow \;\;p + q \in W
    (kp)(0)=k.p(0)=k  0=0kpW  \left( {kp} \right)\left( 0 \right) = k.p\left( 0 \right) = k\;0 = 0 \Rightarrow kp \in W\;

    Demostramos que WW es un subespacio de P2{P_2}.

    177

    Ejemplo 6

    Consideremos el conjunto W={ARn×n/  A=At}W = \left\{ {A \in {\mathbb{R}^{n \times n}}/\;A = {A^t}} \right\}.
    Es decir, el conjunto de matrices simétricas de n×nn \times n.

    Se cumple (a) porque la matriz nula pertenece a WW.
    Se cumple (b) pues si A,BWA,B \in W entonces (A+B)t=At+Bt={(A + B)^t} = {A^t} + {B^t} = =A+B,= A + B, luego (A+B)W.{(A + B)} \in W.
    Se cumple (c) pues si AWA \in W entonces (kA)t=kAt=kA{\left( {kA} \right)^t} = k{A^t} = kA, luego (kA)W.\left( {kA} \right) \in W.
    Demostramos que el conjunto de matrices simétricas de n×nn\times n es un subespacio de Rn×n\mathbb{R}^{n \times n}.

    Ejemplo 7

    Consideremos el conjunto W={AR2×2/det(A)=0}W = \left\{ {A \in {\mathbb{R}^{2 \times 2}}/\det \left( A \right) = 0} \right\}.
    ¿Es un subespacio de R2×2{\mathbb{R}^{2 \times 2}}?

    Se cumple (a) porque la matriz nula pertenece a WW.
    En general det(A+B)det(A)+det(B)\det \left( {A + B} \right) \ne \det \left( A \right) + \det \left( B \right), entonces podría ocurrir que A,BWA,B \in W pero que A+BA + B no esté en WW. Por ejemplo

    A=(1313),  B=(0025),  A+B=(1332) A =\begin{pmatrix} 1 & -3 \\ -1 & 3\end{pmatrix},\; B =\begin{pmatrix} 0 & 0 \\ -2 & 5\end{pmatrix},\; A + B =\begin{pmatrix} 1 & -3 \\ -3 & 2\end{pmatrix} A+BWA + B \notin W entonces no se cumple (b).

    WW no es un subespacio de R2×2{\mathbb{R}^{2 \times 2}}.

    178

    3.1.5   Resumen de los subespacios de R2{\mathbb{R}^2} y R3{\mathbb{R}^3}

    Después de estos ejemplos podemos resumir cuales son los diferentes tipos de subespacios de R2  {\mathbb{R}^2\;} y   R3{\;\mathbb{R}^3}:

    Subespacios de R2\mathbf{R^2}Subespacios de R3\mathbf{R^3}
    {(0; 0)}{(0; 0; 0)}
    Rectas que pasan por el origen.Rectas que pasan por el origen.
    Planos que pasan por el origen.
    R2{\mathbb{R}^2} (como subespacio de sí mismo.)R3{\mathbb{R}^3} (como subespacio de sí mismo.)

    No hay ninguna otra clase de subespacios en R2  {\mathbb{R}^2\;} y   R3{\;\mathbb{R}^3}.

    179

    3.2   Conjunto generador. L.I. L.D. Bases

    3.2.1   Combinación lineal

    Definición

    Sean v1,  v2,  ,  vr,  w\large{v_1},\;{v_2},\; \ldots ,\;{v_r},\;w vectores de un espacio vectorial VV. Se dice que el vector w  \large w\; es una combinación lineal de los vectores v1,  v2,    ,  vr  \large{v_1},\;{v_2},\; \ldots \;,\;{v_r}\; si se puede expresar como sigue:

    w=k1v1+k2v2+  +krvr  \large w = {k_1}{v_1} + {k_2}{v_2} + \; \ldots + {k_r}{v_r}\;

    donde k1,  k2,  ,  kr\large{k_1},\;{k_2},\; \ldots ,\;{k_r} son escalares.

    Observación: Nosotros estamos trabajando con espacios vectoriales reales, o sea que los escalares son números reales.

    180

    3.2.2   Conjunto generador

    Sea {  v1  ,  v2,    vr}\left\{ {\;{v_1}\;,\;{v_2},\; \ldots \;{v_r}} \right\} un conjunto de vectores de un espacio vectorial VV.

    Si todo vector de VV puede expresarse como combinación lineal de v1,  v2,    ,  vr{v_1},\;{v_2},\; \ldots \;,\;{v_r}, entonces se dice que {  v1,  v2,    ,  vr  }\left\{ {\;{v_1},\;{v_2},\; \ldots \;,\;{v_r}\;} \right\} es un conjunto generador de VV o también que v1,  v2,    ,  vr{v_1},\;{v_2},\; \ldots \;,\;{v_r} generan VV.

    3.2.3   Subespacio generado

    Sean v1,  v2,  ,  vr{v_1},\;{v_2},\; \ldots ,\;{v_r} vectores de un espacio vectorial VV.
    1. El vector nulo puede expresarse como combinación lineal de dichos vectores:

    0V  =  0  v1    +  0  v2  +    +  0vr{0_V}\; = \;0\;{v_1}\;\; + \;0\;{v_2}\; + \; \ldots \; + \;0{v_r}

    181

    2. Si sumamos dos combinaciones lineales de los vectores dados, obtenemos otra combinación lineal:

    (a1v1  +    +  arvr)+  (b1v1    +    +  brvr)=( {{a_1}{v_1}\; + \; \ldots \; + \;{a_r}{v_r}} ) + \;( {{b_1}{v_1}\;\; + \; \ldots \; + \;{b_r}{v_r}}) =

    =(a1+  b1)v1++(  ar+br)vr=( {{a_1} + \;{b_1}} ){v_1} + \ldots + ( {\;{a_r} + {b_r}} ){v_r}

    3. Si multiplicamos un escalar k por una combinación lineal de los vectores dados, obtenemos una combinación lineal de dichos vectores:

    k(  a1  v1  +    +  ar  vr)  =  (ka1)v1  +    +  (k  ar)vrk\left( {\;{a_1}\;{v_1}\; + \; \ldots \; + \;{a_r}\;{v_r}} \right)\; = \;\left( {k{a_1}} \right){v_1}\; + \; \ldots \; + \;\left( {k\;{a_r}} \right){v_r}

    Estas tres condiciones permiten afirmar que el conjunto de todas las combinaciones lineales de v1,  v2,  ,  vr{v_1},\;{v_2},\; \ldots ,\;{v_r} es un subespacio de VV.
    Entonces:

    Dados los vectores v1,  v2,  ,  vr{v_1},\;{v_2},\; \ldots ,\;{v_r} en VV, llamamos subespacio generado por v1,  v2,  ,  vr{v_1},\;{v_2},\; \ldots ,\;{v_r} al conjunto de todas las combinaciones lineales de estos vectores. Lo denotamos con la expresión gen{v1,  v2,  ,  vr}gen\left\{{{v_1},\;{v_2},\; \ldots ,\;{v_r}} \right\}
    gen{v1,v2,,  vr}={vV:v=  α1v1+α2v2+  +αrvr, gen\{{v_1},{v_2}, \ldots ,\;{v_r} \} = \{ v \in V: v = \;{\alpha _1}{v_1} + {\alpha _2}{v_2} + \; \ldots + {\alpha _r}{v_r}, con   αiR  }\; {\alpha _i} \in R\; \} subespacio de VV.

    182

    3.2.4   Independencia lineal y dependencia lineal

    En los ejemplos 1 y 2 de Conjunto generador vimos que los conjuntos {(1,1),  (1,1)}   \{(1,1),\;(1, - 1)\}\; y   {(1,1),  (1,1),  (2,0)}\; \{(1,1),\;(1, - 1),\;(2,0)\} generan R2{\mathbb{R}^2}. Si tuviéramos que elegir uno de ellos como generador de R2{\mathbb{R}^2}, ¿por cuál nos inclinaríamos?
    El problema de encontrar los conjuntos generadores más «pequeños» para un e. v., depende de la noción de independencia lineal, que presentamos en esta sección. Si A={v1,v2,,vr}A = \left\{ {{v_1},{v_2}, \ldots ,{v_r}} \right\} es un conjunto de vectores de un espacio vectorial VV, entonces:

    α1v1+α2  v2++αr  vr=0V{\alpha _1}{v_1} + {\alpha _2}\;{v_2} + \ldots + {\alpha _r}\;{v_r} = {0_V}

    tiene al menos la solución trivial:   α1=α2==αr=0{\;\alpha _1} = {\alpha _2} = \ldots = {\alpha _r} = 0
    Si ésta es la única solución , entonces se dice que AA es un conjunto linealmente independiente.
    Si hay otras soluciones (además de la trivial) entonces A es un conjunto linealmente dependiente.

    183

    Una forma alternativa de caracterizar la dependencia lineal es la siguiente:
    Un conjunto de vectores {  v1,v2,  ,vr  }\left\{ {\;{v_1},{v_2},\; \ldots ,{v_r}\;} \right\} de un espacio vectorial VV es linealmente dependiente si y sólo si al menos uno de los vectores puede expresarse como combinación lineal de los demás. {*}

    Demostración:

    \Rightarrow Si el conjunto {  v1,v2,  ,vr  }\left\{ {\;{v_1},{v_2},\; \ldots ,{v_r}\;} \right\} es linealmente dependiente, la ecuación

    α1v1+α2v2+  +  αrvr=0V{\alpha _1}{v_1} + {\alpha _2}{v_2} + \; \ldots + \;{\alpha _r}{v_r} = {0_V} admite otras soluciones además de la trivial. O sea, existe una combinación lineal donde al menos uno de los escalares es distinto de cero, que da el vector nulo.
    Supongamos que α10{\alpha _1} \ne 0. Entonces resulta:

    v1=(α2α1)v2++(αrα1)vr{v_1} = \left( { - \frac{{{\alpha _2}}}{{{\alpha _1}}}} \right){v_2} + \ldots + \left( { - \frac{{{\alpha _r}}}{{{\alpha _1}}}} \right){v_r}

    Por lo tanto, el vector v1{v_1} es combinación lineal de los demás.

    \Leftarrow Sabemos que uno de los vectores puede expresarse como combinación de los demás. Sin perder generalidad, supongamos que:

    v1=k2v2++krvr{\color{red}{v_1}} = {k_2}{v_2} + \ldots + {k_r}{v_r}

    184

    Entonces:

    1v1+k2v2++krvr=0V - 1{\color{red}{v_1}} + {k_2}{v_2} + \ldots + {k_r}{v_r} = {0_V}

    Existe una combinación lineal no trivial (al menos uno de los escalares es distinto de cero) que es igual al vector nulo. Por lo tanto, el conjunto {  v1,v2,  ,vr  }\left\{ {\;{v_1},{v_2},\; \ldots ,{v_r}\;} \right\} es linealmente dependiente, como queríamos demostrar.

    3.2.5   Propiedades

    Propiedad 1

    Un conjunto formado por un solo vector, ¿es linealmente independiente (LI) o dependiente (LD)?

    Planteamos la combinación lineal: α  v=0V\alpha \;v = {0_V}

    185

    Si v=0Vv = {0_V}, α\alpha puede tomar cualquier valor. Por lo tanto: {0V}\left\{ {{0_V}} \right\} es LD.

    Si v0Vv \ne {0_V} , la única solución es α=0\alpha = 0. Por lo tanto: {v}\left\{ v \right\} es LI.

    Propiedad 2

    Si un conjunto de vectores contiene al vector nulo, entonces es linealmente dependiente (LD).

    Demostración:
    Sea A={v1,  v2,  ,  vr,  0V}VA = \left\{ {{v_1},\;{v_2},\; \ldots ,\;{v_r},\;{0_V}} \right\} \subset V

    Entonces se tiene que:

    0v1+0v2++0vr+1  0V=0V0{v_1} + 0{v_2} + \ldots + 0{v_r} + 1\;{0_V} = {0_V}

    Probamos que existe una combinación lineal con escalares no todos nulos, que da el vector nulo. Por lo tanto, AA es linealmente dependiente.

    3.2.6   Interpretación geométrica

    Observación previa: En lo que sigue consideramos los vectores colocados a partir del origen de coordenadas.

    De {*} resulta que dos vectores v1{v_1} y v2{v_2} son linealmente dependientes (LD) si y sólo si uno de ellos es un múltiplo escalar del otro.

    186

    Podemos afirmar entonces que dos vectores en R2{\mathbb{R}^2} o R3{\mathbb{R}^3} son LD si y sólo si están sobre la misma recta que pasa por el origen (vectores paralelos).

    187

    En R3{\mathbb{R}^3}, tres vectores v1,  v2,  v3{v_1},\;{v_2},\;{v_3} son LD si y sólo si están situados en el mismo plano que pasa por el origen (vectores coplanares).


    3.2.7   Base y dimensión de un espacio vectorial

    Habíamos visto que los conjuntos B={(1,  1),  (1,  1)}B = \left\{ {\left( {1,\;1} \right),\;\left( {1,\; - 1} \right)} \right\} y C={(1,  1),  (1,  1),  (2,  0)}C = \left\{ {\left( {1,\;1} \right),\;\left( {1,\; - 1} \right),\;\left( {2,\;0} \right)} \right\} generan R2.{\mathbb{R}^2}. ¿Cuál es la diferencia entre ellos?

    188

    BB es un conjunto linealmente independiente, en cambio CC es linealmente dependiente porque (2,  0)=(1,  1)+(1,  1)\left( {2,\;0} \right) = \left( {1,\;1} \right) + \left( {1,\; -1} \right).

    Un conjunto de vectores B  ={v1  ,  v2  ,  ,  vn}B\; = \left\{ {{v_1}\;,\;{v_2}\;,\; \ldots ,\;{v_n}} \right\} de un espacio vectorial VV se denomina base de V  V\;si y sólo si:

    i. BB es linealmente independiente;
    ii. BB genera a VV.

    Teniendo en cuenta esta definición, B={(1,  1),  (1,  1)}B = \left\{ {\left( {1,\;1} \right),\;\left( {1,\; - 1} \right)} \right\} es una base de R2{\mathbb{R}^2}. Otra base de R2{\mathbb{R}^2} muy usual es la que contiene a los versores canónicos: E={(1,  0)  ,  (0,  1)}E = \left\{ {\left( {1,\;0} \right)\;,\;\left( {0,\;1} \right)} \right\}.

    Puede demostrarse que:

    Si B  =  {v1,  v2,  ,  vn}B\; = \;\left\{ {{v_1},\;{v_2},\; \ldots ,\;{v_n}} \right\} es una base de VV, entonces todo conjunto con más de n vectores es linealmente dependiente.

    De acuerdo con esta propiedad, podemos deducir una característica común a toda base de un espacio vectorial:

    Sean B={v1,v2,  ,  vn}B = \left\{ {{v_1},{v_2},\; \ldots ,\;{v_n}} \right\} y B={w1,  w2,  ,  wq}B' = \left\{ {{w_1},\;{w_2},\; \ldots ,\;{w_q}} \right\} dos bases del espacio vectorial VV.

    Como BB es una base, todo conjunto de más de nn vectores es LD. Pero BB' es LI, entonces:
    qn          [1]q \le n\;\;\;\;\;\left[ 1 \right]

    189

    Como BB' es una base, todo conjunto de más de qq vectores es LD. Pero BB es LI, entonces:

    nq          [2]n \le q\;\;\;\;\;\left[ 2 \right]

    De [1] y [2] se deduce que q=nq = n.

    En consecuencia:

    Si B={v1,v2,,  vn}B = \left\{ {{v_1},{v_2}, \ldots ,\;{v_n}} \right\} es una base de VV, cualquier otra base de VV tiene nn vectores. Esto permite definir el concepto de dimensión.

    La dimensión de un espacio vectorial VV es la cantidad de vectores que componen una base de VV.

    Si B={v1,v2,,  vn}B = \left\{ {{v_1},{v_2}, \ldots ,\;{v_n}} \right\} es una base de VV, la dimensión de VV es nn y lo indicamos como dim(V)  =  ndim\left( V \right)\; = \;n.

    Si no existe una base de VV formada por un conjunto finito de vectores, se dice que VV es un espacio de dimensión infinita. Un ejemplo es el espacio de todos los polinomios (de cualquier grado).

    Como el vector nulo es linealmente dependiente, el espacio {0V}\left\{ {{0_V}} \right\} no tiene base. A este espacio compuesto únicamente por el vector nulo, se le asigna dimensión cero:

    dim(  {0V}  )=0\dim \left( {\;\left\{ {{0_V}} \right\}\;} \right) = 0

    190

    Para determinar la dimensión de un espacio vectorial, es suficiente hallar una base de dicho espacio. Veamos qué dimensión tienen los espacios vectoriales con los cuales trabajaremos:

    En R2{\mathbb{R}^2} conocemos la base canónica:

    E2={(1,0),(0,1)}  dim(R2)=2{E_2} = \left\{ {\left( {1,0} \right),\left( {0,1} \right)} \right\}\; \Rightarrow \dim \left( {{\mathbb{R}^2}} \right) = 2

    En R3{\mathbb{R}^3} la base canónica es:

    E3={(1,0,0),(0,1,0),(0,0,1)}dim(R3)=3{E_3} = \left\{ {\left( {1,0,0} \right),\left( {0,1,0} \right),\left( {0,0,1} \right)} \right\} \Rightarrow \dim \left( {{\mathbb{R}^3}} \right) = 3

    Análogamente en R4{\mathbb{R}^4}:

    E4={(1,0,0,0),(0,1,0,0),(0,0,1,0),(0,0,0,1)}{E_4} = \left\{ {\left( {1,0,0,0} \right),\left( {0,1,0,0} \right),\left( {0,0,1,0} \right),\left( {0,0,0,1} \right)} \right\} \Rightarrow dim(R4)=4 \Rightarrow \dim \left( {{\mathbb{R}^4}} \right) = 4

    De acuerdo con el número de vectores que componen cada una de estas bases, podemos afirmar que:

    dim(Rn)=n\dim \left( {{\mathbb{R}^n}} \right) = n

    Para determinar la dimensión de los espacios de matrices, consideremos por ejemplo V=R3×2V = {\mathbb{R}^{3 \times 2}}.
    Cualquier matriz de 3×23 \times 2 puede expresarse como sigue:

    191
    (abcdef)=a(100000)+b(010000)+c(001000)+\begin{pmatrix} a & b \\ c & d \\ e & f\end{pmatrix} = a\begin{pmatrix} 1 & 0 \\ 0 & 0 \\ 0 & 0\end{pmatrix}+ b\begin{pmatrix} 0 & 1 \\ 0 & 0 \\ 0 & 0\end{pmatrix} + c\begin{pmatrix} 0 & 0 \\ 1 & 0 \\ 0 & 0\end{pmatrix} + +  d(000100)+e(000010)+f(000001) +\; d\begin{pmatrix} 0 & 0 \\ 0 & 1 \\ 0 & 0\end{pmatrix} + e\begin{pmatrix} 0 & 0 \\ 0 & 0 \\ 1 & 0\end{pmatrix} + f\begin{pmatrix} 0 & 0 \\ 0 & 0 \\ 0 & 1\end{pmatrix}

    Observemos la similitud con R6{\mathbb{R}^6} , sólo cambia el formato.

    Encontramos seis matrices linealmente independientes que generan V=R3×2V = {\mathbb{R}^{3 \times 2}}. Es decir, encontramos una base (llamada base canónica) de este espacio y por lo tanto:
    dim(R3×2)=3×2=6\dim \left( {{\mathbb{R}^{3×2}}} \right) = 3×2 = 6.

    Generalizando, podemos afirmar que:

    dim(Rmxn)=mxn\dim \left( {{\mathbb{R}^{mxn}}} \right) = mxn

    Busquemos la dimensión de los espacios de polinomios. Consideremos por ejemplo V=P2  V = {P_{2\;}}. Cualquier polinomio de P2{P_2} puede expresarse cómo sigue:

      p(x)=a0+a1x+a2x2=a01+a1x+a2x2\;p\left( x \right) = {a_0} + {a_1}x + {a_2}{x^2} = {a_0}\color{red}{1} \color{black} + {a_1} \color{red}{x}\color{black} + {a_2}{\color{red}{x^2}}

    El conjunto {1,x,x2}\left\{ {1,x,{x^2}} \right\} genera P2{P_2} y además es linealmente independiente. Hemos obtenido una base (llamada canónica) de P2{P_2}, y por lo tanto dim(P2)=3\dim \left( {{P_2}} \right) = 3.

    192

    Análogamente:

    {1,x,x2,x3}    base  canoˊnica  de  P3\left\{ {1,x,{x^2},{x^3}} \right\}\;\;base\;canó nica\;de\;{P_3}

    {1,x,x2,x3,,  xn}  base  canoˊnica  de  Pn\left\{ {1,x,{x^2},{x^3}, \ldots ,\;{x^n}} \right\}\;base\;canó nica\;de\;{P_n}

    Entonces:

    dim(Pn)=n+1\dim \left( {{P_n}} \right) = n + 1


    3.2.8   Propiedades relacionadas con la dimensión

    Si dim(V)=n,\dim \left( V \right) = n,\, puede afirmarse que:

    1) Todo conjunto de nn vectores linealmente independientes en VV es una base.

    193

    2) Todo conjunto de nn vectores que genere VV es una base.

    3) Todo conjunto de más de nn vectores en el espacio vectorial VV es linealmente dependiente.

    4) Todo conjunto linealmente independiente en VV puede extenderse a una base.

    3.2.9   Coordenadas de un vector respecto de una base

    Propiedad: Si B={v1,v2,,vn}B = \left\{ {{v_1},{v_2}, \ldots ,{v_n}} \right\} es una base del espacio vectorial VV, todo vector de V  V\;puede expresarse de forma única como combinación lineal de los vectores de BB.

    194

    Esta propiedad permite definir coordenadas de un vector respecto de una base.

    Sea B={v1,v2,,vn}B = \left\{ {{v_1},{v_2}, \ldots ,{v_n}} \right\} base de VV.

    Para cada uVu \in V, existen únicos escalares α1,α2,,αnR{\alpha _1},{\alpha _2}, \ldots ,{\alpha _n} \in \mathbb{R} tales que: u=α1v1+α2v2++αnvnu = {\alpha _1}{v_1} + {\alpha _2}{v_2} + \ldots + {\alpha _n}{v_n}

    Estos escalares se denominan coordenadas del vector uu respecto de la base BB.

    Indicaremos las coordenadas mediante la siguiente notación:

    [u]B=(α1α2αn) [u]_{B}= \begin{pmatrix} \alpha_{1} \\ \alpha_{2} \\ \vdots \\ \alpha_{n} \end{pmatrix}


    Ejercicios resueltos
    195

    3.2.10   Base y dimensión de un subespacio vectorial

    Recordemos que un subespacio es un espacio vectorial en sí mismo, por lo tanto podemos hallar una base y su dimensión.

    Si S es un subespacio de V, entonces: dim(S)dim(V){\rm{dim}}\left( S \right) \le {\rm{dim}}\left( V \right).

    Veamos cuáles son las dimensiones de los distintos tipos de subespacios de R3{\mathbb{R}^3}:

    • S={(0,0,0)}S = \left\{ {\left( {0,0,0} \right)} \right\} no tiene base y como habíamos dicho, se le asigna dimensión 0.

      dim({0V})  =0\dim \left( {\left\{ {{0_V}} \right\}} \right)\; = 0

    • Rectas que pasan por el origen,
    • Planos que pasan por el origen y
    • R3{\mathbb{R}^3}\, (Sabemos que R3{\mathbb{R}^3}\, tiene dimensión 3.)

    Consideremos un plano que pase por el origen, por ejemplo:

    π:    x+3y2z=0\pi :\;\;x + 3y - 2z = 0

    x=3y+2zx = - 3y + 2z

    (x,y,z)π(x,y,z)=(3y+2z,y,z)=\left( {x,y,z} \right) \in \pi \Leftrightarrow \left( {x,y,z} \right) = \left( {-3y + 2z,y,z} \right) = =(3y,y,0)+(2z,0,z)=y(3,1,0)+z(2,0,1) = \left( { - 3y,y,0} \right) + \left( {2z,0,z} \right) = y\left( { -3,1,0} \right) + z\left( {2,0,1} \right)

    196

    Esto quiere decir que cualquier vector en ese plano se puede escribir como combinación lineal de (3,1,0)\left( { - 3,1,0} \right) y (2,0,1)\left( {2,0,1} \right). Cómo son LI:

    {(3,1,0),(2,0,1)}  es  una  base  de  S1\left\{ {\left( { - 3,1,0} \right),\left( {2,0,1} \right)} \right\}\;es\;una\;base\;de\;{S_1}

    dim(S1)=2\dim \left( {{S_1}} \right) = 2

    Planos que pasan por el origen son subespacios de dimensión 2.

    Ahora consideremos el subespacio:

    S2={(x,y,z)R3    x+y=0      ,        xyz=0}{S_2} = \left\{ {\left( {x,y,z} \right) \in {\mathbb{R}^3}|\;\;x + y = 0\;\;\;,\;\;\;\;x - y - z = 0} \right\}

    r:{π1:x+y=0π2:xyz=0;    r=π1π2 r: \begin{cases} \pi_{1}: x + y = 0 \\ \pi_{2}: x - y - z = 0 \end{cases}; \;\; r = \pi_{1} \cap \pi_{2}

    La intersección de dos planos no paralelos es una recta. ¿Cómo podemos encontrar una base de una recta?

    y=x  y = - x\;

          x(x)z=0      z=2x\; \Rightarrow \;\;x -\left( { - x} \right) - z = 0\; \Rightarrow \;\;z = 2x

    Si llamamos x=tx = t , resulta:

    (x,y,z)=(t,t,2t)=t(1,1,2)\left( {x,y,z} \right) = \left( {t, - t,2t} \right) = t\left( {1, - 1,2} \right)

    197

    Observamos que todos los vectores de la recta pueden expresarse como combinación lineal del vector director (1,1,2)\left( {1, - 1,2} \right), que además es LI. Por lo tanto, {(1,1,2)}\left\{ {\left( {1, - 1,2} \right)} \right\} es una base de este subespacio.

    Rectas que pasan por el origen son subespacios de dimensión 1.

    En los ejemplos anteriores observamos cómo disminuye la dimensión de un subespacio a medida que agregamos ecuaciones, tal como se muestra en el siguiente cuadro:

    Número mínimo de ecuaciones que definen SSDimensión del subespacioObjeto geométricoEjemplo
    03 R3\mathbf{R^3}
    12 Planos por el origen z=0 \displaystyle z = 0
    21 Rectas por el origen {y=0z=0 \begin{cases} y = 0 \\ z = 0 \end{cases}
    30 {(0;0;0)} \Set{(0; 0; 0)} {x=0y=0z=0 \begin{cases} x = 0 \\ y = 0 \\ z = 0 \end{cases}
    198

      Ejercicios y problemas

    1) Indicar si A={(1,1,1),  (1,3,2)}A = \left\{ {\left( {1,1,1} \right),\;\left( {1,3,2} \right)} \right\} es una base de
    W={(x,y,z)R3:  x+y2z=0}W = \left\{ {\left( {x,y,z} \right) \in {R^3}:\; x + y - 2z = 0} \right\}, justificando la respuesta.
    Posible estrategia de resolución: hacer el producto vectorial de los generadores de AA y verificar que el resultado es proporcional al vector normal al plano WW.

    2) Considere el subespacio de R4R^4

    T1={(x1,x2,x3,x4)R4  x1+x4=0    x2x4=0}{T_1} = \left\{ {\left( {{x_1},{x_2},{x_3},{x_4}} \right) \in {R^4}|\;{x_1} + {x_4} = 0\; \wedge \;{x_2} - {x_4} = 0} \right\}

    Halle base y dimensión.
    Teniendo en cuenta las ecuaciones que definen T1{T_1}, ¿podría anticipar su dimensión?

    3) Hallar una base y la dimensión del subespacio de matrices simétricas de 2×2:

    T2={AR2×2      A=At}{T_2} = \left\{ {A \in {R^{2 \times 2}}\;|\;\;A = {A^t}} \right\}

    ¿Cómo es una matriz simétrica de 2×2?

    A=(abbd) A = \begin{pmatrix} a & b \\ b & d \end{pmatrix}

    199

    (abbd)=a(1000)+b(0110)+d(0001)\begin{pmatrix} a & b \\ b & d \end{pmatrix} = a\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}+ b\begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} + d\begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix}

    Las tres matrices halladas son L.I. y generan el subespacio de matrices simétricas, por lo tanto hemos encontrado una base de dicho subespacio.

    4) Hallar base y dimensión de

    S={p(x)P2  a02a1+3a2=0}S = \left\{ {p\left( x \right) \in {P_2}|\;{a_0} - 2{a_1} + 3{a_2} = 0} \right\}

    5) Hallar una base y la dimensión del subespacio de matrices antisimétricas de 2×2 y del de 3×3.

    200

    3.2.11   Bases de subespacios definidos por generadores

    Hasta ahora hemos buscado bases de subespacios definidos por ecuaciones. ¿Qué ocurre cuando el subespacio está definido por sus generadores?

    Veamos el siguiente ejemplo:

    Hallar una base y la dimensión de
    S = gen{(1,1,2),(1,1,0),  (0,1,1)}R3 gen\left\{ {\left( {1,1,2} \right),\left( {1, -1,0} \right),\;\left( {0,1,1} \right)} \right\} \subset {R^3}.

    En este caso, por la definición de S sabemos que {(1,1,2),(1,1,0),  (0,1,1)}\left\{ {\left( {1,1,2} \right),\left( {1, -1,0} \right),\;\left( {0,1,1} \right)} \right\} es un conjunto generador de S. Para determinar si es una base, tendremos que analizar la independencia lineal:

    - Si son vectores LI, entonces son base del subespacio.

    - Si son LD, tendremos que extraer una base eliminando los vectores «que sobren».

    En el caso específico de 3 vectores en R3, podemos utilizar el determinante (que es el producto mixto). Como en este caso el determinante da cero, los vectores son coplanares y por lo tanto L.D. Tenemos que extraer una base eliminando alguno de los vectores, por ejemplo:

    • {(1,1,2),(1,1,0)}\left\{ {\left( {1,1,2} \right),\left( {1, - 1,0} \right)} \right\} es una base de S
    • {(1,1,0),  (0,1,1)}\left\{ {\left( {1, - 1,0} \right),\;\left( {0,1,1} \right)} \right\} es otra base de S.
    201

    Podemos afirmar que la dimensión de S es 2.

    Una forma práctica de extraer bases es armar una matriz con los vectores dados y llevarla a la forma escalonada:

    Una matriz es escalonada (por filas) si satisface las siguientes propiedades:

    1. Las filas nulas (todos sus elementos son ceros) se encuentran en la parte inferior.

    2. En cada fila no nula, el primer elemento distinto de cero (pivote) está a la derecha del pivote de la fila superior.

    Una matriz cualquiera puede llevarse a la forma escalonada aplicando operaciones elementales entre sus filas. Por ejemplo consideremos la matriz que armamos con los generadores de SS:

    Se denomina rango de una matriz al número de filas LI que tiene la matriz. Veremos en el siguiente capítulo la importancia de este concepto en el estudio de sistemas de ecuaciones lineales.

    Puede demostrarse que:

    1. Si se realizan operaciones elementales entre las filas de una matriz, el rango se conserva.

    2. Las filas no nulas de una matriz escalonada son LI.

    202

    Por lo tanto, para determinar el rango de una matriz se aplican operaciones elementales para obtener una matriz escalonada y se cuentan las filas no nulas.

    En el ejemplo, la matriz escalonada tiene rango 2, por lo tanto la matriz que armamos con los generadores de S tiene rango 2. Esto significa que de los tres generadores de S hay dos linealmente independientes. ¿Cuál es entonces la dimensión de S?
    dim(S)=2  \dim \left( S \right) = 2\;

    Este método también permite obtener bases: las filas no nulas de la última matriz son otra base de S, ya que fueron obtenidas como combinaciones lineales de los vectores de S:

    {(1,1,2),  (0,2,2)}  otra  base  de  S\left\{ {\left( {1,1,2} \right),\;\left( {0,2,2} \right)} \right\}\;otra\;base\;de\;S

      Ejemplo

    Dado el conjunto:

    A={1+x;1x2;2+3x+kx2}P2A = \left\{ {1 + x;1 -{x^2};2 + 3x + k{x^2}} \right\} \subset {P_2}

    Hallar todos los valores de kk para que AA genere un subespacio de dimensión 2

    Para el kk hallado encontrar las ecuaciones del subespacio generado por AA.

    203
    Resolución

    Sabemos que dim(P2)=3\dim \left( {{P_2}} \right) = 3 , entonces todo conjunto de 3 vectores LI en P2{P_2} es base de P2{P_2}. Como se pide que la dimensión del subespacio sea 2, debemos hallar kk de modo que los vectores sean LD.

    Los dos primeros vectores de AA son LI, entonces se trata de analizar para qué valores de kk el tercer vector es combinación lineal de los anteriores:

    2+3x+kx2=α(1+x)+β(1x2)2 + 3x + k{x^2} = \alpha \left( {1 + x} \right) + \beta \left( {1 - {x^2}} \right)

    {2=α+β3=αk=βα=3,β=1,k=1\begin{cases} 2 = \alpha + \beta \\ 3 = \alpha \\ k = -\beta \end{cases} \Rightarrow \alpha = 3,\: \beta = -1,\: k = 1

    Para k=1k = 1, el polinomio 2+3x+kx22 + 3x + k{x^2} es combinación lineal de 1+x  1 + x \; y   1x2.   \;1 - {x^2}.\; Entonces:

    gen{1+x;1x2;2+3x+x2}=gen{1+x;1x2}gen\left\{ {1 + x;1 - {x^2};2 + 3x + {x^2}} \right\} = gen\left\{ {1 + x;1 - {x^2}} \right\}

    {1+x;1x2}  es  base  de  P2\left\{ {1 + x;1 - {x^2}} \right\}\;es\;base\;de\;{P_2}

    Ahora busquemos la ecuación del subespacio. Tomamos un polinomio genérico, p(x)=a0+a1x+a2x2p\left( x \right) = {a_0} + {a_1}x + {a_2}{x^2} y lo escribimos como combinación lineal de los vectores de la base:

    a0+a1x+a2x2=α(1+x)+β(1x2){a_0} + {a_1}x + {a_2}{x^2} = \alpha \left( {1 + x} \right) + \beta \left( {1 - {x^2}} \right)

    204

    {α+β=a0β=a2α=a1α=a1,β=a2,  a0=a1a2 \begin{cases} \alpha + \beta = {a_0} \\ - \beta = {a_2} \\ \alpha = {a_1} \end{cases} \large \Rightarrow \alpha = {a_1},\: \beta = - {a_2},\: \Rightarrow \;{a_0} = {a_1} - {a_2}

    Entonces esa es la ecuación que define al subespacio:

    S=gen(A)={pP2        a0a1+a2=0}S = gen\left( A \right) = \left\{ {p \in {P_2}\;\;|\;\;{a_0} - {a_1} + {a_2} = 0} \right\}

    Como verificación, puede comprobarse que los dos vectores de la base verifican la ecuación obtenida.


    Video con ejercicio resuelto

    Video 3.1. Ejercicio Álgebra CBC A-62 (video de unamunoenlinea en YouTube, Licencia Atribución de Creative Commons)

    205

    3.3   Operaciones con subespacios

    3.3.1   Intersección

    Sean SS y TT subespacios del mismo espacio vectorial VV. Definimos la intersección como sigue:

    ST={vV  :  vS        vT}  S \cap T = \left\{ {v \in V\;:\;v \in S\;\; \wedge \;\;v \in T} \right\}\;

    Propiedad

    STS \cap T es subespacio de VV

      Ejemplo 1

    Sean los subespacios de R3{\mathbb{R}^3}:

    S={(x,y,z)  x3z=0}S = \left\{ {\left( {x,y,z} \right)|\;x - 3z = 0} \right\}

    T={(x,y,z)  x+yz=0}T = \left\{ {\left( {x,y,z} \right)|\;x + y - z = 0} \right\}

    Hallar STS \cap T.

    Resolución

    Por definición STS \cap T es un conjunto que estará formado por los vectores que pertenezcan a SS y a TT. Es decir aquellos vectores que satisfagan las ecuaciones de SS y las de TT:

    ST={(x,y,z)  x3z=0          x+yz=0}S \cap T = \left\{ {\left( {x,y,z} \right)|\;x - 3z = 0\;\; \wedge \;\;\;x + y - z = 0} \right\}

    206

    Se trata de una recta definida como intersección de dos planos. Una base de la recta es un vector director.

    Geométricamente podemos buscar el vector director como el producto vectorial de los vectores normales a los planos:

    v=n1×n2=ijk103111=(3,2,1) v = n_1 \times n_2 = \begin{vmatrix} i & j & k \\ 1 & 0 & -3 \\ 1 & 1& -1 \end{vmatrix} = \left( {3, - 2,1} \right)

    ST={(x,y,z)=λ(3,2,1),  λR}S \cap T = \left\{ {\left( {x,y,z} \right) = \lambda \left( {3, - 2,1} \right),\;\lambda \in \mathbb{R}} \right\}

    Entonces {(3,2,1)}\left\{ {\left( {3, - 2,1} \right)} \right\} es una base de STS \cap T.

    Otra forma de resolverlo es buscar la solución del sistema de ecuaciones:

    {x3z=0x+yz=0  {x=3zy=zx=2z   \begin{cases} x - 3z = 0 \\ x + y - z = 0 \end{cases} \; \Rightarrow \begin{cases} x = 3z \\ y = z -x = -2z \end{cases} \; \Rightarrow (3z,2z,z),  zR. \Rightarrow {(3z, - 2z, z)}, \; \forall z \in \mathbb{R}.

    (x,y,z)=z(3,2,1)\left( {x,y,z} \right) = z\left( {3, - 2,1} \right)

    Y entonces otra vez llegamos a que {(3,2,1)}\left\{ {\left( {3, - 2,1} \right)} \right\} es una base de STS \cap T.

    207

      Ejemplo 2

    Sean los subespacios de R2×2{\mathbb{R}^{2 \times 2}}:

    S={AR2×2      A=At}S = \left\{ {A \in {\mathbb{R}^{2 \times 2}}\;|\;\;A = {A^t}} \right\}

    T=gen  {(1021),  (1010)}T = gen\; \bigg \{\begin{pmatrix} 1 & 0 \\ 2 & -1 \end{pmatrix}, \;\begin{pmatrix} 1 & 0 \\ 1 & 0 \end{pmatrix}\bigg\} Hallar STS \cap T.

    Video con ejercicio resuelto

    Video 3.2. Ejercicio Álgebra CBC A-62 (video de unamunoenlinea en YouTube, Licencia Atribución de Creative Commons)

    208

    3.3.2   Suma de subespacios

    Dados S  y  T  S\;y\;T\; subespacios de VV, se define la suma como sigue:

    S+T={vV      v=v1+v2  ,      con    v1S,    v2T}S + T = \left\{ {v \in V\;|\;\;v = {v_1} + {v_2}\;,\;\;\;con\;\;{v_1} \in S,\;\;{v_2} \in T} \right\} suma de subespacios

    Propiedad: S+TS + T es un subespacio del espacio vectorial VV.

    Si conocemos conjuntos generadores de SS y de TT, podemos hallar generadores de la suma:

    S=gen{v1,v2,,vq}S = gen\left\{ {{v_1},{v_2}, \ldots ,{v_q}} \right\} y T=gen{w1,w2,,wr}T = gen\left\{ {{w_1},{w_2}, \ldots ,{w_r}} \right\}

    S+T=gen{v1,v2,vq,w1,w2,,wr} \Rightarrow S + T = gen\left\{ {{v_1},{v_2}, \ldots {v_q},{w_1},{w_2}, \ldots ,{w_r}} \right\}

    Para hallar la suma es usual buscar las bases de SS y TT. Como las bases son conjuntos generadores LI, si conocemos una base de cada subespacio podremos obtener un conjunto generador de la suma:

    Dadas las bases:

    BS={v1,v2,,vq}{B_S} = \left\{ {{v_1},{v_2}, \ldots ,{v_q}} \right\} y BT={w1,w2,,wr}{B_T} = \left\{ {{w_1},{w_2}, \ldots ,{w_r}} \right\}

    Resulta:

    {v1,,vq,w1,,wr}\left\{ {{v_1}, \ldots ,{v_q},{w_1}, \ldots ,{w_r}} \right\} conjunto generador de la suma.

    Observación: Se obtiene así un conjunto generador de la suma pero no siempre es linealmente independiente.

    209
    • Si es LI, encontramos una base de la suma.
    • Si es LD, podemos extraer una base de la suma eliminando los vectores «que sobran».

    3.3.3   Suma directa

    La suma de dos subespacios es directa si y sólo si la intersección de los subespacios es el vector nulo.

    S+T  es  directa    ST={0V}S + T\;es\;directa\; \Leftrightarrow \;S \cap T = \left\{ {{0_V}} \right\}

    Cuando la suma es directa se escribe:

    STS \oplus T


    210

    3.3.4   Teorema de la dimensión de la suma

    Si S1{S_1} y S2{S_2} son subespacios de un espacio vectorial VV (de dimensión finita), entonces:

    dim(S1+S2)=dim(S1)+dim(S2)dim(S1S2)\dim \left( {{S_1} + {S_2}} \right) = \dim \left( {{S_1}} \right) + \dim \left( {{S_2}} \right) - dim\left( {{S_1} \cap {S_2}} \right)

    En el caso particular de que la suma sea directa, como
    S1S2={0V}{S_1} \cap {S_2} = \left\{ {{0_V}} \right\} , resulta:

    dim(S1S2)=dim(S1)+dim(S2)\dim \left( {{S_1} \oplus {S_2}} \right) = \dim \left( {{S_1}} \right) + \dim \left( {{S_2}} \right)

      Ejemplo

    Dados los subespacios de P2{P_2}:

    S1={pP2      p(0)=0}{S_1} = \left\{ {p \in {P_2}\;|\;\;p\left( 0 \right) = 0} \right\}

    S2={pP2      p(1)=0}{S_2} = \left\{ {p \in {P_2}\;|\;\;p\left( 1 \right) = 0} \right\}

    Hallar bases de ambos subespacios y de la intersección.

    Resolución

    Hallemos una base de S1{S_1}:

    p(0)=a0+a10+a202=0a0=0p\left( 0 \right) = {a_0} + {a_1}0 + {a_2}{0^2} = 0 \Rightarrow {a_0} = 0

    211

    Entonces son los polinomios de la forma:

    a1x+a2x2{a_1}x + {a_2}{x^2}

    Luego una base de S1{S_1} es:

    BS1={x,x2}dim(S1)=2{B_{{S_1}}} = \left\{ {x,{x^2}} \right\} \Rightarrow \dim \left( {{S_1}} \right) = 2

    Hallemos una base de S2{S_2}:

    p(1)=a0+a11+a212=0a0+a1+a2=0p\left( 1 \right) = {a_0} + {a_1}1 + {a_2}{1^2} = 0 \Rightarrow {a_0} + {a_1} + {a_2} = 0

    Entonces son los polinomios de la forma:

    (a1a2)+a1x+a2x2=a1(1+x)+a2(1+x2)\left( { -{a_1} - {a_2}} \right) + {a_1}x + {a_2}{x^2} = {a_1}\left( { - 1 + x} \right) + {a_2}\left( { - 1 + {x^2}} \right)

    BS2={1+x    ,      1+x2}    dim(S2)=2{B_{{S_2}}} = \left\{ { - 1 + x\;\;,\;\;\; - 1 + {x^2}} \right\}\; \Rightarrow \;\dim \left( {{S_2}} \right) = 2

    Para buscar S1S2{S_1} \cap {S_2} debemos plantear que se cumplan las ecuaciones de S1{S_1} y también las de S2{S_2}:

    a0=0      a0+a1+a2=0    a0=0        a1=a2{a_0} = 0\; \wedge \;\;{a_0} + {a_1} + {a_2} = 0 \Rightarrow \;\;{a_0} = 0\;\; \wedge \;\;{a_1} = - {a_2}

    Los polinomios serán de la forma:

    a1xa1x2=a1(x2x){a_1}x -{a_1}{x^2} = -{a_1}\left( {{x^2} - x} \right)

    Luego:

    BS1S2={x2x}    dim(S1S2)=1{B_{{S_1} \cap {S_2}}} = \left\{ {{x^2} - x} \right\}\; \Rightarrow \;\dim \left( {{S_1} \cap {S_2}} \right) = 1

    212

    Nótese que como conocemos las dimensiones de S1{S_1}, S2{S_2} y S1S2{S_1} \cap {S_2}, podemos calcular la dimensión de S1+S2{S_1} + {S_2}:

    dim(S1+S2)=2+21=3\dim \left( {{S_1} + {S_2}} \right) = 2 + 2 - 1 = 3

    Pero el único subespacio de P2{P_2} con dimensión 3 es P2{P_2}. Luego: S1+S2=P2{S_1} + {S_2} = {P_2}.

      Ejercicio para el lector

    Dados los subespacios de R2×2{\mathbb{R}^{2 \times 2}}:

    W1:{AR2×2        A=At}{W_1}:\left\{ {A \in {\mathbb{R}^{2 \times 2}}\;\;|\;\;A = {A^t}} \right\}

    W2:{AR2×2      A=At}{W_2}:\left\{ {A \in {\mathbb{R}^{2 \times 2}}\;|\;\;A = - {A^t}} \right\}

    a) Hallar bases de W1{W_1} y W2{W_2}

    b) Obtener W1W2{W_1} \cap {W_2}.

    c) Sin hallar W1+W2{W_1} + {W_2} analizar la validez de la siguiente afirmación:

    W1W2=R2×2{W_1} \oplus {W_2} = {\mathbb{R}^{2 \times 2}}

    d) Proponer una base de R2×2  {\mathbb{R}^{2 \times 2}}\; formada por matrices simétricas y antisimétricas.

    213

    3.3.5   Producto interno

    En el primer capítulo vimos producto escalar entre vectores y sus aplicaciones a la Geometría. En esta sección nos proponemos generalizar esta operación a otros espacios vectoriales definiendo la noción general de producto interno a partir de las propiedades del producto escalar.

    Definición: Un producto interno en un espacio vectorial real VV es una operación que asigna a cada par de vectores uu y vv de VV un número real u.vu.v tal que se verifican las siguientes propiedades u,  v,  wV \forall {u,\;v,\;w} \in V y   αR   \;\forall \alpha \in \mathbb{R}\; :

    1. uv=vuu \cdot v = v \cdot u
    2. u(v+w)=(uv)+(uw)u \cdot \left( {v + w} \right) = \left( {u \cdot v} \right) + \left( {u \cdot w} \right)
    3. αuv=α(uv)\alpha u \cdot v = \alpha \left( {u \cdot v} \right)
    4. u.u0            u.u=0                u=0Vu.u \ge 0\;\;\; \wedge \;\;\;u.u = 0\;\;\;\; \Leftrightarrow \;\;\;\;u = {0_V}

    Es posible definir así distintos productos internos en cualquier espacio vectorial (mientras se verifiquen estas propiedades). En este libro sólo trabajaremos con el producto interno canónico en Rn,  {\mathbb{R}^n},\; que es la extensión del producto escalar:

    (x1,x2,,xn).(y1,y2,,yn)=x1y1+x2y2+xnyn\left( {{x_1},{x_2}, \ldots ,{x_n}} \right).\left( {{y_1},{y_2}, \ldots ,{y_n}} \right) = {x_1}{y_1} + {x_2}{y_2} + \ldots {x_n}{y_n}

    producto interno canónico en Rn\mathbb{R}^n

    214

    Esta definición nos permite extender el concepto de ortogonalidad a Rn{\mathbb{R}^n}:

    uv          u.v=0u \bot v\;\;\; \Leftrightarrow \;\;u.v = 0\: condición de ortogonalidad

      Ejemplo

    Realicemos el producto interno de los vectores de R4{\mathbb{R}^4}:

    u=(1,2,3,4)                  v=(1,0,1,1)u = \left( {1,2,3,4} \right)\;\;\;\;\;\;\;\;\;v = \left( {1,0,1, - 1} \right)

    u.v=1.1+2.0+3.1+4.(1)=0u.v = 1.1 + 2.0 + 3.1 + 4.\left( { - 1} \right) = 0

    Como u.v=0u.v = 0 entonces uu y vv son ortogonales.  \;


    3.3.6   Complemento ortogonal de un subespacio

    Sea SS subespacio de un espacio vectorial VV con producto interno.

    El complemento ortogonal de SS, que denotamos como S{S^ \bot}, es el conjunto de vectores de VV que son ortogonales a cada uno de los vectores de SS:

    S={  vV    vw=0        wS}{S^ \bot } = \left\{ {\;v \in V|\;\;v \cdot w = 0\;\;\;\;\forall w \in S} \right\}
    complemento ortogonal de SS

    Propiedad:  S{S^ \bot } es un subespacio de VV.

    215

    3.3.7   Propiedades del complemento ortogonal

    Sea VV un espacio vectorial de dimensión finita, con producto interno, y sea SS un subespacio de VV. Entonces se verifican las siguientes propiedades:

    1. (S)=S{\left( {{S^ \bot }} \right)^ \bot } = S
    2. V={0V}  {V^ \bot } = \left\{ {{0_V}} \right\}\; y   {0V}=V\;\left\{ {{0_V}} \right\}^ \bot = V
    3. SS={0V}S \cap {S^ \bot } = \left\{ {{0_V}} \right\}
    4. S+S=VS + {S^ \bot } = V

    Esta última propiedad significa que cualquier vector de VV puede expresarse como suma de un vector de S más otro de S{S^ \bot }.

    216

    Ilustramos con un ejemplo geométrico de R3{\mathbb{R}^3}:

    De las propiedades 3 y 4 se deduce:

    SS=VS \oplus {S^ \bot } = V

    Y por lo tanto:

    dim(S)+dim(S)=dim(V)\dim \left( S \right) + \dim \left( {{S^ \bot }} \right) = \dim \left( V \right)

    La unión de una base de SS con una base de S{S^ \bot } es base de V.V. Esto se aplica para extender una base de S  S\;a una base de VV, como muestra el siguiente ejemplo.

    217

    Ejemplo

    Sea S={(x1,x2,x3,x4)R4  x1+x4=0,  x1x2+3x4=0}S = \left\{ {\left( {{x_1},{x_2},{x_3},{x_4}} \right) \in {\mathbb{R}^4}|\;{x_1} + {x_4} = 0,\;{x_1} - {x_2} + 3{x_4} = 0} \right\}

    Hallar una base de SS y extenderla a una base de R4{\mathbb{R}^4}.

    Resolución

    Buscamos una base de S, por ejemplo:

    BS={  (1,2,0,1),  (0,0,1,0)}{B_S} = \left\{ {\;\left( {1, - 2,0, - 1} \right),\;\left( {0,0,1,0} \right)} \right\}

    Como dim(S)=2,  dim\left( S \right) = 2,\; podemos anticipar que: dim(S)=42=2\,dim({S^ \bot }) = 4 - 2 = 2

    A partir de la base de S,  S,\; obtenemos las ecuaciones de S{S^ \bot }:

    S:{x12x2x4=0x3=0 S^ \bot: \begin{cases} {x_1} - 2{x_2} - {x_4} = 0 \\ {x_3} = 0 \end{cases}

    Y hallamos una base de S.  {S^ \bot }.\; por ejemplo:

    BS  ={(2,1,0,0),  ((0,1,0,2)}{B_{{S^{ \bot \;}}}} = \{\left( {2,1,0,0} \right),\;(\left( {0,1,0, -2} \right)\}

    Entonces uniendo las bases de SS y S{S^ \bot } resulta:

    B={(1,2,0,1),(0,0,1,0),(2,1,0,0),(0,1,0,2)}B = \left\{ {\left( {1, - 2,0, - 1} \right),\left( {0,0,1,0} \right),\left( {2,1,0,0} \right),\left( {0,1,0, - 2} \right)} \right\}
    que es base de R4{\mathbb{R}^4}

    218

    Videos con ejercicios resueltos

    Video 3.3. Ejercicio Álgebra y Geometría Analítica UTN FRBA (video de elaboración propia)

    Video 3.4. Ejercicio Álgebra y Geometría Analítica UTN FRBA (video de elaboración propia)

    219
    Capítulo IV

    Tres espacios de una matriz

          Muchos problemas de ingeniería, por complicados que sean, se pueden reducir a álgebra lineal:

    AX=B    oˊ    AX=λX    oˊ    AXB AX = B \;\; ó \;\; AX = \lambda X \;\; ó \;\; AX \approx B

    Si las matrices son cuadradas, el primer problema AX=BAX = B tiene una única solución cuando las columnas de AA son independientes. El segundo problema AX=λX,AX = \lambda X, cuando tiene solución para λ\lambda y para XX produce vectores independientes (llamados vectores propios). Para resolver el tercer problema, AXB,AX \approx B, suele recurrirse a métodos numéricos.
    Se puede observar que: A=(112123134)A =\begin{pmatrix} 1 & 1 & 2 \\ 1 & 2 & 3 \\ 1 & 3 & 4\end{pmatrix} no tiene columnas independientes. Si sumamos las columnas 1 y 2 obtenemos la columna 3. Un resultado sorprendente e importante del álgebra lineal establece que las filas tampoco son independientes.
    Este capitulo analiza La interacción de las columnas y las filas. A continuación se enumeran cuatro de los conceptos más importantes:

    1.   El espacio columna (todas las combinaciones de las columnas).
    2.   El espacio fila (todas las combinaciones de las filas).
    3.   El rango (el número de columnas o filas independientes).
    4.   Eliminación (la mejor manera de encontrar el rango de una matriz).

    223

    4.1   Espacios fila y columna de A. Rango

    En esta presentación, resumimos la teoría necesaria para resolver la ejercitación:

    Método para hallar el rango de una matriz

    Recordemos que:

    1.  Si se realizan operaciones elementales entre las filas de una matriz, el rango se conserva.

    2.  Las filas no nulas de una matriz escalonada son LI.

    Por lo tanto, para determinar el rango de una matriz se aplican operaciones elementales para obtener una matriz escalonada y se cuentan las filas no nulas.

    224

    Compatibilidad y rango

    Consideremos el siguiente sistema:

    {x1+x2           +x4=3x3x4=2x1+x2x3+2x4=0\begin{cases} \begin{aligned}{x_1} + {x_2} \:\:\:\:\:\:\:\:\:\:\: + {x_4} = 3 \\ {x_3} - {x_4} = 2 \\ {x_1} + {x_2} - {x_3}+ 2{x_4} = 0 \end{aligned} \end{cases}

    Es un sistema de 3 ecuaciones con 4 incógnitas (3 x 4) cuya matriz de coeficientes es la AA del ejemplo 4 de los ejemplos precedentes.

    El sistema puede expresarse en forma matricial como: AX=BA\,X = B

    225

    con A=(110100111112)A = \begin{pmatrix} 1 & 1 & 0 & 1\\ 0 & 0 & 1 & -1 \\ 1 & 1 & -1 & 2\end{pmatrix}\: y B=(320)\:B = \begin{pmatrix} 3 \\ 2 \\ 0\end{pmatrix}

    Este sistema, ¿tendrá solución? ¿Cuántas soluciones tendrá?

    Veamos cómo escribir el sistema en función de las columnas de la matriz AA:

    x1(101)+x2(101)+x3(011)+x4(112)=(320) {x_1} \begin{pmatrix} 1 \\ 0 \\ 1\end{pmatrix} + {x_2}\begin{pmatrix} 1 \\ 0 \\ 1\end{pmatrix} + {x_3}\begin{pmatrix} 0 \\ 1 \\ -1\end{pmatrix} + {x_4}\begin{pmatrix} 1 \\ -1 \\ 2\end{pmatrix} = \begin{pmatrix} 3 \\ 2 \\ 0\end{pmatrix}

    Queda una combinación lineal de las columnas de AA igualada al vector de los términos independientes:

    x1  A1+x2  A2+x3  A3+x4  A4=B{x_1}\;{A_1} + {x_2}\;{A_2} + {x_3}\;{A_3} + {x_4}\;{A_4} = B

    ¿Cuándo tiene solución el sistema? Cuando podemos encontrar valores para x1,x2,x3,x4{x_1},{x_2},{x_3},{x_4} que satisfagan la igualdad. Estos valores pueden ser únicos o no.

    En otras palabras:

    El sistema AX=BAX = B es compatible si y sólo si BB es combinación lineal de las columnas de A.A.

    226

    Decir que BB es combinación lineal de las columnas de AA significa que BB está en el subespacio generado por las columnas de AA.
    O sea:

    El  sistema    AX=B    es  compatible          BCol(A)El\;sistema\;\;AX = B\;\;es\;compatible\;\;\; \Leftrightarrow \;\;B \in Col\left( A \right)

    Consideremos la matriz ampliada del sistema, que se obtiene agregando la columna de los términos independientes:

    A=(A1      A2      A3      A4        B)A' = \left( {{A_1}\;\;\;{A_2}\;\;\;{A_3}\;\;\;{A_4}\;\;\;\;B} \right)

    ¿Qué valores puede tomar rg(A)?rg\left( {A'} \right)?

    El rango de AA' dependerá de si BB es combinación lineal o no de las columnas de AA:

    • rg(A)=rg(A)=2BCol(A)rg\left( {A'} \right) = rg\left( A \right) = 2 \Leftrightarrow B \in Col\left( A \right)
    • rg(A)=rg(A)+1=3BCol(A)rg\left( {A'} \right) = rg\left( A \right) + 1 = 3 \Leftrightarrow B \notin Col\left( A \right)

    Estamos en condiciones de enunciar un teorema central sobre la compatibilidad de sistemas de ecuaciones lineales.

    Un sistema de ecuaciones lineales AX=BAX = B es compatible si y sólo si el rango de la matriz ampliada es igual al rango de AA.

    El  sistema    AX=B  es  compatiblerg(A)=rg(A)El\;sistema\;\;AX = B\;es\;compatible \Leftrightarrow rg\left( A \right) = rg\left( {A'} \right)
    227

    Oprimí el botón en la parte inferior derecha del siguiente recurso y podrás trabajar con distintos espacios columna.

    Objeto interactivo de construcción propia
    228

    4.2   Espacio nulo de A.

    Dada una matriz AA, las soluciones de AX=0AX=0 constituyen un espacio vectorial, que llamaremos el espacio nulo de AA.
    Encontremos el Espacio Nulo de la matriz A:

    A=(111231142) A = \begin{pmatrix} -1 & 1 & 1 \\ 2 & 3 & 1 \\ 1 & 4 & 2\end{pmatrix} Debemos encontrar que matrices XX "anulan" a la matriz AA. Es decir, el espacio solución del sistema AX=0.AX=0. La matriz ampliada del sistema es:

    A=(111023101420) A' = \begin{pmatrix} -1 & 1 & 1 & 0 \\ 2 & 3 & 1 & 0\\ 1 & 4 & 2 & 0 \end{pmatrix}

    Resolviendo el sistema, la solución es

    S={X=(xyz)=z(235)    con  zR} S = \begin{Bmatrix} X = \begin{pmatrix} x \\ y \\ z \end{pmatrix} = z \begin{pmatrix} 2 \\ -3 \\ 5 \end{pmatrix} \;\; con \; z \in \mathbb{R} \end{Bmatrix}

    (Espacio nulo de AA)

    229

    En el siguiente applet podemos observar la interpretación geométrica: todos estos vectores XX (en rojo) forman una recta, al ser multiplicados por AA se anulan.

    Objeto interactivo de construcción propia

    Más ejemplos

    Puedes oprimír el botón situado en la parte inferior derecha del recurso de la página siguiente. En la actividad, primero elige el rango de la matriz con la que quieras trabajar y podrás encontrar el espacio nulo asociado a la misma. En el applet, a dicho espacio nulo se lo indica como el "núcleo" de AA (Terminología que tendrá significado cuano estudiemos transformaciones lineales.)

    230
    Objeto interactivo de construcción propia

    Por último, en el siguiente ejemplo podemos cambiar la matriz con la que trabajamos y graficar los tres espacios asociados (columna, fila y nulo) juntos o por separado.

    El applet muestra los espacios de una matriz M de 3x3. El espacio fila se muestra en azul, el espacio columna se muestra en naranja y el espacio nulo se muestra en verde. Puedes editar la matriz para ver cómo cambian estos espacios y puedes alternar la visibilidad de los tres espacios haciendo clic en las casillas de verificación de la derecha. Algunas observaciones notables:

    • Los espacios fila y columna siempre tienen la misma dimensión, igual al rango de M.
    • La dimensión de los espacios fila/columna y la dimensión del espacio nulo siempre suman 3, de acuerdo con el Teorema de nulidad-rango.
    • El espacio de filas siempre es ortogonal al espacio nulo de M.
    231

    Aquí hay algunas matrices interesantes para probar:

    (121242121) \begin{pmatrix} 1 & 2 & -1 \\ 2 & 4 & -2 \\ 1& 2 & -1 \end{pmatrix} Antes de escribirla en el applet, mirá las columnas de esta matriz: ¿qué observas? ¿Puedes adivinar cómo se verá el espacio columna?
    Otra: (124217103) \begin{pmatrix} 1 & 2 & -4 \\ 2 & -1 & 7 \\ 1 & 0 & 3 \end{pmatrix}

    ¿Cuál es el espacio nulo en este caso? (Pista: recuerda que uno de los subespacios 'triviales' consta solo del vector nulo).

    La última: (000000000) \begin{pmatrix} 0 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix} ¿Cuáles serán los espacio fila y columna en este caso?

    Recordá que podés introducir la matriz que desees.

    232
    Objeto interactivo de construcción propia
    233

    4.3 Relación entre las soluciones de AX=BAX = B
    y AX=0.  AX = 0.\; Variables libres.

    4.3.1   Espacio solución de un sistema homogéneo

    Consideremos un sistema homogéneo AX=0,  AX = 0,\; con ARmxnA \in {\mathbb{R}^{mxn}}.

    Sea Sh={XRnx1  AX=  0}{S_h} = \left\{ {X \in {\mathbb{R}^{nx1}}| \; AX = \;0} \right\}, el conjunto solución del sistema.

    1) X=0  X = 0\; pertenece a Sh{S_h} (solución trivial del sistema).

    2) Sean X1{X_1} y X2{X_2} soluciones del sistema.

    A(X1+X2)=  AX1+AX2=0    X1+X2A\left( {{X_1} + {X_2}} \right) = \;A{X_1} + A{X_2} = 0\; \Rightarrow \;{X_1} + {X_2} también es solución.

    3) Sea XX una solución, y kk un número real:

    A(kX)  =k(AX)=k0=0A\left( {kX} \right)\; = k\left( {AX} \right) = k0 = 0 kXkX también es solución.

    Por lo tanto:

    El conjunto solución de un sistema homogéneo AX  =  0AX\; = \;0 (con ARmxnA \in {\mathbb{R}^{mxn}}) es un subespacio de Rnx1{\mathbb{R}^{nx1}}.

    Si el sistema (homogéneo) es compatible determinado, la única solución es la trivial.  En ese caso,   dim(Sh)  =  0.\;dim\left( {{S_h}} \right)\; = \;0.

    ¿De qué dependerá la dimensión de Sh{S_h} para una matriz AA dada?

    234
    Propiedad

    Sea Sh{S_h} el espacio solución del sistema homogéneo AX=0AX = 0, con ARmxn.A \in {\mathbb{R}^{mxn}}.\: Puede demostrarse que:

    dim(Sh)  =  nrg(A)dim\left( {{S_h}} \right)\; = \;n - rg\left( A \right)

    Nota: Recordamos que el conjunto solución de un sistema homogéneo AX=0AX = 0 se denomina espacio nulo de la matriz AA.

    Ejercicio

    Dado el sistema homogéneo cuya matriz de coeficientes es

    (11101011036k) \begin{pmatrix} 1 & 1 & 1 & 0 \\ 1 & 0 & -1 & 1 \\ 0 & 3 & 6 & k \end{pmatrix}

    a) Determinar la dimensión del espacio solución y el rango de AA de acuerdo con los valores de kk.

    b) Para k=0  k = 0\; hallar una base del espacio solución.

    Te pedimos que compruebes en este ejercicio la propiedad mencionada previamente.

    235

    4.3.2   Relación entre las soluciones de AX=BAX = B y AX=0AX = 0

    Consideremos los siguientes sistemas de ecuaciones:

    Sistema de ecuaciones:

    {2xy+z=4x+2yz=3 \begin{cases} \begin{aligned} 2x - y + z = 4 \\ -x +2y - z = 3 \\ \end{aligned} \end{cases}

    Sistema homogéneo asociado

    {2xy+z=0x+2yz=0 \begin{cases} \begin{aligned} 2x - y + z = 0 \\ -x +2y - z = 0 \\ \end{aligned} \end{cases}

    Resolvemos primero el sistema de ecuaciones original:

    {2xy+z=4x+2yz=3        {6+4y2zy+z=4x=2yz3         \begin{cases} \begin{aligned} 2x - y + z = 4 \\ -x +2y - z = 3 \\ \end{aligned} \end{cases} \;\; \Rightarrow \;\; \begin{cases} \begin{aligned} -6 +4y -2z - y + z = 4 \\ x = 2y - z - 3 \\ \end{aligned} \end{cases} \;\; \Rightarrow \;\;         {z=10+3yx=y+7 \;\; \Rightarrow \;\; \begin{cases} \begin{aligned} z = -10 + 3y \\ x = -y + 7 \\ \end{aligned} \end{cases}

    Luego el conjunto solución es:

    (x,y,z)=(y+7,      y,  10+3y)=y(1,1,3)+(7,0,10)\left( {x,y,z} \right) = \left( { - y + 7,\;\;\;y,\; - 10 + 3y} \right) = y\left( { - 1,1,3} \right) + \left( {7,0, - 10} \right)

    S={(x,y,z)R3    (x,y,z)=λ(1,1,3)+(7,0,10)}S = \left\{ {\left( {x,y,z} \right) \in {\mathbb{R}^3}|\;\;\left( {x,y,z} \right) = \lambda \left( { - 1,1,3} \right) + \left( {7,0, - 10} \right)} \right\}

    Geométricamente, el conjunto solución es una recta que no pasa por el origen.

    236

    Resolvamos ahora el sistema homogéneo asociado:

    {2xy+z=0x+2yz=0        {4y2zy+z=0x=2yz         \begin{cases} \begin{aligned} 2x - y + z = 0 \\ -x +2y - z = 0 \\ \end{aligned} \end{cases} \;\; \Rightarrow \;\; \begin{cases} \begin{aligned} 4y -2z - y + z = 0 \\ x = 2y - z \\ \end{aligned} \end{cases} \;\; \Rightarrow \;\;         {z=3yx=y \;\; \Rightarrow \;\; \begin{cases} \begin{aligned} z = 3y \\ x = - y \\ \end{aligned} \end{cases}

    Luego el conjunto solución es:

    (x,y,z)=(y,    y,  3y)=y(1,1,3)\left( {x,y,z} \right) = \left( { - y,\;\;y,\; 3y} \right) = y\left( { - 1,1,3} \right)

    Sh={(x,y,z)  (x,y,z)=λ(1,1,3)}{S_h} = \left\{ {\left( {x,y,z} \right)|\;\left( {x,y,z} \right) = \lambda \left( { - 1,1,3} \right)} \right\}

    Se trata de una recta paralela a la anterior, que pasa por el origen.

    Notemos que tal como establecimos previamente, el conjunto solución de un sistema homogéneo es un subespacio.

    ¿Qué relación existe entre los conjuntos solución del sistema original y de su homogéneo asociado?

    En la gráfica de la siguiente página se muestran las soluciones de los sistemas. En naranja la recta solución del sistema no homogéneo, y en violeta la recta solución del sistema homogéneo:

    237

    Son dos rectas paralelas: la recta que es solución del sistema homogéneo pasa por el origen y la recta que es solución del sistema no homogéneo no pasa por el origen.

    Observamos que:

    S=(x,  y,  z)=λ  (1,  1,  3)Sh+(7,  0,  10)Xp S = {(x,\;y,\;z) = \lambda \; \underbrace{(-1, \;1, \;3)}_{S_{h}} + \underbrace{(7,\; 0, \; -10)}_{X_{p}}}

        S=Sh+Xp \Rightarrow \;\;S = {S_h} + {X_p}

    A continuación veremos un conjunto de propiedades que permiten generalizar este resultado para un sistema de ecuaciones lineales AX=BAX = B con ARm×n,    BRm×1A \in {\mathbb{R}^{m \times n}},\;\;B \in {\mathbb{R}^{m \times 1}}:

    238
    Propiedad 1

    La diferencia de dos soluciones particulares es solución del sistema homogéneo asociado.

    Si X1{X_1} y X2{X_2} son dos soluciones particulares del sistema AX=BAX = B, entonces se cumple que:

    AX1=B        AX2=B A{X_1} = B \;\; \land \;\; A{X_2} = B

    Restando miembro a miembro:

    AX1AX2=0    A  (X1X2)=0 A{X_1} - A{X_2} = 0 \; \Rightarrow \; A \;({{X_1} - {X_2}}) = 0

    Por lo tanto X1X2  {X_1} - {X_2}\; es solución de AX=0AX = 0

    Propiedad 2

    La suma de una solución particular de AX=BAX = B y una solución del homogéneo asociado, es solución de AX=B.AX = B.

    Sean:

    • Xp{X_p} una solución particular de AX=BAX = B
    • Xh{X_h} una solución del sistema homogéneo asociado.

    Queremos probar que (Xp+Xh)\left( {{X_p} + {X_h}} \right) es solución de AX=BAX = B

    A(Xp+Xh)=AXp+AXh=B+0=BA\left( {{X_p} + {X_h}} \right) = A{X_p} + A{X_h} = B + 0 = B

    239
    Propiedad 3

    En la propiedad anterior habíamos probado que (Xp+Xh)\left( {{X_p} + {X_h}} \right) es solución de AX=B.  AX = B.\; Recíprocamente:

    Cualquier solución del sistema AX=BAX = B puede expresarse como Xp+Xh,  {X_p} + {X_h},\; siendo Xp  {X_p}\; una solución particular del sistema y Xh  {X_h}\; una solución del homogéneo asociado.

    Sea XX una solución del sistema AX=B,  AX = B,\; entonces:

    X=(XXp)+XpX = \left( {X - {X_p}} \right) + {X_p}

    Por propiedad 1,   (XXp)\;\left( {X - {X_p}} \right) es solución del homogéneo asociado, por lo tanto:

    X=Xh+XpX = {X_h} + {X_p}

    A partir de las propiedades 2 y 3 se deduce que:

    El conjunto solución de un sistema AX=BAX = B puede expresarse como suma de una solución particular y la solución del sistema homogéneo asociado.

    S=Xp+Sh={Xp+Xh  XhSh} S = {X_p} + {S_h} = \begin{Bmatrix} {{X_p} + {X_h}|\;{X_h} \in {S_h}} \end{Bmatrix}
    240

    4.4 Número de variables libres de un sistema

    Recordemos que:

    El sistema AX=BAX = B es compatible \Leftrightarrow BCol(A)B \in Col\left( A \right) \Leftrightarrow rg(A)=rg(A) \Leftrightarrow rg\left( A \right) = rg\left( {A'} \right)

    Supongamos que BCol(A)B \in Col\left( A \right), es decir que el sistema es compatible. ¿Cuántas soluciones tiene?

    Sea   rg(A)=rg(A)=r\;rg\left( A \right) = rg\left( {A'} \right) = r

    Si nn es el número de incógnitas, puede demostrarse que:

    241
    • Si r=nr = n, el sistema tiene solución única (SCD)
    • Si r<nr < n, el sistema admite infinitas soluciones (SCI), con nrn - r variables libres

    Observación: El número de variables libres es igual a la dimensión del espacio solución del sistema homogéneo asociado.

      Ejemplo

    Analizar cada uno de los siguientes sistemas de ecuaciones, e indicar el número de variables libres.

    a) {x+y+z=0xy   =1z=3 \quad \quad \quad \quad \quad \quad \quad \begin{cases} \begin{aligned} x + y + z = 0 \\ x - y \quad \: \: \: = 1 \\ z = 3 \end{aligned} \end{cases}

    b) {x+y+z=0xy   =12x+z=1 \quad \quad \quad \quad \quad \quad \quad \begin{cases} \begin{aligned} x + y + z = 0 \\ x - y \quad \: \: \: = 1 \\ 2x \quad \: + z = 1 \end{aligned} \end{cases}

    c) {x1+x2+x3+x4=1x2  x4=0x1       +x3+2x4=1 \quad \quad \quad \quad \quad \quad \quad \begin{cases} \begin{aligned} {x_1} + {x_2} + {x_3} + {x_4} = 1 \\ \quad {x_2} \quad \quad\; - {x_4} = 0 \\ {x_1} \:\:\:\:\:\:\: + {x_3} + 2{x_4} = 1 \end{aligned} \end{cases}

    242

      Resolución

    Ítem a

    {x+y+z=0xy   =1z=3 \begin{cases} \begin{aligned} x + y + z = 0 \\ x - y \quad \: \: \: = 1 \\ z = 3 \end{aligned} \end{cases}

    Realicemos el análisis a partir del rango de AA y de AA':

    (111011010033)F2F2F1(111002110033)\begin{pmatrix} 1 & 1 & 1 & 0 \\ 1 & -1 & 0 & 1 \\ 0& 0 & 3 & 3\end{pmatrix} \: F_2 \to F_2 - F_1 \: \begin{pmatrix} 1 & 1 & 1 & 0 \\ 0 & -2 & -1 & 1 \\ 0& 0 & 3 & 3\end{pmatrix}

    Como rg(A)=rg(A)=n=3rg\left( A \right) = rg\left( {A'} \right) = n = 3 entonces el sistema es SCD. No hay variables libres.

    Ítem b

    {x+y+z=0xy   =12x+z=1 \begin{cases} \begin{aligned} x + y + z = 0 \\ x - y \quad \: \: \: = 1 \\ 2x \quad \: + z = 1 \end{aligned} \end{cases}

    Escalonamos AA' para hallar rg(A)rg(A) y rg(A)rg(A'):

    (111011012011)F2F2F1F3F32F1(111002110211)\begin{pmatrix} 1 & 1 & 1 & 0 \\ 1 & -1 & 0 & 1 \\ 2 & 0 & 1 & 1\end{pmatrix} \: \begin{matrix} F_2 \to F_2 - F_1 \\ F_3 \to F_3 - 2F_1\end{matrix} \: \begin{pmatrix} 1 & 1 & 1 & 0 \\ 0 & -2 & -1 & 1 \\ 0 & -2 & -1 & 1 \end{pmatrix}
    243
    (111002110211)F3F3F2(111002110000) \begin{pmatrix} 1 & 1 & 1 & 0 \\ 0 & -2 & -1 & 1 \\ 0 & -2 & -1 & 1\end{pmatrix} F_3 \to F_3 - F_2 \begin{pmatrix} 1 & 1 & 1 & 0 \\ 0 & -2 & -1 & 1 \\ 0 & 0 & 0 & 0 \end{pmatrix}

    Al escalonar AA' se elimina una de las filas. Luego:

    rg(A)=rgA=2   rg(A) = rg{A'} = 2 \; y   n=3      \; n = 3 \; \Rightarrow \;\; SCI con una variable libre.

    Le proponemos al lector hallar el conjunto solución, comprobar que el sistema es un SCI con una variable libre, e indicar una base de Sh{S_h}.

    Ítem c

    {x1+x2+x3+x4=1x2  x4=0x1       +x3+2x4=1 \begin{cases} \begin{aligned} {x_1} + {x_2} + {x_3} + {x_4} = 1 \\ \quad {x_2} \quad \quad\; - {x_4} = 0 \\ {x_1} \:\:\:\:\:\:\: + {x_3} + 2{x_4} = 1 \end{aligned} \end{cases}

    Escalonamos AA' para hallar rg(A)rg(A) y rg(A)rg(A'):

    (111110101010121)F2F2F1(111110101001010)\begin{pmatrix} 1 & 1 & 1 & 1 & 1 \\ 0 & 1 & 0 & -1 & 0 \\ 1 & 0 & 1 & 2 & 1 \end{pmatrix} F_2 \to F_2 - F_1 \begin{pmatrix} 1 & 1 & 1 & 1 & 1 \\ 0 & 1 & 0 & -1 & 0 \\ 0 & -1 & 0 & 1 & 0 \end{pmatrix} (111110101001010)F3F3+F2(111110101000000)\begin{pmatrix} 1 & 1 & 1 & 1 & 1 \\ 0 & 1 & 0 & -1 & 0 \\ 0 & -1 & 0 & 1 & 0 \end{pmatrix} F_3 \to F_3 + F_2 \begin{pmatrix} 1 & 1 & 1 & 1 & 1 \\ 0 & 1 & 0 & -1 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{pmatrix}

    Al escalonar AA' se elimina una de las filas. Luego:

    244

    rg(A)=rgA=2   rg(A) = rg{A'} = 2 \; y   n=4  \; n = 4 \; \Rightarrow SCI con dos variables libres.

    Le proponemos al lector hallar el conjunto solución, comprobar que el sistema es un SCI con dos variables libres, e indicar una base de Sh{S_h}.


    Resumen

    Sea el sistema de ecuaciones AX=B,  AX = B,\; con ARm×n, A \in {\mathbb{R}^{m \times n}},
    XRn×1,  BRm×1,   X \in {\mathbb{R}^{n \times 1}},\;B \in {\mathbb{R}^{m \times 1}},\; cuya matriz ampliada es:

    A=(A1        A2      A3          An        B)A' = \left( {{A_1}\;\;\;\;{A_2}\;\;\;{A_3}\;\; \ldots \;\;\;{A_n}\;\;\;\;B} \right)

    El siguiente esquema resume el análisis de compatibilidad que realizamos previamente:

    245

    En el caso particular de un sistema homogéneo AX=0:AX = 0:

    246

    Acerca de los objetos interactivos de esta obra

    Todos los objetos interactivos del libro hechos en Geogebra han sido diseñados por el autor.

    Los objetos interactivos hechos en Descartesjs utilizados en el capítulo 2 han sido diseñados y creados originalmente por la profesora María José García Cebrian [1] y adaptados a este libro por el autor; con la excepción del interactivo de la página 156, cuya autora es la profesora Melissa Méndez Servin [3] y el objeto interactivo de la página 146 cuyos autores originales son John Jairo García Mora, Sonia Jaquelliny Moreno Jiménez y Margarita Patiño Jaramillo.[4] Este último objeto fue modificado por quien esto escribe para adaptarlo a esta publicación.

    Los objetos interactivos hechos en Descartesjs utilizados en el capítulo 1 han sido diseñados y creados originalmente por las profesoras María José García Cebrian y Esperanza Álvarez Sáiz[2] y adaptados a este libro por el autor.

    247

    Bibliografía recomendada

    • Lay, D. (2012). Algebra Lineal y sus Aplicaciones.
      Ed. Pearson.
    • Kozak, A., Pastorelli, S. y Vardanega, P. (2007). Nociones de Geometría Analítica y Álgebra Lineal. McGraw-Hill, Buenos Aires.
    • Merino González, L. y Santos Aláez E. (2010). Álgebra lineal con metodos elementales. Ediciones PARANINFO.
    • Lipschutz S. (1992). Álgebra lineal. Mc Graw – Hill, Madrid.
    • Noble B. y Daniel J. W. (1989). Álgebra lineal aplicada.
      Ed. Prentice Hall.
    • Nakos G., Joyner D. (1999). Algebra Lineal con Aplicaciones. Thomson Ed.
    • Poole D. (2011) Álgebra lineal, Una introducción moderna. Cengage Learning Editores.
    249