Uge 8: Riemann-Integralet i n-D#

Nøglebegreber#

  • Riemann-integralet for skalar-funktioner af n variable

  • Riemann-integralet for vektor-funktioner

  • Transformationssætningen: Koordinatskifte i \(\mathbb{R}^n\)

  • Jacobi-determinanten

  • Typiske koordinater:

    • I \(\mathbb{R}^2\): Kartesiske og polære koordinater

    • I \(\mathbb{R}^3\): Kartesiske, sfæriske/kugle, cylinder/semi-polære koordinater

Forberedelse og pensum#

  • Læsepensum: Resten af kapitel 6

  • Python demo


Opgaver – Store Dag#

1: Planintegraler over rektangler. Håndregning#

Spørgsmål a#

Betragt området \(B=\left\lbrace (x,y) \bigm| 0\leq x\leq 2 \wedge -1\leq y\leq 0\right\rbrace\) i \(\mathbb{R}^2\). Udregn planintegralet

\[\begin{equation*} \int_B (x^2y^2+x) \mathrm{d}\pmb{x} \end{equation*}\]

via formlen for dobbelt integraler over (akseparallelle) rektangler.

Spørgsmål b#

Vi vil udregne samme planintegral igen, men på en umiddelbart mere kompliceret måde, nemlig via Transformationssætningen for integraler over \(\mathbb{R}^2\).

Spørgsmål c#

Udregn planintegralet

\[\begin{equation*} \int_B \frac{y}{1+xy} \;\mathrm{d}\pmb{x}, \quad\text{hvor}\quad B=\left\lbrace (x,y) \mid 0\leq x\leq 1 \, \wedge \, 0\leq y\leq 1\right\rbrace \end{equation*}\]

2: Polære koordinater. Håndregning#

En funktion \(f:\mathbb{R}^2 \to \mathbb{R}\) er givet ved

\[\begin{equation*} f(x,y)=x^2-y^2 \end{equation*}\]

For et givet punkt \(\pmb{x}=(x,y)\) i planen betegner \(r = \Vert \pmb{x} \Vert\) punktets afstand til origo \((0,0)\). Tilsvarende betegner \(\theta\) vinkel mellem \(x\)-aksen og punktets stedvektor, regnet med fortegn i omdrejningsretningen mod uret. En punktmængde \(B\) er i polære koordinater beskrevet som de punkter for hvilket

\[\begin{equation*} 0\leq r \leq a \, \text{ og } \, -\frac{\pi}{4} \leq \theta \leq \frac{\pi}{2}, \end{equation*}\]

hvor \(a\) er et vilkårligt positivt reelt tal.

Spørgsmål a#

Lav en skitse af \(B\), og bestem arealet af \(B\) både ved integration og elementær geometrisk betragtning.

Spørgsmål b#

Bestem planintegralet \(\int_B f(x,y) \;\mathrm{d}\pmb{x}\).

3: Volumen af et parallellotop#

Et parallellotop \(P\) i \(\mathbb{R}^n\) “udspændt” af vektorerne \(\pmb{a}_1, \pmb{a}_2, \dots, \pmb{a}_n\) defineres ved:

\[\begin{equation*} P = \left\{ \pmb{y} \in \mathbb{R}^n \mid \, \pmb{y} = A\pmb{x}, \quad \text{hvor } x_i \in [0,1] \text{ for $i=1,2,\dots, n$} \right\} \end{equation*}\]

hvor \(A = [\pmb{a}_1 | \pmb{a}_2 | \cdots | \pmb{a}_n]\) er \(n \times n\) matricen hvis \(i\)’te søjle er \(\pmb{a}_i\). Punktmængden kan kort skrives \(P=A([0,1]^n)\).

Det kan vises med værktøjer kun fra Matematik 1a (specielt karakteriseringen af determinanten) at det \(n\)-dimensionale volumen af \(P\) er:

\[\begin{equation*} \mathrm{vol}_n(P) = |\mathrm{det}(A)| \end{equation*}\]

