Barycentriska koordinater

Från testwiki
Hoppa till navigering Hoppa till sök

Inom geometri betecknar barycentriska koordinater (från grekiska βαρύς, barys, "tung" och κέντρον, kentron, "centrum") en uppsättning av n+1 tal, vilka anger en punkts läge i förhållande till en n-dimensionell simplex (sträcka, triangel, tetraeder, etcetera) i det n-dimensionella rummet genom att ange relativa vikter som, om de placeras i hörnen på denna simplex, gör punkten till simplexens geometriska tyngdpunkt. I allmänhet avses läget av en punkt i planet i förhållande till en triangel.[1][2] De skall inte förväxlas med begreppet "barycentrum" som används inom astronomi för att ange den gemensamma tyngdpunkten för en uppsättning himlakroppar (även om begreppet är närbesläktat). Barycentriska koordinater infördes av August Ferdinand Möbius 1827 i Der Barycentrische Calcul.[3]

Barycentriska koordinater skrivs vanligtvis separerade av kolon (exempelvis α:β:γ för en punkt i planet i förhållande till en triangel i samma plan).

Om alla koordinaterna är större än noll ligger punkten innanför simplexens begränsningar och är en eller flera koordinater noll ligger punkten på begränsningarna. Alla koordinater kan inte vara noll. Är någon koordinat negativ ligger punkten utanför simplexen (det motsvarar att en "negativ vikt", eller en "lyftkraft", måste placeras i hörnet). Någon koordinat måste ha ett positivt värde.

De barycentriska koordinaterna är relativa, vilket innebär att endast deras inbördes förhållanden spelar roll: 2:1:0 är detsamma som 4:2:0 eller 100:50:0.

Med absoluta barycentriska koordinater menas att koordinaterna normerats så att deras summa blir lika med ett. För att normera koordinaterna delar man dem med deras summa. Exempelvis om koordinaterna 2:1:0 divideras med summan av dem (2+1+0=3) får vi de absoluta barycentriska koordinaterna 23:13:0.[4]

Inom astronomi används termen barycentriskt koordinatsystem för att ange ett koordinatsystem (sfäriskt eller kartesiskt) med origo i systemets tyngdpunkt (exempelvis solsystemets tyngdpunkt).

Barycentriska koordinater i en dimension

Figur 1. Endimensionella barycentriska koordinater är väldigt triviala, men mönstret följer med när antalet dimensioner ökar.

Endimensionella barycentriska koordinater beskriver läget av en punkt på en linje i förhållande till en given sträcka på linjen. Låt oss kalla sträckans ändpunkter för A och B som i figur 1. De barycentriska koordinaterna anger då hur mycket massa vi skall placera i (eller hur mycket kraft vi skall applicera på) A respektive B för att tyngdpunkten skall befinna sig i P relativt sett. Placerar vi all massa i A ligger tyngdpunkten i A och P ligger alltså i A och har, exempelvis, koordinaterna 1:0 (och motsvarande för B såklart). För att ange de barycentriska koordinaterna för P skall vi alltså beräkna två krafter applicerade i A respektive B (vilket ger oss koordinaterna FA:FB) som ger ett motverkande vridmoment i förhållande till P, det vill säga att FAAP=FBBP=FBPB, vilket ger oss koordinaterna FA:FAAPPB. Eftersom endast koordinaternas relativa värden är av intresse kan vi multiplicera dem med PBFA, vilket ger oss de likvärdiga koordinaterna PB:AP. Dessa kan "normeras" genom att dividera dem med AP+PB=AB varvid deras summa blir lika med ett. Vi ser att respektive koordinat är proportionell mot det "riktade" avståndet från den andra ändpunkten. I det fall man anger de fakiska avstånden PB:AP talar man om de homogena barycentriska koordinaterna med avseende på sträckan AB.[5]

Barycentriska koordinater i två dimensioner

