Gráfica de la función Omega de Wright en el eje real.
En matemáticas , la función Omega de Wright , también llamada función de Wright , es una función que está definida por la función W de Lambert :
ω
(
z
)
=
W
⌈
I
m
(
z
)
−
π
2
π
⌉
(
e
z
)
{\displaystyle \omega (z)=W_{\left\lceil {\frac {\mathrm {Im} (z)-\pi }{2\pi }}\right\rceil }(e^{z})}
Donde
I
m
(
z
)
{\displaystyle \mathrm {Im} (z)}
es la parte imagniaria de z,
⌈
z
⌉
{\displaystyle \lceil z\rceil }
es la función techo de z y
K
(
z
)
=
⌈
I
m
(
z
)
−
π
2
π
⌉
{\displaystyle \mathrm {K} (z)=\left\lceil {\frac {\mathrm {Im} (z)-\pi }{2\pi }}\right\rceil }
es un "número desenrollado" de z. Anteriormente
K
(
z
)
{\displaystyle \mathrm {K} (z)}
estaba definida como
K
(
z
)
=
⌊
π
−
I
m
(
z
)
2
π
⌋
{\displaystyle \mathrm {K} (z)=\left\lfloor {\frac {\pi -\mathrm {Im} (z)}{2\pi }}\right\rfloor }
,[ 1] pero posteriormente se cambió a la versión actual porque evita tener que usar de manera reiterada el signo menos (-) en las fórmulas relacionadas.[ 2]
Esta función satisface la relación
W
k
(
z
)
=
ω
(
ln
(
z
)
+
2
π
i
k
)
{\displaystyle W_{k}(z)=\omega (\ln(z)+2\pi ik)}
.
Donde ln(z) indica al logaritmo principal
ln
(
z
)
=
ln
|
z
|
+
i
arg
(
z
)
{\displaystyle \ln(z)=\ln |z|+i\arg(z)}
(arg(z) es el argumento principal ).
ω
:
C
→
C
{\displaystyle \omega :\mathbb {C} \to \mathbb {C} }
es epiyectiva en
C
/
{\displaystyle \mathbb {C} /}
{
0
{\displaystyle 0}
} . "/" se refiere a la diferencia de conjuntos .
ω
:
C
→
C
{\displaystyle \omega :\mathbb {C} \to \mathbb {C} }
es inyectiva si
ω
(
z
)
≠
−
1
{\displaystyle \omega (z)\neq -1}
, ya que
z
=
−
1
±
i
π
→
ω
(
z
)
=
−
1
{\displaystyle z=-1\pm i\pi \to \omega (z)=-1}
.
ω
(
z
)
{\displaystyle \omega (z)}
es continua
∀
z
≠
q
±
i
π
,
q
≤
−
1
{\displaystyle \forall z\neq q\pm i\pi ,\quad q\leq -1}
.
Sea
K
(
z
)
=
⌈
I
m
(
z
)
−
π
2
π
⌉
{\displaystyle \mathrm {K} (z)=\left\lceil {\frac {\mathrm {Im} (z)-\pi }{2\pi }}\right\rceil }
. Entonces
K
(
z
)
=
K
(
ω
(
z
)
)
+
ln
(
ω
(
z
)
)
⟺
z
=
ln
(
ω
(
z
)
)
+
ω
(
z
)
{\displaystyle \mathrm {K} (z)=\mathrm {K} (\omega (z))+\ln(\omega (z))\iff z=\ln(\omega (z))+\omega (z)}
.
Derivación e integración[ editar ]
La función Omega de Wright satisface la ecuación diferencial :
(
1
+
ω
)
d
ω
d
z
=
ω
{\displaystyle (1+\omega ){\frac {d\omega }{dz}}=\omega }
Así:
d
ω
d
z
=
ω
(
z
)
1
+
ω
(
z
)
{\displaystyle {\frac {d\omega }{dz}}={\frac {\omega (z)}{1+\omega (z)}}}
Donde la función sea analítica .
La integral de las funciones de la forma
ω
(
z
)
n
{\displaystyle \omega (z)^{n}}
, viene dado por:
∫
ω
n
d
z
=
{
ω
n
+
1
−
1
n
+
1
+
ω
n
n
si
n
≠
−
1
ln
(
ω
)
−
1
ω
si
n
=
−
1
{\displaystyle \int \omega ^{n}\,dz={\begin{cases}{\frac {\omega ^{n+1}-1}{n+1}}+{\frac {\omega ^{n}}{n}}&{\mbox{si }}n\neq -1\\\ln(\omega )-{\frac {1}{\omega }}&{\mbox{si }}n=-1\end{cases}}}
Una consecuencia de las restricciones sobre la inyectividad y epiyectividad de
ω
{\displaystyle \omega }
, es que la relación inversa
ω
−
1
(
u
)
{\displaystyle \omega ^{-1}(u)}
es función inversa para todos los valores de u menos en 0 y -1. Por lo tanto la función inversa se define como:
ω
−
1
(
u
)
=
{
u
+
ln
(
u
)
−
2
π
i
−
∞
<
u
<
−
1
u
+
ln
(
u
)
en cualquier otro caso
∀
(
u
≠
0
∧
u
≠
−
1
)
{\displaystyle \omega ^{-1}(u)={\begin{cases}u+\ln(u)-2\pi i&-\infty <u<-1\\u+\ln(u)&{\text{en cualquier otro caso}}\quad \forall \ (u\neq 0\quad \land \quad u\neq -1)\end{cases}}}
Nótese que para u = -1,
ω
−
1
{\displaystyle \omega ^{-1}}
toma dos valores,
−
1
±
i
π
{\displaystyle -1\pm i\pi }
(por el hecho de que se podría evaluar de cualquiera de las dos formas de la definición). No obstante, al igual que con otras funciones multivaluadas, se puede adoptar un convenio para denotar un valor principal. Así para términos prácticos,
ω
−
1
(
−
1
)
{\displaystyle \omega ^{-1}(-1)}
solo tomaría un valor -aun no decidido-.
ω
(
0
)
=
W
0
(
1
)
≈
0.56714
ω
(
1
)
=
1
ω
(
1
+
e
)
=
e
ω
(
i
(
1
+
1
2
π
)
)
=
i
ω
(
−
1
+
i
π
)
=
ω
(
−
1
−
i
π
)
=
−
1
ω
(
−
1
3
+
ln
(
1
3
)
+
i
π
)
=
W
0
(
−
1
3
e
−
1
3
)
=
−
1
3
ω
(
−
1
3
+
ln
(
1
3
)
−
i
π
)
=
W
−
1
(
−
1
3
e
−
1
3
)
≈
−
2.23714
ω
(
−
2
+
ln
(
2
)
+
i
π
)
=
W
0
(
−
2
e
−
2
)
≈
−
0.40637
ω
(
−
2
+
ln
(
2
)
−
i
π
)
=
W
−
1
(
−
2
e
−
2
)
=
−
2
{\displaystyle {\begin{array}{lll}\omega (0)&=W_{0}(1)&\approx 0.56714\\\omega (1)&=1&\\\omega (1+e)&=e&\\\omega \left(i\left(1+{\frac {1}{2}}\pi \right)\right)&=i&\\\omega (-1+i\pi )&=\omega (-1-i\pi )&=-1\\\omega (-{\frac {1}{3}}+\ln \left({\frac {1}{3}}\right)+i\pi )&=W_{0}\left(-{\frac {1}{3}}e^{-{\frac {1}{3}}}\right)&=-{\frac {1}{3}}\\\omega (-{\frac {1}{3}}+\ln \left({\frac {1}{3}}\right)-i\pi )&=W_{-1}\left(-{\frac {1}{3}}e^{-{\frac {1}{3}}}\right)&\approx -2.23714\\\omega (-2+\ln(2)+i\pi )&=W_{0}\left(-2e^{-2}\right)&\approx -0.40637\\\omega (-2+\ln(2)-i\pi )&=W_{-1}\left(-2e^{-2}\right)&=-2\\\end{array}}}
La ecuación x + ln(x) = z, tiene por solución:
x
=
{
ω
(
z
)
z
≠
q
±
i
π
,
q
≤
−
1
ω
(
z
)
,
ω
(
z
−
2
π
i
)
z
=
q
+
i
π
,
q
≤
−
1
{\displaystyle x={\begin{cases}\omega (z)&z\neq q\pm i\pi ,&q\leq -1\\\omega (z),\omega (z-2\pi i)&z=q+i\pi ,&q\leq -1\\\end{cases}}}
Para valores
z
=
q
−
i
π
,
q
≤
−
1
{\displaystyle z=q-i\pi ,\quad q\leq -1}
, la ecuación no tiene solución. Cabe destacar que de los dos valores z que satisfacen
ω
(
z
)
=
−
1
{\displaystyle \omega (z)=-1}
, solo
z
=
−
1
+
i
π
{\displaystyle z=-1+i\pi }
cumple
ω
(
z
)
+
ln
(
ω
(
z
)
)
=
z
{\displaystyle \omega (z)+\ln(\omega (z))=z}
, lo cual motiva a que
ω
−
1
(
−
1
)
{\displaystyle \omega ^{-1}(-1)}
-como valor principal- represente a
−
1
+
i
π
{\displaystyle -1+i\pi }
por sobre
−
1
−
i
π
{\displaystyle -1-i\pi }
.
Por otra parte, las dos soluciones complejas -no reales- (conjugadas entre sí) de la ecuación x = ln(x) se pueden representar mediante esta función:
e
−
ω
(
−
i
π
)
=
e
−
W
−
1
(
−
1
)
≈
0.31813
+
1.33723
i
e
−
ω
(
i
π
)
=
e
−
W
0
(
−
1
)
≈
0.31813
−
1.33723
i
{\displaystyle {\begin{array}{lll}e^{-\omega (-i\pi )}&=e^{-W_{-1}(-1)}&\approx 0.31813+1.33723i\\e^{-\omega (i\pi )}&=e^{-W_{0}(-1)}&\approx 0.31813-1.33723i\\\end{array}}}
A continuación se mostrará un ejemplo de cómo resolver ecuaciones usando esta función.
Ejemplo 1
3
x
+
4
ln
(
x
)
=
2
ln
(
x
)
=
1
2
−
3
4
x
x
=
e
1
2
e
−
3
4
x
3
4
x
e
3
4
=
3
4
e
1
2
{\displaystyle {\begin{aligned}3x+4\ln(x)&=2\\\ln(x)&={\frac {1}{2}}-{\frac {3}{4}}x\\x&=e^{\frac {1}{2}}e^{-{\frac {3}{4}}}x\\{\frac {3}{4}}xe^{\frac {3}{4}}&={\frac {3}{4}}e^{\frac {1}{2}}\\\end{aligned}}}
Aplicando
Y
=
X
e
X
⟺
X
=
W
(
Y
)
{\displaystyle Y=Xe^{X}\;\Longleftrightarrow \;X=W(Y)}
3
4
x
=
W
k
(
3
4
e
1
2
)
x
=
4
3
W
k
(
3
4
e
1
2
)
x
=
4
3
ω
(
ln
(
3
4
e
1
2
)
+
2
π
i
k
)
{\displaystyle {\begin{aligned}{\frac {3}{4}}x&=W_{k}\left({\frac {3}{4}}e^{\frac {1}{2}}\right)\\x&={\frac {4}{3}}W_{k}\left({\frac {3}{4}}e^{\frac {1}{2}}\right)\\x&={\frac {4}{3}}\omega \left(\ln \left({\frac {3}{4}}e^{\frac {1}{2}}\right)+2\pi ik\right)\end{aligned}}}
Lo anterior por la primera propiedad. Ahora como 2 es un número positivo y tanto x como ln(x) son inyectivas, la (única) solución buscada es un número real. Sabiendo esto, se puede deducir que k = 0 puesto que
3
4
e
1
2
>
1
{\displaystyle {\frac {3}{4}}e^{\frac {1}{2}}>1}
. Así solo k = 0 cumple lo pedido.
Finalmente la solución es:
x
=
4
3
ω
(
ln
(
3
4
e
1
2
)
)
≈
0.86304
{\displaystyle x={\frac {4}{3}}\omega \left(\ln \left({\frac {3}{4}}e^{\frac {1}{2}}\right)\right)\approx 0.86304}
Función Omega de Wright en los números reales[ editar ]
En los números reales se puede observar mediante la definición, que Im(z) = 0. Por ello, K(z) = 0. Así la función se define como:
ω
(
x
)
=
W
0
(
e
x
)
{\displaystyle \omega (x)=W_{0}(e^{x})}
Con esto la función está definida para todo
R
{\displaystyle \mathbb {R} }
, pero su recorrido es
R
+
{\displaystyle \mathbb {R} ^{+}}
.
Esta función es biyectiva en
R
→
R
+
{\displaystyle \mathbb {R} \to \mathbb {R} ^{+}}
, por lo tanto tiene función inversa, la cual es:
ω
−
1
(
x
)
=
x
+
ln
(
x
)
{\displaystyle \omega ^{-1}(x)=x+\ln(x)}
La cual se deduce por la definición de la función inversa.
Demostración
La inversa de
W
0
(
x
)
:
[
−
1
e
,
+
∞
)
→
[
−
1
,
+
∞
)
{\displaystyle W_{0}(x):\left[-{\frac {1}{e}},+\infty \right)\to \left[-1,+\infty \right)}
es
x
e
x
{\displaystyle xe^{x}}
.
Así que de manera análoga se tiene que al resolver
ω
(
x
)
=
W
0
(
e
x
)
=
y
{\displaystyle \omega (x)=W_{0}(e^{x})=y}
, ocupando lo anterior se tiene (modificando el dominio y recorrido, que es igual al codominio en este caso):
(1 )
e
x
=
y
e
y
{\displaystyle e^{x}=ye^{y}}
Así aplicando ln() en (1):
x
=
ln
(
y
e
y
)
=
ln
(
y
)
+
ln
(
e
y
)
=
ln
(
y
)
+
y
{\displaystyle x=\ln(ye^{y})=\ln(y)+\ln(e^{y})=\ln(y)+y}
Finalmente reemplazando x en y, se tiene que
ω
−
1
(
x
)
=
ln
(
x
)
+
x
{\displaystyle \omega ^{-1}(x)=\ln(x)+x}
Gráficos de la función Omega de Wright en el plano complejo
z = Re(ω(x + i y ))
z = Im(ω(x + i y ))
ω(x + i y )
↑ R.M. Corless, y D.J. Jeffrey, "The unwinding number" , 'Sigsam Bulletin', vol 30, no. 2, junio 1996, 28-35
↑ R.M. Corless, y D.J. Jeffrey, "The Wright ω function" .
, 'AISC (Artificial Intelligence, Automated Reasoning, and Symbolic Computation)', vol. 2385, 2002, 76-89