Функция Райта Омега
wrightOmega( вычисляет функцию Райта omega x)x. Если z является матрицей, wrightOmega действует элементарно на z.
Вычислите функцию омега Райта для этих чисел. Поскольку эти числа не являются символическими объектами, получаются результаты с плавающей запятой:
wrightOmega(1/2)
ans =
0.7662wrightOmega(pi)
ans =
2.3061wrightOmega(-1+i*pi)ans = -1.0000 + 0.0000
Вычислите функцию Райта omega для чисел, преобразованных в символические объекты. Для большинства символических (точных) чисел, wrightOmega возвращает неразрешенные символьные вызовы:
wrightOmega(sym(1/2))
ans = wrightOmega(1/2)
wrightOmega(sym(pi))
ans = wrightOmega(pi)
Для некоторых точных чисел, wrightOmega имеет особые значения:
wrightOmega(-1+i*sym(pi))
ans =
-1Вычислите функцию Райта omega для x и sin(x) + x*exp(x). Для символьных переменных и выражений: wrightOmega возвращает неразрешенные символьные вызовы:
syms x wrightOmega(x) wrightOmega(sin(x) + x*exp(x))
ans = wrightOmega(x) ans = wrightOmega(sin(x) + x*exp(x))
Теперь вычислите производные этих выражений:
diff(wrightOmega(x), x, 2) diff(wrightOmega(sin(x) + x*exp(x)), x)
ans = wrightOmega(x)/(wrightOmega(x) + 1)^2 -... wrightOmega(x)^2/(wrightOmega(x) + 1)^3 ans = (wrightOmega(sin(x) + x*exp(x))*(cos(x) +... exp(x) + x*exp(x)))/(wrightOmega(sin(x) + x*exp(x)) + 1)
Вычисление омега-функции Райта для элементов матрицы M и вектор V:
M = [0 pi; 1/3 -pi]; V = sym([0; -1+i*pi]); wrightOmega(M) wrightOmega(V)
ans =
0.5671 2.3061
0.6959 0.0415
ans =
lambertw(0, 1)
-1[1] Корлесс, Р. М. и Д. Дж. Джеффри. «Функция Райта омега». Искусственный интеллект, автоматизированные рассуждения и символические вычисления (J. Calmet, B. Benhamou, O. Caprotti, L. Henocque, и V. Sorge, eds.). Берлин: Springer-Verlag, 2002, стр. 76-89.