Figur 2.
Absoluta barycentriska koordinater för vissa punkter i en liksidig respektive rätvinklig triangel.

Tvådimensionella barycentriska koordinater beskriver läget av en punkt P i planet i förhållande till en triangel ABC i samma plan (se figur 2). Genom att placera tre "vikter" i de tre triangelhörnen (eller applicera tre krafter på hörnen) skall vi "balansera" triangeln i P.

Vi börjar med att placera all vikt FP i P. Därefter balanserar vi vikten längs linjen AL (i enlighet med resonemanget för endimensionella koordinater) så att vi placerar FA=FPPLAL i A och FL=FPAPAL i L, varefter vi balanserar vikten i L längs BC genom att flytta FB=FLLCBC till B och FC=FLBLBC till C.

Detta innebär att L fortfarande är tyngdpunkt på BC, vilket innebär att triangeln fortfarande är balanserad i P och att all vikt befinner sig fördelad på A, B eller C.

Vi noterar av det ovanstående att FBFC=LCBL. På samma sätt kan vi visa att FAFB=NBAN och FAFC=MCAM.[6]

Detta leder till att BPC:APC:APB (i det fall arean av en deltriangel ligger helt utanför ABC är dess area negativ.) är barycentiska koordinter för P, eftersom

FBFC=LCBL=CLABLA=CLPBLP=CLACLPBLABLP=APCAPB

och på samma sätt är FAFB=BPCAPC och FAFC=BPCAPB.

Detta ger oss att

FA:FB:FC=FA:FAAPCBPC:FAAPBBPC=

=BPCFA(FA:FAAPCBPC:FAAPBBPC)==BPC:APC:APB.

Dessa kallas de homogena barycentriska koordinaterna relativt triangeln ABC.[7][8] Vi kan normera dem genom att dividera var och en av dem med ABC så att deras summa blir lika med ett, vilket ger oss de absoluta barycentriska koordinaterna (vilka även kallas areal coordinates, "areella koordinater", på engelska[9]):

BPCABC:APCABC:APBABC

Barycentriska koordinater i tre eller flera dimensioner

Figur 3.

Det tredimensionella fallet: Analogt med det tvådimensionella fallet (i vilket vikterna i triangelhörnen är proportionella mot respektive "motstående deltriangels" area) är vikten i respektive tetraederhörn proportionell mot volymen av dess "motstående deltetraeders" volym.

Resonemang
Betrakta tetraedern ABCD och punkten G i figur 3. Placera all vikt i G och fördela sedan denna vikt på A och O så att G fortfarande är tyngdpunkt. Den tilldelade vikten i O fördelas sedan på hörnen i triangeln BCD (se ovan under två dimensioner) så att denna triangel är balanserad i O och varvid tetraedern fortfarande balanserar i G. För areorna av deltrianglarna BCO, BDO och CDO gäller
BCOμD=BDOμC=BCOμB
där μX betecknar vikten i hörnet X. Dessa deltrianglar utgör baserna för tetraedrar med det fjärde hörnet i A, vilka sålunda har volymer som är proportionella mot bastriangelns yta (och därmed mot vikten i det "motstående triangelhörnet"). På samma sätt är volymerna av de tre tetraedrarna med de tre deltrianglarna som basytor och det fjärde hörnet i G proportionella mot vikten i respektive "motstående hörn". Detta ger
ABCOGBCOμD=ABDOGBDOμC=ACDOGCDOμB
ABCGμD=ABDGμC=ACDGμB(=BCDGμA).
Det vill säga att vikterna i tetraederhörnen är proportionella mot volymerna av de "motstående deltetraedrarna", vilka har den till hörnet motstående tetraedersidan som basyta och det fjärde hörnet i G (detta gäller såklart även för μA och BCDG, vilket ju visas på samma sätt genom att utgå från ett annat hörn och dess motstående tetraedersida).
Punkten G har alltså de barycentriska koordinaterna μA:μB:μC:μD i förhållande till tetraedern ABCD. De absoluta barycentriska koordinaterna för G är BCDGABCD:ACDGABCD:ABDGABCD:ABCGABCD, eftersom deras summa ju är lika med ett.

