Mathematics · Boek 5 · Bachelor Year 3

Universitaire wiskunde — Bachelor jaar 3

Universitaire wiskunde — Bachelor jaar 3 · Bachelor Year 3

7Volledige ruimten: Baire, Ascoli, Stone–Weierstrass

Volledigheid — elke Cauchyrij convergeert — is de eigenschap waarmee de analyse objecten voortbrengt: vaste punten van contracties, sommen van reeksen, oplossingen van vergelijkingen verkregen als limiet. Dit hoofdstuk brengt de drie grote bestaansmachines van de metrische theorie samen. De stelling van Baire laat zien dat een volledige ruimte geen aftelbare vereniging van verwaarloosbare stukken kan zijn, en tovert objecten tevoorschijn (continue nergens differentieerbare functies!) uit een redenering die niet veel meer dan tellen is. Arzelà–Ascoli bepaalt de compacte deelverzamelingen van C(K)\mathcal C(K) en is het compactheidswerkpaard van de analyse — de weekendopgave bewijst er de bestaansstelling van Peano voor differentiaalvergelijkingen mee. Stone–Weierstrass toont aan dat veeltermen, en veel meer, dicht liggen in C(K)\mathcal C(K): benaderen wordt een algebraïsche controle. Onderweg construeren we vervollediging en bewijzen we de voortzettingsstelling voor uniform continue afbeeldingen, het dagelijks brood van Hoofdstukken 12, 13 en 14.

7.1 Volledige ruimten, vervollediging, voortzetting

Definitie 7.1

Een metrische ruimte heet volledig als elke Cauchyrij convergeert (bachelorjaar 2: Rn\R^n is volledig; C([0,1])\mathcal C(\intcc01) met \norm\cdot_\infty is volledig). Een gesloten deelverzameling van een volledige ruimte is volledig; en een volledige deelverzameling van een willekeurige metrische ruimte is gesloten.

Bewijs. Voor de twee uitspraken: een Cauchyrij uit de gesloten FF convergeert in XX, en haar limiet, adherent aan FF, ligt in FF; en een in XX convergerende rij uit een volledige AA is Cauchy, convergeert dus in AA, en limieten zijn uniek.

Stelling 7.2 (Voortzetting van uniform continue afbeeldingen)

Zij DXD \subseteq X dicht, YY volledig, en f ⁣:DYf \colon D \to Y uniform continu. Dan zet ff zich op precies één manier voort tot een continue fˉ ⁣:XY\bar f \colon X \to Y, en fˉ\bar f is uniform continu.

