Archivo del Autor: Blanca Radillo

Álgebra Lineal I: Aplicaciones del teorema espectral, bases ortogonales y más propiedades de transformaciones lineales

Por Blanca Radillo

Introducción

Hoy es la última clase del curso. Ha sido un semestre difícil para todas y todos. El quedarnos en casa, obligados a buscar alternativas digitales que sean de fácil acceso para la mayoría de las personas, aprender a realizar toda nuestra rutina diaria en un mismo espacio; sin dudarlo, un semestre lleno de retos que de una u otra manera, haciendo prueba y error, hemos aprendido a sobrellevar.

El día de hoy terminaremos con el tema de teoría espectral. Veremos algunos problemas donde usaremos las técnicas de búsqueda de eigenvalores y eigenvectores, así como aplicaciones de uno de los teoremas más importante: el Teorema Espectral.

Matrices simétricas, matrices diagonalizables

En entradas anteriores hemos discutido sobre qué condiciones me garantizan que una matriz $A$ es diagonalizable. No volveremos a repetir cuál es la definición de matriz diagonalizable ya que en múltiples ocasiones lo hicimos.

Sabemos que una matriz simétrica en $M_{n} (R)$ siempre es diagonalizable, gracias al teorema espectral, pero el siguiente problema nos ilustra que si cambiamos de campo $F$ , no tenemos la garantía de que las matrices simétricas en $M_{n} (F)$ también lo sean.

Problema 1. Demuestra que la matriz simétrica con coeficientes complejos

$A = (\begin{matrix} 1 & i \\ i & - 1 \end{matrix})$

no es diagonalizable.

Solución. Por la primera proposición de la clase «Eigenvalores y eigenvectores de transformaciones y matrices», si $A$ fuese diagonalizable, es decir, que existe una matriz invertible $P$ y una diagonal $D$ tal que $A = P^{- 1} D P$ , entonces $A$ y $D$ tienen los mismos eigenvalores. Entonces, encontremos los eigenvalores de $A$ : buscamos $λ \in C$ tal que $det (λ I - A) = 0$ ,

$\begin{aligned} det (λ I - A) & = | \begin{array}{c} λ - 1 & i \\ i & λ + 1 \end{array} | \\ = (λ - 1) (λ + 1) - i^{2} = λ^{2} - 1 + 1 \\ = λ^{2} = 0. \end{aligned}$

Por lo tanto, el eigenvalor con multiplicidad 2 de $A$ (y también el eigenvalor de $D$ ) es $λ = 0$ . Si $D$ es de la forma

$D = (\begin{matrix} a & 0 \\ 0 & b \end{matrix})$ ,

es fácil ver (y calcular) que sus eigenvalores son $a$ y $b$ , pero por lo anterior, podemos concluir que $a = b = 0$ , y por lo tanto $D$ es la matriz cero. Si fuese así, $A = P^{- 1} D P = 0$ , contradiciendo la definición de $A$ .

Problema 2. Sea $A$ una matriz simétrica con entradas reales y supongamos que $A^{k} = I$ para algún entero positivo $k$ . Prueba que $A^{2} = I$ .

Solución. Dado que $A$ es simétrica y con entradas reales, todos sus eigenvalores son reales. Más aún son $k$ -raíces de la unidad, entonces deben ser $\pm 1$ . Esto implica que todos los eigenvalores de $A^{2}$ son iguales a 1. Dado que $A^{2}$ también es simétrica, es diagonalizable y, dado que sus eigenvalores son iguales a 1, por lo tanto $A^{2} = I$ .

Más propiedades de transformaciones lineales y bases ortogonales

En otras clases como Cálculo, Análisis, hablamos de funciones continuas, discontinuas, acotadas, divergentes; mientras que en este curso nos hemos enfocado únicamente en la propiedad de linealidad de las transformaciones. Si bien no es interés de este curso, podemos adelantar que, bajo ciertas condiciones del espacio $V$ , podemos tener una equivalencia entre continuidad y acotamiento de una transformación.

Decimos que la norma de una transformación está definida como

$‖ T ‖ = sup_{x \in V ∖ 0} \frac{‖ T (x) ‖}{‖ x ‖}$ .

Por ende, decimos que una transformación es acotada si su norma es acotada, $‖ T ‖ < \infty$ .

Problema 1. Sea $V$ un espacio euclideano y sea $T$ una transformación lineal simétrica en $V$ . Sean $λ_{1}, \dots, λ_{n}$ los eigenvalores de $T$ . Prueba que

$sup_{x \in V ∖ 0} \frac{‖ T (x) ‖}{‖ x ‖} = max_{1 \leq i \leq n} | λ_{i} | .$

Solución. Renumerando a los eigenvalores, podemos decir que $max_{i} | λ_{i} | = | λ_{n} |$ . Sea $e_{1}, \dots, e_{n}$ una base ortonormal de $V$ tal que $T (e_{i}) = λ_{i} e_{i}$ para todo $i$ . Si $x \in V ∖ 0$ , podemos escribirlo como $x = x_{1} e_{1} + \dots + x_{n} e_{n}$ para algunos reales $x_{i}$ . Entonces, por linealidad de $T$ ,

$T (x) = \sum_{i = 1}^{n} λ_{i} x_{i} e_{i} .$

Dado que $| λ_{i} | \leq | λ_{n} |$ para toda $i$ , tenemos que

$\frac{‖ T (x) ‖}{‖ x ‖} = \sqrt{\frac{\sum_{i = 1}^{n} λ_{i}^{2} x_{i}^{2}}{\sum_{i = 1}^{n} x_{i}^{2}}} \leq | λ_{n} |,$

por lo tanto

$\begin{aligned} max_{1 \leq i \leq n} | λ_{i} | & = | λ_{n} | = \frac{‖ T (e_{n}) ‖}{‖ e_{n} ‖} \\ \leq sup_{x \in V ∖ 0} \frac{‖ T (x) ‖}{‖ x ‖} \\ \leq | λ_{n} | = max_{1 \leq i \leq n} | λ_{i} | . \end{aligned}$

Obteniendo lo que queremos.

Para finalizar, no olvidemos que una matriz es diagonalizable si y sólo si el espacio tiene una base de eigenvectores, y que está íntimamente relacionado con el teorema espectral.

Problema 2. Encuentra una base ortogonal consistente con los eigenvectores de la matriz

$A = \frac{1}{7} (\begin{matrix} - 2 & 6 & - 3 \\ 6 & 3 & 2 \\ - 3 & 2 & 6 \end{matrix}) .$

Solución. Para encontrar los eigenvectores, primero encontrar los eigenvalores y, después, para cada eigenvalor, encontrar el/los eigenvectores correspondientes.

Calculemos:

$\begin{aligned} 0 & = det (λ I - A) = | \begin{array}{c} λ + 2 / 7 & - 6 / 7 & 3 / 7 \\ - 6 / 7 & λ - 3 / 7 & - 2 / 7 \\ 3 / 7 & - 2 / 7 & λ - 6 / 7 \end{array} | \\ = λ^{3} - λ^{2} - λ + 1 \\ = (λ - 1) (λ^{2} - 1), \end{aligned}$

entonces los eigenvalores de $A$ son $1, - 1$ , ( $λ = 1$ tiene multiplicidad 2).

