W функция Ламберта

редактировать

Многозначная «функция» График y = W (x) для действительного x < 6 and y>−4. Верхняя ветвь (синяя) с y ≥ −1 - это график функции W 0 (главная ветвь), нижняя ветвь (пурпурная) с y ≤ −1 - график функции W <103.>−1. Минимальное значение x равно {−1 / e, −1}

. В математике используется функция Ламберта W, также называемая омега-функция или логарифм произведения, является многозначной функцией, а именно ветвями обратной связи функции f (w) = we, где w - любое комплексное число, а e - экспоненциальная функция .

Для каждого целого числа k существует одна ветвь, обозначенная W k (z), которая является комплексной функцией одного комплексного аргумента. W 0 известен как основная ветвь. Эти функции являются следующим своим: если z и w - любые комплексные числа, то

wew = z {\ displaystyle we ^ {w} = z}{ \ displaystyle we ^ {w} = z}

выполняется тогда и только тогда, когда

w = W k (z) для некоторого целого k. {\ displaystyle w = W_ {k} (z) \ \ {\ text {для некоторого целого числа}} k.}{\ displaystyle w = W_ {k} (z) \ \ {\ text {для некоторого целого числа}} k.}

При работе только с действительными числами две ветви W 0 и W −1 достаточно: для действительных чисел x и y уравнение

yey = x {\ displaystyle ye ^ {y} = x}{\ displaystyle ye ^ {y} = x}

может быть решено относительно y, только если x ≥ −1 / e; мы получаем y = W 0 (x), если x ≥ 0, и два значения y = W 0 (x) и y = W −1 ( x), если −1 / e ≤ x < 0.

Отношение Ламберта W не может быть выражено в терминах элементарных функций. Это полезно в комбинаторике, например, при перечислении деревьев. Его можно использовать для решения различных соотношений, включающих экспоненты (например, максимумы распределений Планка, Бозе-Эйнштейна и Ферми-Дирака ), а также решение дифференциальных уравнений с запаздыванием, например y ′ (t) = ay (t - 1). В биохимии и, в частности, в кинетике ферментов, раствор в открытой для анализа форме кинетики времени-изменения кинетики Михаэлиса - Ментен описывается с точки зрения W-функция Ламберта.

Основная ветвь функции Ламберта W на комплексной плоскости. Обратите внимание на отрезок ветви вдоль отрицательной действительной оси, заканчивающийся на -1 / e. На этом рисунке оттенок точки z определяется абсолютным аргументом W (z), а яркость - большим значением W (z).
модуля главной ветви функции Ламберта W, окрашенный в соответствии с arg W (z)

.

Содержание
  • 1 Терминология
  • 2 История
  • 3 Элементарные свойства, ветви и диапазон
    • 3.1 Обратное
  • 4 Исчисление
    • 4.1 Производственное
    • 4.2 Первообразное
  • 5 Асимптотические разложения
    • 5.1 Целочисленные и комплексные степени
  • 6 Идентичности
  • 7 Специальные значения
  • 8 Представления
  • 9 Другие формулы
    • 9.1 Определенные интегралы
    • 9.2 Неопределенные интегралы
  • 10 Приложения
    • 10.1 Решение уравнений
    • 10.2 Вязкие потоки
    • 10.3 Нейровизуализация
    • 10.4 Химическая инженерия
    • 10.5 Материалы
    • 10.6 Пористая среда
    • 10.7 Числа Бернулли и род Тодда
    • 10.8 Статистика
    • 10.9 Точные решения уравнения Шредингера
    • 10.10 Точные решения вакуумных уравнений Эйнштейна
    • 10.11 Резонансы дельты -шель l потенциал
    • 10.12 Термодинамическое равновесие
    • 10.13 Соответствие AdS / CFT
    • 10.14 Эпидемиология
    • 10.15 Определение времени полета снаряда
  • 11 Обобщения
  • 12 Графики
  • 13 Численная оценка
  • 14 Программное обеспечение
  • 15 См. Также
  • 16 Примечания
  • 17 Ссылки
  • 18 Внешние ссылки
Терминология

Функция Ламберта W названа в честь Иоганна Генриха Ламберта. Основная ветвь W 0 ена обозначена Wp в Цифровой библиотеке математических функций, а ветвь W -1 здесь обозначена Wm.

Выбранное здесь обозначение (с W 0 и W −1) следует канонической ссылке на функцию Ламберта W Корлессом, Гоннетом, Хэром, Джеффри и Кнут.

Имя «логарифм произведений» можно понимать так: поскольку обратная функция f (w) = e называется логарифмом, это имеет смысл для возврата обратной функции продукта мы как «логарифм продукта». Это связано с омега-константой, которая равна W 0 (1).

История

Ламберт впервые рассмотрел родственное трансцендентное уравнение Ламберта в 1758 году, что привело к статье Леонарда Эйлера в 1783 году, в котором обсуждался особый случай нас.

Функция, которую рассматривал Ламберт, была

x = x m + q. {\ displaystyle x = x ^ {m} + q.}{\ displaystyle x = x ^ {m} + q.}

Эйлер преобразовал это уравнение в форму

x a - x b = (a - b) c x a + b. {\ displaystyle x ^ {a} -x ^ {b} = (a-b) cx ^ {a + b}.}{\ displaystyle x ^ {a} -x ^ {b} = (ab) cx ^ {a + b}.}

Оба автора получили решение своих уравнений.

После того, как Эйлер решил это уравнение, он рассмотрел случай a = b. Взяв пределы, он вывел уравнение

ln ⁡ x = c x a. {\ displaystyle \ ln x = cx ^ {a}.}{\ displaystyle \ ln х = cx ^ {a}.}

Затем он положил a = 1 и получил решение сходящегося ряда для полученного уравнения, выразив x через c.

После взятия производных по x и некоторых манипуляций получается стандартный вид функции Ламберта.

В 1993 году, когда было объявлено, что W-функция Ламберта обеспечивает точное решение квантово-механической двухъямной модели дельта-функции Дирака для равных зарядов - фундаментальные проблемы физики - Корлесс и разработчики системы компьютерной алгебры Maple провели поиск в библиотеке и представили, что эта функция является повсеместной по своей природе.

Другой пример, где эта функция встречается, находится в Michaelis - Menten кинетика.

Хотя фольклор знал, что функция Ламберта W не может быть выражена в терминах элементарных (лиувиллевских) функций, первое опубликованное доказательство появилось только в 2008 году.

Элементарные свойства, ветви и диапазон
Диапазон функции W с отображением всех ветвей. Черные кривые (включая действительную ось) образуют изображение обратной оси, оранжевые кривые - изображение мнимой оси. Фиолетовая кривая - это изображение небольшого круга вокруг точки z = 0; красные кривые - это изображение маленького круга вокруг точки z = −1 / e. График мнимой части W [n, x + iy] для ветвей n = -2, -1,0,1, 2. График аналогичен графику многозначной функции комплексного логарифма , за исключением того, что расстояние между листами не является постоянным, а соединение основного листа отличается

Существует счетное множество ветвей функции W, обозначенных на W k (z) для целого k; W 0 (z) - главная или основная ветвь. W 0 (z) определен для всех комплексных чисел z, тогда как W k (z) с k ≠ 0 определен для всех ненулевых z. Имеем W 0 (0) = 0 и limz → 0 W k (z) = −∞ для всех k ≠ 0.

Точка ветвления для главной ветвь находится в точке z = −1 / e, с разрезом ветви, который продолжается до −∞ вдоль отрицательной действительной оси. Этот разрез ветвей отделяет главную ветвь от двух ветвей W -1 и W 1. Во всех ветвях W k с k ≠ 0 имеется точка ветвления в точке z = 0 и ветвь, разрезанная вдоль всей отрицательной действительной оси.

Функции W k (z), k ∈ Z все инъективны и их диапазоны не пересекаются. Образом всей многозначной функции W является комплексная плоскость. Изображение действительной оси представляет собой объединение действительной и квадратисы Гиппия, параметрической кривой w = −t cot t + it.

Обратный

Области комплексной плоскости, для которых W (n, z e z) = z {\ displaystyle W (n, ze ^ {z}) = z}{\ displaystyle W (n, ze ^ {z}) = z} . Более темные границы того же цвета включаются в более светлую область того же цвета. Точка в {-1,0} входит как в область n = - 1 {\ displaystyle n = -1}n = -1 (синий), так и в область n = 0 {\ displaystyle n = 0}n = 0 (серая) область. Горизонтальные линии сетки кратны π.

