Mathematics · Boek 5 · Bachelor Year 3

Universitaire wiskunde — Bachelor jaar 3

Universitaire wiskunde — Bachelor jaar 3 · Bachelor Year 3

10De Lebesgue-integraal

De integraal van Riemann snijdt het domein in kleine intervallen; die van Lebesgue snijdt het beeld: om ff te integreren, meet je de verzamelingen {f>t}\{f > t\}. De wijziging oogt onschuldig en is revolutionair. Limieten en integralen, in de theorie van Riemann eeuwig kibbelend (uniforme convergentie vereist!), worden verzoend door drie convergentiestellingen — monotone convergentie, Fatou, gedomineerde convergentie — waarvan de hypothesen bijna gênant zwak zijn. Dit hoofdstuk bouwt de integraal op een willekeurige maatruimte (X,A,μ)(X, \mathcal A, \mu), bewijst de drie stellingen, beslecht de precieze verhouding met de integraal van Riemann (een begrensde functie is riemann-integreerbaar dan en slechts dan als ze bijna overal continu is) en industrialiseert het differentiëren van integralen met een parameter — de techniek waarmee de weekendopgave 0sinxx ⁣dx\int_0^\infty\frac{\sin x}x\,\dd x en Rex2 ⁣dx\int_\R \eu^{-x^2}\dd x berekent.

10.1 Meetbare functies

Definitie 10.1

Zij (X,A)(X, \mathcal A) en (Y,B)(Y, \mathcal B) meetbare ruimten. Dan heet f ⁣:XYf \colon X \to Y meetbaar als f1(B)Af^{-1}(B) \in \mathcal A voor elke BBB \in \mathcal B. Voor reële (of [,+][-\infty,+\infty]-waardige) functies draagt Y=RY = \R haar Borel-σ\sigma-algebra, en volstaat het f1((t,+))={f>t}Af^{-1}(\intoo t{+\infty}) = \{f > t\} \in \mathcal A voor alle tRt \in \R na te gaan: de goede verzamelingen {B:f1(B)A}\{B : f^{-1}(B) \in \mathcal A\} vormen een σ\sigma-algebra (originelen verdragen zich met verzamelingsbewerkingen) die de voortbrengende halfrechten bevat (Definitie 9.2, Methode 9.17).

Propositie 10.2

(a) Samenstellingen van meetbare afbeeldingen zijn meetbaar; continue afbeeldingen zijn borel-meetbaar. (b) Zijn f,g ⁣:XRf, g \colon X \to \R meetbaar, dan ook f+gf + g, fgfg, max(f,g)\max(f,g), f\abs f en λf\lambda f. (c) Zijn de (fn)(f_n) meetbaar met waarden in [,+][-\infty, +\infty], dan zijn supnfn\sup_nf_n, infnfn\inf_nf_n, lim supfn\limsup f_n en lim inffn\liminf f_n meetbaar; en convergeert fnff_n \to f puntsgewijs, dan is ff meetbaar.

Bewijs. (a) (gf)1(B)=f1(g1(B))(g\circ f)^{-1}(B) = f^{-1}(g^{-1}(B)); en de continuïteit geeft de meetbaarheid via de voortbrengende open verzamelingen (Probleem 9.1, vraag 10, in algemene vorm). (b) (f,g) ⁣:XR2(f, g) \colon X \to \R^2 is meetbaar voor de Borel-σ\sigma-algebra van R2\R^2 — ga het na op de open blokken, die voortbrengen (de open verzamelingen van R2\R^2 zijn aftelbare verenigingen van rationale blokken): (f,g)1(U×V)=f1(U)g1(V)(f,g)^{-1}(U\times V) = f^{-1}(U)\cap g^{-1}(V) — en ++, ×\times en max\max zijn continu van R2\R^2 naar R\R: stel samen. (c) {supfn>t}=n{fn>t}\{\sup f_n > t\} = \bigcup_n\{f_n > t\}; inf=sup()\inf = -\sup(-); lim sup=infNsupnN\limsup = \inf_N\sup_{n \geq N}; en een puntsgewijze limiet is haar eigen lim sup\limsup.

Definitie 10.3

Een elementaire functie is een meetbare functie met eindig veel waarden: s=i=1nci1Ais = \sum_{i=1}^n c_i\,\mathbf 1_{A_i} met AiAA_i \in \mathcal A disjunct en ci0c_i \geq 0 (voor de niet-negatieve theorie). Haar integraal is

s ⁣dμ=iciμ(Ai)[0,+]\int s\,\dd\mu = \sum_i c_i\,\mu(A_i) \in [0, +\infty]

(met de afspraak 0=00\cdot\infty = 0); de waarde hangt niet van de voorstelling af (verfijn twee partities).

Stelling 10.4 (Benadering door elementaire functies)

Elke meetbare f ⁣:X[0,+]f \colon X \to [0, +\infty] is de puntsgewijze limiet van een stijgende rij elementaire functies:

sn=k=1n2nk12n1{k12nf<k2n}+n1{fn}f.s_n = \sum_{k=1}^{n2^n} \frac{k-1}{2^n}\, \mathbf 1_{\{\frac{k-1}{2^n} \leq f < \frac k{2^n}\}} + n\,\mathbf 1_{\{f \geq n\}} \nearrow f .

Bewijs. Elke sns_n is elementair (de verzamelingen zijn originelen van borelverzamelingen). Monotonie: van nn naar n+1n+1 gaan splitst elk dyadisch niveau in tweeën en verlaagt de toegekende waarde nooit (een punt met k12nf(x)<k2n\frac{k-1}{2^n} \leq f(x) < \frac k{2^n} krijgt 2k22n+1\frac{2k-2}{2^{n+1}} of 2k12n+1\frac{2k-1}{2^{n+1}}, allebei k12n\geq \frac{k-1}{2^n}; en de afkapwaarde nn stijgt mee). Convergentie: is f(x)<f(x) < \infty, dan is voor n>f(x)n > f(x) het verschil f(x)sn(x)2nf(x) - s_n(x) \leq 2^{-n}; en is f(x)=f(x) = \infty, dan is sn(x)=ns_n(x) = n \to \infty.

10.2 De integraal en de convergentiestellingen

Definitie 10.5

Voor meetbare f0f \geq 0:

f ⁣dμ=sup{s ⁣dμ:s elementair, 0sf}[0,+].\int f \,\dd\mu = \sup\Bigl\{\int s\,\dd\mu : s \text{ elementair},\ 0 \leq s \leq f\Bigr\} \in [0, +\infty].

Ze is per constructie monotoon in ff, en zet het elementaire geval voort (voor elementaire ff wordt het supremum in ff bereikt: vergelijk elementaire integralen via gemeenschappelijke verfijningen).

Stelling 10.6 (Monotone convergentie, Beppo Levi)

Is 0fnf0 \leq f_n \nearrow f puntsgewijs (meetbaar), dan is

fn ⁣dμf ⁣dμ.\int f_n\,\dd\mu \nearrow \int f\,\dd\mu .

Bewijs. ff is meetbaar (Propositie 10.2(c)) en fn\int f_n stijgt naar een zekere LfL \leq \int f (monotonie). Leg omgekeerd een elementaire s=ci1Aifs = \sum c_i\mathbf 1_{A_i} \leq f vast en een θ(0,1)\theta \in (0,1); de verzamelingen En={fnθs}E_n = \{f_n \geq \theta s\} zijn meetbaar en stijgen naar XX (waar s(x)>0s(x) > 0: f(x)s(x)>θs(x)f(x) \geq s(x) > \theta s(x), dus vanaf zeker moment fn(x)θs(x)f_n(x) \geq \theta s(x); en waar s(x)=0s(x) = 0: triviaal). Dan is

fnEnθs ⁣dμ=θiciμ(AiEn)nθiciμ(Ai)=θs\int f_n \geq \int_{E_n}\theta s\,\dd\mu = \theta\sum_i c_i\,\mu(A_i \cap E_n) \xrightarrow[n\to\infty]{} \theta\sum_ic_i\,\mu(A_i) = \theta\int s

wegens de continuïteit van onderen (Propositie 9.6(c)). Dus LθsL \geq \theta\int s voor alle θ<1\theta < 1 en alle elementaire sfs \leq f: LfL \geq \int f.

Gevolg 10.7

Voor meetbare f,g0f, g \geq 0 en c0c \geq 0 geldt (f+g)=f+g\int(f + g) = \int f + \int g en cf=cf\int cf = c\int f; en voor een reeks niet-negatieve meetbare functies is nfn=nfn\int\sum_nf_n = \sum_n\int f_n.

Bewijs. Voor elementaire functies is de additiviteit een berekening op een gemeenschappelijke verfijning. In het algemeen neem je snfs_n \nearrow f en tngt_n \nearrow g (Stelling 10.4): dan is sn+tnf+gs_n + t_n \nearrow f + g, en de monotone convergentie draagt de additiviteit over op de limiet. De uitspraak over de reeks is de monotone convergentie toegepast op de partiële sommen.

Stelling 10.8 (Lemma van Fatou)

Voor meetbare fn0f_n \geq 0 geldt

lim infnfn ⁣dμ    lim infnfn ⁣dμ.\int \liminf_n f_n \,\dd\mu \;\leq\; \liminf_n \int f_n\,\dd\mu .

Bewijs. Zij gN=infnNfng_N = \inf_{n\geq N}f_n: meetbaar, met 0gNlim inffn0 \leq g_N \nearrow \liminf f_n en gNfng_N \leq f_n voor elke nNn \geq N, zodat gNinfnNfn\int g_N \leq \inf_{n \geq N}\int f_n. Pas de monotone convergentie toe op het linkerlid: lim inffn=limNgNlimNinfnNfn=lim inffn\int\liminf f_n = \lim_N\int g_N \leq \lim_N\inf_{n\geq N}\int f_n = \liminf\int f_n.

Definitie 10.9

Een meetbare f ⁣:XRf \colon X \to \R (of C\C) heet integreerbaar als f ⁣dμ<\int\abs f\,\dd\mu < \infty; dan is f=f+f\int f = \int f^+ - \int f^- (positief en negatief deel; reëel en imaginair deel in het complexe geval). De integraal is lineair op de integreerbare functies (ontbind en zet weer samen op de positieve delen; het complexe geval herleidt tot het reële) en voldoet aan ff\abs{\int f} \leq \int\abs f (reëel geval: ±f=(±f)f\pm\int f = \int(\pm f) \leq \int\abs f; complex geval: vermenigvuldig met een constante van modulus 11 om de integraal reëel te maken). Een eigenschap geldt bijna overal als ze alleen op een μ\mu-nulverzameling faalt; en ff op een nulverzameling wijzigen verandert geen enkele integraal (het verschil wordt gedomineerd door 1N\infty\cdot\mathbf 1_N, met integraal 00).

Stelling 10.10 (Gedomineerde convergentie)

Zij fnff_n \to f bijna overal, met fng\abs{f_n} \leq g bijna overal voor een vaste integreerbare gg. Dan is ff integreerbaar en

fn ⁣dμf ⁣dμ,en zelfsfnf ⁣dμ0.\int f_n\,\dd\mu \longrightarrow \int f\,\dd\mu, \qquad\text{en zelfs}\quad \int\abs{f_n - f}\,\dd\mu \to 0 .