Ahora, hay que encontrar los vectores $v = (x, y, z)$ tal que $A v = λ v$ , para todo eigenvalor $λ$ .

Si $λ = - 1$ ,

$(λ I - A) v = \frac{1}{7} (\begin{matrix} - 5 & - 6 & 3 \\ - 6 & - 10 & - 2 \\ 3 & - 2 & - 13 \end{matrix}) v = 0,$

reduciendo, obtenemos que $v = (3 α, - 2 α, α)$ para todo $α \in R$ .

Si $λ = 1$ , resolviendo de la misma manera $(λ I - A) v = (I - A) v = 0$ , tenemos que $v = (β, γ, - 3 β + 2 γ)$ para todo $β, γ$ . Entonces el conjunto de eigenvectores es

$B = {v_{1} = (3, - 2, 1), v_{2} = (1, 0, - 3), v_{3} = (0, 1, 2)} .$

Es fácil ver que el conjunto $B$ es linealmente independiente, más aún $dim (R^{3}) = 3 = | B |$ , por lo tanto, $B$ es la base consistente con los eigenvectores de $A$ .

$△$

Agradecemos su esfuerzo por llegar hasta el final a pesar de todas las adversidades. Esperamos pronto volver a ser sus profesores/ayudantes. Mucha suerte en la última parcial, es el último esfuerzo. Pero también les deseamos mucho éxito en su proyecto de vida. ¡Gracias!

Entradas relacionadas

Ir a Álgebra Lineal I
Entrada anterior del curso: Teorema espectral para matrices simétricas reales

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104721 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM»

Álgebra Lineal I: Problemas de determinantes y ecuaciones lineales

Por Blanca Radillo

2 respuestas

Introducción

En esta entrada, realizaremos problemas que nos ayudarán a repasar el tema visto el pasado lunes, sobre soluciones de sistemas lineales, Teorema de Rouché-Capelli y la regla de Cramer.

Problemas de ecuaciones lineales

Una de las maneras más usuales para demostrar que un conjunto de vectores es linealmente independientes es probar que tomamos una combinación lineal de éstos tal que es igual a 0, sólo es posible si todos los coeficientes son igual a cero. Pero como ya lo hemos visto anteriormente en diversos problemas, algunas veces ésto nos genera un sistema de ecuaciones que puede ser difícil y/o tardado resolver.

Por ello, otra manera de demostrar independencia lineal es ilustrada con el siguiente problema.

Problema 1. Considera los vectores

$v_{1} = (1, x, 0, 1), v_{2} = (0, 1, 2, 1), v_{3} = (1, 1, 1, 1)$

en $R^{4}$ . Prueba que para cualquier elección de $x \in R$ , los vectores $v_{1}, v_{2}, v_{3}$ son linealmente independientes.

Solución. Sea $A$ la matriz cuyas columnas son $v_{1}, v_{2}, v_{3}$ , es decir,

$A = (\begin{matrix} 1 & 0 & 1 \\ x & 1 & 1 \\ 0 & 2 & 1 \\ 1 & 1 & 1 \end{matrix}) .$

Sabemos que $v_{1}, v_{2}, v_{3}$ son linealmente independiente si y sólo si $dim(span (v_{1}, v_{2}, v_{3})) = 3$ , ya que $rank (A) = 3$ , y eso es equivalente (por la clase del lunes) a demostrar que $A$ tiene una submatriz de $3 \times 3$ invertible.

Notemos que si borramos el segundo renglón, obtenemos la submatriz cuyo determinante es

$| \begin{matrix} 1 & 0 & 1 \\ 0 & 2 & 1 \\ 1 & 1 & 1 \end{matrix} | = - 1,$

lo que implica que es invertible, y por lo tanto $v_{1}, v_{2}, v_{3}$ son vectores linealmente independientes.

En este curso, los ejemplos usualmente utilizan espacios vectoriales sobre $R$ o sobre $C$ . Como $R \subset C$ , es natural preguntarnos si los resultados obtenidos en los problemas trabajados en $R$ se cumplen en $C$ . En este caso particular, si las soluciones de una matriz en $M_{m, n} (R)$ son soluciones de la misma matriz pero vista como elemento en $M_{m, n} (C)$ . El siguiente teorema nos da el resultado a esta pregunta.

Teorema. Sea $A \in M_{m, n} (F)$ y sea $F_{1}$ un campo contenido en $F$ . Consideremos el sistema lineal $A X = 0$ . Si el sistema tiene una solución no trivial en $F_{1}^{n}$ , entonces tiene una solución no trivial en $F^{n}$ .

Demostración. Dado que el sistema tiene una solución no trivial en $F_{1}^{n}$ , $r := rank (A) < n$ vista como elemento en $M_{m, n} (F_{1})$ . Por el primer teorema visto en la clase del lunes, el rango es el tamaño de la submatriz cuadrada más grande que sea invertible, y eso es independiente si se ve a $A$ como elemento de $M_{m, n} (F_{1})$ o de $M_{m, n} (F)$ . Y por el teorema de Rouché-Capelli, el conjunto de soluciones al sistema es un subespacio de $F^{n}$ de dimensión $n - r > 0$ . Por lo tanto, el sistema $A X = 0$ tiene una solución no trivial en $F^{n}$ .

A continuación, se mostrarán dos ejemplos de la búsqueda de soluciones a sistemas lineales donde usaremos todas las técnicas aprendidas a lo largo de esta semana.

Problema. 2 Sea $S_{a}$ el siguiente sistema lineal:

$\begin{matrix} x - 2 y + z = 1 \\ 3 x + 2 y - 2 z = 2 \\ 2 x - y + a z = 3 \end{matrix} .$

Encuentra los valores de $a$ para los cuales el sistema no tiene solución, tiene exactamente una solución y tiene un número infinito de soluciones.

Solución. El sistema lo podemos escribir como $A X = b$ donde

$A = (\begin{matrix} 1 & - 2 & 1 \\ 3 & 2 & - 2 \\ 2 & - 1 & a \end{matrix}) y b = (\begin{matrix} 1 \\ 2 \\ 3 \end{matrix}) .$

Notemos que

$| \begin{matrix} 1 & - 2 & 1 \\ 3 & 2 & - 2 \\ 2 & - 1 & a \end{matrix} | = 8 a - 1,$

entonces si $a \neq 1 / 8$ , $A$ es invertible, y por lo tanto $rank (A) = 3$ , mientras que si $a = 1 / 8$ , $A$ no es invertible y $rank (A) = 2$ ya que la submatriz es invertible

$| \begin{matrix} 1 & - 2 \\ 3 & 2 \end{matrix} | = 8.$

Además, si la matriz $(A, b)$ es igual a

$(\begin{matrix} 1 & - 2 & 1 & 1 \\ 3 & 2 & - 2 & 2 \\ 2 & - 1 & a & 3 \end{matrix}),$

quitando la tercera columna, obtenemos una submatriz invertible (ejercicio). Por lo tanto, $rank (A, b) = 3$ .

Aplicando el Teorema de Rouché-Capelli, para $a = 1 / 8$ , el sistema $A X = b$ no tiene soluciones. También podemos concluir que como $rank (A) = 3$ para todo $a \neq 1 / 8$ , el sistema tiene exactamente una solución. (Y $A X = b$ nunca tiene infinitas soluciones).

