Archivo de la etiqueta: fourier

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

Por Blanca Radillo

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=WW2, para todo vV podemos definir su proyección en W, que denotamos πW(v), como el único elemento en W tal que vπW(v)W2.

Debido a las discusiones sobre bases ortogonales, no es difícil ver que si w,u=0 para todo wW, entonces uW2. Como consecuencia de esto, tenemos el siguiente resultado:

Teorema. Sea V un espacio vectorial sobre R con producto interior ,, y sea W un subespacio de V de dimensión finita. Sea v1,,vn una base ortogonal de W. Entonces para todo vV tenemos que

πW(v)=i=1nv,vivi2vi.

Demostración. Escribimos v como v=πW(v)+u con uW2. Por la observación previa al teorema, u,vi=0 para todo i. Además existen a1,,an tales que πW(v)=a1v1++anvn. Entonces

0=u,vi=v,viπW(v),vi=v,vij=1najvj,vi=v,viaivi,vi,

porque v1,,vn es una base ortogonal. Por lo tanto, para todo i, obtenemos

ai=v,vivi2.

◻

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)=minxWxv.

Dado que xW, xπW(v)W, y por definición de proyección vπW(v)W2, entonces

xv2=(xπW(v))+(πW(v)v)2=xπW(v)2+2xπW(v),πW(v)v+πW(v)v2=xπW(v)2+πW(v)v2πW(v)v2.

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 ,, y sea W un subespacio de V de dimensión finita. Sea v1,,vn una base ortonormal de W. Entonces para todo vV tenemos que

πW(v)=i=1nv,vivi,

y

d(v,W)2=vi=1nv,vivi2=v2i=1nv,vi2.

En particular

i=1nv,vi2v2.

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

Demostración. Por el teorema anterior y dado que v1,,vn es una base ortonormal, obtenemos la primera ecuación. Ahora, por Pitágoras,

d(v,W)2=vπW(v)2=v2πW(v)2.

Por otro lado, tenemos que

πW(v)2=i=1nv,vivi2=i,j=1nv,vivi,v,vjvj=i,j=1nv,viv,vjvi,vj=i=1nv,vi2.

Por lo tanto, se cumple la igualdad de la distancia. Finalmente como d(v,W)20, 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=01xp(x)q(x)dx.

Aplica el algoritmo de Gram-Schmidt a los vectores 1,x,x2.

Solución. Es fácil ver que ese sí es un producto interior en V (tarea moral). Nombremos v1=1,v2=x,v3=x2. Entonces

e1=v1v1=2v1=2,

ya que v12=01xdx=12.

Sea z2=v2v2,e1e1. Calculando, v2,e1=012x2dx=23. Entonces z2=x232=x23. Esto implica que

e2=z2z2=6(x23)=6x4.

Finalmente, sea z3=v3v3,e1e1v3,e2e2. Haciendo los cálculos obtenemos que

z3=x2(24)2(15)(6x4)

z3=x265x+310.

Por lo tanto

e3=z3z3=106(x265x+310).

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:RR definida como f(0)=f(π)=0, f(x)=1xπ en el intervalo (π,0), y f(x)=1xπ 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

n=11n2=π26,n=01(2n+1)2=π28.

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

a0(f)=1πππf(x)dx,

ak(f)=1πππf(x)cos(kx)dx,

bk(f)=1πππf(x)sen(kx)dx.

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 k1, a0=0,ak=0,bk=2kπ.

Entonces por el teorema de Plancherel,

k=14k2π2=1πππf2(x)dx=1π(π0(1+xπ)2dx+0π(1xπ)2dx)=23,

teniendo que k=11k2=23π24=π26.

Ahora para obtener la otra identidad de Euler, notemos que

n=01(2n+1)2=n=11n2n=11(2n)2=π26π246=π28.

Entradas relacionadas

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: Bases ortogonales

Por Leonardo Ignacio Martínez Sandoval

Introducción

Como ya discutimos en las entradas anteriores, si tenemos un espacio vectorial V con producto interior, entonces podemos definir varias nociones geométricas en V, como ángulos, norma y distancia. Ahora vamos a definir una noción muy útil en álgebra lineal: la de bases ortogonales. Para ello, combinaremos las nociones de bases y producto interior.

Las bases ortogonales no sólo tienen aplicaciones en álgebra lineal. También son el punto de partida de muchos conceptos matemáticos avanzados. Un primer ejemplo es el análisis de Fourier, que estudia cómo aproximar funciones mediante funciones trigonométricas y que tiene aplicaciones en el mundo real en análisis de señales. Otro ejemplo es la vasta teoría de polinomios ortogonales, con aplicaciones en el mundo real en aproximación e integración numérica.