(For den interesserede studerende: et sådant bevis findes her https://textbooks.math.gatech.edu/ila/determinants-volumes.html)

I \(\mathbb{R}^2\) er et parallellotop det velkendte pallellogram, og \(\mathrm{vol}_n(P)\) er arealet af \(P\), mens man i \(\mathbb{R}^3\) genfinder parallelepipedummet og det almindelige volumen.

Spørgsmål a#

Vis \(\mathrm{vol}_n(P) = |\mathrm{det}(A)|\) ved hjælp af Transformationssætningen for integraler over \(\mathbb{R}^n\).

I resten af opgaven ønsker vi at undersøge udsagnet \(\mathrm{vol}_n(P) = |\mathrm{det}(A)|\) uden brug af integrationsteknikker.

Spørgsmål b#

Lad \(n=2\). Vælg to lineært uafhængige vektorer \(\pmb{a}_1, \pmb{a}_2\) i \(\mathbb{R}^2\). Du kan med fordel vælge \(\pmb{a}_1 \in \mathrm{span}(\pmb{e}_1)\). Udregn (ved simple geometriske overvejelser) arealet af parallellogrammet “udspændt” af de to vektorer. Udregn også \(|\mathrm{det}(A)|\) og sammenlign størrelserne.

Spørgsmål c#

Lad \(n=2\) og lad nu \(\pmb{a}_1, \pmb{a}_2\) være vilkårlige, men lineært uafhængige vektorer i \(\mathbb{R}^2\). Kan du bevise formlen \(\mathrm{areal}(P) = |\mathrm{det}(A)|\), hvor \(P\) er parallellogrammet “udspændt” af de to vektorer. Du kan antage (hvorfor?) at \(\pmb{a}_1 \in \mathrm{span}(\pmb{e}_1)\) hvis det gør dit argument lettere.

Spørgsmål d#

Lad \(n=3\). Vælg tre lineært uafhængige vektorer \(\pmb{a}_1, \pmb{a}_2, \pmb{a}_3\) i \(\mathbb{R}^3\). Du kan med fordel vælge \(\pmb{a}_1, \pmb{a}_2 \in \mathrm{span}(\pmb{e}_1, \pmb{e}_2)\). Udregn (ved simple geometriske overvejelser) volumen af parallelepipedummet “udspændt” af de tre vektorer. Udregn også \(|\mathrm{det}(A)|\) og sammenlign størrelserne.

Spørgsmål e (frivillig/ekstra)#

Lad \(n=3\) og lad nu \(\pmb{a}_1, \pmb{a}_2, \pmb{a}_3\) være vilkårlige, men lineært uafhængige vektorer i \(\mathbb{R}^3\). Kan du bevise formlen \(\mathrm{areal}(P) = |\mathrm{det}(A)|\), hvor \(P\) er parallelepipedummet “udspændt” af de tre vektorer. Du kan antage (hvorfor?) at \(\pmb{a}_1, \pmb{a}_2 \in \mathrm{span}(\pmb{e}_1, \pmb{e}_2)\), hvis det gør dit argument lettere.

4: Planintegral med parametrisering I#

I \((x,y)\)-planen er der givet punktet \(P_0=(1,2)\) og punktmængden

\[\begin{equation*} C=\left\lbrace (x,y)\Big\vert \frac 32\leq y \leq \frac 52 \wedge 0\leq x\leq \frac 12 y^2\right\rbrace \end{equation*}\]

Spørgsmål a#

Lav en foreløbig skitse af \(C\) og angiv en parameterfremstilling \(\pmb{r}(u,v)\) for \(C\) med passende intervaller for \(u\) og \(v\), dvs angiv \(\Gamma\)\(\pmb{r}(\Gamma)=C\). Argumentér for at den valgte parametrisering er injektiv (hvis den valgte parametrisering ikke er injektiv, skal du finde en ny).

Spørgsmål b#

Bestem de to parameterværdier \(u_0\) og \(v_0\) således at \(\pmb{r}(u_0,v_0)=P_0\). Lav en illustration af \(C\) (fx med Sympy) hvor du udfra \(P_0\) afsætter tangentvektorerne \(\pmb{r}'_u(u_0,v_0)\) og \(\pmb{r}'_v(u_0,v_0)\). Bestem arealet af det parallelogram som udspændes af tangentvektorerne, jf. Opgave 3: Volumen af et parallellotop.

Spørgsmål c#

Bestem den til \(\pmb{r}(u,v)\) hørende Jacobi-determinant og argumenter for at de to søjle-vektorer i Jacobi-matricen er lineært uafhængige for alle \((u,v) \in \Gamma\). Udregn Jacobi-determinanten i punktet \((u_0,v_0)\).

Spørgsmål d#

Udregn planintegralet:

\[\begin{equation*} \int_C \frac{1}{y^2+x} \mathrm{d}\pmb{x} \end{equation*}\]

via Transformationssætningen for integraler over \(\mathbb{R}^2\). Du skal argumentere for at Transformationssætningen kan bruges.

5: Planintegral med parametrisering II#

Vi skal bestemme planintegralet

\[\begin{equation*} \int_B 2xy\,\mathrm{d} \pmb{x} \quad\text{hvor}\quad B=\left\lbrace (x,y) \mid 0\leq x \, \wedge \, 0\leq y, x+y\leq 1\right\rbrace \end{equation*}\]

Følg de nedenstående trin.

Spørgsmål a#

Skitsér først området \(B\). Bestem derefter en parameterfremstilling for \(B\).

Spørgsmål b#

Bestem den Jacobi-determinant, som svarer til denne parametrisering. Er Jacobi-determinanten forskellig fra nul på det indre af parameterområdet (som kræves af Transformationssætningen)?

Spørgsmål c#

Bestem nu det ønskede integral.

6: Partiel integration og substitution i to variable#

Spørgsmål a#

Bestem \(\displaystyle{\int_0^{\frac{\pi}{2}}\left(\int_0^{\frac{\pi}{2}} u\cos(u+v)\mathrm{d}u\right)\mathrm{d}v.}\)

Spørgsmål b#

Bestem \(\displaystyle{\int_0^1\left(\int_0^1 \frac{v}{(uv+1)^2}\mathrm{d}u\right)\mathrm{d}v.}\)

7: Et tripelintegral#

Bestem trippelintegralet

\[\begin{equation*} \displaystyle{\int_1^2\int_1^2\int_1^2 \frac{xy}{z} \mathrm dx\mathrm dy\mathrm dz.}\ \end{equation*}\]

Opgaver – Lille Dag#

Bemærk: På engelsk kaldes \(|\det{\pmb{J}_{\pmb{r}}(\pmb{u})}|\) ofte for “the Jacobian”. På dansk vil man nok kalde den “Jacobianten” eller bruge den fulde beskrivelse: “absolut-værdien af determinanten af Jacobi-matricen”. På DTU kaldes den undertiden for “Jacobi-funktionen”, men dette er ikke en betegnelse der bruger uden for DTU - og man tales således ikke om “the Jacobi function” i engelsk sproget litteratur.

1: Parametriseret rumligt område. Håndregning.#

Et område \(B\) i \((x,y,z)\)-rummet er givet ved parameterfremstillingen

\[\begin{equation*} \pmb{r}(u,v,w)=\big(\frac{1}{2}u^2-v^2,-uv,w\big),\quad u\in \left[ 0,2\right],v\in \left[ 0,2\right],w\in \left[ 0,2\right] \end{equation*}\]

Spørgsmål a#

I \(B\) er der givet punktet

\[\begin{equation*} \pmb{x}_0=\pmb{r}(1,1,1) \end{equation*}\]

Find \(\pmb{x}_0\). Afsat ud fra \(\pmb{x}_0\) udspænder tangentvektorerne \(\pmb{r}_u'(1,1,1),\pmb{r}_v'(1,1,1)\) og \(\pmb{r}_w'(1,1,1)\) et parallelepipedum \(P\), jf. Opgave 3: Volumen af et parallellotop. Bestem volumen af dette parallelepipedum. Illustrér evt. med Sympy.

Spørgsmål b#

Bestem absolut-værdien af Jacobi-determinanten der hører til \(\pmb{r}\). Evaluer den i \(\pmb{x}_0\).

Spørgsmål c#

Bestem voluminet af \(B\).

2: Massefordelinger i \((x,y)\)-planen#

Betragt punktmængderne i \(\mathbb{R}^2\) givet ved:

\[\begin{equation*} B=\left\lbrace (x,y)\vert 1\leq x\leq 2 \, \wedge \, 0\leq y\leq x^3\right\rbrace \end{equation*}\]

og (igen)

\[\begin{equation*} C=\left\lbrace (x,y)\Big\vert \frac 32\leq y \leq \frac 52 \wedge 0\leq x\leq \frac 12 y^2\right\rbrace. \end{equation*}\]

Vi opfatter \(f(x,y)\) som massetætheden/massefylden (kg/m\(^2\)) i punktet \((x,y)\).

Spørgsmål a#

Antag massetætheden er konstant \(f(x,y)=1\) for \((x,y)\in B\). Bestem massen og massemidtpunktet af \(B\).

Spørgsmål b#

Antag massetætheden er \(f(x,y)=x^2\) for \((x,y)\in B\). Bestem massen og massemidtpunktet af \(B\).

Spørgsmål c#

Antag massetætheden er konstant \(f(x,y)=1\) for \((x,y)\in C\). Bestem massen og massemidtpunktet af \(C\).

Spørgsmål d#

Antag massetætheden er \(f(x,y)=x^2\) for \((x,y)\in C\). Bestem massen og massemidtpunktet af \(C\).

3: Kugleformede områder i rummet#

Betragt det rumlige område \(\pmb{r}(\Gamma)\) givet ved

\[\begin{equation*} \pmb{r}(u,v,w)=\big(u\sin(v)\cos(w),u\sin(v)\sin(w),u\cos(v)\big), \quad (u,v,w) \in \Gamma, \end{equation*}\]

hvor \(\Gamma = [a,b] \times [c,d] \times [e,f] \subset [0, \infty[ \times [0,\pi] \times [0,2\pi]\). Vi betragter altså følgende parameterværdier: \(u\in [a,b],v\in [c,d],w\in [e,f]\).

Spørgsmål a#

Hvilken betydning har parametrene?

Spørgsmål b#

Lad \(A\) være det område der er bestemt ved valget:

\[\begin{equation*} a=1,b=3,c=\frac{\pi}{4},d=\frac{\pi}{3},e=0,f=\frac{3\pi}{4} \end{equation*}\]

og \(B\) ved valget

\[\begin{equation*} a=2,b=4,c=\frac{\pi}{4},d=\frac{\pi}{2},e=-\frac{\pi}{4},f=\frac{\pi}{4} \end{equation*}\]

Beskriv i ord hvert af områderne \(A\), \(B\) og \(A\cap B\), og bestem deres volumen.

Spørgsmål c#

Find integralerne

\[\begin{equation*} \int_A x_1 \, \mathrm{d}\pmb{x}, \quad \int_Bx_1 \, \mathrm{d}\pmb{x} \quad \text{og} \quad \int_{A\cap B}x_1 \, \mathrm{d}\pmb{x} \end{equation*}\]

4: Et uegentligt integral i planen#

Lad \(B\) være enhedskvadratet \([0,1]^2\). Vi vil undersøge det uegentlige planintegral

\[\begin{equation*} I := \int_B \frac{1}{x_2-x_1-1} \mathrm{d}\pmb{x} \end{equation*}\]

Integranden \(f(x_1,x_2)=\frac{1}{x_2-x_1-1}\) er ikke Riemann-integrabel over \(B\), da \(f\) ikke er defineret i punktet \((x_1,x_2)=(0,1)\). Vi ønsker at finde ud af om vi alligevel kan tillægge integralet en værdi via grænseovergang.

Spørgsmål a#

Find de punkter i planen \((x,y)\) hvor \(f(x_1,x_2)\) ikke er defineret. Find værdimængden af \(f\) som funktion på \(B \setminus \{(0,1)\}\).

Spørgsmål b#

Lad \(B_a = [a,1] \times [0,1]\) for et fast \(a \in [0,1]\). Lav en skitse af \(B_a\) og lav en parametrisering af \(B_a\). Udregn Jacobi-determinanten af parametriseringen.

Spørgsmål c#

Udregn Riemann-integralet

\[\begin{equation*} I_a := \int_{B_a} \frac{1}{x_2-x_1-1} \mathrm{d}\pmb{x} \end{equation*}\]

for hvert \(a \in ]0,1]\).

Spørgsmål d#

Udregn grænsen af \(I_a\) for \(a \to 0\).

Spørgsmål e#

Lad \(B_b = [0,1] \times [0,b]\). Definer \(I_b := \int_{B_b} \frac{1}{x_2-x_1-1} \mathrm{d}\pmb{x} \). Find \(\lim_{b \to 1} I_b\) og sammenlign med overstående.