Mathematics · Book 5 · Bachelor Year 3

Universitaire wiskunde — Bachelor jaar 3

Universitaire wiskunde — Bachelor jaar 3 · Bachelor Year 3

14De Fouriertransformatie

Fourierreeksen ontbinden periodieke signalen in discrete harmonischen; de Fouriertransformatie doet hetzelfde voor signalen op de hele lijn, met een continuum van frequenties. Ze zet differentiatie om in vermenigvuldiging, convolutie in producten, en Gaussen in Gaussen — de redenen waarom ze differentiaalvergelijkingen oplost, signaalverwerking aandrijft, en de centrale limietstelling in Hoofdstuk 23 zal bewijzen. Dit hoofdstuk ontwikkelt de L1L^1-theorie (Riemann–Lebesgue, inversie, injectiviteit), de Schwartz-klasse waar de transformatie een perfecte bijectie is, en de L2L^2-theorie (Plancherel: de transformatie is, op een constante na, een unitaire operator), met twee pronkstuktoepassingen: de warmtevergelijking, van begin tot eind opgelost in de weekendopgave, en de somformule van Poisson. Conventie:

f^(ξ)=Rf(x)eiξx ⁣dx.\hat f(\xi) = \int_\R f(x)\,\eu^{-\iu\xi x}\,\dd x .

14.1 De transformatie op L1L^1

Propositie 14.1

Voor fL1(R)f \in L^1(\R): f^\hat f is welgedefinieerd, begrensd (f^f1\norm{\hat f}_\infty \leq \norm f_1), continu, en:

  1. τaf^(ξ)=eiaξf^(ξ)\widehat{\tau_af}(\xi) = \eu^{-\iu a\xi}\hat f(\xi) en eiaxf^(ξ)=f^(ξa)\widehat{\eu^{\iu ax}f}(\xi) = \hat f(\xi - a);
  2. f(/λ)^(ξ)=λf^(λξ)\widehat{f(\cdot/\lambda)}(\xi) = \lambda\hat f(\lambda\xi) voor λ>0\lambda > 0;
  3. als xfL1xf \in L^1, dan is f^\hat f C1\mathcal C^1 met (f^)(ξ)=(ix)f^(ξ)(\hat f)'(\xi) = \widehat{(-\iu x)f}(\xi);
  4. als fC1f \in \mathcal C^1 met fL1f' \in L^1 (en f0f \to 0 in ±\pm\infty, hier automatisch), dan f^(ξ)=iξf^(ξ)\widehat{f'}(\xi) = \iu\xi\hat f(\xi);
  5. fg^=f^g^\widehat{f * g} = \hat f\,\hat g voor f,gL1f, g \in L^1.

Bewijs. Begrensdheid: f^f\abs{\hat f} \leq \int\abs f. Continuïteit: DCT met dominator f\abs f (Stelling 10.14). (1), (2): substituties (Stelling 11.10). (3): differentiatie onder de integraal, dominator xf\abs{xf} (Stelling 10.15). (4): eerst, f(x)=f(0)+0xff(x) = f(0) + \int_0^xf' heeft een limiet in ±\pm\infty (fL1f' \in L^1), die 00 moet zijn (fL1f \in L^1); integreer dan partieel op [A,A][-A, A] en laat AA \to \infty. (5): Fubini, legitiem omdat (x,y)f(xy)g(y)eiξx(x,y)\mapsto f(x - y)g(y)\eu^{-\iu\xi x} absoluut integreerbaar is (Stelling 11.9):

fg^(ξ)=f(xy)g(y)eiξ(xy)eiξy ⁣dx ⁣dy=f^(ξ)g^(ξ).\widehat{f*g}(\xi) = \iint f(x - y)g(y)\eu^{-\iu\xi(x - y)} \eu^{-\iu\xi y}\dd x\,\dd y = \hat f(\xi)\,\hat g(\xi).

Voorbeeld 14.2

De Gauss: voor a>0a > 0,

eax2^(ξ)=πa  eξ2/4a:\widehat{\eu^{-ax^2}}(\xi) = \sqrt{\frac\pi a}\;\eu^{-\xi^2/4a} :