En estas entradas de bases ortogonales tomaremos espacios vectoriales sobre R con un producto interior ,.

Conjuntos ortogonales y ortonormales

Comenzamos con la siguiente definición. Recuerda que V es un espacio vectorial sobre R con producto interior, así que induce una norma .

Definición. Sea S un conjunto de vectores en V. Decimos que S es

  • Ortogonal si cualquier par de vectores distintos de S es ortogonal, es decir, si para todo v,w en S, con vw se tiene que v,w=0.
  • Ortonormal si es ortogonal, y además todo vector de S tiene norma 1.

En otras palabras, S es ortonormal si para todo v en S se tiene v,v=1 y para v y w en S distintos se tiene v,w=0.

Ejemplo. Si tomamos a Rn con el producto punto, entonces la base canónica es un conjunto ortonormal pues, en efecto, eiei=1 y para ij se tiene eiej=0.

Todo conjunto de un sólo elemento es ortogonal, pues no hay nada que probar. Otro conjunto ortonormal en R2 es el conjunto que sólo tiene al vector (35,45), pues este es un vector de norma 1.

Los vectores (1,1,0), (1,1,0) y (0,0,1) forman otro conjunto ortogonal en R3, pues en efecto
(1,1,0)(1,1,0)=11=0(1,1,0)(0,0,1)=0(0,0,1)(1,1,0)=0.

Sin embargo, este no es un conjunto ortonormal, pues la norma de (1,1,0) es 21. Si normalizamos a cada vector, es decir, si lo dividimos entre su norma, entonces obtenemos los vectores ortonormales (1/2,1/2,0), (1/2,1/2,0) y (0,0,1).

Propiedades de conjuntos ortogonales y ortonormales

Todo conjunto ortogonal de vectores no nulos se puede normalizar como en el ejemplo de la sección anterior para obtener un conjunto ortonormal. Es decir, si S es un conjunto de vectores distintos de 0, entonces S={vv:vS} es un conjunto ortonormal.

Una propiedad fundamental de los conjuntos ortonormales de vectores es que son linealmente independientes. Se puede probar algo un poco más general.

Proposición. Si S es un conjunto ortogonal de vectores no nulos, entonces los elementos de V son linealmente independientes.

Demostración. Tomemos v1,,vn elementos de S y supongamos que existen α1,,αn escalares tales que v:=i=1nαivi=0.

Tomemos un índice j en 1,,n y hagamos el producto interior v,vj. Por un lado, como v=0, este produto es 0. Por otro lado, por linealidad es i=1nαivi,vj.

Cuando ij, el sumando correspondiente es igual a 0. De este modo, el único sumando no cero es cuando i=j, el cual es αjvj,vj. De estos argumentos, deducimos que αjvj,vj=0. Como los vectores son no nulos, se tiene que vj,vj0. Así, αj=0 para todo j=1,,n, lo cual muestra que los vectores son linealmente independientes.

◻

Como cada elemento de un conjunto ortonormal tiene norma 1, entonces no puede ser nulo, así que como corolario de la proposición anterior, todo conjunto ortonormal es linealmente independiente. Otro corolario es el siguiente.

Corolario. En un espacio Euclideano de dimensión d, los conjuntos ortogonales sin vectores nulos tienen a lo más d elementos.

Bases ortogonales y ortonormales

Cuando una base de un espacio vectorial es ortogonal (o bien, ortonormal), pasan varias cosas buenas. Esto amerita una definición por separado.

Definición. Sea S un conjunto de vectores en V. Decimos que S es

  • Una base ortogonal si S es una base de V y es un conjunto ortogonal.
  • Una base ortonormal si S una base de V y es un conjunto ortonormal.

Ejemplo. En Rn la base canónica es una base ortonormal.

En R2 el conjunto S={(2,3),(9,6)} es un conjunto ortogonal. Además, se puede verificar fácilmente que son dos vectores linealmente independientes. De este modo, S es una base ortogonal.

Sin embargo, S no es una base ortonormal pues el primero de ellos tiene norma 22+32=13. Si quisiéramos convertir a S en una base ortonormal, podemos normalizar a cada uno de sus elementos.

En la sección anterior vimos que los conjuntos ortonormales son linealmente independientes. Otro corolario de este resultado es lo siguiente.

Corolario. En un espacio Euclideano de dimensión n, un conjunto ortonormal de n vectores es una base ortonormal.

La importancia de las bases ortogonales yace en que dada una base ortonormal B y un vector v, podemos encontrar varias propiedades de v en términos de B fácilmente. Por ejemplo, veremos más adelante que:

  • Las coordenadas de v con respecto a la base B son sencillas.
  • Hay una fórmula simple para la norma de v en términos de sus coordenadas en la base B.
  • Si B es una base de un subespacio W de V, entonces es fácil encontrar la distancia de v a W.