$△$

Problema 3. Sean $a, b, c$ números reales dados. Resuelve el sistema lineal

$\begin{matrix} (b + c) x + b y + c z = 1 \\ a x + (a + c) y + c z = 1 \\ a x + b y + (a + b) z = 1 \end{matrix} .$

Solución. La matriz del sistema es

$A = (\begin{matrix} b + c & b & c \\ a & a + c & c \\ a & b & a + b \end{matrix}) .$

No es difícil ver que $det (A) = 4 a b c$ . Si $a b c \neq 0$ , usando la regla de Cramer, la única solución al sistema está dada por

$x = \frac{| \begin{matrix} 1 & b & c \\ 1 & a + c & c \\ 1 & b & a + b \end{matrix} |}{4 a b c}, y = \frac{| \begin{matrix} b + c & 1 & c \\ a & 1 & c \\ a & 1 & a + b \end{matrix} |}{4 a b c}$

$y = \frac{| \begin{matrix} b + c & b & 1 \\ a & a + c & 1 \\ a & b & 1 \end{matrix} |}{4 a b c},$

resolviendo los determinantes obtenemos que

$x = \frac{a^{2} - (b - c)^{2}}{4 a b c}, y = \frac{b^{2} - (a - c)^{2}}{4 a b c}, z = \frac{c^{2} - (a - b)^{2}}{4 a b c} .$

Ahora, si $a b c = 0$ , entonces $A$ no es invertible ( $rank (A) < 3$ ). El sistema es consistente si y sólo si $rank (A) = rank (A, b)$ .

Sin pérdida de generalidad, decimos que $a = 0$ (pues $a b c = 0$ ). Esto reduce el sistema a

$\begin{matrix} (b + c) x + b y + c z = 1 \\ c (y + z) = 1 \\ b (y + z) = 1 \end{matrix} .$

El sistema es consistente si $b = c$ y distintos de cero. En este caso, tenemos que $b (2 x + y + z) = 1$ y $b (y + z) = 1$ , implicando $x = 0$ , $y + z = 1 / b$ . De manera similar, obtenemos las posibles soluciones si $b = 0$ o si $c = 0$ .

Resumiendo:

Si $a b c \neq 0$ , el sistema tiene una solución única dada por la regla de Cramer.
Si tenemos alguno de los siguientes tres casos: caso 1) $a = 0$ y $b = c \neq 0$ ; caso 2) $b = 0$ y $a = c \neq 0$ ; caso 3) $c = 0$ y $a = b \neq 0$ , tenemos infinitas soluciones descritas como, para todo $w \in R$ : caso 1) $(0, w, 1 / b - w)$ ; caso 2) $(w, 0, 1 / a - w)$ ; caso 3) $(w, 1 / a - w, 0)$ .
Si no se cumplen ninguno de las cuatro condiciones anteriores para $a, b, c$ , el sistema no es consistente.

$△$

Entradas relacionadas

Ir a Álgebra Lineal I
Entrada anterior del curso: Determinantes en sistemas de ecuaciones lineales y regla de Cramer
Siguiente entrada del curso: Eigenvectores y eigenvalores de transformaciones y matrices

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104721 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM»

Álgebra Lineal I: Problemas de definición y propiedades de determinantes

Por Blanca Radillo

3 respuestas

Introducción

En esta entrada haremos una serie de problemas que nos ayudarán como repaso de los temas vistos durante las últimas dos semanas. Mostraremos algunas propiedades bastante interesantes acerca de las transformaciones alternantes y antisimétricas, así como de la transformación estrella de esta semana: el determinante.

Problemas de transformaciones antisimétricas

En la entrada del miércoles 6 de mayo, hablábamos sobre la equivalencia entre transformaciones alternantes y antisimétricas, justo resaltamos que ésto no es cierto si el campo $F$ es $Z_{2}$ , y el siguiente ejemplo lo expone:

Ejemplo. Sea $f : Z_{2} \times Z_{2} \to Z_{2}$ definido como $f (x, y) = x y$ . Claramente $f$ es bilineal, pero no es alternate ya que $f (1, 1) = 1 \neq 0$ . Por otro lado, $f$ es antisimétrica, porque $f (x, y) + f (y, x) = x y + y x = 2 x y = 0$ .

$△$

De manera natural surge la pregunta: ¿cómo podemos construir una transformación $d$ -lineal antisimétrica o alternante? El siguiente problema muestra un camino para obtener una transformación antisimétrica dada un mapeo $d$ -lineal $f$ .

Problema. Sea $f : V^{d} \to W$ una transformación $d$ -lineal. Demuestra que

$A (f) := \sum_{σ \in S_{d}} sign (σ) σ (f)$

es un mapeo $d$ -lineal antisimétrico.

Solución. Es fácil ver que $A (f)$ es una transformación $d$ -lineal, dado que $A (f)$ es una combinación lineal de mapeos $d$ -lineales. Queremos probar que, para $τ \in S_{d}$ , $τ (A (f)) = sign (τ) A (f)$ . Notemos que

$\begin{aligned} τ (A (f)) & = \sum_{σ \in S_{d}} sign (σ) τ (σ (f)) \\ = \sum_{σ \in S_{d}} sign (σ) (τ σ) (f) . \end{aligned}$

Usando el hecho que $sign (τ) sign (σ) = sign (τ σ)$ y que ${τ σ : σ \in S_{d}} = S_{d}$ , obtenemos que

$\begin{aligned} sign (τ) τ (A (f)) & = \sum_{σ \in S_{d}} sign (τ σ) (τ σ) (f) \\ = \sum_{η \in S_{d}} sign (η) (η) (f) = A (f) . \end{aligned}$

Por lo tanto, $τ (A (f)) = sign (τ) A (f)$ .

Problemas de determinantes

Ahora continuando con la discusiones del determinante, sabemos que éste es una forma $n$ -lineal alternante, y además que cualquier otra forma $n$ -lineal alternante varía de $det (b_{1}, \dots, b_{n})$ únicamente por un factor multiplicativo. Otro resultado interesante ese teorema es el siguiente:

Problema 1. Sea $V$ un espacio vectorial sobre $F$ de dimensión finita. Sea $e_{1}, \dots, e_{n}$ una base de $V$ y sea $T : V \to V$ una transformación lineal. Demuestra que para todo $v_{1}, \dots, v_{n} \in V$ tenemos que

$\sum_{i = 1}^{n} det (v_{1}, \dots, v_{i - 1}, T (v_{i}), v_{i + 1}, \dots, v_{n}) = Tr (T) \cdot det (v_{1}, \dots, v_{n}),$

donde todos los determinantes están calculados en la base canónica y $Tr (T)$ es la traza de la matriz de $T$ (con respecto a la base canónica).

Solución. Definimos el mapeo $ϕ : V^{n} \to F$ como

$ϕ (v_{1}, \dots, v_{n}) = \sum_{i = 1}^{n} det (v_{1}, \dots, v_{i - 1}, T (v_{i}), v_{i + 1}, \dots, v_{n}) .$

Esta transformación es la suma de transformaciones $n$ -lineales, por lo tanto $ϕ$ es $n$ -lineal. Más aún, es alternante, ya que si asumimos, por ejemplo, que $v_{1} = v_{2}$ , entonces

$\begin{aligned} ϕ (v_{1}, v_{1}, v_{3}, \dots, v_{n}) & = det (T (v_{1}), v_{1}, v_{3}, \dots, v_{n}) + det (v_{1}, T (v_{1}), v_{3}, \dots, v_{n}) \\ + \sum_{i = 3}^{n} det (v_{1}, v_{1}, \dots, v_{i - 1}, T (v_{i}), v_{i + 1}, \dots, v_{n}) \\ = det (T (v_{1}), v_{1}, v_{3}, \dots, v_{n}) + det (v_{1}, T (v_{1}), v_{3}, \dots, v_{n}) \\ = det (T (v_{1}), v_{1}, v_{3}, \dots, v_{n}) - det (T (v_{1}), v_{1}, v_{3}, \dots, v_{n}) \\ = 0, \end{aligned}$

debido a que el determinante es antisimétrico.

Por el último teorema visto en la clase del viernes pasado, existe escalar $α$ tal que

$ϕ (v_{1}, \dots, v_{n}) = α det (v_{1}, \dots, v_{n})$

para todo $v_{1}, \dots, v_{n}$ . Sea $A = [a_{i j}]$ la matriz de $T$ con respecto a la base canónica. Si tomamos $v_{1} = e_{1}, \dots, v_{n} = e_{n}$ , por el mismo teorema tenemos que

$\begin{aligned} α & = ϕ (e_{1}, \dots, e_{n}) \\ = \sum_{i = 1}^{n} det (e_{1}, \dots, e_{i - 1}, \sum_{j = 1}^{n} a_{j i} e_{j}, e_{i + 1}, \dots, e_{n}) \\ = \sum_{i = 1}^{n} \sum_{j = 1}^{n} a_{j i} det (e_{1}, \dots, e_{i - 1}, e_{j}, e_{i + 1}, \dots, e_{n}) \\ = \sum_{i = 1}^{n} a_{i i} = Tr (T) . \end{aligned}$

Por lo tanto, obtenemos lo que queremos.

Por último, los siguientes dos problemas nos ilustran como podemos obtener información de las matrices de manera fácil y «bonita», usando algunas propiedades de los determinantes vistas en la sesión del martes pasado.

Problema 2. Sea $n$ un número impar y sean $A, B \in M_{n} (R)$ matrices tal que $A^{2} + B^{2} = 0_{n}$ . Prueba que la matriz $A B - B A$ no es invertible.

Solución. Notemos que

$(A + i B) (A - i B) = A^{2} + B^{2} + i (B A - A B) = i (B A - A B) .$

Por la propiedad del determinante de un producto, tenemos que

$det (A + i B) det (A - i B) = i^{n} det (B A - A B) .$

Suponemos que $A B - B A$ es invertible, entonces $det (B A - A B) \neq 0$ . Además sabemos que

$det (A - i B) = det (\overset{―}{A + i B}) = \overset{―}{det (A + i B)},$

esto implica que $| det (A + i B) |^{2} = i^{n} det (B A - A B) .$ Como consecuencia, $i^{n}$ es un número real, contradiciendo al hecho que $n$ es impar. Por lo tanto $det (B A - A B) = 0$ .

Problema 3. Para $1 \leq i, j \leq n$ , definimos $a_{i j}$ como el número de divisores positivos en común de $i$ y $j$ y definimos $b_{i j}$ igual a 1 si $j$ divide $i$ e igual a 0 si no.

Probar que $A = B \cdot^{t} B$ , donde $A = [a_{i j}]$ y $B = [b_{i j}]$ .
¿Qué podemos decir de la forma de $B$ ?
Calcula $det (A)$ .

Solución. 1) Fijando $i, j$ tenemos que

$det (B \cdot^{t} B)_{i j} = \sum {k = 1}^{n} b_{i k} b_{j k} .$

Notemos que $b_{i k} b_{j k}$ no es cero ( $b_{i j}, b_{j k} = 1$ ) si y sólo si $k$ divide a $i$ y a $j$ , esto implica que la cantidad de términos de la suma no ceros corresponde exactamente con la cantidad de los divisores en común que tengan $i$ y $j$ . Por lo tanto $det (B \cdot^{t} B)_{i j} = a_{i j}$ .

2) Si $i < j$ , no es posible que $j$ divida a $i$ . Entonces $b_{i j} = 0$ para todo $i < j$ , esto significa que $B$ es, al menos, triangular inferior. Un dato más que podemos asegurar es que $b_{i i} = 1$ para toda $i$ , por lo tanto, al menos, todos los términos de la diagonal de $B$ son iguales a 1.

3) Dada la propiedad multiplicativa del determinante, dado que $det (B) = det (^{t} B)$ y usando el inciso (1), tenemos que $det (A) = det (B \cdot^{t} B) = (det B)^{2} .$ Pero por el inciso (2), $det B = 1$ , concluimos que $det A = 1$ .

$△$

Entradas relacionadas

Ir a Álgebra Lineal I
Entrada anterior del curso: Propiedades de determinantes
Siguiente entrada del curso: Técnicas básicas de cálculo de determinantes

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104721 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM»

Álgebra Lineal I: Problemas de bases ortogonales, Fourier y proceso de Gram-Schmidt

Por Blanca Radillo

2 respuestas

Introducción

Durante las últimas clases hemos visto problemas y teoremas que nos demuestran que las bases ortogonales son extremadamente útiles en la práctica, ya que podemos calcular fácilmente varias propiedades una vez que tengamos a nuestra disposición una base ortogonal del espacio que nos interesa. Veamos más problemas de bases ortogonales y otros resultados que nos permitirán reforzar estas ideas.

Problemas resueltos de bases ortogonales y proyecciones

Para continuar con este tema, veremos que las bases ortogonales nos permiten encontrar de manera sencilla la proyección de un vector sobre un subespacio. Primero, recordemos que si $V = W \oplus W_{2}$ , para todo $v \in V$ podemos definir su proyección en $W$ , que denotamos $π_{W} (v)$ , como el único elemento en $W$ tal que $v - π_{W} (v) \in W_{2}$ .

Debido a las discusiones sobre bases ortogonales, no es difícil ver que si $⟨ w, u ⟩ = 0$ para todo $w \in W$ , entonces $u \in W_{2}$ . Como consecuencia de esto, tenemos el siguiente resultado:

Teorema. Sea $V$ un espacio vectorial sobre $R$ con producto interior $⟨ \cdot, \cdot ⟩$ , y sea $W$ un subespacio de $V$ de dimensión finita. Sea $v_{1}, \dots, v_{n}$ una base ortogonal de $W$ . Entonces para todo $v \in V$ tenemos que

$π_{W} (v) = \sum_{i = 1}^{n} \frac{⟨ v, v_{i} ⟩}{{‖ v_{i} ‖}^{2}} v_{i} .$

Demostración. Escribimos $v$ como $v = π_{W} (v) + u$ con $u \in W_{2}$ . Por la observación previa al teorema, $⟨ u, v_{i} ⟩ = 0$ para todo $i$ . Además existen $a_{1}, \dots, a_{n}$ tales que $π_{W} (v) = a_{1} v_{1} + \dots + a_{n} v_{n}$ . Entonces