via Oefening 10.7 (de ODE F=ξ2FF' = -\frac\xi{2}F truc, herschaald), of via (3): g=eax2^g = \widehat{\eu^{-ax^2}} voldoet aan g(ξ)=ξ2ag(ξ)g'(\xi) = -\frac{\xi}{2a}g(\xi) (partieel integreren), g(0)=π/ag(0) = \sqrt{\pi/a}. Gaussen zijn vaste punten van de transformatie op schaling na — de diepe reden waarom ze de centrale limietstelling regeren.

Stelling 14.3 (Riemann–Lebesgue)

Voor fL1(R)f \in L^1(\R): f^(ξ)0\hat f(\xi) \to 0 als ξ\abs\xi \to \infty. Dus f^ ⁣:L1C0(R)\widehat{\phantom f} \colon L^1 \to \mathcal C_0(\R) (continue functies die in oneindig verdwijnen).

Bewijs. Voor een indicator van een interval, 1[a,b]^(ξ)=eiaξeibξiξ0\widehat{\mathbf 1_{\intcc ab}}(\xi) = \frac{\eu^{-\iu a\xi} - \eu^{-\iu b\xi}}{\iu\xi} \to 0; dus voor stapfuncties. Stapfuncties zijn dicht in L1L^1 (Stelling 12.6(1) plus benadering van eindig-maatverzamelingen door eindige unies van intervallen, Oefening 9.7), en de transformatie is \norm\cdot_\infty-1\norm\cdot_1 continu: voor fs1<ε\norm{f - s}_1 < \varepsilon, lim supξf^(ξ)ε\limsup_{\abs\xi\to\infty}\abs{\hat f(\xi)} \leq \varepsilon.

14.2 Inversie en injectiviteit

Lemma 14.4 (Vermenigvuldigingsformule)

Voor f,gL1(R)f, g \in L^1(\R): f^g=fg^\displaystyle\int \hat f\,g = \int f\,\hat g.

Bewijs. Beide zijden zijn gelijk aan f(x)g(ξ)eixξ ⁣dx ⁣dξ\iint f(x)g(\xi)\eu^{-\iu x\xi}\dd x\,\dd\xi (Tonelli–Fubini: de dubbele integraal van de absolute waarde is f1g1\norm f_1\norm g_1).

Stelling 14.5 (Inversie)

Zij fL1(R)f \in L^1(\R).

  1. (Gaussiaanse sommeerbaarheid) Voor elke xx,

    (fgε)(x)=12πRf^(ξ)eεξ2eixξ ⁣dξ,waar gε(y)=12πεey2/4ε,(f * g_\varepsilon)(x) = \frac1{2\pi}\int_\R \hat f(\xi)\, \eu^{-\varepsilon\xi^2}\,\eu^{\iu x\xi}\,\dd\xi, \qquad\text{waar } g_\varepsilon(y) = \frac{1}{2\sqrt{\pi\varepsilon}}\, \eu^{-y^2/4\varepsilon},

    en fgεff * g_\varepsilon \to f in L1L^1 als ε0\varepsilon \to 0.

  2. Als bovendien f^L1\hat f \in L^1, dan voor bijna elke xx

    f(x)=12πRf^(ξ)eixξ ⁣dξ,f(x) = \frac{1}{2\pi}\int_\R \hat f(\xi)\,\eu^{\iu x\xi}\,\dd\xi ,

    en ff heeft een continue representant.

  3. (Injectiviteit) Als f^=0\hat f = 0 dan f=0f = 0 a.e.

Bewijs. (1) Fixeer xx en pas Lemma 14.4 toe op ff en g(ξ)=12πeεξ2eixξg(\xi) = \frac1{2\pi}\eu^{-\varepsilon\xi^2}\eu^{\iu x\xi}: via Voorbeeld 14.2 (met de modulatieregel),

g^(y)=12ππεe(yx)2/4ε=gε(xy),\hat g(y) = \frac1{2\pi}\sqrt{\frac\pi\varepsilon}\, \eu^{-(y - x)^2/4\varepsilon} = g_\varepsilon(x - y),

dus 12πf^(ξ)eεξ2eixξ ⁣dξ=f(y)gε(xy) ⁣dy=(fgε)(x)\frac1{2\pi}\int\hat f(\xi)\eu^{-\varepsilon\xi^2} \eu^{\iu x\xi}\dd\xi = \int f(y)g_\varepsilon(x - y)\dd y = (f*g_\varepsilon)(x). De gεg_\varepsilon zijn een benaderende eenheid: gε0g_\varepsilon \geq 0, gε=1\int g_\varepsilon = 1 (Gauss-integraal), concentrerend in 00; het bewijs van Stelling 12.9(2) geldt woordelijk (alleen gε=1\int g_\varepsilon = 1 en concentratie werden gebruikt: voor de staart, y>δgε0\int_{\abs y > \delta}g_\varepsilon \to 0): fgεf10\norm{f * g_\varepsilon - f}_1 \to 0.

(2) Als f^L1\hat f \in L^1: de rechterkant van (1) convergeert, via DCT (dominator f^\abs{\hat f}), naar 12πf^(ξ)eixξ ⁣dξ\frac1{2\pi}\int\hat f(\xi)\eu^{\iu x\xi}\dd\xi voor elke xx, en deze limietfunctie is continu (DCT opnieuw). Anderzijds fgεff * g_\varepsilon \to f in L1L^1, dus langs een deelrij a.e. (Stelling 12.4): de twee limieten stemmen a.e. overeen.

(3) f^=0\hat f = 0 maakt de rechterkant van (1) verdwijnen: fgε=0f * g_\varepsilon = 0 voor alle ε\varepsilon, en fgεff * g_\varepsilon \to f in L1L^1: f=0f = 0 a.e.

14.3 De Schwartz-klasse

Definitie 14.6

De Schwartz-klasse S(R)\mathcal S(\R) bestaat uit de C\mathcal C^\infty-functies ff met supxxmf(n)(x)<\sup_x\abs{x^m f^{(n)}(x)} < \infty voor alle m,n0m, n \geq 0 (alle afgeleiden vervallen sneller dan elke macht). Voorbeelden: eax2\eu^{-ax^2}, Cc\mathcal C_c^\infty. Duidelijk SLp\mathcal S \subseteq L^p voor elke pp (grens door C(1+x2)1C(1 + x^2)^{-1}), en S\mathcal S is stabiel onder afgeleiden, vermenigvuldiging met polynomen, en producten.

Stelling 14.7

De Fouriertransformatie beeldt S(R)\mathcal S(\R) bijectief op zichzelf af, met inverse gˇ(x)=12πg(ξ)eixξ ⁣dξ\check g(x) = \frac1{2\pi}\int g(\xi)\eu^{\iu x\xi}\dd\xi.

Bewijs. Zij fSf \in \mathcal S. Itereren van Propositie 14.1(3), (f^)(n)=(ix)nf^(\hat f)^{(n)} = \widehat{(-\iu x)^nf} (elke xkfL1x^kf \in L^1); itereren van (4) met h=(ix)nfSh = (-\iu x)^nf \in \mathcal S (al wier afgeleiden integreerbaar zijn), (iξ)mh^=h(m)^(\iu\xi)^m\hat h = \widehat{h^{(m)}}. Combineren,

ξm(f^)(n)(ξ)=h(m)^(ξ)((ix)nf)(m)1<\abs{\xi^m\,(\hat f)^{(n)}(\xi)} = \bigl|\widehat{\,h^{(m)}}(\xi)\bigr| \leq \bigl\|\bigl((-\iu x)^nf\bigr)^{(m)}\bigr\|_1 < \infty

uniform in ξ\xi: f^S\hat f \in \mathcal S. Omdat f^L1\hat f \in L^1, geldt inversie (Stelling 14.5(2)) overal (beide zijden continu): f^ˇ=f\check{\hat f} = f, en symmetrisch gˇ^=g\widehat{\check g} = g (de check-transformatie is g12πg^()g \mapsto \frac1{2\pi}\hat g(-\cdot), opnieuw S\mathcal S bewarend): bijectie.

14.4 Plancherel en L2L^2

Stelling 14.8 (Plancherel)

Voor fL1L2(R)f \in L^1 \cap L^2(\R):

f^22=2πf22.\norm{\hat f}_2^2 = 2\pi\,\norm f_2^2 .

Bijgevolg strekt f^\widehat{\phantom f} zich uniek uit tot een continue lineaire afbeelding F ⁣:L2(R)L2(R)\mathcal F \colon L^2(\R) \to L^2(\R) met Ff2=2πf2\norm{\mathcal Ff}_2 = \sqrt{2\pi}\norm f_2; F\mathcal F is bijectief, met F1=12πFσ\mathcal F^{-1} = \frac1{2\pi}\,\mathcal F\circ\sigma waar σf=f()\sigma f = f(-\cdot), en ze bewaart inwendige producten op de factor 2π2\pi na.

Bewijs. Zij fL1L2f \in L^1\cap L^2 en h=ff~h = f * \tilde f met f~(x)=f(x)\tilde f(x) = \overline{f(-x)}. Dan hL1h \in L^1 (Stelling 11.9), hh is continu en begrensd (Oefening 12.6: f,f~L2f, \tilde f \in L^2), h(0)=ffˉ=f22h(0) = \int f\bar f = \norm f_2^2, en h^=f^f~^=f^f^=f^20\hat h = \hat f\,\widehat{\tilde f} = \hat f\,\overline{\hat f} = \abs{\hat f}^2 \geq 0 (bereken f~^=f^\widehat{\tilde f} = \overline{\hat f}). Pas Stelling 14.5(1) toe op hh in x=0x = 0:

(hgε)(0)=12πh^(ξ)eεξ2 ⁣dξ.(h * g_\varepsilon)(0) = \frac1{2\pi}\int \hat h(\xi)\,\eu^{-\varepsilon\xi^2}\dd\xi .

Als ε0\varepsilon \to 0: de linkerkant gaat naar h(0)h(0) (hh continu begrensd: (hgε)(0)h(0)=(h(y)h(0))gε(y) ⁣dy0(h*g_\varepsilon)(0) - h(0) = \int(h(-y) - h(0))g_\varepsilon(y)\dd y \to 0 door splitsen klein/groot yy); de rechterkant stijgt naar 12πh^\frac1{2\pi}\int\hat h via MCT (h^0\hat h \geq 0). Dus 12πf^2=f22\frac1{2\pi}\int\abs{\hat f}^2 = \norm f_2^2, eindig of niet a priori — en eindig, bewijzend zowel lidmaatschap als de identiteit.

Uitbreiding: L1L2CcL^1\cap L^2 \supseteq \mathcal C_c is dicht in L2L^2 (Stelling 12.6); de transformatie is 2π\sqrt{2\pi}-isometrisch daar, dus strekt uniek uit tot een isometrie-op-constante-na F\mathcal F op L2L^2 (Stelling 7.2). Inversie voor S\mathcal S (Stelling 14.7) transporteert via dezelfde dichtheid (beide zijden L2L^2-continu): F(12πF(σf))=f\mathcal F\bigl(\frac1{2\pi}\mathcal F(\sigma f)\bigr) = f op S\mathcal S, dus op L2L^2: bijectiviteit. Inwendige producten: polarisatie uit de normidentiteit.

Stelling 14.9 (Somformule van Poisson)

Zij fS(R)f \in \mathcal S(\R) (continue ff met f+f^C(1+)2\abs{f} + \abs{\hat f} \leq C(1 + \abs\cdot)^{-2} volstaat). Dan

nZf(n)  =  kZf^(2πk).\sum_{n\in\Z} f(n) \;=\; \sum_{k\in\Z}\hat f(2\pi k) .

Bewijs. Zij F(x)=nZf(x+n)F(x) = \sum_{n\in\Z}f(x + n): de reeks convergeert normaal op compacta (verval van ff), dus FF is continu, en ze is 11-periodiek. Haar Fouriercoëfficiënten (periode 11: ck(F)=01F(t)e2iπkt ⁣dtc_k(F) = \int_0^1F(t)\eu^{-2\iu\pi kt}\dd t):

ck(F)=n01f(t+n)e2iπkt ⁣dt=Rf(t)e2iπkt ⁣dt=f^(2πk)c_k(F) = \sum_n\int_0^1 f(t + n)\,\eu^{-2\iu\pi kt}\dd t = \int_\R f(t)\,\eu^{-2\iu\pi kt}\dd t = \hat f(2\pi k)

(normale convergentie rechtvaardigt de verwisseling; de fase is 11-periodiek). De reeks kck(F)\sum_k\abs{c_k(F)} convergeert (verval van f^\hat f), dus de Fourierreeks van FF convergeert normaal; haar som is een continue functie met dezelfde Fouriercoëfficiënten als FF, dus gelijk aan FF (injectiviteit op de cirkel: het verschil heeft nule coëfficiënten, en Stelling 13.9 geeft nul in L2L^2, dus overal door continuïteit). Evalueer in x=0x = 0.

Voorbeeld 14.10 (De theta-identiteit)

Poisson toepassen op f(x)=eπtx2f(x) = \eu^{-\pi tx^2} (t>0t > 0), wier transform f^(ξ)=t1/2eξ2/4πt\hat f(\xi) = t^{-1/2}\eu^{-\xi^2/4\pi t} is (Voorbeeld 14.2 met a=πta = \pi t):

nZeπn2t=1tkZeπk2/t:\sum_{n\in\Z}\eu^{-\pi n^2t} = \frac1{\sqrt t}\sum_{k\in\Z}\eu^{-\pi k^2/t} :

de functionaalvergelijking van Jacobi’s theta-functie, sleutel tot de functionaalvergelijking van Riemanns ζ\zeta — en een spectaculaire numerieke accelerator: voor tt klein convergeert de linkerkant langzaam, de rechterkant bliksemsnel.

Methode 14.11

Werkbereiken: L1L^1 — transform puntsgewijs gedefinieerd, inversie heeft f^L1\hat f \in L^1 nodig; S\mathcal S — alles is legaal, bewijs hier eerst; L2L^2 — transform gedefinieerd via dichtheid (niet via de integraal!), perfecte symmetrie, Parseval-boekhouding. Om een transform te berekenen: reduceer tot de tabel (indicator, exponentiaal, Gauss) via de regels van Propositie 14.1; om een identiteit te bewijzen: vestig ze op S\mathcal S (of Cc\mathcal C_c^\infty) en strek uit via dichtheid en continuïteit (Methode 12.13); om een lineaire PDE of ODE met constante coëfficiënten op te lossen: transformeer, deel, inverteer.

De warmtekern g_t(x) = 12√π t\, -x2/4t op drie tijden: totale massa 1 voor altijd, hoogte t-1/2, breedte √ t. Convolueren van beginwaarden met deze spreidende Gauss is de hele inhoud van de weekendopgave; in frequentie leest hetzelfde plaatje g_t( ) = -t 2 — hoge frequenties sterven eerst, en die asymmetrie is de pijl van de tijd.
De warmtekern gt(x)=12πtex2/4tg_t(x) = \frac1{2\sqrt{\pi t}}\,\eu^{-x^2/4t} op drie tijden: totale massa 11 voor altijd, hoogte t1/2\sim t^{-1/2}, breedte t\sim \sqrt t. Convolueren van beginwaarden met deze spreidende Gauss is de hele inhoud van de weekendopgave; in frequentie leest hetzelfde plaatje g^t(ξ)=etξ2\hat g_t(\xi) = \eu^{-t\xi^2} — hoge frequenties sterven eerst, en die asymmetrie is de pijl van de tijd.

14.5 Oefeningen

Oefening 14.1

Bereken de Fouriertransformaties van: 1[a,a]\mathbf 1_{\intcc{-a}a}; eax\eu^{-a\abs x} (a>0a > 0); de tentfunctie max(0,1x)\max(0, 1 - \abs x); 1x2+a2\frac1{x^2 + a^2} (gebruik inversie op de tweede). Noteer de opkomende tabel.

Oplossing

Oplossing van Oefening 14.1.

1[a,a]^(ξ)=aaeiξx ⁣dx=2sin(aξ)ξ\widehat{\mathbf 1_{\intcc{-a}a}}(\xi) = \int_{-a}^a\eu^{-\iu\xi x}\dd x = \frac{2\sin(a\xi)}{\xi} (waarde 2a2a in 00). eax^(ξ)=0e(a+iξ)x+e(aiξ)x ⁣dx=1a+iξ+1aiξ=2aa2+ξ2\widehat{\eu^{-a\abs x}}(\xi) = \int_0^\infty\eu^{-(a + \iu\xi)x} + \eu^{-(a - \iu\xi)x}\,\dd x = \frac1{a + \iu\xi} + \frac1{a - \iu\xi} = \frac{2a}{a^2 + \xi^2}. Tent: max(0,1x)=1[1/2,1/2]1[1/2,1/2]\max(0, 1 - \abs x) = \mathbf 1_{\intcc{-1/2}{1/2}} * \mathbf 1_{\intcc{-1/2}{1/2}}, dus haar transform is (2sin(ξ/2)ξ)2=(sin(ξ/2)ξ/2)2\bigl(\frac{2\sin(\xi/2)}\xi\bigr)^2 = \bigl(\frac{\sin(\xi/2)}{\xi/2}\bigr)^2. Laatste: 2aa2+ξ2L1\frac{2a}{a^2+\xi^2} \in L^1, dus inversie (Stelling 14.5(2)) toegepast op eax\eu^{-a\abs x} geeft, na hernoemen van variabelen,

(1x2+a2)^(ξ)=πaeaξ.\widehat{\Bigl(\frac1{x^2 + a^2}\Bigr)}(\xi) = \frac{\pi}{a}\,\eu^{-a\abs\xi} .

Oefening 14.2

Zij fL1f \in L^1. Druk uit in termen van f^\hat f de transformaties van: f(xa)f(x - a), f(x)cos(bx)f(x)\cos(bx), f(ax+b)f(ax + b), f(x)\overline{f(-x)}, (ff)(x)(f * f)(x). Verifieer elke regel op de Gauss.

Oplossing

Oplossing van Oefening 14.2.

Uit Propositie 14.1: f(a)^=eiaξf^(ξ)\widehat{f(\cdot - a)} = \eu^{-\iu a\xi}\hat f(\xi); fcos(b)^=12(f^(ξb)+f^(ξ+b))\widehat{f\cos(b\cdot)} = \frac12\bigl(\hat f(\xi - b) + \hat f(\xi + b)\bigr); f(a+b)^(ξ)=1aeibξ/af^(ξ/a)\widehat{f(a\cdot + b)}(\xi) = \frac1a\,\eu^{\iu b\xi/a}\,\hat f(\xi/a) (a>0a > 0); f()^=f^\widehat{\overline{f(-\cdot)}} = \overline{\hat f}; ff^=f^2\widehat{f * f} = \hat f^2. Op de Gauss (ex2^=πeξ2/4\widehat{\eu^{-x^2}} = \sqrt\pi\eu^{-\xi^2/4}) is elke regel een eenregelige controle — b.v. e(xa)2\eu^{-(x-a)^2} heeft transform πeiaξeξ2/4\sqrt\pi\,\eu^{-\iu a\xi}\eu^{-\xi^2/4}, die de directe berekening (maak het kwadraat volledig) bevestigt.

Oefening 14.3 ★★

(a) Toon dat 1[1,1]1[1,1]\mathbf 1_{\intcc{-1}1} * \mathbf 1_{\intcc{-1}1} transform (2sinξξ)2\bigl(\frac{2\sin\xi}\xi \bigr)^2 heeft, en concludeer R(sinξξ)2 ⁣dξ=π\int_\R\bigl(\frac{\sin\xi}\xi\bigr)^2\dd\xi = \pi via Plancherel — of via inversie in 00. Vergelijk Probleem 10.1. (b) Bereken R ⁣dx(x2+1)2\int_\R\frac{\dd x}{(x^2+1)^2} via Plancherel toegepast op ex\eu^{-\abs x}.

Oplossing

Oplossing van Oefening 14.3.

(a) h=1[1,1]1[1,1]h = \mathbf 1_{\intcc{-1}1}*\mathbf 1_{\intcc{-1}1} heeft h^=(2sinξξ)2L1\hat h = \bigl(\frac{2\sin\xi}\xi\bigr)^2 \in L^1; inversie in x=0x = 0, waar h(0)=λ([1,1][1,1])=2h(0) = \lambda(\intcc{-1}1\cap\intcc{-1}1) = 2:

2=12πR(2sinξξ)2 ⁣dξ  R(sinξξ)2 ⁣dξ=π,2 = \frac1{2\pi}\int_\R\Bigl(\frac{2\sin\xi}\xi\Bigr)^2 \dd\xi \ \Longrightarrow\ \int_\R\Bigl(\frac{\sin\xi}\xi\Bigr)^2\dd\xi = \pi ,

consistent met 0sin2ξ2=π2\int_0^\infty\frac{\sin^2}{\xi^2} = \frac\pi2 (Probleem 10.1).

(b) Plancherel voor f=exf = \eu^{-\abs x}: f^2=2πf2\int\abs{\hat f}^2 = 2\pi\int\abs f^2 leest 4 ⁣dξ(1+ξ2)2=2πe2x ⁣dx=2π\int\frac{4\,\dd\xi}{(1 + \xi^2)^2} = 2\pi\int\eu^{-2\abs x}\dd x = 2\pi: R ⁣dξ(1+ξ2)2=π2\int_\R\frac{\dd\xi}{(1+\xi^2)^2} = \frac\pi2.

Oefening 14.4 ★★

(Warmtekern-algebra) Met gt(x)=12πtex2/4tg_t(x) = \frac1{2\sqrt{\pi t}}\eu^{-x^2/4t}: (a) verifieer g^t(ξ)=etξ2\hat g_t(\xi) = \eu^{-t\xi^2}; (b) concludeer de semigroepwet gtgs=gt+sg_t * g_s = g_{t+s} zonder enige integraalberekening; (c) toon gt1=1\norm{g_t}_1 = 1 en gt22=(8πt)1/2\norm{g_t}_2^2 = (8\pi t)^{-1/2}.

Oplossing

Oplossing van Oefening 14.4.

(a) gt(x)=12πtex2/4tg_t(x) = \frac1{2\sqrt{\pi t}}\eu^{-x^2/4t}: via Voorbeeld 14.2 met a=14ta = \frac1{4t}, g^t(ξ)=12πt4πtetξ2=etξ2\hat g_t(\xi) = \frac1{2\sqrt{\pi t}}\sqrt{4\pi t}\,\eu^{-t\xi^2} = \eu^{-t\xi^2}. (b) gtgs^=etξ2esξ2=gt+s^\widehat{g_t * g_s} = \eu^{-t\xi^2}\eu^{-s\xi^2} = \widehat{g_{t+s}}, en de transform is injectief op L1L^1 (Stelling 14.5(3)): gtgs=gt+sg_t * g_s = g_{t+s}. (c) gt1=1\norm{g_t}_1 = 1 (Gauss-integraal); gt22=14πtex2/2t ⁣dx=2πt4πt=18πt\norm{g_t}_2^2 = \frac1{4\pi t}\int\eu^{-x^2/2t}\dd x = \frac{\sqrt{2\pi t}}{4\pi t} = \frac1{\sqrt{8\pi t}}.

Oefening 14.5 ★★

Toon dat als fL1f \in L^1 even en reëel is, f^\hat f even en reëel is; als ff oneven en reëel is, is f^\hat f oneven en zuiver imaginair. Wat berekent f^(0)\hat f(0)? Concludeer dat f0f \geq 0 f^=f^(0)=f\norm{\hat f}_\infty = \hat f(0) = \int f forceert, en interpreteer voor kansdichtheden (Hoofdstuk 23: een karakteristieke functie heeft modulus 1\leq 1, bereikt in 00).

Oplossing

Oplossing van Oefening 14.5.

Voor reële even ff: f^(ξ)=fcos(ξx) ⁣dx\hat f(\xi) = \int f\cos(\xi x)\dd x (het sinusdeel cancelleert): reëel en even. Oneven: f^(ξ)=ifsin(ξx)\hat f(\xi) = -\iu\int f\sin(\xi x): oneven, zuiver imaginair. f^(0)=f\hat f(0) = \int f: de totale massa. Als f0f \geq 0: f^(ξ)f=f=f^(0)\abs{\hat f(\xi)} \leq \int\abs f = \int f = \hat f(0), dus het sup wordt bereikt in 00. Voor een kansdichtheid is f^(ξ)\hat f(-\xi) de karakteristieke functie van Hoofdstuk 23: modulus 1\leq 1 overal, =1= 1 in de oorsprong.

Oefening 14.6 ★★★

(Niet-surjectiviteit) Toon dat f^ ⁣:L1C0\widehat{\phantom f}\colon L^1 \to \mathcal C_0 injectief en continu is, maar niet surjectief, in drie stappen. (i) Injectiviteit (Stelling 14.5) en continuïteit (f^f1\norm{\hat f}_\infty \leq \norm f_1), en C0\mathcal C_0 is een Banachruimte (gesloten in \norm\cdot_\infty). (ii) Als de afbeelding surjectief was, zou ze bijectief zijn, en de open-afbeeldingsstelling (Stelling 8.12) zou een constante CC geven met f1Cf^\norm f_1 \leq C\norm{\hat f}_\infty voor alle fL1f \in L^1. (iii) Weerspreek dit met fn(x)=sinxxsin(x/n)x/nf_n(x) = \frac{\sin x}{x}\cdot\frac{\sin(x/n)}{x/n}: haar transform is (op constanten na) de convolutie 1[1,1]1[1/n,1/n]\mathbf 1_{\intcc{-1}1} * \mathbf 1_{\intcc{-1/n}{1/n}}-type trapezium — toon f^nπ\norm{\hat f_n}_\infty \leq \pi uniform, terwijl fn1clnn\norm{f_n}_1 \geq c\ln n door de bogen van sinxx\frac{\abs{\sin x}}x op [1,n][1, n] te tellen (waar de tweede factor van onderen begrensd is), zoals in Stelling 8.11.

Oplossing

Oplossing van Oefening 14.6.

(i) Injectiviteit is Stelling 14.5(3); continuïteit is f^f1\norm{\hat f}_\infty \leq \norm f_1 (met waarden in C0\mathcal C_0 via Riemann–Lebesgue); C0\mathcal C_0 is gesloten in de sup-norm (uniforme limieten van verdwijnende-in-oneindig-functies verdwijnen in oneindig): Banach.

(ii) Een continue bijectie tussen Banachruimten heeft continue inverse (Stelling 8.12): er zou CC zijn met f1Cf^\norm f_1 \leq C\norm{\hat f}_\infty.

(iii) Zij fn(x)=sinxxsin(x/n)x/nf_n(x) = \frac{\sin x}x\cdot\frac{\sin(x/n)}{x/n}: een product van twee L2L^2-functies, en O(x2)O(x^{-2}) in oneindig, dus fnL1L2f_n \in L^1\cap L^2. Omdat (sin(ax)ax)\bigl(\frac{\sin(ax)}{ax}\bigr) L2L^2-transform πa1[a,a]\frac\pi a\mathbf 1_{\intcc{-a}a} heeft, geeft de productformule gh^=12πg^h^\widehat{gh} = \frac1{2\pi}\hat g * \hat h (geldig voor g,hL2g, h \in L^2 met ghL1gh \in L^1; controleer op Schwartz-functies via Fubini en strek uit via L2L^2-continuïteit van beide zijden via Plancherel)

f^n=12π(π1[1,1])(πn1[1/n,1/n]):\hat f_n = \frac1{2\pi}\,\bigl(\pi\mathbf 1_{\intcc{-1}1}\bigr) * \bigl(\pi n\,\mathbf 1_{\intcc{-1/n}{1/n}}\bigr):

een trapezium van hoogte πn22n=π\frac{\pi n}2\cdot\frac2n = \pi: f^n=π\norm{\hat f_n}_\infty = \pi voor elke nn. Maar op [1,n][1, n], sin(x/n)x/nsin1>0\frac{\sin(x/n)}{x/n} \geq \sin 1 > 0, dus

fn1sin11nsinxx ⁣dxclnn\norm{f_n}_1 \geq \sin 1\int_1^n\frac{\abs{\sin x}}x\dd x \geq c\ln n

(bogentelling, zoals in Stelling 8.11). De grens fn1Cπ\norm{f_n}_1 \leq C\pi faalt voor grote nn: niet surjectief. (Het beeld is een dichte — via Stone–Weierstrass-type argumenten — maar echte deelruimte van C0\mathcal C_0.)

Oefening 14.7 ★★

(Gladheid \leftrightarrow verval-woordenboek) Bewijs: fL1f \in L^1 met f^(ξ)=O(ξk1δ)\hat f(\xi) = O(\abs\xi^{-k-1-\delta}) voor zekere δ>0\delta > 0 impliceert dat ff een Ck\mathcal C^k-representant heeft. Omgekeerd fCckf \in \mathcal C^k_c impliceert f^(ξ)=O(ξk)\hat f(\xi) = O(\abs\xi^{-k}). Illustreer beide richtingen op de tentfunctie.

Oplossing

Oplossing van Oefening 14.7.

Als f^(ξ)=O(ξk1δ)\hat f(\xi) = O(\abs\xi^{-k-1-\delta}): dan ξjf^L1\xi^j\hat f \in L^1 voor 0jk0 \leq j \leq k (integreerbaar in oneindig via het verval, lokaal via continuïteit van f^\hat f). Inversie (Stelling 14.5(2)) representeert ff a.e. door x12πf^(ξ)eixξ ⁣dξx \mapsto \frac1{2\pi}\int\hat f(\xi)\eu^{\iu x\xi}\dd\xi, en differentiatie onder de integraal (dominators ξjf^\abs{\xi^j\hat f}) maakt deze representant Ck\mathcal C^k. Omgekeerd voor fCckf \in \mathcal C_c^k: itereren van Propositie 14.1(4), (iξ)kf^=f(k)^(\iu\xi)^k\hat f = \widehat{f^{(k)}}, dus f^f(k)1ξk\abs{\hat f} \leq \norm{f^{(k)}}_1\abs\xi^{-k}. Tentfunctie: continu met compacte drager (k=0k = 0: transform begrensd), en haar transform ξ2=O(ξ011)\sim \xi^{-2} = O(\abs\xi^{-0-1-1}) geeft, via de eerste richting, een C0\mathcal C^0-representant — beide scherp: de tent is niet C1\mathcal C^1, en haar transform vervalt niet sneller dan ξ2\xi^{-2}.

Oefening 14.8 ★★★

(Ongelijkheid van Heisenberg) Voor fS(R)f \in \mathcal S(\R) reëel met f2=1\norm f_2 = 1, bewijs

(x2f(x)2 ⁣dx)(12πξ2f^(ξ)2 ⁣dξ)    14,\Bigl(\int x^2f(x)^2\dd x\Bigr)\cdot \Bigl(\frac1{2\pi}\int \xi^2\abs{\hat f(\xi)}^2\dd\xi\Bigr) \;\geq\; \frac14 ,

met gelijkheid voor Gaussen. (Schrijf 1=f2=x(f2)1 = \int f^2 = -\int x\,(f^2)' partieel, begrens door Cauchy–Schwarz, en converteer f2\norm{f'}_2 via Plancherel.) Interpretatie: een signaal en zijn spectrum kunnen niet beide geconcentreerd zijn.

