6 svar
127 visningar
naytte 7419 – Moderator
Postad: 9 aug 20:36 Redigerad: 9 aug 20:42

Är detta en rimlig intuition bakom Jacobianen?

Halloj!

Jag håller på att befatta mig med variabelbyten i dubbel- och trippelintegraler och håller på att försöka förstå intuitivt varför man behöver jacobianen när man genomför variabelbyten. Jag undrar om det jag presenterar nedan är en rimlig motivering till att Jacobianen behövs. Betrakta följande delmängd av 2\mathbb{R}^2:

S:={(x,y)2:x2+y2<4}\displaystyle S:=\{ (x,y)\in\mathbb{R}^2:x^2+y^2<4 \}

Denna delmängd definerar en öppen disk med radie 22 i xyxy-planet. Det följer att vi har:

SdA=4π\displaystyle \iint_SdA = 4\pi

Betrakta nu följande transformation T:S2T: S \to \mathbb{R}^2:

u,v=Tx,y:=x2,y2\displaystyle \left(u,v\right)=T\left(x,y\right):=\left(\frac{x}{2},\frac{y}{2}\right)

Bilden av hela området SS under TT är då:

T(S)=:{T(x):xS}={(u,v)2:u2+v2<1}\displaystyle T(S)=:\{ T(x):x\in S \}=\{ (u,v)\in\mathbb{R}^2:u^2+v^2< 1 \}

Om vi alltså naivt skulle försöka beräkna vår integral med de nya koordinaterna utan förbehåll skulle vi kanske beräkna:

T(S)dudv=π\displaystyle \iint_{T(S)}dudv=\pi

Nu har vi alltså fått en area som är fyra gånger för liten! Faktum är att om vi lägger till denna faktor har vi exakt:

SdA=T(S)4dA=4π\displaystyle \iint_SdA=\iint_{T(S)}4dA=4\pi

Detta motiverar att vi behöver någon typ av "skalningsfaktor" när vi genomför koordinatbyten och byter "värld" (vi går från xyxy-världen till uvuv-världen). När vi transformerade området så "tryckte vi ihop" punkterna ur SS till ett mindre område (ett område med mindre Lebesguemått) T(S)T(S). Är detta en rimlig (extremt heuristisk) motivering till att man behöver Jacobianen, fastän vi inte explicit har talat om denna ännu? Vi behöver någon typ av konverteringsfaktor mellan xyxy- och uvuv-planen.

LuMa07 495
Postad: 9 aug 21:08 Redigerad: 9 aug 21:16

Det stämmer bra.

I kursen i linjär algebra har du säkert stött på linjära avbildningar/transformationer i planet. När man har en figur i planet med en viss area och denna figur utsätts för en linjär transformation TT, så kommer figurens area förändras i enlighet med determinanten av transformationens avbildningsmatris AA, d.v.s.

area(T(S))=|detA|area(S)\text{area}(T(S)) = |\det A| \, \text{area}(S)

Determinanten funkar alltså som en skalningsfaktor för areor. (Tecknet av determinanten avgör om avbildningen bevarar orientering eller vänder på den, så tecknet är oväsentligt för figurens area.)

Om man nu kollar på koordinatvariabler för den linjära transformationen TT, d.v.s.

