11 votos

Problema de producto infinito mediante la iteración de una función: $ \exp(x) = x \cdot f^{\circ 1}(x)\cdot f^{\circ 2}(x) \cdot \ldots $

[actualización]: He hecho la pregunta más precisa, más general y he añadido una pregunta de seguimiento

Considerando la iteración de funciones (con enfoque en la exponenciación iterada) estoy buscando, si alguna función que quiero iterar puede -espero que con alguna ventaja- ser expresada por iteraciones de una función -por así decirlo- "más básica".

Fo ejemplo asumo una función $f(x)$ tal que

$ \qquad \displaystyle \exp(x) = x \cdot f^{\circ 1}(x)\cdot f^{\circ 2}(x)\cdot f^{\circ 3}(x)\cdots $
$\qquad \qquad$ (donde la anotación del círculo significa iteración, y $f^{\circ 0}=x, f^{\circ 1}(x)=f(x)$ )

y primero pregunto: ¿qué hace esta función $f(x)$ ¿se ve así?

Lo que hago es esta sustitución: $$ \small \begin{array} {lrll} 1.& \exp(x) & = &x & \cdot f^{\circ 1}(x) & \cdot f^{\circ 2}(x) & \cdot f^{\circ 3}(x) & \cdots \\ 2.& \exp(f(x))&= && f^{\circ 1}(x) & \cdot f^{\circ 2}(x) & \cdot f^{\circ 3}(x) & \cdots \\ \\ \\ 3.& {\exp(f(x))\over \exp(x) } & = & \frac 1x \\ \\ & \exp(f(x)) & = & &{ \exp(x) \over x} \\ \\ \\ 4. & f(x)&=& x & - \log(x) \end{array} $$ $ \qquad \qquad $ (De 1. y 4. Sé, que x está ahora restringido a $x \gt 0$ )

Como señalan algunos comentarios, la construcción de la función $f(x)$ está infradeterminado, por lo que en el paso $3.$ el numerador y el denominador pueden tener un factor común $c$ de tal manera que tendremos
$$ \small \begin{array} {lrll} 1a.& \exp(x)\cdot c & = &x & \cdot f^{\circ 1}(x) & \cdot f^{\circ 2}(x) & \cdot f^{\circ 3}(x) & \cdots \\ 2a.& \exp(f(x)) \cdot c&= && f^{\circ 1}(x) & \cdot f^{\circ 2}(x) & \cdot f^{\circ 3}(x) & \cdots \end{array} $$

Si hago ahora el cálculo con algún ejemplo $x$ por $$ y = x \cdot f^{\circ 1}(x)\cdot f^{\circ 2}(x)\cdot f^{\circ 3}(x)\cdots $$ Recibo por todo lo probado $x>0$ el resultado $$ y = \exp(x) / \exp(1) $$ de manera que, efectivamente, un cofactor $c$ se produce y que es precisamente $1 / \exp(1)$

Q1: ¿De dónde procede este factor adicional en la evaluación empírica? ¿Dónde se ha perdido alguna información crucial en los pasos analíticos anteriores?

La pregunta puede ser más precisa:
Q2: ¿Cómo determina el cálculo empírico, que el cofactor $c$ es sólo $1/\exp(1)$ ?

(Los comentarios de @Eric Wong abordan esta cuestión, pero aún no la he explicitado)


Al revisar mi propia pregunta, es posible una generalización en el sentido de que puedo utilizar cualquier base $b$ con $\log(b)\ge 1$ tal que

$ \qquad \displaystyle {b^x \over b} = x \cdot f_b^{\circ 1}(x)\cdot f_b^{\circ 2}(x)\cdot f_b^{\circ 3}(x)\cdots $

y la constante $c$ sale siempre igual a $1/b$ de tal manera que también podríamos escribir
$ \qquad \displaystyle b^{x-1} = x \cdot f_b^{\circ 1}(x)\cdot f_b^{\circ 2}(x)\cdot f_b^{\circ 3}(x)\cdots $

Q3: ¿El comentario de Eric puede ser más explícito para que sirva para todas las bases?


Observo que para $\eta < \log(b) <1$ con $\eta \lt 0.39996 $ la secuencia de iterados de las funciones $f_b(x)$ cualquiera de los dos enfoques $1$ alternando desde abajo y desde arriba o no se acercan $1$ en absoluto, sino que se acercan a distintos puntos de acumulación... Esta observación es importante aquí, porque para tales bases $b$ la fórmula del producto anterior no funciona correctamente, porque a veces no se consigue la convergencia a un único punto fijo. Pero debido al comentario de @did he movido esa pregunta a un hilo conductor separado

Un fragmento de código utilizando Pari/GP:

f(x) = x-log(x)  \\ define the function 

x0=1.5
     \\  = 1.50000000000
[tmp=x0,pr=1]              \\ initialize
for(k=1,64,pr *= tmp;tmp = f(tmp));   pr   \\ compute 64 terms, show result
      \\ = 1.64872127070

exp(x0)        \\ show expected value
       \\ = 4.48168907034

pr*exp(1)      \\ show, how it matches
       \\  = 4.48168907034

Este es un ejemplo que muestra el tipo de convergencia; utilizo $x_0=1.5$ y una precisión interna de 200 dígitos decimales. Entonces obtenemos los términos del producto parcial como $$ \small \begin{array} {r|r} x_k=f^{\circ k}(x) & (x_k-1) \\ \hline 1.50000000000 & 0.500000000000 \\ 1.09453489189 & 0.0945348918918 \\ 1.00420537512 & 0.00420537512103 \\ 1.00000881788 & 0.00000881787694501 \\ 1.00000000004 & 3.88772483656E-11 \\ 1.00000000000 & 7.55720220223E-22 \\ 1.00000000000 & 2.85556525627E-43 \\ 1.00000000000 & 4.07712646640E-86 \\ 1.00000000000 & 8.31148011150E-172 \\ 1.00000000000 & 1.020640763E-202 \\ 1.00000000000 & 1.020640763E-202 \\ \cdots & \cdots \end{array} $$

