6 svar
382 visningar
AlvinB behöver inte mer hjälp
AlvinB 4014
Postad: 2 aug 2019 13:58 Redigerad: 2 aug 2019 14:31

Integral av ln(ln(x))

Hej!

Jag sitter och klurar på integralen:

01lnlnx dx\displaystyle\int_0^1\ln\left(\ln\left(x\right)\right)\ dx

Jag är medveten om att ln(x)\ln(x) ger negativa värden på intervallet [0,1][0,1] och att den yttre logaritmen då kommer vara komplexvärd. Jag tar därför principalvärdet av logaritmen.

I försök att beräkna detta använde jag partialintegration för att få fram en primitiv funktion:

lnlnx dx=xlnlnx-1ln(x) dx=xlnlnx-lix+C\displaystyle\int\ln\left(\ln\left(x\right)\right)\ dx=x\ln\left(\ln\left(x\right)\right)-\int\frac{1}{\ln(x)}\ dx=x\ln\left(\ln\left(x\right)\right)-\text{li}\left(x\right)+C

där li(x)\text{li}(x) betecknar logaritmintegralen:

lix=0x1ln(t) dt\displaystyle\text{li}\left(x\right)=\int_0^x\frac{1}{\ln(t)}\ dt

(med Cauchy-principalvärdet av integralen på grund av diskontinuiteten vid t=1t=1)

Denna primitiva funktion ger:

01lnlnx dx=[xlnlnx-lix]01=limx1-(xlnlnx-lix)-limx0+(xlnlnx-lix)=\displaystyle\int_0^1\ln\left(\ln\left(x\right)\right)\ dx=[x\ln\left(\ln\left(x\right)\right)-\text{li}\left(x\right)]_0^1=\lim_{x\to1^-}(x\ln\left(\ln\left(x\right)\right)-\text{li}\left(x\right))-\lim_{x\to0^+}(x\ln\left(\ln\left(x\right)\right)-\text{li}\left(x\right))=

=limx1-(xlnlnx-lix)-0=limx1-xlnlnx-lix=\lim_{x\to1^-}(x\ln\left(\ln\left(x\right)\right)-\text{li}\left(x\right))-0=\lim_{x\to1^-}x\ln\left(\ln\left(x\right)\right)-\text{li}\left(x\right)

men detta gränsvärde lyckas jag inte beräkna, hur jag än bär mig åt. Någon som har några idéer?

tomast80 4245
Postad: 2 aug 2019 15:13 Redigerad: 2 aug 2019 15:17

Man borde kunna beräkna gränsvärdet med någon form av serieutveckling, men gissar att det blir rätt bökigt.

För övrigt noterade jag att den snarlika integralen:

01ln(ln(1x))dx=-γ\displaystyle \int_0^1 \ln (\ln (\frac{1}{x})) dx=-\gamma

Kan det vara till hjälp?

AlvinB 4014
Postad: 2 aug 2019 15:19 Redigerad: 2 aug 2019 15:21
tomast80 skrev:

Man borde kunna beräkna gränsvärdet med någon form av serieutveckling, men gissar att det blir rätt bökigt.

För övrigt noterade jag att det snarlika integralen:

01ln(ln(1x))dx=-γ\int_0^1 \ln (\ln (\frac{1}{x})) dx=-\gamma

Kan det vara till hjälp?

Jo, jag är hyfsat säker på att min integral har värdet:

01lnlnx dx=-γ+iπ\displaystyle\int_0^1\ln\left(\ln\left(x\right)\right)\ dx=-\gamma+i\pi

men jag är mer intresserad av hur man kommer fram till värdet än vad det faktiskt är. Har du någon idé om var man kan börja i de bökiga serieberäkningarna?

EDIT: Hur beräknar du förresten din reellvärda integral?

tomast80 4245
Postad: 2 aug 2019 16:50

Jag förstår. Har inte räknat ut den själv utan tagit sambandet för givet. Har tyvärr ingen konkret idé på hur man går vidare.

Hittade dock en hel del matnyttig info här:

http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.231.4513&rep=rep1&type=pdf

AlvinB 4014
Postad: 2 aug 2019 22:27 Redigerad: 3 aug 2019 00:40

Jag tror jag knäckte din integral nu. Med substitutionen t=-ln(x)t=-\ln(x) blir dx=-e-t dtdx=-e^{-t}\ dt och integralen fås till:

01lnln1x dx=01ln-lnx dx=-0e-tlnt dt=0e-tlnt dt\displaystyle\int_0^1\ln\left(\ln\left(\frac{1}{x}\right)\right)\ dx=\int_0^1\ln\left(-\ln\left(x\right)\right)\ dx=-\int_\infty^0 e^{-t}\ln\left(t\right)\ dt=\int_0^\infty e^{-t}\ln\left(t\right)\ dt

Denna integral kan sedan beräknas genom att utgå från gammafunktionen

Γz=0xz-1e-x dx\displaystyle\Gamma\left(z\right)=\int_0^\infty x^{z-1}e^{-x}\ dx

