91ÊÓÆ”

De flesta matematiska problem som studeras inom teknik, eller fysik, Ă€r vĂ€lstĂ€llda. Det betyder att de har en entydig lösning som dessutom beror kontinuerligt pĂ„ indata. InnehĂ„ller data ett litet mĂ€tfel sĂ„ Ă€ndras lösningen inte sĂ€rskilt mycket. Är istĂ€llet ett problem illa-stĂ€llt betyder det ofta att lösningen helt kan Ă€ndra karaktĂ€r pĂ„ grund av ett mycket litet mĂ€tfel. Eftersom mĂ€tfel i praktiken Ă€r oundvikliga mĂ„ste lösningsmetoden konstrueras pĂ„ ett sĂ„dant sĂ€tt att man kan vara sĂ€ker pĂ„ att den lösning som berĂ€knar Ă€r vettig.

Figur 1: Den tÀnkta mÀtsituationen dÀr vi söker temperaturen pÄ en yta och mÀtningar Àr tillgÀngliga inne i materialet.


Figur 2: Temperaturfördelning under markytan i Tibet berÀknad utifrÄn mÀtningar pÄ ytan. Vi ser Àven markytan och en modell över markens egenskaper.

En tillĂ€mpning Ă€r att man ibland vill uppskatta en tidsberoende temperatur pĂ„ ytan av en kropp. Är ytan inte tillgĂ€nglig för direkta mĂ€tningar kan man istĂ€llet mĂ€ta nĂ€ra ytan och berĂ€kna den sökta temperaturen genom att lösa ett illa-stĂ€llt problem.

MÀtsituationen kan beskrivas som ett Cauchyproblem för vÀrmeledningsekvationen dÀr mÀtningar Àr tillgÀngliga vid x=1, och temperaturen söks vid x=0, se Figur 1. TillÀmpningar finns inom flera vetenskapliga sammanhang. Inom geologi vill man uppskatta temperaturen flera kilometer under markytan genom att anvÀnda mÀtningar pÄ, eller Ätminstone nÀra, ytan. HÀr löser man ett stationÀrt vÀrmeledningsproblem, med en icke-linjÀr modell, med hjÀlp av en effektiv numerisk metod som ger stabilitet i berÀkningarna. Liknande tillÀmpningar finns inom stÄlindustrin dÀr man vill studera i detalj vad som hÀnder dÄ man exempelvis vÀrmebehandlar stÄlplattor och placerar mÀtpunkter inuti materialet. 

Liknande tillÀmpningar finns för Helmholtz ekvation, som beskriver vÄgutbredning, dÀr man vill lokalisera en ljudkÀlla givet mÀtningar en bit bort. Cauchyproblemet för Helmholtz ekvation Àr illa-stÀllt och speciella metoder krÀvs. HÀr studerar jag, i samarbete med andra forskare, en klass av metoder som bygger pÄ att en viss kvadratisk form kan göras positivt definit genom att en artificiell rand introduceras, se Figur 3. Vi fÄr dÄ ett sÀtt att introducera en skalÀrprodukt, som Àr naturlig för problemet, och som gör det möjligt att implementera effektiva iterativa regulariseringsmetoder.


Figur 3: Vi söker lösningen till Helmholtz ekvation pÄ en del av randen givet Cauchydata pÄ en annan del av randen. Den kvadratiska formen blir positivt definit om vi inför artificiella inre rÀnder enligt exemplen.

Programverktyg

Programverktyg för att lösa the Sideways Heat Equation

Flera MATLAB-rutiner har skrivits för att lösa the Sideways Heat Equation. Dessa rutiner Àr skrivna för att testa idéer, och mestadels för anvÀndning inom vÄr forskningsgrupp. Men det mesta av koden Àr ganska lÀtt att förstÄ och online-dokumentationen Àr sÄdan att koden ska kunna anvÀndas med ett minimum av problem. I synnerhet finns det flera demonstrationer som syftar till att visa hur programmen ska anvÀndas. Det finns ocksÄ ett exempel dÀr verktygen tillÀmpas pÄ ett problem med faktiska uppmÀtta data.

  • En fullstĂ€ndig lista över alla program som Ă€r skrivna för att lösa the Sideways Heat Equation finns tillgĂ€nglig. [Contents.m]
  • KĂ€llkoden för alla program. Endast tillgĂ€nglig som komprimerat tararkiv. [Shetools.tar.gz]

Kontakt