Archivo de la categoría: Matemáticas

Posts de matemáticas, la ciencia más cercana a las artes.

Geometría Analítica I: Diagonalización ortogonal de matrices simétricas

Por Paola Lizeth Rojas Salazar

Introducción

Anteriormente, estudiamos los vectores y valores propios de las matrices simétricas, en esta entrada vamos a usar que ya sabemos muchas cosas sobre el comportamiento respecto al producto interior, para hablar sobre la diagonalización ortogonal de matrices simétricas, cuyo procedimiento inicia resolviendo su polinomio característico.

Teoremas importantes

Antes de ver el proceso para la diagonalización ortogonal de matrices simétricas, vamos a enunciar un lema y un teorema que van a justificar la «receta» a seguir para esta diagonalización.

Lema 4.12: Considera una matriz simétrica A. Si λ1,u y λ2,v, son pares propios de A con λ1λ2, entonces u y v son ortogonales.

Demostración

Sabemos que:

(1)λ1(uv)=(λ1u)v=Auv=uAv=u(λ2v)=λ2(uv)

Esto implica que (λ1λ2)(uv)=0

Y λ1λ2, entonces uv=0.

Con lo que hemos terminado la demostración.

Teorema 4.13: Considera una matriz simétrica de 2×2, A. Entonces existe una rotación BO(2) tal que BTAB es diagonal de la siguiente forma:

(2)(λ100λ2)

Con λ1 y λ2, los valores propios de A.

Demostración

Por las entradas anteriores, las siguientes implicaciones son ciertas, puedes comprobarlo tú mismo con facilidad.

Como A es simétrica de 2×2, entonces A tiene valores propios λ1,λ2R.

Caso 1 λ1=λ2

Entonces A es diagonal y puede tomarse a B como la matriz identidad que es rotación en O(2).

Caso 2 λ1λ2

Consideramos a u,v, los vectores propios correspondientes a λ1 y λ2. Observa que u es diferente al vector cero.

Sabemos que u y v son ortogonales, entonces v es paralelo a uT que también es vector propio correspondiente a λ2.

Considera B=1|u|(u,uT), donde se puede comprobar fácilmente que B es la matriz de una rotación y que cumple que BTAB es diagonal.

«Receta»

Ingredientes

  1. Una matriz simétrica A=AT de 2×2

Procedimiento

  1. Resolver su polinomio característico con det(AλI).
  2. Encontrar u0 tal que (Aλ1I)u=0.
  3. Declarar B=1|u|(u,uT).
  4. La matriz diagonal, con entradas λ1 y λ2, estará dada por BTAB.

Tarea moral

  1. Termina de escribir la demostración del Teorema 4.13.
  2. Demuestra que, si una matriz A cualquiera, tiene dos valores propios distintos, entonces existe una matriz BGl(2) tal que B1AB es diagonal.
  3. Encuentra la matriz B de una rotación que diagonalice las siguientes matrices simétricas: Además, calcula BTAB:
    • (3)A=(1111)
    • (4)A=(612121)
    • (5)A=(7662)

Más adelante…

Avanza a las siguientes entradas, en las que usaremos estos conocimientos para dar dos nuevas formas de clasificación de las curvas.

Geometría Analítica I: Valores y vectores propios

Por Paola Lizeth Rojas Salazar

Introducción

Mucho hemos hablado, a lo largo de las entradas anteriores, sobre las soluciones de las ecuaciones resultantes de la forma det(AIλ)=0 con A matriz y x vector; sin embargo, aún no nos hemos dedicado a resolver este tipo de ecuaciones. En esta entrada, hablaremos de las soluciones de estas ecuaciones que se llaman valores propios y que tienen un vector propio asociado.

¿Qué son?

Si tienes una matriz cuadrada A, vamos a decir que un vector v es vector propio de A con valor propio λR, si v0 y Av=λv. A la pareja λ,v, la vamos a llamar pareja propia de A y el principal problema de estudio de esta entrada será encontrar a estas parejas.

Implicaciones importantes

Lema 4.7: Si u y v son vectores propios de A con el mismo valor λR, entonces cualquier combinación lineal no trivial de ellos también es vector propio de A.