och derivera båda led med avseende på zz:

Γ'z=ddz(0xz-1e-x dx)\displaystyle\Gamma\prime\left(z\right)=\frac{d}{dz}(\int_0^\infty x^{z-1}e^{-x}\ dx)

Γ'z=0zxz-1e-x dx\displaystyle\Gamma\prime\left(z\right)=\int_0^\infty\frac{\partial}{\partial z}\left(x^{z-1}e^{-x}\right)\ dx

Γ'z=0xz-1lnxe-x dx\displaystyle\Gamma\prime\left(z\right)=\int_0^\infty x^{z-1}\ln\left(x\right)e^{-x}\ dx

Insättning av z=1z=1 ger därefter:

Γ'1=0x1-1lnxe-x dx=0e-xlnx dx\displaystyle\Gamma\prime\left(1\right)=\int_0^\infty x^{1-1}\ln\left(x\right)e^{-x}\ dx=\int_0^\infty e^{-x}\ln\left(x\right)\ dx

Då gammafunktionens derivata är lika med Γ'(z)=Γ(z)ψ(z)\Gamma\prime(z)=\Gamma(z)\psi(z) blir Γ'(1)=Γ(1)ψ(1)=ψ(1)=-γ\Gamma\prime(1)=\Gamma(1)\psi(1)=\psi(1)=-\gamma, och alltså får vi:

0e-xlnx dx=ψ1=-γ\displaystyle\int_0^\infty e^{-x}\ln\left(x\right)\ dx=\psi\left(1\right)=-\gamma

Förmodligen är det även möjligt att med hjälp av satsen om monoton konvergens bevisa att integralen har värdet -γ-\gamma.

EDIT: Lyckades även visa integralen utan hänvisning till digammafunktionen ψ(z)\psi(z):

Visa spoiler

Vi vill beräkna integralen:

I=0e-xlnx dx\displaystyle I=\int_0^\infty e^{-x}\ln\left(x\right)\ dx

Med hjälp av gränsvärdet:

limn(1-xn)n\lim_{n\to\infty}(1-\dfrac{x}{n})^n

kan vi skapa funktionsföljden:

fnx={(1-xn)n ifall xn0 annars        f_n\left(x\right)=\{\begin{matrix}(1-\frac{x}{n})^n\ \text{ifall}\ x\leq n\\0\ \text{annars}\ \ \ \ \ \ \ \ \end{matrix}

och vi ser att fn(x)exf_n(x)\nearrow e^x. Satsen om monoton konvergens ger då:

0e-xlnx dx=limn0fnxlnx dx=\displaystyle\int_0^\infty e^{-x}\ln\left(x\right)\ dx=\lim_{n\to\infty}\int_0^\infty f_n\left(x\right)\ln\left(x\right)\ dx=

fn(x)=0f_n(x)=0 för alla x>nx>n kan vi skriva integralen som:

=limn0nfnxlnx dx=limn0n1-xnnlnx dx=\displaystyle=\lim_{n\to\infty}\int_0^n f_n\left(x\right)\ln\left(x\right)\ dx=\lim_{n\to\infty}\int_0^n\left(1-\frac{x}{n}\right)^n\ln\left(x\right)\ dx=

Med substitutionen x=nux=nu, u=x/nu=x/n, dx=n dudx=n\ du erhålls:

=limnn011-unln(nu du=limnn011-unlnn+lnu du=\displaystyle=\lim_{n\to\infty}n\int_0^1\left(1-u\right)^n\ln(\left(nu\right)\ du=\lim_{n\to\infty}n\int_0^1\left(1-u\right)^n\left(\ln\left(n\right)+\ln\left(u\right)\right)\ du=

Substitutionen t=1-ut=1-u, u=1-tu=1-t, du=-dtdu=-dt leder nu till:

=limn-n10tnlnn+ln1-t dt=limnn(01tnlnn dt+01tn1-t dt)=\displaystyle=\lim_{n\to\infty}-n\int_1^0 t^n\left(\ln\left(n\right)+\ln\left(1-t\right)\right)\ dt=\lim_{n\to\infty}n(\int_0^1 t^n\ln\left(n\right)\ dt+\int_0^1 t^n\left(1-t\right)\ dt)=

Maclaurinutvecklingen av ln(1+x)\ln(1+x) lyder:

ln1+x=n=1-1n+1·xnn\displaystyle\ln\left(1+x\right)=\sum_{n=1}^\infty \left(-1\right)^{n+1}\cdot\frac{x^n}{n}

Byts xx mot -x-x erhålls:

ln1-x=n=1-1n+1·-xnn=-n=1xnn\displaystyle\ln\left(1-x\right)=\sum_{n=1}^\infty \left(-1\right)^{n+1}\cdot\frac{\left(-x\right)^n}{n}=-\sum_{n=1}^\infty\frac{x^n}{n}

Sätts detta nu in i den högra integralen fås:

I=limnnlnn01tn dt-n01tnk=1tkk dt=limnnlnn[tn+1n+1]01-nk=101tn+kk dt=\displaystyle I=\lim_{n\to\infty}n\ln\left(n\right)\int_0^1 t^n\ dt-n\int_0^1 t^n\sum_{k=1}^\infty\frac{t^k}{k}\ dt=\lim_{n\to\infty} n\ln\left(n\right)[\frac{t^{n+1}}{n+1}]_0^1-n\sum_{k=1}^\infty\int_0^1\frac{t^{n+k}}{k}\ dt=

=limnnn+1·lnn-nk=1[tn+k+1k(n+k+1)]01=limnnn+1·lnn-nk=11k(n+k+1)\displaystyle=\lim_{n\to\infty}\frac{n}{n+1}\cdot\ln\left(n\right)-n\sum_{k=1}^\infty[\frac{t^{n+k+1}}{k(n+k+1)}]_0^1=\lim_{n\to\infty}\frac{n}{n+1}\cdot\ln\left(n\right)-n\sum_{k=1}^\infty\frac{1}{k(n+k+1)}

Partialbråksuppdelning ger:

1k(n+k+1)=Ak+Bn+k+1\dfrac{1}{k(n+k+1)}=\dfrac{A}{k}+\dfrac{B}{n+k+1}

Med Heavisides handpåläggningsmetod fås A=1/(n+1)A=1/(n+1) och B=-1/(n+1)B=-1/(n+1). Detta ger:

I=limnnn+1·lnn-nn+1k=11k-1n+k+1=\displaystyle I=\lim_{n\to\infty}\frac{n}{n+1}\cdot\ln\left(n\right)-\frac{n}{n+1}\sum_{k=1}^\infty\frac{1}{k}-\frac{1}{n+k+1}=

Analyserar man nu summan som bildas till höger inser man att de vänstra termerna bildar summan 1+1/2+1/3+...1+1/2+1/3+... och de högra termerna -1/(n+2)-1/(n+3)-1/(n+4)-...-1/(n+2)-1/(n+3)-1/(n+4)-.... Således kommer alla termer efter k=n+1k=n+1 att ta ut tidigare termer, vilket gör att summan i själva verket bara är lika med 1+2+3+...1/n+1/(n+1)1+2+3+...1/n+1/(n+1):

limnnn+1(lnn-k=1n+11k)=limnnn+1(lnn-k=1n1k-1n+1\displaystyle\lim_{n\to\infty}\frac{n}{n+1}(\ln\left(n\right)-\sum_{k=1}^{n+1}\frac{1}{k})=\lim_{n\to\infty}\frac{n}{n+1}(\ln\left(n\right)-\sum_{k=1}^n\frac{1}{k}-\frac{1}{n+1}

Beräknas nu gränsvärdet erhålls slutligen:

I=limnnn+1limn(lnn-k=1n1k-1n+1)=1·(limn-(k=1n1k-lnn)-limn1n+1)=-γ-0=-γ\displaystyle I=\lim_{n\to\infty}\left(\frac{n}{n+1}\right)\lim_{n\to\infty}(\ln\left(n\right)-\sum_{k=1}^n\frac{1}{k}-\frac{1}{n+1})=1\cdot(\lim_{n\to\infty}-(\sum{k=1}^n\frac{1}{k}-\ln\left(n\right))-\lim_{n\to\infty}\frac{1}{n+1})=-\gamma-0=-\gamma

V.S.V.

AlvinB 4014
Postad: 2 aug 2019 22:49 Redigerad: 2 aug 2019 23:49

Nu tror jag även att bitarna föll på plats för min egen integral. Samma substitution, t=-ln(x)t=-\ln(x), ger:

01lnlnx dx=-0e-tln-t dt=0e-tln-t dt=\displaystyle\int_0^1\ln\left(\ln\left(x\right)\right)\ dx=-\int_\infty^0 e^{-t}\ln\left(-t\right)\ dt=\int_0^\infty e^{-t}\ln\left(-t\right)\ dt=

Då vi tar principalvärdet av logaritmen blir ln(-t)=ln(t)+iπ\ln(-t)=\ln(t)+i\pi, och vips så har vi trollat fram e-xln(x)e^{-x}\ln(x)-integralen från det förra inlägget!

=0e-tlnt+iπ dt=0e-tlnt dt+0e-t·iπ dt=-γ+iπ[-e-t]0=-γ+iπ0--1=-γ+iπ\displaystyle=\int_0^\infty e^{-t}\left(\ln\left(t\right)+i\pi\right)\ dt=\int_0^\infty e^{-t}\ln\left(t\right)\ dt+\int_0^\infty e^{-t}\cdot i\pi\ dt=-\gamma+i\pi[-e^{-t}]_0^\infty=-\gamma+i\pi\left(0-\left(-1\right)\right)=-\gamma+i\pi

Det verkar som tipset med den andra integralen gjorde susen. Jag får tacka så mycket. :-)

tomast80 4245
Postad: 3 aug 2019 05:49

Mycket vackert, AlvinB!

Kul att kunna bidra!

Svara
Close