Bewijs. Laat een nulverzameling weg om de hypothesen puntsgewijs te maken. Uit fg\abs f \leq g volgt dat ff integreerbaar is. De functies hn=2gfnf0h_n = 2g - \abs{f_n - f} \geq 0 voldoen aan lim infhn=2g\liminf h_n = 2g; Fatou geeft

2glim inf(2gfnf)=2glim supfnf,\int 2g \leq \liminf\int\bigl(2g - \abs{f_n - f}\bigr) = \int 2g - \limsup\int\abs{f_n - f},

dus lim supfnf0\limsup\int\abs{f_n - f} \leq 0 (het aftrekken is geoorloofd: 2g<\int 2g < \infty). Ten slotte is fnffnf0\abs{\int f_n - \int f} \leq \int\abs{f_n - f} \to 0.

Methode 10.11

Sta je voor limnfn\lim_n\int f_n, probeer dan op volgorde: (1) is de rij monotoon (of is het een reeks met niet-negatieve termen)? Dan de monotone convergentie, zonder enige integreerbaarheidseis. (2) Is er één integreerbare dominant gfng \geq \abs{f_n}, gevonden met grove afschattingen (neem het “supn\sup_n” van de schattingen)? Dan de gedomineerde convergentie. (3) Geen dominant en geen monotonie? Dan begrenst Fatou nog steeds één kant, en kan de gelijkheid werkelijk falen: de ontsnappende bult fn=n1(0,1/n)f_n = n\mathbf 1_{\intoo0{1/n}} heeft fn=1\int f_n = 1 terwijl fn0f_n \to 0 bijna overal. Dominatie is precies wat verbiedt dat massa naar oneindig ontsnapt, verticaal of horizontaal.

10.3 Riemann tegenover Lebesgue

Stelling 10.12 (Criterium van Lebesgue)

Zij f ⁣:[a,b]Rf \colon \intcc ab \to \R begrensd. Dan is ff riemann-integreerbaar dan en slechts dan als ff λ\lambda-bijna overal continu is; in dat geval is ff lebesgue-integreerbaar en vallen de twee integralen samen.

Bewijs. Zij voor een verdeling σ=(a=x0<<xN=b)\sigma = (a = x_0 < \dots < x_N = b) de trapfuncties LσL_\sigma en UσU_\sigma op elke (xi1,xi)\intoo{x_{i-1}}{x_i} gelijk aan mi=inf[xi1,xi]fm_i = \inf_{[x_{i-1}, x_i]}f respectievelijk Mi=supM_i = \sup; de darbouxsommen zijn hun integralen (Riemann en Lebesgue vallen op trapfuncties samen en geven beide miΔxi\sum m_i\Delta x_i). Neem een rij verdelingen σn\sigma_n, elk een verfijning van de vorige, met maaswijdte 0\to 0 en met darbouxsommen die naar de onderste en bovenste darbouxintegraal van ff convergeren. Door de verfijningen is LσnL_{\sigma_n} niet-dalend en UσnU_{\sigma_n} niet-stijgend, puntsgewijs buiten de aftelbare verzameling DD van alle deelpunten; noem de limieten \ell en uu (meetbaar, Propositie 10.2). Voor xDx \notin D, met In(x)I_n(x) het open σn\sigma_n-interval dat xx bevat: (x)=supninfIn(x)f\ell(x) = \sup_n\inf_{I_n(x)}f en u(x)=infnsupIn(x)fu(x) = \inf_n\sup_{I_n(x)}f; en omdat de maaswijdten naar 00 krimpen, zijn dat de onderste en bovenste omhullenden van ff in xxu(x)(x)u(x) - \ell(x) is de oscillatie van ff in xx — zodat (x)=u(x)\ell(x) = u(x) precies wanneer ff in xx continu is. Met de monotone en de gedomineerde convergentie (begrensd, eindig interval):

[a,b] ⁣dλ=limnLσn=f,[a,b]u ⁣dλ=f.\int_{\intcc ab}\ell\,\dd\lambda = \lim_n\int L_{\sigma_n} = \underline{\int}f, \qquad \int_{\intcc ab}u\,\dd\lambda = \overline{\int}f .

Dus: ff riemann-integreerbaar     \iff f=f\underline\int f = \overline\int f     \iff (u)=0\int(u - \ell) = 0     \iff u=u = \ell bijna overal (u0u - \ell \geq 0; Oefening 10.5)     \iff ff bijna overal continu. In dat geval is fu\ell \leq f \leq u met =u\ell = u bijna overal: ff is bijna overal gelijk aan de meetbare \ell en dus lebesgue-meetbaar (volledigheid van λ\lambda), met f ⁣dλ= ⁣dλ=f=abf\int f\,\dd\lambda = \int\ell\,\dd\lambda = \underline\int f = \int_a^bf.

Voorbeeld 10.13

1Q\mathbf 1_\Q is nergens continu: niet riemann-integreerbaar — maar voor Lebesgue triviaal: 1Q ⁣dλ=λ(Q)=0\int\mathbf 1_\Q\,\dd\lambda = \lambda(\Q) = 0. De functie van Thomae (1q\frac1q in de rationale getallen pq\frac pq en 00 elders) is precies in de irrationale punten continu: riemann-integreerbaar met integraal 00. En oneigenlijke riemannintegralen zijn een ander begrip: 0sinxx ⁣dx\int_0^\infty\frac{\sin x}x\,\dd x convergeert als limiet van 0A\int_0^A (de weekendopgave berekent hem als π2\frac\pi2), maar sinxxL1((0,+))\frac{\sin x}x \notin L^1(\intoo0{+\infty}): de absolute integraal divergeert als de harmonische reeks (Oefening 10.6). De theorie van Lebesgue ruilt voorwaardelijke convergentie in voor robuuste limietstellingen.

10.4 Integralen met parameters

Overal is (X,A,μ)(X, \mathcal A, \mu) een maatruimte, TT een metrische ruimte (de parameter) en f ⁣:T×XCf \colon T \times X \to \C met f(t,)f(t, \cdot) integreerbaar voor elke tt; stel F(t)=Xf(t,x) ⁣dμ(x)F(t) = \int_X f(t, x)\,\dd\mu(x).

Stelling 10.14 (Continuïteit)

Stel dat tf(t,x)t \mapsto f(t,x) voor bijna elke xx continu is in t0t_0, en dat er een integreerbare gg is met f(t,x)g(x)\abs{f(t,x)} \leq g(x) voor alle tt in een omgeving van t0t_0 en bijna elke xx. Dan is FF continu in t0t_0.

Bewijs. Voor elke rij tnt0t_n \to t_0 is f(tn,)f(t0,)f(t_n, \cdot) \to f(t_0, \cdot) bijna overal, gedomineerd door gg: de gedomineerde convergentie geeft F(tn)F(t0)F(t_n) \to F(t_0); en in metrische ruimten volstaat de rijcontinuïteit (Opmerking 6.8).

Stelling 10.15 (Differentiëren onder de integraal)

Zij TT een open interval van R\R. Stel dat tf(t,x)t \mapsto f(t,x) voor bijna elke xx differentieerbaar is op TT, met

ft(t,x)g(x)voor alle tT en bijna elke x,\Bigl|\frac{\partial f}{\partial t}(t, x)\Bigr| \leq g(x) \quad \text{voor alle } t \in T \text{ en bijna elke } x,

met gg integreerbaar. Dan is FF differentieerbaar op TT met F(t)=Xft(t,x) ⁣dμ(x)F'(t) = \int_X \frac{\partial f}{\partial t}(t, x)\,\dd\mu(x).

Bewijs. Leg tt vast en zij hn0h_n \to 0: de differentiequotiënten

φn(x)=f(t+hn,x)f(t,x)hnft(t,x)bijna overal,\varphi_n(x) = \frac{f(t + h_n, x) - f(t, x)}{h_n} \longrightarrow \frac{\partial f}{\partial t}(t,x) \quad\text{bijna overal},

en de middelwaardeongelijkheid begrenst φn(x)supstf(s,x)g(x)\abs{\varphi_n(x)} \leq \sup_{s}\abs{\partial_tf(s,x)} \leq g(x): de gedomineerde convergentie is van toepassing, en F(t+hn)F(t)hn=φntf(t,)\frac{F(t + h_n) - F(t)}{h_n} = \int\varphi_n \to \int\partial_t f(t, \cdot).

Voorbeeld 10.16 (De gammafunctie)

Stel voor t>0t > 0

Γ(t)=0+xt1ex ⁣dx.\Gamma(t) = \int_0^{+\infty} x^{t-1}\eu^{-x}\,\dd x .

De integraal convergeert: nabij 00 is xt1x^{t-1} integreerbaar (t>0t > 0); en in het oneindige is xt1exCex/2x^{t-1}\eu^{-x} \leq C\eu^{-x/2}. Partiële integratie (op [ε,A][\varepsilon, A], met daarna limieten via de monotone convergentie) geeft de functionaalvergelijking Γ(t+1)=tΓ(t)\Gamma(t + 1) = t\,\Gamma(t), waaruit Γ(n+1)=n!\Gamma(n+1) = n!: de faculteit geïnterpoleerd. Op elke [a,b](0,+)\intcc ab \subseteq \intoo0{+\infty} wordt t(xt1ex)=lnxxt1ex\partial_t\bigl(x^{t-1}\eu^{-x}\bigr) = \ln x\cdot x^{t-1}\eu^{-x} gedomineerd door lnx(xa1+xb1)ex\abs{\ln x}(x^{a-1} + x^{b-1})\eu^{-x}, integreerbaar: dus is Γ\Gamma C1\mathcal C^1 en, per inductie, C\mathcal C^\infty, met Γ(k)(t)=0(lnx)kxt1ex ⁣dx\Gamma^{(k)}(t) = \int_0^\infty(\ln x)^kx^{t-1}\eu^{-x}\dd x. De waarde Γ(12)=π\Gamma(\frac12) = \sqrt\pi is de gaussische integraal in vermomming (Probleem 10.1).

De integrand xx: de oneigenlijke integraal ∈t_0∈fty convergeert door de afwisselende opheffing tussen de bogen, maar de oppervlakten | | van de bogen gedragen zich als 2π k — een harmonische reeks: xx ∉ L1. Lebesgue-integreerbaarheid is absolute integreerbaarheid.
De integrand sinxx\frac{\sin x}x: de oneigenlijke integraal 0\int_0^\infty convergeert door de afwisselende opheffing tussen de bogen, maar de oppervlakten \abs{\cdot} van de bogen gedragen zich als 2πk\frac2{\pi k} — een harmonische reeks: sinxxL1\frac{\sin x}x \notin L^1. Lebesgue-integreerbaarheid is absolute integreerbaarheid.

10.5 Oefeningen

Oefening 10.1

(a) Toon aan dat een monotone functie RR\R \to \R borel-meetbaar is, en dat een afgeleide (van een overal differentieerbare functie) borel-meetbaar is. (b) Toon aan dat f ⁣:XRf \colon X \to \R meetbaar is dan en slechts dan als {f>q}A\{f > q\} \in \mathcal A voor elke rationale qq.

Oplossing

Oplossing van Oefening 10.1.

(a) Is ff niet-dalend, dan is {f>t}\{f > t\} gelijk aan \varnothing, R\R of een halfrechte (a,+)\intoo a{+\infty} of [a,+)\intco a{+\infty}: in elk geval borel; en voor niet-stijgende ff net zo. Een afgeleide: f(x)=limnn(f(x+1n)f(x))f'(x) = \lim_n n\bigl(f(x + \frac1n) - f(x)\bigr) is een puntsgewijze limiet van continue (en dus meetbare) functies: Propositie 10.2(c).