График вышеупомянутой области в плоскости плоскости, где простая обратная 'W (n, zez) = z {\ displaystyle W (n, ze ^ {z}) = z}{\ displaystyle W (n, ze ^ {z}) = z} верно. f = ze означает, что существует n такое, что z = W (n, f) = W (n, zez) {\ displaystyle z = W (n, f) = W (n, ze ^ {z})))}{\ displaystyle z = W (n, f) = W (n, ze ^ {z})} , где n будет зависеть от значений z. Значение целого числа n резко изменится, когда ze станет означать на отрезке ветви W (n, zez) {\ displaystyle W (n, ze ^ {z})}{ \ displaystyle W (n, ze ^ {z})} , что будет означать, что ze ≤ 0, за исключением n = 0 {\ displaystyle n = 0}n = 0 , где ze ≤ -1 / e.

Определите z = x + i y {\ displaystyle z = x + iy}z = x + iy , где x и y действительны. Выражая e в полярных координатах, видно, что:

zez = (x + iy) ex (cos ⁡ (y) + i sin ⁡ (y)) = ex (x cos ⁡ (y) - y sin ⁡ (y)) + iex (x грех ⁡ (y) + y соз ⁡ (y)) {\ displaystyle {\ begin {align} ze ^ {z} = (x + iy) e ^ {x} (\ cos (y) + i \ sin (y)) \\ = e ^ {x} (x \ cos (y) -y \ sin (y)) + ie ^ {x} (x \ sin (y) + y \ cos (y)) \\\ конец {выровнен}}}{\ displaystyle {\ begin {align} ze ^ {z} = (x + iy) e ^ {x} (\ cos (y) + i \ sin (y)) \\ = e ^ {x} (x \ cos (y) -y \ sin (y)) + ie ^ {x} (x \ sin (y) + y \ cos (y)) \ \\ конец {выровнен}}}

Для n ≠ 0 {\ displaystyle n \ neq 0}n \ neq 0 , отрезок ветви для W [n, zez] {\ displaystyle W [n, ze ^ {z}]}{\ displaystyle W [n, ze ^ {z}]} будет неположительной действительной осью, так что:

x sin ⁡ (y) + y cos ⁡ (y) = 0 ⇒ x = - y / загар ⁡ (y) {\ displaystyle x \ sin (y) + y \ cos (y) = 0 \ Rightarrow x = -y / \ tan (y)}{\ displaystyle x \ sin (y) + y \ cos (y) = 0 \ Rightarrow x = -y / \ tan (y)}

и

(Икс соз ⁡ ( Y) - Y грех ⁡ (Y)) ex ≤ 0 {\ Displaystyle (x \ cos (y) -y \ sin (y)) e ^ {x} \ Leq 0}{\ displaystyle (x \ cos (y) -y \ sin (y)) e ^ {x} \ leq 0}

