16 svar
408 visningar
Qetsiyah 6567 – Livehjälpare
Postad: 13 aug 2020 21:02

Flervariabelanalys: föredrar ni ∇ eller grad, div respektive rot?

Jag råkade veta om nablasymbolen och dess betydelse innan jag började plugga flervarren, jag tycker att det är konstigt att böjers undviker att använda den, den är smidig.

Något mindre pedagogiskt (man får minnas det helt enkelt), men det sparar plats och tid!

AlvinB 4014
Postad: 13 aug 2020 21:06

Jag gillar att använda \nabla överallt där det är möjligt. Låter man =(/x,/y,/z)\nabla=(\partial /\partial x,\partial/\partial y,\partial/\partial z) ser man ju direkt på uttrycken f\nabla f, ×F\nabla\times\mathbf{F} och ·F\nabla\cdot\mathbf{F} hur de beräknas, något som man inte gör om man skriver, grad\text{grad}, rot\text{rot} och div\text{div}.

Freewheeling 220 – Fd. Medlem
Postad: 13 aug 2020 21:19

Föredrar också nabla, mycket p.g.a. alla smidiga samband som involverar operatorn, se här: https://en.wikipedia.org/wiki/Vector_calculus_identities

Mycket användbart när jag pluggade elektromagnetisk fältteori minns jag :)

Qetsiyah 6567 – Livehjälpare
Postad: 13 aug 2020 21:24 Redigerad: 13 aug 2020 21:25

Åhhh den där wikipediasidan har jag sett en gång. Jag blir imponerad och förundrad, men förstår inte.

Jag minns inte i vilken bok eller hemsida jag läste det, men de påstod att nabla främst användes vid europeiska universitet. Jag kan inte uttrycka mig om påståendes sanningshalt, men... Nää?

--------------

Ett annat sätt att störa examinatorn är att skriva nabla framför funktionen, den kommuterar ju hehe. Inte kryssprodukten då!

AlvinB 4014
Postad: 13 aug 2020 21:36
Qetsiyah skrev:

Ett annat sätt att störa examinatorn är att skriva nabla framför funktionen, den kommuterar ju hehe. Inte kryssprodukten då!

Nja, det skulle jag faktiskt inte säga att det är korrekt. Det är ju faktiskt så att operatorerna /x\partial/\partial x o.s.v. inte kommuterar. Operatorn borde stå till vänster.

Studera till exempel följande exempel, där u\mathbf{u} är ett vektorfält och ff är ett skalärfält:

(·u)f(\nabla\cdot\mathbf{u})f

och

(u·)f(\mathbf{u}\cdot\nabla)f

Dessa två uttryck är inte lika med varandra. Ser du varför?

Qetsiyah 6567 – Livehjälpare
Postad: 13 aug 2020 21:41 Redigerad: 13 aug 2020 21:54

Är de inte? 

=f(n=1?unxn)

u·=·u och f är en skalär? Eller? Har min hjärna lämnat kontoret såhär tidigt på kvällen?

EDIT: jag kollade wikipedia, jag förstår. Isåfall betyder ·u ingenting alls eftersom det inte finns någon funktion framför att derivera. Detta visste jag inte.

EDIT2: Men skalärprodukten måste ju kommutera, vad spelar det för roll om nablas egenskaper?

EDIT3: Isåfall är nabla inte en vektor, eller så är detta inte en skalärprodukt

AlvinB 4014
Postad: 13 aug 2020 21:57 Redigerad: 13 aug 2020 21:59

Differentialoperatorn deriverar bara saker som står till höger om den.

fx\dfrac{\partial f}{\partial x} är alltså derivatan av ff med avseende på xx, medan i fallet fxf\dfrac{\partial}{\partial x} "väntar" differentialoperatorn fortfarande på att en funktion ska sättas till höger om den så att den kan deriveras.

I mina exempel resulterar det i följande (nu låtsas jag att u\mathbf{u} och ff är tredimensionella för enkelhetens skull):

