Многозначная «функция»
График 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 - любые комплексные числа, то
выполняется тогда и только тогда, когда
При работе только с действительными числами две ветви W 0 и W −1 достаточно: для действительных чисел x и y уравнение
может быть решено относительно 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 Элементарные свойства, ветви и диапазон
- 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 году, в котором обсуждался особый случай нас.
Функция, которую рассматривал Ламберт, была
Эйлер преобразовал это уравнение в форму
Оба автора получили решение своих уравнений.
После того, как Эйлер решил это уравнение, он рассмотрел случай a = b. Взяв пределы, он вывел уравнение
Затем он положил 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.
Обратный
Области комплексной плоскости, для которых
. Более темные границы того же цвета включаются в более светлую область того же цвета. Точка в {-1,0} входит как в область
(синий), так и в область
(серая) область. Горизонтальные линии сетки кратны π.
График вышеупомянутой области в плоскости плоскости, где простая обратная 'верно. f = ze означает, что существует n такое, что , где n будет зависеть от значений z. Значение целого числа n резко изменится, когда ze станет означать на отрезке ветви , что будет означать, что ze ≤ 0, за исключением , где ze ≤ -1 / e.
Определите , где x и y действительны. Выражая e в полярных координатах, видно, что:
Для , отрезок ветви для будет неположительной действительной осью, так что:
и
Для , отрезок ветви для будет действительной осью с
Внутри границ, ограниченных указанным выше, не будет прерывистых изменений в , и эти области будут указывать, где функция W просто обратима: т.е. .
Исчисление
Производная
С помощью неявного дифференцирования можно показать, что все ветви W удовлетворяет дифференциальному уравнению
( W не дифференцируемо для z = −1 / e.) Как следствие, мы получаем формулу для производной W:
Используя тождество e = z / W (z), мы получаем эквивалентную формулу:
В начале координат мы имеем
Первообразная
Функция W (x) и многие выражения, включающие W (x), могут быть интегрированы с помощью подстановки w = W (x), т.е. x = мы:
(Последнее уравнение часто встречается в литературе, но не выполняется при x = 0). Одним из следствий этого (используя тот факт, что W 0 (e) = 1) является тождеством
Асимптотические разложения
Ряд Тейлора W 0 около 0 можно найти с помощью теоремы об обращении Лагранжа и определяется как
Радиус сходимости равенство 1 / e, что можно увидеть с помощью теста отношения . Функция, определяемая этим рядом, может быть расширена до голоморфной функции , определенные на всех комплексных числах с отрезком ветви вдоль интервала (−∞, −1 / e] ; эта голоморфная функция определяет главная ветвь функции Ламберта W.
Для больших значений x W 0 асимптотичен
где L 1 = ln x, L 2 = ln ln x, и [. l + 1 ]- неотрицательное число Стирлинга первого рода. Сохраняя только первые два членских разложения,
Другая действительная ветвь, W −1, определенная в интервале [−1 / e, 0), имеет приближение того же вида, когда x стремится к нулю, в этом случае L 1 = ln (−x) и L 2 = ln (−ln ( −x)).
Показано, что имеет место следующая оценка (верхняя оценка только для x ≥ e):
В 2013 году было доказано, что ветвь W −1 может быть ограниченным следующим образом:
Целочисленные и комплексные степени
Целочисленные вложения степени W 0 Разные степени W 0 ряд Тейлора (или Лорана ) в нуле:
В общем, для r ∈ ℤ Формула обращения Лагранжа дает
который, в общем, является рядо м Лорана порядка r. Эквивалентно, последний может быть записан в форме Taylor ex расширение степеней W 0 (x) / x:
это происходит для любого r ∈ ℂ и | х | < 1/e.
Identities
График W j (x e), где синий - для j = 0, а красный - для j = -1. Диагональная линия представляет интервалы, где W j (xe) = x
Из определения следует несколько тождеств:
Обратите внимание, что, поскольку 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) - другое решение.
Некоторые другие тождества:
- (который может быть расширен на другие n и x, если правильная ветка выбрано).
Подставляя −ln x в определение:
С повторяющейся экспонентой Эйлера h (x):
Специальные значения
Для любого ненулевого алгебраическое число x, W (x) трансцендентным номером. В самом деле, если W (x) равно нулю, то x также должно быть равен нулю, а если W (x) ненулевой и алгебраический, то по теореме Линдемана - Вейерштрасса, e должен быть трансцендентным, что означает, что x = W (x) также должен быть трансцендентным.
Следующие специальные значения главной ветви:
- (постоянная омега ).
Представления
Основныевь функц ии Лам берта может быть представлен интегралом благодаря Пуассону:
В более широкой области −1 / e ≤ x ≤ e значительно более простое представление находится с помощью Mező:
Следующая непрерывная дробь представление также справедливо для главной ветви:
Также, если | W (z) | < 1:
В свою очередь, если | W (z) |>е, то
.
Другие формулы
Определенные интегралы
Существует несколько полезных формул определенного интеграла, включающих основные ветвь функции W, в том числе следующие:
Первую идентичность можно найти, записав интеграл Гаусса в полярных координат.
Вторая идентичность можно получить, сделав замену u = W (x), которая дает
Таким образом,
Третья идентичность может быть получена из второго, если замена u = x, и первое также может быть получено из третьего пути замены z = 1 / √2 tan x.
За исключением z вдоль отрезка ветви (−∞, −1 / e] (где интеграл не сходится), главную ветвь функции Ламберта W можно вычислить с помощью следующего интеграла:
, где два интегральных выражения эквивалентны из-за симметрии подынтегрального выражения.
Неопределенные интегралы
Приложения
Решение уравнений
Функция Ламберта W используется для решений, в которых неизвестные количества встречаются как в основании, так и в экспоненте, или как внутри, так и вне логарифма. Стратегия состоит в том, чтобы преобразовать такое уравнение в одну из форм ze = w, а затем решить относительно z. используя функцию W.
Например, уравнение
(где x - неизвестное действительное число) может быть решено следующим образом: переписав его как
Это последнее уравнение имеет желаемую формулу, а решения для действительного x::
и, таким образом:
. Как правило, решение для
это:
где a, b и c - комплексные константы, где b и c не равны нулю, а функция W имеет любой целочисленный порядок.
Вязкие потоки
Фронты и отложения зернистых и селевых потоков, а также фронты вязких флюидов в природных явлениях и в лабораторных экспериментах можно описать с помощью омега-функции Ламберта-Эйлера следующим образом:
где H (x) - высота селевого потока, x - положение канала ниже по потоку, L - параметр единой модели, состоящий из нескольких физических и геометрических параметров потока, высоты потока и градиента гидравлического давления.
В потоке в трубе функция W Ламберта является частью явной формуки уравнения Коулбрука для нахождения коэффициента трения Дарси. Этот коэффициент используется для определения падения давления на прямом участке трубы, когда поток турбулентный.
Нейровизуализация
Функция Ламберта использовалась в области нейровизуализации для связи церебрального кровотока и потребления кислорода в вокселе мозга до соответствующего сигнала, зависящего от уровня оксигенации крови (жирный шрифт).
Химическая инженерия
Функция Ламберта использовалась в области химической инженерии для моделирования пористой толщины электродной пленки в суперконденсаторе на основе стеклоуглерода для электрохимического накопления энергии. Функция Ламберта W оказалась точным решением для процесса термической активации в газовой фазе, когда рост углеродной пленки и горения одной и той же пленки конкурируют друг с другом.
Материаловедение
Ламберт Функция использовалась в области роста эпитаксиальной пленки для определения критической дислокации начальной толщины пленки. Это расчетная толщина эпитаксиальной пленки, согласно термодинамическим принципам в пленке будут развиваться кристаллографические дислокации, чтобы минимизировать запасенную в пленках упругую энергию. Перед применением метода Ламберта W для решения этой задачи необходимо определить критическую толщину решения неявного уравнения. Ламберт W с легкостью превращает его в явное уравнение для аналитической обработки.
Пористая среда
Функция Ламберта W использовалась в области течения жидкости в пористой среде для моделирования наклона граница раздела двух гравитационно разделенных жидкостей в однородном наклонном пористом слое с постоянным падением и толщиной, где более тяжелая жидкость, закачиваемая в нижний конец, вытесняет более легкую жидкость, которая производит той же скоростью, из верхнего конца. Основная ветвь решения соответствует стабильному смещению, а ветвь −1 применяется, если смещение нестабильно, когда более тяжелая жидкость течет под более легкой.
Числа Бернулли и род Тодда
Уравнение (связано с производственными функциями чисел Бернулли и рода Тодда ):
можно решить с помощью двух ветвей W 0 и W −1:
Это приложение показывает, что разность ветвлений функции W может быть для решения других трансцендентных уравнений.
Статистика
Центроид гистограмм, определенно симметризованной дивергенции Кульбака - Лейблера (также называемой дивергенцией Джефф), имеет замкнутую форму с использованием функций Ламберта W.
Точные решения уравнения Уравнение Шредингера
Функция W Ламберта появляется в квантово-механическом потенциале, который дает пятое - после гармонического осциллятора плюс центробежный, кулоновского плюс обратный квадрат, Морса и - точное решение стационарного одномерного уравнения Шредингера в терминах конфлюэнтные гипергеометрические функции. Потенциал задается как
Особенность решения состоит в том, что каждый из двух фундаментальных решений, составляющих общее решение уравнения Шредингера, задается комбинацией двух конфлюэнтных гипергеометрических функций аргумента, пропорционального
Функция W Ламберта также представляет в точном решении для энергии связанного состояния одномерное уравнение Шредингера с двойным дельта-потенциалом.
Точные решения уравнений вакуума Эйнштейна
В решении уравнений вакуума Эйнштейна с метрикой Шварцшильда требуется функция W для перехода от координат Эддингтона - Финкельштейна к координатам Шварцшильда. По этой причине он также появляется при построении координат Крускала – Секереса.
Резонансы потенциала дельта-оболочки
s-волновые резонансы потенциала дельта-оболочки можно записать точно в терминах функции Ламберта W.
Термодинамическое равновесие
Если в реакции участвуют реагенты и продукты, имеющие теплоемкость, постоянную с температурой, то константа равновесия K подчиняется
для некоторых констант a, b, и c. Когда c (равно ΔC p / R) не равно нулю, мы можем найти значение или значения T, где K равно заданному значению, следующим образом, где мы используем L для ln T.
Если a и c имеют одинаковые знак будет либо два решения, либо ни одного (или одно, если аргумент W точно равен −1 / e). (Верхнее решение может не иметь отношения.) Если они имеют противоположные знаки, будет одно решение.
Соответствие AdS / CFT
Классические поправки конечного размера к дисперсионным соотношениям пиков и могут быть выражены в терминах функций Ламберта W.
Эпидемиология
В пределе t → ∞ модели SIR соотношение восприимчивых и выздоровевших людей имеет решение в терминах функций Ламберта W.
Определение времени полета снаряда
Общее время полета снаряда, которое испытывает сопротивление воздуха, в точной форме с использованием функций W Ламберта, пропорциональное его скорость , может быть определено.
Обобщения
Стандартная функция Ламберта W выражает точные решения трансцендентных алгебраических уравнений (по x) в форме:
| | (1) |
где a 0, c и r - действительные константы. Решение:
Обобщения функции Ламберта W включает:
- Приложение к общей теории относительности и квантовой механике (квантовой гравитации ) в более низких измерениях, фактически ссылка (неизвестно до 2007 г.) между этими двумя областями, где правая часть (1) заменена квадратичным многочленом от x:
| | (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:
| | (3) |
- где r i и s i - различные действительные константы, а x является функцией собственной энергии и межъядерного расстояния R. Уравнение (3) с его специальными случаями, выраженными в (1) и (2), относится к большому классу дифференциальные уравнения с запаздыванием. Г. Понятие «ложной производной» Х. Харди обеспечивает точные множественные корни для частных случаев (3).
Применения W-функции Ламберта в фундаментальных физических проблемах не исчерпывает даже для стандартного случая, выраженного в (1), как это было недавно замечено в области атомной, молекулярной и оптической физики.
Графики
Числовая оценка
Функция W может быть аппроксимирована использованием метода Ньютона с последовательными приближениями к w = W (z) (поэтому z = we)
Функция W также может быть аппроксимирована использованием метода Галлея,
приведено у Corless et al. для вычислений W.
Программное обеспечение
Функция Lambert W реализована как LambertWв Maple, lambertwв GP (и glambertWв 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 доступно на домашней странице Иштвана Мезо.
См. также
Примечания
Ссылки
- Corless, R.; Gonnet, G.; Заяц, Д.; Джеффри, Д.; Кнут, Дональд (1996). «О функции W Ламберта» (PDF). Успехи в вычислительной математике. 5 : 329–359. DOI : 10.1007 / BF02124750. ISSN 1019-7168. S2CID 29028411. Архивировано из оригинального (PDF) 14 декабря 2010 г. Проверено 10 марта 2007 г.
- Шапо-Блондо, Ф.; Монир, А. (2002). «Оценка W-функций Ламберта и ее применение для генерации обобщенного гауссовского шума с показателем 1/2» (PDF). IEEE Trans. Сигнальный процесс. 50 (9). DOI : 10.1109 / TSP.2002.801912. Архивировано из оригинального (PDF) 28 марта 2012 г. Проверено 10 марта 2004 г.
- Фрэнсис; и другие. (2000). «Количественная общая теория периодического дыхания». Тираж. 102 (18): 2214–21. CiteSeerX 10.1.1.505.7194. doi : 10.1161 / 01.cir.102.18.2214. PMID 11056095. S2CID 14410926.(функция Ламберта используется для решения дифференциальной динамики при заболеваниях человека.)
- Hayes, B. (2005). «Почему W?» (PDF). Американский ученый. 93 (2): 104–108. doi : 10.1511 / 2005.2.104.
- Рой, Р.; Олвер, Ф. В. Дж. (2010), «W-функция Ламберта», в Олвер, Фрэнк У. Дж. ; Лозье, Даниэль М.; Бойсверт, Рональд Ф.; Кларк, Чарльз У. (ред.), Справочник NIST по математическим функциям, Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248
- Стюарт, Шон М. (2005). «Новая элементарная функция для наших учебных программ?» (PDF). Австралийский старший математический журнал. 19 (2): 8–26. ISSN 0819-4564. ЭРИК EJ720055. Краткое содержание.
- Веберик, Д., «Развлечение с функцией Ламберта W (x)» arXiv: 1003.1628 (2010) ; Веберик, Д. (2012). «W-функция Ламберта для приложений в физике». Компьютерная физика. 183 (12): 2622–2628. arXiv : 1209.0735. Bibcode : 2012CoPhC.183.2622V. doi : 10.1016 / j.cpc.2012.07.008. S2CID 315088.
- Чатзигеоргиу, И. (2013). «Границы функции Ламберта и их применение для анализа сбоев взаимодействия пользователей». Письма связи IEEE. 17 (8): 1505–1508. arXiv : 1601.04895. DOI : 10.1109 / LCOMM.2013.070113.130972. S2CID 10062685.
Внешние ссылки
| Викискладе есть медиафайлы, связанные с функцией Ламберта W. |