Fler dimensioner: Ökar vi på antalet dimensioner förfar vi på samma sätt. För fyra dimensioner balanserar vi först mellan ett hörn och skärningspunkten (för linjen genom hörnet och G) med den till hörnet motstående tetraedern, varefter vi balanserar denna tetraeders tilldelade vikt på dess fyra hörn enligt ovan. Vi kan fortsätta att öka på med en dimension i taget på samma sätt. Detta innebär att vi för en n-dimensionell simplex A1A2A3...An och en punkt G får

A2A3...AnGμA1=A1A3...AnGμA2=....
Punkten G har de barycentriska koordinaterna μA1:μA2:μA3:...:μAn i förhållande till simplexen A1A2A3...An.

Ortsvektorer och kartesiska koordinater

Figur 4.

I figur 4 visas en triangel ABC och en punkt P som har de absoluta barycentriska koordinaterna μA:μB:μC. Vi noterar att, i enlighet med endimensionella koodinater ovan, APAB=μBAB och APAC=μCAC. Därför har vi också att AP=APAB+APAC=μBAB+μCAC. Om nu triangelhörnen har ortsvektorerna OA, OB respektive OC i förhållande till en punkt O får vi (i sista ledet utnyttjas att de absoluta koordinaternas summa är ett):

OP=OA+AP=OA+μBAB+μCAC=
=OA+μB(OBOA)+OA+μC(OCOA)==(1μBμC)OA+μBOB+μCOC==μAOA+μBOB+μCOC

Från barycentriska till kartesiska koordinater

Uttrycket för punktens ortsvektor ger direkt att, om triangelhörnen har de kartesiska koordinaterna A=(xA,yA), B=(xB,yB) respektive C=(xC,yC), så är de kartesiska koordinaterna för P:

P=(xP,yP)=(μAxA+μBxB+μCxC, μAyA+μByB+μCyC)

Speciellt märker vi att om triangelhörnen har de kartesiska koordinaterna (0,0), (1,0) och (0,1) så har punkten (x,y) de absoluta barycentriska koordinaterna 1xy: x: y.

Från kartesiska till barycentriska koordinater

Om vi skriver om uttrycken för punktens kartesiska koordinater och utnyttjar att μC=1μAμB får vi två ekvationer med två obekanta:

(xAxC)μA+(xBxC)μB+xCxP=0
(yAyC)μA+(yByC)μB+yCyP=0

vilka har lösningen

μA=(yByC)(xPxC)+(xCxB)(yPyC)(yByC)(xAxC)+(xCxB)(yAyC) och
μB=(yCyA)(xPxC)+(xAxC)(yPyC)(yByC)(xAxC)+(xCxB)(yAyC).
(ur vilka vi även får μC=1μAμB)

Trilinjära koordinater

Låt a, b och c beteckna längden av de motstående sidorna till triangelhörnen. En punkt med de trilinjära koordinaterna x:y:z har då de barycentriska koordinaterna ax:by:cz.[10]

Omvänt har därför en punkt med de barycentriska koordinaterna α:β:γ de trilinjära koodinaterna αa:βb:γc.

De homogena barycentriska koordinaterna α:β:γ motsvaras av de exakta trilinjära koordinaterna a:b:c=2αa:2βb:2γc.[8]

Bevis
α:β:γ=|BCP|:|ACP|:|ABP|=aa2:bb2:cc2=aa:bb:cc=ax:by:cy
Om α:β:γ är homogena är de alltså lika med aa2:bb2:cc2, och omvänt är de exakta trilinjära koordinaterna a:b:c=2αa:2βb:2γc.


Barycentriska koordinater för vissa punkter