Для n = 0 {\ displaystyle n = 0}n = 0 , отрезок ветви для W [n, zez] {\ displaystyle W [n, ze ^ {z}]}{\ displaystyle W [n, ze ^ {z}]} будет действительной осью с - ∞ < z ≤ − 1 / e {\displaystyle -\infty {\ displaystyle - \ infty <z \ leq -1 / е} , так что неравенство будет имет ь вид:

(x соз ⁡ (y) - y sin ⁡ (y)) ex ≤ - 1 / e {\ dis playstyle (x \ cos (y) -y \ sin (y)) e ^ {x } \ leq -1 / e}{\ displaystyle (x \ cos (y) -y \ sin (y)) e ^ {x} \ leq -1 / e}

Внутри границ, ограниченных указанным выше, не будет прерывистых изменений в W (n, zez) {\ displaystyle W (n, ze ^ {z})}{ \ displaystyle W (n, ze ^ {z})} , и эти области будут указывать, где функция W просто обратима: т.е. W (n, zez) = z {\ displaystyle W (n, ze ^ {z}) = z}{\ displaystyle W (n, ze ^ {z}) = z} .

Исчисление

Производная

С помощью неявного дифференцирования можно показать, что все ветви W удовлетворяет дифференциальному уравнению

z (1 + W) d W dz = W для z - 1 e. {\ displaystyle z (1 + W) {\ frac {dW} {dz}} = W \ quad {\ text {for}} z \ neq - {\ frac {1} {e}}.}{\ displaystyle z (1 + W) {\ frac {dW} {dz}} = W \ quad {\ text {for}} z \ neq - {\ frac {1} {e}}.}

( W не дифференцируемо для z = −1 / e.) Как следствие, мы получаем формулу для производной W:

d W dz = W (z) z (1 + W (z)) для z ∉ {0, - 1 e}. {\ displaystyle {\ frac {dW} {dz}} = {\ frac {W (z)} {z (1 + W (z))}} \ quad {\ text {for}} z \ not \ in \ left \ {0, - {\ frac {1} {e}} \ right \}.}{\ displaystyle {\ frac {dW} {dz}} = {\ frac {W (z)} {z (1 + W (z))}} \ quad {\ text {for}} z \ not \ in \ left \ {0, - {\ frac {1} {e}} \ right \}.}

Используя тождество e = z / W (z), мы получаем эквивалентную формулу:

d W dz = 1 z + e W (z) для z ≠ - 1 e. {\ displaystyle {\ frac {dW} {dz}} = {\ frac {1} {z + e ^ {W (z)}}} \ quad {\ text {for}} z \ neq - {\ frac { 1} {e}}.}{\ displaystyle {\ frac {dW} {dz}} = {\ frac {1} {z + e ^ {W (z)}}} \ quad {\ text {for}} z \ neq - {\ frac {1} {e}}.}

В начале координат мы имеем

W 0 ′ (0) = 1. {\ displaystyle W '_ {0} (0) = 1.}{\displaystyle W'_{0}(0)=1.}

Первообразная

Функция W (x) и многие выражения, включающие W (x), могут быть интегрированы с помощью подстановки w = W (x), т.е. x = мы:

∫ W (x) dx = x W (x) - x + e W (x) + C = x (W (x) - 1 + 1 W (x)) + C. {\ Displaystyle {\ begin {align} \ int W (x) \, dx = xW (x) -x + e ^ {W (x)} + C \\ = x \ left (W (x) -1 + {\ frac {1} {W (x)}} \ right) + C. \ end {align}}}{\ displaystyle {\ begin {align} \ int W (x) \, dx = xW (x) -x + e ^ {W (x)} + C \\ = x \ left (W (x) -1 + {\ frac {1} {W (x)}} \ right) + C. \ End {align}} }

(Последнее уравнение часто встречается в литературе, но не выполняется при x = 0). Одним из следствий этого (используя тот факт, что W 0 (e) = 1) является тождеством

∫ 0 e W 0 (x) dx = e - 1. {\ displaystyle \ int _ { 0} ^ {e} W_ {0} (x) \, dx = e-1.}{\ displaystyle \ int _ { 0} ^ {e} W_ {0} (x) \, dx = e-1.}
Асимптотические разложения

Ряд Тейлора W 0 около 0 можно найти с помощью теоремы об обращении Лагранжа и определяется как

W 0 (x) = ∑ n = 1 ∞ (- n) n - 1 n! Икс N знак равно Икс - Икс 2 + 3 2 Икс 3-8 3 Икс 4 + 125 24 Х 5 - ⋯. {\ displaystyle W_ {0} (x) = \ sum _ {n = 1} ^ {\ infty} {\ frac {(-n) ^ {n-1}} {n!}} x ^ {n} = xx ^ {2} + {\ tfrac {3} {2}} x ^ {3} - {\ tfrac {8} {3}} x ^ {4} + {\ tfrac {125} {24}} x ^ {5} - \ cdots.}{\ displaystyle W_ {0} (x) = \ sum _ {n = 1} ^ {\ infty} {\ frac {(-n) ^ { n-1}} {n!}} x ^ {n} = xx ^ {2} + {\ tfrac {3} {2}} x ^ {3} - {\ tfrac {8} {3}} x ^ {4} + {\ tfrac {125} {24}} x ^ {5} - \ cdots.}

Радиус сходимости равенство 1 / e, что можно увидеть с помощью теста отношения . Функция, определяемая этим рядом, может быть расширена до голоморфной функции , определенные на всех комплексных числах с отрезком ветви вдоль интервала (−∞, −1 / e] ; эта голоморфная функция определяет главная ветвь функции Ламберта W.

Для больших значений x W 0 асимптотичен

W 0 (x) = L 1 - L 2 + L 2 L 1 + L 2 (- 2 + L 2) 2 L 1 2 + L 2 (6-9 L 2 + 2 L 2 2) 6 L 1 3 + L 2 (- 12 + 36 L 2 - 22 L 2 2 + 3 L 2 3) 12 L 1 4 + ⋯ = L 1 - L 2 + ∑ l = 0 ∞ ∑ m = 1 ∞ (- 1) l [l + ml + 1] ì! L 1 - l - м L 2 м, {\ displaystyle {\ begin {выровнено} W_ {0} (x) = L_ {1} -L_ {2} + {\ frac {L_ {2}} {L_ {1}}} + {\ frac {L_ {2} \ left (-2 + L_ {2} \ right)} {2L_ {1} ^ {2}}} + {\ frac {L_ {2} \ left (6-9L_ { 2} + 2L_ {2} ^ {2} \ right)} {6L_ {1} ^ {3}}} + {\ frac {L_ {2} \ left (-12+ 36L_ {2} -22L_ {2} ^ {2} + 3L_ {2} ^ {3} \ right)} {12L_ {1} ^ {4}}} + \ cdots \\ [5pt] = L_ {1} -L_ {2} + \ sum _ {l = 0} ^ {\ infty} \ sum _ {m = 1} ^ {\ inft y} {\ frac {(-1) ^ {l} \ left [{\ begin {smallmatrix} l + m \\ l + 1 \ end {smallmatrix}} \ right]} {m!}} L_ {1} ^ {- lm} L_ {2} ^ {m}, \ end {align}}}{\ displaystyle {\ begin {align} W_ {0} (x) = L_ {1} -L_ {2} + {\ frac {L_ { 2}} {L_ {1}}} + {\ frac {L_ {2} \ left (-2 + L_ {2} \ right)} {2L_ {1} ^ {2}}} + {\ frac {L_ {2} \ left (6-9L_ {2} + 2L_ {2} ^ {2} \ right)} {6L_ {1} ^ {3}}} + {\ fra c {L_ {2} \ left (- 12 + 36L_ {2} -22L_ {2} ^ {2} + 3L_ {2} ^ {3} \ right)} {12L_ {1} ^ {4}}} + \ cdots \\ [5pt] = L_ {1} -L_ {2} + \ sum _ {l = 0} ^ {\ infty} \ sum _ {m = 1} ^ {\ infty} {\ frac {(-1) ^ {l} \ left [ {\ begin {smallmatrix} l + m \\ l + 1 \ end {smallmatrix}} \ right]} {m!}} L_ {1} ^ {- lm} L_ {2} ^ {m}, \ end { align}}}

где L 1 = ln x, L 2 = ln ln x, и [. l + 1 ]- неотрицательное число Стирлинга первого рода. Сохраняя только первые два членских разложения,

W 0 (x) = ln ⁡ x - ln ⁡ ln ⁡ x + o (1). {\ displaystyle W_ {0} (x) = \ ln x- \ ln \ ln x + o (1).}{\ displaystyle W_ {0} (x) = \ ln x- \ ln \ ln х + о (1).}

Другая действительная ветвь, W −1, определенная в интервале [−1 / e, 0), имеет приближение того же вида, когда x стремится к нулю, в этом случае L 1 = ln (−x) и L 2 = ln (−ln ( −x)).

Показано, что имеет место следующая оценка (верхняя оценка только для x ≥ e):

ln ⁡ x - ln ⁡ ln ⁡ x + ln ⁡ ln ⁡ Икс 2 ln ⁡ x ≤ W 0 (x) ≤ ln ⁡ x - ln ⁡ ln ⁡ x + ee - 1 ln ⁡ ln ⁡ x ln ⁡ x. {\ Displaystyle \ пер х- \ пер \ пер Икс + {\ гидроразрыва {\ пер \ пер х} {2 \ пер х}} \ leq W_ {0} (х) \ Leq \ пер х- \ пер \ пер х + {\ frac {e} {e-1}} {\ frac {\ ln \ ln x} {\ ln x}}.}{\ displaystyle \ ln x- \ ln \ ln x + {\ frac {\ ln \ ln x} {2 \ ln x}} \ leq W_ {0} ( x) \ leq \ ln x- \ ln \ ln x + {\ frac {e} {e-1}} {\ frac {\ ln \ ln x} {\ ln x}}.}

В 2013 году было доказано, что ветвь W −1 может быть ограниченным следующим образом:

- 1 - 2 u - u < W − 1 ( − e − u − 1) < − 1 − 2 u − 2 3 u for u>0. {\ displaystyle -1 - {\ sqrt {2u}} - u 0.}{\displaystyle -1-{\sqrt {2u}}-u<W_{-1}\left(-e^{-u-1}\right)<-1-{\sqrt {2u}}-{\tfrac {2}{3}}u\quad {\text{for }}u>0.}

Целочисленные и комплексные степени

Целочисленные вложения степени W 0 Разные степени W 0 ряд Тейлора (или Лорана ) в нуле:

W 0 (x) 2 = ∑ n = 2 ∞ - 2 (- n) n - 3 (n - 2)! xn знак равно x 2 - 2 x 3 + 4 x 4 - 25 3 x 5 + 18 x 6 - ⋯. {\ displaystyle W_ {0} (x) ^ {2} = \ sum _ {n = 2} ^ {\ infty } {\ frac {-2 \ left (-n \ right) ^ {n-3}} {(n-2)!}} x ^ {n} = x ^ {2} -2x ^ {3} + 4x ^ {4} - {\ tfrac {25} {3}} x ^ {5} + 18x ^ {6} - \ cdots.}{\ displaystyle W_ {0} (x) ^ {2} = \ sum _ {n = 2} ^ {\ infty} {\ frac {-2 \ left (-n \ right) ^ {n-3}} {(n-2)!}} x ^ {n} = x ^ {2} -2x ^ {3} + 4x ^ {4} - {\ tfrac {25} {3}} x ^ {5} + 18x ^ {6} - \ cdots.}

В общем, для r ∈ ℤ Формула обращения Лагранжа дает

W 0 (x) r = ∑ n = r ∞ - r (- n) n - r - 1 (n - r)! Xn, {\ displaystyle W_ {0} (x) ^ {r} = \ sum _ {n = r} ^ {\ infty} {\ frac {-r \ left (-n \ right) ^ {nr-1}} {(nr)!}} x ^ { n},}{\ displaystyle W_ {0} (x) ^ {r} = \ sum _ {n = r} ^ {\ infty} {\ frac {-r \ left (-n \ right) ^ {nr-1}} {(nr)!}} x ^ {n},}

который, в общем, является рядо м Лорана порядка r. Эквивалентно, последний может быть записан в форме Taylor ex расширение степеней W 0 (x) / x:

(W 0 (x) x) r = e - r W 0 (x) = ∑ n = 0 ∞ г (п + г) п - 1 п! (- х) n, {\ displaystyle \ left ({\ frac {W_ {0} (x)} {x}} \ right) ^ {r} = e ^ {- rW_ {0} (x)} = \ sum _ {n = 0} ^ {\ infty} {\ frac {r \ left (n + r \ right) ^ {n-1}} {n!}} \ left (-x \ right) ^ {n},}{\ displaystyle \ left ({\ frac {W_ {0} (x)} {x}} \ right) ^ {r} = e ^ {- rW_ {0} (x)} = \ sum _ {n = 0} ^ {\ infty} {\ frac {r \ left (n + r \ right) ^ {n-1}} {n!}} \ left (-x \ right) ^ {n},}

это происходит для любого r ∈ ℂ и | х | < 1/e.

Identities
График W j (x e), где синий - для j = 0, а красный - для j = -1. Диагональная линия представляет интервалы, где W j (xe) = x

Из определения следует несколько тождеств:

W 0 (xex) = x для x ≥ - 1, W - 1 (xex) = x для x ≤ - 1. {\ displaystyle {\ begin {align} W_ {0} (xe ^ {x}) = x {\ text {for}} x \ geq -1, \\ W_ { -1} (xe ^ {x}) = x {\ text {for}} x \ leq -1. \ End {align}}}{\ displaystyle {\ begin {align} W_ {0} (xe ^ {x}) = x {\ text {for}} x \ geq -1, \\ W _ {- 1} (xe ^ {x}) = x {\ text {для}} x \ leq -1. \ end {align}}}

Обратите внимание, что, поскольку f (x) = xe не injective, не всегда верно, что W (f (x)) = x, как и в случае с обратными тригонометрическими функциями . Для фиксированного x < 0 and x ≠ −1, the equation xe = ye has two solutions in y, one of which is of course y = x. Then, for i = 0 and x < −1, as well as for i = −1 and x ∈ (−1, 0), y = Wi (xe) - другое решение.

Некоторые другие тождества:

W (x) ⋅ e W (x) = x, следовательно: e W (x) = x W (x), e - W (x) = W (x) x, en W (x) = (x W (x)) n. {\ Displaystyle {\ begin {align} W (x) \ cdot e ^ {W (x)} = x, \ quad {\ text {, следовательно:}} \\ [5pt] e ^ {W (x)} = {\ frac {x} {W (x)}}, \ qquad e ^ {- W (x)} = {\ frac {W (x)} {x}}, \ qquad e ^ {nW ( x)} = \ left ({\ frac {x} {W (x)}} \ right) ^ {n}. \ end {align}}}{\ displaystyle {\ begin {align} W (x) \ cdot e ^ {W (x)} = x, \ quad {\ text {поэтому: }} \\ [5pt] e ^ {W (x)} = {\ frac {x} {W (x)}}, \ qquad e ^ {- W (x)} = {\ frac {W (x)} {x}}, \ qquad e ^ {nW (x)} = \ left ({\ frac {x} {W (x)}} \ right) ^ {n}. \ end {align}}}
ln ⁡ W 0 (x) = ln ⁡ x - W 0 (x) для x>0. {\ displaystyle \ ln W_ {0} (x) = \ ln x-W_ {0} (x) \ quad {\ text {for}} x>0.}{\displaystyle \ln W_{0}(x)=\ln x-W_{0}(x)\quad {\text{for }}x>0.}
W 0 (x ln ⁡ x) = пер ⁡ Икс и е W 0 (Икс пер ⁡ Икс) = Икс для Икс ≥ 1 е. {\ Displaystyle W_ {0} \ left (x \ ln x \ right) = \ ln x \ quad {\ text {и} } \ quad e ^ {W_ {0} \ left (x \ ln x \ right)} = x \ quad {\ text {for}} x \ geq {\ frac {1} {e}}.}{\ displaystyle W_ {0} \ left (x \ ln x \ right) = \ ln x \ quad {\ text {и}} \ quad e ^ {W_ { 0} \ left (x \ ln x \ right)} = x \ quad {\ text {for}} x \ geq {\ frac {1} {e}}.}
W - 1 (Икс пер ⁡ Икс) знак равно пер ⁡ Икс и е W - 1 (Икс пер ⁡ Икс) = Икс для Икс ≤ 1 е. {\ Displaystyle W _ {- 1} \ влево (х \ пер х \ вправо) = \ ln x \ quad {\ text {and}} \ quad e ^ {W _ {- 1} \ left (x \ ln x \ right)} = x \ quad {\ text {for}} x \ leq {\ frac {1} {e}}.}{\ displaystyle W _ {- 1} \ left (x \ ln x \ right) = \ ln x \ quad {\ text {and}} \ quad e ^ {W _ {- 1 } \ left (x \ ln x \ right)} = x \ quad {\ text {for}} x \ leq {\ frac {1} {e}}.}
W (x) = ln ⁡ x W (x) для x ≥ - 1 e, W (nxn W (x) n - 1) = n W ( x) для n, x>0 {\ displaystyle {\ begin {align} W (x) = \ ln {\ frac {x} {W (x)}} {\ text {for}} x \ geq - {\ frac {1} {e}}, \\ [5pt] W \ left ( {\ frac {nx ^ {n}} {W \ left (x \ right) ^ {n-1}}} \ right) = nW (x) {\ text {for}} n, x>0 \ end {align}}}{\displaystyle {\begin{aligned}W(x)=\ln {\frac {x}{W(x)}}{\text{for }}x\geq -{\frac {1}{e}},\\[5pt]W\left({\frac {nx^{n}}{W\left(x\right)^{n-1}}}\right)=nW(x){\text{for }}n,x>0 \ end {align}}}
(который может быть расширен на другие n и x, если правильная ветка выбрано).
W (x) + W (y) = W (xy (1 W (x) + 1 W (y))) для x, y>0. {\ Displaystyle W (x) + W (y) = W \ left (xy \ left ({\ frac {1} {W (x)}} + {\ frac {1} {W (y)}} \ справа))) \ right) \ quad {\ text {for}} x, y>0.}{\displaystyle W(x)+W(y)=W\left(xy\left({\frac {1}{W(x)}}+{\frac {1}{W(y)}}\right)\right)\quad {\text{for }}x,y>0.}