4voto

Did Puntos 1

Digamos que $f$ resuelve $(\ast)$ para $g$ si $g(x)=x\cdot f(x)\cdot f\circ f(x)\cdot f\circ f\circ f(x)\cdots$ en el sentido de que el producto en el lado derecho converge y que su valor es $g(x)$ para cada positivo $x$ .

Los cálculos en el post y en algunos comentarios muestran que, si $f$ resuelve $(\ast)$ para $g:x\mapsto c\mathrm e^x$ para algún positivo $c$ entonces $f:x\mapsto x-\log x$ y que $f:x\mapsto x-\log x$ efectivamente resuelve $(\ast)$ para $g:x\mapsto\mathrm e^{x-1}$ y no para cualquier otro $g:x\mapsto c\mathrm e^x$ .

En general, hay que tener en cuenta que si $f$ resuelve $(\ast)$ para $g$ entonces $g(f(x))=f(x)\cdot f\circ f(x)\cdot f\circ f\circ f(x)\cdots$ por lo que $$ x\cdot g(f(x))=g(x). $$ Suponiendo que $g$ es invertible, por ejemplo porque $g$ es creciente, consideramos la función $Tg$ definido por $$ Tg(x)=g^{-1}(g(x)/x). $$ Entonces la única solución posible de $(\ast)$ es $f=Tg$ En particular $f$ es único y tal que $f(1)=1$ cuando existe. Así, para $(\ast)$ para tener soluciones, hay que suponer que $g(1)=1$ .

Supongamos además que $g'(1)$ existe, por lo tanto, $g(1+\varepsilon)=1+g'(1)\varepsilon+o(\varepsilon)$ cuando $\varepsilon\to0$ . Entonces $Tg(1+\varepsilon)=1+\eta$ con $g'(1)\eta+o(\eta)=(g'(1)-1)\varepsilon+o(\varepsilon)$ . Si $g'(1)\lt1/2$ el punto $1$ es repulsivo para $f$ por lo que el producto en el lado derecho de $(\ast)$ diverge y $(\ast)$ no tiene solución. Si $g'(1)\gt1/2$ la relación $(g'(1)-1)/g'(1)$ está en $(-1,1)$ por lo que el producto infinito converge y $(\ast)$ tiene la solución única $Tg$ .

En el caso concreto de que $g(x)=\mathrm e^{x-1}$ entonces $g'(1)=1$ por lo que $(\ast)$ tiene una solución única $f=Tg$ y uno sabe que $Tg:x\mapsto x-\log x$ .

En el caso concreto de que $g(x)=b^{x-1}$ para algún positivo $b$ entonces $g'(1)=\log b$ por lo que $(\ast)$ tiene una solución única $f=Tg$ cuando $b\gt\sqrt{\mathrm e}$ y uno sabe que $Tg:x\mapsto x-\log_b x$ y $(\ast)$ no tiene solución cuando $b\lt\sqrt{\mathrm e}$ .

En el régimen límite $b=\sqrt{e}$ un análisis más cuidadoso de las secuencias $(x_n)$ definido por $1+x_{n+1}=Tg(1+x_n)$ muestra que $(x_n)$ es finalmente alternativo (decreciente en módulo y alternativamente positivo y negativo para cada $n$ suficientemente grande), por lo que el producto infinito $\prod\limits_n(1+x_n)$ converge y $(\ast)$ tiene una solución única $f=Tg:x\mapsto x-2\log x$ también en este caso.

1voto

Fernando Puntos 11

La constante c de (1a) se anula en el paso (3). La constante c debe calcularse después de obtener f(x). Una forma es numéricamente. Otra forma es algebraicamente, pero aún no sé cómo.

0voto

Jorrit Reedijk Puntos 129

No es una respuesta completa, pero posiblemente sea un comienzo... Podemos observar que para $x=0$ la lhs es distinta de cero, por lo que en la rhs el producto infinito es $0 \cdot f(0) \cdot \ldots$ y el $f^{\circ k}(0)$ debe tener en cuenta ese cofactor cero. Ahora bien, si tenemos $f(x) = x - \log(x)$ entonces esto es efectivamente infinito en cero - pero esta opción podría ser la/una razón del problema.

Si reformulamos la función f(x) de forma que $$ \exp(x) = (x+1) \cdot f^{ \circ 1}(x+1) \cdot \ldots $$ entonces en x=0 tenemos en el lado derecho un producto de todos los unos, que es lo que "nos gustaría" y podríamos reescribir $$ \exp(x-1) = {\exp(x) \over e } = x \cdot f^{ \circ 1}(x) \cdot \ldots $$ y obtener la relación, que en realidad encontramos empíricamente.
Sin embargo, creo que ese argumento debe hacerse más fuerte y formal; en lo que respecta a mis pensamientos, es simplemente una insinuación.
[actualización]: Además, cambiando la función de esta manera, tal que $g(x) = (x+1) + \log(1+x)$ se estrella con el producto infinito, se vuelve muy divergente...

Además, esto me recuerda a un efecto de aspecto similar, a saber, la función gamma incompleta, en la que la relación de los gammas incompletos en x y en x+1 se comportan como la propia gamma, pero tienen un término constante adicional -por ejemplo $\frac 1e$ -, que también se expresa mediante alguna integral. Sin embargo - esto es sólo una vaga sensación de estar relacionado de alguna manera, no puedo realmente formalizar ...

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X