Den geometriska tyngdpunkten

Triangelns geometriska tyngdpunkt har barycentriska koordinater

1:1:1
Bevis
Beviset följer direkt ur att den geometriska tyngdpunkten är medianernas skärningspunkt och att medianerna delar triangeln i sex likstora trianglar.

Den inskrivna cirkelns medelpunkt

Den inskrivna cirkelns medelpunkt har barycentriska koordinater som kan skrivas som

BC : AC : AB
Bevis
Den inskrivna cirkelns medelpunkt P har samma avstånd till triangelns sidor, dess radie r. Radien är lika med deltrianglarnas höjd och de barycentriska koordinaterna är därför
BCP:ACP:ABP=rBC2:rAC2:rAB2=BC : AC : AB

De vidskrivna cirklarnas medelpunkter

Medelpunkterna för de vidskrivna cirklarna till BC, AC och AB har de barycentriska koordinaterna

|BC|:|AC|:|AB|,
|BC|:|AC|:|AB| respektive
|BC|:|AC|:|AB|.
Bevis
Eftersom medelpunkten för den vidskrivna cirkeln, med radien r, till BC har de trilinjära koordinaterna
r:r:r=1:1:1
har den de barycentriska koordinaterna
|BC|:|AC|:|AB|
i enlighet med avsnittet om trilinjära koordinater ovan. Motsvarande gäller de båda övriga vidskrivna cirklarnas medelpunkter.

Den omskrivna cirkelns medelpunkt

Figur 5.

Den omskrivna cirkelns medelpunkt har barycentriska koordinater som kan skrivas som

BCcosBAC : ACcosCBA : ABcosBCA
Bevis
Den omskrivna cirkelns medelpunkt P har samma avstånd till triangelhörnen, cirkelns radie r. Låt F vara fotpunkt åt PBC (figur 5). Vi har då
BCP=BFP+CFP=2BFP=BFFP=
 =BFrcosBPF=BFrcosBPC2=
 =BFrcosBAC=BC2rcosBAC
I näst sista ledet utnyttjas randvinkelsatsen. (BC är ju en korda i den omskrivna cirkeln och vinkeln i medelpunkten är enligt denna sats dubbelt så stor som vinkeln i en punkt på omkretsen.)
Samma resonemang för de båda andra triangelsidorna ger oss barycentriska koordinater
BC2rcosBAC : AC2rcosCBA : AB2rcosBCA=
 =BCcosBAC : ACcosCBA : ABcosBCA

Ortocentrum

Ortocentrum har barycentriska koordinater som kan skrivas som

BCcosBAC:ACcosCBA:ABcosACB

eller

BCcosCBAcosACB:ACcosACBcosBAC:BCcosBACcosCBA

Det andra uttrycket erhålls genom multiplikation med cosBACcosCBAcosACB

Bevis
Vi utnyttjar förhållandet mellan barycentriska koordinater och trilinjära koordinater (se avsnittet ovan).
Den omskrivna cirkelns medelpunkt har barycentriska koordinater (se ovan)
 BCcosBAC : ACcosCBA : ABcosBCA
och därför trilinjära koordinater
cosBAC:cosCBA:cosACB
Ortocentrum är isogonalkonjugat till den omskrivna cirkelns medelpunkt och dess trilinjära koordinater är därför
1cosBAC:1cosCBA:1cosACB
vilka motsvarar de barycentriska koordinaterna
BCcosBAC:ACcosCBA:ABcosACB

Symmedianpunkten

Symmedianpunkten har de barycentriska koordinaterna

BC2:AC2:AB2
Bevis
För beviset utnyttjar vi sambandet mellan barycentriska och trilinjära koordinater och att symmedianpunkten är isogonalkonjugat till den geometriska tyngdpunkten.
Den geometriska tyngdpunkten har barycentriska koordinater
1:1:1
vilket motsvarar de trilinjära koordinaterna
1BC:1AC:1AB
Symmedianpunkten har, eftersom den är isogonalkonjugat till den geometriska tyngdpunkten, de trilinjära koordinaterna
BC:AC:AB
vilka motsvarar de barycentriska koordinaterna
BC2:AC2:AB2