Подставляя −ln x в определение:

W 0 (- ln ⁡ xx) = - е. {\ Displaystyle {\ begin {align} W_ {0} \ left (- {\ frac {\ ln x} {x}} \ right) = - \ ln x {\ text {for}} 0 е. \ End {выравнивание}}}{\displaystyle {\begin{aligned}W_{0}\left(-{\frac {\ln x}{x}}\right)=-\ln x{\text{for }}0<x\leq e,\\[5pt]W_{-1}\left(-{\frac {\ln x}{x}}\right)=-\ln x{\text{for }}x>е. \ End {выравнивание}}}

С повторяющейся экспонентой Эйлера h (x):

h (x) = e - W (- ln ⁡ x) = W (- пер ⁡ Икс) - пер ⁡ Икс для Икс ≠ 1. {\ Displaystyle {\ begin {Выровнено} h (x) = e ^ {- W (- \ ln x)} \\ = {\ frac { W (- \ ln x)} {- \ ln x}} \ quad {\ text {for}} x \ ne q 1. \ end {align}}}{\ displaystyle {\ begin {align} h (x) = e ^ {- W (- \ ln x)} \\ = {\ frac { W (- \ ln x)} {- \ ln x}} \ quad {\ text {for}} x \ neq 1. \ end {align}}}
Специальные значения

Для любого ненулевого алгебраическое число x, W (x) трансцендентным номером. В самом деле, если W (x) равно нулю, то x также должно быть равен нулю, а если W (x) ненулевой и алгебраический, то по теореме Линдемана - Вейерштрасса, e должен быть трансцендентным, что означает, что x = W (x) также должен быть трансцендентным.

Следующие специальные значения главной ветви:

W (- π 2) = i π 2. {\ displaystyle W \ left (- {\ frac {\ pi} {2}} \ right) = { \ frac {i \ pi} {2}}.}{\ displaystyle W \ left (- {\ frac {\ pi} {2}} \ right) = {\ frac {i \ pi} {2}}. }
W (- 1 e) = - 1. {\ displaystyle W \ left (- {\ frac {1} {e}} \ right) = - 1.}{\ displaystyle W \ left (- {\ frac {1} {e}} \ right) = - 1.}
W (- ln ⁡ aa) = - ln ⁡ a (1 e ≤ a ≤ e). {\ displaystyle W \ left (- {\ frac {\ ln a} {a}} \ right) = - \ ln a \ quad \ left ({\ frac {1} {e}} \ leq a \ leq e \ справа).}{\ displaystyle W \ left (- {\ frac {\ ln a} {a} } \ right) = - \ ln a \ quad \ left ({\ frac {1} {e}} \ leq a \ leq e \ right).}
W (2 пер ⁡ 2) = пер ⁡ 2. {\ displaystyle W \ left (2 \ ln 2 \ right) = \ ln 2.}{\ displaystyle W \ left (2 \ ln 2 \ right) = \ ln 2.}
W (0) = 0. { \ Displaystyle W (0) = 0.}{\ displaystyle W (0) = 0.}
W (1) = Ω = (∫ - ∞ ∞ dt (et - t) 2 + π 2) - 1-1 ≈ 0,56714329… {\ displaystyle W ( 1) = \ Omega = \ left (\ int _ {- \ infty} ^ {\ infty} {\ frac {dt} {\ left (e ^ {t} -t \ right) ^ {2} + \ pi ^ {2}}} \ right) ^ {- 1} -1 \ приблизительно 0,56714329 \ ldots}{\ displaystyle W (1) = \ Omega = \ left (\ int _ {- \ infty} ^ {\ infty} {\ frac {dt} {\ left (e ^ {t} -t \ right) ^ {2} + \ pi ^ {2}}} \ справа) ^ {- 1} -1 \ приблизительно 0,56714329 \ ldots} (постоянная омега ).
W (1) = e - W (1) = пер ⁡ (1 W (1)) знак равно - пер ⁡ W (1). {\ Displaystyle W (1) = e ^ {- W (1)} = \ ln \ left ({\ frac {1} {W (1)}} \ right) = - \ ln W (1).}{\ Displaystyle W (1) = е ^ {- W (1)} = \ ln \ left ({\ frac {1} {W (1)}} \ right) = - \ ln W (1). }
W (e) = 1. {\ displaystyle W (e) = 1.}{\ displaystyle W (e) = 1}
W (e 1 + e) ​​Знак равно е. {\ Displaystyle W \ left (e ^ {1 + e} \ right) = e.}{\ displaystyle W \ left (e ^ {1 + e} \ right) = e.}
W (- 1) ≈ - 0,31813 + 1,33723 i. {\ Displaystyle W (- 1) \ приблизительно -0,31813 + 1,33723i.}{\ displaystyle W (-1) \ приблизительно -0,31813 + 1, 33723i.}
Представления

Основныевь функц ии Лам берта может быть представлен интегралом благодаря Пуассону:

- π 2 W (- x) = ∫ 0 π sin ⁡ (3 2 т) - х e cos ⁡ t sin ⁡ (5 2 t - sin ⁡ t) 1 - 2 xe cos ⁡ t cos ⁡ (t - sin ⁡ t) + x 2 e 2 cos ⁡ t sin ⁡ (1 2 t) dt для | х | < 1 e. {\displaystyle -{\frac {\pi }{2}}W(-x)=\int _{0}^{\pi }{\frac {\sin \left({\tfrac {3}{2}}t\right)-xe^{\cos t}\sin \left({\tfrac {5}{2}}t-\sin t\right)}{1-2xe^{\cos t}\cos(t-\sin t)+x^{2}e^{2\cos t}}}\sin \left({\tfrac {1}{2}}t\right)\,dt\quad {\text{for }}|x|<{\frac {1}{e}}.}{\ displaystyle - {\ frac {\ pi} {2}} W (-x) = \ int _ {0} ^ {\ pi} {\ frac {\ sin \ left ({\ tfrac {3} {2}} t \ right) -xe ^ {\ cos t} \ sin \ left ({\ tfrac {5} {2}} t- \ sin t \ right)} {1- 2xe ^ {\ cos t} \ cos (t- \ sin t) + x ^ {2} e ^ {2 \ cos t}}} \ sin \ left ({\ tfrac {1} {2}} t \ right) \, dt \ quad {\ text {for}} | х | <{\ frac {1} {e}}.}

В более широкой области −1 / e ≤ x ≤ e значительно более простое представление находится с помощью Mező:

W (x) = 1 π Re ⁡ ∫ 0 π ln ⁡ (eeit - xe - iteeit - xeit) dt. {\ displaystyle W (x) = {\ frac {1} {\ pi}} \ operatorname {Re} \ int _ {0} ^ {\ pi} \ ln \ left ({\ frac {e ^ {e ^ { it}} - xe ^ {- it}} {e ^ {e ^ {it}} - xe ^ {it}}} \ right) \, dt.}{\ displaystyle W (x) = {\ frac {1} {\ pi}} \ operatorname {Re} \ int _ {0} ^ {\ pi} \ ln \ left ({\ frac {e ^ {e ^ {it}} - xe ^ {- it}} {e ^ {e ^ {it}} - xe ^ {it}}} \ right) \, dt.}

Следующая непрерывная дробь представление также справедливо для главной ветви:

W (x) = x 1 + x 1 + x 2 + 5 x 3 + 17 x 10 + 133 x 17 + 1927 x 190 + 13582711 x 94423 + ⋱. {\ Displaystyle W (x) = {\ cfrac {x} {1 + {\ cfrac {x} {1 + {\ cfrac {x} {2 + {\ cfrac {5x} {3 + {\ cfrac {17x}))) {10 + {\ cfrac {133x} {17 + {\ cfrac {1927x} {190 + {\ cfrac {13582711x} {94423+ \ ddots}}}}}}}}}}}}}}}}. }{\ displaystyle W (x) = {\ cfrac {x} {1 + {\ cfrac {x} {1 + {\ cfrac {x} {2 + {\ cfrac {5x} {3 + { \ cf) rac {17x} {10 + {\ cfrac {133x} {17 + {\ cfrac {1927x} {190 + {\ cfrac {13582711x} {94423+ \ ddots}}}}}}}}}}}} }}}}}}.}

Также, если | W (z) | < 1:

W (x) = x exp ⁡ x exp ⁡ x ⋱. {\ displaystyle W (x) = {\ cfrac {x} {\ exp {\ cfrac {x} {\ exp {\ cfrac {x} {\ ddots}}}}}}.}{\ displaystyle W (x) = {\ cfrac {x} {\ exp {\ cfrac {x} {\ exp {\ cfrac {x} {\ ddots}}}}}}.}

В свою очередь, если | W (z) |>е, то

W (x) = ln ⁡ x ln ⁡ x ln ⁡ x ⋱. {\ displaystyle W (x) = \ ln {\ cfrac {x} {\ ln {\ cfrac {x} {\ ln {\ cfrac {x} {\ ddots}}}}}}.}{\ displaystyle W (x) = \ ln {\ cfrac {x} {\ ln {\ cfrac {x} {\ ln {\ cfrac {x} {\ ddots}}}}}}.}

.

Другие формулы

Определенные интегралы

Существует несколько полезных формул определенного интеграла, включающих основные ветвь функции W, в том числе следующие:

∫ 0 π W (2 cot 2 ⁡ x) sec 2 ⁡ xdx = 4 π. ∫ 0 ∞ W (Икс) Икс Икс d Икс знак равно 2 2 π. ∫ 0 ∞ W (1 x 2) d x = 2 π. {\ displaystyle {\ begin {align} \ int _ {0} ^ {\ pi} W \ left (2 \ cot ^ {2} x \ right) \ sec ^ {2} x \, dx = 4 {\ sqrt {\ pi}}. \\ [5pt] \ int _ {0} ^ {\ infty} {\ frac {W (x)} {x {\ sqrt {x}}}} \, dx = 2 {\ sqrt {2 \ pi} }. \\ [5pt] \ int _ {0} ^ {\ infty} W \ left ({\ frac {1} {x ^ {2}}} \ right) \, dx = {\ sqrt {2 \ pi} }. \ end {align}}}{\ displaystyle {\ begin {align} \ int _ {0} ^ {\ pi} W \ left (2 \ cot ^ {2} x \ right) \ sec ^ {2} x \, dx = 4 {\ sqrt {\ pi}}. \\ [5pt] \ int _ {0} ^ {\ infty} {\ frac {W (x)} {x {\ sqrt {x}}}} \, dx = 2 {\ sqrt {2 \ pi} }. \\ [5pt] \ int _ {0} ^ {\ infty} W \ left ({\ frac {1} {x ^ {2}}} \ right) \, dx = {\ sqrt {2 \ pi} }. \ End {align}}}

Первую идентичность можно найти, записав интеграл Гаусса в полярных координат.

Вторая идентичность можно получить, сделав замену u = W (x), которая дает

х = ueu, dxdu = (u + 1) eu. {\ displaystyle {\ begin {align} x = ue ^ {u}, \\ [5pt] {\ frac {dx} {du}} = (u + 1) e ^ {u}. \ end {выравнивается}}}{\ displaystyle {\ begin {align} x = ue ^ {u}, \ \ [5pt] {\ frac {dx} {du}} = (u + 1) e ^ {u}. \ end {align}}}

Таким образом,

