람베르트 W 함수

람베르트 W함수에서 넘어옴

역사 raw
대문 랜덤 문서 최근 토론
소인수의 수를 세는 오메가 함수에 대한 내용은 소인수 계량 함수 문서를 참고해 주세요
1. 개요
1.1. 미적분1.2. 알려진 값
2. 활용
2.1. 방정식의 해 구하기 12.2. 방정식의 해 구하기 2
3. 관련 문서

1. 개요[편집]

람베르트 W\boldsymbol W 함수(Lambert W\boldsymbol W function)는 특수함수의 하나로, 오메가 함수(Omega function) 또는 곱 로그(Product logarithm)[1]라고도 한다.

함수의 정의에 앞서 우선 다음과 같은 지수함수를 정의해 보자.
y=xexy = xe^x
여기서 xx와 yy를 서로 바꾸어 역함수
x=yeyx = ye^y
를 얻는데, 이를 만족하는 yy를 y≡W(x)y \equiv W(x)로 정의하고, 람베르트 WW 함수라 한다.[2] 즉, W(x)eW(x)=xW(x)e^{W(x)}=x이다.

이 함수는 절대로 초등함수로 나타낼 수 없는데, 동일한 수로 곱과 지수가 엉켜 있는 형태이기 때문.

게다가 y=xexy = xe^x가 x=−1x = -1에서 극솟값 및 최솟값 −e−1-e^{-1}을 나타내므로, 람베르트 WW 함수는 기본적으로 음함수이다. 그래서 양함수로 나타내기 위해 y=−1y = -1을 기점으로 W−1(x)W_{-1}(x)[3]와 W0(x)W_0(x)[4]로 쪼개서 나타낸다. 즉 W(x)W(x)로 이르는 것은 실제로는
W−1(x)1R(W−1(x))∪W0(x)1R(W0(x))\dfrac{W_{-1}(x)}{{\bf1}_{\mathbb R}(W_{-1}(x))} \cup \dfrac{W_0(x)}{{\bf1}_{\mathbb R}(W_0(x))}[5]
인 셈이다.

아래는 람베르트 WW 함수의 그래프이다. 위의 설명과 같이 두 영역으로 나뉘어 나타난다.

파일:나무_람베르트W함수_그래프_NEW.png

W0(x)W_0(x)은 매클로린 전개를 이용하여 무한급수로 나타낼 수 있는데 다음과 같다.
W0(x)=∑n=1∞(−n)n−1n!xn=x−x2+32x3−83x4+12524x5−⋯\begin{aligned} W_{0}(x) &= \sum_{n=1}^\infty \frac{(-n)^{n-1}}{n!}x^n \\ &= x - x^2 + \frac32x^3 - \frac83x^4 + \frac{125}{24}x^5 - \cdots\end{aligned}

한편, xex=1xe^x = 1, 즉 W0(1)W_0(1)을 오메가 상수라고 하며 Ω\Omega로 나타낸다. 위 무한급수 식에 x=1x=1을 대입하면 얻을 수 있고, 구체적인 값은 약 0.567143이다.

일반화된 버전으로 Wn(x)W_n(x)가 있는데, nn이 −1-1, 00이 아닌 경우 무조건 복소수 값을 띤다. 심지어 nn이 −1-1, 00인 경우에도 상술한 범위[6]를 벗어나면 복소수가 된다.

1.1. 미적분[편집]

이 함수의 미분은 x=W(x)eW(x)x = W(x)e^{W(x)}를 이용해서 구할 수 있으며, 음함수의 미분법을 사용하여 양변을 미분하면
1=W′(x)eW(x)+W(x)eW(x)W′(x)=W′(x)eW(x)(W(x)+1)\begin{aligned} 1 &= W'(x)e^{W(x)} + W(x) e^{W(x)}W'(x) \\ &= W'(x)e^{W(x)}(W(x)+1) \end{aligned}
로 쓸 수 있다. 그런데 x=W(x)eW(x)x = W(x) e^{W(x)}으로부터
eW(x)=xW(x)e^{W(x)} = \dfrac x{W(x)}
이므로
1=W′(x)eW(x)(W(x)+1)=W′(x)xW(x)(W(x)+1)W(x)=xW′(x)(W(x)+1)\begin{aligned} 1 &= W'(x)e^{W(x)}(W(x)+1) \\ &= W'(x)\dfrac x{W(x)}(W(x)+1) \\ W(x) &= xW'(x)(W(x)+1) \end{aligned}
으로 쓸 수 있음에 따라
W′(x)=W(x)x(W(x)+1)W'(x)=\dfrac{W(x)}{x(W(x)+1)}
을 얻는다. 단, x=−e−1x=-e^{ -1 }에서는 미분 가능하지 않으며, x=0x=0에서는 치환 전의 식에 따라 W′(0)=1W'(0)=1을 얻을 수 있다.