·uf=(uxx+uyy+uzz)f=uxxf+uyyf+uzzf\left(\nabla\cdot\mathbf{u}\right)f=(\dfrac{\partial\mathbf{u}_x}{\partial x}+\dfrac{\partial\mathbf{u}_y}{\partial y}+\dfrac{\partial\mathbf{u}_z}{\partial z})f=\dfrac{\partial\mathbf{u}_x}{\partial x}f+\dfrac{\partial\mathbf{u}_y}{\partial y}f+\dfrac{\partial\mathbf{u}_z}{\partial z}f

u·f=(uxx+uyy+uzz)f=uxfx+uyfy+uzfz\left(\mathbf{u}\cdot\nabla\right)f=(\mathbf{u}_x\dfrac{\partial}{\partial x}+\mathbf{u}_y\dfrac{\partial}{\partial y}+\mathbf{u}_z\dfrac{\partial}{\partial z})f=\mathbf{u}_x\dfrac{\partial f}{\partial x}+\mathbf{u}_y\dfrac{\partial f}{\partial y}+\mathbf{u}_z\dfrac{\partial f}{\partial z}

I det första fallet slutar det med att vi deriverar komponenterna i u\mathbf{u} medan i det andra fallet slutar det med vi att vi deriverar ff, eftersom deriveringsoperatorn bara deriverar det som står till höger om den.

Nablaoperatorn är av denna inledning inte kommutativ. Det gäller att skilja på operatorer och tal. Differentialoperatorn är inget vanligt tal, vilket gör att nablaoperatorn inte heller är en vanlig vektor. Därför kan skalärprodukter med nabla inte förväntas ha samma egenskaper som en vanlig skalärprodukt.

Qetsiyah 6567 – Livehjälpare
Postad: 13 aug 2020 22:05 Redigerad: 13 aug 2020 22:11

Så vilken (eller båda) av mina påståenden gäller? 1) nabla är inte en vektor 2) pricken betecknar inte skalärprodukt.

Det kanske inte är en intressant frågeställning, och nabla och pricken är bara konvention/räkneregel?

EDIT: är (u·) isåfall en operator?

AlvinB 4014
Postad: 13 aug 2020 22:21
Qetsiyah skrev:

Så vilken (eller båda) av mina påståenden gäller? 1) nabla är inte en vektor 2) pricken betecknar inte skalärprodukt.

Det kanske inte är en intressant frågeställning, och nabla och pricken är bara konvention/räkneregel?

Ettan är definitivt sann. Nabla är inte en vektor, utan en vektoroperator precis som differentialoperatorn inte är ett tal, utan en operator. Sen är frågan om man kan utföra skalärprodukt med något som inte är en vektor. Jag skulle säga nej, eftersom skalärprodukten då slutar kommutera, så vi får väl säga att tvåan också är sann och kalla pricken för något annat.

Man skulle så klart som du säger kunna se alltsammans som en fiffig minnesregel, men jag tycker faktiskt att operatorer är såpass rigoröst definierade att man kan använda dem "på riktigt". Man får bara vara uppmärksam på när deras räkneregler skiljer sig från de vanliga.

Freewheeling 220 – Fd. Medlem
Postad: 13 aug 2020 22:21 Redigerad: 13 aug 2020 22:24
Qetsiyah skrev:

Så vilken (eller båda) av mina påståenden gäller? 1) nabla är inte en vektor 2) pricken betecknar inte skalärprodukt.

Det kanske inte är en intressant frågeställning, och nabla och pricken är bara konvention/räkneregel?

EDIT: är (u·) isåfall en operator?

Jag skulle säga att båda dina påståenden är korrekta. \nabla är en operator och u·\mathbf{u} \cdot \nabla är också en operator, som båda kan ta en skalär funktion som input. Nabla och pricken är som du säger bara en konvention/räkneregel. Att man använder just samma notation som man gör för för skalärprodukten är p.g.a. hur suggestiv den är, då man förstår precis hur man ska "räkna" med operatorn, som AlvinB visar.