Mejor aún, las bases ortonormales siempre existen.

Teorema. Todo espacio Euclideano tiene una base ortonormal.

Es decir, sin importar qué espacio vectorial real de dimensión finita tomemos, y sin importar qué producto punto le pongamos, podemos dar una base ortogonal. De hecho, veremos un resultado un poco más fuerte, que nos dará un procedimiento para encontrar dicha base, incluso imponiendo restricciones adicionales.

Ejemplo de bases ortogonales en polinomios

Ejemplo. Tomemos Rn[x] el espacio de polinomios de grado a lo más n con coeficientes reales. Además, tomemos números reales distintos x0,,xn. A partir de estos reales podemos definir la operación P,Q=j=0nP(xj)Q(xj), la cual es claramente bilineal y simétrica.

Tenemos que P,P es una suma de cuadrados, y por lo tanto es no negativa. Además, si P,P=0, es porque j=0nP(xj)2=0, y como estamos trabajando en R esto implica que cada sumando debe ser cero. Pero las igualdades P(x0)==P(xn)=0 dicen que los n+1 reales distintos xi son raíces de P, y como P es de grado a lo más n, tenemos que P es el polinomio 0. En resumen, , es un producto interior en Rn[x]. Vamos a dar una base ortogonal con respecto a este producto interior.

Para i=0,,n, consideremos los polinomios Li(x)=0kn,kixxkxixk. Observa que Lj(xj)=1 y si ji, tenemos Li(xj)=0. Afirmamos que B={Lj:j=0,,n+1} es una base ortonormal de Rn[x] con el producto interior que definimos. Como consiste de n+1 polinomios y dim(Rn[x])=n+1, basta con que veamos que es un conjunto ortonormal.

Primero, notemos que
Li,Li=j=0nLi(xj)2=Li(xi)2=1,

de modo que cada Li tiene norma 1.

Luego, notemos que si ij, entonces Li(xk)Lj(xk)=0 pues xk no puede ser simultáneamente xi y xj. De este modo,

Li,Lj=k=0nLi(xk)Lj(xk)=0.

Con esto mostramos que cada par de polinomios distintos es ortogonal. Esto termina la demostración de que B es base ortonormal.

◻

Ejemplo de conjuntos ortogonales en funciones periódicas

Ejemplo. Consideremos V el conjunto de funciones f:RR continuas y periódicas de periodo 2π. Definimos f,g=ππf(x)g(x)dx. Se puede mostrar que , así definido es un producto interior en V.

Para cada entero positivo n, definimos
Cn(x)=cos(nx)πSn(x)=sin(nx)π.

Además, definimos C0(x)=12π. Afirmamos que F:={Cn:n0}{Sn:n1} es un conjunto ortonormal de vectores. Mostremos esto.

Para empezar, notamos que C02=ππ12πdx=1.

Luego, tenemos que para n1 que
Cn2=ππ1πcos2(nx)dx=ππ1+cos(2nx)2πdx=1,

ya que para todo entero m0 se tiene que ππcos(mx)dx=0. De manera similar, usando la identidad sin2(nx)=1cos(nx)2, se puede ver que la norma de Sn es 1.

Para ver que las parejas de elementos distintas son ortogonales, tenemos varios casos. Si tomamos n1, el resultado para C0,Cn ó C0,Sn se deduce de que
ππcos(mx)dx=ππsin(mx)dx=0 para todo entero m0.

Si tomamos dos Ci’s distintos, dos Sis distintos o un Ci y un Si, el resultado se deduce de las fórmulas «producto a suma» de las funciones trigonométricas.

◻

Más adelante…

En esta entrada combinamos las nociones de bases y el producto interior, estudiadas en entradas anteriores, para definir a las bases ortogonales. Vimos algunas propiedades de conjuntos ortogonales y ortonormales, para extenderlos a bases ortogonales y ortonormales. Vimos unos ejemplos de bases ortogonales de los polinomios y otros ejemplos de conjuntos ortogonales en funciones periódicas.

En la siguiente entrada veremos aplicaciones de estos conceptos, culminando en una descomposición de Fourier.

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.

  • Encuentra un conjunto ortogonal de vectores en R4 tal que ninguna de las entradas de ninguno de sus vectores sea igual a 0.
  • Escribe las demostraciones de los corolarios enunciados en esta entrada.
  • Muestra que , definido en el ejemplo de funciones periódicas es un producto interior.
  • Termina de mostrar que la familia F del ejemplo de funciones periódicas es ortonormal. Sugerencia: Usa identidades de suma y resta de ángulos para poner el producto de senos (o cosenos o mixto) como una suma de senos y/o cosenos.

Entradas relacionadas

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»