∫ 0 ∞ W (x) xxdx = ∫ 0 ∞ uueuueu (u + 1) eudu = ∫ 0 ∞ u + 1 ueudu = ∫ 0 ∞ u + 1 u 1 eudu = ∫ 0 ∞ u 1 2 e - u 2 du + ∫ 0 ∞ u - 1 2 e - u 2 du = 2 ∫ 0 ∞ (2 w) 1 2 e - wdw + 2 ∫ 0 ∞ (2 w) - 1 2 e - wdw (u = 2 w) = 2 2 ∫ 0 ∞ w 1 2 e - wdw + 2 ∫ 0 ∞ w - 1 2 e - wdw = 2 2 ⋅ Γ (3 2) + 2 ⋅ Γ (1 2) = 2 2 (1 2 π) + 2 (π) = 2 2 π. {\ displaystyle {\ begin {align} \ int _ {0} ^ {\ infty} {\ frac {W (x)} {x {\ sqrt {x}}}} \, dx = \ int _ {0 } ^ {\ infty} {\ frac {u} {ue ^ {u} {\ sqrt {ue ^ {u}}}}} (u + 1) e ^ {u} \, du \\ [5pt] = \ int _ {0} ^ {\ infty} {\ frac {u + 1} {\ sqrt {ue ^ {u}}}} du \\ [5pt] = \ int _ {0} ^ {\ infty } {\ frac {u + 1} {\ sqrt {u}}} {\ frac {1} {\ sqrt {e ^ {u}}}} du \\ [5pt] = \ int _ {0} ^ {\ infty} u ^ {\ tfrac {1} {2}} e ^ {- {\ frac {u} {2}}} du + \ int _ {0} ^ {\ infty} u ^ {- {\ tfrac {1} {2}}} e ^ {- {\ frac {u} {2}}} du \\ [5pt] = 2 \ int _ {0} ^ {\ infty} (2w) ^ {\ tfrac {1} {2}} e ^ {- w} \, dw + 2 \ int _ {0} ^ {\ infty} (2w) ^ {- {\ tfrac {1} {2}}} e ^ { - w} \, dw \ quad (u = 2w) \\ [5pt] = 2 {\ sqrt {2}} \ int _ {0} ^ {\ infty} w ^ {\ tfrac {1} {2 }} e ^ {- w} \, dw + {\ sqrt {2}} \ int _ {0} ^ {\ infty} w ^ {- {\ tfrac {1} {2}}} e ^ {- w } \, dw \\ [5pt] = 2 {\ sqrt {2}} \ cdot \ Gamma \ left ({\ tfrac {3} {2}} \ right) + {\ sqrt {2}} \ cdot \ Гамма \ влево ({\ tfrac {1} {2}} \ right) \\ [5pt] = 2 {\ sqrt {2}} \ left ({\ tfrac {1} {2}} {\ sqrt {\ pi}} \ right) + {\ sqrt {2}} \ left ({\ sqrt {\ pi}} \ right) \\ [5pt] = 2 {\ sqrt {2 \ pi}}. \ end {align}}}{\ displaystyle {\ begin {align} \ int _ {0} ^ {\ infty} {\ frac {W (x)} {x {\ sqrt {x}}}} \, dx = \ int _ {0} ^ {\ infty} {\ frac {u} {ue ^ {u} {\ sqrt {ue ^ {u}}}}} (u + 1) e ^ {u} \, du \\ [5pt] = \ int _ {0} ^ {\ infty} {\ frac {u + 1} {\ sqrt {ue ^ {u}}} } ду \\ [5pt] = \ int _ {0} ^ {\ infty} {\ frac {u + 1} {\ sqrt {u}}} {\ frac {1} {\ sqrt {e ^ {u }}}} du \\ [5pt] = \ int _ {0} ^ {\ infty} u ^ {\ tfrac {1} {2}} e ^ {- {\ frac {u} {2}}} du + \ int _ {0} ^ {\ infty} u ^ {- {\ tfrac {1} {2}}} e ^ {- {\ frac {u} {2}}} du \\ [5pt] = 2 \ int _ {0} ^ {\ infty} (2w) ^ {\ tfrac {1} {2}} e ^ {- w} \, dw + 2 \ int _ {0} ^ {\ infty} ( 2w) ^ {- {\ tfrac {1} {2}}} e ^ {- w} \, dw \ quad (u = 2w) \\ [5pt] = 2 {\ sqrt {2}} \ int _ {0} ^ {\ infty} w ^ {\ tfrac {1} {2}} e ^ {- w} \, dw + {\ sqrt {2}} \ int _ {0} ^ {\ infty} w ^ {- {\ tfra c {1} {2}}} e ^ {- w} \, dw \\ [5pt] = 2 {\ sqrt {2}} \ cdot \ Gamma \ left ({\ tfrac {3} {2}} \ right) + {\ sqrt {2}} \ cdot \ Gamma \ left ({\ tfrac {1} {2}} \ right) \\ [5pt] = 2 {\ sqrt {2}} \ left ({ \ tfrac {1} {2}} {\ sqrt {\ pi}} \ right) + {\ sqrt {2}} \ left ({\ sqrt {\ pi}} \ right) \\ [5pt] = 2 {\ sqrt {2 \ pi}}. \ end {align}}}

Третья идентичность может быть получена из второго, если замена u = x, и первое также может быть получено из третьего пути замены z = 1 / √2 tan x.

За исключением z вдоль отрезка ветви (−∞, −1 / e] (где интеграл не сходится), главную ветвь функции Ламберта W можно вычислить с помощью следующего интеграла:

W (z) = z 2 π ∫ - π π (1 - ν детская кроватка ⁡ ν) 2 + ν 2 z + ν csc ⁡ ν e - ν детская кроватка ⁡ ν d ν = z π ∫ 0 π (1 - ν детская кроватка ⁡ ν) 2 + ν 2 z + ν csc ⁡ ν е - ν детская кроватка ⁡ ν d ν, {\ displaystyle {\ begin {align} W (z) = {\ frac {z} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} {\ frac {\ left (1- \ nu \ cot \ nu \ right) ^ {2} + \ nu ^ {2}} {z + \ nu \ csc \ nu e ^ { - \ nu \ cot \ nu}}} \, d \ nu \\ [5pt] = {\ frac {z} {\ pi}} \ int _ {0} ^ {\ pi} {\ frac {\ left (1- \ nu \ cot \ nu \ right) ^ {2} + \ nu ^ {2}} {z + \ nu \ csc \ nu e ^ {- \ nu \ cot \ nu}}} \, d \ nu, \ end {align}}}{\ displaystyle {\ begin {align} W (z) = {\ frac {z} {2 \ pi}} \ int _ {- \ pi} ^ {\ pi} {\ frac {\ left (1- \ nu \ cot \ nu \ right) ^ {2} + \ nu ^ {2}} {z + \ nu \ csc \ nu e ^ {- \ nu \ cot \ nu}}} \, d \ nu \\ [5pt] = {\ frac {z} {\ pi }} \ int _ {0} ^ {\ pi} {\ frac {\ left (1- \ nu \ cot \ nu \ right) ^ {2} + \ nu ^ {2}} {z + \ nu \ csc \ nu e ^ {- \ nu \ cot \ nu}}} \, d \ ню, \ конец {выровнено}}}

, где два интегральных выражения эквивалентны из-за симметрии подынтегрального выражения.

Неопределенные интегралы

∫ W (x) xdx = 1 2 (1 + W ( x)) 2 + C. ∫ W (A e B x) dx = 1 2 B (1 + W (A e B x)) 2 + C. {\ Displaystyle {\ begin {align} \ int {\ frac {W (x)} {x}} dx = {\ tfrac {1} {2}} {\ bigl (} 1 + W (x) {\ bigr)} ^ {2} + C. \ \ [5pt] \ int W \ left (Ae ^ {Bx} \ right) dx = {\ frac {1} {2B}} {\ bigl (} 1 + W \ left (Ae ^ {Bx} \ right) {\ bigr)} ^ {2} + C. \ end {align}}}{\ displaystyle {\ begin {align} \ int {\ frac {W (x)} {x}} dx = {\ tfrac {1} {2}} {\ bigl (} 1 + W (x) {\ bigr)} ^ {2} + C. \\ [5pt] \ int W \ left (Ae ^ {Bx} \ right) dx = {\ frac {1} {2B }} {\ bigl (} 1 + W \ left (Ae ^ {Bx} \ right) {\ bigr)} ^ {2} + C. \ end {выровнено }}}
Приложения

Решение уравнений

Функция Ламберта W используется для решений, в которых неизвестные количества встречаются как в основании, так и в экспоненте, или как внутри, так и вне логарифма. Стратегия состоит в том, чтобы преобразовать такое уравнение в одну из форм ze = w, а затем решить относительно z. используя функцию W.

Например, уравнение

3 x = 2 x + 2 {\ displaystyle 3 ^ {x} = 2x + 2}{\ displaystyle 3 ^ {x} = 2x + 2}

(где x - неизвестное действительное число) может быть решено следующим образом: переписав его как

(x + 1) 3 - x = 1 2 ⇔ (- x - 1) 3 - x - 1 = - 1 6 ⇔ (ln ⁡ 3) (- x - 1) e ( пер ⁡ 3) (- Икс - 1) = - пер ⁡ 3 6 {\ Displaystyle {\ begin {Выровнено} (x + 1) \ 3 ^ {- x} = {\ frac {1} {2}} \ \\ Leftrightarrow \ (- x-1) \ 3 ^ {- x-1} = - {\ frac {1} {6}} \\\ Leftrightarrow \ (\ ln 3) (- x-1) \ e ^ {(\ ln 3) (- x-1)} = - {\ frac {\ ln 3} {6}} \ end {align}}}{\ displaystyle {\ begin {align} ( x + 1) \ 3 ^ {- x} = {\ frac {1} {2}} \\\ Стрелка влево \ (-x-1) \ 3 ^ {- x-1} = - {\ frac {1 } {6}} \\\ Leftrightarrow \ (\ ln 3) (- x-1) \ e ^ {(\ ln 3) (- x-1)} = - {\ frac {\ ln 3} {6 }} \ end {align}}}

Это последнее уравнение имеет желаемую формулу, а решения для действительного x::

(ln 3) (- x - 1) = W 0 (- ln 3 6) или (ln ⁡ 3) (- x - 1) = W - 1 (- ln ⁡ 3 6) { \ Displaystyle (\ ln 3) (- x-1) = W_ {0} \ left ({\ frac {- \ ln 3} {6}} \ right) \ \ \ {\ textrm {или}} \ \ \ (\ ln 3) (- x-1) = W _ {- 1} \ left ({\ frac {- \ ln 3} {6}} \ right)}{\ Displaystyle (\ ln 3) (- x-1) = W_ {0} \ left ({\ frac {- \ ln 3} {6}} \ right) \ \ \ {\ textrm {или}} \ \ \ (\ ln 3) ( - x-1) = W _ {- 1} \ left ({\ frac {- \ ln 3} {6}} \ right)}

и, таким образом:

x = - 1 - W 0 (- ln ⁡ 3 6) ln ⁡ 3 = - 0,7901 1... или x = - 1 - W - 1 (- ln ⁡ 3 6) ln ⁡ 3 = 1,44456... {\ displaystyle x = -1- {\ frac {W_ {0} \ left (- {\ frac {\ ln 3} {6}} \ right)} {\ ln 3}} = - 0,79011... \ \ {\ textrm {или}} \ \ x = -1 - {\ frac { W _ {- 1} \ left (- {\ frac {\ ln 3} {6}} \ right)} {\ ln 3}} = 1.44456...}{\ displaystyle x = -1 - {\ frac {W_ {0} \ left (- {\ frac {\ ln 3} {6}} \ right)} {\ ln 3} } = - 0,79011... \ \ {\ textrm {или}} \ \ x = -1 - {\ frac {W _ {- 1} \ left (- {\ frac {\ ln 3} {6}} \ справа)} {\ ln 3}} = 1.44456...}

. Как правило, решение для

x = a + beansx {\ displaystyle x = a + b \, e ^ {cx}}{\ displaystyle x = a + b \, e ^ {cx}}

это:

x = a - 1 c W (- bceac) {\ displaystyle x = a - {\ frac { 1} {c}} W (-bc \, e ^ {ac})}{\ displaystyle x = a - {\ frac {1} {c}} W (-bc \, e ^ {ac})}

где a, b и c - комплексные константы, где b и c не равны нулю, а функция W имеет любой целочисленный порядок.

Вязкие потоки

Фронты и отложения зернистых и селевых потоков, а также фронты вязких флюидов в природных явлениях и в лабораторных экспериментах можно описать с помощью омега-функции Ламберта-Эйлера следующим образом:

ЧАС (Икс) знак равно 1 + W ((ЧАС (0) - 1) е (ЧАС (0) - 1) - х L), {\ Displaystyle Н (х) = 1 + W \ влево ((Н (0) -1) e ^ {(H (0) -1) - {\ frac {x} {L}}} \ right),}{\ displaystyle H (x) = 1 + W \ слева ((H (0) -1) e ^ {(H (0) -1) - {\ frac {x} {L}}} \ right),}

где H (x) - высота селевого потока, x - положение канала ниже по потоку, L - параметр единой модели, состоящий из нескольких физических и геометрических параметров потока, высоты потока и градиента гидравлического давления.

В потоке в трубе функция W Ламберта является частью явной формуки уравнения Коулбрука для нахождения коэффициента трения Дарси. Этот коэффициент используется для определения падения давления на прямом участке трубы, когда поток турбулентный.

Нейровизуализация

Функция Ламберта использовалась в области нейровизуализации для связи церебрального кровотока и потребления кислорода в вокселе мозга до соответствующего сигнала, зависящего от уровня оксигенации крови (жирный шрифт).

Химическая инженерия

Функция Ламберта использовалась в области химической инженерии для моделирования пористой толщины электродной пленки в суперконденсаторе на основе стеклоуглерода для электрохимического накопления энергии. Функция Ламберта W оказалась точным решением для процесса термической активации в газовой фазе, когда рост углеродной пленки и горения одной и той же пленки конкурируют друг с другом.

Материаловедение

Ламберт Функция использовалась в области роста эпитаксиальной пленки для определения критической дислокации начальной толщины пленки. Это расчетная толщина эпитаксиальной пленки, согласно термодинамическим принципам в пленке будут развиваться кристаллографические дислокации, чтобы минимизировать запасенную в пленках упругую энергию. Перед применением метода Ламберта W для решения этой задачи необходимо определить критическую толщину решения неявного уравнения. Ламберт W с легкостью превращает его в явное уравнение для аналитической обработки.

Пористая среда

Функция Ламберта W использовалась в области течения жидкости в пористой среде для моделирования наклона граница раздела двух гравитационно разделенных жидкостей в однородном наклонном пористом слое с постоянным падением и толщиной, где более тяжелая жидкость, закачиваемая в нижний конец, вытесняет более легкую жидкость, которая производит той же скоростью, из верхнего конца. Основная ветвь решения соответствует стабильному смещению, а ветвь −1 применяется, если смещение нестабильно, когда более тяжелая жидкость течет под более легкой.

Числа Бернулли и род Тодда

Уравнение (связано с производственными функциями чисел Бернулли и рода Тодда ):

Y = X 1 - e X {\ displaystyle Y = {\ frac {X} {1- e ^ {X}}}}Y = {\ frac {X} {1-e ^ {X}}}

можно решить с помощью двух ветвей W 0 и W −1:

X (Y) = {W - 1 (Y e Y) - W 0 (Y e Y) = Y - W 0 (Y e Y) для Y < − 1, W 0 ( Y e Y) − W − 1 ( Y e Y) = Y − W − 1 ( Y e Y) for − 1 < Y < 0. {\displaystyle X(Y)={\begin{cases}W_{-1}\left(Ye^{Y}\right)-W_{0}\left(Ye^{Y}\right)=Y-W_{0}\left(Ye^{Y}\right){\text{for }}Y<-1,\\W_{0}\left(Ye^{Y}\right)-W_{-1}\left(Ye^{Y}\right)=Y-W_{-1}\left(Ye^{Y}\right){\text{for }}-1{\ displaystyle X (Y) = {\ begin {cases} W _ {- 1} \ left (Ye ^ {Y} \ right) - W_ {0} \ left (Ye ^ {Y} \ right) = Y-W_ {0} \ left (Ye ^ {Y} \ right) {\ text {for}} Y <-1, \\ W_ { 0} \ left (Ye ^ {Y} \ right) -W _ {- 1} \ left (Ye ^ {Y} \ right) = YW _ {- 1} \ left (Ye ^ {Y} \ right) {\ text {for}} - 1 <Y <0. \ end {case}}}

Это приложение показывает, что разность ветвлений функции W может быть для решения других трансцендентных уравнений.

Статистика

Центроид гистограмм, определенно симметризованной дивергенции Кульбака - Лейблера (также называемой дивергенцией Джефф), имеет замкнутую форму с использованием функций Ламберта W.

Точные решения уравнения Уравнение Шредингера

Функция W Ламберта появляется в квантово-механическом потенциале, который дает пятое - после гармонического осциллятора плюс центробежный, кулоновского плюс обратный квадрат, Морса и - точное решение стационарного одномерного уравнения Шредингера в терминах конфлюэнтные гипергеометрические функции. Потенциал задается как

V = V 0 1 + W (e - x σ). {\ displaystyle V = {\ frac {V_ {0}} {1 + W \ left (e ^ {- {\ frac {x} {\ sigma}}}}} \ right)}}.}{\ displaystyle V = {\ frac {V_ {0}} {1 + W \ left (e ^ {- {\ frac {x} {\ sigma}}} \ right)}}.}

Особенность решения состоит в том, что каждый из двух фундаментальных решений, составляющих общее решение уравнения Шредингера, задается комбинацией двух конфлюэнтных гипергеометрических функций аргумента, пропорционального

z = W (e - x σ). {\ displaystyle z = W \ left (e ^ {- {\ frac {x} {\ sigma}}}} \ right).}{\ displaystyle z = W \ left (e ^ {- {\ frac {x} {\ sigma}}} \ right).}

Функция W Ламберта также представляет в точном решении для энергии связанного состояния одномерное уравнение Шредингера с двойным дельта-потенциалом.

Точные решения уравнений вакуума Эйнштейна

В решении уравнений вакуума Эйнштейна с метрикой Шварцшильда требуется функция W для перехода от координат Эддингтона - Финкельштейна к координатам Шварцшильда. По этой причине он также появляется при построении координат Крускала – Секереса.

Резонансы потенциала дельта-оболочки

s-волновые резонансы потенциала дельта-оболочки можно записать точно в терминах функции Ламберта W.

Термодинамическое равновесие

Если в реакции участвуют реагенты и продукты, имеющие теплоемкость, постоянную с температурой, то константа равновесия K подчиняется

ln ⁡ K = a T + b + c ln ⁡ T {\ displaystyle \ ln K = {\ frac {a} {T}} + b + c \ ln T}{\ displaystyle \ ln K = {\ frac {a} {T}} + б + с \ пер T}

для некоторых констант a, b, и c. Когда c (равно ΔC p / R) не равно нулю, мы можем найти значение или значения T, где K равно заданному значению, следующим образом, где мы используем L для ln T.