(b) {f>t}=qQ,q>t{f>q}\{f > t\} = \bigcup_{q \in \Q,\, q > t}\{f > q\}: zijn de rationale niveaus meetbaar, dan alle niveaus, en de halfrechten brengen B(R)\mathcal B(\R) voort.

Oefening 10.2

Bereken, met volledige verantwoording:

limn0+cosx(1+x/n)n ⁣dx,limn01nxn11+x ⁣dx.\lim_{n\to\infty}\int_0^{+\infty} \frac{\cos x}{(1 + x/n)^{n}}\,\dd x, \qquad \lim_{n\to\infty}\int_0^1 \frac{n\,x^{n-1}}{1 + x}\,\dd x .

(Voor de tweede: substitueer u=xnu = x^n vóór het domineren.)

Oplossing

Oplossing van Oefening 10.2.

Eerste: (1+x/n)nex(1 + x/n)^n \nearrow \eu^x voor x0x \geq 0, dus gaat de integrand puntsgewijs naar excosx\eu^{-x}\cos x; en voor n2n \geq 2 is (1+x/n)n(1+x/2)2(1 + x/n)^n \geq (1 + x/2)^2, wat de integreerbare dominant (1+x/2)2(1 + x/2)^{-2} geeft. Gedomineerde convergentie:

limn0cosx(1+x/n)n ⁣dx=0excosx ⁣dx=Re0e(1i)x ⁣dx=Re11i=12.\lim_n\int_0^\infty\frac{\cos x}{(1 + x/n)^n}\dd x = \int_0^\infty \eu^{-x}\cos x\,\dd x = \operatorname{Re}\int_0^\infty\eu^{-(1 - \iu)x}\dd x = \operatorname{Re}\frac{1}{1 - \iu} = \frac12 .

Tweede: substitueer u=xnu = x^n (een C1\mathcal C^1-bijectie van (0,1)\intoo01):

01nxn11+x ⁣dx=01 ⁣du1+u1/n01 ⁣du2=12,\int_0^1\frac{nx^{n-1}}{1 + x}\dd x = \int_0^1\frac{\dd u}{1 + u^{1/n}} \longrightarrow \int_0^1\frac{\dd u}{2} = \frac12,

via de gedomineerde convergentie: voor u(0,1)u \in \intoo01 is u1/n1u^{1/n} \to 1, en de integrand is begrensd door 11 op een ruimte van eindige maat.

Oefening 10.3 ★★

(a) Geef een strikte ongelijkheid in het lemma van Fatou. (b) Geef op drie manieren een fn0f_n \to 0 puntsgewijs met fn=1\int f_n = 1: ontsnapping in hoogte, in breedte, en naar oneindig. Welke ene hypothese van de gedomineerde convergentie schendt elk van hen? (c) Toon aan dat men in het lemma van Fatou de lim inf\liminf aan geen van beide kanten door lim sup\limsup kan vervangen.

Oplossing

Oplossing van Oefening 10.3.

(a) fn=n1(0,1/n)f_n = n\,\mathbf 1_{\intoo0{1/n}}: puntsgewijs is lim inffn=0\liminf f_n = 0, en fn=1\int f_n = 1: dus 0<10 < 1.

(b) Hoogte: n1(0,1/n)n\mathbf 1_{\intoo0{1/n}}; breedte: 1n1(0,n)\frac1n\mathbf 1_{\intoo0n}; verschuiving: 1(n,n+1)\mathbf 1_{\intoo n{n+1}}. Alle gaan puntsgewijs naar 00 met =1\int = 1. In elk geval faalt de hypothese van de dominatie: supnfn\sup_nf_n is respectievelijk 1/x\approx 1/x nabij 00, ongeveer een niet-integreerbaar constant profiel, en van het type 1(1,)\mathbf 1_{\intoo1\infty} — nooit integreerbaar.

(c) “lim supfnlim supfn\int\limsup f_n \geq \limsup\int f_n” faalt voor de verschuivende bult: linkerlid 00, rechterlid 11. En “lim suplim sup\limsup\int \leq \int\limsup” is dezelfde uitspraak. Ook Fatou voor de lim sup\limsup met omgekeerde \leq (“omgekeerde Fatou”) vergt een dominant — dezelfde bult is het tegenvoorbeeld.

Oefening 10.4 ★★

(a) Toon aan dat 0+xex1 ⁣dx=n11n2=π26\displaystyle\int_0^{+\infty}\frac{x}{\eu^x - 1}\,\dd x = \sum_{n\geq1}\frac1{n^2} = \frac{\pi^2}6 (ontwikkel 1ex1\frac1{\eu^x - 1} in een meetkundige reeks en integreer term voor term — welke stelling staat dat toe?). (b) (De droom van de tweedejaars) Toon aan dat 01xx ⁣dx=n1nn\displaystyle\int_0^1 x^{-x}\,\dd x = \sum_{n\geq1}n^{-n}. (Schrijf xx=exlnx=k(xlnx)kk!x^{-x} = \eu^{-x\ln x} = \sum_k\frac{(-x\ln x)^k}{k!} en bereken 01(xlnx)k ⁣dx\int_0^1(-x\ln x)^k\dd x door x=eu/(k+1)x = \eu^{-u/(k+1)} te substitueren en Γ\Gamma te herkennen.)

Oplossing

Oplossing van Oefening 10.4.

(a) Voor x>0x > 0 is 1ex1=ex1ex=n1enx\frac1{\eu^x - 1} = \frac{\eu^{-x}}{1 - \eu^{-x}} = \sum_{n\geq1}\eu^{-nx}, dus xex1=n1xenx\frac{x}{\eu^x - 1} = \sum_{n\geq1}x\eu^{-nx}, een reeks van niet-negatieve meetbare functies: Gevolg 10.7 staat term voor term integreren toe:

0x ⁣dxex1=n10xenx ⁣dx=n11n2=π26\int_0^\infty\frac{x\,\dd x}{\eu^x - 1} = \sum_{n\geq1}\int_0^\infty x\eu^{-nx}\dd x = \sum_{n\geq1}\frac1{n^2} = \frac{\pi^2}6

(0xenx ⁣dx=n2\int_0^\infty x\eu^{-nx}\dd x = n^{-2} met partiële integratie; het Bazelprobleem uit het volume van bachelorjaar 2, of Oefening 13.5 verderop).

(b) Op (0,1)\intoo01 is xlnx0-x\ln x \geq 0, dus is xx=exlnx=k(xlnx)kk!x^{-x} = \eu^{-x\ln x} = \sum_k\frac{(-x\ln x)^k}{k!} een reeks met niet-negatieve termen: opnieuw verwisselen. Substitutie x=eu/(k+1)x = \eu^{-u/(k+1)} geeft

01(xlnx)k ⁣dx=0(uk+1)kekuk+1  euk+1k+1 ⁣du=1(k+1)k+10ukeu ⁣du=k!(k+1)k+1.\int_0^1(-x\ln x)^k\dd x = \int_0^\infty\Bigl(\frac{u}{k+1}\Bigr)^{k} \eu^{-\frac{ku}{k+1}}\;\frac{\eu^{-\frac u{k+1}}}{k+1}\,\dd u = \frac{1}{(k+1)^{k+1}}\int_0^\infty u^k\eu^{-u}\dd u = \frac{k!}{(k+1)^{k+1}} .

Bijgevolg is 01xx ⁣dx=k01(k+1)k+1=n1nn\int_0^1x^{-x}\dd x = \sum_{k\geq0}\frac{1}{(k+1) ^{k+1}} = \sum_{n\geq1}n^{-n}: de droom van de tweedejaars, streng bewezen.

Oefening 10.5 ★★

(a) Toon aan dat een meetbare f0f \geq 0 met f ⁣dμ=0\int f\,\dd\mu = 0 voldoet aan f=0f = 0 bijna overal. (Beschouw {f1/n}\{f \geq 1/n\} en de ongelijkheid van Markov: μ({fa})1af\mu(\{f \geq a\}) \leq \frac1a\int f — bewijs die.) (b) Toon aan dat een integreerbare ff bijna overal eindig is. (c) Toon aan dat uit Af ⁣dμ=0\int_A f\,\dd\mu = 0 voor elke meetbare AA volgt dat f=0f = 0 bijna overal.

Oplossing

Oplossing van Oefening 10.5.

(a) Markov: a1{fa}fa\,\mathbf 1_{\{f \geq a\}} \leq f; integreren geeft μ({fa})1af\mu(\{f \geq a\}) \leq \frac1a\int f. Is f=0\int f = 0, dan is μ({f1n})=0\mu(\{f \geq \frac1n\}) = 0 voor elke nn, en is {f>0}=n{f1n}\{f > 0\} = \bigcup_n\{f \geq \frac1n\} een nulverzameling.

(b) μ({f=})μ({fn})1nf0\mu(\{\abs f = \infty\}) \leq \mu(\{\abs f \geq n\}) \leq \frac1n\int\abs f \to 0.

(c) Neem A={f>0}A = \{f > 0\}: dan is f+ ⁣dμ=Af ⁣dμ=0\int f^+\dd\mu = \int_Af\,\dd\mu = 0, dus f+=0f^+ = 0 bijna overal volgens (a); en net zo f=0f^- = 0 bijna overal.

Oefening 10.6 ★★

(a) Pas Stelling 10.12 toe om de riemann-integreerbaarheid te beslissen van: 1Q\mathbf 1_\Q; de functie van Thomae; en 1K\mathbf 1_K voor KK een dikke cantorverzameling (Oefening 9.5). (b) Toon aan dat 1+sinxx ⁣dx=+\int_1^{+\infty}\abs{\frac{\sin x}x}\,\dd x = +\infty, terwijl limA1Asinxx ⁣dx\lim_{A\to\infty}\int_1^A\frac{\sin x}x\,\dd x bestaat (partiële integratie): oneigenlijke convergentie zonder integreerbaarheid.

Oplossing

Oplossing van Oefening 10.6.

(a) 1Q\mathbf 1_\Q: overal discontinu, dus niet riemann-integreerbaar (Stelling 10.12); haar lebesgue-integraal is λ(Q)=0\lambda(\Q) = 0. Thomae: continu in elk irrationaal punt (bij gegeven ε\varepsilon hebben slechts eindig veel rationale getallen in [0,1]\intcc01 een noemer 1/ε\leq 1/\varepsilon; mijd ze met een kleine omgeving) en discontinu in de rationale punten (dichtheid van de irrationale getallen): dus bijna overal continu, riemann-integreerbaar, met integraal 00 (ze verdwijnt bijna overal). 1K\mathbf 1_K met KK een dikke cantorverzameling: de discontinuïteitsverzameling is K=K\partial K = K (gesloten met leeg inwendige), van maat 12>0\frac12 > 0: dus niet riemann-integreerbaar — en toch lebesgue-integreerbaar met integraal λ(K)=12\lambda(K) = \frac12.

(b) kπ(k+1)πsinxx ⁣dx1(k+1)πkπ(k+1)πsinx ⁣dx=2(k+1)π\int_{k\pi}^{(k+1)\pi}\frac{\abs{\sin x}}x\dd x \geq \frac1{(k+1)\pi}\int_{k\pi}^{(k+1)\pi}\abs{\sin x}\dd x = \frac{2}{(k+1)\pi}: de reeks divergeert. Convergentie van de oneigenlijke integraal: voor A>πA > \pi is