$\begin{aligned} 0 & = ⟨ u, v_{i} ⟩ = ⟨ v, v_{i} ⟩ - ⟨ π_{W} (v), v_{i} ⟩ \\ = ⟨ v, v_{i} ⟩ - \sum_{j = 1}^{n} a_{j} ⟨ v_{j}, v_{i} ⟩ \\ = ⟨ v, v_{i} ⟩ - a_{i} ⟨ v_{i}, v_{i} ⟩, \end{aligned}$

porque $v_{1}, \dots, v_{n}$ es una base ortogonal. Por lo tanto, para todo $i$ , obtenemos

$a_{i} = \frac{⟨ v, v_{i} ⟩}{{‖ v_{i} ‖}^{2}} .$

Distancia de un vector a un subespacio y desigualdad de Bessel

En la clase de ayer, vimos la definición de distancia entre dos vectores. También se puede definir la distancia entre un vector y un subconjunto como la distancia entre el vector y el vector «más cercano» del subconjunto, en símbolos:

$d (v, W) = min_{x \in W} ‖ x - v ‖ .$

Dado que $x \in W$ , $x - π_{W} (v) \in W$ , y por definición de proyección $v - π_{W} (v) \in W_{2}$ , entonces

$\begin{aligned} {‖ x - v ‖}^{2} & = {‖ (x - π_{W} (v)) + (π_{W} (v) - v) ‖}^{2} \\ = {‖ x - π_{W} (v) ‖}^{2} + 2 ⟨ x - π_{W} (v), π_{W} (v) - v ⟩ + {‖ π_{W} (v) - v ‖}^{2} \\ = {‖ x - π_{W} (v) ‖}^{2} + {‖ π_{W} (v) - v ‖}^{2} \\ \geq {‖ π_{W} (v) - v ‖}^{2} . \end{aligned}$

Y dado que la proyección pertenece a $W$ , la desigualdad anterior muestra que la proyección es precisamente el vector en $W$ con el que $v$ alcanza la distancia a $W$ . En conclusión, $d (v, W) = ‖ π_{W} (v) - v ‖ .$

Teorema. Sea $V$ un espacio vectorial sobre $R$ con producto interior $⟨ \cdot, \cdot ⟩$ , y sea $W$ un subespacio de $V$ de dimensión finita. Sea $v_{1}, \dots, v_{n}$ una base ortonormal de $W$ . Entonces para todo $v \in V$ tenemos que

$π_{W} (v) = \sum_{i = 1}^{n} ⟨ v, v_{i} ⟩ v_{i},$

$\begin{aligned} d (v, W)^{2} & = {‖ v - \sum_{i = 1}^{n} ⟨ v, v_{i} ⟩ v_{i} ‖}^{2} \\ = {‖ v ‖}^{2} - \sum_{i = 1}^{n} ⟨ v, v_{i} ⟩^{2} . \end{aligned}$

En particular

$\sum_{i = 1}^{n} ⟨ v, v_{i} ⟩^{2} \leq {‖ v ‖}^{2} .$

A esta última desigualdad se le conoce como desigualdad de Bessel.

Demostración. Por el teorema anterior y dado que $v_{1}, \dots, v_{n}$ es una base ortonormal, obtenemos la primera ecuación. Ahora, por Pitágoras,

$d (v, W)^{2} = {‖ v - π_{W} (v) ‖}^{2} = {‖ v ‖}^{2} - {‖ π_{W} (v) ‖}^{2} .$

Por otro lado, tenemos que

$\begin{aligned} {‖ π_{W} (v) ‖}^{2} & = {‖ \sum_{i = 1}^{n} ⟨ v, v_{i} ⟩ v_{i} ‖}^{2} \\ = \sum_{i, j = 1}^{n} ⟨ ⟨ v, v_{i} ⟩ v_{i}, ⟨ v, v_{j} ⟩ v_{j} ⟩ \\ = \sum_{i, j = 1}^{n} ⟨ v, v_{i} ⟩ ⟨ v, v_{j} ⟩ ⟨ v_{i}, v_{j} ⟩ \\ = \sum_{i = 1}^{n} ⟨ v, v_{i} ⟩^{2} . \end{aligned}$

Por lo tanto, se cumple la igualdad de la distancia. Finalmente como $d (v, W)^{2} \geq 0$ , inmediatamente tenemos la desigualdad de Bessel.

Veamos ahora dos problemas más en los que usamos la teoría de bases ortonormales.

Aplicación del proceso de Gram-Schmidt

Primero, veremos un ejemplo más del uso del proceso de Gram-Schmidt.

Problema. Consideremos $V$ como el espacio vectorial de polinomios en $[0, 1]$ de grado a lo más $2$ , con producto interior definido por $⟨ p, q ⟩ = \int_{0}^{1} x p (x) q (x) d x .$

Aplica el algoritmo de Gram-Schmidt a los vectores $1, x, x^{2}$ .

Solución. Es fácil ver que ese sí es un producto interior en $V$ (tarea moral). Nombremos $v_{1} = 1, v_{2} = x, v_{3} = x^{2}$ . Entonces

$e_{1} = \frac{v_{1}}{‖ v_{1} ‖} = \sqrt{2} v_{1} = \sqrt{2},$

ya que ${‖ v_{1} ‖}^{2} = \int_{0}^{1} x d x = \frac{1}{2} .$

Sea $z_{2} = v_{2} - ⟨ v_{2}, e_{1} ⟩ e_{1}$ . Calculando, $⟨ v_{2}, e_{1} ⟩ = \int_{0}^{1} \sqrt{2} x^{2} d x = \frac{\sqrt{2}}{3} .$ Entonces $z_{2} = x - \frac{\sqrt{2}}{3} \sqrt{2} = x - \frac{2}{3} .$ Esto implica que

$e_{2} = \frac{z_{2}}{‖ z_{2} ‖} = 6 (x - \frac{2}{3}) = 6 x - 4.$

Finalmente, sea $z_{3} = v_{3} - ⟨ v_{3}, e_{1} ⟩ e_{1} - ⟨ v_{3}, e_{2} ⟩ e_{2}$ . Haciendo los cálculos obtenemos que

$z_{3} = x^{2} - (\frac{\sqrt{2}}{4}) \sqrt{2} - (\frac{1}{5}) (6 x - 4)$

$z_{3} = x^{2} - \frac{6}{5} x + \frac{3}{10} .$

Por lo tanto

$e_{3} = \frac{z_{3}}{‖ z_{3} ‖} = 10 \sqrt{6} (x^{2} - \frac{6}{5} x + \frac{3}{10}) .$

$△$

El teorema de Plancherel y una fórmula con $π$

Finalmente, en este ejemplo, usaremos técnicas de la descomposición de Fourier para solucionar un problema bonito de series.

Problema. Consideremos la función $2 π -$ periódica $f : R \to R$ definida como $f (0) = f (π) = 0,$ $f (x) = - 1 - \frac{x}{π}$ en el intervalo $(- π, 0)$ , y $f (x) = 1 - \frac{x}{π}$ en el intervalo $(0, π)$ .