Freewheeling 220 – Fd. Medlem
Postad: 13 aug 2020 22:28
Qetsiyah skrev:

Jag minns inte i vilken bok eller hemsida jag läste det, men de påstod att nabla främst användes vid europeiska universitet. Jag kan inte uttrycka mig om påståendes sanningshalt, men... Nää?

Hmm, i vår lärobok i den s.k. "TET:en" (kurs i elektromagnetisk fältteori) så användes nabla-operatorn väldigt friskt i alla fall och det var en amerikansk bok. Det är ju i och för sig en fysikkurs :)

Qetsiyah 6567 – Livehjälpare
Postad: 14 aug 2020 00:52

Ja, pricken och speciellt ×är VÄLDIGT suggestiva om hur man räknar med de. När jag clickar mig vidare på "vector operator" i wikipedia finns inte många andra exempel på det (det fanns laplacianen och d'Alembert operatorn).

Jag får nog nöja mig med att säga att nabla är en operator som uppvisar vektorliknande egenskaper, men är inte en vektor.

Jroth 1191 – Fd. Medlem
Postad: 14 aug 2020 03:19 Redigerad: 14 aug 2020 03:25

Låt oss säga att vektorn T\nabla T är gradienten till ett temperaturfält T i ett område.

Då talar gradienten om i vilken riktning temperaturen ändrar sig som mest. Längden av vektorn talar om hur stor förändringen är.

 I en viss given punkt, oavsett i vilket koordinatsystem vi mäter, borde fältet ändra sig som mest i en given riktning (och med en given storlek). Temperaturen vet ju inte om i vilket koordinatsystem vi valt att studera systemet. Fältet borde alltså ge oss samma resultat, oavsett koordinatsystem.

Men hur uppför sig egentligen den "vektorlika" vektoroperatorn ∇ i andra koordinatsystem?

Och hur kan vi konstruera vektorstorheten T\nabla T så att den blir oberoende av valet av koordinatsystem?

Qetsiyah 6567 – Livehjälpare
Postad: 14 aug 2020 22:47 Redigerad: 14 aug 2020 22:47

Ja... det är det jag ska komma till! Jag har köpt boken vektoranalys. Det verkar vara viktigt, olika basvektorer och koordinatsystem. Jag hoppas att det blir ett annat tillgängligare sätt att förstå mig på tensorn också, men det är långt bak i boken.

Om du vill så avslöja svaren på dina frågor själv! Jag sitter verkligen inte på svaren.

Jroth 1191 – Fd. Medlem
Postad: 16 aug 2020 14:15 Redigerad: 16 aug 2020 14:44

Ja, det har med vektoranalys, basvektorer och tensorer att göra :)

Antag att vi lägger ut två alternativa koordinatsystem i rummet. I det vänstra är det 1m mellan strecken. I det högra är det 2m mellan strecken. Vi kallar det vänstra för xyxy-systemet och det högra för x'y'x^{'}y^{'}-systemet. En transformation mellan systemen ser då ut så här:

x=2x'x=2x^{'}

y=2y'y=2y^{'}

I det vänstra systemet är gradienten 

T=nTxnen=Txi^+Tyj^\displaystyle \nabla T=\sum_n\, \frac{\partial T}{\partial x^n}\mathbf{e} _n=\frac{\partial T}{\partial x}\mathbf{\hat{i}}+\frac{\partial T}{\partial y}\mathbf{\hat{j}}

I det högra systemet är gradienten:

T'=nT'xne'n=T'x'i^'+T'y'j^'\displaystyle \nabla {T^{'}}=\sum_n\, \frac{\partial T^{'}}{\partial x^n} \mathbf{e^{'}}_n=\frac{\partial T^{'}}{\partial x^{'}}\mathbf{\hat{i}}^{'}+\frac{\partial T^{'}}{\partial y^{'}}\mathbf{\hat{j}}^{'}

Basvektorerna är tydligen dubbelt så långa; för att garantera en invariant storlek på gradienten måste därför derivatorna halveras. Gör de det?

Tyvärr inte! I det högra systemet ändrar sig saker istället dubbelt så fort. Om man tänker på det gäller ungefär:

ΔF'Δx'=2ΔFΔx\frac{\Delta F^{'}}{\Delta x^{'}}=2\frac{\Delta F}{\Delta x}

Dessvärre verkar alltså storleken av gradienten i det högra systemet vara fyra gånger så stor.

En vanlig vektor har vett nog att under transformation röra sig åt motsatt håll (kontra) en förändring av basvektorernas längd. Men gradienten tycks vara motvalls! Den rör sig MED en förändring och gör saken VÄRRE (co-moving). Den tycks komma från ett annat universum, the upside-down.

Hur lagar vi det då? För att svara på den frågan behöver man egentligen en ganska omfattande matematisk apparat, men syftet med det här inlägget är bara att ge en känsla för hur det hänger ihop.

Det visar sig (i kursen linjär algebra) att en transformation av en" vanlig" vektor vv ges av en matrismultiplikation med  transformationens jakobian (J),

v=Jv'v=Jv^{'}

Vill man kan man skriva matrismultiplikationen som en summa (för varje komponent j):

vj=xjxh'vh'=Jh'jvh'v^j=\frac{\partial x^j}{\partial x^{h^'}}v^{h^'}=J^j_{h^'}v^{h^'}

Här använder vi summakonventionen över indexet h'h^{'}.

Vektorer som uppför sig så kallar vi kontravarianta

Den andra sorten, scary creatures from the upside-down, kallar vi således kovarianta. De transformerar enligt

