Operaciones con números complejos

Vamos a hacer un ejercicio de números complejos muy completo para alumnos de Bachillerato y también para que alumnos de otros niveles más avanzados recuerden algunos procesos.

Calcular:

\sqrt[3]{{\left( {\frac{{1 + \sqrt 3 i}}{{1 - i}}} \right)^2 }}

Si z={\frac{{1 + \sqrt 3 i}}{{1 - i}}} =>  \frac{{1 + \sqrt 3 i}}{{1 - i}} \cdot \frac{{1 + i}}{{1 + i}} = \frac{{1 + i + \sqrt 3 i + \sqrt 3 i^2 }}{{1 - i^2 }} = \frac{{1 - \sqrt 3 }}{2} + \frac{{1 + \sqrt 3 }}{2}i

Pasamos a forma polar:

Módulo: \left| z \right| = \sqrt {\left( {\frac{{1 - \sqrt 3 }}{2}} \right)^2  + \left( {\frac{{1 + \sqrt 3 }}{2}} \right)^2 }  = 1

Argumento: \theta = arctg\frac{{\frac{{1 + \sqrt 3 }}{2}}}{{\frac{{1 - \sqrt 3 }}{2}}} = 75

De forma que tenemos \sqrt[3]{{\left( {1_{75} } \right)^2 }} = \sqrt[3]{{1^2 _{75 \cdot 2} }} = \sqrt[3]{{1_{150} }} = \left\{ \begin{array}{l}\sqrt[3]{1}_{\frac{{150 + 0 \cdot 360}}{3}}  = 1_{50}  \\ \sqrt[3]{1}_{\frac{{150 + 1 \cdot 360}}{3}}  = 1_{170}  \\ \sqrt[3]{1}_{\frac{{150 + 2 \cdot 360}}{3}} = 1_{290}  \\ \end{array} \right.

Si queremos podemos volver a pasar las soluciones a forma binómica mediante la fórmula z_\theta   = \left| z \right|\cos \theta  + \left| z \right|sen\theta i

Resuelve un problema de optimización

Iremos publicando problemas para que los lectores del Blog intenten hallar su solución (pueden publicarla en los comentarios):

Calcula las dimensiones de un cilindro de volumen máximo inscrito en un cono de revolución.

Sorry, the GeoGebra Applet could not be started. Please make sure that Java 1.4.2 (or later) is installed and active in your browser (Click here to install Java now)

Cáculo de resíduos (difíciles) 2- ejemplos.

Hola a todos. Lo prometido es deuda!! Estoy aquí, de nuevo, para daros algunos ejemplos de aplicación del método que os dí en mi primer artículo. En efecto,  veréis que, en muchos casos, el cálculo de resíduos correspondientes a polos de orden mayores a la unidad se vuelve muy engorroso mediante la fórmula (1) general que podeis encontrar en mi primer artículo. Empecemos pues.

Sea la integral impropia:

\displaystyle\int_{-\infty}^{+\infty}\displaystyle\frac{\cos{mx}}{(x^2+a^2)^2}\,dx\quad con\quad a,m>0

Para resolverla gracias al teorema de los resíduos, extendemos la función integrando al plano complejo, definiendo así:

h(z)=\displaystyle\frac{{\rm e}^{imz}}{(z^2+a^2)^2}\quad con\quad a,m\in\mathbb{R}^+

En este artículo, mi intención no es la de resolver esta integral (quien tenga dudas de cómo se hace o simplemente quiera saberlo, no dudeis en preguntar) sino que quería contextualizar el cálculo que sigue. Efectivamente, el teorema de los resíduos nos lleva a calcular el resíduo de h en el polo z_0=ia. Esto es:

\boxed{Res(h(z),z_0)}

Ante todo, es fácil ver que z_0=ia es un polo doble (de orden 2) de h. Veamos cómo sería la resolución de este problema con la fórmula general. En primer lugar, tendremos que calcular la deriada de (z-z_0)h(z)=\displaystyle\frac{{\rm e}^{imz}}{(z+ia)^2} que, de por sí, lleva algo de cálculo en el cual hay que estar atento para no equivocarse:

\displaystyle\frac{d}{dz}\{(z-z_0)h(z)\}=\displaystyle\frac{im{\rm e}^{imz}(z+ia)^2-2(z+ia){\rm e}^{imz}}{(z+ia)^4}

\displaystyle\frac{d}{dz}\{(z-z_0)h(z)\}={\rm e}^{imz}\displaystyle\frac{im(z+ia)^2-2(z+ia)}{(z+ia)^4}

\displaystyle\frac{d}{dz}\{(z-z_0)h(z)\}=\displaystyle\frac{{\rm e}^{imz}}{(z+ia)^3}(im(z+ia)-2)

\boxed{\displaystyle\frac{d}{dz}\{(z-z_0)h(z)\}=\displaystyle\frac{{\rm e}^{imz}}{(z+ia)^3}(imz-ma-2)}

Ahora, ya que en este caso 1/(n-1)!=1, falta coger el límite z\to z_0:

\displaystyle\lim_{z\to ia}\displaystyle\frac{{\rm e}^{imz}}{(z+ia)^3}(imz-ma-2)=\displaystyle\frac{{\rm e}^{-ma}}{(2ia)^3}(-ma-ma-2)

\displaystyle\lim_{z\to ia}\displaystyle\frac{{\rm e}^{imz}}{(z+ia)^3}(imz-ma-2)=\displaystyle\frac{{\rm e}^{-ma}}{-8ia^3}(-2ma-2)

\displaystyle\lim_{z\to ia}\displaystyle\frac{{\rm e}^{imz}}{(z+ia)^3}(imz-ma-2)=\displaystyle\frac{{\rm e}^{-ma}}{4ia^3}(ma+1)

\Longrightarrow\boxed{Res(h,ia)=\displaystyle\frac{{\rm e}^{-ma}}{4ia^3}(ma+1)}

Podeis ver que este cálculo no es ni muy largo ni muy difícil (luego pondré algunos más complicados) pero el método que os propuse en mi artículo es mucho más rápido y eficaz. Para entender lo que sigue, supongo que no será necesario deciros que teneis que haber entendido el planteamiento de mi post anterior. Usaré exactamente las misma notaciones. En este caso, tenemos que:

f(z)={\rm e}^{imz}

g_r(z)=(z+ia)^2

EL primer paso es el de desarrollar estas dos funciones analíticas en un disco centrado en z_0=ia en serie de potencias. Además, sabemos que, del resultado de la división según las potencias crecientes de estas dos funciones, sólo necestiremos el coeficiente del termino en \xi^{n-1}=\xi (ya que n=2 es el orden del polo), así que no necesitamos muchos términos de los desarrollos  (obviamente pondremos todos los de g_r por ser un polinomio de desarrollo finito):

f(z)={\rm e}^{imz}={\rm e}^{-ma}(1+im\xi+\cdots)

g_r(z)=(z+ia)^2=-4a^2+4ia\xi+\xi^2

Finalmente, procedemos a la división según potencias crecientes (que no es más que una división euclídea de polinomios con dichos polinomios ordenados de potencias menores a mayores) parándonos en el coeficiente del término en \xi (sólo tenemos que hacer 2 pasos de la división!!!!) obteniendo:

1+im\xi+\cdots\,:\,-4a^2+4ia\xi+\xi^2=\displaystyle\frac{1}{4a^2}-\displaystyle\frac{i}{4a^3}(ma+1)\xi

De donde, multiplicando por el factor {\rm e}^{-ma} que nos hemos dejado para hacer la división, obtenemos el resultado deseado:

\boxed{Res(h,ia)=\displaystyle\frac{{\rm e}^{-ma}}{4ia^3}(ma+1)}

Por lo tanto, hay bastante menos cálculo con éste método que con la fórmula general. Además, cuando tengáis práctica, lo que acabo de hacer se vuelve muy fácil.

Mañana os pondré más ejemplos. Que os vaya bien.

Alex.

P.D: No dudéis en pedir que exponga un ejemplo en particular.

Cálculo de resíduos (difíciles)

Hola.  Recuerdo que, cuando estudié análisis de la variable compleja, enunciamos un teorema fantástico y poderoso: El teorema de los resíduos.  Me parecía increible que el valor de una integral, a lo largo de un camino cerrado homotópico a un punto, sólo dependiera de las contribuciones de las singularidades del integrando (los resíduos) abarcadas por el camino. No obstante, así es y así se demuestra. No entraré en el detalle de la demostración de este teorema (hay mucha información en la web) ni en sus diversas aplicaciones. Lo que me interesa desarrollar a continuación es un método de cálculo de los resíduos de una función en un caso concreto de singularidades. En efecto,  dichas singularidades pueden ser aísladas o no (la singularidad de una función no es necesariamente un punto pero puede ser un segmento lo que no impide tener un resíduo), y en el caso de que sean aisladas pueden ser evitables, esenciales o ser polos. Nos centraremos, en todo lo que sigue, en las singularidades aisladas. Además, descartaremos las singularidades evitables ya que tienen resíduo nulo (el desarrollo en serie da Laurent no tiene potencias negativas). También, dejaremos de lado a las singularidades esenciales que precisan hallar explícitamente el desarrollo de Laurent para obtener el resíduo. Así pues, nos ceñiremos al estudio de los resíduos de funciones que presentan polos. Para ello, recordemos rigurosamente lo que es un polo y un resíduo:

  • Polo de orden n: Consideremos el desarrollo de Laurent de una función f en z_0:

 f(z)=\displaystyle\sum_{-\infty}^{+\infty}a_n(z-z_0)^n

z_0 es un polo de orden n \in\mathbb{N} de f si los coeficientes del desarrollo de Laurent cumplen que:

a_{-n}\neq 0   y   a_k=0   si   k<-n

Es decir que el desarrollo de Laurent tiene potencias negativas en un número finito (sino la singularidad sería esencial) y el orden del polo viene dado por el mínimo entero n tal que a_{-n}\neq 0.

  • Resíduo de f en una singularidad aislada: En el desarrollo de Laurent que hemos considerado, llamamos al coeficiente a_{-1} el resíduo de f en z_0

Ahora bien, consideremos una función compleja de la forma:

h(z)=\displaystyle\frac{f(z)}{g(z)}

con f y g analíticas. El conjunto de singularidades S de h es, por lo tanto, el conjunto de los ceros de g. Esto es:

S=\{z\in\mathbb{C}|g(z)=0\}

S es por lo tanto un conjunto de puntos aislados y son los polos de la función h. Para este tipo de singularidades, existe una fórmula con la cual podemos siempre calcular los resíduos:

\boxed{Res(h,z_0)=\displaystyle\lim_{z\to z_0}\displaystyle\frac{1}{(n-1)!}\displaystyle\frac{d^{n-1}}{dz^{n-1}} \{(z-z_0)^nh(z)\}}          (1)

donde n es el orden del polo. Veamos, pues, cual es la mejor opción para calcular resíduos según el orden del polo considerado

  • Polo simple (n=1) Cuando n=1 (polo simple), esta fórmula suele ser útil porqué es simple. En efecto, se tiene que:

Res(f,z_o)=\displaystyle\lim_{z\to z_0}(z-z_0)h(z)

Además, en nuestro caso particular  h(z)=\displaystyle\frac{f(z)}{g(z)} tenemos que:

Res(f,z_o)=\displaystyle\lim_{z\to z_0}(z-z_0)h(z)=\displaystyle\lim_{z\to z_0}\displaystyle\frac{f(z)}{\frac{g(z)}{z-z_0}}

Pero, nos fijamos en que g(z_0)=0 y \displaystyle\lim_{z\to z_0}\displaystyle\frac{g(z)-g(z_0)}{z-z_0}=g'(z_0) para obtener que:

\boxed{Res(f,z_o)=\displaystyle\frac{f(z_0)}{g'(z_0)}}

  • Polos de orden superior a uno: En este caso, la fórmula (1) no es, en general, muy útil. En efecto, en primer lugar se tendría que encontrar la derivada (no trivial en general)

\displaystyle\frac{d^{n-1}}{dz^{n-1}} \{(z-z_0)^nh(z)\}

que da una indeterminación, y luego calcular el límite (Hôpital, desarrollos limitados…). Parece pues que ésta no es la manera más inteligente de calcular el resíduo. Creo, en mi opinión, que, en la mayoría de los textos, es un error pedagógico colocar esta fórmula como “remedio milagro” para la búsqueda de resíduos correspondientes a polos. En efecto, el estudiante la memoriza (o la anota en su formulario) pensando que le servirá siempre (por ser completamente general) y se olvida entender realmente  la esencia de lo que esta haciendo. Luego, si en el exámen le sale un polo de orden 4, por ejemplo, se encuentra con unos cálculos tediosos para los cuales no está preparado y, por ende, pierde un tiempo valioso con muchas probabilidades de equivocarse o rendirse delante la complejidad. Por ello, la esencia de mi artículo toma todo su sentido en este preciso instante. Intentaré proporcionaros un método fiable, rápido e intutivo de como encontrar los resíduos de una función en polos de orden mayores a la unidad.

La verdad es que, a los que teneis que usar el teorema de los resíduos, os recomiendo encarecidamente que cojais lapiz y papel y os apuntéis lo que sigue (no es complicado pero no basta con leerlo). No os preocupéis si, al principio, os cuesta un poco entenderlo. No deja de ser un método así que requiere un poco de práctica. Además, el tiempo que invertáis  en practicar estará, de sobra, compensado por la rápidez con la cual seréis capaces de hallar los resíduos de un función que presenta polos (que, en general, son la mayoría en las clases de variable compleja). Ánimo y a estudiar jejeje.

Nuestro propósito es el de calcular Res(h,z_0) siendo z_0 un polo de orden n de h. Entonces, esto quiere decir que z_0 es un cero de orden (multiplicidad) n de g. Al ser, g analítica, le corresponde un desarrollo en serie como sigue:

g(z)=\displaystyle\sum_{k=0}^{\infty}a_k(z-z_0)^k

Pero todos los términos hasta el rango n son nulos, por lo que podemos escribir que:

\boxed{g(z)=\displaystyle\sum_{k=n}^{\infty}a_k(z-z_0)^k=(z-z_0)^n\underbrace{\{a_n+a_{n+1}(z-z_0)+a_{n+2}(z-z_0)^2+\cdots\}}_{g_r(z)}}

con a_n\neq 0 (sino z_0 sería un cero de orden n+1 de g) y donde hemos introducido g_r que llamaremos función reducida de g. Esta función reducida de g es analítica en un disco centrado en z_0 y no se anula en dicho punto, por lo que podemos definir la función reducida de h como:

h_r(z)=\displaystyle\frac{f(z)}{g_r(z)}

Está claro, pues, que h_r es analítica en el disco de analiticidad de g_r centrado en z_0.  ¿Porqué hemos intoducido estas funciones reducidas? Pues porqué, ahora, podemos hallar la serie de Laurent de h obteniendo la serie entera (de Taylor con potencias positivas exclusivamente) de h_r, y,como consecuencia, podremos hallar el resíduo que buscamos, que corresponde al coeficiente del término en (z-z_0)^{-1} de la serie de Laurent de h, a partir del desarrollo en serie de h_r. Y ¿a qué coeficientes de dicha serie corresponde el resíduo? Pues al término en (z-z_0)^{n-1} del desarrollo de  h_r. En efecto, tenemos que:

\boxed{h(z)=\displaystyle\frac{f(z)}{(z-z_0)^n g_r (z)}=\displaystyle\frac{h_r(z)}{(z-z_0)^n}}

Finalmente, bastará con hallar el desarrollo en serie entera de h_r hasta el orden n-1 para hallar Res(h,z_0).

Ahora, nos podemos fijar en que, a partir de este resultado, volvemos a encontrar la fórmula (1). En efecto, el coeficiente del desarrollo de Taylor de h_r correspondiente a la potencia n-1-ésima es, por definición:

\displaystyle\frac{1}{(n-1)!}\displaystyle\frac{d^{n-1}}{dz^{n-1}}\{h_r(z_0)\}

Pero, h_r(z)=(z-z_0)^nh(z). Entonces, este coeficiente (que es el resíduo) se puede encontrar mediante el límite de la fórmula (1). No obstante, a la práctica, será mucho más rápido y eficiente hallar el desarrollo en serie de Taylor de h_r en z_0 mediante la división de potencias crecientes de la serie de f respecto de g_r. Esto es, tenemos los siguientes desarrollos en serie:

f(z)=b_0+b_1\xi+b_2\xi^2+\cdots\\g_r(z)=a_n+a_{n+1}\xi+a_{n+2}\xi^2+\cdots

donde los a_k son los coeficientes del desarrollo de g y donde hemos hecho \xi=z-z_0. Planteamos la división en potencias crecientes:

b_0+b_1\xi+b_2\xi^2+\cdots\quad |\underline{ a_n+a_{n+1}\xi+a_{n+2}\xi^2+\cdots}\Longrightarrow\displaystyle\frac{b_0}{a_n}+\cdots

Cuando lleguemos a la potencia n-1, tendremos el resíduo deseado. El proceso de la división en potencias crecientes parece ser engorroso, no obstante, a la práctica, lo es muchísimo menos que el cálculo de la derivada y del límite planteado en la fórmula (1) .

Espero que os sirva. En los días próximos, incluiré en este post algunos ejemplos de la potencia de este método frente a la fórmula general.

Que os vaya muy bien!!

Alex.

Ecuaciones diofánticas

Una ecuación diofántica es una ecuación algebraica con soluciones en los números enteros. Su nombre lo deben al matemático Diofanto de Alejandría.

Comenzamos estudiando la ecuación diofántica lineal de dos incógnitas:

Ax + By = C\,\,con\,\,A,\,\,B\,\,y\,\,C \in {\rm Z}

Esta ecuación tiene solución si y solo si mcd(A,B) divide a C. En este caso la ecuación tiene infinitas soluciones:

\left\{ \begin{gathered}x = x_0  + \lambda \frac{B}{d} \hfill \\y = y_0  - \lambda \frac{A}{d} \hfill \\ \end{gathered}  \right. donde \left( {x_0 ,y_0 } \right) es una solución particular de la ecuación, d=mcd(A,B) y \lambda un parámetro entero.

Para hallar la solción particular usaremos la Identidad de Bezout junto con el Algortimo de Euclides.

Veamos un ejemplo:

30x + 12y = 1200

mcd(30,12)=6/1200 luego la ecuación diofántica tiene solución 6 = 30 - 12 \cdot 2 \Rightarrow 200 \cdot 6 = 200 \cdot 30 - 400 \cdot 12

Por tanto la solución particular es (200,-400) y la solución general será: \left( {200 + 2t, - 400 - 5t} \right)

Problema clásico del mono y los cocos:

Cinco hombres y un mono naufragan en una isla desierta. Los hombres pasan todo el primer día recogiendo cocos. Por la noche, uno de ellos despierta y, desconfiado, decide separar su parte. Dividió los cocos en cinco montones, y como sobraba un coco, se lo dio al mono. Poco después un segundo náufrago se despierta y hace lo mismo. Al dividir los cocos en cinco montones volvió a sobrar un coco y también se lo dio al mono. Uno tras otro, el tercero, cuarto y quinto náufragos hacen lo mismo. Por la mañana, al día siguiente, dividieron los cocos en cinco montones sin que sobrara ninguno. ¿Cuántos se habían recolectado inicialmente?

\begin{gathered}x=5x_1 +1\hfill\\x-x_1-1=4x_1=5x_2+1\hfill\\x-x_1 - x_2 - 2 = 4x_2 = 5x_3 + 1 \hfill \\x - x_1 - x_2 - x_3 - 3 = 4x_3 = 5x_4 + 1 \hfill \\x - x_1 - x_2 - x_3 - x_4 - 4 = 4x_4 = 5x_5 + 1 \hfill \\x - x_1 - x_2 - x_3 - x_4 - x_5 - 5 = 4x_5 = 5y \hfill \\\end{gathered}

Sustituyendo empezando por la última y hacia la primera:

4 \cdot 4x_4 = 4 \cdot 5x_5 + 4 \cdot 1 \Rightarrow 16x_4 = 5 \cdot 4x_5 + 4 = 5 \cdot 5y + 4 64x_3 = 5\left( {25y + 4} \right) + 16 = 125y + 36 256x_2 = 5\left( {125y + 36} \right) + 64 = 625y + 244
1024x_1 = 5\left( {625y + 244} \right) + 256 = 3125y + 1476
1024x = 5\left( {3125y + 1476} \right) + 1024 = 15625y + 8404

Obtenemos finalmente la ecuación:

1024x - 15625y = 8404

Nota: continuaremos el post resolviendo esta última ecuación diofántica, haciendo algunas demostraciones de lo visto hasta ahora, y por último nos enfrentaremos a la famosa ecuación pitagórica.