- a = (b - ln ⁡ K) T + c T ln ⁡ T = (b - ln ⁡ K) e L + c L e L - ac = (b - ln ⁡ K c + L) e L - aceb - ln ⁡ K c = (L + b - ln ⁡ K c) e L + b - ln ⁡ K c L = W (- aceb - ln ⁡ K c) + ln ⁡ K - bc T = exp ⁡ (W (- aceb - ln ⁡ K c) + ln ⁡ K - bc). {\ displaystyle {\ begin {align} -a = (b- \ ln K) T + cT \ ln T \\ = (b- \ ln K) e ^ {L} + cLe ^ {L} \\ [ 5pt] - {\ frac {a} {c}} = \ left ({\ frac {b- \ ln K} {c}} + L \ right) e ^ {L} \\ [5pt] - {\ frac {a} {c}} e ^ {\ frac {b- \ ln K} {c}} = \ left (L + {\ frac {b- \ ln K} {c}} \ right) e ^ { L + {\ frac {b- \ ln K} {c}}} \\ [5pt] L = W \ left (- {\ frac {a} {c}} e ^ {\ frac {b- \ ln K}) {c}} \ right) + {\ frac {\ ln Kb} {c}} \\ [5pt] T = \ exp \ left (W \ left (- {\ frac {a} {c}} e ^ { \ frac {b- \ ln K} {c}} \ right) + {\ frac {\ ln Kb} {c}} \ right). \ end {align}}}{\ displaystyle {\ begin {выровнено} -a = (b- \ ln K) T + cT \ ln T \\ = (b- \ ln К) е ^ {L} + cLe ^ {L} \\ [5pt] - {\ frac {a} {c}} = \ left ({\ frac {b- \ ln K} {c}} + L \ right) e ^ {L} \\ [5pt] - {\ frac {a} {c}} e ^ {\ frac {b- \ ln K} {c}} = \ left (L + {\ frac {b- \ ln K} {c} } \ right) e ^ {L + {\ frac {b - \ ln K} {c}}} \\ [5pt] L = W \ left (- {\ frac {a} {c}} e ^ { \ frac {b- \ ln K} {c}} \ right) + {\ frac {\ ln Kb} {c}} \\ [5pt] T = \ exp \ left (W \ left (- {\ frac {a} {c}} e ^ {\ frac {b- \ ln K} {c}} \ right) + {\ frac {\ ln Kb} {c}} \ right). \ end {align}}}

Если a и c имеют одинаковые знак будет либо два решения, либо ни одного (или одно, если аргумент W точно равен −1 / e). (Верхнее решение может не иметь отношения.) Если они имеют противоположные знаки, будет одно решение.

Соответствие AdS / CFT

Классические поправки конечного размера к дисперсионным соотношениям пиков и могут быть выражены в терминах функций Ламберта W.

Эпидемиология

В пределе t → ∞ модели SIR соотношение восприимчивых и выздоровевших людей имеет решение в терминах функций Ламберта W.

Определение времени полета снаряда

Общее время полета снаряда, которое испытывает сопротивление воздуха, в точной форме с использованием функций W Ламберта, пропорциональное его скорость , может быть определено.

Обобщения

Стандартная функция Ламберта W выражает точные решения трансцендентных алгебраических уравнений (по x) в форме:

e - cx = a 0 (x - r) {\ displaystyle e ^ { - cx} = a_ {0} (xr)}{\ displaystyle e ^ {- cx} = a_ {0} (xr)}

(1)

где a 0, c и r - действительные константы. Решение:

x = r + 1 c W (c e - c r a 0). {\ displaystyle x = r + {\ frac {1} {c}} W \ left ({\ frac {c \, e ^ {- cr}} {a_ {0}}} \ right).}{\ displaystyle x = r + {\ frac {1} {c}} W \ left ({\ frac {c \, e ^ {-cr}} {a_ {0}}} \ right).}

Обобщения функции Ламберта W включает:

e - cx = a 0 (x - r 1) (x - r 2), {\ displaystyle e ^ { - cx} = a_ {0} \ left (x-r_ {1} \ right) \ left (x-r_ {2} \ right),}{\ displaystyle e ^ {- cx} = a_ {0} \ left (x-r_ {1} \ right) \ left (x- r_ {2} \ справа),}

(2)

где r 1 и r 2 - действительные различные константы, корни квадратичного многочлена. Здесь решение - это функция с одним аргументом x, но такие термины, как r i и a 0, являются действующими функциями. В этом отношении обобщение напоминает функцию гипергеометрическую и функцию G Мейера, но принадлежит к другому классу функций. Когда r 1 = r 2, обе стороны (2) могут быть разложены на множители и уменьшены до (1), и, таким образом, решение сводится к стандартному решению W-функции.. Уравнение (2) выражает уравнение, определяющее поле дилатон, из которого выводится метрика R = T или линейной задачи двух тел гравитации в 1 + 1 измерениях (одно пространственное измерение и одно временное измерение измерение измерения для случая неравных масс покоя, а также собственных сил квантово-механической двухъямной модели дельта-функции Дирака для неравных зарядов в одном измерении.
  • Аналитический анализ решения собственных энергий частного случая квантово-механической задачи трех тел, а именно (трехмерной) молекулы-иона водорода. Здесь правая часть (1) заменена отношением многочленов бесконечного порядка по x:
e - cx = a 0 ∏ i = 1 ∞ (x - ri) ∏ i = 1 ∞ (x - си) {\ displaystyle e ^ {- cx} = a_ {0} {\ frac {\ displaystyle \ prod _ {i = 1} ^ {\ infty} (x-r_ {i})} {\ displaystyle \ prod _ {i = 1} ^ { \ infty} (x-s_ {i})}}}{\ displaystyle e ^ {- cx} = a_ {0} {\ frac {\ Displaystyle \ prod _ {я = 1} ^ {\ infty} (x-r_ {i})} {\ displaystyle \ prod _ {i = 1} ^ {\ infty} (x-s_ {i})} }}

(3)

где r i и s i - различные действительные константы, а x является функцией собственной энергии и межъядерного расстояния R. Уравнение (3) с его специальными случаями, выраженными в (1) и (2), относится к большому классу дифференциальные уравнения с запаздыванием. Г. Понятие «ложной производной» Х. Харди обеспечивает точные множественные корни для частных случаев (3).

Применения W-функции Ламберта в фундаментальных физических проблемах не исчерпывает даже для стандартного случая, выраженного в (1), как это было недавно замечено в области атомной, молекулярной и оптической физики.

Графики
Числовая оценка

Функция W может быть аппроксимирована использованием метода Ньютона с последовательными приближениями к w = W (z) (поэтому z = we)

wj + 1 = wj - wjewj - zewj + wjewj. {\ Displaystyle w_ { j + 1} = w_ {j} - {\ frac {w_ {j} e ^ {w_ {j}} - z} {e ^ {w_ {j}} + w_ {j} e ^ {w_ {j} }}}.}w_ {j + 1} = w_ {j} - {\ frac {w_ {j} e ^ {w_ { j}} - z} {e ^ {w_ {j}} + w_ {j} e ^ {w_ {j}}}}.

Функция W также может быть аппроксимирована использованием метода Галлея,

wj + 1 = wj - wjewj - zewj (wj + 1) - (wj + 2) (wjewj - z) 2 wj + 2 {\ displaystyle w_ {j + 1} = w_ {j} - {\ frac {w_ {j} e ^ {w_ {j}} - z} {e ^ {w_ {j}} \ left (w_ { j} +1 \ right) - {\ dfrac {\ left (w_ {j} +2 \ right) \ left (w_ {j} e ^ {w_ {j}} - z \ right)} {2w_ {j} + 2}}}}}{\ displaystyle w_ {j + 1} = w_ {j} - {\ frac {w_ {j} e ^ {w_ {j}} - z} {e ^ {w_ {j}} \ left (w_ {j} +1 \ right) - {\ dfrac {\ left (w_ {j} +2 \ right) \ left (w_ {j} e ^ {w_ {j}} - z \ right)} {2w_ {j} +2}}}}}

приведено у Corless et al. для вычислений W.

Программное обеспечение

Функция Lambert W реализована как LambertWв Maple, lambertwв GPglambertWв PARI ), lambertwв Matlab, также lambertwв Octave с пакетом specfun, как lambert_wв Maxima, как ProductLog(с скрытым псевдонимом LambertW) в Mathematica, как lambertwв пакете специальных функций Python scipy, как LambertWв модуле Perl ntheoryи как gsl_sf_lambert_W0, gsl_sf_lambert_Wm1функции в разделе специальные функции Научной библиотеки GNU (GSL). В библиотеках Boost C ++ используются следующие вызовы: lambert_w0, lambert_wm1, lambert_w0_primeи lambert_wm1_prime. В R функция Ламберта W реализована как функции lambertW0и lambertWm1в lamWпакете.

Код C ++ для всех ветвей сложной функции Ламберта W доступно на домашней странице Иштвана Мезо.

См. также
Примечания
Ссылки
Внешние ссылки
Викискладе есть медиафайлы, связанные с функцией Ламберта W.
Последняя правка сделана 2021-05-26 12:21:09
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте