venres, 29 de marzo de 2024

Fraccións continuas teito e as montañas de Galicia

por Andrés Ventas

Fraccións continuas teito

Unha fracción continua teito ($fct$) é unha fracción continua na que usamos a función teito para obter os seus coeficientes. Os seus converxentes $\frac{p_i}{q_i}$ obtéñense mediante unha recorrencia de resta,

$p_{-1}=1, \ p_{0}=c_0, \ p_i=c_i p_{i-1} - p_{i-2}$.

$q_{-1}=0, \ q_{0}=1, \ q_i=c_i q_{i-1} - q_{i-2}$.

Para o cálculo da fracción continua, se $x$ é racional usamos o algoritmo de Euclides coa función teito e se o número e irracional usamos unha iteración sobre o inverso de $x$ restando en cada paso o valor teito do resultado.

Hai varias notacións.É moi cómoda a reducida, choendo os coeficientes entre símbolos da función teito $\lceil c_0, c_1, c_2, c_3, \cdots \rceil$.

A notación tradicional é \[ \mathrm{x}= c_0-\cfrac{1}{c_1-\cfrac{1}{c_2-\cfrac{1}{c_3- \cdots\vphantom{\cfrac{1}{1}} }}} \]

E a maiores hai outra notación semireducida $c_0 - \dfrac{1}{c_1 \ -} \ \dfrac{1}{c_2 \ - } \ \dfrac{1}{c_3 \ - } \ \dfrac{}{\ldots} $,

Vexamos dous exemplos de cálculo:

Con Euclides teito, calculamos a $fct$ de $\dfrac{93}{20}$ e os seus converxentes,

Cociente teitoresto
9320$c_0=5$-7
207$c_1=3$-1
71$c_2=7$0
537
$p_i$ 1 51493
$q_i$ 0 1320

Así $\dfrac{93}{20}=\lceil 5, 3, 7 \rceil = 5-\cfrac{1}{3-\cfrac{1}{7\vphantom{\cfrac{1}{1}} }}$. (Como fracción continua regular $\dfrac{93}{20}=[4, 1, 1, 1, 6 ]).$

Agora con inverso iterativo, calculamos $\sqrt{3} \approx 1.7320$ e os seus converxentes,

$c_i = \lceil x_{i-1} \rceil$ $r_i=c_i - x_{i-1}$$x_i = 1/r_i $
$x_{-1} \approx 1.7320$
$c_0=2$0.26803.7320
$c_0=4$0.26803.7320
$c_0=4$0.26803.7320
$\ldots$$\ldots$ $\ldots$
2444$\cdots$
$p_i$ 1 272697$\cdots$
$q_i$ 0 141556$\cdots$

Así temos $\sqrt{3}=\lceil 2, 4, 4, \cdots \rceil$,

sendo $\sqrt{3} \approx 1.7320$ e o cuarto converxente $\dfrac{97}{56} \approx 1.7321$.

Comparten moitas propiedades coas fraccións continuas regulares, mais evidentemente teñen as súas pecularidades. No libro de Sergey Khrushchev, Orthogonal Polynomials and Continued Fractions temos un minucioso tratado.

As $fct$ converxen sempre polo lado superior.

Nalgúns casos as $fct$ dan o menor número de pasos, aínda que normalmente son peores en número de pasos. De feito non cumpren coa constante de Khinchin, porque se atoan nos coeficientes, fundamentalmente no $2$.

No 1657 Brouncker atopa, seguindo outro camiño, unha fracción continua teito para solucionar a ecuación de Pell $x^2 - 2y^2 =1$, onde as solucións son certos converxentes da $fct=\lceil 6, 6, 6, \ldots \rceil$, con recorrencia $x_n=6x_{n-1}-x_{n-2}$ e por tanto a solución consiste en resolver a ecuación $x=6-\dfrac{1}{x}$, que dá $x=3+2\sqrt{2}$. (Lema 2.21 do libro de Khrushchev).

Se se permiten coeficientes en $\mathbb{Q}$ ou mesmo en $\mathbb{C}$ os resultados son máis interesantes aínda, serían funcións teito xeneralizadas que van en paralelo coas fraccións continuas xeneralizadas.

Fraccións continuas mixtas

Na fracción continua mixta ($fcm$) usamos a función teito ou chan en cada paso en función do menor residuo. A notación sería igual que a da fracción continua regular mais levando un superindice "$\mathbf{-}$" sobre os coefcientes que se obteñan coa función teito $[ c_0, c_1^{\mathbf{-}}, c_2^{\mathbf{-}}, c_3, \cdots ]$.

Para obter os seus converxentes $\dfrac{p_i}{q_i}$ temos recorrencia de suma ou resta en función do superíndice do coeficiente $c_{i-1}$, restamos cando $c_{i-1}^{\mathbf{-}}$ ten superíndice negativo.

$p_{-1}=1, p_{0}=c_0, p_i=c_i p_{i-1} \pm p_{i-2}$.

$q_{-1}=0, q_{0}=1, q_i=c_i q_{i-1} \pm q_{i-2}$.

(Signo negativo se $c_{i-1}$ ten superíndice $\mathbf{-}$).

Este tipo de fracción dá o menor número de pasos no algoritmo de Euclides e igualmente dá o menor número de termos na fracción continua. Outra vez vexamos dous exemplos.

Con Euclides mixto, calculamos a $fcm$ de $\frac{2114}{61}$, e os seus converxentes,

Cociente teito ou chanresto
211461$c_0=35^{\mathbf{-}}$-21
6121$c_1=3^{\mathbf{-}}$-2
212$c_2=10$ 1
21$c_3=2$ 0
$35^{\mathbf{-}}$$3^{\mathbf{-}}$ 102
$p_i$ 1 3510410052114
$q_i$ 0 132961

Por tanto $\dfrac{2114}{61}=[ 35^{\mathbf{-}}, 3^{\mathbf{-}}, 10, 2]$.

Xa temos unha moi boa cousa positiva das $fcm$ pois obtemos o $mcd$ nun menor número de pasos e por tanto tamén o multiplicativo modular inverso.

Constantes das montañas de Galicia

Chámanse medias ou constantes metálicas os valores das fraccións continuas que teñen todos os seus coeficientes iguais, a máis famosa precisamente a proporción aurea $\varphi=\frac{1+\sqrt{5}}{2}=1.618033\ldots$ con fracción continua regular $[1,1,1,1, \ldots]$.

Así que imos aproveitar para definir, usando as $fct$ de forma similar, as constantes das montañas de Galicia. Temos que $fct=\lceil 3, 3, 3, \ldots \rceil = \frac{3+\sqrt{5}}{2}=2.618033\ldots$ que como vemos é $\varphi + 1$ e a esta constante ímoslle chamar Pena Trevinca pois o nome évos ben acaído para ese conxunto de treses.

É fácil demostrar que o valor para calquera $fct$ cos coeficientes repetidos é $$\lceil c, c, c, c, \ldots \rceil=\dfrac{c+\sqrt{c^2-4}}{2}$$ por tanto podemos asignar unha constante a cadanseu monte:

O caso da $fct=\lceil 2, 2, 2, \ldots \rceil$ é particular pois ten como solución unha constante enteira $1$, así que imos deixar esa constante para o Monte da Guía en Vigo que ten só $100$ metros.

$fct$ ConstanteMonteAltura
$ \lceil 2, 2, 2, \ldots \rceil$ $\frac{2+\sqrt{0}}{2}=1$ Monte da Guía (Vigo) 100 metros
$ \lceil 3, 3, 3, \ldots \rceil$ $\frac{3+\sqrt{5}}{2}=2.618033$ Pena Trevinca (Serra do Eixo) 2127 metros
$ \lceil 4, 4, 4, \ldots \rceil$ $\frac{4+\sqrt{12}}{2}=2+\sqrt{3}=3.732050$ Cuíña (Serra dos Ancares) 1987 metros
$ \lceil 5, 5, 5, \ldots \rceil$ $\frac{5+\sqrt{21}}{2}=4.791287$ Manzaneda (Serra da Queixa) 1781 metros
$ \lceil 6, 6, 6, \ldots \rceil$ (Brouncker) $\frac{6+\sqrt{32}}{2}=3+2\sqrt{2}=5.828427$ Formigueiros (Serra do Courel) 1639 metros
$ \lceil 7, 7, 7, \ldots \rceil$ $\frac{7+\sqrt{45}}{2}=6.854101$ O Turrieiro (Serra da Enciña da Lastra) 1612 metros
$ \lceil 8, 8, 8, \ldots \rceil$ $\frac{8+\sqrt{60}}{2}=4+\sqrt{15}=7.872983$ Monte Faro (Serra do Faro) 1187 metros
$ \lceil 9, 9, 9, \ldots \rceil$ $\frac{9+\sqrt{77}}{2}=8.887482$ O Cadramón (Serra do Xistral) 1056 metros

Bibliografia

  • Sergey Khrushchev, Orthogonal Polynomials and Continued Fractions
  • Wikipedia, Fracción continua
  • Wikipedia Fraccións continuas xeneralizadas.
  • Wikipedia, Constante de Khinchin
  • Wiwipedia, Metallic Means
  • luns, 11 de marzo de 2024

    Un regalo da xeometría inversiva

    Nunca estudei a xeometría inversiva. O máis cerca que estiven diso foi nas clases de Álxebra II, no segundo cuso da carreira, cando traballamos a razón dobre. Aquel achegamento , facendo honra á denominación da materia, foi puramente alxébrico. Nesas clases nunca debuxamos unha circunferencia. Así, cando vexo algúns apuntes sobre ese tema ando sempre ás atoutiñadas, todo me sorprende. 

    Nunha entrada anterior xa se explicaba en que consistía a inversión dun plano mediante unha circunferencia. Recórdoo de seguido. Trátase dunha transformación do plano (agás un punto) en si mesmo. Dada unha circunferencia, (que chamaremos circunferencia inversiva) de centro O (centro de inversión) e raio R construiremos a inversión así:


    Se $P$ é un punto do círculo de centro $O$ e raio $R=OT$ trazamos a semirecta $OP$ e a súa perpendicular polo punto $P$. Esta perpendicular cortará en dous puntos á circunferencia. A tanxente nun destes puntos cortará a semirecta $OP$ nun punto $P'$ que será a inversión de $P$

    No caso de que $P$ fique fóra do círculo a obtención de $P'$ sería semellante. Desde $P$ trazamos unha das tanxentes á circunferencia $PT$. Despois trázase a perpendicular a $OP$ por $T$ e obtemos $P'$

    En calquera caso  o inverso dun punto $P$ respecto dunha circunferencia de centro $O$ e raio $R$ é outro punto $P'$ na semirecta $OP$ tal que $OP\cdot OP'=R^{2}$ Son bastante evidentes as seguintes propiedades:

    • O inverso do inverso é o propio punto
    • O inverso dun punto interior á circunferencia fica fóra da mesma e viceversa.
    • Os únicos puntos auto-inversos son os da circunferencia

    Na entrada á que facía referencia anteriormente deducíranse algunhas propiedades da inversión usando un método un tanto estrafalario (mediante o uso da proxección estereográfica). Alí establecimos os seguintes resultados:

    • A inversión conserva os ángulos
    • A inversa dunha recta que pasa por O é a propia recta (sempre que omitamos o propio punto O, que é o único que non ten imaxe)
    • A inversión transforma circunferencias que non pasan por O en circunferencias
    • A inversión transforma circunferencias que pasan por O en rectas.

    Imos demostrar a última propiedade mediante un procedemento máis estándar. Recóllo do libro "Regreso a la geometría" de H.S.M Coxeter e S.L. Greitzer (La tortuga de Aquiles 1993).

    Consideremos unha recta $a$ que non pase por $O$. Tracemos desde $O$ a perpendicular a $a$, $=OA$. Sexa $A'$ a inversa de $A$. Debuxemos a circunferencia $\alpha $ de diámetro $OA$. Desde un punto $P\in a$ trazamos o segmento $OP$ que cortará a $\alpha $ nun punto $P'$. 

    Os triángulos $\triangle OP'A'$ e $\triangle OPA$ son semellantes xa que comparte o ángulo en O e ambos teñen ademais un ángulo recto. De aí que tamén $$\frac{OP}{OA}=\frac{OA'}{OP'}\Rightarrow OP\cdot OP'=OA\cdot OA'=R^{2}$$

    Entón $P'$ é o inverso de P.  Recíprocamente calquera punto da circunferencia $\alpha$ invértese noutro da recta $a$


    A lonxitude dun segmento invertido

    Seguindo o libro de Coxeter e Greitzer abordaremos agora un teorema que explica como se modifica a distancia mediante a inversión. 

    Fórmula da lonxitude dun segmento invertido. Se unha circunferencia $\omega$ de centro $O$ e raio $R$ inverte os puntos $A$ e $B$ en $A'$ e $B'$, as distancias verifican a seguinte relación $$A'B'=\frac{R^{2}\cdot AB}{OA\cdot OB}$$


    Sexa $A'$ inverso de $A$: $OA\cdot OA'=R^{2}$

    Sexa $B'$ inverso de $B$: $OB\cdot OB'=R^{2}$

    Entón $OA\cdot OA'=OB\cdot OB'$ polo que e  $$\frac{OA}{OB}=\frac{OB'}{OA'}$$

    Daquela os triángulos $\triangle OAB$ e $\triangle OA'B'$ son semellantes pois o ángulo en O é común e os lados que o determinan son proporcionais (aplicamos o chamado criterio LAL). Polo tanto temos tamén que $$\frac{A'B'}{AB}=\frac{OA'}{OB}=\frac{OA\cdot OA'}{OA\cdot OB}=\frac{R^{2}}{OA\cdot OB}$$

    Despexando, obtemos a igualdade prometida: $$A'B'=\frac{R^{2}\cdot AB}{OA\cdot OB}$$

    Un regalo: o teorema de Ptolomeo

    Hai moitas demostracións do teorema de Ptolomeo. Algunhas teñen base trigonométrica, outras usan a semellanza de triángulos. No libro que vimos mencionando demóstrase este teorema usando a recta de Simson-Wallace (para este tema ver estas outras entradas neste mesmo blogue [1] e [2]). Resulta que tamén hai unha demostración inversiva, ademais é directa e simple.

    Antes de nada lembremos que un cuadrilátero cíclico é aquel que ten os catro vértices nunha mesma circunferencia. Sabendo isto podemos enunciar o

    Teorema de Ptolomeo. Se ABCD é un cuadrilátero cíclico, a suma dos produtos dos pares de lados opostos é igual ao produto das diagonais, isto é: $AC\cdot BD=AD\cdot BC+AB\cdot BC$

    Sexa $\alpha$ a circunferencia pola que pasan os catro vértices do cuadrilátero. Tomando o vértice $A$ como centro, construímos outra circunferencia $\omega$ que conteña ao cuadrilátero. Agora, como os outros vértices, $B$, $C$ e $D$ están en $\alpha$, unha circunferencia que pasa polo centro de inversión, se os invertemos, as súas imaxes $B'$, $C'$ e $D'$ ficarán todas nunha recta. Velaí que $$B'D'=B'C'+C'D'$$

    Fagamos agora uso da fórmula, dada anteriormente,  que nos dá a distancia dun segmento invertido:$$\frac{R^{2}\cdot BD}{AB\cdot AD}=\frac{R^{2}\cdot BC}{AB\cdot AC}+\frac{R^{2}\cdot CD}{AC\cdot AD}$$

    Eliminando $R^{2}$ e sacando denominadores facendo uso de que o seu mínimo común múltiplo é $AB\cdot AC\cdot AD$:

    $$AB\cdot BD=AD\cdot BC+AB\cdot CD$$

    O recíproco do teorema de Ptolomeo tamén se verifica, isto é:

    Recíproco do teorema de Ptolomeo. Se ABCD é un cuadrilátero tal que: $AC\cdot BD=AD\cdot BC+AB\cdot BC$, entón ABCD é un cuadrilátero cíclico.
    Por outra banda, dados catro puntos calquera, sempre se verificaría a desigualdade, coñecida tamén como desigualdade de Ptolomeo: $$AB\cdot BD\leq AD\cdot BC+AB\cdot CD$$
    Ademais a igualdade que nos ofrece o teorema de Ptolomeo caracteriza os cuadriláteros cíclicos.


    Como colofón citaremos a seguinte aplicación do teorema de Ptolomeo. Consideremos o pentágono regular de lado unidade. As súas diagonais son evidentemente iguais (son a base de triángulos isósceles de lado 1); poñamos que miden $\varphi$ Desbotando por exemplo o vértice superior quédanos un cuadrilátero cícliclo.

    Se lle aplicamos a este cuadrilátero o teorema de Ptolomeo teremos que $$\varphi \cdot \varphi =1\cdot \varphi +1$$

    $$\varphi ^{2}=\varphi +1$$

    Esta é unha ecuación doada de resolver, ou de recoñecer.