Simpsons regelkalkylator
Beräkna kurvlängd med Simpsons regel när exakt integration är opraktisk eller onödig.
Ursprunglig kurva och vald båge
Guldpunkter markerar mittpunkten i varje delintervall.
Integrand g(x) och mittpunktsrektanglar
Hovra över en rektangel för att se dess delintervall, mittpunkt, integrandvärde och bidrag.
Detaljerad beräkningstabell för mittpunktsregeln
Konvergensjämförelse
Skillnader mellan på varandra följande förfiningar är en praktisk felindikator, inte en garanterad felgräns.
Beräkningsfel
Sektionsnavigering
Vad är Simpsons regel?
Simpsons regel är en metod för numerisk integration som approximerar värdet av en bestämd integral med hjälp av andragradspolynom. I stället för att uppskatta arean under en kurva med rektanglar eller trapetser anpassas en parabel genom varje grupp om tre på varandra följande stödpunkter, och den exakta arean under varje parabel beräknas. Summan av dessa parabelareor ger approximationen.
Metoden är uppkallad efter Thomas Simpson, en engelsk matematiker från 1700-talet, även om grundidén förekom tidigare i andras arbete. Dess styrka kommer från en enkel observation: en parabel kan följa en jämn kurva mycket närmare än vad en rät linje kan. När funktionen som integreras är jämn ligger de paraboliska bitarna tätt mot den sanna kurvan, och den resulterande approximationen blir anmärkningsvärt noggrann även med ett måttligt antal delintervall.
För att förstå varför tre punkter spelar roll, minns att en rät linje bestäms fullständigt av två punkter, medan en parabel (ett andragradspolynom) behöver tre. Simpsons regel utnyttjar detta genom att dela integrationsintervallet i ett jämnt antal delintervall, gruppera dem i par och anpassa en parabel över varje par. Varje par spänner över två delintervall och tre punkter: vänster ändpunkt, mittpunkten och höger ändpunkt.
Den klassiska formen kallas Simpsons 1/3-regel eftersom faktorn 1/3 förekommer i formeln. Det finns också en Simpsons 3/8-regel som använder tredjegradspolynom över grupper om tre delintervall, men 1/3-regeln är standardmetoden och den version som implementeras i de flesta kalkylatorer, inklusive den på den här sidan.
Simpsons regel tillhör familjen Newton-Cotes formler, som approximerar integraler genom att utvärdera integranden i jämnt fördelade punkter och kombinera dessa värden med fasta vikter. Det som skiljer Simpsons regel är dess viktmönster: 1, 4, 2, 4, 2, ..., 4, 1. Ändpunkterna får vikt 1, inre punkter på udda positioner får vikt 4 och inre punkter på jämna positioner får vikt 2. Detta alternerande mönster är metodens kännetecken och källan till dess noggrannhet.
Förstå båglängd i matematik
Längden av ett rakt linjesegment är enkel: mät avståndet mellan dess ändpunkter. En kurva är annorlunda. Den böjer och vrider sig, så ingen enskild rak mätning fångar dess sanna utsträckning. För att hitta längden av en kurva bryter matematiker den i små bitar, approximerar varje bit som rak och summerar resultaten.
Föreställ dig att du går längs en slingrande stig. Om du tar mycket små steg är varje steg nästan rakt, även om stigen som helhet kurvar. Summera alla de små raka stegen så får du en bra uppskattning av den totala tillryggalagda sträckan. Gör stegen mindre och mindre, och uppskattningen konvergerar mot den exakta längden. Detta är kärnidén bakom båglängd.
För en funktion y = f(x) på ett intervall [a, b], betrakta ett litet horisontellt steg dx. Över det steget stiger funktionen med ungefär dy = f'(x) dx, där f'(x) är derivatan (lutningen). Det lilla kurvsegmentet har längden som ges av Pythagoras sats: ds = √((dx)2 + (dy)2) = √(1 + (f'(x))2) dx. Att summera alla dessa infinitesimala bitar genom integration ger båglängdsformeln:
L = ∫ab √(1 + (f'(x))2) dx
Varje del av denna formel har en tydlig betydelse. Symbolen ∫ab säger att vi ska ackumulera bidrag från x = a till x = b. Uttrycket f'(x) är derivatan, som mäter hur brant kurvan stiger i varje punkt. Att kvadrera den, addera 1 och ta kvadratroten omvandlar lutningen till den lokala töjningsfaktorn: hur mycket längre kurvan är än en platt horisontell linje på den platsen. Där kurvan är platt är f'(x) = 0 och faktorn är 1. Där kurvan är brant växer faktorn.
Denna formel kräver att funktionen är tillräckligt jämn: f måste ha en kontinuerlig derivata på intervallet, och integranden √(1 + (f'(x))2) måste vara integrerbar. För de flesta funktioner som förekommer i analyskurser och ingenjörsarbete håller dessa villkor. Problem uppstår med skarpa hörn, vertikala tangenter eller diskontinuiteter, där derivatan beter sig illa.
Svårigheten är att denna integral ofta inte har någon elementär primitiv funktion. Även för enkla funktioner som y = x2 involverar den primitiva funktionen inversa hyperboliska funktioner. För många realistiska kurvor finns det ingen sluten primitiv funktion alls. Det är precis här numeriska metoder som Simpsons regel blir nödvändiga: de uppskattar integralen direkt, utan att behöva en primitiv funktion.
Hur Simpsons regel beräknar båglängd
Att tillämpa Simpsons regel på båglängd innebär en tydlig följd av steg. Den centrala insikten är att Simpsons regel inte tillämpas på den ursprungliga funktionen f(x). Den tillämpas på båglängdsintegranden g(x) = √(1 + (f'(x))2). Denna distinktion är kritisk och en vanlig källa till fel.
Här är hela processen:
- Börja med funktionen. Identifiera y = f(x) och bekräfta att den är deriverbar på ditt intervall.
- Beräkna derivatan. Hitta f'(x) analytiskt eller verifiera att den är tillgänglig.
- Bygg integranden. Bilda g(x) = √(1 + (f'(x))2). Detta är funktionen som Simpsons regel kommer att integrera.
- Identifiera intervallet. Bekräfta den nedre gränsen a och den övre gränsen b för det kurvsegment du vill mäta.
- Välj delintervall. Välj ett jämnt antal n delintervall. Antalet måste vara jämnt eftersom Simpsons regel arbetar i par.
- Beräkna steglängden. Beräkna h = (b - a) / n, bredden på varje delintervall.
- Hitta stödpunkterna. Punkterna är x0 = a, x1 = a + h, x2 = a + 2h, och så vidare upp till xn = b.
- Utvärdera integranden. Beräkna g(x0), g(x1), ..., g(xn) i varje stödpunkt.
- Tillämpa vikterna. Multiplicera med Simpsons mönster: 1 för ändpunkter, 4 för inre punkter med udda index, 2 för inre punkter med jämnt index.
- Kombinera och skala. Addera de viktade värdena och multiplicera med h/3 för att få båglängdsuppskattningen.
Den sammansatta Simpsons 1/3-regelns formel för båglängd är:
L ≈ (h/3)[g(x0) + 4g(x1) + 2g(x2) + 4g(x3) + ... + 4g(xn-1) + g(xn)]
där h = (b - a)/n och g(x) = √(1 + (f'(x))2).
Kravet att n ska vara jämnt kommer från parstrukturen. Varje parabel spänner över två delintervall (tre punkter). Om n vore udda skulle det sista delintervallet sakna partner, och metoden skulle bryta samman. När du av någon anledning behöver ett udda antal paneler används i stället trapetsregeln eller en blandad ansats.
Lägg noga märke till viktmönstret. Det första och sista värdet (g(x0) och g(xn)) förekommer vardera en gång med vikt 1. Inre punkter med udda index (x1, x3, x5, ...) får vikt 4 eftersom var och en tjänar som mittpunkt för sin paraboliska panel, där parabelns inflytande är starkast. Inre punkter med jämna index (x2, x4, ...) får vikt 2 eftersom var och en tjänar som en delad ändpunkt mellan två angränsande paneler.
Simpsons regels formel förklarad steg för steg
Den sammansatta Simpsons regels formel har flera komponenter som var och en spelar en specifik roll. Att förstå dem individuellt gör metoden mycket mindre skrämmande.
Steglängd (h). Detta är bredden på varje delintervall, beräknad som h = (b - a)/n. Mindre h betyder fler paneler och generellt bättre noggrannhet, till priset av fler funktionsutvärderingar. För båglängd styr h hur fint kurvan samplas.
Integrationsgränser (a och b). Dessa definierar det kurvsegment som mäts. Valet av intervall påverkar resultatet dramatiskt: att förlänga b längre längs kurvan lägger till mer längd. Dubbelkolla alltid att a < b och att båda ligger inom funktionens definitionsmängd.
Stödpunkter. De n+1 punkterna x0 till xn är jämnt fördelade över [a, b]. Simpsons regel utvärderar integranden i var och en av dessa punkter. Ingen punkt hoppas över.
Ändpunktsbidrag. Värdena g(x0) och g(xn) ingår vardera en gång i summan, med vikt 1. De förankrar approximationen vid gränserna.
Bidrag med udda index. Punkterna x1, x3, x5, och så vidare får vikt 4. Dessa är mittpunkterna i varje parabolisk panel. Vikten 4 återspeglar att mittpunkten har störst inflytande på arean under en parabel.
Inre bidrag med jämnt index. Punkterna x2, x4, x6, och så vidare (exklusive ändpunkterna) får vikt 2. Var och en av dessa punkter delas som ändpunkt mellan två angränsande paneler, så dess totala bidrag fördubblas.
Viktad summering. Multiplicera varje funktionsvärde med dess vikt, addera allt och multiplicera sedan med h/3. Faktorn h/3 skalar den viktade summan till rätt enheter och storlek.
Följande tabell visar viktmönstret för n = 6:
| Punkt | x-värde | Vikt | Roll |
|---|---|---|---|
| x0 | a | 1 | Vänster ändpunkt |
| x1 | a + h | 4 | Panel 1 mittpunkt |
| x2 | a + 2h | 2 | Delad ändpunkt |
| x3 | a + 3h | 4 | Panel 2 mittpunkt |
| x4 | a + 4h | 2 | Delad ändpunkt |
| x5 | a + 5h | 4 | Panel 3 mittpunkt |
| x6 | b | 1 | Höger ändpunkt |
Så använder du Simpsons regels kalkylator
Kalkylatorn på den här sidan uppskattar båglängd med Simpsons regel i flera inmatningslägen. Varje läge hanterar ett sätt att beskriva en kurva. Att förstå vilket läge som passar ditt problem är det första steget.
Funktionsläge. Ange en formel som x^2, sin(x) eller exp(x), plus intervallgränserna a och b. Kalkylatorn deriverar numeriskt eller använder strukturen du anger, bygger båglängdsintegranden och tillämpar Simpsons regel. Detta är det vanligaste läget för analysproblem.
Radie- och vinkelläge. För cirkelbågar, ange radien och centralvinkeln. Båglängden följer direkt från geometrin, och Simpsons regel behövs inte för själva cirkelfallet, men läget är användbart för jämförelse.
Parametriskt läge. Ange x(t) och y(t) som funktioner av en parameter t, plus parameterintervallet. Kalkylatorn använder den parametriska båglängdsformeln med Simpsons regel tillämpad på √((dx/dt)2 + (dy/dt)2).
Polärt läge. Ange r(θ) och vinkelintervallet. Simpsons regel integrerar √(r2 + (dr/dθ)2) över vinkelområdet.
Punktläge. Ange en lista med koordinatpar. Kalkylatorn summerar rätlinjeavstånden mellan på varandra följande punkter. Detta är en polygonal approximation snarare än Simpsons regel, användbar när du har uppmätta data i stället för en formel.
3D-parametriskt läge. Ange x(t), y(t) och z(t) för en rymdkurva. Simpsons regel integrerar √((dx/dt)2 + (dy/dt)2 + (dz/dt)2).
Så här kör du en båglängdsberäkning med Simpsons regel i funktionsläge:
- Ange funktionen. Skriv formeln med standardnotation, till exempel
x^2ellersin(x). - Ställ in den nedre gränsen. Ange startvärdet a för x för kurvsegmentet.
- Ställ in den övre gränsen. Ange slutvärdet b för x. Bekräfta att a < b.
- Välj delintervall. Välj ett jämnt tal n. Börja med ett måttligt värde som 20 eller 50.
- Bekräfta metoden. Den här sidan låser den numeriska metoden till Simpsons regel, så inget val behövs.
- Kör beräkningen. Verktyget utvärderar integranden i n+1 punkter och returnerar uppskattningen.
- Kontrollera konvergensen. Kör om med ett större n (till exempel det dubbla) och jämför. Om resultatet knappt ändras är det tillförlitligt.
För ett realistiskt exempel, prova f(x) = x2 på [0, 1] med n = 20, sedan n = 40. De två resultaten bör stämma överens på flera decimaler, vilket bekräftar att uppskattningen är stabil.
Räknat exempel: Beräkna båglängden av y = x<sup>2</sup> med Simpsons regel
Detta avsnitt går igenom ett komplett exempel med varje mellansteg visat. Funktionen är y = x2 på intervallet [0, 1], med n = 4 delintervall. Alla värden är beräknade oberoende och kan verifieras med valfri kalkylator.
Steg 1: Funktionen. Vi börjar med y = x2, en enkel parabel.
Steg 2: Derivatan. Derivering ger y' = 2x.
Steg 3: Båglängdsintegranden. Insättning i formeln ger: g(x) = √(1 + (2x)2) = √(1 + 4x2).
Steg 4: Integralen. Den exakta båglängden är L = ∫01 √(1 + 4x2) dx. Denna integral har en sluten form som involverar inversa hyperboliska funktioner, och dess värde är ungefär 1.4789428575. Vi ska se hur nära Simpsons regel kommer.
Steg 5: Steglängd. Med n = 4 på [0, 1] är h = (1 - 0)/4 = 0.25.
Steg 6: Stödpunkter. De fem jämnt fördelade punkterna är:
- x0 = 0.00
- x1 = 0.25
- x2 = 0.50
- x3 = 0.75
- x4 = 1.00
Steg 7: Funktionsvärden. Utvärdera g(x) = √(1 + 4x2) i varje punkt:
- g(0.00) = √(1 + 0) = 1.0000000000
- g(0.25) = √(1 + 0.25) = √1.25 ≈ 1.1180339887
- g(0.50) = √(1 + 1) = √2 ≈ 1.4142135624
- g(0.75) = √(1 + 2.25) = √3.25 ≈ 1.8027756377
- g(1.00) = √(1 + 4) = √5 ≈ 2.2360679775
Steg 8: Tillämpa vikter. Simpsons mönster för n = 4 är 1, 4, 2, 4, 1:
- 1 × g(x0) = 1 × 1.0000000000 = 1.0000000000
- 4 × g(x1) = 4 × 1.1180339887 = 4.4721359550
- 2 × g(x2) = 2 × 1.4142135624 = 2.8284271247
- 4 × g(x3) = 4 × 1.8027756377 = 7.2111025508
- 1 × g(x4) = 1 × 2.2360679775 = 2.2360679775
Steg 9: Viktad summa. Addition ger 1.0000000000 + 4.4721359550 + 2.8284271247 + 7.2111025508 + 2.2360679775 = 17.7477336082.
Steg 10: Slutresultat. Multiplicera med h/3 = 0.25/3 ≈ 0.0833333333: L ≈ 0.0833333333 × 17.7477336082 ≈ 1.4789778007.
Jämfört med det exakta värdet 1.4789428575 är felet cirka 0.0000349, eller ungefär 3.5 × 10-5. Med endast 4 delintervall (5 funktionsutvärderingar) uppnår Simpsons regel redan 4 decimalers noggrannhet. Detta är metodens styrka: snabb konvergens på jämna integrander.
Som jämförelse ger en fördubbling till n = 6 värdet 1.4789405084 (fel ≈ 2.4 × 10-6), och även det grova n = 2 ger 1.4821537045 (fel ≈ 0.0032). Varje fördubbling av n minskar felet dramatiskt, vilket är karakteristiskt för Simpsons konvergens av fjärde ordningen.
Förstå noggrannheten hos Simpsons regel
Varje numerisk metod ger en approximation, inte ett exakt svar. Skillnaden mellan approximationen och det sanna värdet kallas felet. Att förstå hur detta fel beter sig hjälper dig att välja lämpliga inställningar och lita på dina resultat.
Simpsons regel är en metod av fjärde ordningen, vilket betyder att felet krymper proportionellt mot h4 (eller motsvarande 1/n4) när steglängden minskar. I praktiska termer minskar en fördubbling av antalet delintervall typiskt felet med en faktor på cirka 16. Det är därför det räknade exemplet ovan förbättrades så dramatiskt från n = 2 till n = 4 till n = 6.
Den klassiska felgränsen för den sammansatta Simpsons regeln är:
|Es| ≤ M(b - a)5 / (180n4)
Här är M en övre gräns för absolutbeloppet av integranden g(x):s fjärde derivata över intervallet [a, b]. Det vill säga M ≥ |g(4)(x)| för alla x i [a, b]. Den fjärde derivatan mäter hur snabbt krökningen hos själva integranden förändras. När g är mycket jämn (dess fjärde derivata är liten) garanterar gränsen ett litet fel. När g slingrar sig snabbt är M stort och gränsen mindre betryggande.
Flera viktiga förbehåll gäller. För det första kräver denna gräns att g har en kontinuerlig fjärde derivata på intervallet. Om integranden har en knick, en spets eller en diskontinuitet i en högre derivata gäller inte gränsen. För det andra är det ofta svårt att beräkna M i praktiken, så gränsen är mer användbar för teoretisk förståelse än för vardaglig feluppskattning. För det tredje kan gränsen vara mycket pessimistisk: faktiska fel är ofta mycket mindre än vad formeln förutspår.
I praktiken är den mest tillförlitliga noggrannhetskontrollen konvergenstestning. Beräkna båglängden med n delintervall, sedan med 2n, sedan med 4n. Om resultaten stämmer överens på det antal decimaler du behöver är uppskattningen pålitlig. Om de fortsätter att ändras väsentligt, öka n ytterligare eller undersök om integranden har problematiskt beteende (skarpa svängar, nästan singulariteter, snabb oscillation).
Den ursprungliga funktionen f:s jämnhet spelar in indirekt genom g. Om f har ett skarpt hörn hoppar f' diskontinuerligt, och g = √(1 + (f')2) ärver en knick. Simpsons regel fungerar fortfarande, men konvergensen saktar ner. I sådana fall återställer en uppdelning av intervallet vid problempunkten och separat tillämpning av Simpsons regel på varje jämn bit den snabba konvergensen.
Simpsons regel mot trapetsregeln
Trapetsregeln är Simpsons regels närmaste släkting och vanligaste alternativ. Båda är Newton-Cotes formler med jämnt fördelade punkter, men de skiljer sig fundamentalt i hur de approximerar integranden mellan stödpunkterna.
Trapetsregeln förbinder på varandra följande stödpunkter med räta linjer och bildar trapetser. Simpsons regel anpassar parabler genom grupper om tre punkter. Eftersom en parabel kan kurva sig för att följa integranden medan en rät linje inte kan, fångar Simpsons regel krökning som trapetsregeln missar. Denna geometriska skillnad översätts direkt till noggrannhet: Simpsons regel är av fjärde ordningen (fel proportionellt mot 1/n4) medan trapetsregeln är av andra ordningen (fel proportionellt mot 1/n2).
I praktiska termer betyder detta att Simpsons regel ofta uppnår med n = 20 vad trapetsregeln behöver n = 200 eller mer för att matcha, åtminstone för jämna integrander. Kostnaden är att Simpsons regel kräver ett jämnt n och något noggrannare bokföring med de alternerande vikterna. Trapetsregeln fungerar med valfritt n och har ett enklare viktmönster (1, 2, 2, ..., 2, 1).
Ingen av metoderna är universellt överlägsen. Trapetsregeln har en överraskande fördel för periodiska funktioner integrerade över hela perioder, där den kan konvergera exponentiellt snabbt. Simpsons regel kan ha svårt med funktioner som inte är tillräckligt jämna för att dess fjärde ordningens teori ska gälla. För grova eller brusiga data är den enklare trapetsregeln ibland mer robust.
| Kriterium | Simpsons regel | Trapetsregeln |
|---|---|---|
| Approximationsform | Parabelsegment | Rätlinjesegment |
| Konvergensordning | Fjärde ordningen (1/n4) | Andra ordningen (1/n2) |
| Krav på delintervall | n måste vara jämnt | Valfritt positivt n |
| Viktmönster | 1, 4, 2, 4, ..., 4, 1 | 1, 2, 2, ..., 2, 1 |
| Funktionsutvärderingar | n + 1 | n + 1 |
| Bäst för | Jämna integrander som kräver hög noggrannhet | Enkla fall, periodiska funktioner, grova data |
| Nödvändig jämnhet | Kontinuerlig fjärde derivata för full teori | Kontinuerlig andra derivata för full teori |
Ett sunt arbetssätt är att använda Simpsons regel som primär metod och korskontrollera med trapetsregeln vid högre upplösning när integrandens beteende är osäkert. Om båda metoderna är överens är tilliten hög. Om de skiljer sig väsentligt behöver integranden sannolikt fler delintervall eller intervalluppdelning.
Olika matematiska ansatser till båglängd
Kurvor kan beskrivas i flera koordinatsystem, och var och en har sin egen båglängdsformel. Simpsons regel gäller för alla, eftersom varje formel är en bestämd integral av en jämn integrand.
Kartesisk båglängd
För y = f(x) på [a, b] är formeln L = ∫ab √(1 + (dy/dx)2) dx. Detta är standardformen som lärs ut i analyskurser. Simpsons regel integrerar g(x) = √(1 + (f'(x))2) som beskrivs genomgående i denna guide. Den fungerar bra när f är deriverbar med en kontinuerlig derivata.
Parametrisk båglängd
För en kurva beskriven av x = x(t) och y = y(t) med t från a till b är formeln L = ∫ab √((dx/dt)2 + (dy/dt)2) dt. Parametrisk form hanterar kurvor som inte klarar vertikalinjetestet, såsom cirklar, ellipser och cykloider, som inte kan skrivas som y = f(x). Simpsons regel integrerar den parametriska hastighetsfunktionen √((x'(t))2 + (y'(t))2), som representerar hur snabbt en punkt rör sig längs kurvan när t ökar.
Polär båglängd
För r = r(θ) med θ från α till β är formeln L = ∫αβ √(r2 + (dr/dθ)2) dθ. Polära koordinater beskriver spiraler, kardioider och roskurvor naturligt. Integranden kombinerar det radiella avståndet r med den vinkelmässiga förändringshastigheten dr/dθ. Simpsons regel tillämpas direkt och integrerar med avseende på θ.
Tredimensionell båglängd
För en rymdkurva x = x(t), y = y(t), z = z(t) utvidgas formeln naturligt: L = ∫ab √((dx/dt)2 + (dy/dt)2 + (dz/dt)2) dt. Varje koordinat bidrar med sin kvadrerade hastighetskomponent. Detta är väsentligt för spiralbanor, 3D-banor och kurvor i datorgrafik. Simpsons regel hanterar 3D-integranden på samma sätt som 2D-varianterna.
I varje fall är strategin identisk: identifiera integranden, bekräfta att den är jämn på intervallet, välj ett jämnt n och tillämpa den viktade summan. Kalkylatorn på den här sidan stöder funktions-, parametriskt, polärt och 3D-parametriskt läge, vart och ett kopplat genom Simpsons regel på lämpligt sätt.
Praktiska tillämpningar av Simpsons regel
Simpsons regel för båglängd förekommer i många fält där krökta avstånd måste uppskattas från matematiska modeller eller uppmätta data.
Maskinteknik. Kamprofiler, kugghjulständer och styrskenor designas som jämna matematiska kurvor. Att uppskatta deras längd hjälper till att bestämma materialbehov, kontaktsträckor och slitage mönster. Simpsons regel ger snabba numeriska kontroller under designiterationscykeln.
Väg- och transportgeometri. Motorvägskurvor, järnvägsövergångar och banlayouts följer noggrant designade geometriska profiler (klotoider, cirkelbågar, paraboliska vertikalkurvor). Att beräkna den sanna längden av dessa sträckningar spelar roll för kostnadsuppskattning, skyltplacering och byggplanering.
CAD och datorstödd design. Designprogramvara representerar kurvor som splines och Beziersegment. Numerisk integration uppskattar längden av dessa sammansatta kurvor för tillverkning, där skärbanor och materialåtgång beror på noggranna längddata.
Robotik och rörelseplanering. En robotarm eller ett autonomt fordon följer en planerad bana genom rymden. Den totala banlängden påverkar restid, energiförbrukning och uppdragsplanering. Simpsons regel integrerar hastighetsprofilen längs parametriserade banor.
Vetenskapliga beräkningar. Forskare som modellerar fysikaliska fenomen, från partikelbanor till fältlinjer, behöver rutinmässigt kurvlängder. När de styrande ekvationerna inte ger någon sluten båglängd är numerisk integration standardansatsen.
Byggnadsteknik. Bågformade broar, krökta balkar och skalkonstruktioner har längder och ytmått härledda från deras geometriska definitioner. Noggrann båglängdsberäkning stöder lastanalys och materialavdrag.
Matematisk modellering och forskning. I ren och tillämpad matematik förekommer båglängd i differentialgeometri, variationsproblem och optimering. Numeriska uppskattningar ger intuition och verifiering för teoretiska resultat.
I alla dessa sammanhang värdesätts Simpsons regel eftersom den ger hög noggrannhet med relativt få funktionsutvärderingar, vilket gör den effektiv för både interaktiva verktyg och storskaliga beräkningar.
Vanliga misstag vid användning av Simpsons regel
Flera återkommande fel snubblar studenter och praktiker. Medvetenhet om dessa fallgropar förhindrar de flesta dåliga resultat.
Använda ett udda antal delintervall. Simpsons 1/3-regel kräver fundamentalt ett jämnt n. Med ett udda n bryter parstrukturen samman. Verifiera alltid att n är jämnt innan du beräknar. Om ditt problem antyder ett udda antal, avrunda uppåt till nästa jämna tal.
Fel derivata. Integranden använder f'(x), inte f(x). Ett vanligt misstag är att derivera felaktigt (till exempel att skriva derivatan av x2 som x i stället för 2x) eller att glömma kedjeregeln för sammansatta funktioner. Dubbelkolla varje derivata.
Tillämpa Simpsons på f i stället för g. Besläktat med ovanstående tillämpar vissa användare vikterna direkt på den ursprungliga funktionens värden f(xi) snarare än på båglängdsintegranden g(xi) = √(1 + (f'(xi))2). Detta beräknar en helt fel integral.
Felaktiga intervallgränser. Att ange a > b, eller gränser utanför funktionens definitionsmängd, ger meningslösa resultat. För funktioner med begränsade definitionsmängder (som √(x) som kräver x ≥ 0), verifiera att intervallet är giltigt.
För få delintervall. Även om Simpsons regel konvergerar snabbt kan n = 2 eller n = 4 vara otillräckligt för starkt krökta integrander. Utför alltid en konvergenskontroll med större n.
Grad-radian-förvirring. När f involverar trigonometriska funktioner, bekräfta om kalkylatorn förväntar sig radianer eller grader. Att blanda dem korrumperar tyst derivatan och resultatet. De flesta matematiska verktyg använder radianer.
Anta exakthet. Simpsons regel är exakt endast för polynom upp till grad 3. För allmänna integrander är resultatet en approximation. Att behandla den som exakt utan konvergensbevis är ett misstag.
Ignorera singulariteter. Om f' exploderar inom [a, b] (till exempel nära en vertikal tangent) blir integranden g stor och Simpsons regel kan konvergera långsamt eller otillförlitligt. Dela intervallet eller använd specialiserade metoder.
När ska du använda Simpsons regel?
Simpsons regel är ett utmärkt standardval för numerisk båglängdsberäkning, men att veta när den glänser och när man ska se sig om efter annat gör dig till en mer effektiv problemlösare.
Använd Simpsons regel när:
- Integranden g(x) = √(1 + (f'(x))2) är jämn på intervallet, utan knickar eller singulariteter.
- Ingen elementär primitiv funktion finns, eller att hitta en är mer möda än en numerisk uppskattning motiverar.
- Du behöver god noggrannhet med ett måttligt antal funktionsutvärderingar.
- Du beräknar båglängd för kursarbete, ingenjörskontroller eller designvalidering.
- Kurvan beskrivs parametriskt eller i polära koordinater och den resulterande integranden är jämn.
Överväg alternativ eller extra försiktighet när:
- Integranden har en diskontinuitet eller ett skarpt hörn inom intervallet. Dela intervallet vid problempunkten först.
- Funktionen oscillerar snabbt. Du kan behöva mycket stort n, eller en metod designad för oscillerande integraler.
- Du bara har brusiga uppmätta data snarare än en jämn formel. Trapetsregeln eller splinebaserade ansatser kan vara lämpligare.
- Extremt hög precision krävs. Gaussisk kvadratur eller adaptiv Simpsons regel (som förfinar delintervallen automatiskt där det behövs) kan vara bättre.
- Intervallet är oändligt eller halv-oändligt. Transformera integralen till ett ändligt intervall först, eller använd specialiserad kvadratur.
För de flesta analysstudenter och arbetande ingenjörer ger Simpsons regel med en konvergenskontroll (jämför n, 2n och 4n) all noggrannhet som behövs. Den är enkel att implementera, allmänt tillgänglig i programvara och teoretiskt väl förstådd.
Relaterade verktyg
Vanliga frågor om Simpsons regel
Vad används Simpsons regel till? +
Simpsons regel approximerar bestämda integraler numeriskt. För båglängd uppskattar den L = ∫<sub>a</sub><sup>b</sup> √(1 + (f'(x))<sup>2</sup>) dx genom att utvärdera integranden i jämnt fördelade punkter och kombinera dem med viktmönstret 1, 4, 2, 4, ..., 4, 1. Den används flitigt när ingen elementär primitiv funktion finns eller när en snabb numerisk uppskattning räcker.
Hur beräknar Simpsons regel båglängd? +
Den tillämpar Simpsons viktningsformel på båglängdsintegranden g(x) = √(1 + (f'(x))<sup>2</sup>), inte på den ursprungliga funktionen. Intervallet [a, b] delas i ett jämnt antal n delintervall med bredden h = (b - a)/n. Uppskattningen är L ≈ (h/3)[g(x<sub>0</sub>) + 4g(x<sub>1</sub>) + 2g(x<sub>2</sub>) + ... + 4g(x<sub>n-1</sub>) + g(x<sub>n</sub>)].
Varför kräver Simpsons regel ett jämnt antal delintervall? +
Varje parabolisk panel i Simpsons regel spänner över två delintervall (tre punkter). Metoden grupperar delintervallen i par och anpassar en parabel över varje par. Om n vore udda skulle det sista delintervallet sakna partner, och parstrukturen skulle bryta samman. Ett jämnt n säkerställer att varje delintervall tillhör exakt en parabolisk panel.
Vad är skillnaden mellan Simpsons 1/3-regel och trapetsregeln? +
Simpsons 1/3-regel anpassar parabler genom grupper om tre punkter och uppnår noggrannhet av fjärde ordningen (fel proportionellt mot 1/n<sup>4</sup>). Trapetsregeln förbinder punkter med räta linjer och uppnår endast noggrannhet av andra ordningen (fel proportionellt mot 1/n<sup>2</sup>). Simpsons regel behöver typiskt långt färre delintervall för samma noggrannhet på jämna integrander, men den kräver ett jämnt n.
Hur beräknar man steglängden i Simpsons regel? +
Steglängden är h = (b - a)/n, där a och b är den nedre och övre integrationsgränsen och n är det (jämna) antalet delintervall. Till exempel, på [0, 1] med n = 4 är h = 0.25. Mindre h betyder fler paneler och generellt bättre noggrannhet.
Är Simpsons regel mer noggrann än trapetsregeln? +
För jämna integrander, ja, typiskt med stor marginal. Simpsons konvergens av fjärde ordningen betyder att en fördubbling av n minskar felet med ungefär 16x, mot 4x för trapetsregeln. Denna fördel förutsätter dock tillräcklig jämnhet. För grova, brusiga eller periodiska integrander kan trapetsregeln ibland prestera jämförbart eller bättre.
Hur många delintervall ska jag använda? +
Börja med ett måttligt jämnt tal som 20 eller 50, fördubbla det sedan och jämför resultaten. Om uppskattningen knappt ändras har den konvergerat och du kan lita på den. Om den ändras väsentligt, fortsätt att fördubbla. För de jämna integrander som är typiska i båglängdsproblem är n = 50 till 200 vanligtvis mer än tillräckligt.
Kan Simpsons regel beräkna längden av en parametrisk kurva? +
Ja. För x = x(t), y = y(t) på [a, b], tillämpa Simpsons regel på den parametriska hastighetsintegranden √((dx/dt)<sup>2</sup> + (dy/dt)<sup>2</sup>). Metoden är identisk: välj ett jämnt n, utvärdera i jämnt fördelade t-värden, tillämpa vikterna 1-4-2-...-4-1 och multiplicera med h/3.
Kan Simpsons regel användas för tredimensionella kurvor? +
Ja. För en rymdkurva x(t), y(t), z(t) blir båglängdsintegranden √((dx/dt)<sup>2</sup> + (dy/dt)<sup>2</sup> + (dz/dt)<sup>2</sup>). Simpsons regel integrerar denna precis som i 2D-fallet. Kalkylatorn på den här sidan inkluderar ett 3D-parametriskt läge.
Vad är felformeln för Simpsons regel? +
Den sammansatta Simpsons regels felgräns är |E<sub>s</sub>| ≤ M(b - a)<sup>5</sup>/(180n<sup>4</sup>), där M är en övre gräns för |g<sup>(4)</sup>(x)|, absolutbeloppet av integrandens fjärde derivata, över [a, b]. Detta kräver att g är fyra gånger kontinuerligt deriverbar. I praktiken är konvergenstestning (jämföra n, 2n, 4n) mer användbar än att beräkna M.
Ger Simpsons regel ett exakt svar? +
Endast för integrander som är polynom av grad 3 eller lägre. För allmänna jämna funktioner ger den en approximation vars fel krymper som 1/n<sup>4</sup>. Det räknade exemplet i denna guide (y = x<sup>2</sup> på [0,1]) uppnår ett fel under 0.00004 med endast n = 4, vilket visar hur snabbt approximationen blir praktiskt taget exakt.
Hur kan jag kontrollera om mitt numeriska resultat är tillförlitligt? +
Använd konvergenstestning: beräkna med n, sedan 2n, sedan 4n delintervall. Om resultaten stämmer överens på dina krävda decimaler är uppskattningen tillförlitlig. Som en andra kontroll, jämför mot trapetsregeln vid hög upplösning. Överensstämmelse mellan två olika metoder tyder starkt på korrekthet.