Demostración

Si tenemos dos parejas de A, u, v que cumplen Au=λu y Av=λv, entonces, para cualquier par de coeficientes α,βR se tiene que, la combinación lineal de u y v con estos vectores, se cumple:

(6)A(αu+βv)=α(Au)+β(Av)=α(λu)+β(λv)=λ(αu+βv)

Lo que significa que, si αu+βv0, entonces es vector propio de A con valor propio λ.

Con lo que hemos terminado la demostración.

El siguiente lema es muy importante para realizar los cálculos.

Lema 4.8: Para cualquier matriz A cuadrada, se cumple: λ es un valor propio de A si y solo si, det(AλI)=0

Demostración

Sabemos que, como λ es valor propio de A, entonces tiene a su vector propio v correspondiente que cumple que: Av=λv.

Avλv=0

(AλI)v=0

Y, como v0, entonces:

det(AλI)=0

Fin de la demostración.

Observa que, con los conocimientos que tenemos hasta el momento, ya puedes demostrar fácilmente el siguiente Lema y Corolario.

Lema 4.9: Si A es una matriz simétrica de 2×2, entonces A tiene dos valores propios λ1 y λ2 en R.

Corolario 4.10: Sea A una matriz simétrica de 2×2, con valores propios λ1 y λ2. Entonces, sus valores propios coinciden (λ1=λ2) si y solo si, A=λ1I. En este caso, cualquier vector v0 es vector propio.

Corolario 4.11: Considera una matriz simétrica A de 2×2. Entonces, existe una base u,vR2, donde u y v son vectores propios de A.

Demostración

Por el Lema 4.9, sabemos que A tiene dos valores propios λ1 y λ2.

Caso 1, λ1=λ2

Por el Corolario 4.10, cualquier base u,vR2 funciona.

Caso 2, λ1λ2

Por la definición de valor propio, existen u,v, vectores distintos de 0 que corresponden a los valores propios λ1 y λ2 respectivamente.

Estos vectores no pueden ser paralelos porque por el Lema 4.8, esto implicaría que λ1=λ2.

Entonces, u y v forman una base de R2.

Terminamos la demostración.

Ejemplo

Calculemos los valores y vectores propios de la siguiente matriz simétrica:

(7)A=(2221)

Valores propios

Recordemos que, para encontrar los valores propios, debemos resolver su polinomio característico que está dado por det(AλI):

(8)det(AλI)=det((2221)(λ00λ))

Si continuamos con el desarrollo de la expresión anterior, comprueba que llegamos al siguiente polinomio característico:

(9)det(AλI)=λ2λ6

Para resolver el polinomio anterior, debemos igualarlo a 0, de donde vamos a obtener que, las raíces del polinomio son: λ1=1+1+242=3 y λ2=152=2

Vectores propios

Para encontrar los vectores propios correspondientes a λ1 y λ2, debemos encontrar una solución no trivial para los sistemas (AλiI)x=0 con i=1,2

Para λ1=3

(10)(A3I)x=(1224)(xy)=(x+2y2x4y)=(00)

De donde obtenemos el siguiente sistema:

(11)x+2y=02x4y=0

Donde, una de sus soluciones no triviales es uT=(2,1)

Para λ2=2

(12)(A(2)I)x=(4221)(xy)=(4x+2y2x+y)=(00)

De donde obtenemos el siguiente sistema:

(13)4x+2y=02x+y=0

Donde, una de sus soluciones no triviales es vT=(1,2)

Observa que estos vectores u y v son ortogonales, ¿será coincidencia? Lo veremos más adelante.

Tarea moral

  1. Comprueba que, para los vectores propios obtenidos en los sistemas de ecuaciones (6) y (7), se cumple que Au=3u y que Av=2v.
  2. Demuestra, con un argumento algebraico y uno geométrico, que la matriz (14)(0110) no tiene vectores propios.
  3. Demuestra que la matriz (15)(abba) no tiene vectores propios para b0.
  4. Usa el Lema 4.9 para demostrar el Corolario 4.10.
  5. Demuestra el Lema 4.9. Hint: usa que, al ser A matriz simétrica, entonces A=AT, después, expresa a A de la siguiente forma y desarrolla:

(16)A=(abbc)

Más adelante…

En la siguiente entrada, concluiremos nuestro estudio de los valores y vectores propios, analizando la diagonalización ortogonal de matrices simétricas.

Cálculo Diferencial e Integral I: La regla de la cadena

Por Juan Manuel Naranjo Jurado

Introducción

Anteriormente revisamos, entre otras cosas, cómo derivar la suma, el producto y el cociente de funciones. La siguiente operación a analizar es la composición de funciones, tema del cual tratará esta entrada.

Demostración de la regla de la cadena

Teorema. Sean A, BR, g:AR, f:BR y x0A tales que

  1. Para todo xA, g(x)B.
  2. g es derivable en x0, es decir limxx0g(x)g(x0)xx0=g(x0).
  3. f es derivable en g(x0), es decir limtx0f(t)f(g(x0))tg(x0)=f(g(x0)).

Entonces fg es derivable en x0, además (fg)(x0)=f(g(x0))g(x0).

Demostración.

Para realizar esta demostración haremos uso de una función auxiliar de la que probaremos propiedades específicas.

ρ(t)={f(t)f(g(x0))tg(x0)f(g(x0)), si tg(x0)0, si t=g(x0).

Podemos observar que la función ρ está «inspirada» en la definición de derivada de f en el punto g(x0). Procederemos a puntualizar 5 observaciones de nuestra función auxiliar.

  1. Como f:BR, entonces ρ:BR.
  2. El límite de ρ en g(x0) es cero, puesto que
    limtg(x0)ρ(t)=limtg(x0)(f(t)f(g(x0))tg(x0)f(g(x0)))=limtg(x0)f(t)f(g(x0))tg(x0)limtg(x0)f(g(x0))=f(g(x0))f(g(x0)), por el supuesto 3=0.

    limtg(x0)ρ(t)=0.
  3. ρ es continua en g(x0), puesto que limtg(x0)ρ(t)=0=ρ(g(x0)).
  4. Para todo tB, se sigue de la definición de ρ que f(t)f(g(x0))=(ρ(t)+f(g(x0))(tg(x0)).
  5. Por el supuesto 2, g es derivable en x0 lo que implica que también es continua en tal punto, además por la observación 3, sabemos que ρ es continua en g(x0). Por tanto, se tiene que
    limxx0ρ(g(x))=ρ(g(x0))=0.ρg es continua en x0.

Ahora que establecimos las 5 observaciones, estamos listos para calcular la derivada de la composición:

(fg)(x0)=limxx0f(g(x))f(g(x0))xx0=limxx0(ρ(g(x))+f(g(x0)))(g(x)g(x0))xx0, por la obs 4=limxx0((ρ(g(x))+f(g(x0)))g(x)g(x0)xx0)=limxx0(ρ(g(x))+f(g(x0)))limxx0g(x)g(x0)xx0=(0+f(g(x0)))g(x0), por la obs 3 y el supuesto 2=f(g(x0))g(x0).

(fg)(x0)=f(g(x0))g(x0).

◻

Aplicando la regla de la cadena

A continuación revisaremos algunos ejemplos donde aplicaremos la proposición anterior. La idea general de los ejercicios será expresar una función en términos de la composición de otras dos.

Ejemplo 1. Encuentra la derivada de la función F(x)=(3x+1)2.

Notemos que podemos ver a F como la composición de las siguientes dos funciones
f(x)=x2,g(x)=3x+1.

Así, F(x)=f(g(x)). Y empleando la regla de la cadena se tiene que

F(x)=f(g(x))g(x)=2g(x)g(x)=2(3x+1)(3)=6(3x+1)=18x+6.

Ejemplo 2. Deriva la función F(x)=x2+1x3+3.

Definimos las funciones

f(x)=x con derivada f(x)=12x y g(x)=x2+1x3+3 con derivada
g(x)=(x3+3)(2x)(x2+1)(3x2)(x3+3)2=x43x2+6x(x3+3)2.

Con lo anterior, se tiene que F(x)=f(g(x)), y empleando la regla de la cadena tenemos

F(x)=f(g(x))g(x)=12g(x)x43x2+6x(x3+3)2=12x2+1x3+3x43x2+6x(x3+3)2=x3+32x2+1x43x2+6x(x3+3)2.

Más adelante…

En las siguientes entradas haremos un resumen de las «reglas de derivación» que hemos visto hasta ahora y probaremos algunas más; particularmente se hará la revisión de las derivadas para las funciones trigonométricas.

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.

  • Teorema de Carathéodory. Sea f definida en un intervalo A y sea aA. Entonces f es derivable en a si y solo si existe una función ρ en A que es continua en a y satisface:
    f(x)f(a)=ρ(x)(xa) para xA.
    En este caso, se tiene que ρ(a)=f(a).
  • Deriva la función f(x)=52x+x2.
  • Si f:RR es derivable en x0 y f(x0)=0. Prueba que g(x):=|f(x)| es derivable en x0 si y solo si f(x0)=0.
  • Determina en dónde es derivable cada una de las siguientes funciones de RR y encuentra la derivada:
    • f(x)=|x|+|x+1|.
    • g(x)=2x+|x|.
    • h(x)=x|x|.

Entradas relacionadas

Agradecimientos

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

Geometría Analítica I: Encontrar el centro y los ejes de una cónica

Por Paola Lizeth Rojas Salazar

Introducción

En esta entrada, continuaremos con el estudio de las cónicas, pero en esta ocasión, vamos a encontrar su centro y ejes, a partir de dos grupos de isometrías que ya son familiares para nosotros, las rotaciones y traslaciones y usando otro tema que ya ha sido estudiado con anterioridad, la equivalencia de polinomios y reducción de términos lineales y cuadráticos.

Encontrando el centro de las traslaciones

Para cualquier vector hR2, consideremos la traslación g(x)=x+h y veamos cómo se escribe el polinomio Pg con P(x)=xAx+kx+f:

(17)(Pg)(x)=P(x+h)=(x+h)A(x+h)+k(x+h)+f

Factorizando y desarrollando un poco la expresión anterior, obtenemos:

(18)(Pg)(x)=xAx+xAh+hAx+kx+(hAh+kh+f)

Donde hAh+kh+f=P(h)

Lo que nos lleva, finalmente, a:

(19)(Pg)(x)=xAx+xAh+hAx+kx+P(h)

La pregunta ahora es, ¿hay una forma de encontrar el centro de la traslación g a partir de esta expresión? Lo que, por muy extraño que parezca, es cierto, pero, ¿cómo?

Enunciemos unos lemas que nos ayudarán a encontrar la respuesta a la pregunta anterior.

Lema 4.4: Dadas A y B dos matrices que se puedan multiplicar, se cumple que (AB)T=BTAT

Lema 4.5: Si tenemos una matriz simétrica A (recordemos que una matriz A es simétrica si A=AT), entonces, para todo par de vectores x,y en R, se cumple que xAy=Axy

Demostración

Sean x,y vectores, recordemos que x=xT y que y=yT, por esto y el lema anterior, tenemos que:

(20)xAy=xTAy=(xTAy)T=(Ay)T(xT)T=yTATx=yAx=Axy

Ahora sí podemos encontrar el centro de la traslación considerando:

(21)(Pg)(x)=xAx+xAh+hAx+kx+P(h)

Ya que, considerando el lema anterior, podemos simplificar (Pg)(x) de la siguiente manera:

(22)(Pg)(x)=xAx+xAh+hAx+kx+P(h)=2(Ahx)+kx+P(h)

Y, finalmente:

(23)(Pg)(x)=(2Ah+k)x+P(h)

Donde (2Ah+k)x es la parte lineal de esta composición por lo que, si podemos encontrar una hR2 que cumpla que 2Ah+k=0, entonces habremos encontrado una traslación que no contenga la parte lineal del polinomio. Si esta h existe, es el centro de la traslación (en el caso de este capítulo, estaremos hablando de traslaciones de cónicas).

Lo anterior lo podemos resumir en el siguiente lema:

Lema 4.6: Sea P(x)=xAx+kx+f un polinomio cuadrático (es decir que A=AT) tal que det(A)0. Si definimos c:=A1k, c es el centro de la curva asociada al polinomio P, C(P) donde:

(24)P(x+c)=xAx+P(c)

Como buena conclusión de este apartado, observa que las traslaciones afectan la parte lineal de los polinomios cuadráticos.

Encontrando los ejes de las rotaciones

Ahora considera la rotación g(x)=Bx con B en el general lineal de R2, es decir, BGl(2) y P el polinomio cuadrático general. Entonces:

(25)(Pg)(x)=P(Bx)=(Bx)A(Bx)+k(Bx)+f

Si desarrollamos y simplificamos esta expresión, obtenemos:

(26)(Pg)(x)=x(BTAB)x+(BTk)x+f

La pregunta en este caso es, ¿existe una forma de encontrar los ejes de la rotación a partir de esta expresión? La respuesta es sí.

A diferencia de las traslaciones, en las que se afectaba la parte lineal, para las rotaciones nos vamos a enfocar en la parte cuadrática. Debemos encontrar una manera de simplificar la expresión BTAB.

Considera a B como matriz ortogonal (BO(2)), esto implica que BTAB=B1AB que es la matriz que expresa la función A en la base de las columnas de B.

Finalmente, toma a u,v columnas de B que forman una base ortonormal y que A alarga estas columnas en factores λ,μ, es decir, que Au=λu y Av=μv. Entonces, las siguientes igualdades se cumplen:

(27)A=(λ00μ)

(28)BTAB=(u,v)TA(u,v)=(uTAuuTAvvTAuvTAv)

Si desarrollamos esta última igualdad, obtenemos:

(29)BTAB=(λ00μ)

Si encontramos una matriz B que cumpla (3), podemos eliminar el término mixto del polinomio P y acercarnos a los polinomios canónicos.

Tarea moral

  1. Demuestra el Lema 4.4.
  2. Demuestra que, para A,B,C matrices que se pueden multiplicar, se tiene que: (ABC)T=CTBTAT
  3. Encuentra el centro, si es que tienen, de las curvas asociadas a los siguientes polinomios:
    • xy3x2y2,
    • x2+2y26x+4y+3
    • 9x24xy+6y258x+24y+59

Más adelante…

Continuaremos con el estudio de la equivalencia y reducción de polinomios, con valores y vectores propios.

Cálculo Diferencial e Integral I: Rectas tangente y normal a una curva

Por Karen González Cárdenas

Introducción

En la unidad anterior vimos la teoría relacionada a las funciones derivables. A lo largo de esta última parte del curso, veremos una serie de aplicaciones de la derivada en distintos ámbitos. Esperamos que te parezcan interesantes los ejemplos que aquí expondremos y la relación del Cálculo en problemáticas de otras áreas. Comenzaremos con obtener la recta tangente y normal de una función en un punto dado.

¿Qué dice la geometría?

Recordemos algunos conceptos geométricos para entrar en contexto:
Decimos que una recta T es tangente si toca a una curva en un sólo punto. Y que una recta N es normal si es perpendicular a la recta tangente en el punto de tangencia.

  • T es la recta tangente en el punto p
  • N es la recta normal en p

En los cursos de geometría probablemente te encontraste con la siguiente ecuación para definir a una recta:
yy1=m(xx1)
ésta es conocida como la forma punto-pendiente.


Vemos que gráficamente estamos considerando un punto (x1,y1) sobre la recta y decimos que un punto cualquiera (x,y) se encuentra también sobre la recta si cumple la igualdad anterior.

Recordando…

A principios de la unidad pasada vimos que una función f es derivable en un punto x0 si existe el siguiente límite:
limxx0f(x)f(x0)xx0=f(x0).

Y que además la interpretación geométrica de dicho límite es justo la pendiente de la recta tangente a la gráfica de nuestra función f en un (x0,f(x0)).
Con ayuda de este concepto y la definición vista en la sección anterior, vemos que la recta que pasa por el punto (x0,f(x0)) y que es tangente a la gráfica sería:
yy1=m(xx1)yf(x0)=f(x0)(xx0)y=f(x0)(xx0)+f(x0)
donde m=f(x0) y consideramos (x1,y1)=(x0,f(x0)).

Definición de la recta tangente

Motivados por lo anterior tenemos la siguiente definición:
Definición (recta tangente): Sea f una función derivable en un punto x0. Definimos a la recta tangente a la gráfica de f en el punto (x0,f(x0)) como:
T(x)=f(x0)(xx0)+f(x0).

Esta definición es la que estaremos usando en todos los ejercicios de esta entrada por lo que recomendamos tenerla presente. Pasaremos ahora a definir la recta normal a la gráfica de f en el punto (x0,f(x0)).

Definición de la recta normal

Como ya vimos que geométricamente la recta normal es perpendicular a la recta tangente, modificaremos la pendiente a la definición anterior tomando m=1f(x0) con f(x0)0 :
Definición (recta normal): Tomando f una función derivable en un punto x0. Definimos a la recta normal a la gráfica de f en el punto (x0,f(x0)) con la ecuación:
N(x)=1f(x0)(xx0)+f(x0).

Con ambas rectas definidas pasaremos a resolver algunos ejercicios.

Ejemplo 1

Encuentra la recta tangente y normal de la función:
f(x)=x3+2x2x+2
en el punto (4,94).
Solución:
Comenzaremos por obtener la derivada de f(x) haciendo uso de las reglas de derivación:
f(x)=3x2+4x1.

Para obtener la pendiente en el punto indicado debemos sustituir x=4, así:
f(4)=3(4)2+4(4)1=48+161=63

Ahora comenzamos sustituyendo lo anterior en la definición de recta tangente:
T(x)=63(x4)+94=63x252+94=63x158
T(x)=63x158.

Finalmente sustituyendo en la definición de la recta normal:
N(x)=163(x4)+94=x63+463+94=x63+592663
N(x)=x63+592663.

Ejemplo 2

Encuentra la recta tangente y normal con x0=2 de la función:
f(x)=3x25x+6.
Solución:
Comenzamos por sustituir x0=2 para obtener el punto p por donde pasarán ambas rectas:
f(2)=3(2)25(2)+6=1210+6=8
p=(2,8).
Ahora pasemos a obtener la pendiente derivando la función y sustituyendo x0=2:
f(x)=6x5f(2)=6(2)5=125=7.

Procedamos a sustituir en las definiciones para la tangente y la normal:
T(x)=7(x2)+8N(x)=17(x2)+8=7x14+8=x7+27+8=7x6=x7+587
Así concluimos que:
T(x)=7x6N(x)=x7+587

Ejemplo 3

Hallar la recta tangente y normal de la función:
f(x)=x
en el punto p=(9,3).
Solución:
Procederemos a derivar la función haciendo uso de la Regla de la cadena:
f(x)=12(x)121(1)=12(x)12=12x

Obtenemos la pendiente al sustituir x0=9:
f(9)=12(9)=129=16

Ahora hallamos la recta tangente y normal sustituyendo f(9)=16:
T(x)=16(x(9))+3N(x)=116(x(9))+3=16(x+9)+3=6(x+9)+3=x632+3=6x+54+3=x6+32=6x+57

Por lo que finalmente tenemos:
T(x)=x6+32N(x)=6x+57

Más adelante

En la siguiente entrada veremos cómo encontrar máximos y mínimos de una función. Por lo tanto, definiremos dichos conceptos y probaremos algunos resultados que nos brindarán los criterios necesarios, haciendo uso de la derivada, para identificarlos.

Tarea moral

Encuentra la recta tangente y normal en cada uno de los incisos:

  • f(x)=2x3+3x2+4x2 con x0=2.
  • f(x)=x33x en p=(2,2).
  • f(x)=4x2 en p=(2,16).
  • f(x)=sen(π2x) en p=(π3,12).
  • f(x)=x+1x1 en p=(2,3).

Entradas relacionadas

Agradecimientos

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