Oplossing

Oplossing van Oefening 14.8.

Partieel integreren (fSf \in \mathcal S; randtermen verdwijnen):

1=f2=[xf2]x(f2)=2xff2xf2f2.1 = \int f^2 = \bigl[xf^2\bigr]_{-\infty}^{\infty} - \int x\,(f^2)' = -2\int xff' \leq 2\,\norm{xf}_2\,\norm{f'}_2 .

Plancherel en f^=iξf^\widehat{f'} = \iu\xi\hat f: f22=12πξ2f^2\norm{f'}_2^2 = \frac1{2\pi}\int\xi^2\abs{\hat f}^2. Kwadrateer de weergave:

14xf2212πξ2f^2 ⁣dξ.\frac14 \leq \norm{xf}_2^2\cdot\frac1{2\pi} \int\xi^2\abs{\hat f}^2\dd\xi .

Gelijkheid vereist gelijkheid in Cauchy–Schwarz: f=λxff' = \lambda xf met λ<0\lambda < 0 (integreerbaarheid), d.w.z. f(x)=ceλx2/2f(x) = c\,\eu^{\lambda x^2/2}: Gaussen. Een signaal geconcentreerd in xx (kleine xf2\norm{xf}_2) moet uitgespreid spectrum hebben, en omgekeerd: het onzekerheidsprincipe.

Oefening 14.9 ★★

Rechtvaardig Voorbeeld 14.10 in detail (hypothesen van Poisson voor de Gauss), en gebruik de identiteit om nZeπn2\sum_{n\in\Z}\eu^{-\pi n^2} tot zes decimalen met drie termen te evalueren. Hoeveel termen van de definiërende reeks zou dezelfde nauwkeurigheid vereisen bij t=102t = 10^{-2}, versus de getransformeerde reeks?

Oplossing

Oplossing van Oefening 14.9.

De Gauss f(x)=eπtx2f(x) = \eu^{-\pi tx^2} is Schwartz, dus Stelling 14.9 geldt, en f^(ξ)=t1/2eξ2/4πt\hat f(\xi) = t^{-1/2}\eu^{-\xi^2/4\pi t}; in ξ=2πk\xi = 2\pi k wordt de rechterkant t1/2eπk2/tt^{-1/2}\eu^{-\pi k^2/t}: de theta-identiteit. Bij t=1t = 1:

nZeπn2=1+2eπ+2e4π+1+0.0864278+0.0000070=1.0864348,\sum_{n\in\Z}\eu^{-\pi n^2} = 1 + 2\eu^{-\pi} + 2\eu^{-4\pi} + \cdots \approx 1 + 0.0864278 + 0.0000070 = 1.0864348,

nauwkeurig tot 66 decimalen met drie termen (e9π51013\eu^{-9\pi} \approx 5\cdot10^{-13}). Bij t=102t = 10^{-2}: de definiërende reeks heeft eπn2/100<107\eu^{-\pi n^2/100} < 10^{-7} nodig, d.w.z. n23n \gtrsim 23 — ongeveer 4747 termen — terwijl de getransformeerde reeks 10ke100πk210\sum_k\eu^{-100\pi k^2} is, waar al de k=1k = 1-term 10136\sim 10^{-136} is: één term volstaat.

Oefening 14.10 ★★

(Bandbegrensde functies) Zij fL2(R)f \in L^2(\R) met Ff\mathcal Ff gedragen in [π,π]\intcc{-\pi}\pi. Toon dat ff een representant heeft die zich uitstrekt waarvan elke waarde herwinbaar is uit samples: bewijs de Shannon-interpolatie in de gehele getallen,

f(x)=nZf(n)sin(π(xn))π(xn)in L2,f(x) = \sum_{n\in\Z} f(n)\, \frac{\sin\bigl(\pi(x - n)\bigr)}{\pi(x - n)} \quad\text{in } L^2,

door Ff\mathcal Ff te ontwikkelen in de Fourier-basis van L2([π,π])L^2(\intcc{-\pi}\pi) (Stelling 13.9) en term-voor-term terug te transformeren.

Oplossing

Oplossing van Oefening 14.10.

FfL2([π,π])L1([π,π])\mathcal Ff \in L^2(\intcc{-\pi}\pi) \subseteq L^1(\intcc{-\pi}\pi) (eindige maat), dus inversie geeft de continue representant f(x)=12πππFf(ξ)eixξ ⁣dξf(x) = \frac1{2\pi}\int_{-\pi}^\pi\mathcal Ff(\xi)\eu^{\iu x\xi}\dd\xi, met

f(n)=12πππFf(ξ)einξ ⁣dξ=en,Fff(n) = \frac1{2\pi}\int_{-\pi}^{\pi}\mathcal Ff(\xi)\,\eu^{\iu n\xi}\dd\xi = \langle e_{-n}, \mathcal Ff\rangle

in de notatie van Stelling 13.9. Ontwikkelen in die Hilbert-basis: Ff=nf(n)einξ\mathcal Ff = \sum_nf(n)\,\eu^{-\iu n\xi} in L2([π,π])L^2(\intcc{-\pi}\pi). Pas de L2L^2-continue F1\mathcal F^{-1} term-voor-term toe:

F1(1[π,π]einξ)(x)=12πππeiξ(xn) ⁣dξ=sin(π(xn))π(xn),\mathcal F^{-1}\bigl(\mathbf 1_{\intcc{-\pi}\pi}\eu^{-\iu n\xi}\bigr)(x) = \frac1{2\pi}\int_{-\pi}^{\pi}\eu^{\iu\xi(x - n)}\dd\xi = \frac{\sin\bigl(\pi(x-n)\bigr)}{\pi(x - n)} ,

gevend f=nf(n)sinc(n)f = \sum_nf(n)\operatorname{sinc}(\cdot - n) in L2L^2: een bandbegrensd signaal is bepaald door haar geheeltallige samples — Shannons samplingstelling.

Oefening 14.11 ★★

(De transform als operator van orde vier) Op S(R)\mathcal S(\R), zij Ff=f^\mathcal F f = \hat f. (a) Met de inversieformule, toon (F2f)(x)=2πf(x)(\mathcal F^2f)(x) = 2\pi\,f(-x), en concludeer F4=(2π)2id\mathcal F^4 = (2\pi)^2\, \mathrm{id}. (b) Concludeer dat elke eigenwaarde van F\mathcal F op S\mathcal S tot {±2π,±i2π}\{\pm\sqrt{2\pi}, \pm\iu\sqrt{2\pi}\} behoort, en exposeer een eigenfunctie voor +2π+\sqrt{2\pi} (welke functie van dit hoofdstuk is evenredig met haar eigen transform?). (c) Toon dat even functies F2f=2πf\mathcal F^2f = 2\pi f en oneven F2f=2πf\mathcal F^2f = -2\pi f voldoen; produceer een eigenfunctie voor de eigenwaarde i2π-\iu\sqrt{2\pi} uit xex2/2x\eu^{-x^2/2} door haar transform te berekenen (differentiëer de transform van de Gauss).

Oplossing

Oplossing van Oefening 14.11.

(a) Inversie op S\mathcal S: f(x)=12πf^(ξ)eixξ ⁣dξ=12π(Ff^)(x)f(x) = \frac1{2\pi}\int\hat f(\xi)\eu^{\iu x\xi}\dd\xi = \frac1{2\pi}(\mathcal F\hat f)(-x), d.w.z. (F2f)(x)=2πf(x)(\mathcal F^2f)(x) = 2\pi f(-x). Tweemaal toepassen: F4f=2πF2f()=(2π)2f\mathcal F^4f = 2\pi\,\mathcal F^2f(-\cdot) = (2\pi)^2f.

(b) Als Ff=λf\mathcal Ff = \lambda f met f0f \neq 0: (2π)2f=F4f=λ4f(2\pi)^2f = \mathcal F^4f = \lambda^4f, dus λ4=(2π)2\lambda^4 = (2\pi)^2: λ{±2π,±i2π}\lambda \in \{\pm\sqrt{2\pi}, \pm\iu\sqrt{2\pi}\}. De Gauss g(x)=ex2/2g(x) = \eu^{-x^2/2} heeft g^=2πg\hat g = \sqrt{2\pi}\,g (Voorbeeld 14.2 bij a=12a = \frac12): eigenfunctie voor +2π+\sqrt{2\pi}.

(c) F2f=2πf()\mathcal F^2f = 2\pi f(-\cdot) is gelijk aan ±2πf\pm2\pi f volgens pariteit. Voor h(x)=xex2/2h(x) = x\eu^{-x^2/2}: differentiëren van g^(ξ)=2πeξ2/2\hat g(\xi) = \sqrt{2\pi}\eu^{-\xi^2/2} met de regel xf^=i ⁣d ⁣dξf^\widehat{xf} = \iu\frac{\dd}{\dd\xi}\hat f:

h^(ξ)=i ⁣d ⁣dξ(2πeξ2/2)=i2πξeξ2/2=i2πh(ξ):\hat h(\xi) = \iu\,\frac{\dd}{\dd\xi}\bigl(\sqrt{2\pi} \eu^{-\xi^2/2}\bigr) = -\iu\sqrt{2\pi}\,\xi\eu^{-\xi^2/2} = -\iu\sqrt{2\pi}\,h(\xi) :

een eigenfunctie voor i2π-\iu\sqrt{2\pi}. (De Hermite-functies zetten het patroon voort, cyclerend door de vier eigenwaarden — de discrete Fourier-klok.)

Oefening 14.12 ★★

(Autocorrelatie en lemma van Wiener) Voor fL2(R)f \in L^2(\R) definieer f~(x)=f(x)\tilde f(x) = \overline{f(-x)} en de autocorrelatie Af=ff~A_f = f * \tilde f. (a) Toon dat AfA_f een begrensde continue functie is met Af(0)=f22Af(x)A_f(0) = \norm f_2^2 \geq \abs{A_f(x)} voor alle xx (Oefening 12.6 en Cauchy–Schwarz). (b) Toon, eerst voor fL1L2f \in L^1\cap L^2, dat Af^=f^20\widehat{A_f} = \abs{\hat f\,}^2 \geq 0: de autocorrelatie heeft niet-negatieve transform — spectra van autocorrelaties zijn vermogensspectra. (c) Concludeer de identiteit Rf^(ξ)2eixξ ⁣dξ=2πAf(x)\int_\R\abs{\hat f(\xi)}^2\eu^{\iu x\xi}\,\dd\xi = 2\pi A_f(x) (inversie; rechtvaardig toepasbaarheid wanneer f^L2\hat f \in L^2 f^2L1\abs{\hat f}^2 \in L^1 heeft), en evalueer ze voor f=1[1/2,1/2]f = \mathbf 1_{\intcc{-1/2}{1/2}} in x=0x = 0: herwin R(sinuu)2 ⁣du=π\int_\R\bigl(\frac{\sin u}u\bigr)^2\dd u = \pi.

Oplossing

Oplossing van Oefening 14.12.

(a) f~L2\tilde f \in L^2 met f~2=f2\norm{\tilde f}_2 = \norm f_2; Oefening 12.6 (geconjugeerde exponenten p=q=2p = q = 2) maakt Af=ff~A_f = f * \tilde f begrensd en uniform continu, met

Af(x)=f(y)f(yx) ⁣dy,Af(0)=f22,Af(x)f2f(x)2=Af(0)A_f(x) = \int f(y)\,\overline{f(y - x)}\,\dd y, \qquad A_f(0) = \norm f_2^2, \qquad \abs{A_f(x)} \leq \norm f_2\,\norm{f(\cdot - x)}_2 = A_f(0)

via Cauchy–Schwarz.

(b) Voor fL1L2f \in L^1\cap L^2: f~L1\tilde f \in L^1 ook, en de convolutiestelling geeft Af^=f^f~^\widehat{A_f} = \hat f\, \widehat{\tilde f}; berekenend, f~^(ξ)=f(x)eiξx ⁣dx=f(u)eiξu ⁣du=f^(ξ)\widehat{\tilde f}(\xi) = \int\overline{f(-x)}\eu^{-\iu\xi x}\dd x = \overline{\int f(u)\eu^{-\iu\xi u}\dd u} = \overline{\hat f(\xi)}: Af^=f^20\widehat{A_f} = \abs{\hat f}^2 \geq 0.

(c) Wanneer f^2L1\abs{\hat f}^2 \in L^1, geldt inversie voor de continue AfA_f (haar transform is integreerbaar; Stelling 14.5):

Af(x)=12πf^(ξ)2eixξ ⁣dξ.A_f(x) = \frac1{2\pi}\int\abs{\hat f(\xi)}^2 \eu^{\iu x\xi}\,\dd\xi .

Voor f=1[1/2,1/2]f = \mathbf 1_{\intcc{-1/2}{1/2}}: f^(ξ)=2sin(ξ/2)ξ=sin(ξ/2)ξ/2\hat f(\xi) = \frac{2\sin(\xi/2)}\xi = \frac{\sin(\xi/2)}{\xi/2}, en in x=0x = 0:

1=f22=12πR(sin(ξ/2)ξ/2)2 ⁣dξ=12π2R(sinuu)2 ⁣du1 = \norm f_2^2 = \frac1{2\pi}\int_\R \Bigl(\frac{\sin(\xi/2)}{\xi/2}\Bigr)^2\dd\xi = \frac1{2\pi}\cdot2\int_\R\Bigl(\frac{\sin u}u\Bigr)^2\dd u

(ξ=2u\xi = 2u), d.w.z. R(sinuu)2 ⁣du=π\int_\R\bigl(\frac{\sin u}u\bigr)^2\dd u = \pi — Plancherels favoriete integraal, herwonnen via autocorrelatie.

14.6 Probleem: de warmtevergelijking op de lijn

Probleem 14.1

Weekendopgave — tu=xx2u\partial_tu = \partial^2_{xx}u, van begin tot eind opgelost

Warmte spreidt; de vergelijking tu=xx2u\partial_tu = \partial_{xx}^2u zegt dat haar dichtheid diffundeert met een snelheid gegeven door de lokale kromming van het temperatuurprofiel. We lossen het Cauchy-probleem op R\R op — gegeven ff, vind u(t,x)u(t, x) voor t>0t > 0 met u(0,)=fu(0, \cdot) = f — bewijzen de opmerkelijke eigenschappen van de oplossing, en zien waarom tijd niet omkeerbaar is. Doorheen is gt(x)=12πtex2/4tg_t(x) = \frac1{2\sqrt{\pi t}}\eu^{-x^2/4t} de warmtekern en u(t,)=gtfu(t, \cdot) = g_t * f.

Deel I — De kern afleiden. Werk eerst formeel: veronderstel u(t,)Su(t, \cdot) \in \mathcal S lost de vergelijking op, en laat u^(t,ξ)\hat u(t, \xi) de transform in xx zijn.

  1. Toon (formeel) tu^=ξ2u^\partial_t\hat u = -\xi^2\hat u, dus u^(t,ξ)=etξ2f^(ξ)\hat u(t,\xi) = \eu^{-t\xi^2}\hat f(\xi), en herken u(t)=gtfu(t) = g_t * f (Oefening 14.4). Dit motiveert de definitie van uu; alles wordt nu direct bewezen, voor fCb(R)f \in \mathcal C_b(\R) (begrensd continu) of fLpf \in L^p.

Deel II — Verificatie.

  1. Toon dat voor t>0t > 0, u(t,x)=gt(xy)f(y) ⁣dyu(t, x) = \int g_t(x-y)f(y)\dd y welgedefinieerd is voor fCbf \in \mathcal C_b, en dat uu C\mathcal C^\infty is in (t,x)(t, x) op (0,)×R\intoo0\infty\times\R (differentiëer onder de integraal; domineer Gauss-afgeleiden lokaal uniform in (t,x)(t,x)).
  2. Verifieer tgt=xx2gt\partial_tg_t = \partial^2_{xx}g_t door directe berekening, en concludeer tu=xx2u\partial_tu = \partial^2_{xx}u voor t>0t > 0.
  3. (Beginvoorwaarde) Toon dat voor fCbf \in \mathcal C_b, u(t,x)f(x)u(t, x) \to f(x) als t0+t \to 0^+, uniform op compacta (benaderende eenheid: splits yδ\abs y \leq \delta, y>δ\abs y > \delta); voor fLpf \in L^p (p<p < \infty), toon u(t)fp0\norm{u(t) - f}_p \to 0.
  4. (Instant gladmaking) Concludeer: zelfs voor louter begrensd continue ff is de oplossing C\mathcal C^\infty voor elke t>0t > 0 — warmte wist onmiddellijk ruwheid. Bereken u(t,)u(t, \cdot) expliciet voor f=1(0,)f = \mathbf 1_{\intoo0\infty} (een errorfunctie) en schets haar profiel voor drie waarden van tt.

Deel III — Structurele eigenschappen.

  1. (Positiviteit en vergelijking) Als f0f \geq 0 dan u>0u > 0 voor alle t>0t > 0 (strikt, tenzij f=0f = 0 a.e.); als f1f2f_1 \leq f_2 dan u1u2u_1 \leq u_2. Een koude plek warmt onmiddellijk: commentarieer.
  2. (Behoud) Voor fL1f \in L^1: u(t,x) ⁣dx=f\int u(t, x)\dd x = \int f voor alle tt (Tonelli) — totale warmte is behouden.
  3. (Dissipatie) Voor fL1L2f \in L^1\cap L^2, toon via Plancherel dat tu(t)2t \mapsto \norm{u(t)}_2 niet-stijgend is, strikt tenzij f=0f = 0, en bereken haar limiet als tt \to \infty. Toon bovendien u(t)f12πt0\norm{u(t)}_\infty \leq \frac{\norm f_1}{2\sqrt{\pi t}} \to 0: warmte spreidt en sterft.
  4. (Uniciteit, L2L^2-klasse) Zij uu een oplossing met u(t)L2u(t) \in L^2 voor alle tt, uC1((0,),L2)u \in \mathcal C^1(\intoo0\infty, L^2) in de natuurlijke zin en u(t)fu(t) \to f in L2L^2 als t0t\to0; toelatend dat de transform haar omzet in tu^=ξ2u^\partial_t\hat u = -\xi^2\hat u puntsgewijs a.e. in ξ\xi voor a.e. tt (gerechtvaardigd door testen tegen Cc\mathcal C_c^\infty in ξ\xi — schets dit), toon u^(t,ξ)=etξ2f^(ξ)\hat u(t,\xi) = \eu^{-t\xi^2}\hat f(\xi), dus uniciteit in deze klasse.

Deel IV — De pijl van de tijd.

  1. Toon dat het achterwaartse probleem slecht gesteld is: opdat de oplossing in tijd s-s (s>0s > 0) bestaat met data ff in tijd 00 — d.w.z. opdat f=gshf = g_s * h een oplossing hL2h \in L^2 heeft — is nodig dat esξ2f^(ξ)L2\eu^{s\xi^2}\hat f(\xi) \in L^2: een extreme vervalvoorwaarde op f^\hat f. Exposeer een expliciete gladde fL2f \in L^2 waarvoor geen achterwaartse oplossing op enig tijdsinterval bestaat: neem de functie met f^(ξ)=eξ\hat f(\xi) = \eu^{-\abs\xi} — identificeer ff (Oefening 14.1) en toon esξ2eξL2\eu^{s\xi^2}\eu^{-\abs\xi} \notin L^2 voor elke s>0s > 0.
  2. (Gladmaking vs. informatie) Leg in een korte alinea, met vragen 5, 9 en 10, uit waarom de warmtesemigroep (fgtf)t0(f \mapsto g_t * f)_{t\geq0} injectief maar niet surjectief is op L2L^2, en waarom dit de onomkeerbaarheid van diffusie uitdrukt.

Deel V — Shannons samplingstelling. Een functie fL2(R)f \in L^2(\R) is bandbegrensd tot Ω\Omega als f^=0\hat f = 0 a.e. buiten [Ω,Ω]\intcc{-\Omega}\Omega; schrijf PWΩPW_\Omega (Paley–Wiener-ruimte) voor deze.

  1. Toon dat elke fPWΩf \in PW_\Omega a.e. overeenstemt met de C\mathcal C^\infty-functie 12πΩΩf^(ξ)eixξ ⁣dξ\frac1{2\pi}\int_{-\Omega}^{\Omega}\hat f(\xi)\eu^{\iu x\xi}\,\dd\xi (rechtvaardig gladheid en de identificatie), met alle afgeleiden begrensd: bandbeperking is een extreme vorm van regulariteit. Vanaf nu duidt ff deze representant aan.
  2. Ontwikkel f^L2([Ω,Ω])\hat f \in L^2(\intcc{-\Omega}\Omega) in de Fourier-basis van dat interval en identificeer de coëfficiënten als samples van ff:

    f^(ξ)=πΩnZf(nπΩ)einπξ/Ωin L2([Ω,Ω]).\hat f(\xi) = \frac\pi\Omega\sum_{n\in\Z} f\Bigl(\frac{n\pi}\Omega\Bigr)\, \eu^{-\iu n\pi\xi/\Omega} \quad\text{in } L^2(\intcc{-\Omega}\Omega) .
  3. Concludeer de samplingstelling: voor fPWΩf \in PW_\Omega,

    f(x)=nZf(nπΩ)sinc(Ωxnπ),sinct=sintt,f(x) = \sum_{n\in\Z}f\Bigl(\frac{n\pi} \Omega\Bigr)\,\operatorname{sinc}(\Omega x - n\pi), \qquad \operatorname{sinc}t = \frac{\sin t}t,

    met convergentie in L2(R)L^2(\R) en uniform op R\R (injecteer de reeks van vraag 13 in de inversieformule en bereken de elementaire integraal): een bandbegrensd signaal is volledig bepaald door haar waarden op een rooster van stap π/Ω\pi/\Omega — de Nyquist-snelheid.

  4. Toon dat de functies xsinc(Ωxnπ)x \mapsto \operatorname{sinc}(\Omega x - n\pi), nZn \in \Z, een orthogonale familie in L2(R)L^2(\R) vormen met constante norm π/Ω\sqrt{\pi/\Omega}, en concludeer de energie-identiteit f22=πΩnf(nπ/Ω)2\norm f_2^2 = \frac\pi\Omega\sum_n\abs{f(n\pi/\Omega)}^2.
  5. (Aliasing) Exposeer een niet-nul gPW2Ωg \in PW_{2\Omega} die verdwijnt in elk samplepunt nπΩ\frac{n\pi}\Omega (beschouw g(x)=sin(Ωx)sinc(Ωx)g(x) = \sin(\Omega x)\operatorname{sinc}(\Omega x) en controleer haar band): samplen onder de Nyquist-snelheid verliest informatie — twee verschillende signalen kunnen alle samples delen: het stroboscopische wagenwiel-effect, gemathematiseerd.
  6. (Vrijheidsgraden) Met vragen 14–15, rechtvaardig de ingenieursregel: een signaal bandbegrensd tot Ω\Omega waarvan de energie in wezen door een tijdsvenster van lengte TT wordt gedragen, wordt beschreven door ongeveer ΩTπ\frac{\Omega T}\pi reële samples — maak “in wezen” precies via de energie-identiteit en de staart nπ/Ω>T/2\sum_{\abs{n\pi/\Omega} > T/2}.
  7. (Consistentiecontroles) Verifieer de samplingstelling met de hand op twee leden van PWΩPW_\Omega: (a) f=sinc(Ω)f = \operatorname{sinc}(\Omega\,\cdot), wier samples δn0\delta_{n0} zijn; (b) f(x)=cos(ωx)sinc(εx)f(x) = \cos(\omega x) \operatorname{sinc}(\varepsilon x)-type smalleband- signalen — preciezer, toon dat voor fPWΩf \in PW_{\Omega'} met Ω<Ω\Omega' < \Omega, de Ω\Omega-snelheidsreeks ff ook reconstrueert (oversampling is onschuldig), door PWΩPWΩPW_ {\Omega'} \subseteq PW_\Omega in te bedden.

Deel VI — Onzekerheid, nog tweemaal. Heisenbergs ongelijkheid (Oefening 14.8) begrens hoe geconcentreerd ff en f^\hat f gezamenlijk kunnen zijn; hier zijn haar alles-of-niets-zus en haar exacte verzadiging.

  1. Zij fL1f \in L^1 met suppf[A,A]\operatorname{supp}f \subseteq \intcc{-A}A. Toon dat f^\hat f de som is van een overal-convergente machtreeks:

    f^(ξ)=k0(iξ)kk!mk,mk=AAxkf(x) ⁣dx,mkAkf1\hat f(\xi) = \sum_{k\geq0}\frac{(-\iu\xi)^k}{k!} \,m_k, \qquad m_k = \int_{-A}^{A}x^kf(x)\,\dd x, \quad \abs{m_k} \leq A^k\norm f_1

    (ontwikkel eiξx\eu^{-\iu\xi x} en rechtvaardig de verwisseling door normale convergentie): de transform is reëel-analytisch, met oneindige straal van convergentie in elk punt.

  2. Concludeer de drager-dichotomie: een reëel-analytische functie die verdwijnt op een niet-leeg open interval verdwijnt identiek (de verzameling waar alle afgeleiden verdwijnen is open en gesloten — werk het Taylor-argument uit); concludeer dat geen niet-nul ff zowel ff als f^\hat f compact gedragen heeft, en dat PWΩPW_\Omega geen niet-nul compact-gedragen functie bevat — bandbegrensde signalen duren eeuwig, en tijdsbegrensde signalen lekken in alle frequenties.
  3. (Verzadiging van Heisenberg) Op de Gauss-familie f=eax2f = \eu^{-ax^2}, bereken beide concentratiefactoren en verifieer dat het genormaliseerde product (x2f2)(12πξ2f^2)/f24\bigl(\int x^2\abs f^2\bigr)\bigl(\frac1{2\pi}\int \xi^2\abs{\hat f}^2\bigr)\big/\norm f_2^4 gelijk is aan 14\frac14 voor elke aa — de gelijkheidsfamilie van Oefening 14.8 in levenden lijve; leg via een schalingsargument uit waarom het product constant moet zijn langs de familie.
  4. Leg de fysische lezing uit (positie/impulsdichtheden van een kwantumtoestand; \hbar in de normalisatie geeft σxσp2\sigma_x\sigma_p \geq \frac\hbar2), en verbind over het hoofdstuk: instant gladmaking (Deel II), onomkeerbaarheid (Deel IV), sampling (Deel V), Heisenberg en de drager-dichotomie zijn vijf uitdrukkingen van één wet — het gedrag van f^\hat f in oneindig wetgeeft wat ff overal mag doen.

Deel VII — De algebra van de kern, en één oplosbaar voorbeeld.

  1. (Semigroep) Bewijs de Chapman–Kolmogorov- identiteit gtgs=gt+sg_t * g_s = g_{t+s} voor t,s>0t, s > 0 (via de convolutiestelling en injectiviteit van de transform op L1L^1), en concludeer u(t+s)=gsu(t)u(t + s) = g_s * u(t): evolueren voor tijd t+st + s is evolueren voor tt, dan voor ss. Verscherp de dissipatie van vraag 8: schrijvend u(t)22=12πe2tξ2f^(ξ)2 ⁣dξ\norm{u(t)}_2^2 = \frac1{2\pi}\int \eu^{-2t\xi^2}\abs{\hat f(\xi)}^2\dd\xi, toon via Cauchy–Schwarz dat

    tlnu(t)2convex is op (0,+):t \longmapsto \ln\,\norm{u(t)}_2 \quad\text{convex is op } \intoo0{+\infty} :

    de L2L^2-energie daalt niet louter, ze daalt log-convex.

  2. (Waar de warmte naartoe gaat) Zij f0f \geq 0, fL1f \in L^1, met x2f(x) ⁣dx<\int x^2f(x)\dd x < \infty. Toon dat voor alle t>0t > 0

    Rxu(t,x) ⁣dx=Rxf(x) ⁣dx,Rx2u(t,x) ⁣dx=Rx2f(x) ⁣dx+2tRf:\int_\R x\,u(t, x)\,\dd x = \int_\R x f(x)\,\dd x, \qquad \int_\R x^2u(t, x)\,\dd x = \int_\R x^2f(x)\,\dd x + 2t\int_\R f :

    het zwaartepunt van de warmte beweegt nooit, en de variantie groeit lineair in de tijd — de diffusieve schaling x2tx \sim \sqrt{2t}, te herlezen wanneer Brownse beweging verschijnt in Hoofdstuk 22. (Bereken de eerste twee momenten van gtg_t en gebruik Tonelli op de convolutie.)

  3. (De Gauss, van begin tot eind opgelost) Voor f(x)=ex2f(x) = \eu^{-x^2}, vestig de gesloten vorm

    u(t,x)=11+4texp(x21+4t),u(t, x) = \frac1{\sqrt{1 + 4t}}\, \exp\Bigl(-\frac{x^2}{1 + 4t}\Bigr),

    en verifieer erop, met de hand: de vergelijking tu=xx2u\partial_tu = \partial^2_{xx}u; behoud u(t)=π\int u(t) = \sqrt\pi; de dissipatiewet u(t)2=(π/2)1/4(1+4t)1/4\norm{u(t)}_2 = (\pi/2)^{1/4}(1 + 4t)^{-1/4} (vergelijk haar t1/4t^{-1/4}-verval met het t1/2t^{-1/2} sup-norm-verval van vraag 8); en de exacte variantiegroei van vraag 24. Bij t=6t = 6: de piek is gedaald tot 15\frac15 van haar beginhoogte terwijl het profiel vijf keer breder is — dezelfde warmte, uitgespreid.

Oplossing

Oplossing van Probleem 14.1.

1. Transformeren van de vergelijking in xx (formeel): tu^(t,ξ)=xx2u^=(iξ)2u^=ξ2u^\partial_t\hat u(t,\xi) = \widehat{\partial^2_{xx}u} = (\iu\xi)^2\hat u = -\xi^2\hat u, een ODE in tt voor elke frequentie: u^(t,ξ)=etξ2f^(ξ)\hat u(t,\xi) = \eu^{-t\xi^2}\hat f(\xi). Omdat etξ2=g^t\eu^{-t\xi^2} = \hat g_t (Oefening 14.4), is het product de transform van gtfg_t * f.

2. u(t,x)fgt=f\abs{u(t,x)} \leq \norm f_\infty\int g_t = \norm f_\infty: welgedefinieerd. Op [t0,T]×[A,A][t_0, T]\times[-A, A]: elke gemengde afgeleide tmxngt(xy)\partial^m_t\partial^n_xg_t(x - y) is een polynoom in (xy)(x - y) en t1t^{-1} maal e(xy)2/4t\eu^{-(x-y)^2/4t}, begrensd voor y2A\abs y \geq 2A door C(1+y2)Ne(yA)2/4TC\,(1 + y^2)^N\eu^{-(\abs y - A)^2/4T}, een integreerbare dominator onafhankelijk van (t,x)(t, x) in het venster (en begrensd voor y2A\abs y \leq 2A): herhaalde differentiatie onder de integraal (Stelling 10.15) geldt: uC((0,)×R)u \in \mathcal C^\infty(\intoo0\infty\times\R).

3. Met gt(x)=12πtex2/4tg_t(x) = \frac1{2\sqrt{\pi t}}\eu^{-x^2/4t}:

tgt=gt(x24t212t)=xx2gt\partial_tg_t = g_t\Bigl(\frac{x^2}{4t^2} - \frac1{2t}\Bigr) = \partial^2_{xx}g_t

(differentiëer tweemaal in xx: xgt=x2tgt\partial_xg_t = -\frac x{2t}g_t, xx2gt=(x24t212t)gt\partial^2_{xx}g_t = \bigl(\frac{x^2}{4t^2} - \frac1{2t}\bigr)g_t). Via vraag 2 gaan de afgeleiden onder de integraal: tu=xx2u\partial_tu = \partial^2_{xx}u.

4. u(t,x)f(x)=gt(y)(f(xy)f(x)) ⁣dyu(t,x) - f(x) = \int g_t(y)\bigl(f(x - y) - f(x)\bigr)\dd y. Gegeven een compact KK en ε\varepsilon: uniforme continuïteit van ff op een omgeving van KK geeft δ\delta met f(xy)f(x)<ε\abs{f(x-y) - f(x)} < \varepsilon voor xKx \in K, yδ\abs y \leq \delta; de staart draagt 2fy>δgt(y) ⁣dy=2fP\leq 2\norm f_\infty\int_{\abs y > \delta}g_t(y)\dd y = 2\norm f_\infty\,\P-massa voorbij δ\delta, wat 2πδ/2tez2 ⁣dz0\frac2{\sqrt\pi}\int_{\delta/2\sqrt t}^\infty\eu^{-z^2}\dd z \to 0 is als t0t \to 0. Voor fLpf \in L^p: u(t)fpgt(y)τyffp ⁣dy\norm{u(t) - f}_p \leq \int g_t(y)\norm{\tau_yf - f}_p\dd y (Minkowski/Jensen zoals in Stelling 12.9), splits op dezelfde manier met Stelling 12.6(3).

5. Instant gladmaking is vraag 2 (u(t)u(t) is C\mathcal C^\infty voor t>0t > 0 zonder gladheid van ff te gebruiken). Voor f=1(0,)f = \mathbf 1_{\intoo0\infty}:

u(t,x)=0gt(xy) ⁣dy=1πx/2tez2 ⁣dz=12(1+erf(x2t)),erf(s)=2π0sez2 ⁣dz:u(t, x) = \int_0^\infty g_t(x - y)\dd y = \frac1{\sqrt\pi}\int_{-x/2\sqrt t}^{\infty}\eu^{-z^2}\dd z = \frac12\Bigl(1 + \operatorname{erf}\Bigl(\frac{x}{2\sqrt t}\Bigr)\Bigr), \qquad \operatorname{erf}(s) = \frac2{\sqrt\pi}\int_0^s\eu^{-z^2}\dd z :

een gladgemaakte stap waarvan de overgangszone als t\sqrt t verbreedt (profielen bij t1<t2<t3t_1 < t_2 < t_3: steeds vlakkere hellingen door (0,12)(0, \frac12)).

6. De integrand gt(xy)f(y)g_t(x-y)f(y) is 0\geq 0 en de kern is strikt positief: u(t,x)=0u(t,x) = 0 zou f=0f = 0 a.e. forceren. Monotonie in ff is monotonie van de integraal. Een plek waar f=0f = 0 op een interval heeft nog steeds u(t,)>0u(t, \cdot) > 0 daar voor elke t>0t > 0: warmte propageert met oneindige snelheid (elke positiviteit ergens wordt overal onmiddellijk gevoeld).

7. Tonelli (gt(xy)f(y)g_t(x-y)\abs{f(y)} is integreerbaar op R2\R^2): u(t,x) ⁣dx=f(y)(gt(xy) ⁣dx) ⁣dy=f\int u(t,x)\dd x = \int f(y)\bigl(\int g_t(x - y)\dd x\bigr)\dd y = \int f.

8. Plancherel: 2πu(t)22=e2tξ2f^(ξ)2 ⁣dξ2\pi\norm{u(t)}_2^2 = \int\eu^{-2t\xi^2}\abs{\hat f(\xi)}^2\dd\xi, niet-stijgend in tt (puntsgewijs), strikt tenzij f^=0\hat f = 0 a.e. (=f=0= f = 0), met limiet 00 als tt\to\infty via DCT. En u(t)gtf1=f12πt0\norm{u(t)}_\infty \leq \norm{g_t}_\infty\norm f_1 = \frac{\norm f_1}{2\sqrt{\pi t}} \to 0.

9. Voor φCc\varphi \in \mathcal C_c^\infty, tφ,u^(t)t \mapsto \langle\varphi, \hat u(t)\rangle is C1\mathcal C^1 met afgeleide φ,tu^=φ,xxu^=ξ2φ\langle\varphi, \partial_t\hat u\rangle = \langle\varphi, \widehat{\partial_{xx}u}\rangle = \langle\xi^2\varphi\dots\rangle — precies, xx2u^=ξ2u^\widehat{\partial^2_{xx}u} = -\xi^2\hat u transporteert de vergelijking. Dan voor a.e. ξ\xi heeft de absoluut continue functie tetξ2u^(t,ξ)t \mapsto \eu^{t\xi^2}\hat u(t,\xi) afgeleide etξ2(ξ2u^+tu^)=0\eu^{t\xi^2}(\xi^2\hat u + \partial_t\hat u) = 0 in de geïntegreerde zin: ze is constant, en t0t \to 0 latend (u^(t)f^\hat u(t) \to \hat f in L2L^2, a.e. langs een deelrij): u^(t,ξ)=etξ2f^(ξ)\hat u(t, \xi) = \eu^{-t\xi^2}\hat f(\xi) a.e. Twee oplossingen in de klasse hebben dezelfde transform: ze zijn gelijk.

10. f=gshf = g_s * h met hL2h \in L^2 forceert f^=esξ2h^\hat f = \eu^{-s\xi^2}\hat h, d.w.z. h^=esξ2f^L2\hat h = \eu^{s\xi^2}\hat f \in L^2. Neem f^(ξ)=eξ\hat f(\xi) = \eu^{-\abs\xi}: dan f(x)=1π11+x2f(x) = \frac1\pi\cdot\frac1{1 + x^2} (Oefening 14.1, inversie), een perfect gladde L2L^2-functie; maar e2sξ22ξ\eu^{2s\xi^2 - 2\abs\xi} \to \infty: esξ2f^L2\eu^{s\xi^2}\hat f \notin L^2 voor elke s>0s > 0. Het Cauchy-profiel is nooit het resultaat van eerdere diffusie.

11. De warmtesemigroep vermenigvuldigt transforms met etξ2\eu^{-t\xi^2}, die nergens verdwijnt: injectief — formeel wordt geen informatie vernietigd. Maar haar bereik bestaat uit functies wier transforms als etξ2\eu^{-t\xi^2} vervallen: een piepkleine, dichte-maar-echte deelruimte van L2L^2 (vraag 10 toont dat zelfs uitstekende functies erbuiten liggen). Inverteren zou de frequentie ξ\xi versterken met etξ2\eu^{t\xi^2}: onbegrensd, dus instabiel tegen elke verstoring. Diffusie is onomkeerbaar niet omdat de afbeelding vergeet, maar omdat haar inverse niet continu kan zijn — een pijl van de tijd gemaakt van functionaalanalyse.

12. f^L2([Ω,Ω])L1\hat f \in L^2(\intcc{-\Omega}\Omega) \subseteq L^1 (Cauchy–Schwarz op een begrensd interval), dus F(x)=12πΩΩf^(ξ)eixξ ⁣dξF(x) = \frac1{2\pi}\int_{-\Omega}^\Omega\hat f(\xi)\eu^{\iu x\xi}\dd\xi is overal gedefinieerd, en differentiatie onder de integraal (gedomineerd door Ωkf^L1\Omega^k\abs{\hat f} \in L^1 op de band) maakt ze C\mathcal C^\infty met F(k)Ωk2πf^L1\abs{F^{(k)}} \leq \frac{\Omega^k}{2\pi}\norm{\hat f}_{L^1} overal. En F=fF = f a.e.: beide zijden hebben dezelfde transform, en de transform is injectief op L2L^2 (Stelling 14.8 en haar L2L^2-uitbreiding).

13. De exponentialen ξeinπξ/Ω\xi \mapsto \eu^{-\iu n\pi\xi/\Omega}, nZn \in \Z, vormen een Hilbert-basis van L2([Ω,Ω])L^2(\intcc{-\Omega}\Omega) (Stelling 13.9, herschaald). De coëfficiënt van f^\hat f langs de nn-de is

12ΩΩΩf^(ξ)einπξ/Ω ⁣dξ=2π2Ω12πΩΩf^(ξ)ei(nπ/Ω)ξ ⁣dξ=πΩf(nπΩ),\frac1{2\Omega}\int_{-\Omega}^\Omega\hat f(\xi)\, \eu^{\iu n\pi\xi/\Omega}\dd\xi = \frac{2\pi}{2\Omega}\cdot \frac1{2\pi}\int_{-\Omega}^{\Omega}\hat f(\xi)\, \eu^{\iu(n\pi/\Omega)\xi}\dd\xi = \frac\pi\Omega\,f\Bigl(\frac{n\pi}\Omega\Bigr),

via de formule van vraag 12 in x=nπΩx = \frac{n\pi}\Omega: de gestelde ontwikkeling geldt in L2L^2 van de band.

14. Injecteer de ontwikkeling in de inversieformule van vraag 12; de verwisseling van som en integraal is de continuïteit van de L2L^2-paring tegen 12πeixξ1ξΩ\frac1{2\pi} \eu^{\iu x\xi}\mathbf 1_{\abs\xi\leq\Omega} (van L2L^2-norm 2Ω2π\frac{\sqrt{2\Omega}}{2\pi}, onafhankelijk van xx — vandaar uniformiteit):

f(x)=nf(nπΩ)12ΩΩΩei(xnπ/Ω)ξ ⁣dξ=nf(nπΩ)sinc(Ωxnπ),f(x) = \sum_nf\Bigl(\frac{n\pi}\Omega\Bigr)\cdot \frac1{2\Omega}\int_{-\Omega}^\Omega \eu^{\iu(x - n\pi/\Omega)\xi}\dd\xi = \sum_nf\Bigl(\frac{n\pi}\Omega\Bigr) \operatorname{sinc}(\Omega x - n\pi),

omdat 12ΩΩΩeiuξ ⁣dξ=sin(Ωu)Ωu\frac1{2\Omega}\int_{-\Omega}^\Omega\eu^{\iu u\xi}\dd\xi = \frac{\sin(\Omega u)}{\Omega u}.

15. De berekening van vraag 14 achterstevoren lezend, is de transform van sn=sinc(Ωnπ)s_n = \operatorname{sinc}(\Omega\cdot - n\pi) s^n=πΩeinπξ/Ω1[Ω,Ω]\hat s_n = \frac\pi\Omega\,\eu^{-\iu n\pi\xi/\Omega}\,\mathbf 1_{\intcc{-\Omega}\Omega}. Plancherel:

sn,sm=12π(πΩ)2ΩΩei(nm)πξ/Ω ⁣dξ=πΩδnm:\langle s_n, s_m\rangle = \frac1{2\pi} \Bigl(\frac\pi\Omega\Bigr)^2\int_{-\Omega}^\Omega \eu^{\iu(n-m)\pi\xi/\Omega}\dd\xi = \frac\pi\Omega\,\delta_{nm} :

een orthogonale familie van constante norm π/Ω\sqrt{\pi/\Omega}. Normen nemen in de ontwikkeling van vraag 14: f22=πΩnf(nπ/Ω)2\norm f_2^2 = \frac\pi\Omega\sum_n\abs{f(n\pi/\Omega)}^2.

16. g(x)=sin(Ωx)sinc(Ωx)=sin2(Ωx)Ωxg(x) = \sin(\Omega x)\operatorname{sinc} (\Omega x) = \frac{\sin^2(\Omega x)}{\Omega x} verdwijnt in elk roosterpunt nπΩ\frac{n\pi}\Omega (inclusief 00, via de limiet) en is niet identiek nul. Haar band: schrijf g=12i(eiΩxeiΩx)sinc(Ωx)g = \frac1{2\iu}\bigl(\eu^{\iu\Omega x} - \eu^{-\iu\Omega x}\bigr)\operatorname{sinc}(\Omega x); modulatie door e±iΩx\eu^{\pm\iu\Omega x} verschuift de transform met Ω\mp\Omega, dus g^\hat g is gedragen in [2Ω,2Ω]\intcc{-2\Omega}{2\Omega} (inderdaad in de unie van twee verschoven banden): gPW2Ωg \in PW_{2\Omega}, onzichtbaar voor Ω\Omega-snelheidssampling — aliasing in levenden lijve.

17. Via vraag 15 dragen de samples de energie democratisch: f2=πΩf(nπ/Ω)2\norm f^2 = \frac\pi\Omega\sum \abs{f(n\pi/\Omega)}^2. Als de energie van het signaal buiten het tijdsvenster [T/2,T/2]\intcc{-T/2}{T/2} ε2f2\leq \varepsilon^2\norm f^2 is, voldoen de samples buiten het venster (op randtermen na gecontroleerd door de uniforme grens van vraag 12) πΩnπ/Ω>T/2f(nπ/Ω)2f1x>T/22ε2f2\frac\pi\Omega\sum_{\abs{n\pi/\Omega} > T/2} \abs{f(n\pi/\Omega)}^2 \approx \norm{f\,\mathbf 1_{\abs x > T/2}}^2 \leq \varepsilon^2\norm f^2: de samplingreeks inkorten tot de ΩTπ\approx \frac{\Omega T}\pi in-venster indices reconstrueert ff tot relatieve fout ε\approx\varepsilon. Vandaar telt het tijd–bandbreedte- product ΩTπ\frac{\Omega T}{\pi} de effectieve reële vrijheidsgraden van het signaal — de regel achter elk audioformaat.

18. (a) sinc(Ωx)\operatorname{sinc}(\Omega x) heeft samples f(nπ/Ω)=sinc(nπ)=δn0f(n\pi/\Omega) = \operatorname{sinc}(n\pi) = \delta_{n0}: de reeks reduceert tot haar n=0n = 0-term, sinc(Ωx)\operatorname{sinc}(\Omega x) — de stelling reproduceert haar eigen kern. (b) Als f^\hat f gedragen is in [Ω,Ω][Ω,Ω]\intcc{-\Omega'}{\Omega'} \subseteq \intcc{-\Omega}\Omega, loopt elke stap van vragen 13–14 woordelijk met de grotere band Ω\Omega (de ontwikkeling van f^\hat f op het grotere interval is nog steeds legitiem): sneller samplen dan je eigen Nyquist-snelheid verandert niets in de reconstructie — oversampling is onschuldig, en in de praktijk nuttig (sneller-vervallende reconstructiekernen kunnen dan gebruikt worden).

19. Ontwikkel eiξx=k(iξx)kk!\eu^{-\iu\xi x} = \sum_k\frac{(-\iu\xi x)^k}{k!} binnen de integraal; op [A,A]\intcc{-A}A convergeert de reeks normaal (kξkAkk!fL1\sum_k\frac{\abs{\xi}^kA^k}{k!}\abs f \in L^1), dus term-voor-term-integratie is legitiem:

f^(ξ)=k0(iξ)kk!mk,mkAkf1.\hat f(\xi) = \sum_{k\geq0}\frac{(-\iu\xi)^k}{k!}m_k, \qquad \abs{m_k} \leq A^k\norm f_1 .

De grens maakt de reeks convergent voor elke complexe ξ\xi; rond elk punt ξ0\xi_0 geeft hergroeperen (absolute convergentie) een machtreeks in ξξ0\xi - \xi_0: f^\hat f is reëel-analytisch met oneindige straal overal.

20. Zij gg reëel-analytisch op R\R (Taylorreeks convergerend naar gg nabij elk punt) en Z={ξ:g(k)(ξ)=0 k}Z = \{\xi : g^{(k)}(\xi) = 0\ \forall k\}. ZZ is gesloten (doorsnede van gesloten verzamelingen); ze is open, omdat in ξ0Z\xi_0 \in Z de lokale Taylorontwikkeling van gg de nulreeks is, dus gg verdwijnt identiek nabij ξ0\xi_0, samen met alle afgeleiden. Als gg verdwijnt op een interval, ZZ \neq \varnothing; via samenhang van R\R, Z=RZ = \R: g0g \equiv 0. Nu als f0f \neq 0 compacte drager samen met f^\hat f had: vraag 19 maakt f^\hat f reëel-analytisch, verdwijnend buiten een compact, dus op intervallen: f^0\hat f \equiv 0, dus f=0f = 0 a.e. via injectiviteit — tegenspraak. Evenzo kan een niet-nul fPWΩf \in PW_\Omega niet compact gedragen zijn (verwissel de rollen van ff en f^\hat f via inversie): bandbegrensde signalen sterven nooit, tijdsbegrensde signalen bezetten onbegrensd spectrum.

21. Voor f=eax2f = \eu^{-ax^2}: f22=π2a\norm f_2^2 = \sqrt{\frac\pi{2a}} en x2f2=14aπ2a\int x^2\abs f^2 = \frac1{4a}\sqrt{\frac{\pi}{2a}} (Gauss-tweede moment); f^=πaeξ2/4a\hat f = \sqrt{\frac\pi a}\,\eu^{-\xi^2/4a} (Voorbeeld 14.2) en

12πξ2f^2 ⁣dξ=12ππaξ2eξ2/2a ⁣dξ=12aa2πa=2πa2.\frac1{2\pi}\int\xi^2\abs{\hat f}^2\dd\xi = \frac1{2\pi}\cdot\frac\pi a\int\xi^2 \eu^{-\xi^2/2a}\dd\xi = \frac1{2a}\cdot a\sqrt{2\pi a} = \frac{\sqrt{2\pi a}}2 .

Genormaliseerd product: 14aπ2a2πa2/π2a=14\frac1{4a}\sqrt{\frac\pi{2a}}\cdot \frac{\sqrt{2\pi a}}2\big/\frac{\pi}{2a} = \frac14, onafhankelijk van aa. Schaling verklaart de constantheid: ff vervangen door f(λ)f(\lambda\cdot) vermenigvuldigt x2f2/f2\int x^2\abs f^2/\norm f^2 met λ2\lambda^{-2} en 12πξ2f^2/f2\frac1{2\pi}\int\xi^2\abs{\hat f}^2/\norm f^2 met λ2\lambda^{2}: het product is dilatie-invariant, en de Gaussen vormen één dilatiebaan.

22. Met f2\abs f^2 de positiesdichtheid en 12πf^2\frac1{2\pi}\abs{\hat f}^2 de impulsdichtheid van een kwantumtoestand (fysische eenheden voegen \hbar in), leest Oefening 14.8 σxσp2\sigma_x\sigma_p \geq \frac\hbar2: geen toestand is scherp in beide observabelen. Over het hoofdstuk draagt één wet vijf pakken: warmte gladmaakt onmiddellijk omdat etξ2\eu^{-t\xi^2} hoge frequenties vernietigt (Deel II); de stroom kan niet achteruit lopen omdat ze herstellen onbegrensd is (Deel IV); een bandbegrensd signaal is star genoeg om op een aftelbaar rooster te leven (Deel V); geen functie slaat Heisenbergs vloer; en geen functie is compact gedragen aan beide zijden van de transform (vragen 19–20). Wat f^\hat f in oneindig doet, regelt wat ff overal mag doen.

23. Zowel gtg_t als gsg_s liggen in L1L^1 met gt^(ξ)=etξ2\widehat{g_t}(\xi) = \eu^{-t\xi^2} (berekening van vraag 1), dus de convolutiestelling geeft gtgs^=etξ2esξ2=e(t+s)ξ2=gt+s^\widehat{g_t * g_s} = \eu^{-t\xi^2}\eu^{-s\xi^2} = \eu^{-(t+s)\xi^2} = \widehat{g_{t+s}}; twee L1L^1-functies met dezelfde transform stemmen a.e. overeen (injectiviteit, via de inversiestelling — hier zijn beide zijden continu, dus ze stemmen overal overeen): gtgs=gt+sg_t * g_s = g_{t+s}. Bijgevolg u(t+s)=gt+sf=gs(gtf)=gsu(t)u(t + s) = g_{t+s} * f = g_s * (g_t * f) = g_s * u(t) (associativiteit van convolutie, Tonelli). Log-convexiteit: zij N(t)=u(t)22=12πe2tξ2f^2 ⁣dξN(t) = \norm{u(t)}_2^2 = \frac1{2\pi}\int \eu^{-2t\xi^2}\abs{\hat f}^2\dd\xi (Plancherel, vraag 8). Voor t=t1+t22t = \frac{t_1 + t_2}2, schrijf

e2tξ2f^2=(e2t1ξ2f^2)1/2(e2t2ξ2f^2)1/2,\eu^{-2t\xi^2}\abs{\hat f}^2 = \Bigl(\eu^{-2t_1\xi^2}\abs{\hat f}^2\Bigr)^{1/2} \Bigl(\eu^{-2t_2\xi^2}\abs{\hat f}^2\Bigr)^{1/2},

en Cauchy–Schwarz geeft N(t1+t22)N(t1)N(t2)N\bigl(\frac{t_1+t_2}2\bigr) \leq \sqrt{N(t_1)\,N(t_2)}: lnN\ln N is middenpuntconvex, en continu zijnde (gedomineerde convergentie in tt), convex; zo is lnu(t)2=12lnN(t)\ln\norm{u(t)}_2 = \frac12\ln N(t). Verval met een convexe logaritme: de warmtestroom kan energie niet in een uitbarsting verliezen en dan stagneren.

24. De momenten van de kern: gt=1\int g_t = 1 (vraag 7 met f=gsf = g_s, of direct de Gauss-integraal), xgt(x) ⁣dx=0\int x\,g_t(x)\dd x = 0 (oneven integrand), en, substituerend x=2tvx = 2\sqrt t\,v,

Rx2gt(x) ⁣dx=4tπRv2ev2 ⁣dv=2t.\int_\R x^2g_t(x)\,\dd x = \frac{4t}{\sqrt\pi}\int_\R v^2\eu^{-v^2}\dd v = 2t .

Substituerend x=z+yx = z + y in de convolutie en noteren (z+y)2gt(z)f(y) ⁣dz ⁣dy<\iint(\abs z + \abs y)^2g_t(z)f(y)\,\dd z\,\dd y < \infty (elk van zkgt\int\abs z^kg_t, ykf\int\abs y^kf is eindig voor k2k \leq 2, met y1+y22\abs y \leq \frac{1 + y^2}2), Fubini en Tonelli gelden voor de momentintegralen hieronder:

xu(t,x) ⁣dx=(z+y)gt(z)f(y) ⁣dz ⁣dy=0 ⁣ ⁣f+1 ⁣ ⁣yf(y) ⁣dy,\int x\,u(t,x)\dd x = \iint (z + y)\,g_t(z)f(y)\,\dd z\,\dd y = 0\cdot\!\int\! f + 1\cdot\!\int\! yf(y)\dd y,

wat de eerste claim is; en

(z+y)2gt(z)f(y) ⁣dz ⁣dy=2tf+20 ⁣ ⁣yf+y2f(y) ⁣dy,\iint (z+y)^2g_t(z)f(y)\,\dd z\,\dd y = 2t\int f + 2\cdot0\cdot\!\int\! yf + \int y^2f(y)\dd y ,

de tweede. Gemiddelden tellen op, varianties tellen op, en de kern draagt gemiddelde 00 en variantie 2t2t bij: na tijd tt heeft de warmte zich over een breedte van orde 2t\sqrt{2t} verspreid — afstand groeit als de vierkantswortel van de tijd, de handtekening van diffusie (en van de Brownse paden van Hoofdstuk 22).

25. Transformzijde: f^(ξ)=πeξ2/4\hat f(\xi) = \sqrt\pi\,\eu^{-\xi^2/4}, dus u^(t,ξ)=πe(t+14)ξ2\hat u(t,\xi) = \sqrt\pi\,\eu^{-(t + \frac14)\xi^2}, wat de transform is van (1+4t)1/2exp(x2/(1+4t))(1 + 4t)^{-1/2}\exp\bigl(-x^2/(1+4t)\bigr) (het Gauss-woordenboek eax2π/aeξ2/4a\eu^{-ax^2} \mapsto \sqrt{\pi/a}\,\eu^{-\xi^2/4a} met a=11+4ta = \frac1{1+4t}): de gesloten vorm. Directe controle, met σ=1+4t\sigma = 1 + 4t:

tu=σ1/2ex2/σ(2σ+4x2σ2)=xx2u,\partial_tu = \sigma^{-1/2}\eu^{-x^2/\sigma} \Bigl(-\frac2\sigma + \frac{4x^2}{\sigma^2}\Bigr) = \partial^2_{xx}u ,

beide zijden berekend uit xu=2xσu\partial_xu = -\frac{2x}\sigma\,u. Behoud: u(t)=σ1/2πσ=π\int u(t) = \sigma^{-1/2}\sqrt{\pi\sigma} = \sqrt\pi voor alle tt. Dissipatie:

u(t)22=1σe2x2/σ ⁣dx=1σπσ2=π2(1+4t)1/2,\norm{u(t)}_2^2 = \frac1\sigma\int\eu^{-2x^2/\sigma}\dd x = \frac1\sigma\sqrt{\frac{\pi\sigma}2} = \sqrt{\frac\pi2}\,(1+4t)^{-1/2},

dus u(t)2=(π/2)1/4(1+4t)1/4\norm{u(t)}_2 = (\pi/2)^{1/4}(1+4t)^{-1/4}, niet-stijgend met convexe logaritme (vraag 23); de L2L^2-norm vervalt als t1/4t^{-1/4}, precies de helft van de t1/2t^{-1/2}-exponent van u(t)\norm{u(t)}_\infty — consistent met u22uu1\norm u_2^2 \leq \norm u_\infty\norm u_1 en behoud van u1\norm u_1. Variantie: x2u(t)=1σσ3/2π2=π2(1+4t)=x2f+2tπ\int x^2u(t) = \frac1\sigma\cdot\frac{\sigma^{3/2}\sqrt\pi}2 = \frac{\sqrt\pi}2\,(1 + 4t) = \int x^2f + 2t\sqrt\pi, zoals vraag 24 voorspelt (x2f=π2\int x^2f = \frac{\sqrt\pi}2, f=π\int f = \sqrt\pi). Bij t=6t = 6: σ=25\sigma = 25, piekhoogte u(6,0)=15u(6, 0) = \frac15 versus u(0,0)=1u(0,0) = 1, breedteschaal σ=5\sqrt\sigma = 5 maal de beginwaarde, en u=π1.7725\int u = \sqrt\pi \approx 1.7725 doorheen: de vlek is vijf keer lager, vijf keer breder, en er ontbreekt geen calorie.