HH/ITE/BN Matematisk modellering och Mathematica 1 Något om Matematisk modellering och Mathematica Bertil Nilsson 2015-08-15 t Vtct qin t cin t qut t cut t c 0 c0 t 0, Τ , V t A k1 B k2 k3 C A t B t k1 A k2 B C t k2 B k3 C k1 A k3 C 2 Matematisk modellering och Mathematica HH/ITE/BN ť Förord På följande sidor presenteras en elementär "streetwise guide" till matematisk modellering med med något lite användning av Mathematica. Framställningen är fåordig, fri från pedanteri men i någon mening fullständig. Det man väsentligen behöver veta om begrepp, terminologi, beteckningar och teori för att modellera och lösa problem i framtida kurser och yrkesliv som ingenjör, naturvetare eller lärare klarläggs och typiska exempel ges. ť Introduktion Matematiken har för de flesta människor en undanskymd och inte sällan ifrågasatt tillvaro. Vi har ju datorer!! I själva verkat har den en allt viktigare, om än lite osynlig, betydelse i våra liv. I en konkurrensutsatt verksamhet gäller det för en modern ingenjör att minska ledtider och spara in på kostsamma provserier. Därför är all produktutveckling numera helt beroende av att kunna simulera sina produkter i dator under utvecklingsfasen. Framgångsrik utveckling av exempelvis bilar, datorer, mobiltelefoner, vindkraftverk och medicinsk utrustning vore otänkbart utan stöd av avancerad matematik! Detsamma gäller för mera "osynliga saker", såsom effektiv schemaläggning av flygtrafik (tåg, sopbilar, elnät, mobilnät...) eller sökmotorer (Google) på www. Att modellera verkliga problem med hjälp av matematik kräver erfarenhet och ett öppet sinnelag. Till sin natur är denna verksamhet oftast helt skild från den matematikundervisning man möter i skolan och upplevs av naturliga skäl som svår. Inte minst på grund av att man ska hämta kunskap och inspiration från flera grundläggande ämnen, inte bara matematik även om detta kommer att bli själva språket. Det gäller att "se" på verkligheten genom glasögon som identifierar och tydliggör fenomen som kan kläs med matematiska begrepp som derivata, integral och differentialekvation. Med detta betraktelsesätt tillägnar sig en modern ingenjör en klar konkurrensfördel eftersom modellering och simulering är ett absolut krav för att klara de allt kortare ledtider som krävs för att utveckla nya produkter. Modellering är inte något speciellt för matematik, mer än att den blir lite mer precis, utan något som vi ägnar oss åt jämt och ständigt, undermedvetet eller ej. Det kanske handlar om en mental modell över gruppdynamiken i en klass, vilken i sin tur kan påverka vårt beteende eller hur ska jag kryssa mellan bilarna på parkeringsplatsen utanför Maxi för att minimera avståndet till dörren. Beroende på väder väljer vi kanske helt olika vägar när vi ska förflytta oss inom Högskoleområdet, osv. Vi sammanfattar redan nu några punkter om Matematisk modellering Handlar om att kunna omsätta matematiska kunskaper i praktiken för att erhålla användbara lösningar på problem hämtade från verkligheten. Man måste bland annat kunna avgöra vilken och hur en matematisk teori kan användas i ett specifikt problem. Ett bra sätt att systematiskt struktuera och analysera ett verkligt problem. Konstruktion av modellen leder ofta till ökad kunskap om det verkliga problemet. Matematisk analys av modellen kan ge ökad insikt om det verkliga problemets egenskaper. Med modellen kan vi simulera, eller härma verkligheten under hypotetiska scenarier, t.ex. utföra experiment som inte låter sig göras i verkligheten på grund av att det vore för dyrt, tidsödande eller riskabelt, eller helt enkelt omöjligt. Är ett tvärvetenskapligt ämne där man använder matematiken som ett verktyg i tillämpningar, som inte alltid vid första anblicken ter sig matematiska. Är ingen exakt vetenskap med entydiga korrekta svar. Olika angreppssätt kan ge olika resultat. Lär man sig enklast (endast?) genom att själv praktisera det. Matematisk modellering är en mycket kreativ sysselsättning! ť Vad är en matematisk modell? Begreppet matematisk modell har blivit ett populärt inslag i den tillämpade matematiken, ungefär från 1960-talet och framåt. Modelltänkandet erövrar alltfler områden där matematiken används: naturvetenskap, teknik, ekonomi, samhällsvetenskap osv. En orsak till den här utvecklingen är de enorma möjligheter som numera finns att med datorers hjälp göra snabba numeriska beräkningar. Det går att arbeta med komplicerade matematiska konstruktioner och ända få fram användbara resultat. Före datorernas tid var man hänvisad till analytiskt lösbara formuleringar eller tidsödande numeriska beräkningar. ť Modellbegreppet Ordet modell kommer från det latinska ordet modulus (litet mått) och betyder förebild eller mönster. En modell avbildar eller föreställer något annat. Detta som avbildas kallar vi det verkliga eller riktiga systemet. Den fysiska modellen är ett eller flera objekt som konstrueras för att efterlikna ett system. Modelljärnvägen avbildar "riktiga" tågvagnar, lok, spår, växlar, signalsystem osv. Ibland görs modellen före det riktiga systemet. En fartygsmodell där skrovets form är noga definierad byggs traditionellt före fartyget och används vid byggandet för att ge varje del dess rätta form. Vid modellförsök använder man skalriktiga modeller för att undersöka fysikaliska förhållanden. HH/ITE/BN Matematisk modellering och Mathematica 3 När en fysisk modell skapas nöjer man sig med att avbilda vissa egenskaper hos systemet och bortse från andra. Vad som är intressant eller ointressant avgörs av syftet med modellen. Studerar man luftströmningen runt en flygplanskropp behöver man inte montera in stolar i flygplansmodellen! När en modell av en damm konstrueras för att studera belastningar använder man andra material än de som den riktiga dammen senare byggs av. Den store filosofen och matematiken Rene Descartes (1596–1650. Dog i Stockholm som drottning Kristinas personlige lärare i matematik.) introducerade och utvecklade begreppet matematisk modell. I början hade Descartes en dröm om att hitta en universell metod, en matematisk metod förstås, för att beskriva varje problem som kunde komma från verkligheten och sedan lösa problemet. Även om hans dröm aldrig förverkligats, så finns det mängder av problem som kan lösas i denna anda. I korthet gick den ut på följande: Reducera problemet till ett matematiskt problem. Reducera det matematiska problemet till ett algebraiskt problem. Reducera det algebraiska problemet till att lösa en ekvation. ť Matematisk modell En matematisk modell "avbildar" eller beskriver ett system med hjälp av matematiska begrepp och storheter. Exempel: är lätta att hitta genom historien Människor har i alla tider haft en önskan att förstå dygnets växling mellan dag och natt, årstidernas regelbundna återkomst, sol- och månförmörkelser osv. Många matematiska modeller för solsystemet har utvecklats under historiens gång. På medeltiden föreställde man sig jorden som fast och orörlig. Runt jorden rörde sig solen och planeterna på bestämda sfärer (den "ptolemeiska" världsbilden). Den Kopernikanska modellen placerar i stället solen i medelpunkten med jorden och de övriga planeterna i cirkulära banor runt den. Kepler var den som kom fram till att planetbanorna beskriver ellipser. Dagens astronomer arbetar med förfinade modeller där nya planeter, kometer och asteroider ingår. Modellerna av solsystemet är formulerade i geometriska eller matematiskt analytiska termer och är därför matematiska modeller. Den allmänna gravitationslagen (formulerad av Newton 1684) kallas en naturlag. Den ingår i den klassiska matematiska modellen för partiklars och kroppars dynamik. Modellen har visat sig vara enormt användbar för att beskriva kroppars rörelse i verkligheten. Men liksom andra modeller har den ett begränsat giltighetsområde. Den kan t.ex. inte användas för att beskriva vad som händer inuti atomer. Hormonet insulin balanserar sockerkoncentrationen i blodet. Man kan matematiskt beskriva hur halterna av insulin och socker i blodet påverkar varandra. ť Varför gör man modeller? För att få svar på frågor om hur en system fungerar kan man genomföra experiment med systemet. Men det finns många tillfällen då ett experiment med det riktiga systemet inte är möjligt. Skälen kan vara olika: det kan vara omöjligt: Det går inte att experimentera med solsystemet eller med väder och vind. Inte heller kan man experimentera med en system som ännu inte finns utan ska konstrueras. det kan bli för dyrt: Att experimentera med reglersystemet i en stor processindustri kan kosta hundratusentals kronor i förstörd produktion. det kan vara farligt. Provflygningar med nykonstruerade flygplan skulle vara mycket riskabla utan föregående noggranna beräkningar i matematiska modeller och studier av fysiska modeller av planet. Att testa fram dosering av mediciner på människor utan att utsätta personer för risker är bara möjligt inom snäva gränser. att få bättre förståelse av systemet. att sammanfatta teorier om systemet. att förmedla generell kunskap om systemet. att strukturera eller formalisera kunskap och tankar för diskussion och kritik. att kvantifiera skeenden i mer komplexa sammanhang. att studera ett planerat eller hypotetiskt system när ett verkligt inte finns. att göra projektioner för framtiden utifrån befintlig kunskap. att simulera experiment som är för svåra, farliga, tidskrävande eller dyra i verkligheten för att hinna med i konkurrensen. att använda som redskap i ett prognos- eller förvaltningssammanhang. Vi sammanfattar: Matematisk modellering är en nödvändighet inom all modern utveckling! I en konkurrensutsatt verksamhet gäller det för en modern ingenjör att minska ledtider och spara in på kostsamma provserier. 4 Matematisk modellering och Mathematica HH/ITE/BN ť Modellen ger ny kunskap I många fall är man alltså hänvisad till studier av modeller av systemet. Hur kommer det sig då att man kan få reda på något nytt något som man inte redan visste när modellen konstruerades - genom att undersöka en modell? Den matematiska modellen är uppbyggd på ett antal grundläggande förutsättningar som väljs så att de beskriver centrala egenskaper hos systemet. Egenskaperna är experimentellt eller erfarenhetsmässigt verifierade. Sådana egenskaper kan t.ex. vara naturlagar (som i själva verket också de är matematiska modeller som är väl undersökta och vars giltighet man noga känner sedan tidigare). Det kan också handla om egenskaper man har iakttagit hos det speciella systemet ifråga. Utifrån dessa egenskaper, som formuleras i matematiska termer, och andra antaganden eller förenklingar man kan behöva göra, drar man slutsatser. Då hanterar man sina matematiska objekt och använder logikens lagar enligt vad som gäller inom matematiken. De matematiska teorierna är i sin tur baserade på vissa axiom och logiska regler. Det betyder att de slutsatser man kommer fram till är trovärdiga logiska konsekvenser av de förutsättningar man grundat modellen på. Det gäller åtminstone om de förenklingar man gjort på vägen inte varit för grova. Den nya kunskap om systemet som modellen ger oss är alltså i princip logiska konsekvenser av egenskaper som vi själva valt ut och som vi anser karaktärisera systemet. De konsekvenserna kan emellertid vara svåra att lista ut utan hjälp av matematiken. Matematikens formalism är bekväm och många teorier grundligt genomarbetade. Därför blir matematiken ett mycket kraftfullt verktyg. ť Vilka frågor kan modellen besvara? De frågor man ställer kan röra sig om vitt skilda saker. Exempel: på frågor vars svar man kan söka genom att konstruera och studera matematiska modeller: Varför skakar min bil just vid vissa hastigheter och inte andra? Vilket jakttryck tål en viss älgstam? Hur hög ska en skorsten byggas för att röken inte ska smutsa ned luften i närheten av anläggningen? Vilket väder får vi i morgon? Hur ska insulindosen väljas för en viss patient? Hur ska jag optimera ett visst transportsystem med tanke på energiåtgång, tidskrav, ekonomi? I hur många steg ska en raket byggas? Är vi på väg mot ett varmare klimat på jorden på grund av de koldioxidutsläpp människan orsakat eller är vi på väg mot en ny istid? Hur ska ett effektivt datorprogram för lösning av linjära ekvationssystem se ut? Kommer sälstammen utanför svenska västkusten att kunna återhämta sig? Vilken effekt på skatteplanerandet får förslaget till ny aktiebeskattning? Typen av frågor går att klassificera. Ett sätt att göra det är följande: Förståelse. Vi söker förstå hur ett system fungerar - kanske av ren nyfikenhet. Förutsägelse. Vi önskar kunna förutsäga hur systemet kommer att bete sig i framtiden. Styrning. Vi har eller vill skaffa oss metoder att påverka systemet så att det fungerar på ett förutbestämt sätt. Konstruktion. Vi är ute efter att konstruera ett visst system vars egenskaper vi önskar bestämma i förväg eller som vi vill göra optimala i någon mening. De fyra typerna bildar nivåer i en hierarki. Förståelse är grundläggande och krävs för att kunna göra förutsägelser. Både kontroll och konstruktion bygger på de två andra. För en ingenjör är det naturligt att arbeta med de två sista typerna av frågor. Grundforskning inom naturvetenskapen handlar ofta om den första typen medan tillämpad teknikforskning mest rör sig med frågor av de andra tre typerna. ť Klassificering av matematiska modeller En modell kallas mekanistisk om den beskriver ett orsaks-samband. En icke-mekanistisk modell kan bestå av empiriskt baserade samband mellan parametrar och variabler, i det enklaste fallet en kurvanpassning. Ofta får man börja att arbeta med en icke-mekanistisk modell för att senare kunna sluta sig till hur en mekanistisk ("förklarande") modell kan se ut. En modell är statisk om den beskrivs helt och hållet av de momentana (just aktuella) värdena av vissa parametrar eller insignaler. Modellen har inget "minne". Ett dynamiskt systems tillstånd däremot bestäms också av sin förhistoria. I den matematiska beskrivningen kommer tidsderivator in när det gäller dynamiska modeller. Ett dynamiskt system brukar beteckna ett system som kan beskrivs i en matematisk modell där ekvationerna som ingår är ordinära differential- eller differensekvationer. Man skiljer också på kontinuerliga och diskreta modeller. Diskreta modeller innehåller variabler som förändras i diskreta steg (t.ex. antal) eller har bara ändligt många alternativa tillstånd. Kontinuerliga modeller beskrivs av variabler som förändras kontinuerligt eller infinitesimalt. Ett kontinuerligt system modelleras med differentialekvationer. HH/ITE/BN Matematisk modellering och Mathematica 5 En deterministisk modell ger åt varje ingående variabel ett väldefinierat värde. I en stokastisk modell arbetar man med slumpmässighet. En variabel har inte ett bestämt värde utan beskrivs av en sannolikhetsfördelning. När ett system innehåller inslag av slumpmässighet kan det vara nödvändigt att göra en stokastisk modell för att den ska bli realistisk. En komplikation är att deterministiska system kan ha ett stokastiskt utseende. Det gäller framför allt så kallade kaotiska system. Den matematiska definitionen på kaotiskt dynamiska system är att mycket små skillnader i begynnelsevärden medför stora skillnader i slutändan. Exempel: För att bestämma nedböjning och påkänningar på en bro används ofta Eulers balkteori som modell. För samma indata, brospann, pilonhöjd, balktjocklekar osv. får man samma svar varje gång man räknar. Vi har en deterministisk modell. Hur ska man dimensionera antalet öppna kassor på ICA? Med kännedom om hur många personer det kommer per tidsenhet och deras överläggningar angående kön till kassan kan man med slumptal studera hur vinsten blir vid olika konfigurationer. Här får vi olika svar varje gång vi kör modellen, alltså en stokastisk modell. Vädret är ett mycket bra exempel på svårförutsägbart system. Lite skämtsamt brukar man säga att en fjäril som fladdrar i Amazonas kan påverka vädret i Sverige. Riktigt så illa är det nog inte, men principen är riktig. Vi har att göra med ett kaotiskt system. ť Hur gör man en matematisk modell? Att modellera är något man lär sig genom erfarenhet och det är inte lätt att beskriva en metod som fungerar i alla situationer. Ett krav är kännedom om det aktuella systemet. Ett annat är goda matematiska kunskaper. Det går dock att peka på vissa aspekter av modelleringens konst redan med enkla modeller som inte kräver specialistkunskap eller avancerade matematiska metoder. Ofta får man ta fram en första modell som sedan förbättras i flera omgångar. Arbetet kan inte formaliseras helt och hållet men en viss struktur kan urskiljas. I arbetet med en viss modell följer stegen inte alltid efter varandra i den logiska ordning som presenteras här men i princip kan man säga att man vevar runt tills man är nöjd! Verklig värld Länk mellan världarna Matematisk värld Specificera det Gör nödvändiga verkliga problemet antaganden Formulera det matematiska problemet Utvärdera modellen Tolka den matematiska lösningen Lös det matematiska problemet Använd modellen för att kommunicera, förklara, förutsäga, bestämma, planera, börja, sluta… Man ska vara medveten om att alla system uppfattas och "filtreras" genom de sinnen och instrument som vi har till hands när vi ska uttala oss om "sanningen". Vår uppfattning om verkligheten är alltså inte verkligheten utan redan det en tanke-modell. Dessutom är alla system delar av andra system i en hierarki. Det gäller att göra lämpliga avgränsningar. Exempel: Antag att du vill konstruera en modell av hur mängden grodor varierar i en damm. I det fallet kan man anta, med viss säkerhet, att cellprocesser är så snabba att de inte nämnvärt påverkar variationen i grodantal. Man kan också anta att mängden grodor i en damm flera mil därifrån inte heller nämnvärt påverkar den lokala populationen. Vidare kan man nog anta att eventuella evolutionära processer är så långsamma att de inte heller påverkar antalet. På så sätt kan man begränsa detaljrikedomen i modellen, och göra den hanterlig. Hade däremot syftet med modellen varit att beskriva grodornas näringsupptag hade det varit nödvändigt att avgränsa detaljrikedomen på ett annat sätt. Beteendet hos modellen beror helt och hållet på vad man stoppar in, det vill säga de antaganden man gör om systemet. Ändrar man ett antagande kan man ändra utfallet totalt. Därför är det viktigt att noggrannt ange alla antaganden för att utomstående ska kunna bilda sig en uppfattning om var och när en modell är giltig. Om antagandena inte är giltiga för ett visst system kan modellen ge en felaktig bild av den verklighet den avser att avbilda. För att utvärdera modellens giltigheten är det därför nödvändigt att jämföra med verkligheten. Var också observant på orimliga resultat från modellen. Det går att konstruera modeller för vilka system som helst, fast metoden kan variera kraftigt beroende på frågeställning och mellan vetenskapliga discipliner. I vissa discipliner är verbala eller logiska modeller vanligast medan andra discipliner använder nästan enbart matematiska modeller. I princip är det ingen skillnad, men den matematiska modellen har fördelen att vara betydligt mer koncis och transparent. Oftast är avnämaren inte någon matematiker därför gäller det att föra samtalet kring problemställningen i en god resonerande ton och behålla de matematiska reflektionerna som dyker upp för sig själv. 6 Matematisk modellering och Mathematica HH/ITE/BN ť Precisera frågeställningen Den ursprungliga frågan är ofta vag och oprecis. Gör man modellen på uppdrag av någon annan är det särskilt viktigt att ta god tid på sig att precisera frågan. Det kan också bli aktuellt att föreslå en annorlunda infallsvinkel som leder till en helt ny frågeställning. Det är ofta svårt att precisera frågeställningen utan att samtidigt skaffa sig kunskaper om hur systemet fungerar. Exempel: Antag att vi får i uppdrag att göra ett förslag till belysning av en idrottsarena. Vilket är bästa sättet att ordna belysningen? Vad menas med "bästa"? Starkast, jämnast, utan reflexer ... ? Eller menar man kanske billigast? Vilka minimikrav ställs? Ska vi ta fram alternativa modeller till olika kostnad och olika kvalitet? Vill vi ha en kvalitativ modell, det vill säga en modell som reagerar proportionellt på ingående parametrar på samma sätt som den verkliga modellen eller krävs en kvantitativ modell där vi även kräver numerisk relevans? Exempel: Antag att vi ska dimensionera hjulupphängningen på en bil och har tagit fram en modell över hur någon storhet varierar med mönsterdjupet Δ på däcket. Vi vill nu se vad som händer när vi dubblerar vikten på hjulet. Vi förväntar oss ett högre värde på vår studerade storhet. Till vänster har vi en icke kvalitativ modell. Att den blev rätt i första fallet var förmodligen tur eller "fusk". I mitten har vi en kvalitativt riktig modell. En sådan är ofta enkel, snabb och räcker långt för att bedriva utvecklingsarbete över en kopp kaffe på Heathrow. Den högra är kvantitativ och därmed även kvalitativ. En sådan är oftast mer omfattande, kanske svåröverskådlig, och kräver datorberäkning. Ej kvalitativ Kvantitativ Kvalitativ prov, modell 5 prov, modell 5 prov, modell 5 4 4 4 3 3 3 2 2 2 1 1 0 1 2 3 4 5 Δ 1 0 1 2 3 4 5 Δ 0 1 2 3 4 5 Δ Det är mycket viktigt att formulera sig i kvantifierbara termer eftersom modellens relevans så småningom ska utvärderas mot vad vi preciserat här! ť Identifiera faktorer och samband De faktorer i systemet som kan påverka svaret på vår frågeställning kan vara många. Alla sådana faktorer kandiderar till en plats i den matematiska modellen. Man kan dela upp arbetet med att hitta faktorer i två steg. I en "brain-storming"-fas listar man alla tänkbara faktorer. I nästa fas väljer man ut de faktorer som ska finnas med i modellen. Anledningen att utesluta en viss faktor kan vara någon av följande: Dess inverkan är försumbar jämfört med andra faktorer. Vi kan eventuellt ta med den i en förfinad version av modellen senare. Det är omöjligt att förutse dess konsekvenser. Vi får tänka oss att den faktorn inte varierar eller spelar en liten roll. Det är svårt eller invecklat att beskriva hur faktorn påverkar systemet. Det gäller en slumpmässig variation. Vi är av något skäl inte beredda att göra en stokastisk modell eller tror att slumpmässigheten inte påverkar resultatet nämnvärt. Vi får arbeta med medelvärden i stället. Det är lätt att inse att vi kan införa stora begränsningar när vi bestämmer vilka faktorer som inte kommer med i modellen. Det gäller att hålla detta i minnet när utvärderingen av modellen sker så småningom. Är modellens förutsägelser dåliga kan det bero på att vi har försummat en faktor som inte var oviktig. Hur påverkar de olika faktorerna varandra? Om det går att formulera samband bör det göras redan nu. Än så länge kan man nöja sig med kvalitativa beskrivningar. Vrid och vänd på problemet. Förenkla! Idealisera! Vilka storheter kan ha betydelse? Använd kunskap, förnuft, intuition eller chansa! Vad vet jag? Vad söker jag? Vilken informationen finns? Nödvändiga? Onödiga? Motsägande? Verkar det rimligt? Vad har jag sett för teorier om sådant här? Kanske kan man omformulera problemet? Har jag sett något liknande tidigare? Vad skiljer i så fall problemen åt? Vad behöver jag ändra i den för den ska passa in här? Lista, rita bubblor och pilar och försök länka samman! Kanske kan jag hämta en modell från ett helt annat tillämpningsområde. Kom ihåg att en roll matematiken har är att koka ner olika tillämpningar och peka på gemensamma strukturer i grunden som absolut inte verkar uppenbara för det otränade ögat! Exempel: Så vitt skilda fenomen som radioktivt sönderfall, bakterietillväxt, avsvalning av sockerkaka, blandning av föroreingar i en sjö, m.fl. beskrivs av exakt samma differentialekvation. Det är bara dimensionen på ingående storheter och parametrar som varierar. Tänk på att börja enkelt! Om vi studerar modelleringsloopen ovan ser vi att det finns bara en väg ut! Det duger alltså inte att göra något så komplicerat att vi inte kommer runt! Sitt inte i ett hörn och deppa. Framåt till "varje" pris...!! I ett företag duger det inte att HH/ITE/BN Matematisk modellering och Mathematica 7 akademiskt proklamera "Ett intressant matematiskt problem, det löser vi kanske om 100 år!". Då har förmodligen företaget gått i konkurs för länge sedan. Så Einsteins gamla devis duger gott "Make it as simple as possible, but not simpler!" Börja enkelt! Det kanske räcker! Annars modifierar vi i nästa varv av modelleringsloopen ovan! ť Översätt till matematiskt språk De faktorer som ska inga måste kvantifieras för att kunna användas i en matematisk modell. Vi inför alltså ett antal variabler som beskriver faktorer som vi vill studera. Vi kan också introducera vissa hjälpvariabler eller parametrar som gör det möjligt att beskriva egenskaper hos systemet. Parametrar är storheter som vi väljer att hålla konstant under en betraktelse. Ofta gör man en parameterstudie och parametrar och variabler kan i en annan modell byta plats med varann. Var noga med att tänka igenom vilken dimension olika variabler har. Gör dimensionsanalys på alla samband! Den preciserade frågeställningen ska översättas till matematiskt språk. Den kommer att innehålla frågor av typen: Vilket är det maximala värdet av...? Vad blir storleken av en viss variabel vid en viss tidpunkt? Hur lång tid tar det innan funktionens värde blir mindre än...? För stora system med många variabler blir frågeställningen sammansatt av många delfrågor. Samband vi tidigare har funnit ska också översättas till matematiskt språk. För att det ska vara möjligt måste sambanden kvantifieras. Ibland innebär kvantifieringen en stark förenkling. Att studera extremförhållande ger ofta ledtrådar i modellval. Välj så att uppenbara egenskaper i preoblemställningen uppfylls, exempelvis; När tiden går mot oändligheten går koncentrationen mot noll! När massan går mot noll försvinner kraften! Andra samband som nu formuleras är relationer mellan (hjälp)variablerna som kan beskriva geometriska villkor, kontinuitetssamband (massbalans t ex) osv. De samband man sätter upp har ofta formen av ekvationer. Antingen är det vanliga algebraiska ekvationer eller differential- eller differensekvationer. Inte sällan har man ett optimeringsproblem. Ett samband kan vara ett tidigare känt samband från t.ex. fysik. Men om det inte är välkänt så ska man göra noga klart exakt vad sambandet innebär. Hela analysen och resultaten bygger nämligen på dessa samband och de måste vara rimliga, annars blir modellen inte bra. Exempel: Att modellering ger upphov till ekvationer är ganska naturligt eftersom en modell av verkligheten ofta representeras av något slags samband, exempelvis volymen av en konservburk V Πr2 h. Det är ju sällan man väljer r och h och ser vad volymen blir. Det vanliga scenariot i verkligheten är snarare att en ingenjör ställs inför mängder med krav på sin modell, exempelvis "Vilken är den minsta plåtåtgång för en given volym?". Detta i sin tur leder till krav på modellens variabler, i vårt fall r och h. Om flera krav skall beaktas och jämkas samman leder detta sedan naturligt vidare till optimering. Läsa spelet En av svårigheterna är att kunna dechiffrera en verbal beskrivning till en matematisk beskrivning. Vi exemplifierar med en mycket vanlig situation, nämligen någon storhet, säg y, som varierar med tiden t. Då har vi en differentialekvation att vänta och vi söker funktionen y t . Som hjälp vid dechiffreringen är det bra att känna igen vanliga formuleringar som skall tolkas som derivata av den sökta funktionen y Ändring av y per tidsenhet t y Ändringshastigheten av y t eller y' t . eller y ' t . I en uppsjö av tillämpningar är olika typer av proportionalitet inblandade. Vissa ord kan också direkt associeras med matematiska krumelurer, exempelvis y t är proportionell mot t yt kt. k y t är omvänt proportionell mot t yt k . t k Ändring av y t per tidsenhet är proportionell mot y t y' t y' t Ändring av y t per tidsenhet är proportionell mot y t och t y' t y' t k y t t. k Ändring av y t per tidsenhet är proportionell mot y t och skillnaden mellan m och y t y' t ky t . k k m yt y' t ky t m yt . 8 Matematisk modellering och Mathematica HH/ITE/BN För fler exempel och fler tips inom specifika tillämpningar hänvisas till senare avsnittet i "Något om ..." serien längre fram i kursblocket Tillämpad Matematik. Fysikaliska principer Ofta kommer någon fysikalisk lag eller princip till användning. 2 Newtons accelerationslag m y t2 F my F y KR y y F. Energins bevarande Massans oförstörbarhet. massa = densitet (eller koncentration) gånger volym. Arkimedes princip: Då en kropp nedsänkes i en vätska påverkas denna av en lyftkraft som är lika stor som den undan– trängda vätskans tyngd. Dimensionsanalys Se till att det är lika många "äpplen och päron" på båda sidor om – tecknet i ekvationerna! Dimensionsanalys är ett stort och viktigt stöd vid modellering och vid utvärderandet av modellen. Ta för vana att utnyttja denna hjälp under hela arbetet! Ett studium av de storheter som ingår i problemet ger nästan alltid direkta tips som leder till målet. Se vidare i Något om Dimensionsanalys och Mathematica . ť Analysera modellen Nu vidtar det vanliga matematiska arbetet. Det gäller att lösa ekvationerna. Man kan använda analytiska eller numeriska metoder. Båda typerna av metoder är lika viktiga. Innan man tar till numeriska metoder bör man analysera modellen så långt det är möjligt. Att söka numeriska lösningar betyder ofta att man är tvungen att välja värden på vissa variabler eller parametrar. Det betyder att resultaten inte blir lika allmängiltiga som vid en analytisk lösning. Man måste noga fundera på vilka variabelvärden det är lönt att räkna på. Frågeställningen i modellen bestämmer valet av numerisk strategi. När det gäller stokastiska modeller kan man genomföra simuleringar. Det innebär att man med hjälp av en slumptalsgenerator i modellen efterliknar naturens slumpmässiga beteende. Simuleringar genomförs alltid i dator. Begreppet simulering används ibland i en annan betydelse. Då har det inte med slumpmässighet att göra. Simulering får då helt enkelt betyda numerisk lösning av differentialekvationerna i ett kontinuerligt eller dynamiskt system. Man "simulerar" sitt system i den meningen att man efterliknar det i datorn. Man får då ofta nöja sig med att dra slutsatser utifrån mängder med simuleringar. När man lämnat skolan och dess uppriggade ekvationer för att passa handräkning möter man som praktiskt arbetande ingenjör mängder med ekvationer som inte går att lösa exakt, oftast beroende på olinjäritet. Frågorna som dyker upp är "Finns det någon lösning?", "I så fall hur många?", "Hur identifierar jag den eller de som är relevanta?", "Hur hittar jag den som löser mitt problem?". Oftast har en ekvation flera lösningar, varav en del är mer exotiska än andra, "negativa längder", "komplexa massor" och så vidare. Dessa har kanske inget med verkligheten att göra men likväl är modellen korrekt! Man brukar säga att en matematisk modell är rikare än den fysikaliska modell den beskriver. Sådana här frågeställningar kommer ofta i ny dager om de belyses med grafik under arbetets gång! Naturligtvis används moderna och effektiva datorprogram för analysen och visualisering av resultat. Tillägna dig ett arbetssätt där du så ofta som möjligt använder grafiska representationer! ť Tolka resultatet av analysen I nästa fas går vi tillbaka till systemet och tolkar resultaten av den matematiska analysen i det icke-matematiska språk som problemet ursprungligen formulerades i. Den här fasen är intressant och kan ibland bjuda på spännande överraskningar. Det kan vara så att systemet ger utfall som är oväntade. I så fall har vi fått reda på något nytt om vårt system genom att göra en modell för det. En mindre spännande möjlighet är att det finns någon felaktighet i den matematiska analysen, vi kanske har räknat fel. Slutligen kan modellen vara för grov eller uppenbart felaktig när man ser resultatet. Exempel: Man vill av tunn plåt tillverka en cylindrisk konservburk med given volym V . Bestäm radie och höjd i den burk som kräver minst materialåtgång, det vill säga har minst total area. HH/ITE/BN Matematisk modellering och Mathematica 9 Lösningsförslag: Antag att konservburken har höjden h och radien r. Dessa kan nu inte variera fritt oberoende av varandra, de binds samman av att volymen på burken är given V Π r2 h. Sådana här kopplingar brukar kallas för just kopplingsvillkor. Totala arean av burken byggs upp av två lock, 2Al 2 Π r2 , samt mantelarean Am omkrets höjd 2 Π r h. Gör nu inte för mycket för hand, varje sådan insats är en potentiell risk för att introducera fel. Låt Mathematica göra jobbet! Skriv bara ned alla grundsamband. Glöm inte dimensionsanalys! ekv Π r2 h, Atot V Π h r2 , Atot V 2 Al 2 Al Π r2 , Am Am , Al Π r2 , A m A m , Al 2 Π r h 2 Π h r Utnyttja att V är given för att lösa ut Atot som funktion av r. Ta för vana att lösa ut lika många variabler som vi har ekvationer. Även de variabler som inte primärt används vid optimeringen är oftast intressanta att veta värdena på till slut. Så alla som funktioner av r! Amm Solve ekv, Atot , h, Al , Am 2 Π r3 Atot V V ,h r Πr 2V Π r2 , Am , Al 2 First r Innan man börjar på allvar är det utmärkt att pigga upp sig med en bild över situationen för att se om modellen är sund. Visst, tydligt minimum som sig bör, eftersom Atot både då r 0 och r . PlotAtot V2 3 x V1 3 , x, 0.1, 1 , PlotStyle . Amm . r PlotRange 5, 10 , AxesLabel "r V 1 3 ", "Atot V 2 3 Red, " Atot V 2 3 10 9 8 7 6 0.0 0.2 0.4 0.6 0.8 1.0 r V1 3 Bestäm nu det r som minimerar Atot genom att söka nollställe till derivatan, dAdr D Atot 2 Π r3 6Πr r 0. . Amm, r V r2 Solve dAdr 0, r 3 1 r Atot r 3 3 V 23 1 V , r 3 V , r 2Π 3 3 2Π 2Π N r r 0.270963 3 0.469322 V , r 0.541926 3 V , r 0.270963 0.469322 3 V Här är det bara den mittersta lösningen som är reell, de andra två komplexa har inte med saken att göra. Vi har uppenbart minimum, varav slutligen alla variabler vid detta välsignade tillstånd. Amm . r 2 Atot 3 3 2Π V 23 3 22 3 ,h 3 3 V , Al Π Π V2 3 22 3 , Am 2 3 2 Π V 2 3 Alla symboliska resultat måste underkastas dimensionsanalys!! Här har vi för arean Atot höjden h 3 3 22 3 3 3 Π 2 Π V 2 3 1 V 2 3 m3 23 m2 , Ok! 3 V 1 V 1 m3 1 3 m , Ok! Med symboliska svar kan man lätt få en kvalitativ bild av hur modellen påverkas av olika storheter. 10 Matematisk modellering och Mathematica HH/ITE/BN ť Utvärdera modellen Den systematiska utvärderingen är viktig. Utvärderingen sker först och främst mot det punkter som formulerades när frågeställningen mejslades ut. Det är därför så viktigt att detta görs ordentligt i inledningen till modelleringen. Har man möjlighet att jämföra med resultatet av experiment eller observationer av det riktiga systemet eller ett liknande system (kanske en nedskalad modell) är det bra. De experiment man behöver göra för att testa modellen är färre, enklare och billigare än om man hade undersökt systemet och försökt besvara frågan utan modellens hjälp. Har man inte möjlighet att göra experiment får man utgå från de data man har tillgång till. Det bästa är givetvis att testa med andra data än dem man utgick ifrån för att välja sambanden i modellen. Att testa extremvärde är ofta fruktbart! Modellens giltighetsområde ska undersökas. Alla modeller har någon form av begränsning i sin giltighet. Exempelvis är det helt olika modeller om man ska göra en väderprognos över ett dygn eller tio dygn! Så länge variablerna ligger i intervall där vi har empiriska data rör man sig på mycket säkrare mark än när man försöker extrapolera. Det kan tänkas att utvärderingen leder till slutsatsen att modellen har brister. Det gäller då att inte vara alltför fast vid sin skapelse utan acceptera att modellen behöver ändras. Den "misslyckade" modellen ger information om att någon viktig faktor försummats. De kunskaper om systemet man skaffat genom att arbeta med den första modellen har man nytta av i fortsättningen. I själva verket är det paradoxalt nog själva modelleringsprocessen som är intressantast och ger de nya kunskaperna, inte den färdiga modellen. En felkalkyl, noggrannhetsanalys och känslighetssanalys bör göras. De empiriska data som modellen utgår från har en viss felmarginal som påverkar resultatet av analysen. En modell som är väldigt känslig för variationer i variabler och parametrar är förmodligen inte riktigt sund. En numerisk analys inför också nya felkällor. Det gäller förstås också en stokastisk simulering. Hur säkra kan vi vara på vart resultat? Kom ihåg att en modell får inte användas utanför de förutsättningar som bestämdes då problemet specifiserades! ť Förenklingens roll Det är ofta givande att börja med en starkt förenklad modell. Den kan vara lätt att konstruera och snabb att räkna igenom. Den kan visa hur mycket av systemets egenskaper som förklaras av en eller ett par centrala faktorer och den kan slutligen användas som utgångspunkt för nästa bättre modell. Detta kan vara en god hjälp om man råkar ut för det som är så vanligt i problemlösning, man "sitter fast". Ett annat "knep" är att testa vad modellen säger om ytterlighetsfallen. Ofta vet man hur systemet bär sig åt i de fallen (tiden går mot oändligheten, massan går mot noll osv.). Det ger en möjlighet att utvärdera modellen under arbetets gång. ť Kommunicera modellen Om utvärderingen av modellen har varit nöjaktig är det läge att sätta den i arbete! Annars vore det ju ingen mening med att göra den. Så använd den så mycket som möjligt till att kommunicera, förklara, förutsäga, bestämma, planera, börja, sluta...så! Skaffa er fördelar gentemot andra som håller på med liknande produkter men ännu inte upptäckt Matematisk modellering!
© Copyright 2024