Problemas de bases ortogonales: Aplicando el teorema de Plancherel para una fórmula que involucra a pi. — Gráfica de la función $f$ .

Usa el teorema de Plancherel para deducir las identidades de Euler

$\begin{aligned} \sum_{n = 1}^{\infty} \frac{1}{n^{2}} & = \frac{π^{2}}{6}, \\ \sum_{n = 0}^{\infty} \frac{1}{(2 n + 1)^{2}} & = \frac{π^{2}}{8} . \end{aligned}$

Solución. Notemos que no sólo es $2 π -$ periódica, también es una función impar, es decir, $f (- x) = - f (x)$ . Por lo visto en la clase del miércoles pasado tenemos que calcular

$a_{0} (f) = \frac{1}{π} \int_{- π}^{π} f (x) d x,$

$a_{k} (f) = \frac{1}{π} \int_{- π}^{π} f (x) c o s (k x) d x,$

$b_{k} (f) = \frac{1}{π} \int_{- π}^{π} f (x) s e n (k x) d x .$

Para no hacer más larga esta entrada, la obtención de los coeficientes de Fourier se los dejaremos como un buen ejercicio de cálculo. Para hacer las integrales hay que separar la integral en cada uno de los intervalos $[- π, 0]$ y $[0, π]$ y en cada uno de ellos usar integración por partes.

El resultado es que para todo $k \geq 1$ , $a_{0} = 0, a_{k} = 0, b_{k} = \frac{2}{k π} .$

Entonces por el teorema de Plancherel,

$\begin{aligned} \sum_{k = 1}^{\infty} \frac{4}{k^{2} π^{2}} & = \frac{1}{π} \int_{- π}^{π} f^{2} (x) d x \\ = \frac{1}{π} (\int_{- π}^{0} {(1 + \frac{x}{π})}^{2} d x + \int_{0}^{π} {(1 - \frac{x}{π})}^{2} d x) \\ = \frac{2}{3}, \end{aligned}$

teniendo que $\sum_{k = 1}^{\infty} \frac{1}{k^{2}} = \frac{2}{3} \frac{π^{2}}{4} = \frac{π^{2}}{6} .$

Ahora para obtener la otra identidad de Euler, notemos que

$\begin{aligned} \sum_{n = 0}^{\infty} \frac{1}{(2 n + 1)^{2}} & = \sum_{n = 1}^{\infty} \frac{1}{n^{2}} - \sum_{n = 1}^{\infty} \frac{1}{(2 n)^{2}} \\ = \frac{π^{2}}{6} - \frac{π^{2}}{4 \cdot 6} = \frac{π^{2}}{8} . \end{aligned}$

$△$

Entradas relacionadas

Ir a Álgebra Lineal I
Entrada anterior del curso: Proceso de Gram-Schmidt
Siguiente entrada del curso: Transformaciones multilineales

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104721 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM»

Álgebra Lineal I: Proceso de Gram-Schmidt

Por Blanca Radillo

4 respuestas

Introducción

Durante esta semana hemos introducido el concepto de bases ortogonales y ortonormales, así como algunas propiedades especiales. Para poder aplicar los resultados que hemos visto, es necesario insistir en que las bases sean de este tipo (ortonormales). Ahora veremos cómo encontrar bases ortonormales usando algo llamado el proceso de Gram-Schmidt.

Recordando todos los problemas anteriores de este curso, decíamos que una base es un conjunto de vectores linealmente independientes y que el número de vectores coincide con la dimensión del espacio. Pero hasta este momento no nos interesó determinar si las bases eran ortonormales o no. Si nos pusiéramos a ver si lo eran, es probable que muy pocas lo sean. Entonces surgen dos preguntas, ¿será difícil encontrar una base ortonormal de un espacio vectorial? y ¿habrá alguna manera de construir una base ortonormal?

Proceso de Gram-Schmidt

La respuesta a la primera pregunta es «no, no es difícil», y justo la respuesta de la segunda pregunta es la justificación. Dada una base cualquiera del espacio vectorial, podemos construir una base ortonormal de ese mismo espacio gracias al siguiente teorema.

Teorema (Gram-Schmidt). Sean $v_{1}, v_{2}, \dots, v_{d}$ vectores linealmente independientes en un espacio vectorial $V$ sobre $R$ (no necesariamente de dimensión finita), con producto interior $⟨ \cdot, \cdot ⟩$ . Entonces existe una única familia de vectores ortonormales $e_{1}, e_{2}, \dots, e_{d}$ en $V$ con la propiedad de que para todo $k = 1, 2, \dots, d$ , tenemos que

$\begin{aligned} span (e_{1}, e_{2}, \dots, e_{k}) & = span (v_{1}, v_{2}, \dots, v_{k}), y \\ ⟨ e_{k}, v_{k} ⟩ & > 0. \end{aligned}$

Demostración. Lo haremos por inducción sobre $d$ , la cantidad de vectores con la que empezamos.

La base inductiva es cuando $d = 1$ . Tomamos un vector $e_{1} \in span (v_{1})$ , entonces podemos escribirlo como $e_{1} = λ v_{1}$ para cierta $λ$ . Si queremos que $0 < ⟨ e_{1}, v_{1} ⟩ = λ {‖ v_{1} ‖}^{2}$ , entonces $λ > 0$ . Además queremos que $e_{1}$ tenga norma igual a 1, entonces $1 = {‖ e_{1} ‖}^{2} = ⟨ e_{1}, e_{1} ⟩ = λ^{2} {‖ v_{1} ‖}^{2},$ lo cual es posible si $λ = \frac{1}{‖ v_{1} ‖}$ . Como $e_{1}$ es un múltiplo escalar de $v_{1}$ , se tiene que $span (e_{1}) = span (v_{1})$ . Además, la construcción forzó a que $e_{1} = \frac{1}{‖ v_{1} ‖} v_{1}$ sea el único vector que satisface las condiciones del teorema.

Hagamos ahora el paso inductivo. Tomemos un entero $d \geq 2$ , y supongamos que el teorema es cierto para $d - 1$ . Sean $v_{1}, v_{2}, \dots, v_{d}$ vectores en $V$ linelmente independientes. Por hipótesis, sabemos que existe una única familia de vectores ortonormales $e_{1}, \dots, e_{d - 1}$ que satisfacen las condiciones del teorema respecto a la familia $v_{1}, \dots, v_{d - 1}$ . Es suficiente con probar que existe un único vector $e_{d}$ tal que $e_{1}, \dots, e_{d}$ satisface el teorema con respecto a $v_{1}, \dots, v_{d}$ , esto es
$\begin{aligned} ‖ e_{d} ‖ & = 1, \\ ⟨ e_{d}, e_{i} ⟩ & = 0 \forall 1 \leq i \leq d - 1, \\ ⟨ e_{d}, v_{d} ⟩ & > 0, \end{aligned}$

$span (e_{1}, \dots, e_{d}) = span (v_{1}, \dots, v_{d}),$

ya que, por hipótesis, los casos de $k < d$ se cumplen.

La idea para construir $e_{d}$ es tomarlo de $span (v_{1}, \dots, v_{d})$ , expresarlo como combinación lineal de estos y encontrar condiciones necesarias y suficientes sobre los coeficientes de $e_{d}$ para que satisfaga las conclusiones del teorema. Hagamos esto.