Bewijs. Eenduidigheid: twee continue voortzettingen vallen samen op de dichte DD en dus overal (de verzameling {g=h}\{g = h\} is gesloten: het origineel van de gesloten diagonaal onder x(g(x),h(x))x\mapsto(g(x), h(x))). Bestaan: kies voor xXx \in X een rij dnxd_n \to x met dnDd_n \in D. De rij (f(dn))(f(d_n)) is Cauchy: bij gegeven ε\varepsilon levert de uniforme continuïteit een δ\delta met d(u,v)<δd(f(u),f(v))<εd(u,v) < \delta \Rightarrow d(f(u), f(v)) < \varepsilon, en (dn)(d_n) is Cauchy. Definieer fˉ(x)=limf(dn)\bar f(x) = \lim f(d_n); de limiet hangt niet van de gekozen rij af (verweef er twee). fˉ\bar f zet ff voort (constante rijen) en erft de continuïteitsmodulus: is d(x,x)<δd(x, x') < \delta, dan geeft benadering van beide door punten van DD op afstand <δd(x,x)2< \frac{\delta - d(x,x')}2 in de limiet d(fˉ(x),fˉ(x))εd(\bar f(x), \bar f(x')) \leq \varepsilonfˉ\bar f is uniform continu.

Stelling 7.3 (Vervollediging)

Elke metrische ruimte XX bedt zich isometrisch in als dichte deelverzameling van een volledige metrische ruimte X^\hat X, uniek op een isometrie na die XX puntsgewijs vasthoudt: haar vervollediging.

Bewijs. Bestaan. Zij X\mathcal X de verzameling Cauchyrijen van XX, met de pseudo-afstand

D((xn),(yn))=limnd(xn,yn),D\bigl((x_n), (y_n)\bigr) = \lim_n d(x_n, y_n),

waarbij de limiet bestaat omdat d(xn,yn)d(xm,ym)d(xn,xm)+d(yn,ym)\abs{d(x_n, y_n) - d(x_m, y_m)} \leq d(x_n, x_m) + d(y_n, y_m) de reële rij Cauchy maakt. Stel X^=X/\hat X = \mathcal X/{\sim}, waarbij rijen op DD-afstand 00 vereenzelvigd worden; DD daalt af tot een afstand. Bed XX in via de constante rijen: een isometrie met dicht beeld (een Cauchyrij wordt in DD benaderd door de constanten gebouwd op haar eigen termen: D((xn),(xk)const)=limnd(xn,xk)0D\bigl((x_n), (x_k)_{\rm const}\bigr) = \lim_n d(x_n, x_k) \to 0 als kk \to \infty, wegens de Cauchy-eigenschap). Volledigheid van X^\hat X: zij (ξk)(\xi^k) Cauchy in X^\hat X; kies wegens de dichtheid xkXx_k \in X met D(ξk,xk)2kD(\xi^k, x_k) \leq 2^{-k}; dan is (xk)(x_k) Cauchy in XX (driehoeksongelijkheid via de ξ\xi), definieert het een punt ξX^\xi \in \hat X, en is D(ξk,ξ)2k+D(xk,ξ)0D(\xi^k, \xi) \leq 2^{-k} + D(x_k, \xi) \to 0 (de afstand van de constante xkx_k tot de klasse van (xj)j(x_j)_j is limjd(xk,xj)\lim_j d(x_k, x_j), klein voor grote kk).

Eenduidigheid: twee vervolledigingen X^1\hat X_1 en X^2\hat X_2 bevatten XX dicht; de identiteit van XX, een isometrie, is uniform continu en zet zich dus voort tot X^1X^2\hat X_1 \to \hat X_2 (Stelling 7.2), nog steeds een isometrie op een dichte verzameling en dus overal; symmetrisch in de andere richting, en de samenstellingen houden de dichte XX vast: het zijn de identiteiten.

Stelling 7.4 (Vast punt van Banach)

Zij XX volledig en niet-leeg, en f ⁣:XXf \colon X \to X een contractie: d(f(x),f(y))kd(x,y)d(f(x), f(y)) \leq k\,d(x,y) met k<1k < 1. Dan heeft ff precies één vast punt xx^*, en convergeert elke baan ernaartoe, met de expliciete snelheid d(xn,x)kn1kd(x1,x0)d(x_n, x^*) \leq \frac{k^n}{1-k}\,d(x_1, x_0).

Bewijs. (Bachelorjaar 2 bewees dit; we noteren het argument van twee regels opnieuw, zodat dit hoofdstuk op zichzelf staat.) De baan xn+1=f(xn)x_{n+1} = f(x_n) voldoet aan d(xn+1,xn)knd(x1,x0)d(x_{n+1}, x_n) \leq k^nd(x_1, x_0) en is dus Cauchy (meetkundige reeks); haar limiet xx^* is vast (continuïteit van ff) en uniek, want twee vaste punten voldoen aan dkdd \leq k\,d. De snelheid: sommeer de meetkundige staart.

Voorbeeld 7.5 (De identiteit verstoren)

Zij g ⁣:RdRdg \colon \R^d \to \R^d kk-lipschitz met k<1k < 1. Dan is φ=id+g\varphi = \mathrm{id} + g een homeomorfisme van Rd\R^d op Rd\R^d. Injectiviteit, met een kwantitatieve modulus:

φ(x)φ(y)xyg(x)g(y)(1k)xy.\norm{\varphi(x) - \varphi(y)} \geq \norm{x - y} - \norm{g(x) - g(y)} \geq (1 - k)\norm{x - y} .

Surjectiviteit is de vastepuntstelling: φ(x)=y\varphi(x) = y oplossen betekent x=yg(x)x = y - g(x), en xyg(x)x \mapsto y - g(x) is een kk-contractie van het volledige Rd\R^d — voor elke yy bestaat er precies één oplossing x=ψ(y)x = \psi(y). De getoonde ongelijkheid maakt de inverse ψ\psi lipschitz met constante 11k\frac1{1-k}: een homeomorfisme, met expliciete grenzen op beide moduli. Deze onschuldig ogende uitspraak is de motor binnen de inverse-functiestelling (Hoofdstuk 20): nabij een punt waar DfDf inverteerbaar is, is ff een inverteerbare lineaire afbeelding plus een kleine lipschitz-verstoring, en het voorbeeld van vandaag doet de rest. Het kwantificeert ook de numerieke robuustheid: een stelsel dat minder wordt verstoord dan de marge van de inverse, blijft oplosbaar, en de oplossing verschuift met hoogstens 11k\frac{1}{1-k} maal de verstoring.

7.2 De stelling van Baire

Stelling 7.6 (Baire)

In een volledige metrische ruimte is een aftelbare doorsnede van dichte open verzamelingen dicht. Gelijkwaardig: is X=nFnX = \bigcup_{n}F_n met elke FnF_n gesloten, dan heeft een zekere FnF_n een niet-leeg inwendige.

Bewijs. Zij (Un)(U_n) dicht en open, en B0=B(x0,r0)B_0 = B(x_0, r_0) een willekeurige open bal; we zoeken een punt van Un\bigcap U_n in B0B_0. Inductief: omdat UnU_{n} dicht en open is, snijdt het de open bal Bn1B_{n-1} in een open verzameling, die een gesloten bal Bˉ(xn,rn)\bar B(x_n, r_n) bevat met 0<rnrn1/20 < r_n \leq r_{n-1}/2 en Bˉ(xn,rn)Bn1Un\bar B(x_n, r_n) \subseteq B_{n-1}\cap U_n. De middelpunten vormen een Cauchyrij (xmBnx_m \in B_n voor mnm \geq n, en de stralen gaan naar 00); de limiet xx ligt in elke Bˉ(xn,rn)\bar B(x_n, r_n) (geslotenheid), en dus in elke UnU_n en in B0B_0. Voor de tweede vorm: heeft geen enkele FnF_n een inwendige, dan zijn de Un=XFnU_n = X \setminus F_n open en dicht, en ontsnapt een punt van Un\bigcap U_n aan Fn=X\bigcup F_n = X: absurd.

Opmerking 7.7

Woordenschat: een verzameling heet nergens dicht als haar afsluiting een leeg inwendige heeft, en mager (van de eerste categorie) als ze een aftelbare vereniging van nergens dichte verzamelingen is. Baire zegt: een volledige metrische ruimte is niet mager in zichzelf, en het complement van een magere verzameling is dicht. “Mager” is een begrip van kleinheid dat loodrecht op de maat staat (Hoofdstuk 9 zal magere verzamelingen van volle maat opleveren), en baireargumenten bewijzen bestaan door overvloed: om één object zonder eigenschap P aan te wijzen, toon je aan dat de objecten mét P een magere verzameling vormen.

Gevolg 7.8

(a) R\R is overaftelbaar. (b) Q\Q is geen aftelbare doorsnede van open deelverzamelingen van R\R, en een niet-lege volledige metrische ruimte zonder geïsoleerde punten is overaftelbaar.

Bewijs. (a) Was R=x{x}\R = \bigcup_{x}\{x\} over een aftelbare verzameling, dan had een zeker singleton een inwendige. (b) Was Q=nVn\Q = \bigcap_n V_n met VnV_n open (noodzakelijk dicht, want Q\supseteq \Q), dan vormen de VnV_n en de complementen R{q}\R\setminus\{q\} met qQq \in \Q een aftelbare familie dichte open verzamelingen met lege doorsnede — in strijd met Baire. Is XX volledig zonder geïsoleerde punten en aftelbaar, dan toont X=xX{x}X = \bigcup_{x \in X}\{x\} haar als een aftelbare vereniging gesloten verzamelingen met leeg inwendige (geen geïsoleerde punten): opnieuw Baire.

Stelling 7.9 (De monsters van Weierstrass bestaan)

Er bestaan continue functies op [0,1]\intcc01 die in geen enkel punt differentieerbaar zijn. Sterker: de verzameling fC([0,1])f \in \mathcal C(\intcc01) die zelfs maar in één punt een (eindige) afgeleide hebben, is mager in (C([0,1]),)(\mathcal C(\intcc01), \norm\cdot_\infty).

Bewijs. Stel voor n1n \geq 1

Fn={f:x[0,1], h0 met x+h[0,1], f(x+h)f(x)nh}.F_n = \Bigl\{f : \exists x \in \intcc01,\ \forall h \neq 0 \text{ met } x + h \in \intcc01,\ \abs{f(x+h) - f(x)} \leq n\abs h \Bigr\}.

Is ff differentieerbaar in xx, dan is fFnf \in F_n voor zekere nn: het quotiënt f(x+h)f(x)/h\abs{f(x+h)-f(x)}/\abs h is begrensd voor hδ\abs h \leq \delta (differentieerbaarheid: het gaat naar f(x)\abs{f'(x)}) en begrensd door 2f/δ2\norm f_\infty/\delta voor hδ\abs h \geq \delta. Dus bevat Fn\bigcup F_n alle functies die ergens differentieerbaar zijn, en volstaat het aan te tonen dat elke FnF_n gesloten is met leeg inwendige.

Gesloten: zij fkff_k \to f uniform met fkFnf_k \in F_n en getuigen xkxx_k \to x (compactheid, na extractie). Voor hh met x+h[0,1]x + h \in \intcc01: kies hkhh_k \to h met xk+hk[0,1]x_k + h_k \in \intcc01 (bijvoorbeeld hk=h+xxkh_k = h + x - x_k, afgeknipt); dan is f(x+h)f(x)=limfk(xk+hk)fk(xk)limnhk=nh\abs{f(x + h) - f(x)} = \lim \abs{f_k(x_k + h_k) - f_k(x_k)} \leq \lim n\abs{h_k} = n\abs h, met de uniforme convergentie en de continuïteit van ff in de betrokken punten: dus fFnf \in F_n.

Leeg inwendige: gegeven fFnf \in F_n en ε>0\varepsilon > 0 zoeken we een gg met gfε\norm{g - f}_\infty \leq \varepsilon en gFng \notin F_n. Benader eerst ff tot op ε/2\varepsilon/2 door een stuksgewijs affiene φ\varphi (uniforme continuïteit: interpoleer op een fijn rooster), met hellingen begrensd door zekere MM. Tel er een kleine zaagtand bij op: g=φ+ε2sNg = \varphi + \frac\varepsilon2\,s_N, waarbij sN(x)s_N(x) de 1N\frac1N-periodieke zigzag met amplitude 11 en helling ±2N\pm 2N is. In elke xx is er aan één kant een willekeurig kleine hh waarbij de zaagtand op [x,x+h][x, x+h] een helling ±2Nε2=±εN\pm 2N\cdot\frac\varepsilon2 = \pm\varepsilon N bijdraagt: het differentiequotiënt van gg overschrijdt εNM>n\varepsilon N - M > n voor grote NN. Dus is gFng \notin F_n, op willekeurig kleine uniforme afstand ε\leq \varepsilon van ff. Besluit: Fn\bigcup F_n is mager, en volgens Baire is haar complement — bestaande uit nergens differentieerbare functiesdicht in C([0,1])\mathcal C(\intcc01): zulke functies bestaan in overvloed.

7.3 Arzelà–Ascoli

Overal is KK een compacte metrische ruimte en C(K)=C(K,Rd)\mathcal C(K) = \mathcal C(K, \R^d), met f=supKf(x)\norm f_\infty = \sup_K \norm{f(x)}: een volledige ruimte (uniforme limieten van continue functies zijn continu — bachelorjaar 2).

Definitie 7.10

Een familie FC(K)\mathcal F \subseteq \mathcal C(K) heet equicontinu als er voor elke ε>0\varepsilon > 0 een δ>0\delta > 0 is zodanig dat

d(x,y)<δ    f(x)f(y)<εvoor alle fFd(x, y) < \delta \implies \norm{f(x) - f(y)} < \varepsilon \quad\text{voor \emph{alle} } f \in \mathcal F

(één δ\delta voor de hele familie — bijvoorbeeld elke familie met een gemeenschappelijke lipschitzconstante of een gemeenschappelijke höldermodulus), en puntsgewijs begrensd als supff(x)<\sup_{f}\norm{f(x)} < \infty voor elke xx.

Stelling 7.11 (Arzelà–Ascoli)

Een deelverzameling FC(K)\mathcal F \subseteq \mathcal C(K) is relatief compact (heeft compacte afsluiting) dan en slechts dan als ze equicontinu en puntsgewijs begrensd is. In het bijzonder heeft elke equicontinue, puntsgewijs begrensde rij een uniform convergente deelrij.

Bewijs. (\Leftarrow) Zij (fk)(f_k) een rij uit F\mathcal F. Een compacte metrische KK is separabel: voor elke nn overdekken eindig veel ballen met straal 1n\frac1n de ruimte KK (Stelling 6.16); hun middelpunten vormen een aftelbare dichte verzameling D={x1,x2,}D = \{x_1, x_2, \dots\}. Haal met de puntsgewijze begrensdheid en Bolzano–Weierstrass achtereenvolgens deelrijen tevoorschijn die in x1x_1 convergeren, daarna ook in x2x_2, enzovoort, en neem de diagonaaldeelrij (gj)(g_j): die convergeert in elk punt van DD. De equicontinuïteit tilt dat op tot uniform Cauchy: kies bij gegeven ε\varepsilon een δ\delta als in de definitie, overdek KK met eindig veel ballen B(xi,δ)B(x_i, \delta) met xiDx_i \in D (imi \leq m), en kies JJ zo groot dat gj(xi)gj(xi)<ε\norm{g_j(x_i) - g_{j'}(x_i)} < \varepsilon voor j,jJj, j' \geq J en imi \leq m. Voor willekeurige xB(xi,δ)x \in B(x_i, \delta):

gj(x)gj(x)gj(x)gj(xi)+gj(xi)gj(xi)+gj(xi)gj(x)<3ε.\norm{g_j(x) - g_{j'}(x)} \leq \norm{g_j(x) - g_j(x_i)} + \norm{g_j(x_i) - g_{j'}(x_i)} + \norm{g_{j'}(x_i) - g_{j'}(x)} < 3\varepsilon .

Dus is (gj)(g_j) uniform Cauchy en convergeert ze in het volledige C(K)\mathcal C(K). Elke rij uit F\mathcal F heeft dus een convergente deelrij: Fˉ\bar{\mathcal F} is (rijcompact, en dus volgens Stelling 6.16) compact.

(\Rightarrow) Is Fˉ\bar{\mathcal F} compact, dan is de puntsgewijze begrensdheid duidelijk (evalueren is continu). Voor de equicontinuïteit: overdek F\mathcal F met eindig veel ballen B(fi,ε)B(f_i, \varepsilon) van C(K)\mathcal C(K); elke fif_i is uniform continu (Heine, Gevolg 6.17), wat een gemeenschappelijke δ\delta geeft voor imi \leq m; en dan is voor fB(fi,ε)f \in B(f_i, \varepsilon) en d(x,y)<δd(x,y) < \delta: f(x)f(y)2ε+fi(x)fi(y)<3ε\norm{f(x)-f(y)} \leq 2\varepsilon + \norm{f_i(x)-f_i(y)} < 3\varepsilon.

Voorbeeld 7.12

De gesloten eenheidsbal van C([0,1])\mathcal C(\intcc01) is niet compact (fn(x)=xnf_n(x) = x^n heeft geen uniform convergente deelrij: de puntsgewijze limiet is discontinu), en inderdaad is (xn)(x^n) niet equicontinu in 11. Daartegenover is {f:f1, Lip(f)1}\{f : \norm f_\infty \leq 1,\ \operatorname{Lip}(f) \leq 1\} wél compact: begrensd, equicontinu met lipschitzconstante 11, en gesloten. Ascoli verklaart waarom de compactheid in oneindige dimensie faalt (Riesz, bachelorjaar 2) en wat je moet toevoegen om haar te herstellen: een uniforme continuïteitsmodulus.

7.4 Stone–Weierstrass

Lemma 7.13 (Dini)

Zij KK compact en (fn)(f_n) een monotone rij continue reële functies die puntsgewijs naar een continue ff convergeert. Dan is de convergentie uniform.

Bewijs. Zeg fnff_n \uparrow f; stel gn=ffn0g_n = f - f_n \downarrow 0, continu. Bij gegeven ε\varepsilon stijgen de open verzamelingen Un={gn<ε}U_n = \{g_n < \varepsilon\} en overdekken ze KK (puntsgewijze convergentie); haal er een eindige deeloverdekking uit: K=Un0K = U_{n_0} voor zekere n0n_0 (stijgende familie), dat wil zeggen 0gn<ε0 \leq g_n < \varepsilon overal voor nn0n \geq n_0.

Lemma 7.14

Er is een rij veeltermen unu_n met un(t)tu_n(t) \to \sqrt t uniform op [0,1]\intcc01.

Bewijs. Definieer u0=0u_0 = 0 en un+1(t)=un(t)+12(tun(t)2)u_{n+1}(t) = u_n(t) + \frac12\bigl(t - u_n(t)^2\bigr): veeltermen. Per inductie is 0un(t)t0 \leq u_n(t) \leq \sqrt t op [0,1]\intcc01: aangenomen voor nn, is

tun+1(t)=(tun(t))(112(t+un(t)))0,\sqrt t - u_{n+1}(t) = \bigl(\sqrt t - u_n(t)\bigr) \Bigl(1 - \tfrac12\bigl(\sqrt t + u_n(t)\bigr)\Bigr) \geq 0,

want t+un2\sqrt t + u_n \leq 2; en un+1un0u_{n+1} \geq u_n \geq 0. Dus is (un(t))(u_n(t)) niet-dalend en begrensd door t\sqrt t: ze convergeert puntsgewijs, en de limiet (t)\ell(t) voldoet aan =+12(t2)\ell = \ell + \frac12(t - \ell^2), dus (t)=t\ell(t) = \sqrt t, continu. Dini (Lemma 7.13) tilt dat op tot uniform.

Stelling 7.15 (Stone–Weierstrass, reële versie)

Zij KK een compacte (hausdorffse) ruimte en AC(K,R)\mathcal A \subseteq \mathcal C(K, \R) een deelalgebra (gesloten onder sommen, producten en scalaire veelvouden) die de constanten bevat en de punten scheidt (voor xyx \neq y is er een fAf \in \mathcal A met f(x)f(y)f(x) \neq f(y)). Dan ligt A\mathcal A dicht in (C(K,R),)(\mathcal C(K,\R), \norm\cdot_\infty).

Bewijs. Zij Aˉ\bar{\mathcal A} de afsluiting, opnieuw een algebra (producten van uniforme limieten convergeren op begrensde verzamelingen).

Stap 1: Aˉ\bar{\mathcal A} is een tralie, dat wil zeggen gesloten onder max\max en min\min. Omdat max(f,g)=f+g+fg2\max(f,g) = \frac{f + g + \abs{f-g}}2 en net zo voor min\min, volstaat het dat uit fAˉf \in \bar{\mathcal A} volgt fAˉ\abs f \in \bar{\mathcal A}: met M=f>0M = \norm f_\infty > 0 is f=M(f/M)2\abs f = M\sqrt{(f/M)^2}, en Lemma 7.14 geeft veeltermen unu_n met un((f/M)2)f/Mu_n\bigl((f/M)^2\bigr) \to \abs f/M uniform; en veeltermen in leden van de algebra (met constante term: de constanten zitten erin) blijven in Aˉ\bar{\mathcal A}.

Stap 2: interpolatie in twee punten. Voor xyx \neq y en a,bRa, b \in \R is er een gAg \in \mathcal A met g(x)=ag(x) = a en g(y)=bg(y) = b: neem een hh die xx en yy scheidt en stel g=a+(ba)hh(x)h(y)h(x)g = a + (b - a)\frac{h - h(x)}{h(y) - h(x)}.

Stap 3. Zij fC(K)f \in \mathcal C(K) en ε>0\varepsilon > 0. Kies voor elk paar x,yx, y een gx,yAg_{x,y} \in \mathcal A met gx,y(x)=f(x)g_{x,y}(x) = f(x) en gx,y(y)=f(y)g_{x,y}(y) = f(y) (stap 2; voor x=yx = y neem je de constante functie gx,x=f(x)g_{x,x} = f(x)). Leg xx vast: voor elke yy bevat de open verzameling Vy={gx,y<f+ε}V_y = \{g_{x,y} < f + \varepsilon\} het punt yy; de compactheid levert y1,,ymy_1, \dots, y_m met K=VyjK = \bigcup V_{y_j}, en hx=minjgx,yjAˉh_x = \min_j g_{x, y_j} \in \bar{\mathcal A} (stap 1) voldoet overal aan hx<f+εh_x < f + \varepsilon, met hx(x)=f(x)h_x(x) = f(x). Varieer nu xx: Wx={hx>fε}W_x = \{h_x > f - \varepsilon\} is open en bevat xx; haal er x1,,xlx_1, \dots, x_l uit die KK overdekken, en h=maxihxiAˉh = \max_i h_{x_i} \in \bar{\mathcal A} voldoet aan fε<h<f+εf - \varepsilon < h < f + \varepsilon: dus fhε\norm{f - h}_\infty \leq \varepsilon. Bijgevolg is fAˉf \in \bar{\mathcal A}.

Gevolg 7.16

(a) (Weierstrass) De veeltermen liggen dicht in C([a,b],R)\mathcal C(\intcc ab, \R); en de veeltermen in nn veranderlijken liggen dicht in C(K,R)\mathcal C(K, \R) voor compacte KRnK \subseteq \R^n. (b) (Complexe versie) Is AC(K,C)\mathcal A \subseteq \mathcal C(K, \C) een deelalgebra die de constanten bevat, de punten scheidt en gesloten is onder toevoeging, dan ligt ze dicht. (c) (Trigonometrische versie) De trigonometrische veeltermen nNcneint\sum_{\abs n \leq N}c_n\eu^{\iu n t} liggen dicht in de ruimte van continue 2π2\pi-periodieke functies met \norm\cdot_\infty.

Bewijs. (a) De veeltermen vormen een algebra met de constanten; en de coördinaatfuncties scheiden de punten van Rn\R^n. (b) De reële en imaginaire delen f+fˉ2\frac{f + \bar f}2 en ffˉ2i\frac{f - \bar f}{2\iu} van leden van A\mathcal A vormen een reële algebra ARC(K,R)\mathcal A_\R \subseteq \mathcal C(K, \R) met de constanten; ze scheidt de punten (uit f(x)f(y)f(x) \ne f(y) volgt dat Ref\operatorname{Re}f of Imf\operatorname{Im}f scheidt). Pas de reële stelling toe en zet weer samen. (c) Beschouw de 2π2\pi-periodieke functies als C(S1,C)\mathcal C(S^1, \C) (met S1=R/2πZS^1 = \R/2\pi\Z, compact: Oefening 6.5); de algebra voortgebracht door eit\eu^{\iu t}, eit\eu^{-\iu t} en de constanten is gesloten onder toevoeging en scheidt de punten van de cirkel (eit\eu^{\iu t} is er injectief op). Pas (b) toe.

Opmerking 7.17

De trigonometrische versie repareert, en veralgemeent enorm, het gat dat het fourierhoofdstuk van bachelorjaar 2 openliet: de dichtheid van de trigonometrische veeltermen in (C(S1),2)(\mathcal C(S^1), \norm\cdot_2) volgt er a fortiori uit (2\norm\cdot_2 \leq \norm\cdot_\infty op de normeringsconstante na), wat van het fouriersysteem in Hoofdstuk 13 een orthonormale basis zal maken en Parseval eindelijk in volle algemeenheid zal bewijzen.

7.5 Oefeningen

Oefening 7.1

(a) Toon aan dat C([0,1],R)\mathcal C(\intcc01, \R) met f1=01f\norm f_1 = \int_0^1\abs f niet volledig is: de functies fnf_n die gelijk zijn aan 00 op [0,12]\intcc0{\frac12} en aan 11 op [12+1n,1]\intcc{\frac12 + \frac1n}1, en daartussen affien, zijn 1\norm\cdot_1-Cauchy zonder continue limiet. (b) Toon aan dat een genormeerde ruimte waarin elke absoluut convergente reeks convergeert, volledig is. (Haal uit een Cauchyrij een deelrij met xnk+1xnk2k\norm{x_{n_{k+1}} - x_{n_k}} \leq 2^{-k}.)

Oplossing

Oplossing van Oefening 7.1.

(a) Voor mnm \geq n verdwijnt fnfmf_n - f_m buiten een interval van lengte 1n\frac1n en is ze begrensd door 11: fnfm11n0\norm{f_n - f_m}_1 \leq \frac1n \to 0, dus Cauchy. Was fnff_n \to f in 1\norm\cdot_1 met ff continu, dan is voor een vaste α<12\alpha < \frac12: 0αf=lim0αffn+fnlim(ffn1+0)=0\int_0^{\alpha}\abs f = \lim \int_0^\alpha\abs{f - f_n + f_n} \leq \lim\bigl(\norm{f - f_n}_1 + 0\bigr) = 0 (want fn0f_n \equiv 0 daar voor grote nn), dus f0f \equiv 0 op [0,12)[0, \frac12) (continuïteit); en net zo f1f \equiv 1 op (12,1](\frac12, 1]: geen continue functie doet beide. De ruimte is dus onvolledig — de vervollediging is L1L^1, gebouwd in Hoofdstuk 12.

(b) Zij (xn)(x_n) Cauchy; kies n1<n2<n_1 < n_2 < \cdots met xnk+1xnk2k\norm{x_{n_{k+1}} - x_{n_k}} \leq 2^{-k}. De reeks k(xnk+1xnk)\sum_k (x_{n_{k+1}} - x_{n_k}) convergeert absoluut en dus; haar partiële sommen zijn xnk+1xn1x_{n_{k+1}} - x_{n_1}, dus convergeert (xnk)(x_{n_k}), en een Cauchyrij met een convergente deelrij convergeert.

Oefening 7.2

Met Baire: (a) toon aan dat een volledige genormeerde ruimte geen aftelbare (algebraïsche) basis heeft — leid af dat de ruimte van de veeltermen voor geen enkele norm volledig is; (b) toon aan dat als een rij continue fn ⁣:RRf_n \colon \R \to \R puntsgewijs naar ff convergeert, de verzameling continuïteitspunten van ff dicht is. (Voor (b): neem aan of bewijs dat Ωδ={x:oscxf<δ}\Omega_\delta = \{x : \operatorname{osc}_x f < \delta\} open is, en toon aan dat ze dicht is met behulp van FN={x:fn(x)fm(x)δ/3 n,mN}F_{N} = \{x: \abs{f_n(x) - f_m(x)} \leq \delta/3\ \forall n,m \geq N\}, gesloten verzamelingen die R\R overdekken; werk in een willekeurige gesloten bal om Baire toe te passen.)

Oplossing

Oplossing van Oefening 7.2.

(a) Stel dat EE volledig is met algebraïsche basis (en)nN(e_n)_{n\in\N}, en stel Fn=Vect(e1,,en)F_n = \operatorname{Vect}(e_1, \dots, e_n): gesloten (eindigdimensionale deelruimten zijn volledig en dus gesloten — bachelorjaar 2), met leeg inwendige: was B(x,r)FnB(x, r) \subseteq F_n, neem dan vFnv \notin F_n; dan is x+r2vvB(x,r)Fnx + \frac{r}{2\norm v}v \in B(x,r) \setminus F_n, absurd. Maar E=FnE = \bigcup F_n (elke vector is een eindige combinatie): in strijd met Baire (Stelling 7.6). De ruimte R[X]\R[X] heeft de aftelbare basis (Xn)(X^n), dus maakt geen enkele norm haar volledig.

(b) Leg δ>0\delta > 0 en een niet-lege open bal BB vast; we zoeken in BB een punt van Ωδ={x:oscxf<δ}\Omega_\delta = \{x : \operatorname{osc}_xf < \delta\}, waarbij oscxf=infVxdiamf(V)\operatorname{osc}_xf = \inf_{V \ni x}\operatorname{diam} f(V). (Ωδ\Omega_\delta is open: is diamf(V)<δ\operatorname{diam}f(V) < \delta voor een open VxV \ni x, dan heeft elke yVy \in V oscillatie <δ< \delta.) De verzamelingen

FN={x:fn(x)fm(x)δ3  n,mN}F_N = \bigl\{x : \abs{f_n(x) - f_m(x)} \leq \tfrac\delta3\ \ \forall\, n, m \geq N \bigr\}

zijn gesloten (doorsneden van originelen van gesloten verzamelingen) en overdekken R\R (de puntsgewijze convergentie maakt (fn(x))(f_n(x)) Cauchy). Baire toegepast binnen de volledige Bˉ\bar B: een zekere FNBˉF_N \cap \bar B bevat een bal B=B(x0,ρ)B' = B(x_0, \rho). Laat mm \to \infty gaan: fNfδ3\abs{f_N - f} \leq \frac\delta3 op BB'. Krimp wegens de continuïteit van fNf_N in x0x_0 tot een Bx0B'' \ni x_0 waar fNfN(x0)δ3\abs{f_N - f_N(x_0)} \leq \frac\delta3; dan is voor xBx \in B''

f(x)f(x0)ffN(x)+fN(x)fN(x0)+fNf(x0)δ:\abs{f(x) - f(x_0)} \leq \abs{f - f_N}(x) + \abs{f_N(x) - f_N(x_0)} + \abs{f_N - f}(x_0) \leq \delta:

dus diamf(B)2δ\operatorname{diam} f(B'') \leq 2\delta en x0Ω3δBx_0 \in \Omega_{3\delta} \cap B. Elke Ωδ\Omega_\delta is dus open en dicht, en kΩ1/k\bigcap_k\Omega_{1/k} — de verzameling continuïteitspunten — is dicht volgens Baire.

Oefening 7.3 ★★

(a) Toon aan dat f(x)=12(x+ax)f(x) = \frac12\bigl(x + \frac ax\bigr) (met a>1a > 1) een contractie van [a,+)[\sqrt a, +\infty) is en bepaal haar vaste punt — de methode van Heron. Schat het aantal iteraties voor een nauwkeurigheid 101210^{-12} vanaf x0=ax_0 = a, voor a=2a = 2. (b) (Vergelijking van Kepler) Toon aan dat voor 0e<10 \leq e < 1 en mRm \in \R de vergelijking x=m+esinxx = m + e\sin x precies één oplossing heeft, die continu van mm afhangt.

Oplossing

Oplossing van Oefening 7.3.

(a) f(x)=12(x+ax)f(x) = \frac12(x + \frac ax) beeldt [a,)[\sqrt a, \infty) op zichzelf af (ongelijkheid van rekenkundig en meetkundig gemiddelde: f(x)xax=af(x) \geq \sqrt{x \cdot \frac ax} = \sqrt a), en daar is f(x)=12(1ax2)[0,12)f'(x) = \frac12(1 - \frac a{x^2}) \in [0, \frac12): een 12\frac12-lipschitz contractie van een gesloten (dus volledige) verzameling. Vast punt: x=f(x)    x2=ax = f(x) \iff x^2 = a, dus x=ax^* = \sqrt a. Snelheid (Stelling 7.4): d(xn,a)2n+1d(x1,x0)d(x_n, \sqrt a) \leq 2^{-n+1}\,d(x_1, x_0). Voor a=2a = 2 en x0=2x_0 = 2 is d(x1,x0)=12d(x_1, x_0) = \frac12, dus garandeert n=40n = 40 dat 240<10122^{-40} < 10^{-12}. (In werkelijkheid convergeert de methode van Newton kwadratisch: een handvol iteraties volstaat; de contractieschatting is pessimistisch maar gratis.)

(b) fm(x)=m+esinxf_m(x) = m + e\sin x is ee-lipschitz met e<1e < 1 op het volledige R\R: dus precies één vast punt x(m)x(m). Voor twee parameters:

x(m)x(m)=fm(x(m))fm(x(m))mm+ex(m)x(m),\abs{x(m) - x(m')} = \abs{f_m(x(m)) - f_{m'}(x(m'))} \leq \abs{m - m'} + e\abs{x(m) - x(m')},

dus x(m)x(m)mm1e\abs{x(m) - x(m')} \leq \frac{\abs{m - m'}}{1 - e}: zelfs lipschitz in mm.

Oefening 7.4 ★★

(a) Zij XX volledig en f ⁣:XXf \colon X \to X zodanig dat een zekere iteratie fpf^p een contractie is. Toon aan dat ff precies één vast punt heeft. Toepassing: de integraaloperator TT op C([0,a])\mathcal C(\intcc0a), Tf(x)=0xfTf(x) = \int_0^x f, voldoet aan Tpap/p!\norm{T^p} \leq a^p/p! — leid af dat u=g+λTuu = g + \lambda Tu voor elke λ\lambda oplosbaar is. (b) (Edelstein) Zij KK compact en f ⁣:KKf\colon K \to K met d(f(x),f(y))<d(x,y)d(f(x), f(y)) < d(x, y) voor xyx \neq y. Toon aan dat ff precies één vast punt heeft, maar dat het contractietempo verloren kan gaan: op X=[1,+)X = [1, +\infty) (volledig, niet compact) heeft f(x)=x+1xf(x) = x + \frac1x geen vast punt, hoewel de afstanden strikt afnemen.

Oplossing

Oplossing van Oefening 7.4.

(a) Zij xx^* het unieke vaste punt van fpf^p. Dan is fp(f(x))=f(fp(x))=f(x)f^p(f(x^*)) = f(f^p(x^*)) = f(x^*): f(x)f(x^*) is een vast punt van fpf^p, dus f(x)=xf(x^*) = x^*. Een vast punt van ff is er een van fpf^p: de eenduidigheid draagt over. Voor TT: per inductie is Tpf(x)fxp/p!\abs{T^pf(x)} \leq \norm f_\infty x^p/p! (elke integratie voegt een factor xk\frac xk toe), dus Tpap/p!\norm{T^p} \leq a^p/p!. De afbeelding S(u)=g+λTuS(u) = g + \lambda Tu voldoet aan Sp(u)Sp(v)=λpTp(uv)S^p(u) - S^p(v) = \lambda^pT^p(u - v), met norm λpap/p!uv0\leq \abs\lambda^pa^p/p!\,\norm{u - v} \to 0: een zekere SpS^p is een contractie, en SS heeft precies één vast punt: de vergelijking van Volterra u=g+λTuu = g + \lambda Tu is voor elke λ\lambda eenduidig oplosbaar.

(b) φ(x)=d(x,f(x))\varphi(x) = d(x, f(x)) is continu op de compacte KK: ze bereikt haar minimum in een zekere x0x_0. Was f(x0)x0f(x_0) \neq x_0, dan is φ(f(x0))=d(f(x0),f2(x0))<d(x0,f(x0))=minφ\varphi(f(x_0)) = d(f(x_0), f^2(x_0)) < d(x_0, f(x_0)) = \min\varphi: absurd. Eenduidigheid: twee vaste punten xyx \ne y zouden d(x,y)=d(f(x),f(y))<d(x,y)d(x,y) = d(f(x), f(y)) < d(x,y) geven. Zonder compactheid: f(x)=x+1xf(x) = x + \frac1x op [1,)[1, \infty) voldoet voor x<yx < y aan f(y)f(x)=(yx)(11xy)<yxf(y) - f(x) = (y - x)\bigl(1 - \frac1{xy}\bigr) < y - x, en toch is altijd f(x)>xf(x) > x: geen vast punt — een strikte afname van de afstanden is zwakker dan een uniforme contractiefactor.

Oefening 7.5 ★★

(a) Twee continue afbeeldingen naar een hausdorffruimte die op een dichte deelverzameling samenvallen, vallen overal samen; waar is dat in dit hoofdstuk gebruikt? (b) Zij DXD \subseteq X dicht en f ⁣:DYf \colon D \to Y een isometrische bijectie op een dichte deelverzameling van een volledige YY, met XX volledig. Toon aan dat ff zich voortzet tot een isometrische bijectie XYX \to Y. Leid daaruit opnieuw de eenduidigheid van de vervollediging af.

Oplossing

Oplossing van Oefening 7.5.

(a) De verzameling A={g=h}A = \{g = h\} is het origineel van de diagonaal ΔY\Delta_Y onder de continue x(g(x),h(x))x \mapsto (g(x), h(x)); en ΔY\Delta_Y is gesloten omdat YY Hausdorff is (voor y1y2y_1 \neq y_2 geven disjuncte open omgevingen een open blok rond (y1,y2)(y_1, y_2) dat de diagonaal mijdt, zodat het complement van ΔY\Delta_Y open is): dus is AA gesloten, bevat het een dichte verzameling en is het gelijk aan XX. Gebruikt: de eenduidigheid in Stelling 7.2, en daarmee in de eenduidigheid van de vervollediging.

(b) ff is een isometrie en dus uniform continu: ze zet zich voort tot F ⁣:XYF\colon X \to Y (Stelling 7.2), nog steeds isometrisch (de betrekking d(F(x),F(x))=d(x,x)d(F(x), F(x')) = d(x, x') geldt op een dichte verzameling paren en beide leden zijn continu). Net zo zet f1 ⁣:f(D)Xf^{-1}\colon f(D) \to X zich voort tot G ⁣:YXG \colon Y \to X. De samenstelling GFG \circ F is continu en houdt de dichte DD vast: ze is idX\mathrm{id}_X (onderdeel (a)); en symmetrisch is FG=idYF \circ G = \mathrm{id}_Y. Dus is FF een isometrische bijectie. Eenduidigheid van de vervollediging: pas dit toe op D=XD = X, dicht gelegen in twee vervolledigingen.

Oefening 7.6 ★★

Welke van de volgende families zijn equicontinu, puntsgewijs begrensd, relatief compact in C([0,1])\mathcal C(\intcc01)?

{xsin(nx)}n,{xxn}n,\{x \mapsto \sin(nx)\}_{n},\qquad \{x \mapsto x^n\}_n,
{fC1:f1, f1},{f:Lip(f)1, f(0)=0}.\Bigl\{f \in \mathcal C^1 : \norm f_\infty \leq 1,\ \norm{f'}_\infty \leq 1\Bigr\},\qquad \{f : \operatorname{Lip}(f)\leq 1,\ f(0) = 0\}.

Verantwoord elk antwoord met Ascoli of met een tegenvoorbeeldrij.

Oplossing

Oplossing van Oefening 7.6.

{sin(nx)}\{\sin(nx)\}: puntsgewijs begrensd door 11; niet equicontinu: in x=0x = 0 is sin(nπ2n)=1\sin(n\cdot\frac{\pi}{2n}) = 1 met π2n0\frac\pi{2n} \to 0, wat elke gemeenschappelijke δ\delta voor ε=12\varepsilon = \frac12 schendt. Niet relatief compact (de noodzakelijkheid in Ascoli, Stelling 7.11).

{xn}\{x^n\}: begrensd; niet equicontinu in 11: 1(1δ)n11 - (1 - \delta)^n \to 1 als nn \to \infty voor vaste δ\delta. Niet relatief compact — en consistent daarmee is haar puntsgewijze limiet discontinu, zodat geen enkele deelrij uniform convergeert.

{f1,f1}\{\norm f_\infty \leq 1, \norm{f'}_\infty \leq 1\}: de middelwaardeongelijkheid maakt de familie 11-lipschitz en dus equicontinu; begrensd: relatief compact volgens Ascoli. (Niet compact: ze is niet gesloten — uniforme limieten hoeven niet C1\mathcal C^1 te zijn; haar afsluiting bestaat uit de 11-lipschitz functies met norm 1\leq 1.)

{Lip(f)1,f(0)=0}\{\operatorname{Lip}(f) \leq 1, f(0) = 0\}: equicontinu; puntsgewijs begrensd (f(x)x1\abs{f(x)} \leq x \leq 1); en gesloten onder uniforme limieten (de lipschitzongelijkheid en de waarde in 00 gaan over op de limiet): compact.

Oefening 7.7 ★★★

(Compacte integraaloperatoren) Zij kC([0,1]2)k \in \mathcal C(\intcc01^2) en, voor fC([0,1])f \in \mathcal C(\intcc01), Tf(x)=01k(x,y)f(y) ⁣dyTf(x) = \int_0^1k(x,y)f(y)\,\dd y. (a) Toon aan dat TT de eenheidsbal van C([0,1])\mathcal C(\intcc01) naar een equicontinue, uniform begrensde verzameling stuurt; besluit dat TT een compacte operator is: beelden van begrensde verzamelingen zijn relatief compact. (b) Leid af dat voor begrensde (fn)(f_n) de rij (Tfn)(Tf_n) een uniform convergente deelrij heeft, en dat TT geen bijectie met continue inverse kan zijn. (Het beeld van de eenheidsbal zou een compacte omgeving van 00 in C([0,1])\mathcal C(\intcc01) zijn: verboden door de stelling van Riesz uit bachelorjaar 2.)

Oplossing

Oplossing van Oefening 7.7.

(a) Voor f1\norm f_\infty \leq 1 is Tf(x)k\abs{Tf(x)} \leq \norm k_\infty, en

Tf(x)Tf(x)01k(x,y)k(x,y) ⁣dyωk(xx),\abs{Tf(x) - Tf(x')} \leq \int_0^1\abs{k(x,y) - k(x',y)}\,\dd y \leq \omega_k\bigl(\abs{x - x'}\bigr),

waarbij ωk\omega_k een modulus van uniforme continuïteit van kk op het compacte vierkant is (Heine): het beeld van de eenheidsbal is uniform begrensd en equicontinu, en dus relatief compact (Ascoli). Wegens de lineariteit heeft elke begrensde verzameling een relatief compact beeld: TT is een compacte operator.

(b) De uitspraak over deelrijen is de definitie van relatieve compactheid, toegepast op {Tfn}\{Tf_n\}. Was TT bijectief met continue inverse, dan is T(B(0,1))B(0,ε)T(B(0,1)) \supseteq B(0, \varepsilon) voor zekere ε>0\varepsilon > 0 (T1T^{-1} is continu in 00); en dan was de gesloten bal Bˉ(0,ε)\bar B(0,\varepsilon), een gesloten deelverzameling van de compacte T(B(0,1))\overline{T(B(0,1))}, compact — onmogelijk in het oneindigdimensionale C([0,1])\mathcal C(\intcc01) volgens de stelling van Riesz (bachelorjaar 2).

Oefening 7.8 ★★

Bewijs of weerleg, voor continue fn ⁣:[0,1]Rf_n \colon \intcc01 \to \R: (a) uit fn0f_n \downarrow 0 puntsgewijs volgt uniform (Dini — bewijs het opnieuw); (b) hetzelfde zonder monotonie; (c) hetzelfde met monotonie maar met discontinue ff; (d) hetzelfde met monotonie en continue limiet, maar op (0,1)\intoo01.

Oplossing

Oplossing van Oefening 7.8.

(a) Dini: zie Lemma 7.13 — het overdekkingsargument. (b) Onjuist: de bewegende bult fn(x)=max(0,1nx1)f_n(x) = \max(0, 1 - \abs{nx - 1}) gaat puntsgewijs naar 00 (voor x>0x > 0 is fn(x)=0f_n(x) = 0 zodra n>2/xn > 2/x; en fn(0)=0f_n(0) = 0), maar fn=1\norm{f_n}_\infty = 1. (c) Onjuist: fn(x)=xnf_n(x) = x^n daalt naar de discontinue 1{1}\mathbf 1_{\{1\}}, en sup[0,1]fnf=sup[0,1)xn=1\sup_{[0,1]}\abs{f_n - f} = \sup_{[0,1)} x^n = 1. (d) Onjuist: xn0x^n \downarrow 0 puntsgewijs op het niet-compacte (0,1)\intoo01, met sup(0,1)xn=1\sup_{(0,1)}x^n = 1. Elke hypothese van Dini is dus nodig.

Oefening 7.9 ★★

(a) (Momenten leggen vast) Zij fC([0,1])f \in \mathcal C(\intcc01) met 01f(x)xn ⁣dx=0\int_0^1 f(x)\,x^n\,\dd x = 0 voor alle nNn \in \N. Toon aan dat f=0f = 0. (Benader ff uniform door veeltermen en bereken f2\int f^2.) (b) Toon aan dat de even veeltermen dicht liggen in C([0,1])\mathcal C(\intcc01) maar niet in C([1,1])\mathcal C(\intcc{-1}1); welke hypothese van Stone–Weierstrass faalt daar? (c) Ligt de algebra voortgebracht door xeixx \mapsto \eu^{\iu x} alleen (zonder eix\eu^{-\iu x}) dicht in C(S1,C)\mathcal C(S^1, \C)? (Beschouw 02πf(t)eit ⁣dt\int_0^{2\pi}f(t)\,\eu^{\iu t}\dd t.)

Oplossing

Oplossing van Oefening 7.9.

(a) Wegens de lineariteit is 01fP=0\int_0^1 fP = 0 voor elke veelterm PP. Kies PnfP_n \to f uniform (Gevolg 7.16): dan is 01f2=lim01fPn=0\int_0^1 f^2 = \lim\int_0^1 fP_n = 0, en de continue f20f^2 \geq 0 met integraal nul verdwijnt identiek.

(b) Op [0,1]\intcc01 vormen de veeltermen in x2x^2 een algebra met de constanten die de punten scheidt (xx2x \mapsto x^2 is injectief op [0,1]\intcc01): dicht volgens Stone–Weierstrass. Op [1,1]\intcc{-1}1 neemt x2x^2 in ±x\pm x dezelfde waarde aan, en elke veelterm in x2x^2 ook: een uniforme limiet van zulke functies is even. Convergeerden even functies gng_n uniform naar de identiteit, dan was x=limgn(x)=limgn(x)=xx = \lim g_n(x) = \lim g_n(-x) = -x voor alle xx: absurd — niet dicht. De scheidingshypothese faalt bij de paren {x,x}\{x, -x\}.

(c) Nee. Voor ff in de algebra A\mathcal A voortgebracht door de constanten en eit\eu^{\iu t} — lineaire combinaties van eint\eu^{\iu nt} met n0n \geq 0 — is Λ(f)=02πf(t)eit ⁣dt=0\Lambda(f) = \int_0^{2\pi}f(t)\,\eu^{\iu t}\,\dd t = 0 (elke 02πei(n+1)t ⁣dt=0\int_0^{2\pi}\eu^{\iu(n+1)t}\dd t = 0). Λ\Lambda is continu voor \norm\cdot_\infty (Λ(f)2πf\abs{\Lambda(f)}\leq 2\pi\norm f_\infty), dus verdwijnt Λ\Lambda op Aˉ\bar{\mathcal A}; maar Λ(eit)=2π0\Lambda(\eu^{-\iu t}) = 2\pi \neq 0: dus eitAˉ\eu^{-\iu t} \notin \bar{\mathcal A}. (Stone–Weierstrass is niet van toepassing: A\mathcal A is niet gesloten onder toevoeging — en die belemmering is precies degene die de theorie van de holomorfe functies in Hoofdstuk 16 zal systematiseren.)

Oefening 7.10 ★★★

(Uniforme begrensdheid, metrische versie) Zij XX een volledige metrische ruimte en FC(X,R)\mathcal F \subseteq \mathcal C(X, \R) een familie die puntsgewijs begrensd is: supfFf(x)<\sup_{f \in \mathcal F}\abs{f(x)} < \infty voor elke xx. Toon aan dat er een niet-lege open UXU \subseteq X is waarop F\mathcal F uniform begrensd is: supfsupUf<\sup_{f}\sup_U \abs f < \infty. (Beschouw Fn={x:f(x)n f}F_n = \{x : \abs{f(x)} \leq n\ \forall f\}.) Dit is de motor achter Banach–Steinhaus in Hoofdstuk 8.

Oplossing

Oplossing van Oefening 7.10.

Fn=fF{x:f(x)n}F_n = \bigcap_{f \in \mathcal F}\{x : \abs{f(x)} \leq n\} is een doorsnede van gesloten verzamelingen: gesloten. De puntsgewijze begrensdheid geeft X=nFnX = \bigcup_n F_n. Baire (Stelling 7.6) levert een n0n_0 met U=F˚n0U = \mathring F_{n_0} \neq \varnothing: op UU is fn0\abs f \leq n_0 voor elke fFf \in \mathcal F tegelijk.

Oefening 7.11 ★★

(C1\mathcal C^1 heeft haar eigen norm nodig) Beschouw op E=C1([0,1],R)E = \mathcal C^1(\intcc01, \R) de norm fC1=f+f\norm f_{\mathcal C^1} = \norm f_\infty + \norm{f'}_\infty. (a) Toon aan dat (E,C1)(E, \norm\cdot_{\mathcal C^1}) volledig is (voor een C1\mathcal C^1-Cauchyrij geldt fnff_n \to f en fngf_n' \to g uniform; identificeer g=fg = f' door in fn(x)=fn(0)+0xfnf_n(x) = f_n(0) + \int_0^xf_n' tot de limiet over te gaan). (b) Toon aan dat (E,)(E, \norm\cdot_\infty) niet volledig is: geef een uniforme limiet van C1\mathcal C^1-functies die niet differentieerbaar is (bijvoorbeeld gladde benaderingen van x12\abs{x - \tfrac12}). (c) Leid uit (a), (b) en de ideeënkring rond de open-afbeeldingsstelling — of rechtstreeks — af dat geen enkele constante CC voldoet aan fCf\norm{f'}_\infty \leq C\,\norm f_\infty op EE: geef een rij die dat aantoont. Differentiëren is onbegrensd; dat is de afgrond achter Stelling 7.9.

Oplossing

Oplossing van Oefening 7.11.

(a) Een Cauchyrij voor C1\norm\cdot_{\mathcal C^1} is uniform Cauchy, samen met haar afgeleiden: fnff_n \to f en fngf_n' \to g uniform, met f,gf, g continu. Overgaan tot de limiet (de uniforme convergentie laat dat onder de integraal toe) in fn(x)=fn(0)+0xfn(t) ⁣dtf_n(x) = f_n(0) + \int_0^xf_n'(t)\,\dd t geeft f(x)=f(0)+0xgf(x) = f(0) + \int_0^xg: ff is C1\mathcal C^1 met f=gf' = g, en fnfC10\norm{f_n - f}_{\mathcal C^1} \to 0. Volledig.

(b) h(x)=x12h(x) = \abs{x - \frac12} is een uniforme limiet van C1\mathcal C^1-functies, bijvoorbeeld hn(x)=(x12)2+1nh_n(x) = \sqrt{(x - \frac12)^2 + \frac1n} (met hnh1n\abs{h_n - h} \leq \frac1{\sqrt n} via de toegevoegde grootheid), en toch is hEh \notin E: de supnorm op EE is niet volledig — haar vervollediging is C([0,1])\mathcal C(\intcc01).

(c) fn(x)=sin(2πnx)f_n(x) = \sin(2\pi nx) heeft fn=1\norm{f_n}_\infty = 1 en fn=2πn\norm{f_n'}_\infty = 2\pi n \to \infty: er is geen CC. (Conceptueel: was het differentiëren begrensd voor de supnorm, dan waren de twee normen uit (a)–(b) gelijkwaardig, wat (E,)(E, \norm\cdot_\infty) volledig zou maken — in strijd met (b). Die onbegrensdheid is precies waarom generieke continue functies nergens differentieerbaar kunnen zijn, Stelling 7.9.)

Oefening 7.12 ★★★

(Lemma van Croft) Zij f ⁣:(0,)Rf\colon\intoo0\infty\to\R continu en stel dat voor elke x>0x > 0 geldt f(nx)0f(nx) \to 0 als het gehele getal nn \to \infty. Toon aan dat f(t)0f(t) \to 0 als t+t \to +\infty. (Leg ε>0\varepsilon > 0 vast; de verzamelingen FN={x>0:f(nx)ε  nN}F_N = \{x > 0 : \abs{f(nx)} \leq \varepsilon\ \ \forall n \geq N\} zijn gesloten en overdekken (0,)\intoo0\infty; Baire in een zeker interval [a,b]\intcc ab geeft een NN en een deelinterval [a,b]FN\intcc{a'}{b'} \subseteq F_N; en dan overdekken de uitrekkingen [na,nb]\bigl[na', nb'\bigr] met nNn \geq N een hele omgeving van ++\infty zodra n(ba)an(b' - a') \geq a'.) Waar wordt de hypothese “voor elke xx” (en niet alleen voor rationale xx) gebruikt?

Oplossing

Oplossing van Oefening 7.12.

Leg ε>0\varepsilon > 0 vast. Elke FNF_N is een doorsnede over nNn \geq N van originelen van de gesloten [ε,ε]\intcc{-\varepsilon} \varepsilon onder de continue xf(nx)x \mapsto f(nx): gesloten. De hypothese zegt dat elke x>0x > 0 in een zekere FNF_N ligt. Baire, toegepast binnen het volledige [a,b]\intcc ab (met willekeurige 0<a<b0 < a < b), maakt een zekere FNF_N dicht in een deelinterval; en omdat ze gesloten is, bevat ze een interval [a,b]\intcc{a'}{b'} met 0<a<b0 < a' < b'. Dan overlappen voor elke nmax(N,aba)n \geq \max(N, \frac{a'}{b' - a'}) de intervallen [na,nb]\intcc{na'}{nb'} en [(n+1)a,(n+1)b]\intcc{(n+1)a'}{(n+1)b'} elkaar (want nb(n+1)anb' \geq (n+1)a'), dus

nn0[na,nb][n0a,+),\bigcup_{n \geq n_0}\intcc{na'}{nb'} \supseteq \intco{n_0a'}{+\infty},

en elke tn0at \geq n_0a' is t=nxt = nx met nNn \geq N en x[a,b]FNx \in \intcc{a'}{b'} \subseteq F_N: dus f(t)ε\abs{f(t)} \leq \varepsilon. Bijgevolg is lim suptfε\limsup_{t\to\infty}\abs f \leq \varepsilon voor elke ε\varepsilon: f0f \to 0. De volledige hypothese is nodig omdat FNF_N een heel interval aan xx-waarden moet overdekken — met alleen rationale xx is de vereniging van de FNF_N aftelbaar en geeft Baire niets; en inderdaad bestaan er continue tegenvoorbeelden die langs alle rationale stralen verdwijnen maar niet in het oneindige.

7.6 Probleem: de bestaansstelling van Peano

Probleem 7.1

Weekendopgave — bestaan van oplossingen van x=f(t,x)x' = f(t, x) zonder Lipschitz

Cauchy–Lipschitz (bachelorjaar 2; opnieuw bewezen in Hoofdstuk 19) eist dat ff lipschitz is in xx. Peano (1890): de continuïteit van ff levert al het bestaan — zij het niet de eenduidigheid. Wij bewijzen het met eulerpolygonen en Ascoli. Kader: f ⁣:RRdf \colon R \to \R^d is continu op de rechthoek R=[t0a,t0+a]×Bˉ(x0,b)R×RdR = \intcc{t_0 - a}{t_0 + a}\times \bar B(x_0, b) \subseteq \R\times\R^d, M=supRfM = \sup_R\norm f, en

T=min(a, b/M)(is M=0, dan is het probleem triviaal).T = \min\bigl(a,\ b/M\bigr) \qquad (\text{is } M = 0 \text{, dan is het probleem triviaal}).

Deel I — Eulerpolygonen. Verdeel voor n1n \geq 1 het interval [t0,t0+T][t_0, t_0 + T] door tk=t0+kT/nt_k = t_0 + kT/n (0kn0 \leq k \leq n) en definieer φn\varphi_n stuksgewijs affien: φn(t0)=x0\varphi_n(t_0) = x_0 en, op [tk,tk+1][t_k, t_{k+1}],

φn(t)=φn(tk)+(ttk)f(tk,φn(tk)).\varphi_n(t) = \varphi_n(t_k) + (t - t_k)\,f\bigl(t_k, \varphi_n(t_k)\bigr).
  1. Toon met inductie aan dat φn\varphi_n goed gedefinieerd is, met φn(t)x0M(tt0)b\norm{\varphi_n(t) - x_0} \leq M(t - t_0) \leq b op [t0,t0+T][t_0, t_0 + T] — zodat de evaluatiepunten in RR blijven. (Hier komt Tb/MT \leq b/M binnen.)
  2. Toon aan dat elke φn\varphi_n MM-lipschitz is.
  3. Leid uit Arzelà–Ascoli (Stelling 7.11) af dat een deelrij (φnj)(\varphi_{n_j}) uniform op [t0,t0+T][t_0, t_0 + T] naar een zekere φ\varphi convergeert, zelf MM-lipschitz met φ(t0)=x0\varphi(t_0) = x_0.

Deel II — De limiet lost de vergelijking op. Definieer het defect Δn(t)=φn(t)f(t,φn(t))\Delta_n(t) = \varphi_n'(t) - f\bigl(t, \varphi_n(t)\bigr) in de punten buiten het rooster.

  1. Toon aan dat ff uniform continu is op RR, en leid af: voor elke ε>0\varepsilon > 0 is er een n0n_0 zodanig dat voor nn0n \geq n_0 en elke tt buiten het rooster Δn(t)ε\norm{\Delta_n(t)} \leq \varepsilon. (Op (tk,tk+1)(t_k, t_{k+1}) is φn(t)=f(tk,φn(tk))\varphi_n'(t) = f(t_k, \varphi_n(t_k)), en ligt (t,φn(t))(t, \varphi_n(t)) op afstand hoogstens (1+M)T/n(1 + M)\,T/n van (tk,φn(tk))(t_k, \varphi_n(t_k)).)
  2. Stel de integraalvorm op: voor alle tt is

    φn(t)=x0+t0tf(s,φn(s)) ⁣ds+t0tΔn(s) ⁣ds,\varphi_n(t) = x_0 + \int_{t_0}^{t} f\bigl(s, \varphi_n(s)\bigr)\dd s + \int_{t_0}^t \Delta_n(s)\,\dd s,

    waarbij de middelste integrand stuksgewijs continu is.

  3. Ga langs (nj)(n_j) tot de limiet over: toon aan dat f(s,φnj(s))f(s,φ(s))f(s, \varphi_{n_j}(s)) \to f(s, \varphi(s)) uniform (opnieuw de uniforme continuïteit van ff), en besluit

    φ(t)=x0+t0tf(s,φ(s)) ⁣ds.\varphi(t) = x_0 + \int_{t_0}^t f\bigl(s, \varphi(s)\bigr)\dd s .
  4. Leid af dat φ\varphi C1\mathcal C^1 is op [t0,t0+T][t_0, t_0 + T] en x=f(t,x)x' = f(t, x) met x(t0)=x0x(t_0) = x_0 oplost; en breid de constructie uit tot [t0T,t0][t_0 - T, t_0] (tijdomkering). Dat is de stelling van Peano.

Deel III — De eenduidigheid faalt werkelijk. Beschouw x=2xx' = 2\sqrt{\abs x} met x(0)=0x(0) = 0, op R\R.

  1. Ga na dat f(x)=2xf(x) = 2\sqrt{\abs x} continu is maar op geen enkele omgeving van 00 lipschitz.
  2. Ga na dat x0x \equiv 0 en, voor elke c0c \geq 0,

    xc(t)={0tc,(tc)2tc,x_c(t) = \begin{cases} 0 & t \leq c,\\ (t - c)^2 & t \geq c,\end{cases}

    allemaal oplossingen door (0,0)(0,0) zijn: een continuüm van verschillende oplossingen.

  3. Waar loopt het argument met de picarditeratie (het vaste punt van Banach) vast voor deze ff?

Deel IV — Grenzen van de methode.

  1. Toon aan dat de stelling van Peano faalt in oneindige dimensie: we nemen het klassieke voorbeeld van Dieudonné in de ruimte c0c_0 van nulrijen zonder bewijs aan; bewijs in plaats daarvan het eindigdimensionale ingrediënt dat daar faalt: de gesloten eenheidsbal van c0c_0 (supnorm) is niet compact — geef een begrensde rij zonder convergente deelrij, en leg uit welke stap van Deel I breekt.
  2. Vat samen: welke hypothesen geven bestaan? En welke bestaan én eenduidigheid? Formuleer de twee stellingen (Peano; Cauchy–Lipschitz) nauwkeurig naast elkaar.

Deel V — Osgood: eenduidigheid voorbij Lipschitz. Zij ω ⁣:(0,)(0,)\omega\colon\intoo0\infty\to\intoo0\infty continu en niet-dalend, met

01 ⁣drω(r)=+(divergentie in 0),\int_0^1\frac{\dd r}{\omega(r)} = +\infty \qquad(\text{divergentie in } 0),

en stel dat ff op RR voldoet aan f(t,x)f(t,y)ω(xy)\norm{f(t, x) - f(t, y)} \leq \omega\bigl(\norm{x - y}\bigr).

  1. Ga na dat ω(r)=Lr\omega(r) = Lr voldoet (Lipschitz), dat ω(r)=rlog1r\omega(r) = r\log\frac1r (voor kleine rr, continu voortgezet) voldoet hoewel het geen O(r)O(r) is, en dat ω(r)=2r\omega(r) = 2\sqrt r niet voldoet. Bereken de integraal in elk geval.
  2. Zij x1,x2x_1, x_2 oplossingen op [t0,t0+T][t_0, t_0 + T] met x1(t0)=x2(t0)x_1(t_0) = x_2(t_0), en δ(t)=x1(t)x2(t)\delta(t) = \norm{x_1(t) - x_2(t)}. Toon uitsluitend met de integraalvormen aan dat voor t0stt_0 \leq s \leq t

    δ(t)δ(s)+stω(δ(v)) ⁣dv.\delta(t) \leq \delta(s) + \int_s^t\omega\bigl(\delta(v)\bigr)\dd v .
  3. (Stelling van Osgood) Stel δ(t1)>0\delta(t_1) > 0 voor zekere t1t_1, en zij τ=sup{tt1:δ(t)=0}\tau = \sup\{t \leq t_1 : \delta(t) = 0\}. Stel voor s(τ,t1)s \in \intoo\tau{t_1}: u(t)=δ(s)+stω(δ(v)) ⁣dvu(t) = \delta(s) + \int_s^t\omega(\delta(v))\,\dd v. Toon aan dat δu\delta \leq u en u=ω(δ)ω(u)u' = \omega(\delta) \leq \omega(u), en leid af dat

    δ(s)u(t1) ⁣drω(r)t1st1τ.\int_{\delta(s)}^{u(t_1)}\frac{\dd r}{\omega(r)} \leq t_1 - s \leq t_1 - \tau .

    Laat sτs \downarrow \tau gaan en leid een tegenspraak af met de divergentie van de integraal. Besluit: oplossingen door een gemeenschappelijke beginvoorwaarde vallen samen — eenduidigheid onder de voorwaarde van Osgood.

  4. Trek de gevolgen: de eenduidigheid van Cauchy–Lipschitz is het geval ω(r)=Lr\omega(r) = Lr; de vergelijking x=xlog1xx' = x\log\frac1{\abs x} (voortgezet met 00 in 00) heeft eenduidige oplossingen hoewel haar rechterlid in 00 niet lipschitz is; en voor x=2xx' = 2\sqrt{\abs x} is de convergentie van 0 ⁣dr2r\int_0\frac{\dd r}{2\sqrt r} precies wat een oplossing toestaat om 00 in eindige tijd te verlaten — breng de waarde van de integraal 0h2 ⁣dr2r=h\int_0^{h^2}\frac{\dd r}{2\sqrt r} = h in verband met het ontsnappingsgedrag van xcx_c.

Deel VI — Snelheden, schema’s, trechters.

  1. (Integraallemma van Grönwall) Zij e,η0e, \eta \geq 0 continu op [t0,t0+T][t_0, t_0+T] met e(t)t0t(Le(s)+η(s)) ⁣dse(t) \leq \int_{t_0} ^t\bigl(L\,e(s) + \eta(s)\bigr)\dd s voor alle tt. Toon aan dat

    e(t)supηL(eL(tt0)1)e(t) \leq \frac{\sup\eta}{L}\bigl(\eu^{L(t - t_0)} - 1\bigr)

    (stel G(t)=t0t(Le+η)G(t) = \int_{t_0}^t(Le + \eta), merk op dat GLG+supηG' \leq LG + \sup\eta, en differentieer eLtG(t)\eu^{-Lt}G(t)).

  2. (Euler convergeert met een snelheid) Stel nu dat ff LL-lipschitz is in xx en LL'-lipschitz in tt op RR. Bewijs, door de defectgrens van vraag 4 (kwantitatief gemaakt: Δn(L+LM)T/n\norm{\Delta_n} \leq (L' + LM)T/n) met vraag 17 te combineren, dat

    sup[t0,t0+T]φnφ    (L+LM)TneLT1L=O(1n),\sup_{[t_0, t_0+T]}\norm{\varphi_n - \varphi} \;\leq\; \frac{(L' + LM)\,T}{n}\cdot \frac{\eu^{LT} - 1}{L} = O\Bigl(\frac1n\Bigr),

    waarbij φ\varphi de oplossing is: met lipschitzgegevens convergeert de hele rij, met een expliciete snelheid — geen deelrijen nodig. Waarom tilt de eenduidigheid convergentie langs deelrijen ook zonder deze berekening op tot volledige convergentie?

  3. (Het schema kiest) Toon voor x=2xx' = 2\sqrt{\abs x} met x(0)=0x(0) = 0 aan dat elke eulerpolygoon identiek nul is, zodat het schema naar de oplossing x0x \equiv 0 convergeert; maar gestart in x(0)=ε>0x(0) = \varepsilon > 0 convergeert het (eerst nn \to \infty, dan ε0\varepsilon \to 0) naar tt2t \mapsto t^2, een andere oplossing door de oorsprong. De niet-eenduidigheid duikt weer op als gevoeligheid van het numerieke schema voor verstoringen.
  4. Toon aan dat de verzameling S\mathcal S van alle oplossingen van x=f(t,x)x' = f(t,x) met x(t0)=x0x(t_0) = x_0 op [t0,t0+T][t_0, t_0+T] (met waarden in Bˉ(x0,b)\bar B(x_0, b)) niet leeg is (Deel II), uniform MM-lipschitz, en gesloten in (C([t0,t0+T],Rd),)\bigl(\mathcal C([t_0, t_0+T], \R^d), \norm\cdot_\infty\bigr); besluit met Ascoli dat S\mathcal S compact is. (De stelling van Kneser voegt toe dat S\mathcal S samenhangend is; die bewijzen we niet.)
  5. Ga het verschijnsel van Kneser na op het voorbeeld: toon voor x=2xx' = 2\sqrt{\abs x} met x(0)=0x(0) = 0 op [0,1][0, 1] aan dat S={xc:c[0,]}\mathcal S = \{x_c : c \in \intcc0\infty\} met x0x_\infty \equiv 0 (neem voor een willekeurige oplossing c=sup{t:x(t)=0}c = \sup\{t : x(t) = 0\} en integreer (x)=1(\sqrt x)' = 1 op {x>0}\{x > 0\}), dat cxcc \mapsto x_c continu is van [0,]\intcc0\infty (de eenpuntscompactificatie, met c=c = \infty als limiet aangelijmd) naar C([0,1])\mathcal C([0,1]), en besluit dat S\mathcal S inderdaad compact en samenhangend is — een trechter in de vorm van een lijnstuk.
  6. (Bereikbare verzamelingen) Leid uit vraag 20 af dat voor elke vaste tt de bereikbare verzameling S(t)={x(t):xS}\mathcal S(t) = \{x(t) : x \in \mathcal S\} compact is; bereken haar voor het voorbeeld van vraag 21 en ga na dat ze ook samenhangend is: S(t)=[0,t2]\mathcal S(t) = \intcc0{t^2} — elke tussenliggende toestand wordt door een zekere oplossing bereikt.

Deel VII — Aanvullingen: afhankelijkheid, optimaliteit, en een schema met de hand berekend.

  1. (Continue afhankelijkheid) Stel dat ff op RR LL-lipschitz is in xx, en zij x,yx, y twee oplossingen met beginwaarden x0x_0 en y0y_0 in t0t_0. Pas het bewijs van vraag 17 aan de ongelijkheid e(t)x0y0+t0tLe(s) ⁣dse(t) \leq \norm{x_0 - y_0} + \int_{t_0}^tL\,e(s)\dd s aan en toon aan dat

    x(t)y(t)x0y0eL(tt0),\norm{x(t) - y(t)} \leq \norm{x_0 - y_0}\, \eu^{L(t - t_0)},

    en ga op x=Lxx' = Lx na dat de grens bereikt wordt: Grönwall is scherp. Leid opnieuw de eenduidigheid af (x0=y0x_0 = y_0), en dat de stromingsafbeelding x0x(t;x0)x_0 \mapsto x(t; x_0) lipschitz is met constante eLT\eu^{LT}, overal waar ze gedefinieerd is.

  2. (Osgood is optimaal) Zij omgekeerd ω\omega continu, niet-dalend en positief op (0,)\intoo0\infty, met ω(r)0\omega(r) \to 0 als r0+r \to 0^+ maar

    01 ⁣drω(r)<+,\int_0^1\frac{\dd r}{\omega(r)} < +\infty,

    en zet f(x)=ω(x)f(x) = \omega(\abs x) voort met f(0)=0f(0) = 0. Toon aan dat Ω(x)=0x ⁣drω(r)\Omega(x) = \int_0^x\frac{\dd r}{\omega(r)} een stijgende bijectie van (0,1]\intoc0{1} op (0,Ω(1)]\intoc0{\Omega(1)} is, dat haar inverse gg de vergelijking g=ω(g)g' = \omega(g) met g(0+)=0g(0^+) = 0 oplost, en dat gg, met 00 voortgezet voor t0t \leq 0, een C1\mathcal C^1-oplossing van x=f(x)x' = f(x) door (0,0)(0, 0) is die van x0x \equiv 0 verschilt (voor g(0)=0g'(0) = 0: schat g(t)t\frac{g(t)}t af met ω(g(t))\omega(g(t))). Besluit: de divergentiehypothese van vraag 15 is geen gemak maar de exacte grens van de eenduidigheid; vind Deel III terug uit ω(r)=2r\omega(r) = 2\sqrt r en Ω(x)=x\Omega(x) = \sqrt x.

  3. (Euler met de hand berekend) Toon voor x=xx' = x met x(0)=1x(0) = 1 op [0,1][0, 1] aan dat de eulerpolygoon voldoet aan φn(1)=(1+1n)n\varphi_n(1) = \bigl(1 + \frac1n\bigr)^n. Bewijs de ontwikkeling

    (1+1n)n=e(112n+O(1n2)),\Bigl(1 + \frac1n\Bigr)^{n} = \eu\Bigl(1 - \frac1{2n} + O\Bigl(\frac1{n^2}\Bigr) \Bigr),

    zodat de fout in t=1t = 1 gelijk is aan e2n+O(n2)\frac{\eu}{2n} + O(n^{-2}): de snelheid O(1n)O(\frac1n) van vraag 18, met de exacte constante. Controleer numeriek voor n=10n = 10: 1.110=2.593741.1^{10} = 2.59374 tegenover e2.71828\eu \approx 2.71828, een fout 0.124540.12454, te vergelijken met e200.13591\frac{\eu}{20} \approx 0.13591.

Oplossing

Oplossing van Probleem 7.1.

1. Inductie naar kk: is φn(tk)x0M(tkt0)MTb\norm{\varphi_n(t_k) - x_0} \leq M(t_k - t_0) \leq MT \leq b, dan ligt het punt (tk,φn(tk))(t_k, \varphi_n(t_k)) in RR, zodat de helling f(tk,φn(tk))f(t_k, \varphi_n(t_k)) gedefinieerd is met norm M\leq M; en dan is voor t[tk,tk+1]t \in [t_k, t_{k+1}]: φn(t)x0φn(tk)x0+M(ttk)M(tt0)MTb\norm{\varphi_n(t) - x_0} \leq \norm{\varphi_n(t_k) - x_0} + M(t - t_k) \leq M(t - t_0) \leq MT \leq b.

2. Elk affien stuk heeft een helling van norm M\leq M; en een stuksgewijs affiene functie met hellingen begrensd door MM is MM-lipschitz (ketting via de roosterpunten).

3. De familie (φn)(\varphi_n) is puntsgewijs begrensd (waarden in Bˉ(x0,b)\bar B(x_0, b)) en equicontinu (gemeenschappelijke lipschitzconstante MM): Ascoli (Stelling 7.11) haalt er φnjφ\varphi_{n_j} \to \varphi uniform op [t0,t0+T][t_0, t_0 + T] uit. De grenzen gaan over op de limiet: φ\varphi is MM-lipschitz met φ(t0)=x0\varphi(t_0) = x_0.

4. RR is compact en ff continu: dus uniform continu (Heine, Gevolg 6.17); zij δ(ε)\delta(\varepsilon) een modulus. Voor t(tk,tk+1)t \in (t_k, t_{k+1}) buiten het rooster is φn(t)=f(tk,φn(tk))\varphi_n'(t) = f(t_k, \varphi_n(t_k)), en verschillen de twee evaluatiepunten van ff met ttkT/n\abs{t - t_k} \leq T/n in de tijd en φn(t)φn(tk)MT/n\norm{\varphi_n(t) - \varphi_n(t_k)} \leq MT/n in de ruimte. Voor nn0n \geq n_0 met (1+M)T/n0<δ(ε)(1 + M)T/n_0 < \delta(\varepsilon) is dus Δn(t)ε\norm{\Delta_n(t)} \leq \varepsilon.

5. Op elke [tk,tk+1][t_k, t_{k+1}] is φn\varphi_n affien, dus φn(tk+1)φn(tk)=tktk+1φn(s) ⁣ds\varphi_n(t_{k+1}) - \varphi_n(t_k) = \int_{t_k}^{t_{k+1}} \varphi_n'(s)\dd s, want de afgeleide is de constante helling; sommeren over de stukken (en het laatste bij tt afsnijden) geeft φn(t)=x0+t0tφn(s) ⁣ds\varphi_n(t) = x_0 + \int_{t_0}^t\varphi_n'(s)\,\dd s. Schrijven we φn=f(s,φn(s))+Δn(s)\varphi_n' = f(s, \varphi_n(s)) + \Delta_n(s) (stuksgewijs continue integranden, met eindig veel sprongen), dan volgt de getoonde formule.

6. Bij gegeven ε\varepsilon is voor grote jj: φnjφ<δ(ε)\norm{\varphi_{n_j} - \varphi}_\infty < \delta(\varepsilon), dus f(s,φnj(s))f(s,φ(s))ε\norm{f(s, \varphi_{n_j}(s)) - f(s, \varphi(s))} \leq \varepsilon voor alle ss: uniforme convergentie van de integranden, en t0tf(s,φnj(s)) ⁣dst0tf(s,φ(s)) ⁣ds\int_{t_0}^t f(s,\varphi_{n_j}(s))\dd s \to \int_{t_0}^tf(s, \varphi(s))\dd s, uniform in tt. Verder is t0tΔnjTsupΔnj0\norm{\int_{t_0}^t\Delta_{n_j}} \leq T\sup\norm{\Delta_{n_j}} \to 0 (vraag 4). Overgaan tot de limiet in de identiteit van vraag 5 geeft φ(t)=x0+t0tf(s,φ(s)) ⁣ds\varphi(t) = x_0 + \int_{t_0}^tf(s, \varphi(s))\,\dd s.

7. De integrand sf(s,φ(s))s \mapsto f(s, \varphi(s)) is continu, dus is het rechterlid C1\mathcal C^1 in tt met afgeleide f(t,φ(t))f(t, \varphi(t)): φ\varphi lost het beginwaardeprobleem op [t0,t0+T][t_0, t_0 + T] op. Voor de linkerhelft: stel g(t,x)=f(2t0t,x)g(t, x) = -f(2t_0 - t, x), continu op de gespiegelde rechthoek met dezelfde grens MM; een oplossing ψ\psi van y=g(t,y)y' = g(t,y) met y(t0)=x0y(t_0) = x_0 op [t0,t0+T][t_0, t_0 + T] levert φ(t)=ψ(2t0t)\varphi(t) = \psi(2t_0 - t), die de oorspronkelijke vergelijking op [t0T,t0][t_0 - T, t_0] oplost; en de twee helften verlijmen tot een C1\mathcal C^1-oplossing (beide eenzijdige afgeleiden in t0t_0 zijn f(t0,x0)f(t_0, x_0)). — Dat is de stelling van Peano: een continue ff laat door elke beginvoorwaarde een lokale oplossing toe.

8. De continuïteit is duidelijk. Lipschitz nabij 00 faalt: f(h)f(0)=2h\abs{f(h) - f(0)} = 2\sqrt h, en 2hLh2\sqrt h \leq L h is onjuist voor h<4/L2h < 4/L^2.

9. Voor tct \leq c lost xc0x_c \equiv 0 op. Voor tct \geq c is xc(t)=2(tc)=2(tc)2=2xcx_c'(t) = 2(t - c) = 2\sqrt{(t-c)^2} = 2\sqrt{\abs{x_c}}. In t=ct = c zijn beide eenzijdige afgeleiden 00: xcx_c is C1\mathcal C^1 en lost overal op, met xc(0)=0x_c(0) = 0 voor elke c0c \geq 0 — samen met x0x \equiv 0 een continuüm van oplossingen door de oorsprong.

10. De picarditeratie stelt Φ(x)(t)=x0+0tf(x(s)) ⁣ds\Phi(x)(t) = x_0 + \int_0^t f(x(s))\dd s op en heeft Φ(x)Φ(y)kxy\norm{\Phi(x) - \Phi(y)} \leq k\norm{x - y} met k<1k < 1 nodig op een geschikte bal — wat volgt uit een lipschitzgrens op ff, onder de integraal overgedragen. Hier laat \sqrt\cdot nabij 00 geen enkele lipschitzgrens toe, en geen keuze van interval of bal repareert dat. En inderdaad kon geen bewijs van eenduidigheid slagen: de eenduidigheid is onjuist (vraag 9).

11. In c0c_0 voldoen de eenheidsvectoren en=(0,,0,1,0,)e_n = (0, \dots, 0, 1, 0, \dots) aan enem=1\norm{e_n - e_m}_\infty = 1 voor nmn \neq m: geen enkele deelrij is Cauchy, dus is de gesloten eenheidsbal niet compact. De stap die breekt is de extractie (vraag 3): Ascoli voor C([t0,t0+T],E)\mathcal C([t_0, t_0+T], E) vergt dat de waarden leven in een ruimte waar begrensde verzamelingen relatief compact zijn — waar in Rd\R^d (Bolzano–Weierstrass), onjuist in c0c_0; puntsgewijze extractie is er niet meer beschikbaar (en inderdaad heeft het voorbeeld van Dieudonné helemaal geen lokale oplossing).

12. Peano: is ff continu op een omgeving van (t0,x0)(t_0, x_0) in R×Rd\R\times\R^d, dan bestaat er een C1\mathcal C^1-oplossing van x=f(t,x)x' = f(t,x) met x(t0)=x0x(t_0) = x_0 op een zeker [t0T,t0+T][t_0 - T, t_0 + T]. Cauchy–Lipschitz (Hoofdstuk 19): is ff bovendien lokaal lipschitz in de veranderlijke xx, dan is de oplossing uniek (twee oplossingen vallen samen op hun gemeenschappelijke interval) — bestaan én eenduidigheid. Het paar (x=2x, x(0)=0)(x' = 2\sqrt{\abs x},\ x(0) = 0) scheidt de twee stellingen.

13. ω(r)=Lr\omega(r) = Lr: 01 ⁣drLr=+\int_0^1\frac{\dd r}{Lr} = +\infty, dus voldoet het. ω(r)=rlog1r\omega(r) = r\log\frac1r (nabij 00):  ⁣drrlog(1/r)=[loglog1r]+\int\frac{\dd r}{r\log(1/r)} = \bigl[-\log\log\frac1r\bigr] \to +\infty als r0r \to 0, dus voldoet het — en toch is ω(r)/r=log1r\omega(r)/r = \log\frac1r \to \infty: niet lipschitz. ω(r)=2r\omega(r) = 2\sqrt r: 01 ⁣dr2r=[r]01=1<\int_0^1\frac{\dd r}{2\sqrt r} = \bigl[\sqrt r\bigr]_0^1 = 1 < \infty, dus voldoet het niet.

14. Trek de twee integraalvormen van elkaar af:

x1(t)x2(t)=x1(s)x2(s)+st(f(v,x1(v))f(v,x2(v))) ⁣dv,x_1(t) - x_2(t) = x_1(s) - x_2(s) + \int_s^t\bigl(f(v, x_1(v)) - f(v, x_2(v))\bigr)\dd v,

neem normen en gebruik de modulus van Osgood: δ(t)δ(s)+stω(δ(v)) ⁣dv\delta(t) \leq \delta(s) + \int_s^t\omega(\delta(v))\,\dd v.

15. τ\tau is goed gedefinieerd (want δ(t0)=0\delta(t_0) = 0) met δ(τ)=0\delta(\tau) = 0 (continuïteit) en δ>0\delta > 0 op (τ,t1]\intoc\tau{t_1}. Leg s(τ,t1)s \in \intoo\tau{t_1} vast. Dan is uu C1\mathcal C^1 met u(s)=δ(s)>0u(s) = \delta(s) > 0, is uδu \geq \delta op [s,t1][s, t_1] (vraag 14), en is u=ω(δ)ω(u)u' = \omega(\delta) \leq \omega(u) (ω\omega is niet-dalend en u>0u > 0). Delen en integreren geeft

δ(s)u(t1) ⁣drω(r)=st1u(v)ω(u(v)) ⁣dvt1st1τ.\int_{\delta(s)}^{u(t_1)}\frac{\dd r}{\omega(r)} = \int_s^{t_1}\frac{u'(v)}{\omega(u(v))}\,\dd v \leq t_1 - s \leq t_1 - \tau .

Hoewel uu van ss afhangt, geldt u(t1)δ(t1)u(t_1) \geq \delta(t_1) voor elke ss, zodat het linkerlid minstens δ(s)δ(t1) ⁣drω(r)\int_{\delta(s)}^{\delta(t_1)}\frac{\dd r}{\omega(r)} is, en dat gaat naar ++\infty als sτs \downarrow \tau (dan is δ(s)δ(τ)=0\delta(s) \to \delta(\tau) = 0, en de integraal divergeert in 00): het begrensde rechterlid wordt tegengesproken. Bijgevolg is δ0\delta \equiv 0: eenduidigheid.

16. Lipschitz is ω=Lr\omega = Lr: de eenduidigheid is teruggevonden. Voor x=xlog1xx' = x\log\frac1{\abs x} voldoet het rechterlid nabij 00 aan de Osgood-modulus ω(r)=rlog1r\omega(r) = r\log\frac1r (middelwaardeongelijkheid op xxlog1xx \mapsto x\log\frac1x, waarvan de afgeleide log1x1\log\frac1x - 1 onbegrensd is — Lipschitz faalt, Osgood geldt): dus eenduidige oplossingen; en merk op dat x0x \equiv 0 er een is, zodat geen andere oplossing 00 kan raken. Voor ω(r)=2r\omega(r) = 2\sqrt r is 0h2 ⁣dr2r=h\int_0^{h^2}\frac{\dd r}{2\sqrt r} = h — het “Osgood-budget” om van 00 naar hoogte h2h^2 te klimmen is precies de tijd hh, en inderdaad is xc(c+h)=h2x_c(c + h) = h^2: de oplossing besteedt precies de tijd hh aan wat de convergente integraal toestaat. Divergentie van de integraal is de onmogelijkheid om 00 in eindige tijd te verlaten; convergentie is de ontsnappingsroute.

17. G(t)=t0t(Le(s)+η(s)) ⁣dsG(t) = \int_{t_0}^t(Le(s) + \eta(s))\dd s is C1\mathcal C^1 met G=Le+ηLG+supηG' = Le + \eta \leq LG + \sup\eta (wegens de hypothese eGe \leq G). Dan is (eLt(G+supηL))=eLt(GLGsupη)0\bigl(\eu^{-Lt}(G + \tfrac{\sup\eta}L)\bigr)' = \eu^{-Lt}(G' - LG - \sup\eta) \leq 0: de uitdrukking tussen haakjes daalt, dus G(t)+supηLeL(tt0)(G(t0)+supηL)=eL(tt0)supηLG(t) + \frac{\sup\eta}L \leq \eu^{L(t - t_0)}\bigl(G(t_0) + \frac{\sup\eta}L\bigr) = \eu^{L(t-t_0)}\frac{\sup\eta}L, oftewel e(t)G(t)supηL(eL(tt0)1)e(t) \leq G(t) \leq \frac{\sup\eta}L(\eu^{L(t-t_0)} - 1).

18. Kwantitatief defect: op (tk,tk+1)(t_k, t_{k+1}) is Δn(t)=f(tk,φn(tk))f(t,φn(t))\Delta_n(t) = f(t_k, \varphi_n(t_k)) - f(t, \varphi_n(t)) met ttkT/n\abs{t - t_k} \leq T/n en φn(t)φn(tk)MT/n\norm{\varphi_n(t) - \varphi_n(t_k)} \leq MT/n, dus ΔnLT/n+LMT/n\norm{\Delta_n} \leq L'T/n + LMT/n. Trek de integraalidentiteiten voor φn\varphi_n (vraag 5) en φ\varphi van elkaar af en gebruik de lipschitzgrens:

φn(t)φ(t)t0tLφnφ(s) ⁣ds+t0tΔn(s) ⁣ds,\norm{\varphi_n(t) - \varphi(t)} \leq \int_{t_0}^t L\,\norm{\varphi_n - \varphi}(s)\,\dd s + \int_{t_0}^t\norm{\Delta_n(s)}\,\dd s,

en vraag 17 met η=Δn\eta = \norm{\Delta_n} geeft de aangegeven grens O(1/n)O(1/n), waarbij φ\varphi uniek is volgens Cauchy–Lipschitz (of Osgood). Ook zonder snelheden: elke deelrij van de equibegrensde en equi-lipschitz rij (φn)(\varphi_n) heeft een sub-deelrij die (Ascoli plus Deel II) naar een oplossing convergeert, die de eenduidigheid dwingt gelijk te zijn aan φ\varphi; en een rij waarvan elke deelrij een sub-deelrij met dezelfde limiet heeft, convergeert.

19. Uit φn(tk)=0\varphi_n(t_k) = 0 volgt met de helling 20=02\sqrt0 = 0 dat φn(tk+1)=0\varphi_n(t_{k+1}) = 0; per inductie is φn0\varphi_n \equiv 0, convergerend naar de nuloplossing. Vanaf ε>0\varepsilon > 0: op [ε,)[\varepsilon', \infty) met ε<ε\varepsilon' < \varepsilon is de functie 2x2\sqrt x lipschitz, dus is vraag 18 van toepassing en convergeert Euler naar de unieke oplossing door (0,ε)(0, \varepsilon), namelijk x(t)=(t+ε)2x(t) = (t + \sqrt\varepsilon)^2 (controle: x=2(t+ε)=2xx' = 2(t + \sqrt \varepsilon) = 2\sqrt x). Als ε0\varepsilon \to 0 gaat (t+ε)2t2(t + \sqrt\varepsilon)^2 \to t^2 uniform op [0,1][0,1]: de dubbele limiet landt op x0(t)=t2x_0(t) = t^2 en niet op 00. Een willekeurig kleine verstoring van het beginpunt stuurt het schema van de ene oplossing naar de andere: niet-eenduidigheid gelezen als numerieke instabiliteit.

20. Niet leeg: Deel II. Elke oplossing voldoet aan x=f(t,x)M\norm{x'} = \norm{f(t, x)} \leq M: S\mathcal S is uniform MM-lipschitz en uniform begrensd (waarden in Bˉ(x0,b)\bar B(x_0, b)). Gesloten: convergeert xnSx_n \in \mathcal S uniform naar xx, ga dan tot de limiet over in xn(t)=x0+t0tf(s,xn(s)) ⁣dsx_n(t) = x_0 + \int_{t_0}^tf(s, x_n(s))\dd s (de integranden convergeren uniform wegens de uniforme continuïteit van ff op de compacte RR): dus xSx \in \mathcal S. Ascoli: S\mathcal S is een gesloten, begrensde, equicontinue deelverzameling van C([t0,t0+T],Rd)\mathcal C([t_0, t_0+T], \R^d): compact.

21. Elke oplossing is niet-dalend (x=2x0x' = 2\sqrt{\abs x} \geq 0) met x(0)=0x(0) = 0, dus x0x \geq 0. Stel c=sup{t[0,1]:x(t)=0}c = \sup\{t \in [0,1] : x(t) = 0\} (eventueel c=c = \infty als x0x \equiv 0, en dan is x=xx = x_\infty). Voor t>ct > c is x>0x > 0 (monotonie plus de definitie van cc), en daar is (x)=x2x=1(\sqrt x)' = \frac{x'}{2\sqrt x} = 1, dus x(t)=tc\sqrt{x(t)} = t - c (continuïteit in cc): x=xcx = x_c. Continuïteit van cxcc \mapsto x_c: voor c,c[0,1]c, c' \in [0, 1] is sup[0,1](tc)+2(tc)+22cc\sup_{[0,1]}\abs{(t - c)_+^2 - (t - c')_+^2} \leq 2\abs{c - c'} (de afbeelding c(tc)+2c \mapsto (t-c)_+^2 is 22-lipschitz, uniform in t[0,1]t \in [0,1]), en omdat xc0x_c \equiv 0 op [0,1][0,1] voor elke c1c \geq 1, herleidt de familie zich tot S={xc:c[0,1]}\mathcal S = \{x_c : c \in [0,1]\} (met x1=0=xx_1 = 0 = x_\infty). Dus is S\mathcal S het beeld van het compacte samenhangende [0,1][0,1] onder de continue afbeelding cxcc \mapsto x_c: compact en samenhangend. De trechter van oplossingen is een continu lijnstuk dat van t2t^2 (onmiddellijke ontsnapping) tot 00 (eeuwige rust) loopt.

22. De evaluatie evt ⁣:C([0,1])R\operatorname{ev}_t\colon \mathcal C([0,1]) \to \R, xx(t)x \mapsto x(t), is continu (x(t)y(t)xy\abs{x(t) - y(t)} \leq \norm{x - y}_\infty), dus is S(t)=evt(S)\mathcal S(t) = \operatorname{ev}_t(\mathcal S) een continu beeld van een compacte verzameling: compact — en van een samenhangende: samenhangend. Voor het voorbeeld doorloopt xc(t)=(tc)+2x_c(t) = (t - c)_+^2, terwijl cc het interval [0,1][0, 1] doorloopt, continu alle waarden van t2t^2 (bij c=0c = 0) tot 00 (bij ctc \geq t): dus S(t)=[0,t2]\mathcal S(t) = \intcc0{t^2}. Op elk tijdstip is de doorsnede van de trechter een vol lijnstuk: tussen rust en maximale ontsnapping wordt elk compromis door een echte oplossing gerealiseerd.

23. Trek de integraalvormen x(t)=x0+t0tf(s,x(s)) ⁣dsx(t) = x_0 + \int_{t_0}^tf(s, x(s))\dd s en haar tegenhanger voor yy van elkaar af, en stel e(t)=x(t)y(t)e(t) = \norm{x(t) - y(t)} en d=x0y0d = \norm{x_0 - y_0}:

e(t)d+t0tLe(s) ⁣ds=G(t).e(t) \leq d + \int_{t_0}^tL\,e(s)\dd s = G(t).

Dan is G(t0)=dG(t_0) = d en G=LeLGG' = Le \leq LG, dus (eL(tt0)G)0\bigl(\eu^{-L(t - t_0)}G\bigr)' \leq 0 en G(t)deL(tt0)G(t) \leq d\,\eu^{L(t - t_0)}; en bijgevolg e(t)deL(tt0)e(t) \leq d\,\eu^{L(t-t_0)}. Scherpte: voor x=Lxx' = Lx zijn de oplossingen door x0x_0 en y0y_0 gelijk aan x0eL(tt0)x_0\eu^{L(t - t_0)} en y0eL(tt0)y_0\eu^{L(t-t_0)}, met afstand precies deL(tt0)d\,\eu^{L(t-t_0)}. Met d=0d = 0 is e0e \equiv 0: de eenduidigheid van Cauchy–Lipschitz, in twee regels opnieuw afgeleid. En voor vaste t[t0,t0+T]t \in [t_0, t_0 + T] is x(t;x0)x(t;y0)eLTx0y0\norm{x(t; x_0) - x(t; y_0)} \leq \eu^{LT}\norm{x_0 - y_0}: de stroming is lipschitz in de beginvoorwaarde — deterministische afhankelijkheid, tegen een beheerste exponentiële prijs.

24. Op (0,1]\intoc01 is Ω\Omega goed gedefinieerd (de integraal convergeert in 00 per hypothese) en C1\mathcal C^1 met Ω=1ω>0\Omega' = \frac1\omega > 0: een stijgende bijectie op (0,Ω(1)]\intoc0{\Omega(1)}, met Ω(x)0\Omega(x) \to 0 als x0+x \to 0^+. Haar inverse g ⁣:(0,Ω(1)](0,1]g \colon \intoc0{\Omega(1)} \to \intoc01 is C1\mathcal C^1 met

g(t)=1Ω(g(t))=ω(g(t))>0,g(t)t0+0.g'(t) = \frac1{\Omega'(g(t))} = \omega\bigl(g(t)\bigr) > 0, \qquad g(t) \xrightarrow[t \to 0^+]{} 0 .

Zet gg voort met 00 voor t0t \leq 0: de continuïteit is duidelijk, en in t=0t = 0 is voor t>0t > 0

g(t)t=1t0tg(s) ⁣ds=1t0tω(g(s)) ⁣dsω(g(t))t0+0\frac{g(t)}t = \frac1t\int_0^tg'(s)\dd s = \frac1t\int_0^t\omega\bigl(g(s)\bigr)\dd s \leq \omega\bigl(g(t)\bigr) \xrightarrow[t\to0^+]{} 0

(ω\omega is niet-dalend, gg stijgend, en ω(0+)=0\omega(0^+) = 0): dus g(0)=0=f(g(0))g'(0) = 0 = f(g(0)), en g=f(g)g' = f(g) geldt aan beide kanten van 00. Bijgevolg zijn x0x \equiv 0 en gg twee verschillende C1\mathcal C^1-oplossingen door (0,0)(0, 0): convergeert 0 ⁣drω(r)\int_0\frac{\dd r}{\omega(r)}, dan faalt de eenduidigheid — de divergentie uit vraag 15 is precies de grens. Voor ω(r)=2r\omega(r) = 2\sqrt r is Ω(x)=0x ⁣dr2r=x\Omega(x) = \int_0^x\frac{\dd r}{2\sqrt r} = \sqrt x en g(t)=t2g(t) = t^2, en tijdverschuivingen geven de hele familie xcx_c uit Deel III.

25. Met stap h=1nh = \frac1n is φn(tk+1)=φn(tk)+hφn(tk)=(1+h)φn(tk)\varphi_n(t_{k+1}) = \varphi_n(t_k) + h\,\varphi_n(t_k) = (1 + h)\varphi_n(t_k), dus φn(1)=(1+1n)n\varphi_n(1) = (1 + \frac1n)^n na nn stappen. Ontwikkeling:

nlog(1+1n)=n(1n12n2+O(1n3))=112n+O(1n2),n\log\Bigl(1 + \frac1n\Bigr) = n\Bigl(\frac1n - \frac1{2n^2} + O\Bigl(\frac1{n^3}\Bigr) \Bigr) = 1 - \frac1{2n} + O\Bigl(\frac1{n^2}\Bigr),

en na exponentiëren: (1+1n)n=ee1/(2n)+O(n2)=e(112n+O(n2))(1 + \frac1n)^n = \eu\,\eu^{-1/(2n) + O(n^{-2})} = \eu\bigl(1 - \frac1{2n} + O(n^{-2})\bigr). De fout in t=1t = 1 is dus eφn(1)=e2n+O(n2)\eu - \varphi_n(1) = \frac{\eu}{2n} + O(n^{-2}): de O(1n)O(\frac1n) van vraag 18, hier met haar exacte constante e2\frac\eu2. Numeriek voor n=10n = 10: 1.110=2.59374246011.1^{10} = 2.5937424601 (1.12=1.211.1^2 = 1.21, 1.14=1.46411.1^4 = 1.4641, 1.18=2.143588811.1^8 = 2.14358881, maal 1.211.21), en e2.593740.12454\eu - 2.59374 \approx 0.12454 tegenover de asymptotische voorspelling e200.13591\frac{\eu}{20} \approx 0.13591: overeenstemming tot op de correctie O(n2)O(n^{-2}), waarvan de hoofdterm hier de voorspelling naar de waargenomen waarde toe verlaagt.

Begrippen gedefinieerd in dit hoofdstuk

Bekijk alle 395 begrippen in de begrippenlijst