Log-polära koordinater

Från Wikipedia

Log-polära koordinater (förkortning för logaritmiskt polära koordinater) är en form av koordinatsystem i två dimensioner där en punkt identifieras av två tal, ett för logaritmen på avståndet från en fix punkt, och ett för en vinkel. Log-polära koordinater är nära besläktade med polära koordinater, som traditionellt används för att beskriva områden i planet med någon form av rotationssymmetri. Inom områden som harmonisk och komplex analys är log-polära koordinater mer naturliga än polära.

Definition och koordinattransformationer[redigera | redigera wikitext]

Log-polära koordinater i planet består av ett par (ρ,θ) av reella tal, där ρ är logaritmen på avståndet till en given punkt (origo) och θ är vinkeln mellan en referenslinje (x-axeln) och linjen som går genom origo och punkten. Vinkelkoordinaten är densamma som för polära koordinater, medan sambandet för den radiella koordinaten är givet av

.

där r är avståndet till origo. Formlerna för transformering från kartesiska koordinater till log-polära koordinater är

.

och formlerna för transformering från log-polära koordinater till kartesiska är

Genom att använda komplexa tal (x,y) = x+iy ser vi att den sistnämnda koordinattransformationen är given av den komplexa exponentialfunktionen

.

Detta får som konsekvens att grundläggande ekvationer inom harmonisk analys och komplex analys får samma enkla form i log-polära koordinater som i kartesiska. Detta är inte fallet med polära koordinater.

Några ekvationer i log-polära koordinater[redigera | redigera wikitext]

Laplaces ekvation[redigera | redigera wikitext]

Laplaces ekvation i två dimensioner kan i kartesiska koordinater skrivas

medan samma ekvation i polära koordinater har det mer komplicerade utseendet

eller ekvivalent

Från sambandet följer dock att så att Laplaces ekvation i log-polära koordinater har formen

d.v.s. samma enkla form som i kartesiska koordinater. Detta är fallet för alla koordinatsystem där koordinattransformationen till kartesiska koordinater ges av en konform avbildning. För ett område i planet med rotationssymmetri, till exempel en cirkelskiva, är därför log-polära koordinater det naturliga valet då man betraktar Laplaces ekvation.

Cauchy-Riemanns ekvationer[redigera | redigera wikitext]

En snarlik situation uppstår då man betraktar analytiska funktioner. En analytisk funktion skriven i kartesiska koordinater uppfyller Cauchy-Riemanns ekvationer:

Om funktionen istället uttrycks i polär form får Cauchy-Riemanns ekvationer den mer komplicerade formen

Precis som i fallet med Laplaces ekvation, återfås den enkla formen i kartesiska koordinater genom att byta polära mot log-polära koordinater (låt ):

Cauchy-Riemanns ekvationer kan också sammanfattas i en enda ekvation som

Genom att uttrycka och i termer av och kan denna ekvation skrivas i den ekvivalenta formen

Eulers ekvation[redigera | redigera wikitext]

Då man vill lösa Dirichlets problem i ett område med rotationssymmetri tillämpar man vanligen separationsmetoden för partiella differentialekvationer på Laplaces ekvation i polär form. Det innebär att man skriver . Detta separerar Laplaces ekvation till två ordinära differentialekvationer

där är en konstant. Den första har konstanta koefficienter och är därför enkel att lösa. Den andra är ett specialfall av Eulers ekvation:

där är konstanter. Denna löses vanligen med ansatsen . Men genom att använda log-polär radie kan denna ekvation överföras till en ekvation med konstanta koefficienter:

I fallet med Laplaces ekvation gäller att och så att ekvationen för får den enkla formen

Löser man Dirichlets problem i kartesiska koordinater är det just de här två ekvationerna man får för och . Återigen ser vi att den naturliga motsvarigheten till kartesiska koordinater för ett område med rotationssymmetri inte är polära utan log-polära koordinater.

Diskret geometri[redigera | redigera wikitext]

Diskret koordinatsystem i cirkelskivan givet av log-polära koordinater (n=25)
Diskret koordinatsystem i cirkelskivan som enkelt kan uttryckas i log-polära koordinater (n=25)

När man till exempel vill lösa en partiell differentialekvation numeriskt i ett område, behöver man införa ett diskret koordinatsystem i detta område. Om området ifråga är rotationssymmetriskt och man vill ha ett rutnät bestående av rektanglar, är polära koordinater inte särskilt användbara eftersom de inte ger någon bra indelning av en cirkelskiva i rektanglar. Anledningen till detta är att man i cirkelns centrum får trianglar snarare än rektanglar. Detta kan avhjälpas med log-polära koordinater på följande sätt. Dela in planet i ett rutnät av kvadrater med en sidlängd 2/n, där n är ett positivt heltal. Använd sedan den komplexa exponentialfunktionen för att skapa ett log-polärt rutnät i planet. Vänstra halvplanet kommer då avbildas på enhetsskivan, och antalet radier är n. Om rutnätet är parallellt med koordinataxlarna och innehåller origo, får man bilden till höger. Det kan vara än mer fördelaktigt att istället avbilda diagonalerna i kvadraterna i detta rutnät, vilket ger ett diskret koordinatsystem i enhetsskivan bestående av spiraler, se bilden nedan till höger.

Dirichlet-Neumannfunktion[redigera | redigera wikitext]

Ett exempel där det sistnämnda koordinatsystemet är lämpligt är i samband med Dirichlet-Neumannfunktioner. Om det diskreta koordinatsystemet tolkas som en oriktad graf[särskiljning behövs] i enhetsskivan, kan det ses som en modell för ett elektriskt nätverk. På varje båge i grafen anges en konduktans, given av en funktion . Det elektriska nätverket kommer då ge en diskret modell av Dirichletproblemet i enhetsskivan, där Laplaces ekvation motsvaras av Kirchhoffs lag. I noderna på cirkelskivans rand läggs en elektrisk potential (Dirichlet-data), vilket ger en elektrisk ström genom dessa randnoder (Neumann-data). Den linjära funktion som avbildar Dirichlet-data på Neumann-data kallas för en Dirichlet-Neumannfunktion, och beror på nätverkets topologi samt konduktans.

I fallet med den kontinuerliga skivan, gäller att om konduktansen är likformig, säg överallt, så uppfyller Dirichlet-Neumannfunktionen följande ekvation:

För att ge en diskret modell av Dirichletproblemet som återspeglar det kontinuerliga fallet, vore det därför önskvärt att hitta en graf i enhetsskivan, vars (diskreta) Dirichlet-Neumannfunktion har samma egenskap. Medan polära koordinater inte ger något naturligt svar på denna fråga, är det just vad det spiralformade nätverket givet av log-polära koordinater förser oss med. [1]

Bildanalys[redigera | redigera wikitext]

Redan under slutet av 1970-talet gavs tillämpningar inom bildanalys för det diskreta spiralformade koordinatsystemet givet ovan. Att representera en bild digitalt i detta koordinatsystem snarare än i kartesiska koordinater ger beräkningstekniska fördelar vid rotation och zoomning i en bild. Det kan också nämnas att fotoreceptorerna i näthinnan i det mänskliga ögat är fördelade i ett mönster med stora likheter med detta spiralformade koordinatsystem. [2]

Källor[redigera | redigera wikitext]

Noter[redigera | redigera wikitext]

  1. ^ Ingerman, David V., The Square of the Dirichlet-to-Neumann Map Equals Minus Laplacian, arXiv:0806.0653, submitted on June 3 2008.
  2. ^ Weiman, Chaikin, Logarithmic Spiral Grids for Image Processing and Display, Computer Graphics and Image Processing 11, 197–226 (1979).