Sea $e_{d}$ un vector tal que $e_{d} \in span (v_{1}, \dots, v_{d})$ . Por ser linealmente independientes y por hipótesis $span (v_{1}, \dots, v_{d}) = span (e_{1}, \dots, e_{d - 1}) + span (v_{d}),$ entonces podemos escribir $e_{d}$ como

$e_{d} = λ v_{d} + \sum_{i = 1}^{d - 1} a_{i} e_{i}$

para algunos $λ, a_{1}, \dots, a_{d - 1}$ . Si resulta que $λ \neq 0$ , esto también implicará que $span (e_{1}, \dots, e_{d}) = span (v_{1}, \dots, v_{d})$ .

Ahora, dado que $e_{d}$ debe formar una familia ortonormal con el resto de los vectores, para todo $j = 1, \dots, d - 1$ , tenemos que

$\begin{aligned} 0 & = ⟨ e_{d}, e_{j} ⟩ \\ = λ ⟨ v_{d}, e_{j} ⟩ + \sum_{i = 1}^{d - 1} a_{i} ⟨ e_{i}, e_{j} ⟩ \\ = λ ⟨ v_{d}, e_{j} ⟩ + a_{j}, \end{aligned}$

entonces $a_{j} = - λ ⟨ v_{d}, e_{j} ⟩$ . Si logramos mostrar que hay un único $λ$ con el que se pueda satisfacer la conclusión del teorema, el argumento anterior muestra que también hay únicos $a_{1}, \dots, a_{d - 1}$ y por lo tanto que hay un único vector $e_{d}$ que satisface el teorema.

Sustituyendo los coeficientes anteriores, obtenemos que

$e_{d} = λ (v_{d} - \sum_{i = 1}^{d - 1} ⟨ v_{d}, e_{i} ⟩ e_{i}) .$

Notemos que si $z := v_{d} - \sum_{i = 1}^{d - 1} ⟨ v_{d}, e_{i} ⟩ e_{i}$ es cero, $v_{d}$ estaría en $span (e_{1}, \dots, e_{d - 1}) = span (v_{1}, \dots, v_{d - 1}),$ contradiciendo que los vectores $v_{i}$ ’s son linealmente independientes, entonces $z \neq 0$ .

Ahora como queremos que $1 = ‖ e_{d} ‖ = | λ | ‖ z ‖$ , esto implica que $| λ | = \frac{1}{‖ z ‖}$ .

Como además queremos que $⟨ e_{d}, v_{d} ⟩ > 0$ y

$⟨ e_{d}, v_{d} ⟩ = ⟨ e_{d}, \frac{e_{d}}{λ} + \sum_{i = 1}^{d - 1} ⟨ v_{d}, e_{i} ⟩ e_{i} ⟩ = \frac{1}{λ},$

se deduce que $λ$ es único y está determinado por $λ = \frac{1}{‖ z ‖} .$ Por lo tanto existe (y es único) el vector $e_{d}$ que satisface el teorema.

Este proceso de construcción es mejor conocido como el proceso de Gram-Schmidt. La demostración da a la vez un algoritmo que nos permite encontrar bases ortogonales (y de hecho ortonormales). Veremos ejemplos de esto en la siguiente sección. Antes de eso, enunciaremos formalmente una de las conclusiones más importantes del teorema anterior.

Recuerda que un espacio Euclideano es un espacio vectorial de dimensión finita sobre $R$ y con un producto interior. Podemos aplicar el proceso de Gram-Schmidt a cualquier base $v_{1}, \dots, v_{d}$ de un espacio Euclideano $V$ y al final obtendremos una familia $e_{1}, \dots, e_{d}$ de vectores ortonormales. Como sabemos que las familias de vectores ortonormales son linealmente independientes, y tenemos $d$ vectores, concluimos que $e_{1}, \dots, e_{d}$ es una base ortonormal. En resumen, tenemos el siguiente resultado.

Corolario. Todo espacio Euclideano tiene una base ortonormal.

Ejemplos de aplicación del proceso de Gram-Schmidt

A continuación veremos algunos ejemplos que nos ayuden a clarificar más este algoritmo.

Ejemplo 1. Sean $v_{1}, v_{2}, v_{3}$ vectores en $R^{3}$ (con el producto interior estándar) definidos por

$v_{1} = (1, 1, 0), v_{2} = (1, 1, 1), v_{3} = (1, 0, 1)$ .

Es fácil ver que estos vectores son linealmente independientes. Entonces construyamos según el proceso de Gram-Schmidt la familia ortonormal de vectores $e_{1}, e_{2}, e_{3}$ . Tenemos que

$e_{1} = \frac{v_{1}}{‖ v_{1} ‖} = \frac{v_{1}}{\sqrt{2}} = (\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0)$ .

Ahora, tomando $z_{2} = v_{2} - ⟨ v_{2}, e_{1} ⟩ e_{1}$ , tenemos que $e_{2}$ está definido como $\frac{z_{2}}{‖ z_{2} ‖}$ , entonces

$\begin{aligned} z_{2} & = (1, 1, 1) - [(1, 1, 1) \cdot (\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0)] (\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0) \\ = (1, 1, 1) - [\frac{2}{\sqrt{2}}] (\frac{1}{\sqrt{2}}, \frac{1}{\sqrt{2}}, 0) \\ = (1, 1, 1) - (2 / 2, 2 / 2, 0) \\ = (1, 1, 1) - (1, 1, 0) = (0, 0, 1) . \end{aligned}$

Esto implica que $e_{2} = \frac{1}{1} (0, 0, 1) = (0, 0, 1)$ . Finalmente tomando $z_{3} = v_{3} - ⟨ v_{3}, e_{1} ⟩ e_{1} - ⟨ v_{3}, e_{2} ⟩ e_{2}$ , sabemos que $e_{3} = \frac{z_{3}}{‖ z_{3} ‖}$ . Entonces

$\begin{aligned} z_{3} & = v_{3} - ⟨ v_{3}, e_{1} ⟩ e_{1} - ⟨ v_{3}, e_{2} ⟩ e_{2} \\ = (1, 0, 1) - (\frac{1}{2}, \frac{1}{2}, 0) - (0, 0, 1) \\ = (\frac{1}{2}, - \frac{1}{2}, 0) . \end{aligned}$

Por lo tanto

$e_{3} = \frac{1}{\sqrt{1 / 2}} (\frac{1}{2}, - \frac{1}{2}, 0) = (\frac{1}{\sqrt{2}}, \frac{- 1}{\sqrt{2}}, 0) .$

$△$

Ejemplo 2. Sea $V$ el espacio de polinomios en $[0, 1]$ con coeficientes reales de grado a lo más 2, con el producto interior

$⟨ p, q ⟩ = \int_{0}^{1} p (x) q (x) d x .$

Sean $v_{1} = 1$ , $v_{2} = 1 + x$ , $v_{3} = 1 + x^{2}$ vectores en $V$ que claramente son linealmente independientes. Encontraremos los vectores que nos da el proceso de Gram-Schmidt.

Primero calculemos

${‖ v_{1} ‖}^{2} = \int_{0}^{1} 1 d x = 1$ ,