vk=xm'xkvm'=Jkm'vm'v_k=\frac{\partial x^{m^'}}{\partial x^{k}}v_{m^'}=J^{m^'}_kv_{m^'}

Nu visar det sig (med kedjeregeln) att gradienten i det högra systemet kan uttryckas som

h'T'=Txh'=Txkxkxh'=Jh'kkT\nabla_{h'} T^{'}=\frac{\partial T }{\partial x^{h^'} }=\frac{\partial T}{\partial x^k}\frac{\partial x^k}{\partial x^{h^'}}=J^k_{h^'}\nabla_k T

Så komponenterna i gradienten transformerar som en kovariant vektor.

För att teckna gradienten som en fristående vektor måste vi använda dualbasen:

T=Txkek\boxed{\nabla T=\frac{\partial T}{\partial x^k}\mathbf{e}^k}

I ett vanligt  kartesiskt system sammanfaller baserna, dvs ei=ei\mathbf{e}_i=\mathbf{e}^i och gradienten och dess nabla-kusiner uppför sig bedrägligt normalt i uttryck som T\nabla T och ·F\nabla \cdot \mathbf{F}.

Men låt dig inte luras, gradienten är en invånare från the upside-down!

Qetsiyah 6567 – Livehjälpare
Postad: 17 aug 2020 00:55

Okej...

I det högra systemet är basvektorerna dubbelt så långa samtidigt som derivatorna är dubbelt så stora? Alltså gradienten är *4?

Varför skriver du T först sedan F?

En vanlig vektor har vett nog att under transformation röra sig åt motsatt håll (kontra) en förändring av basvektorernas längd. Men gradienten tycks vara motvalls! 

Alltså att om vi vill skriva koordinater på en vektor så blir koordinaterna hälften i högra systemet? Verkar vettigt.

Vad heter denna matematiska apparat?

Kontravariant är så som jag förväntar mig att en vektor ska bete sig, ja. 

Jag har bara snuddat vid Einsteins summakonvention, kan du uttolka dina ekvationer i slutet?

Dualbas? Bas för dualrummet? Är det inte en mängd vektorer? Det du skrivit är inte ens en mängd, vad menas?

Jroth 1191 – Fd. Medlem
Postad: 17 aug 2020 20:28 Redigerad: 17 aug 2020 21:34

Summakonventionen betyder bara att om du hittar två index som heter samma, ett där upp och ett där nere, ska man summera dessa över 1n1\dots n där nn är dimensionen. Låt oss titta på den "nya" gradienten:

T=Txheh\nabla T=\frac{\partial T}{\partial x^h}\mathbf{e}^h

Här förekommer hh två gånger, alltså ska vi summera. Eftersom n=2n=2 i vårt exempel blir det bara två termer. Låt oss studera gradienten i x'y'x^'y^'-systemet:

T(x',y')=T(x',y')x'ex'+T(x',y')y'ey'\nabla T(x',y')=\frac{\partial T(x^{'},y^{'})}{\partial x^{'}}\mathbf{e}^{x^{'}}+\frac{\partial T(x^{'},y^{'})}{\partial y^{'}}\mathbf{e}^{y^{'}}

Utan att närmare specificera rummet E\mathbf{E} påstår vi att alla kontravarianta vektorer bor där. Om {e1en}\{e_1\dots e_n\} är en ordnad bas för E\mathbf{E} finns en ordnad bas för dualrummet E*E^* där de kovarianta vektorerna bor. Denna bas {e1en}\{e^1 \dots e^n\} uppfyller <ej,ei>=δij<e^j,e_i>=\delta^j_i. Vidare gäller för vEv\in\mathbf{E} och αE*\alpha\in\mathbf{E}^*

v=<ei,v>eiv=<e^i,v>e_i samt α=<α,ei>ei\alpha=<\alpha,e_i>e^i

En tangentbas för vårt exempelsystem är ex'=(2,0),ex'=(0,2)e^{'}_x=(2,0),\, e^{'}_x=(0,2)

Om vi låter parbildningen mellan E\mathbf{E} och E*\mathbf{E}^* vara den vanliga skalärprodukten ska tangentbasen tillsammans med dualbasen uppfylla:

ei·ej=δji \mathbf{e}^i\cdot \mathbf{e}_j=\delta^i_j

Därför måste dualbasen (eller normalbasen, eller den kontravarianta basen) i vårt exempel vara

ex'=(12,0),ey'=(0,12)\mathbf{e}^{x'}=(\frac{1}{2},0),\, \mathbf{e}^{y'}=(0,\frac{1}{2})

Ett annat sätt att hitta normalbasen under samma förutsättningar är att bilda gradienten av x''(x,y)x''(x,y) och y''(x,y)y''(x,y) med avseende på x,yx,y.

Notera att man ibland kallar tangenbasen ei\mathbf{e}^i den kontravarianta basen, trots att den är en bas för rummet där de kovarianta vektorerna bor.

Låt oss som avslutning beräkna gradienten i en punkt (2,4) till en funktion T(x,y)=x2+y2T(x,y)=x^2+y^2

I vårt primmade system får vi

T'(x',y')=(2x')2+(2y')2T^{'}(x^{'},y^{'})=(2x^{'})^2+(2y^{'})^2

(x',y')=J(2,4)=(1,2)(x^{'},y^{'})=J(2,4)=(1,2)

T'=8x'ex'+8y'ey'=8ex'+16ey'=(4,8)\nabla T^{'}=8x^{'}\mathbf{e}^{x'}+8y^{'}\mathbf{e}^{y'}=8\mathbf{e}^{x'}+16\mathbf{e}^{y'}=(4,8)

 

Som rolig övning kan du testa att beräkna gradienten i polära koordinater

Tangentbasen för polära koordinater är er=r^,eθ=rθ^e_r=\hat{r},\, e_\theta=r\hat{\theta}

Vad blir dualbasen er,eθe^r,e^\theta? Hur ser T(r,θ)T(r,\theta) ut?

Vad är punkten (2,4)(2,4) uttryckt i polära koordinater? Vad blir gradienten i punkten (2,4)(2,4)? Stämmer det med tidigare resultat?

Svara
Close