πAsinxx ⁣dx=[cosxx]πAπAcosxx2 ⁣dx,\int_\pi^A\frac{\sin x}x\dd x = \Bigl[\frac{-\cos x}x\Bigr]_\pi^A - \int_\pi^A\frac{\cos x}{x^2}\dd x,

en beide termen convergeren als AA \to \infty (1x2\frac1{x^2} is integreerbaar): voorwaardelijke convergentie zonder absolute integreerbaarheid.

Oefening 10.7 ★★

Verantwoord dat F(t)=0+ex2cos(tx) ⁣dxF(t) = \int_0^{+\infty}\eu^{-x^2}\cos(tx)\,\dd x C1\mathcal C^1 is op R\R en voldoet aan F(t)=t2F(t)F'(t) = -\frac t2F(t) (partiële integratie); leid af dat F(t)=F(0)et2/4F(t) = F(0)\,\eu^{-t^2/4}. (Met F(0)=π2F(0) = \frac{\sqrt\pi}2 uit de weekendopgave: de gaussische functie is in wezen haar eigen fouriergetransformeerde — Hoofdstuk 14 zal dat systematiseren.)

Oplossing

Oplossing van Oefening 10.7.

Dominatie: t(ex2cos(tx))=xex2sin(tx)xex2\abs{\partial_t(\eu^{-x^2}\cos(tx))} = \abs{x\eu^{-x^2}\sin(tx)} \leq x\eu^{-x^2}, integreerbaar en onafhankelijk van tt: Stelling 10.15 is dus globaal van toepassing, en

F(t)=0xex2sin(tx) ⁣dx=[12ex2sin(tx)]0t20ex2cos(tx) ⁣dx=t2F(t)F'(t) = -\int_0^\infty x\eu^{-x^2}\sin(tx)\,\dd x = \Bigl[\tfrac12\eu^{-x^2}\sin(tx)\Bigr]_0^\infty - \frac t2\int_0^\infty\eu^{-x^2}\cos(tx)\dd x = -\frac t2F(t)

(partiële integratie met  ⁣dv=xex2 ⁣dx\dd v = x\eu^{-x^2}\dd x). De lineaire differentiaalvergelijking geeft F(t)=F(0)et2/4F(t) = F(0)\eu^{-t^2/4}; en met F(0)=π2F(0) = \frac{\sqrt\pi}2 (Probleem 10.1) reproduceert de gaussische functie zichzelf onder deze cosinustransformatie.

Oefening 10.8 ★★★

(Frullani) Zij 0<a<b0 < a < b. Toon aan dat

0+eaxebxx ⁣dx=lnba,\int_0^{+\infty}\frac{\eu^{-ax} - \eu^{-bx}}{x}\,\dd x = \ln\frac ba,

door de integrand als abext ⁣dt\int_a^b \eu^{-xt}\,\dd t te schrijven en de verwisseling te verantwoorden met de niet-negatieve theorie (Gevolg 10.7 in continue vorm — loop op Tonelli vooruit, of snijd [a,b][a,b] in nn gelijke stukken en ga naar de limiet).

Oplossing

Oplossing van Oefening 10.8.

De integraal convergeert: nabij 00 gaat de integrand naar bab - a (begrensd), en in het oneindige neemt ze af als eax\eu^{-ax}. Leg aa vast en bekijk I(b)=0eaxebxx ⁣dxI(b) = \int_0^\infty\frac{\eu^{-ax} - \eu^{-bx}}x\dd x als functie van b[a,+)b \in \intco a{+\infty}. Voor alle bab \geq a is b(integrand)=ebxeax\abs{\partial_b(\text{integrand})} = \eu^{-bx} \leq \eu^{-ax}, en 0eax ⁣dx=1a<\int_0^\infty\eu^{-ax}\dd x = \frac1a < \infty: een integreerbare dominant. Dus geeft Stelling 10.15 dat I(b)=0ebx ⁣dx=1bI'(b) = \int_0^\infty\eu^{-bx}\dd x = \frac1b, met I(a)=0I(a) = 0:

I(b)=ab ⁣dtt=lnba.I(b) = \int_a^b\frac{\dd t}t = \ln\frac ba .

(Gelijkwaardig, langs de route van de hint: de integrand is abext ⁣dt0\int_a^b\eu^{-xt}\dd t \geq 0 en de verwisseling is het continue analogon van Gevolg 10.7, dat wil zeggen Tonelli — bewezen in Hoofdstuk 11; de parameterroute blijft binnen dit hoofdstuk.)

Oefening 10.9 ★★

Zij f0f \geq 0 meetbaar op (X,A,μ)(X, \mathcal A, \mu). Toon aan dat ν(A)=Af ⁣dμ\nu(A) = \int_A f\,\dd\mu een maat definieert (met dichtheid ff ten opzichte van μ\mu), en dat g ⁣dν=gf ⁣dμ\int g\,\dd\nu = \int gf\,\dd\mu voor alle meetbare g0g \geq 0 (bewijs het voor indicatorfuncties, dan voor elementaire functies, en dan met de monotone convergentie — de standaardmachine).

Oplossing

Oplossing van Oefening 10.9.

ν()=0\nu(\varnothing) = 0; en voor disjuncte (An)(A_n) is f1An=nf1Anf\mathbf 1_{\bigsqcup A_n} = \sum_nf\mathbf 1_{A_n} (puntsgewijs, alle termen 0\geq 0), zodat Gevolg 10.7 de σ\sigma-additiviteit geeft. De formule g ⁣dν=gf ⁣dμ\int g\,\dd\nu = \int gf\,\dd\mu: voor g=1Ag = \mathbf 1_A is ze de definitie van ν\nu; voor elementaire gg volgt ze uit de lineariteit; en voor meetbare g0g \geq 0 neem je elementaire sngs_n \nearrow g (Stelling 10.4): dan is snfgfs_nf \nearrow gf, en gaat de monotone convergentie aan beide kanten tot de limiet over. (Die roltrap “indicator \to elementair \to monotone convergentie” is de standaardmachine van de theorie.)

Oefening 10.10 ★★★

(Een falen in de trant van Weierstrass) Definieer f(t)=0+sin(tx)x(1+x2) ⁣dxf(t) = \int_0^{+\infty}\frac{\sin(tx)}{x(1 + x^2)}\,\dd x. (a) Toon aan dat ff goed gedefinieerd en continu is op R\R, en C1\mathcal C^1 met f(t)=0cos(tx)1+x2 ⁣dxf'(t) = \int_0^\infty\frac{\cos(tx)}{1 + x^2}\dd x voor elke tt — maar dat nogmaals onder de integraal differentiëren ongeoorloofd is. (b) Neem aan dat f(t)=π2etf'(t) = \frac\pi2\eu^{-t} voor t>0t > 0 (bewezen in Hoofdstuk 17): wat is dan 0xsin(tx)1+x2 ⁣dx\int_0^\infty\frac{x\sin(tx)}{1+x^2}\dd x voor t>0t > 0, en waarom bevestigt die formule het falen uit (a)?

Oplossing

Oplossing van Oefening 10.10.

(a) Uit sin(tx)tx\abs{\sin(tx)} \leq \abs tx volgt sin(tx)x(1+x2)t1+x2\abs{\frac{\sin(tx)}{x(1+x^2)}} \leq \frac{\abs t}{1+x^2}: de integraal convergeert, en op tT\abs t \leq T levert de dominant T1+x2\frac{T}{1+x^2} de continuïteit (Stelling 10.14). Differentiëren: t=cos(tx)1+x211+x2\abs{\partial_t} = \abs{\frac{\cos(tx)}{1+x^2}} \leq \frac1{1+x^2}, integreerbaar, dus f(t)=0cos(tx)1+x2 ⁣dxf'(t) = \int_0^\infty\frac{\cos(tx)}{1+x^2}\dd x voor alle tt. Een tweede differentiatie zou het integreren van xsin(tx)1+x2\frac{x\sin(tx)}{1 + x^2} vergen, waarvan de absolute waarde zich in het oneindige gedraagt als sin(tx)x\frac{\abs{\sin(tx)}}x: niet integreerbaar — er is geen dominant en Stelling 10.15 kan geen tweede keer worden toegepast.

(b) Nemen we f(t)=π2etf'(t) = \frac\pi2\eu^{-t} voor t>0t > 0 aan, dan is ff' wegens de oneven aard van ff even, dus f(t)=π2etf'(t) = \frac\pi2\eu^{-\abs t} — en dat is niet differentieerbaar in 00: ff is C1\mathcal C^1 maar niet C2\mathcal C^2, wat bevestigt dat de geblokkeerde tweede differentiatie geen technisch toeval was. Voor t>0t > 0 is de oneigenlijke integraal 0xsin(tx)1+x2 ⁣dx\int_0^\infty\frac{x\sin(tx)}{1+x^2}\dd x gelijk aan f(t)=π2et-f''(t) = \frac\pi2\eu^{-t} (differentieer de aangenomen formule waar dat geoorloofd is, dus op (0,)\intoo0\infty) — een oneigenlijke, niet-lebesguewaarde.

Oefening 10.11 ★★

(Lemma van Scheffé) Zij fn,f0f_n, f \geq 0 integreerbaar met fnff_n \to f bijna overal en fnf\int f_n \to \int f. (a) Toon aan dat fnf0\int\abs{f_n - f} \to 0. (Pas de gedomineerde convergentie toe op gn=(ffn)+fg_n = (f - f_n)^+ \leq f en schrijf fnf=2gn(ffn)\int\abs{f_n - f} = 2\int g_n - \int(f - f_n).) (b) Toon met een voorbeeld aan dat de hypothese fnf\int f_n \to \int f niet gemist kan worden (een schuivende of samentrekkende bult), en dat de conclusie voor fnf_n met wisselend teken zonder controle op de absolute waarde faalt: fn=n1(0,1/n]n1(1/n,0]f_n = n\mathbf 1_{\intoc0{1/n}} - n\mathbf 1_{\intoc{-1/n}0} heeft fn0f_n \to 0 bijna overal en fn=00\int f_n = 0 \to 0, en toch fn=2\int\abs{f_n} = 2. (c) Toepassing (dichtheden): convergeren kansdichtheden pnpp_n \to p bijna overal, dan is automatisch pnp0\int\abs{p_n - p} \to 0: puntsgewijze convergentie van dichtheden is L1L^1-convergentie — een gratis opwaardering van de convergentie.

Oplossing

Oplossing van Oefening 10.11.

(a) Zij gn=(ffn)+g_n = (f - f_n)^+: dan is 0gnf0 \leq g_n \leq f (want fn0f_n \geq 0), gn0g_n \to 0 bijna overal, en ff een integreerbare dominant: dus gn0\int g_n \to 0 (gedomineerde convergentie). En omdat fnf=2(ffn)+(ffn)\abs{f_n - f} = 2(f - f_n)^+ - (f - f_n), is

fnf=2gn(ffn)0+0.\int\abs{f_n - f} = 2\int g_n - \Bigl(\int f - \int f_n\Bigr) \longrightarrow 0 + 0 .