uv=Axy\begin{pmatrix}u\\v\end{pmatrix} = A \begin{pmatrix}x\\y\end{pmatrix} , så innebär det att u(x,y)=A11x+A12yv(x,y)=A21x+A22y\left\{\begin{array}{l}u(x,y)=A_{11}x + A_{12}y\\v(x,y)=A_{21}x + A_{22}y\end{array}\right..

När man beräknar partiella derivator och när man bildar Jacobianen, så får man

ux'=A11u^\prime_x = A_{11}uy'=A12u^\prime_y = A_{12}vx'=A21v^\prime_x = A_{21}vy'=A22v^\prime_y = A_{22}

och därmed

d(u,v)d(x,y)=detA11A12A21A22=detA \dfrac{d(u,v)}{d(x,y)} = \det \begin{pmatrix} A_{11} & A_{12} \\ A_{21} & A_{22} \end{pmatrix} = \det A.

Detta visar att jacobianen för denna linjära transformation är lika med avbildningsmatrisens determinant, vilket enligt kunskaper från linjär algebra funkar som skalningsfaktor för transformationens verkan på arean.

 

Hos linjära transformationer har man samma skalningsfaktor överallt, så jacobianen är konstant. Hos icke-linjära transformationer har man olika skalningsfaktorer på olika ställen, men jacobianen varierar och därmed fångar skalningsfaktorn lokalt.

naytte 7419 – Moderator
Postad: 9 aug 21:31 Redigerad: 9 aug 21:31

Tack för ditt svar! Så motiveringen bakom det lite mer komplexa objektet Jacobianen är att området kan skalas annorlunda beroende på var någonstans i området vi är. Om vi går från t.ex. rektangulära koordinater xx och yy till polära koordinater rr och θ\theta, så kommer areaelementen skalas annorlunda beroende på var i området vi befinner oss, dvs. skalningsfaktorn är en funktion i rr och θ\theta, som vi skriver som J(r,θ)J(r,\theta)?

D4NIEL 3345
Postad: 9 aug 21:54 Redigerad: 9 aug 22:10

Ja, det stämmer, Jakobianen är en lokal beskrivning av transformationen och behöver inte vara konstant. För polära koordinater är den

JTr,θ=(x,y)(r,θ)=rJ_T\left(r,\theta\right)=\left|\frac{\partial (x,y)}{\partial(r,\theta)}\right|=r

Tönk på att få transformationen åt rätt håll, formeln är 

kfydy=kfTxJTxdx\displaystyle \int_{\mathbb{R}^k}f\left(\mathbf{y}\right)\,\mathrm{d}\mathbf{y}=\int_{\mathbb{R}^k}f\left(T\left(\mathbf{x}\right)\right)\left|J_T\left(\mathbf{x}\right)\right|\,\mathrm{d}\mathbf{x}

För din exempeltransformation får vi då

x=2u,  y=2vx=2u, \quad y=2v

Och därmed är skalfaktorn åt det hållet JTu,v=x,yu,v=4\left|J_T\left(u,v\right)\right|=\left|\frac{\partial\left(x,y\right)}{\partial\left(u,v\right)}\right|=4

Notationsmässigt anger man alltså variablerna som vi ska integrera eller derivera med avseende på, vilket kan kännas lite bakvänt i början.

Det är alltså lätt hänt att man istället slarvar och beräknar (u,v)(x,y)\frac{\partial(u,v)}{\partial(x,y)} och det blir inte lättare av att det finns böcker som vänder på notationen.

naytte 7419 – Moderator
Postad: 9 aug 23:11 Redigerad: 9 aug 23:12

Hur kommer det sig att du i båda fallen skriver området som k\mathbb{R}^k? Förändras inte området när vi tillämpar transformationen? T.ex. om ursprungsområdet är SS så går vi till området T(S)T(S) under en transformation TT.

D4NIEL 3345
Postad: 10 aug 00:31 Redigerad: 10 aug 00:34

Ja, jag menade inget särskilt med det mer än att vi går från k\mathbb{R}^k till k\mathbb{R}^k. Som du påpekar tänker vi oss att det gäller en transformation TT som mappar en öppen mängd SkS\subset \mathbb{R}^k till k\mathbb{R}^k med JT0J_T\neq 0 (för alla element i SS). Om ff är kontinuerlig och har kompakt stöd i T(S)T(S) är det då tillåtet att använda formeln.

naytte 7419 – Moderator
Postad: 10 aug 01:04 Redigerad: 10 aug 01:05

Fattar! Så för att försöka sammanfatta intuitionen bakom Jakobianen:

Då vi genomför en koordinattransformation med en avbildning TT och byter mellan två områden SS och T(S)T(S) måste vi ta hänsyn till att ett areaelement i SS inte nödvändigtvis är lika stort som motsvarande areaelement i T(S)T(S). I vissa fall kan förändringen vara likformig (t.ex. om vi avbildar C:={(u,v)2:0u1,0v1}C:=\{ (u,v)\in\mathbb{R}^2:0\le u\le1, 0 \le v\le 1 \}D:={(x,y)2:0x2,0y2}D:=\{ (x,y)\in\mathbb{R}^2:0\le x\le 2, 0 \le y\le 2 \} genom variabelbytet x=2ux=2u samt y=2vy=2v), eller bero på var någonstans i det avbildade området man befinner sig (t.ex. då vi går från rektangulära koordinater till polära koordinater), varför man generellt behöver en funktion JTJ_T, Jakobianen under transformationen TT, som "tar hänsyn till" hur ett areaelement ur SS avbildas på motsvarande areaelement i T(S)T(S).

Är det korrekt tolkat?

Svara
Close