entonces $e_{1} = \frac{v_{1}}{‖ v_{1} ‖} = v_{1} = 1$ . Ahora calculemos $z_{2}$ :

$\begin{aligned} z_{2} & = v_{2} - ⟨ v_{2}, e_{1} ⟩ e_{1} \\ = 1 + x - \int_{0}^{1} (1 + x) d x = 1 + x - (1 + \frac{1}{2}) \\ = x - \frac{1}{2} . \end{aligned}$

Haciendo la integral $\int_{0}^{1} {(x - \frac{1}{2})}^{2} d x$ se obtiene que $‖ z_{2} ‖ = \sqrt{\frac{1}{12}}$ , entonces $e_{2} = \sqrt{12} (x - \frac{1}{2})$ .

Por último, hay que calcular $z_{3}$ así como su norma. Primero,

$\begin{aligned} z_{3} & = v_{3} - ⟨ v_{3}, e_{1} ⟩ e_{1} - ⟨ v_{3}, e_{2} ⟩ e_{2} \\ = (1 + x^{2}) - \int_{0}^{1} (1 + x^{2}) d x - 12 (x - \frac{1}{2}) \int_{0}^{1} (1 + x^{2}) (x - \frac{1}{2}) d x \\ = 1 + x^{2} - (1 + \frac{1}{3}) - 12 (x - \frac{1}{2}) (\frac{1}{12}) \\ = x^{2} - \frac{1}{3} - x + \frac{1}{2} \\ = x^{2} - x + \frac{1}{6}, \end{aligned}$

y luego, con la integral $\int_{0}^{1} {(x^{2} - x + \frac{1}{6})}^{2} d x$ se calcula que $‖ z_{3} ‖ = \frac{1}{6 \sqrt{5}}$ , por lo tanto $e_{3} = 6 \sqrt{5} (x^{2} - x + \frac{1}{6})$ .

$△$

Aunque no es un proceso muy eficiente, nos garantiza que podemos encontrar una base ortonormal para cualquier espacio vectorial (con producto interior). Ya con una base ortonormal, podemos usar la descomposición de Fourier de la cual hablamos la entrada anterior y con ella todas las consecuencias que tiene.

Si quieres ver muchos más ejemplos del proceso en $R^{n}$ , puedes usar una herramienta en línea que te permite ver el proceso paso a paso en el conjunto de vectores que tu elijas. Una posible página es el Gram-Schmid Calculator de eMathHelp.

Más adelante…

En esta última entrada teórica de la unidad 3, vimos el método de Gram-Schmidt para construir una base ortonormal, que es un proceso algorítmico que parte de tener una base de un espacio y al final calcula una base ortonormal. También se vieron algunos ejemplos de la aplicación de este proceso para espacios vectoriales finitos como $R^{3}$ y el espacio de polinomios en [0,1] de grado a lo más 2. Aunque no es una manera muy eficaz para encontrar una base ortonormal, sí te garantiza que lo que construye es una.

En la próxima entrada veremos ejercicios resueltos de los temas que hemos estado estudiando a lo largo de esta semana.

Tarea moral

A continuación hay algunos ejercicios para que practiques los conceptos vistos en esta entrada. Te será de mucha utilidad intentarlos para entender más la teoría vista.

Verifica que con el valor $λ$ que se encontró en la demostración del teorema de Gram-Schmidt en efecto se obtiene un vector $e_{d}$ que satisface todas las conclusiones que se desean.
Revisa que los vectores que se obtuvieron en los ejemplos de aplicación del proceso de Gram-Schmidt en efecto son bases ortogonales de los espacios correspondientes.
Aplica el proceso de Gram-Schmidt a los polinomios $1$ , $x$ , $x^{2}$ en el espacio Euclideano de los polinomios reales de grado a lo más dos y producto interior $⟨ p, q ⟩ = p (0) q (0) + p (1) q (1) + p (2) q (2) .$
Aplica el proceso de Gram-Schmidt a los vectores $\begin{array}{r} (1, 1, 1, 1) \\ (0, 1, 1, 1) \\ (0, 0, 1, 1) \\ (0, 0, 0, 1) \end{array}$ de $R^{4}$ con el producto interior canónico (el producto punto).
Usa el Gram-Schmidt Calculator de eMathHelp para ver paso a paso cómo se aplica el proceso de Gram-Schmidt a los vectores $\begin{array}{r} (1, 2, 1, 1, - 1) \\ (0, 0, 1, 0, 0) \\ (2, 0, 0, 1, 1) \\ (0, 2, 0, 0, 1) \\ (- 3, 0, 0, 1, 0) \end{array}$ de $R^{5}$ .

Entradas relacionadas

Ir a Álgebra Lineal I
Entrada anterior del curso: Bases ortonormales y descomposición de Fourier
Siguiente entrada del curso: Problemas de bases ortogonales, Fourier y proceso de Gram-Schmidt

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104721 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM»

El blog de Leo

Aprendiendo, creando y compartiendo matemáticas

Archivo del Autor: Blanca Radillo

Álgebra Lineal I: Aplicaciones del teorema espectral, bases ortogonales y más propiedades de transformaciones lineales

Introducción

Matrices simétricas, matrices diagonalizables

Más propiedades de transformaciones lineales y bases ortogonales

Entradas relacionadas

Agradecimientos

Álgebra Lineal I: Problemas de determinantes y ecuaciones lineales

Introducción

Problemas de ecuaciones lineales

Entradas relacionadas

Agradecimientos

Álgebra Lineal I: Problemas de definición y propiedades de determinantes

Introducción

Problemas de transformaciones antisimétricas

Problemas de determinantes

Entradas relacionadas

Agradecimientos

Álgebra Lineal I: Problemas de bases ortogonales, Fourier y proceso de Gram-Schmidt

Introducción

Problemas resueltos de bases ortogonales y proyecciones

Distancia de un vector a un subespacio y desigualdad de Bessel

Aplicación del proceso de Gram-Schmidt

El teorema de Plancherel y una fórmula con $π$

Entradas relacionadas

Agradecimientos

Álgebra Lineal I: Proceso de Gram-Schmidt

Introducción

Proceso de Gram-Schmidt

Ejemplos de aplicación del proceso de Gram-Schmidt

Más adelante…

Tarea moral

Entradas relacionadas

Agradecimientos

Introducción

Matrices simétricas, matrices diagonalizables

Más propiedades de transformaciones lineales y bases ortogonales

Entradas relacionadas

Agradecimientos

Introducción

Problemas de ecuaciones lineales

Entradas relacionadas

Agradecimientos

Introducción

Problemas de transformaciones antisimétricas

Problemas de determinantes

Entradas relacionadas

Agradecimientos

Introducción

Problemas resueltos de bases ortogonales y proyecciones

Distancia de un vector a un subespacio y desigualdad de Bessel

Aplicación del proceso de Gram-Schmidt

El teorema de Plancherel y una fórmula con π

Entradas relacionadas

Agradecimientos

Introducción

Proceso de Gram-Schmidt

Ejemplos de aplicación del proceso de Gram-Schmidt

Más adelante…

Tarea moral

Entradas relacionadas

Agradecimientos

El teorema de Plancherel y una fórmula con $π$