(b) De schuivende bult fn=1[n,n+1]f_n = \mathbf 1_{\intcc n{n+1}} heeft fn0f_n \to 0 bijna overal en fn=1↛0\int f_n = 1 \not\to 0: zonder de convergentie van de integralen faalt de L1L^1-convergentie (en de hypothese ook). Het voorbeeld met wisselend teken: fn0f_n \to 0 in elke x0x \neq 0 en fn=0\int f_n = 0, maar fn=2\int\abs{f_n} = 2: voor rijen met wisselend teken gaat de stelling werkelijk over controle van het type fn\abs{f_n}, en de positiviteit werd precies gebruikt in gnfg_n \leq f.

(c) Dichtheden voldoen aan pn=1=p\int p_n = 1 = \int p: de hypothese van (a) geldt automatisch, dus dwingt pnpp_n \to p bijna overal af dat pnpL10\norm{p_n - p}_{L^1} \to 0 — en daarmee de convergentie van de kansen ApnAp\int_Ap_n \to \int_Ap, uniform over alle meetbare AA (A(pnp)pnp1\abs{\int_A(p_n - p)} \leq \norm{p_n - p}_1): Scheffé maakt van puntsgewijze convergentie van dichtheden convergentie van de wetten in totale variatie.

Oefening 10.12 ★★

Klassieke limieten, met volledige verantwoording via de monotone en de gedomineerde convergentie:

(a) limn0n(1xn)nex/2 ⁣dx,(b) limn01nxn11+x ⁣dx,\text{(a)}\ \lim_{n\to\infty}\int_0^n\Bigl(1 - \frac xn\Bigr)^n\eu^{x/2}\,\dd x, \qquad \text{(b)}\ \lim_{n\to\infty}\int_0^1\frac{n\,x^{n-1}}{1 + x}\,\dd x,
(c) limn0 ⁣dx(1+x/n)nx1/n.\text{(c)}\ \lim_{n\to\infty}\int_0^\infty \frac{\dd x}{(1 + x/n)^n\,x^{1/n}} .

(Voor (a): (1x/n)nex(1 - x/n)^n \nearrow \eu^{-x} voor vaste xx — bewijs de monotonie via log\log; voor (b): integreer partieel of substitueer x=u1/nx = u^{1/n} en herken een concentratie bij de rand; voor (c): zoek een integreerbare dominant die voor alle n2n \geq 2 geldt, door bij x=1x = 1 te splitsen.)

Oplossing

Oplossing van Oefening 10.12.