부정적분을 구할 때는 부분적분법을 이용하면 된다.
∫W(x) dx=xW(x)−∫xW′(x) dx=xW(x)−∫W(x)eW(x)W′(x) dx=xW(x)−∫W(x)eW(x) dW(x)=xW(x)−eW(x)(W(x)−1)+C=xW(x)−xW(x)(W(x)−1)+C=x[W(x)−1+1W(x)]+C\begin{aligned} \int W(x)\,{\rm d}x &= xW(x) - \int xW'(x)\,{\rm d}x \\ &= xW(x) - \int W(x)e^{W(x)}W'(x)\,{\rm d}x \\ &= xW(x) - \int W(x)e^{W(x)}\,{\rm d}W(x) \\ &= xW(x) - e^{W(x)}(W(x)-1) + C \\ &= xW(x) - \dfrac x{W(x)}(W(x)-1) + C \\ &= x\left[ W(x)-1+\frac1{W(x)} \right]+C \end{aligned}
중간 과정에서 미분할 때 썼던 몇몇 공식을 이용했다.

1.2. 알려진 값[편집]

  • W(−π2)=iπ2W\left( -\dfrac\pi2 \right) = \dfrac{i\pi}2
  • W(−1e)=−1W\left(-\dfrac1e\right)=-1
  • W(−ln⁡aa)=−ln⁡a(e−1≤a≤e)W\left( -\dfrac{\ln a}a \right) = -\ln a \quad (e^{-1} \le a \le e)
  • W(0)=0W(0)=0
  • W(1)≡ΩW(1) \equiv \Omega: 이 값을 오메가 상수라 한다.
  • W(e)=1W(e)=1

2. 활용[편집]

2.1. 방정식의 해 구하기 1[편집]

이 문단에서는 방정식 [math(x^x=a)]의 해를 람베르트 WW 함수로 구해보도록 하자.

양변에 자연로그를 씌우면,
xln⁡x=ln⁡ax\ln x = \ln a
이다. x=eln⁡xx = e^{\ln x} 이므로 위 식은 eln⁡xln⁡x=ln⁡ae^{\ln x}\ln x = \ln a 로 변하는데, 이 때 람베르트 WW 함수를 양변에 취하면 정의에 따라
ln⁡x=W(ln⁡a)\ln x = W(\ln a)
이고, 식을 정리하면
x=eW(ln⁡a)x = e^{W(\ln a)}
를 얻는다.[7]

위 결과를 본 방정식에 대입하면
xx=(eW(ln⁡a))eW(ln⁡a)=eW(ln⁡a) eW(ln⁡a)=eln⁡a=a\begin{aligned} x^x &= (e^{W(\ln a)})^{e^{W(\ln a)}} \\ &= e^{W(\ln a)\,e^{W(\ln a)}} \\ &= e^{\ln a} \\ &= a \end{aligned}
가 되므로 답과 일치한다.

2.2. 방정식의 해 구하기 2[편집]

이 문단에서는 ax=bx+ca^x = bx+c의 해를 람베르트 WW 함수로 구해보자.

z=bx+cz = bx+c으로 놓으면,
x=z−cbx = \dfrac{z-c}b
가 되고, 대입하면 a(z−c)/b=za^{{(z-c)}/b} = z, 양 변에 ac/ba^{ c/b}를 곱하면
(a1/b)z=ac/bz(a^{1/b})^z = a^{c/b}z
가 된다. 여기서 상수를 각각 a1/b≡pa^{1/b} \equiv p, ac/b≡qa^{ c/b} \equiv q로 치환하자.

t=pz→z=log⁡ptt = p^z \to z = \log_pt로 놓으면,
t=qlog⁡ptt = q\log_pt
가 되고 이 식을 변형하면
pt/q=tp^{t/q} = t
가 된다. 이제 양변에 −1/t-1/t제곱을 취하면
p−1/q=(1t)1/tp^{-1/q} = \left( \dfrac1t \right)^{1/t}
가 된다. 이 때, 1/t≡u1/t \equiv u로 치환해주면,
uu=p−1/qu^u = p^{-1/q}
로, 바로 윗 문단에서 푼 xx=ax^x = a 꼴이 되었다. uu에 대해 풀고 치환 했던 문자들을 정리하면, xx에 대한 해는 다음과 같이 나타낼 수 있다.
x=−[1ln⁡aW(−a−c/bln⁡ab)+cb]x = - \left[ \dfrac1{\ln a}W\left(-\dfrac{a^{- c/b}\ln a}b \right)+\dfrac cb \right]

3. 관련 문서[편집]

[1] 정의가 로그함수(x=eyx = e^y)와 유사하다.[2] 예를 들어서 x=ex=e라고 한다면 e=yeye=ye^y를 만족하는 yy의 값은 y=1y=1이므로 y=W(e)=1y=W(e)=1이 된다.[3] 정의역: [−e−1, 0)[-e^{-1}, \,0)[4] 정의역: [−e−1, ∞)[-e^{-1},\, \infty)[5] 1R{\bf1}_{\mathbb R}은 실수 판별 함수이다. 실수이면 11, 실수가 아닌 경우 00이다. 따라서 함숫값이 실수가 아닐 경우 정의역에서 제외된다.[6] n=−1n=-1인 경우 −e−1≤x<0-e^{-1} \le x < 0, n=0n=0인 경우 x≥−e−1x \ge -e^{-1}[7] 이 때, eW(ln⁡a)e^{W(\ln a)}를 aa의 초제곱근(Super-root)이라고 한다.