Barycentrisk interpolation

En Gouraud-skuggad triangel RGB i vilken punkternas färg (Röd, Grön, Blå) är fP=μRfR+μGfG+μBfB (där fR = 100% röd, 0% grön, 0% blå, etcetera).

För en funktion av två variabler, f(x,y), med de kända värdena fA=f(xA,yA), fB=f(xB,yB) och fC=f(xC,yC) för hörnen i triangeln ABC kan en linjär interpolation av värdet i en punkt P=(xP,yP) med de absoluta barycentriska koordinaterna μA:μB:μC göras enligt fP=f(xP,yP)μAfA+μBfB+μCfC. Barycentrisk interpolation kan enkelt utsträckas till fler dimenensioner.[11] Genom att skapa ett nät av trianglar (ett så kallat "mesh"), eller simplexar av högre dimension, kan beräkningar göras för större områden (ett exempel är beräkningar av isolinjer eller värden för olika platser i ett nät av väderstationer). Mer förfinade interpolationer kan göras med polynomapproximationer i stället för linjära sådana.[12]

Barycentrisk interpolation och generaliseringar av denna till godtyckliga polygoner och polyedrar används inom flera områden, exempelvis finita elementmetoden (FEM), och speciellt noterbart är applikationer inom datorgrafik för exempelvis skuggning och animation.[13]

Se även

Referenser

  1. Clifford A. Pickover, 2015, 250 milstolpar i matematikens historia från Pythagoras till 57:e dimensionen, sid. 222. Mall:ISBN
  2. Mall:MathWorld
  3. August Ferdinand Möbius, 1827, Der Barycentrische Calcul, Verlag von Johann Ambosius Barth, Leipzig.
  4. Marie-Nicole Gras, 2014, Distances between the circumcenter of the extouch triangle and the classical centers, Forum Geometricorum, 14, sid. 52.
  5. Paul Yiu, 2013, Introduction to the Geometry of the Triangle Mall:Wayback, sid. 1. Department of Mathematics, Florida Atlantic University.
  6. Med dessa tre uttryck bevisar man dessutom enkelt Cevas sats: LCBLNBANAMMC=1.
  7. Paul Yiu, 2000, The uses of homogeneous barycentric coordinates in plane euclidean geometry Mall:Wayback i "International Journal of Mathematical Education in Science and Technology" 31, sid. 570
  8. 8,0 8,1 Mall:MathWorld
  9. Mall:MathWorld
  10. Matthew Harvey, 2015, Geometry Illuminated: An Illustrated Introduction to Euclidean and Hyperbolic Plane Geometry, sid. 349. Mall:ISBN
  11. Remi Munos, Andrew Moore, 1998, Barycentric Interpolators for Continuous Space & Time Reinforcement Learning i "Proceedings of the 1998 conference on Advances in neural information processing systems II", sid. 1024– 1030, Cambridge, MA, USA, 1999. MIT Press.
  12. Tatiana V. Voitovich, Stefan Vandewalle, Barycentric Interpolation and Exact Integration Formulas for the Finite Volume Element Method i "Numerical analysis and applied mathematcs, AIP Conference Proceedings, 1048", (Simos, T., Psihoyios, G., Tsitouras, C. (Eds.)). International Conference on Numerical Analysis and Applied Mathematics. Kos, Grekland, 16 - 20 september 2008", sid. 575-579.
  13. Michael S Floater, Jiři Kosinka, 2010, Barycentric interpolation and mappings on smooth convex domains i "Proceedings of the 14th ACM Symposium on Solid and Physical Modeling - SPM '10".

Mall:Koordinater