(a) Op (0,n)\intoo0n stijgt φn(x)=nlog(1xn)\varphi_n(x) = n\log(1 - \frac xn) in nn naar x-x (de afbeelding tlog(1xt)tt \mapsto \frac{\log(1 - xt)}{t} daalt als t=1n0t = \frac1n \downarrow 0; of ontwikkel: φn0\varphi_{n}' \geq 0 in nn via log(1u)+u1u0\log(1-u) + \frac{u}{1-u} \geq 0). Dus (1xn)nex/21x<nex/2(1 - \frac xn)^n\eu^{x/2}\mathbf 1_{x<n} \nearrow \eu^{-x/2}, en de monotone convergentie geeft

limn0n(1xn)nex/2 ⁣dx=0ex/2 ⁣dx=2.\lim_n\int_0^n\Bigl(1 - \frac xn\Bigr)^n\eu^{x/2}\dd x = \int_0^\infty\eu^{-x/2}\dd x = 2 .

(b) Substitueer u=xnu = x^n (dus x=u1/nx = u^{1/n} en nxn1 ⁣dx= ⁣dun x^{n-1}\dd x = \dd u):

01nxn11+x ⁣dx=01 ⁣du1+u1/n.\int_0^1\frac{nx^{n-1}}{1 + x}\dd x = \int_0^1\frac{\dd u}{1 + u^{1/n}} .

Voor u(0,1)u \in \intoo01 is u1/n1u^{1/n} \to 1, dus gaat de integrand naar 12\frac12, gedomineerd door 11: de limiet is 12\frac12 (gedomineerde convergentie). (De massa van nxn1nx^{n-1} concentreert zich bij x=1x = 1, waar 11+x=12\frac1{1+x} = \frac12: de substitutie maakt die concentratie zichtbaar.)

(c) Puntsgewijs is (1+x/n)nex(1 + x/n)^n \nearrow \eu^x en x1/n1x^{1/n} \to 1 (voor x>0x > 0): de integrand gaat naar ex\eu^{-x}. Dominant voor n2n \geq 2: op (0,1]\intoc01 is x1/nx1/2x^{-1/n} \leq x^{-1/2} en (1+x/n)n1(1 + x/n)^{-n} \leq 1, wat de integreerbare grens x1/2x^{-1/2} geeft; en op (1,)\intoo1\infty is x1/n1x^{-1/n} \leq 1 en (1+x/n)n1+(n2)x2n21+x24(1 + x/n)^n \geq 1 + \binom n2\frac{x^2}{n^2} \geq 1 + \frac{x^2}4, wat de integreerbare grens 44+x2\frac{4}{4 + x^2} geeft. Gedomineerde convergentie:

limn0 ⁣dx(1+x/n)nx1/n=0ex ⁣dx=1.\lim_n\int_0^\infty\frac{\dd x}{(1 + x/n)^nx^{1/n}} = \int_0^\infty\eu^{-x}\dd x = 1 .

10.6 Probleem: twee beroemde integralen

Probleem 10.1

Weekendopgave — de gaussische integraal en de integraal van Dirichlet, met parameters alleen

Twee integralen beheersen de toegepaste analyse:

G=+ex2 ⁣dx=π,D=0+sinxx ⁣dx=π2G = \int_{-\infty}^{+\infty}\eu^{-x^2}\dd x = \sqrt\pi, \qquad D = \int_0^{+\infty}\frac{\sin x}{x}\,\dd x = \frac\pi2

(de tweede als oneigenlijke integraal, Voorbeeld 10.13). We bewijzen beide met uitsluitend het gereedschap van dit hoofdstuk.

Deel I — De gaussische integraal. Stel voor t0t \geq 0

A(t)=(0tex2 ⁣dx)2,B(t)=01et2(1+x2)1+x2 ⁣dx.A(t) = \Bigl(\int_0^t\eu^{-x^2}\dd x\Bigr)^{2}, \qquad B(t) = \int_0^1\frac{\eu^{-t^2(1 + x^2)}}{1 + x^2}\,\dd x .
  1. Verantwoord dat AA en BB C1\mathcal C^1 zijn op (0,+)\intoo0{+\infty} en bereken AA' en BB'; toon aan dat A(t)+B(t)=0A'(t) + B'(t) = 0. (Substitueer in BB': u=txu = tx.)
  2. Bereken A(0)+B(0)A(0) + B(0) en limt+(A+B)(t)\lim_{t\to+\infty}(A + B)(t) — verantwoord de limiet onder de integraal in BB.
  3. Besluit dat 0ex2 ⁣dx=π2\int_0^\infty \eu^{-x^2}\dd x = \frac{\sqrt\pi}2, dus G=πG = \sqrt\pi, en leid af dat Γ(12)=π\Gamma(\tfrac12) = \sqrt\pi (substitueer x=u2x = u^2 in Γ(12)\Gamma(\frac12)).

Deel II — De integraal van Dirichlet. Stel voor t0t \geq 0

F(t)=0+etxsinxx ⁣dx.F(t) = \int_0^{+\infty}\eu^{-tx}\,\frac{\sin x}{x}\,\dd x .
  1. Toon aan dat de integraal die F(t)F(t) definieert voor elke t>0t > 0 als lebesgue-integraal convergeert, en voor t=0t = 0 als oneigenlijke integraal; toon aan dat D=limA0Asinxx ⁣dxD = \lim_{A\to\infty}\int_0^A\frac{\sin x}x\dd x bestaat (partiële integratie op [π,A][\pi, A]).
  2. Toon aan dat FF C1\mathcal C^1 is op (0,+)\intoo0{+\infty} met

    F(t)=0+etxsinx ⁣dx=11+t2F'(t) = -\int_0^{+\infty}\eu^{-tx}\sin x\,\dd x = -\frac{1}{1 + t^2}

    (dominatie op [t0,)[t_0, \infty) voor elke t0>0t_0 > 0; de laatste integraal met twee partiële integraties of met complexe exponentiëlen).

  3. Toon aan dat F(t)0F(t) \to 0 als t+t \to +\infty, en leid af dat F(t)=π2arctantF(t) = \frac\pi2 - \arctan t op (0,+)\intoo0{+\infty}.
  4. Het delicate punt: D=limt0+F(t)D = \lim_{t\to0^+}F(t). Bewijs het met uniforme controle op de staart: toon met partiële integratie aan dat voor 0t10 \leq t \leq 1 en AπA \geq \pi

    A+etxsinxx ⁣dxCA\Bigl|\int_A^{+\infty}\eu^{-tx}\frac{\sin x}x\,\dd x\Bigr| \leq \frac{C}{A}

    met CC onafhankelijk van tt (differentieer etxx\frac{\eu^{-tx}}x en begrens cos\abs{\cos} door 11; merk op dat tetx1/x(txetx)t\eu^{-tx} \leq 1/x\cdot(tx\eu^{-tx}) met supu0ueu<1\sup_{u\geq0}u\eu^{-u} < 1); splits daarna F(t)DF(t) - D in [0,A][0, A] (waar de gedomineerde convergentie werkt als t0t \to 0) en [A,)[A, \infty).

  5. Besluit: D=π2D = \frac\pi2.

Deel III — Dividenden.

  1. Bereken 0+sin2xx2 ⁣dx\int_0^{+\infty}\frac{\sin^2x}{x^2}\,\dd x (integreer partieel en herleid tot DD via sin2x=2sinxcosx\sin 2x = 2\sin x\cos x).
  2. Bereken 0+1cosxx2 ⁣dx\int_0^{+\infty}\frac{1 - \cos x}{x^2}\,\dd x, en controleer of de twee resultaten met elkaar stroken.
  3. Bereken voor a>0a > 0 de integralen 0+sin(ax)x ⁣dx\int_0^{+\infty}\frac{\sin(ax)}x\dd x en +eax2 ⁣dx\int_{-\infty}^{+\infty}\eu^{-ax^2}\dd x, en noteer de schaalregels (ze worden de werkpaarden van Hoofdstuk 14).
  4. Leg precies uit waarom DD niet rechtstreeks met de gedomineerde convergentie in t=0t = 0 behandeld kon worden (er is geen integreerbare dominant op [0,1]×[0,)[0,1]\times[0,\infty)), en waarom het splitsen van de staart in vraag 7 het eerlijke alternatief is — dat patroon (“uniforme integreerbaarheid van de staarten”) keert overal in de analyse terug.

Deel IV — De gammafunctie volgens Bohr en Mollerup. De functie Γ\Gamma (Voorbeeld 10.16) voldoet aan Γ(1)=1\Gamma(1) = 1 en Γ(x+1)=xΓ(x)\Gamma(x+1) = x\Gamma(x) — maar dat doen oneindig veel andere functies ook (vermenigvuldig met een willekeurige 11-periodieke golf). Eén convexiteitsvoorwaarde legt Γ\Gamma eenduidig vast, en dan volgen haar diepere identiteiten mechanisch. Een positieve functie ff op een interval heet logaritmisch convex als logf\log f convex is.

  1. Toon aan dat logaritmisch convex convex impliceert, dat producten van logaritmisch convexe functies en hun samenstellingen met affiene afbeeldingen logaritmisch convex zijn, en — via de ongelijkheid van Hölder voor twee functies, uv(up)1/p(vq)1/q\int\abs{uv} \leq \bigl(\int\abs u^p\bigr)^{1/p}\bigl(\int\abs v^q\bigr)^{1/q}, rechtstreeks uit de ongelijkheid van Young bewezen — dat Γ\Gamma logaritmisch convex is op (0,)\intoo0\infty.
  2. (Hellingslemma) Zij gg convex op (0,)\intoo0\infty met g(n+1)g(n)=logng(n+1) - g(n) = \log n voor elk geheel getal n1n \geq 1. Vergelijk voor x(0,1]x \in \intoc01 en n2n \geq 2 de hellingen van gg over [n1,n][n-1, n], [n,n+x][n, n+x] en [n,n+1][n, n+1], en leid af dat

    xlog(n1)    g(n+x)g(n)    xlogn.x\log(n-1) \;\leq\; g(n + x) - g(n) \;\leq\; x\log n .
  3. (Bohr–Mollerup) Zij f>0f > 0 met f(1)=1f(1) = 1, f(x+1)=xf(x)f(x+1) = xf(x) en logf\log f convex. Rol de recursie af tot f(n+x)=x(x+1)(x+n1)f(x)f(n + x) = x(x+1)\cdots(x + n - 1)\,f(x) en f(n)=(n1)!f(n) = (n-1)!, en leid met vraag 14 af dat voor x(0,1]x \in \intoc01

    f(x)=limnn!nxx(x+1)(x+n):f(x) = \lim_{n\to\infty} \frac{n!\,n^x}{x(x+1)\cdots(x+n)} :

    ff is uniek, dus f=Γf = \Gamma, en de limietformule van Gauss geldt (breid met de recursie uit tot alle x>0x > 0).

  4. Definieer de bètafunctie B(x,y)=01tx1(1t)y1 ⁣dtB(x, y) = \int_0^1t^{x-1}(1-t)^{y-1}\,\dd t (met x,y>0x, y > 0). Bewijs de convergentie, de recursie B(x+1,y)=xx+yB(x,y)B(x+1, y) = \frac{x}{x+y}\,B(x, y) (partiële integratie), en B(1,y)=1yB(1, y) = \frac1y.
  5. Toon aan dat xB(x,y)x \mapsto B(x, y) logaritmisch convex is (opnieuw Hölder), en pas Bohr–Mollerup toe op

    f(x)=B(x,y)Γ(x+y)Γ(y)f(x) = \frac{B(x, y)\,\Gamma(x + y)}{\Gamma(y)}

    om de formule van Euler te besluiten: B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x, y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} — zonder ook maar één dubbele integraal.

  6. Bereken B(12,12)B(\frac12, \frac12) rechtstreeks (substitueer t=sin2θt = \sin^2\theta) en leid af dat Γ(12)=π\Gamma(\frac12) = \sqrt\pi: de gaussische integraal van Deel I, uit zuivere convexiteit teruggevonden. Vergelijk de twee bewijzen in telkens één zin.
  7. (Verdubbeling van Legendre) Toon aan dat

    g(x)=2x1πΓ(x2)Γ(x+12)g(x) = \frac{2^{x-1}}{\sqrt\pi}\, \Gamma\Bigl(\frac x2\Bigr) \Gamma\Bigl(\frac{x+1}2\Bigr)

    aan de drie hypothesen van Bohr–Mollerup voldoet, en besluit dat g=Γg = \Gamma, dat wil zeggen Γ(2z)=22z1πΓ(z)Γ(z+12)\Gamma(2z) = \frac{2^{2z-1}}{\sqrt\pi}\,\Gamma(z)\,\Gamma(z + \tfrac12) voor alle z>0z > 0.

  8. Leid de gesloten vorm Γ(n+12)=(2n)!4nn!π\Gamma\bigl(n + \tfrac12\bigr) = \dfrac{(2n)!}{4^n\,n!}\sqrt\pi af, en bewijs met het hellingslemma toegepast op logΓ\log\Gamma rond grote gehele getallen de asymptotiek

    Γ(n+12)Γ(n)n1.\frac{\Gamma(n + \frac12)}{\Gamma(n)\,\sqrt n} \longrightarrow 1 .
  9. Combineer de laatste twee vragen tot de asymptotiek van de centrale binomiaalcoëfficiënt

    (2nn)4nπn,\binom{2n}{n} \sim \frac{4^n}{\sqrt{\pi n}},

    en controleer numeriek voor n=10n = 10 ((2010)=184756\binom{20}{10} = 184756, tegenover 410/10π1870794^{10}/\sqrt{10\pi} \approx 187079: verhouding 0.988\approx 0.988).

  10. (Synthese) De constante π\sqrt\pi is nu verschenen als de gaussische integraal (Deel I), als B(12,12)B(\frac12, \frac12) (vraag 18) en binnen de verdubbeling (vraag 19); en de schatting van de centrale binomiaalcoëfficiënt loopt vooruit op zowel Stirling (de weekendopgave van Hoofdstuk 11) als de Moivre–Laplace. Breng de verbanden in kaart: welke uitspraken zijn met welke gelijkwaardig, en wat draagt elke techniek — differentiëren onder de integraal tegenover convexiteit — bij dat de andere niet kan?

Deel V — Nog drie dividenden.

  1. (Wallis, via de bètafunctie) Substitueer voor p>1p > -1 de variabele t=sin2θt = \sin^2\theta om aan te tonen dat

    Wp=0π/2sinpθ ⁣dθ=12B(p+12,12),W_p = \int_0^{\pi/2}\sin^p\theta\,\dd\theta = \frac12\,B\Bigl(\frac{p+1}2, \frac12\Bigr),

    en leid uit de bètarecursie (vraag 16) af dat Wn+2=n+1n+2WnW_{n+2} = \frac{n+1}{n+2}\,W_n voor gehele n0n \geq 0. Bereken W2nW_{2n} en W2n+1W_{2n+1} in gesloten vorm, toon met insluiting aan dat W2n+1/W2n1W_{2n+1}/W_{2n} \to 1, en besluit met het product van Wallis

    π2=limnk=1n4k24k21.\frac\pi2 = \lim_{n\to\infty} \prod_{k=1}^{n}\frac{4k^2}{4k^2 - 1} .
  2. (De gaussische functie ontmoet een frequentie) Stel voor bRb \in \R

    Φ(b)=+ex2cos(2bx) ⁣dx.\Phi(b) = \int_{-\infty}^{+\infty} \eu^{-x^2}\cos(2bx)\,\dd x .

    Toon aan dat Φ\Phi C1\mathcal C^1 is op R\R, dat een partiële integratie de differentiaalvergelijking Φ(b)=2bΦ(b)\Phi'(b) = -2b\,\Phi(b) oplevert, en besluit

    Φ(b)=πeb2:\Phi(b) = \sqrt\pi\,\eu^{-b^2} :

    de gaussische functie reproduceert zichzelf onder deze transformatie — de ene identiteit waarop Hoofdstuk 14 zal draaien.

  3. (De integraal van Frullani) Toon voor 0<a<b0 < a < b aan dat

    0+eaxebxx ⁣dx=logba,\int_0^{+\infty} \frac{\eu^{-ax} - \eu^{-bx}}{x}\,\dd x = \log\frac ba ,

    door naar de parameter aa te differentiëren (verantwoord de dominatie op elke [a0,+)\intco{a_0}{+\infty} met a0>0a_0 > 0, en bepaal de constante door aba \to b te laten gaan). Waar precies heeft de integrand haar ophefbare singulariteit in x=0x = 0 nodig?

Oplossing

Oplossing van Probleem 10.1.

1. AA is C1\mathcal C^1 wegens de hoofdstelling van de integraalrekening en de kettingregel: A(t)=2et20tex2 ⁣dxA'(t) = 2\eu^{-t^2}\int_0^t\eu^{-x^2}\dd x. Voor BB is t(et2(1+x2)1+x2)=2tet2(1+x2)\partial_t\bigl(\frac{\eu^{-t^2(1+x^2)}}{1+x^2}\bigr) = -2t\,\eu^{-t^2(1+x^2)}, continu en begrensd op [t0,T]×[0,1][t_0, T] \times \intcc01 voor elke 0<t0<T0 < t_0 < T (een begrensde dominant op een ruimte van eindige maat volstaat): BB is dus C1\mathcal C^1 op (0,+)\intoo0{+\infty} met

B(t)=2t01et2(1+x2) ⁣dx=2et201te(tx)2 ⁣dx=2et20teu2 ⁣du=A(t)B'(t) = -2t\int_0^1 \eu^{-t^2(1+x^2)}\dd x = -2\eu^{-t^2}\int_0^1 t\,\eu^{-(tx)^2}\dd x = -2\eu^{-t^2}\int_0^t\eu^{-u^2}\dd u = -A'(t)

(substitutie u=txu = tx).

2. A(0)=0A(0) = 0 en B(0)=01 ⁣dx1+x2=π4B(0) = \int_0^1\frac{\dd x}{1+x^2} = \frac\pi4. A+BA + B heeft op (0,+)\intoo0{+\infty} afgeleide nul en is continu in 00 (BB via de dominant 11+x2\frac1{1+x^2} en Stelling 10.14): dus A+Bπ4A + B \equiv \frac\pi4. Als t+t\to+\infty: 0B(t)et200 \leq B(t) \leq \eu^{-t^2} \to 0, en A(t)(0ex2 ⁣dx)2A(t) \to \bigl(\int_0^\infty\eu^{-x^2}\dd x\bigr)^2 (monotone convergentie van de binnenste integraal).

3. Bijgevolg is (0ex2 ⁣dx)2=π4\bigl(\int_0^\infty\eu^{-x^2}\dd x\bigr)^2 = \frac\pi4, dus 0ex2 ⁣dx=π2\int_0^\infty\eu^{-x^2}\dd x = \frac{\sqrt\pi}2, en wegens de pariteit G=πG = \sqrt\pi. Verder is Γ(12)=0x1/2ex ⁣dx=x=u220eu2 ⁣du=π\Gamma(\frac12) = \int_0^\infty x^{-1/2}\eu^{-x}\dd x \overset{x = u^2}{=} 2\int_0^\infty\eu^{-u^2}\dd u = \sqrt\pi.

4. Voor t>0t > 0 is etxsinxxetx\abs{\eu^{-tx}\frac{\sin x}x} \leq \eu^{-tx}, integreerbaar. Voor t=0t = 0 de oneigenlijke convergentie: op [π,A][\pi, A] is

πAsinxx ⁣dx=[cosxx]πAπAcosxx2 ⁣dx,\int_\pi^A\frac{\sin x}x\dd x = \Bigl[-\frac{\cos x}x\Bigr]_\pi^A - \int_\pi^A\frac{\cos x}{x^2}\dd x,

en beide termen convergeren als AA \to \infty; en nabij 00 zet de integrand zich continu voort met de waarde 11.

5. Op [t0,+)[t_0, +\infty) met t0>0t_0 > 0 is t(etxsinxx)=etxsinxet0x\abs{\partial_t\bigl(\eu^{-tx}\tfrac{\sin x}x\bigr)} = \eu^{-tx}\abs{\sin x} \leq \eu^{-t_0x}, integreerbaar: dus is Stelling 10.15 op elk zulk interval van toepassing, en daarmee op heel (0,+)\intoo0{+\infty}:

F(t)=0etxsinx ⁣dx=Im0e(ti)x ⁣dx=Im1ti=11+t2.F'(t) = -\int_0^\infty\eu^{-tx}\sin x\,\dd x = -\operatorname{Im}\int_0^\infty\eu^{-(t - \iu)x}\dd x = -\operatorname{Im}\frac{1}{t - \iu} = -\frac{1}{1 + t^2}.

6. F(t)0etx ⁣dx=1t0\abs{F(t)} \leq \int_0^\infty\eu^{-tx}\dd x = \frac1t \to 0. Integreren van F=11+t2F' = -\frac1{1+t^2} geeft F(t)=CarctantF(t) = C - \arctan t, en tt \to \infty dwingt C=π2C = \frac\pi2 af: dus F(t)=π2arctantF(t) = \frac\pi2 - \arctan t op (0,+)\intoo0{+\infty}.

7. Integreer partieel op [A,R][A, R] met sinx=(cosx)\sin x = (-\cos x)' en laat RR \to \infty gaan:

Aetxsinxx ⁣dx=cosA  etAAAcosx(tx+1x2)etx ⁣dx.\int_A^{\infty}\eu^{-tx}\frac{\sin x}x\dd x = \frac{\cos A\;\eu^{-tA}}{A} - \int_A^\infty \cos x\,\Bigl(\frac tx + \frac1{x^2}\Bigr)\eu^{-tx}\dd x .

Met cos1\abs{\cos} \leq 1: de eerste term is 1A\leq \frac1A; en de integraal is hoogstens Atetx ⁣dxx+A ⁣dxx21A0tetx ⁣dx+1A=2A\int_A^\infty t\eu^{-tx}\frac{\dd x}x + \int_A^\infty\frac{\dd x}{x^2} \leq \frac1A\int_0^\infty t\eu^{-tx}\dd x + \frac1A = \frac2A. Totaal: 3A\leq \frac 3A, uniform voor t[0,1]t \in [0, 1] (het geval t=0t = 0 inbegrepen). Nu is

F(t)D0A(etx1)sinxx ⁣dx+6A.\abs{F(t) - D} \leq \Bigl|\int_0^A(\eu^{-tx} - 1)\,\frac{\sin x}x\,\dd x\Bigr| + \frac6A .

Op [0,A][0, A] is etx1txtA\abs{\eu^{-tx} - 1} \leq tx \leq tA en sinxx1\abs{\frac{\sin x}x} \leq 1, dus is de eerste term hoogstens tA2tA^2. Kies AA met 6A<ε\frac6A < \varepsilon, en daarna t<ε/A2t < \varepsilon/A^2: dan is F(t)D<2ε\abs{F(t) - D} < 2\varepsilon.

8. Bijgevolg is D=limt0+F(t)=limt0+(π2arctant)=π2D = \lim_{t\to0^+}F(t) = \lim_{t\to0^+}\bigl(\frac\pi2 - \arctan t\bigr) = \frac\pi2.

9. Partiële integratie (u=sin2xu = \sin^2x, v=x2v' = x^{-2}):

0sin2xx2 ⁣dx=[sin2xx]0+02sinxcosxx ⁣dx=0sin2xx ⁣dx=D=π2\int_0^\infty\frac{\sin^2x}{x^2}\dd x = \Bigl[-\frac{\sin^2x}{x}\Bigr]_0^\infty + \int_0^\infty\frac{2\sin x\cos x}{x}\dd x = \int_0^\infty\frac{\sin 2x}{x}\dd x = D = \frac\pi2

(substitueer u=2xu = 2x in de laatste stap; de randtermen verdwijnen: sin2x/x0\sin^2 x/x \to 0 aan beide uiteinden).

10. Partiële integratie (u=1cosxu = 1 - \cos x, v=x2v' = x^{-2}): 01cosxx2 ⁣dx=0sinxx ⁣dx=π2\int_0^\infty\frac{1 - \cos x}{x^2}\dd x = \int_0^\infty\frac{\sin x}x\dd x = \frac\pi2. Consistentie: 1cosx=2sin2x21 - \cos x = 2\sin^2\frac x2, en de substitutie x=2ux = 2u maakt van 2sin2(x/2)x2 ⁣dx\int\frac{2\sin^2(x/2)}{x^2}\dd x de integraal sin2uu2 ⁣du\int\frac{\sin^2u}{u^2}\dd u: de twee berekeningen stroken.

11. 0sin(ax)x ⁣dx=π2\int_0^\infty\frac{\sin(ax)}x\dd x = \frac\pi2 voor elke a>0a > 0 (substitueer u=axu = ax: de integraal is schaalinvariant); en Reax2 ⁣dx=π/a\int_\R\eu^{-ax^2}\dd x = \sqrt{\pi/a} (substitueer u=axu = \sqrt a\,x). Schalen in het argument laat de integraal van Dirichlet onveranderd en deelt de gaussische integraal door a\sqrt a.

12. Een dominant die voor alle t[0,1]t \in [0,1] geldt, moet supt[0,1]etxsinxx=sinxx\sup_{t\in[0,1]}\abs{\eu^{-tx}\frac{\sin x}x} = \abs{\frac{\sin x}x} domineren, en dat is niet integreerbaar (Oefening 10.6): de gedomineerde convergentie kan t=0t = 0 dus niet oversteken. Het alternatief van vraag 7 — staarten uniform klein in de parameter, en het compacte deel met de gedomineerde convergentie — is het standaardpatroon van de “uniforme integreerbaarheid”, en keert terug zodra voorwaardelijke convergentie een limietverwisseling ontmoet.

13. Is g=logfg = \log f convex, dan is f=expgf = \exp\circ g convex (exp is convex en stijgend: f(λx+(1λ)y)eλg(x)+(1λ)g(y)λf(x)+(1λ)f(y)f(\lambda x + (1-\lambda)y) \leq \eu^{\lambda g(x) + (1-\lambda)g(y)} \leq \lambda f(x) + (1-\lambda)f(y), met de laatste stap wegens de convexiteit van exp tussen de punten g(x)g(x) en g(y)g(y)). Producten en affiene substituties: logaritmen maken er sommen en affiene substituties van convexe functies van. Hölder (met 1p+1q=1\frac1p + \frac1q = 1): voor up=vq=1\int\abs u^p = \int\abs v^q = 1 geeft Young dat uvupp+vqq\abs{uv} \leq \frac{\abs u^p}p + \frac{\abs v^q}q; integreren geeft uv1\int\abs{uv} \leq 1, en het algemene geval volgt uit de homogeniteit. Pas het dan voor λ(0,1)\lambda \in \intoo01 met p=1λp = \frac1\lambda toe op de ontbinding

tλx+(1λ)y1et=(tx1et)λ(ty1et)1λ:Γ(λx+(1λ)y)Γ(x)λΓ(y)1λ.t^{\lambda x + (1-\lambda)y - 1}\eu^{-t} = \bigl(t^{x-1}\eu^{-t}\bigr)^{\lambda} \bigl(t^{y-1}\eu^{-t}\bigr)^{1-\lambda} : \qquad \Gamma(\lambda x + (1{-}\lambda)y) \leq \Gamma(x)^\lambda\,\Gamma(y)^{1-\lambda} .

14. Voor een convexe gg stijgt de helling van een koorde met haar eindpunten (de driekoordenongelijkheid). Vergelijking van de koorden over [n1,n][n-1, n], [n,n+x][n, n+x] en [n,n+1][n, n+1] geeft

log(n1)=g(n)g(n1)1g(n+x)g(n)xg(n+1)g(n)1=logn,\log(n-1) = \frac{g(n) - g(n-1)}1 \leq \frac{g(n+x) - g(n)}x \leq \frac{g(n+1) - g(n)}1 = \log n,

en vermenigvuldigen met x>0x > 0 geeft de bewering.

15. Met f(n)=(n1)!f(n) = (n-1)! (recursie vanaf f(1)=1f(1) = 1) en f(n+x)=x(x+1)(x+n1)f(x)f(n + x) = x(x+1)\cdots(x+n-1)\,f(x) luidt vraag 14:

(n1)x(n1)!    x(x+1)(x+n1)f(x)    nx(n1)!.(n-1)^x\,(n-1)! \;\leq\; x(x+1)\cdots(x+n-1)\,f(x) \;\leq\; n^x\,(n-1)! .

De bovengrens herschrijft zich als f(x)n!nxx(x+1)(x+n)x+nnf(x) \leq \frac{n!\,n^x}{x(x+1)\cdots(x+n)}\cdot\frac{x+n}n, en de ondergrens op rang n+1n+1 als f(x)n!nxx(x+1)(x+n)f(x) \geq \frac{n!\,n^x}{x(x+1)\cdots(x+n)}. De correctiefactor x+nn1\frac{x+n}n \to 1: de insluiting dwingt

f(x)=limnn!nxx(x+1)(x+n)(x(0,1])f(x) = \lim_n\frac{n!\,n^x}{x(x+1)\cdots(x+n)} \qquad (x \in \intoc01)

af, een uitdrukking die niet van ff afhangt: eenduidigheid op (0,1]\intoc01, en dus overal wegens de recursie. En omdat Γ\Gamma aan alle drie de hypothesen voldoet (vraag 13), is f=Γf = \Gamma en geldt de formule van Gauss — voor alle x>0x > 0, want beide leden voldoen aan dezelfde recursie.

16. Nabij 00 is de integrand tx1\sim t^{x-1}, integreerbaar precies wanneer x>0x > 0; en nabij 11 symmetrisch met yy. Partiële integratie op [ε,1ε]\intcc\varepsilon{1-\varepsilon} en ε0\varepsilon \to 0 (de randtermen verdwijnen voor x,y>0x, y > 0):

B(x+1,y)=[tx(1t)yy]01+xy01tx1(1t)y ⁣dt=xy(B(x,y)B(x+1,y)),B(x{+}1, y) = \Bigl[-t^x\frac{(1-t)^y}y\Bigr]_0^1 + \frac xy\int_0^1t^{x-1}(1-t)^y\,\dd t = \frac xy\bigl(B(x, y) - B(x{+}1, y)\bigr),

met (1t)y=(1t)y1(1t)(1-t)^y = (1-t)^{y-1}(1 - t); oplossen geeft B(x+1,y)=xx+yB(x,y)B(x+1, y) = \frac{x}{x+y}B(x, y). En B(1,y)=01(1t)y1 ⁣dt=1yB(1, y) = \int_0^1(1-t)^{y-1}\dd t = \frac1y.

17. f(1)=1yΓ(1+y)Γ(y)=1f(1) = \frac1y\cdot\frac{\Gamma(1+y)}{\Gamma(y)} = 1; f(x+1)=xx+yB(x,y)(x+y)Γ(x+y)Γ(y)=xf(x)f(x+1) = \frac{x}{x+y}B(x,y)\cdot\frac{(x+y)\Gamma(x+y)} {\Gamma(y)} = x\,f(x); en ff is logaritmisch convex in xx als product van de logaritmisch convexe B(,y)B(\cdot, y) (Hölder op de ontbinding t(λx1+(1λ)x2)1(1t)y1=()λ()1λt^{(\lambda x_1 + (1-\lambda)x_2)-1}(1-t)^{y-1} = (\cdots)^\lambda(\cdots)^{1-\lambda}, als in vraag 13) en Γ(+y)\Gamma(\cdot + y) (affiene verschuiving). Bohr–Mollerup: f=Γf = \Gamma, dat wil zeggen B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x, y) = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}.

18. Met t=sin2θt = \sin^2\theta is  ⁣dt=2sinθcosθ ⁣dθ\dd t = 2\sin\theta\cos\theta\,\dd\theta en t1/2(1t)1/2=1sinθcosθt^{-1/2}(1-t)^{-1/2} = \frac1{\sin\theta\cos\theta}:

B(12,12)=0π/22 ⁣dθ=π=Γ(12)2Γ(1)Γ(12)=π.B\Bigl(\frac12, \frac12\Bigr) = \int_0^{\pi/2}2\,\dd\theta = \pi = \frac{\Gamma(\frac12)^2}{\Gamma(1)} \quad\Longrightarrow\quad \Gamma\Bigl(\frac12\Bigr) = \sqrt\pi .

Deel I bereikte dezelfde constante door naar een parameter te differentiëren en twee functies naar hun limieten te laten wedijveren; hier maakte de convexiteit alleen het probleem zo stijf dat er nog maar één waarde overbleef. Analyse door beweging tegenover analyse door vorm.

19. g(1)=20πΓ(12)Γ(1)=1g(1) = \frac{2^0}{\sqrt\pi}\Gamma(\frac12)\Gamma(1) = 1. Recursie:

g(x+1)=2xπΓ(x+12)Γ(x2+1)=2xπx2Γ(x2)Γ(x+12)=xg(x).g(x+1) = \frac{2^{x}}{\sqrt\pi}\, \Gamma\Bigl(\frac{x+1}2\Bigr)\Gamma\Bigl(\frac x2 + 1\Bigr) = \frac{2^{x}}{\sqrt\pi}\cdot\frac x2\, \Gamma\Bigl(\frac x2\Bigr)\Gamma\Bigl(\frac{x+1}2\Bigr) = x\,g(x) .

Logaritmische convexiteit: product van e(x1)log2\eu^{(x-1)\log2} (logaritmisch affien) en twee affiene herparametriseringen van de logaritmisch convexe Γ\Gamma. Bohr–Mollerup geeft g=Γg = \Gamma; en met x=2zx = 2z: Γ(2z)=22z1πΓ(z)Γ(z+12)\Gamma(2z) = \frac{2^{2z-1}}{\sqrt\pi}\Gamma(z)\Gamma(z + \frac12).

20. Uit Γ(12)=π\Gamma(\frac12) = \sqrt\pi en de recursie volgt Γ(n+12)=(n12)(n32)12π=(2n1)(2n3)12nπ=(2n)!4nn!π\Gamma(n + \frac12) = (n - \frac12)(n - \frac32)\cdots\frac12\, \sqrt\pi = \frac{(2n-1)(2n-3)\cdots1}{2^n}\sqrt\pi = \frac{(2n)!}{4^nn!}\sqrt\pi (vul het oneven product met de even getallen aan). Asymptotiek: vraag 14 met g=logΓg = \log\Gamma en x=12x = \frac12 geeft n1Γ(n+12)Γ(n)n\sqrt{n-1} \leq \frac{\Gamma(n+\frac12)}{\Gamma(n)} \leq \sqrt n, zodat de verhouding tot n\sqrt n tussen 11n\sqrt{1 - \frac1n} en 11 wordt ingeklemd.

21. Uit vraag 20 volgt (2n)!=4nn!πΓ(n+12)(2n)! = \frac{4^nn!}{\sqrt\pi}\Gamma(n + \tfrac12), dus

(2nn)=(2n)!(n!)2=4nΓ(n+12)π  n!=4nπΓ(n+12)nΓ(n)4nπnn=4nπn.\binom{2n}n = \frac{(2n)!}{(n!)^2} = \frac{4^n\,\Gamma(n+\frac12)}{\sqrt\pi\;n!} = \frac{4^n}{\sqrt\pi}\cdot \frac{\Gamma(n+\frac12)}{n\,\Gamma(n)} \sim \frac{4^n}{\sqrt\pi}\cdot\frac{\sqrt n}{n} = \frac{4^n}{\sqrt{\pi n}} .

Numeriek is 410/10π=1048576/5.6050187078.64^{10}/\sqrt{10\pi} = 1048576/5.6050 \approx 187078.6, tegenover (2010)=184756\binom{20}{10} = 184756: verhouding 0.98760.9876 — de fout is O(1/n)O(1/n), bij n=10n = 10 nog zichtbaar.

22. Gelijkwaardigheden: Γ(12)=πG=π\Gamma(\frac12) = \sqrt\pi \Leftrightarrow G = \sqrt\pi (de substitutie x=u2x = u^2 uit vraag 3) B(12,12)=π\Leftrightarrow B(\frac12, \frac12) = \pi (de formule van Euler); en de verdubbeling in z=nz = n is de gesloten vorm van Γ(n+12)\Gamma(n + \frac12), die op het hellingslemma na is de schatting van de centrale binomiaalcoëfficiënt. De parametertechniek (Delen I–II) berekent limieten van bewegende grootheden en is onmisbaar zodra er een echte vervorming in het spel is (de integraal van Dirichlet heeft geen convexiteitsbewijs); de convexiteitstechniek berekent niets maar verbiedt alles — ze blinkt uit in eenduidigheid en functionaalvergelijkingen (Gauss, Euler, Legendre in drie streken), waar het differentiëren in rekenwerk zou verzuipen. Een volleerd analyticus draagt beide.

23. Met t=sin2θt = \sin^2\theta is  ⁣dt=2sinθcosθ ⁣dθ=2t1/2(1t)1/2 ⁣dθ\dd t = 2\sin\theta\cos\theta\,\dd\theta = 2\,t^{1/2}(1-t)^{1/2}\, \dd\theta, dus

Wp=01tp/2 ⁣dt2t1/2(1t)1/2=1201tp+121(1t)121 ⁣dt=12B(p+12,12).W_p = \int_0^1 t^{p/2}\, \frac{\dd t}{2\,t^{1/2}(1-t)^{1/2}} = \frac12\int_0^1 t^{\frac{p+1}2 - 1}(1-t)^{\frac12 - 1}\dd t = \frac12\,B\Bigl(\frac{p+1}2, \frac12\Bigr).

De bètarecursie met x=n+12x = \frac{n+1}2 en y=12y = \frac12 geeft

Wn+2=12(n+1)/2(n+2)/2B(n+12,12)=n+1n+2Wn.W_{n+2} = \frac12\, \frac{(n+1)/2}{(n+2)/2}\,B\Bigl(\frac{n+1}2, \frac12\Bigr) = \frac{n+1}{n+2}\,W_n .

Vanaf W0=π2W_0 = \frac\pi2 en W1=1W_1 = 1:

W2n=π2k=1n2k12k,W2n+1=k=1n2k2k+1.W_{2n} = \frac\pi2\prod_{k=1}^n\frac{2k-1}{2k}, \qquad W_{2n+1} = \prod_{k=1}^n\frac{2k}{2k+1} .

Omdat sinn+1sinn\sin^{n+1} \leq \sin^n op [0,π/2]\intcc0{\pi/2}, is de rij (Wn)(W_n) niet-stijgend, zodat

1W2n+1W2nW2n+1W2n1=2n2n+11.1 \geq \frac{W_{2n+1}}{W_{2n}} \geq \frac{W_{2n+1}}{W_{2n-1}} = \frac{2n}{2n+1} \longrightarrow 1 .

Maar de gesloten vormen geven

W2n+1W2n=2πk=1n(2k)2(2k1)(2k+1)=2πk=1n4k24k21,\frac{W_{2n+1}}{W_{2n}} = \frac2\pi \prod_{k=1}^n\frac{(2k)^2}{(2k-1)(2k+1)} = \frac2\pi\prod_{k=1}^n\frac{4k^2}{4k^2-1},

en nn \to \infty levert het product van Wallis. (Via de formule van Euler is Wp=π2Γ(p+12)/Γ(p2+1)W_p = \frac{\sqrt\pi}2\, \Gamma(\frac{p+1}2)/\Gamma(\frac p2 + 1): Wallis is de gaussische integraal in weer een ander kostuum.)

24. De afgeleide van de integrand naar bb is 2xex2sin(2bx)-2x\,\eu^{-x^2}\sin(2bx), gedomineerd door 2xex2L1(R)2\abs x\,\eu^{-x^2} \in L^1(\R), uniform in bb: dus is Φ\Phi C1\mathcal C^1 met

Φ(b)=+2xex2sin(2bx) ⁣dx.\Phi'(b) = -\int_{-\infty}^{+\infty} 2x\,\eu^{-x^2}\sin(2bx)\,\dd x .

Partiële integratie met u=sin(2bx)u = \sin(2bx) en  ⁣dv=2xex2 ⁣dx\dd v = -2x\,\eu^{-x^2}\dd x (dus v=ex2v = \eu^{-x^2}) laat de randtermen verdwijnen en geeft

Φ(b)=2b+ex2cos(2bx) ⁣dx=2bΦ(b).\Phi'(b) = -2b\int_{-\infty}^{+\infty} \eu^{-x^2}\cos(2bx)\,\dd x = -2b\,\Phi(b).

Bijgevolg is (Φ(b)eb2)=0\bigl(\Phi(b)\,\eu^{b^2}\bigr)' = 0 en Φ(b)=Φ(0)eb2=πeb2\Phi(b) = \Phi(0)\,\eu^{-b^2} = \sqrt\pi\,\eu^{-b^2} volgens Deel I. Op de normering na zegt dat dat de fouriergetransformeerde van ex2\eu^{-x^2} weer een gaussische functie is — het vaste punt waarop de inversietheorie van Hoofdstuk 14 draait.

25. Voor 0<a<b0 < a < b en x>0x > 0 is

0eaxebxx=abesx ⁣ds(ba)eax,0 \leq \frac{\eu^{-ax} - \eu^{-bx}}{x} = \int_a^b \eu^{-sx}\,\dd s \leq (b - a)\,\eu^{-ax},

dus convergeert de integraal I(a)I(a) (in de zin van Lebesgue); en de puntsgewijze grens laat ook zien dat de integrand zich in x=0x = 0 continu voortzet met de waarde bab - a. Leg bb vast; op [a0,+)\intco{a_0}{+\infty} is de afgeleide van de integrand naar aa gelijk aan eax-\eu^{-ax}, gedomineerd door ea0xL1((0,+))\eu^{-a_0x} \in L^1(\intoo0{+\infty}), dus is II C1\mathcal C^1 op (0,b)\intoo0b met

I(a)=0+eax ⁣dx=1a,en dusI(a)=loga+c.I'(a) = -\int_0^{+\infty}\eu^{-ax}\dd x = -\frac1a, \qquad\text{en dus}\qquad I(a) = -\log a + c .

De tweezijdige grens geeft 0I(a)(ba)/a00 \leq I(a) \leq (b-a)/a \to 0 als aba \to b^-, dus c=logbc = \log b en I(a)=logbaI(a) = \log\frac ba. De ophefbare singulariteit is in 00 nodig: elke term eax/x\eu^{-ax}/x heeft afzonderlijk een divergente (logaritmische) integraal nabij 00, en alleen de opheffing van de eerste orde eaxebx=O(x)\eu^{-ax} - \eu^{-bx} = O(x) maakt het verschil daar integreerbaar; in het oneindige is elke term op zichzelf al onschadelijk.

Begrippen gedefinieerd in dit hoofdstuk

Bekijk alle 395 begrippen in de begrippenlijst