Mittpunktsregelns båglängdskalkylator

Mittpunktsregelns båglängdskalkylator approximerar längden av en kurva genom att dela båglängdsintegralen i rektanglar och sampla varje rektangel i dess mittpunkt. Mata in en funktion, välj ett intervall, välj hur många delintervall som ska användas, så visar kalkylatorn approximationen, den fullständiga mittpunktstabellen, feluppskattningen och ett livediagram över rektanglarna. Följ dessa steg för att få ett korrekt svar.

1. Välj beräkningsläge. Använd Explicit funktion för en kurva skriven som y = f(x), till exempel y = x² eller y = sin(x). Använd Parametrisk kurva för en kurva skriven som x = x(t) och y = y(t), till exempel en cirkel eller en ellips som spåras av en parameter.
2. Mata in funktionen. Skriv funktionen i inmatningsrutan med ^ för potenser, kvadratrotsfunktionen för kvadratrötter och sin(), cos(), exp() och ln() för standardfunktionerna. Kalkylatorn hittar derivatan automatiskt, eller så kan du skriva derivatan själv om du redan känner till den.
3. Ställ in intervallet. Mata in den undre gränsen a och den övre gränsen b. Intervallet [a, b] är den del av kurvan vars längd du vill ha. Om du matar in ett b som är mindre än a returnerar kalkylatorn ändå en positiv längd, eftersom båglängd aldrig är negativ.
4. Välj antalet delintervall n. Använd reglaget eller skriv ett värde. Förinställningarna är 1, 2, 4, 5, 8, 10, 20, 50, 100, 500, 1 000 och 5 000. Ett större n ger smalare rektanglar och vanligen en närmare approximation. Börja med n = 4 för att se uträkningen tydligt och höj sedan n för att se hur svaret stabiliseras.
5. Välj visningsprecision. Välj 2, 4, 6, 8 eller 10 decimaler. Kalkylatorn behåller full precision internt och avrundar bara det visade resultatet, så en ändrad precision ändrar aldrig den underliggande beräkningen.
6. Klicka på Beräkna båglängd. Huvudresultatet är mittpunktsapproximationen L_n, tydligt markerad som en approximation. Under den får du delintervallbredden, den första och den sista mittpunkten, det första och det sista integrandvärdet, summan av integrandvärdena, kordans längd mellan ändpunkterna, referenslängden med hög precision samt de uppskattade absoluta och relativa felen.
7. Studera diagrammen. Det första diagrammet visar den ursprungliga kurvan med den valda bågen markerad och mittpunktssampelpunkterna utmärkta. Det andra diagrammet visar båglängdsintegranden med mittpunktsrektanglarna ritade under den. Rektanglarnas totala area är approximationen, så de två diagrammen får inte blandas ihop: rektanglarna hör hemma under integranden, aldrig under den ursprungliga kurvan.
8. Läs beräkningstabellen. Varje delintervall får en rad med sina gränser, sin mittpunkt, derivatan i mittpunkten, integrandvärdet och rektangelns bidrag. Den sista raden visar summan, som är lika med huvudresultatet.
9. Jämför konvergensen och prova Richardson-extrapolation. Konvergenspanelen visar approximationerna för n, 2n och 4n sida vid sida, med skillnaderna mellan dem. Richardson-alternativet kombinerar två approximationer till en förfinad uppskattning. Båda förklaras i guideavsnitten nedan.
10. Kopiera, exportera eller återställ. Använd Kopiera resultat för svaret, exportera tabellen som CSV för din uträkning eller Återställ för att rensa panelen. Förinställningsknapparna laddar verifierade exempel direkt: parabeln, den räta linjen, den horisontella linjen, kvartscirkeln, kuben och sinuskurvan.

Prova ett exempel. y = x² på [0, 2] med n = 4 ger L₄ = 4.6267234873 enheter, mot en referens på 4.6467837624 enheter. y = 3x + 1 på [0, 4] med n = 4 ger exakt 12.6491106407 enheter. y = 2 på [1, 5] med n = 4 ger exakt 4.0000000000 enheter. Den parametriska kvartscirkeln med radie 5 ger 7.8539816340 enheter. y = x³ på [0, 1] med n = 4 ger 1.5330289642 enheter, och y = sin(x) på [0, π] med n = 8 ger 3.8201978627 enheter.

Viktiga tips. Derivera alltid funktionen först: integranden använder f'(x), inte f(x). Sampla i mittpunkter, inte i ändpunkter, annars beräknar du en annan regel. Behåll vinkelenheten radianer för trigonometriska funktioner. Avrunda inte mittpunkterna eller integrandvärdena mitt i beräkningen, eftersom de små avrundningarna adderas över många rektanglar. Om funktionen är odefinierad i en mittpunkt rapporterar kalkylatorn ett fel i stället för att tyst hoppa över samplet.

Numerisk integration är konsten att uppskatta en bestämd integral när en exakt primitiv funktion är svår att hitta, eller när bara samplade värden finns tillgängliga. Mittpunktsregeln är en av de enklaste och mest tillförlitliga metoderna för numerisk integration. Den uppskattar arean under en kurva genom att täcka den med rektanglar och addera deras areor.

En Riemannsumma delar intervallet [a, b] i n lika stora delintervall, vart och ett med bredden Δx = (b - a) / n, och bygger en rektangel på varje delintervall. Rektangelns höjd är funktionsvärdet i en enda sampelpunkt inne i delintervallet. Vänsterändpunktssummor samplar vid den vänstra kanten, högerändpunktssummor samplar vid den högra kanten och mittpunktsregeln samplar i mitten av varje delintervall.

Att sampla i mitten är den centrala idén. På ett kort delintervall stiger en kurva på ena sidan av mittpunkten ungefär lika mycket som den faller på den andra sidan, så överskattningen på ena halvan tar ungefär ut underskattningen på den andra halvan. Den utjämningen gör att mittpunktsregeln är mer exakt än vänster- eller högerändpunktsreglerna med samma antal rektanglar, och att felet krymper med kvadraten på delintervallbredden.

För en funktion g(x) på [a, b] är den sammansatta mittpunktsregeln:

Här är x_i* mittpunkten i det i:te delintervallet, och summan löper över alla n delintervall. Varje term Δx × g(x_i*) är arean av en rektangel, och summan av alla rektangelareor är approximationen av integralen.

Mittpunktsregeln används långt utanför kalkylens läxor. Ingenjörer använder den för att uppskatta integraler i hållfasthetsanalys och signalbehandling, fysiker använder den i numeriska simuleringar och datorgrafik använder den för att mäta kurvlängder och ytareor. Kombinationen av enkelhet och andra ordningens noggrannhet gör den till ett standardval när en integral saknar en prydlig sluten form.

Varför mittpunkter slår ändpunkter. Ta g(x) = x på [0, 1] med n = 1. Vänstersumman ger 0, högersumman ger 1 och mittpunktssumman ger 0.5, vilket är det exakta svaret. Mittpunkten ser funktionens genomsnittsbeteende över delintervallet, medan en ändpunkt bara ser ena sidan av det. Den fördelen gäller för varje glatt funktion, och det är anledningen till att mittpunktsregeln är standard-Riemannsumman i det mesta numeriska arbetet.

Båglängd är avståndet mätt längs en kurva, snarare än det räta avståndet mellan dess ändpunkter. En väg som slingrar sig genom kullar är längre än den räta linje som dras mellan dess två ändar på en karta. En böjd tråd är längre än det gap den spänner över. Båglängd mäter själva den slingrande vägen.

För att hitta den längden kan man tänka sig att kurvan skärs i många pyttesmå bitar. Varje pytteliten bit är så kort att den nästan är rak, så dess längd är nära den räta linjen mellan dess ändar. Enligt Pythagoras sats har en pytteliten bit som rör sig Δx horisontellt och Δy vertikalt längden nära √(Δx² + Δy²). Att addera alla pyttesmå bitar ger den totala längden, och när bitarna krymper mot noll blir summan en integral.

För en kurva y = f(x) ger ett pyttelitet horisontellt steg dx en vertikal förändring dy = f'(x) dx, så det pyttelilla längdelementet är:

Att addera alla längdelement från a till b ger båglängdsintegralen:

Den formeln förklarar varför derivatan är inblandad. En brant kurva ändrar höjd snabbt, så f'(x) är stort, integranden √(1 + [f'(x)]²) är stor och varje enhet horisontellt avstånd bidrar med mer än en enhet kurvlängd. En flat kurva har f'(x) nära noll, integranden är nära 1 och båglängden är nära det horisontella avståndet b - a.

Båglängd är alltid minst lika lång som kordan. Den räta linjen mellan ändpunkterna är den kortaste vägen mellan dem, så kurvan kan aldrig vara kortare än den. För parabeln y = x² på [0, 2] är kordan från (0, 0) till (2, 4) √(20) = 4.4721359550 enheter, medan den sanna båglängden är 4.6467837624 enheter, längre med 0.1746478074 enheter. För sinuskurvan y = sin(x) på [0, π] är kordan π = 3.1415926536 enheter och bågen 3.8201977890 enheter, längre med 0.6786051354 enheter. Att kontrollera att bågen överstiger kordan är ett snabbt rimlighetstest för varje båglängdssvar.

De flesta båglängdsintegraler saknar elementär primitiv funktion. Parabelintegralen ∫ √(1 + 4x²) dx kan lösas med en substitution, men sinusintegralen ∫ √(1 + cos²x) dx kan inte uttryckas med elementära funktioner alls. Det är precis där mittpunktsregeln gör nytta: den uppskattar integralen numeriskt, utan att någon primitiv funktion behövs.

Beräkningen av båglängd med mittpunktsregeln utgår från den exakta båglängdsintegralen och ersätter integralen med en summa av rektangelareor. Definiera båglängdsintegranden:

Då är den exakta båglängden L = ∫_a^b g(x) dx, och mittpunktsregeln approximerar den integralen. Dela [a, b] i n lika stora delintervall. Delintervallbredden är:

Mittpunkten i det i:te delintervallet, där i räknas från 0, är:

Beräkna integranden i varje mittpunkt, multiplicera varje värde med Δx och addera resultaten:

Varje symbol har en exakt betydelse:

En konkret kontroll. För y = x² på [0, 2] med n = 4 är bredden Δx = (2 - 0) / 4 = 0.5, mittpunkterna är 0.25, 0.75, 1.25 och 1.75, och integrandvärdena är 1.1180339887, 1.8027756377, 2.6925824036 och 3.6400549446. Deras summa är 9.2534469747, och multiplikation med 0.5 ger L₄ = 4.6267234873 enheter. Varje tal i den kedjan visas i kalkylatorns steg-för-steg-panel och i dess beräkningstabell.

För en parametrisk kurva x = x(t), y = y(t) gäller samma idé för parameterintervallet [t₁, t₂], där hastigheten v(t) = √([x'(t)]² + [y'(t)]²) spelar integrandens roll. Den parametriska formeln behandlas i ett eget avsnitt nedan.

Beräkningen med mittpunktsregeln följer tio steg, och kalkylatorn utför dem alla automatiskt. Att gå igenom dem för hand en gång gör verktyget mycket lättare att lita på och använda. Det genomarbetade exemplet nedan använder y = x² på [0, 2] med n = 4, samma exempel som kalkylatorn laddar som sin första förinställning.

1. Identifiera funktionen. f(x) = x².
2. Hitta dess första derivata. f'(x) = 2x.
3. Identifiera integrationsintervallet. a = 0, b = 2.
4. Välj antalet delintervall. n = 4.
5. Beräkna delintervallbredden. Δx = (2 - 0) / 4 = 0.5.
6. Bestäm mittpunkten i varje delintervall. x₀* = 0.25, x₁* = 0.75, x₂* = 1.25, x₃* = 1.75. Varje mittpunkt ligger exakt halvvägs mellan sitt delintervalls gränser: 0.25 är halvvägs mellan 0 och 0.5, 0.75 är halvvägs mellan 0.5 och 1, och så vidare.
7. Beräkna båglängdsintegranden i varje mittpunkt. g(x) = √(1 + 4x²), så g(0.25) = 1.1180339887, g(0.75) = 1.8027756377, g(1.25) = 2.6925824036 och g(1.75) = 3.6400549446.
8. Multiplicera varje värde med delintervallbredden. De fyra rektangelbidragen är 0.5 × 1.1180339887 = 0.5590169944, 0.5 × 1.8027756377 = 0.9013878189, 0.5 × 2.6925824036 = 1.3462912018 och 0.5 × 3.6400549446 = 1.8200274723.
9. Addera alla bidrag. 0.5590169944 + 0.9013878189 + 1.3462912018 + 1.8200274723 = 4.6267234873.
10. Tolka den approximativa båglängden. L₄ = 4.6267234873 enheter. Referensen med hög precision är 4.6467837624 enheter, så denna approximation är kort med 0.0200602751 enheter, ett relativt fel på 0.4317%.

Samma tio steg fungerar för varje funktion. För y = x³ på [0, 1] med n = 4 är derivatan f'(x) = 3x², bredden är Δx = 0.25, mittpunkterna är 0.125, 0.375, 0.625 och 0.875, integrandvärdena är 1.0010980300, 1.0853471867, 1.5405489332 och 2.5051217067, deras summa är 6.1321158566, och multiplikation med 0.25 ger L₄ = 1.5330289642 enheter, mot en referens på 1.5478656547 enheter.

Ett andra genomarbetat exempel med en trigonometrisk funktion. För y = sin(x) på [0, π] med n = 8 är derivatan f'(x) = cos(x), bredden är Δx = π / 8 = 0.3926990817, och de åtta mittpunkterna är π/16, 3π/16, 5π/16, 7π/16, 9π/16, 11π/16, 13π/16 och 15π/16. Integrandvärdena, i ordning, är 1.4006926023, 1.3005159423, 1.1439660326, 1.0188524102, 1.0188524102, 1.1439660326, 1.3005159423 och 1.4006926023. Deras summa är 9.7280539749, och multiplikation med 0.3926990817 ger L₈ = 3.8201978627 enheter, mot en referens på 3.8201977890 enheter, ett absolut fel på bara 0.0000000736 enheter.

Varje mittpunktsrektangel är en liten, ärlig bild av en bit av integralen. Den har tre delar, och varje del har en tydlig betydelse.

För parabeln y = x² på [0, 2] med n = 4 har alla fyra rektanglar bredden 0.5. Deras höjder är 1.1180339887, 1.8027756377, 2.6925824036 och 3.6400549446, och deras areor är 0.5590169944, 0.9013878189, 1.3462912018 och 1.8200274723. Den totala arean är 4.6267234873 enheter, vilket är exakt mittpunktsapproximationen L₄.

Varför summan av areorna approximerar integralen. Integralen ∫_a^b g(x) dx är den exakta arean under integrandkurvan. Varje rektangel ersätter den krökta remsan ovanför ett delintervall med en platt remsa av samma bredd, där mittpunktshöjden används som representativ höjd. Där integranden böjer sig uppåt skär den platta toppen av lite area på ena sidan och lägger till lite på den andra. Eftersom samplet tas i mitten tar de två felen nästan ut varandra, och rektangelarean ligger nära den sanna remsarean. Att addera rektanglarna ger en uppskattning av hela arean, och därmed av båglängden.

Det vanligaste diagrammisstaget. Rektanglarna måste ritas under integrandkurvan g(x) = √(1 + [f'(x)]²), inte under den ursprungliga kurvan f(x). Den ursprungliga kurvan och integranden är olika grafer med olika höjder: för y = x² på [0, 2] går den ursprungliga kurvan från 0 till 4, medan integranden går från 1 till √(17) = 4.1231056256. Rektanglar ritade under den ursprungliga kurvan skulle summeras till arean under f(x), vilket inte har något med båglängd att göra. Kalkylatorn håller de två graferna åtskilda, och guiden visar rektanglarna bara under integranden i diagrammet ovan.

Att läsa en rektangel. Att hovra över eller välja en rektangel i kalkylatorn visar dess delintervallnummer, dess vänstra och högra gränser, dess mittpunkt, derivatan i mittpunkten, integrandvärdet och rektangelns bidrag. Den detaljen gör tabellen användbar för läxor: varje tal i det slutliga svaret kan spåras tillbaka till en rektangel.

Antalet delintervall n är den enskilt viktigaste kontrollen i kalkylatorn, eftersom det avgör hur fina rektanglarna är. Ett litet n betyder några få breda rektanglar. Ett stort n betyder många smala rektanglar som följer integrandkurvan noga.

Vad som händer när n är litet. Med n = 1 är hela intervallet en enda rektangel, och svaret är Δx × g(mittpunkt). För parabeln y = x² på [0, 2] ger n = 1 2 × √(5) = 4.4721359550 enheter, vilket är kort mot det sanna 4.6467837624 enheter med 0.1746478074 enheter. Den enda rektangeln kan inte följa integranden när den stiger från 1 till 4.1231056256, så felet är stort. Med n = 2 förbättras svaret till 4.5764912225 enheter. Med n = 4 är det 4.6267234873 enheter. Med n = 8 är det 4.6417306581 enheter. Med n = 16 är det 4.6455205322 enheter, och med n = 32 är det 4.6464679593 enheter.

Hur felet krymper. För en glatt integrand är felet i mittpunktsregeln proportionellt mot 1 / n². Att dubblera n skär felet till ungefär en fjärdedel. Parabeltalen visar detta tydligt: felet faller från 0.0200602751 enheter vid n = 4 till 0.0050531044 enheter vid n = 8, sedan till 0.0012632302 enheter vid n = 16, sedan till 0.0003158031 enheter vid n = 32. Varje dubblering dividerar felet med ungefär fyra. Det beteendet 1 / n² kallas andra ordningens konvergens, och det är anledningen till att mittpunktsregeln är så effektiv: en måttlig ökning av n ger en stor förbättring av noggrannheten.

Varför förbättringen inte alltid är monoton. Felet krymper vanligen när n växer, men det måste inte krympa vid varje steg för varje funktion. Sinuskurvan visar en mild variant av detta: n = 2 ger 3.8476494905 enheter, vilket överskrider referensen 3.8201977890 enheter, medan n = 4 ger 3.8204517876 enheter, vilket är mycket närmare. Approximationen korsade från ovan det sanna värdet till ett värde strax ovanför det, och felet krympte ändå totalt sett. För funktioner med skarpa böjar, svängningar eller nästan singulariteter kan felet vackla när n ändras, även om den långsiktiga trenden är nedåtgående. Anta aldrig att ett större n alltid ger ett bättre svar i ett enskilt steg; kontrollera konvergenstabellen i stället.

Kostnaden för ett större n. Fler delintervall betyder mer aritmetik: fler derivataberäkningar, fler integrandvärden och fler tabellrader. För n = 5 000 utför kalkylatorn 5 000 beräkningar, vilket är snabbt för en dator men ger en tabell som är för lång att läsa rad för rad, så kalkylatorn sidindelar den. Den praktiska regeln är att höja n tills svaret slutar ändras vid den precision du behöver. Om n = 100 och n = 200 stämmer överens till sex decimaler räcker n = 100, och n = 5 000 tillför inget annat än arbete.

Att välja n i praktiken. Använd n = 4 eller n = 8 medan du lär dig, eftersom tabellen då hålls kort nog att kontrollera för hand. Använd n = 50 eller n = 100 för läxsvar som behöver fyra eller fem korrekta decimaler. Använd n = 1 000 eller mer bara när du behöver mycket hög precision, eller när integranden ändras snabbt över intervallet. Konvergenspanelen visar n, 2n och 4n tillsammans, så att du kan se svaret stabiliseras utan att gissa.

Mittpunktsregeln är en av tre klassiska metoder för numerisk integration, och var och en samplar integranden på olika sätt. Att veta hur de skiljer sig hjälper dig att välja rätt verktyg och förstå varför deras svar inte stämmer överens.

Mittpunktsregeln. Samplar i mitten av varje delintervall och bygger rektanglar. Felet är O(1 / n²), samma ordning som trapetsregeln, men mittpunktsregeln är vanligen ungefär dubbelt så exakt som trapetsregeln för samma n, eftersom mittpunktssamplet tar ut det ledande felet. Den behöver bara n beräkningar, en per delintervall, och den fungerar för varje positivt n.

Trapetsregeln. Samplar i ändpunkterna och förenar dem med räta linjer, vilket bildar trapetser. Den delar mittpunktsregelns felordning O(1 / n²), men felet har vanligen motsatt tecken: för en konvex integrand överskattar trapetsregeln medan mittpunktsregeln underskattar, eller tvärtom. Medelvärdet av de två ger kombinationen trapets-mittpunkt, som är grunden för Simpsons regel. Prova trapetsregelns kalkylator för att jämföra de två metoderna på samma funktion.

Simpsons regel. Anpassar en parabel genom de två ändpunkterna och mittpunkten i varje par av delintervall, vilket kräver ett jämnt antal delintervall. Felet är O(1 / n⁴), så den konvergerar mycket snabbare än de andra två för glatta funktioner. Priset är kravet på jämnt n och fler beräkningar per delintervall. Prova Simpsons regels kalkylator när du behöver hög noggrannhet med få delintervall.

En sifferjämförelse. För parabeln y = x² på [0, 2] med n = 4 ger mittpunktsregeln 4.6267234873 enheter, trapetsregeln ger 4.6870560064 enheter och Simpsons regel ger 4.6502010785 enheter, mot referensen 4.6467837624 enheter. Mittpunktsfelet är 0.0200602751 enheter, trapetsfelet är 0.0402722440 enheter och Simpsonfelet är 0.0034173160 enheter. Simpson vinner här, som den vanligen gör för glatta funktioner, men mittpunktsregeln är enklare, behöver inget krav på jämnt n och är det naturliga valet när bara mittpunktssampel finns tillgängliga.

Vilken regel ska du använda? Använd mittpunktsregeln när du vill ha en enkel, tillförlitlig uppskattning, när du samplar en funktion i jämnt fördelade inre punkter, eller när du undervisar om Riemannsummor. Använd trapetsregeln när ändpunktsvärden är de naturliga data, till exempel mätningar tagna i slutet av varje intervall. Använd Simpsons regel när funktionen är glatt och du behöver hög noggrannhet med ett litet n. För båglängd specifikt har mittpunktsregeln en praktisk fördel: den beräknar aldrig integranden i ändpunkterna, där en funktion kan vara odefinierad eller singulär, så den kan hantera kurvor som de andra två reglerna inte kan börja på.

Varje mittpunktsapproximation skiljer sig från den sanna båglängden med ett visst belopp, och att förstå det felet är det som gör ett tal till ett pålitligt svar. Det finns tre olika sätt att tala om fel, och de får inte blandas ihop.

Absolut fel. Det absoluta felet är storleken på gapet mellan approximationen och det sanna värdet: |L - L_n|. För parabeln med n = 4 är det absoluta felet |4.6467837624 - 4.6267234873| = 0.0200602751 enheter. För sinuskurvan med n = 8 är det |3.8201977890 - 3.8201978627| = 0.0000000736 enheter. Absolut fel mäts i samma enheter som själva längden.

Relativt fel. Det relativa felet dividerar det absoluta felet med det sanna värdets storlek, och anges vanligen i procent: |L - L_n| / |L| × 100. För parabeln med n = 4 är det relativa felet 0.4317%. För kuben med n = 4 är det 0.9585%. För sinuskurvan med n = 8 är det 0.00000193%. Relativt fel är det rättvisa sättet att jämföra noggrannhet mellan problem av olika storlek: ett fel på 0.02 enheter betyder lite för en längd på 100 enheter och mycket för en längd på 0.1 enheter.

Felindikator baserad på förfining. I verkligt arbete är det sanna värdet okänt, så det absoluta felet kan inte beräknas direkt. Jämför i stället två approximationer: skillnaden |M_2n - M_n| visar hur mycket svaret rörde sig när indelningen förfinades. För parabeln är |M₈ - M₄| = |4.6417306581 - 4.6267234873| = 0.0150071707 enheter. Den skillnaden är en praktisk varningssignal: om den är stor rör sig svaret fortfarande och n bör höjas. Det är en indikator, inte en garanterad gräns: det sanna felet kan vara större eller mindre än förfiningsskillnaden, och det är vanligen ungefär en tredjedel av den för mittpunktsregeln, eftersom felet krymper med en faktor fyra vid varje dubblering.

Teoretisk felgräns. För en tillräckligt glatt integrand g(x) uppfyller felet i den sammansatta mittpunktsregeln:

Gränsen behöver maximum av integrandens andraderivata över intervallet, vilket ofta är svårt att hitta exakt. För parabelintegranden g(x) = √(1 + 4x²) är andraderivatan g''(x) = 4 / (1 + 4x²)^3/2, vars största värde på [0, 2] är 4 vid x = 0. Gränsen för n = 4 är därför (8 / (24 × 16)) × 4 = 0.0833333333 enheter, och det verkliga felet 0.0200602751 enheter ligger säkert innanför den. För n = 8 är gränsen 0.0208333333 enheter. Gränsen är ärlig men lös: den överskattar det verkliga felet med en faktor fyra här, vilket är typiskt.

När gränsen inte gäller. Felgränsen förutsätter att integranden är glatt, med en begränsad andraderivata. En funktion med ett hörn, ett hopp eller en singularitet bryter det antagandet, och gränsen finns antingen inte eller är oanvändbart stor. Kalkylatorn hittar inte på någon gräns när den inte kan intyga en. I det fallet visar den förfiningsskillnaden och förklarar dess gränser, vilket är det ärliga svaret.

Varför glatthet spelar roll. Felbeteendet 1 / n² beror på att integranden är tillräckligt glatt för att mittpunktsutjämningen ska fungera. En glatt integrand som √(1 + 4x²) konvergerar rent. En integrand med ett skarpt hörn konvergerar långsammare, och en integrand med en singularitet kan konvergera mycket långsamt eller inte alls. Om en dubblering av n knappt ändrar svaret kan funktionen vara glatt och redan exakt. Om en dubblering av n ändrar svaret mycket är antingen funktionen grov eller n fortfarande för litet.

Richardson-extrapolation är ett enkelt knep som förvandlar två mittpunktsapproximation till en mycket bättre uppskattning. Det använder det kända mönstret i mittpunktsregelns fel för att ta ut den ledande feltermen.

För en glatt integrand har mittpunktsapproximationen med n delintervall ett fel som beter sig som C / n², där C är en konstant som beror på funktionen men inte på n. Att dubblera delintervallen till 2n dividerar därför felet med fyra. Det ger två ekvationer med två obekanta, det sanna värdet L och konstanten C:

Att lösa för L eliminerar C och lämnar den Richardson-extrapolerade uppskattningen:

Genomarbetat exempel. För parabeln y = x² på [0, 2] är M₄ = 4.6267234873 enheter och M₈ = 4.6417306581 enheter. Richardson-uppskattningen är L_R = (4 × 4.6417306581 - 4.6267234873) / 3 = 4.6467330483 enheter. Det sanna värdet är 4.6467837624 enheter, så det extrapolerade felet är bara 0.0000507141 enheter, mycket mindre än det ursprungliga M₄-felet på 0.0200602751 enheter. Extrapolationen köpte noggrannheten hos ett mycket större n med bara beräkningarna för n = 4 och n = 8.

Genomarbetat exempel med sinuskurvan. För y = sin(x) på [0, π] är M₈ = 3.8201978627 enheter och M₁₆ = 3.8201977890 enheter. Richardson-uppskattningen är L_R = (4 × 3.8201977890 - 3.8201978627) / 3 = 3.8201977645 enheter, mot referensen 3.8201977890 enheter, ett fel på bara 0.0000000245 enheter.

Genomarbetat exempel med kuben. För y = x³ på [0, 1] är M₄ = 1.5330289642 enheter och M₈ = 1.5441590533 enheter. Richardson-uppskattningen är L_R = (4 × 1.5441590533 - 1.5330289642) / 3 = 1.5478690830 enheter, mot referensen 1.5478656547 enheter, ett fel på 0.0000034283 enheter. Det ursprungliga M₄-felet var 0.0148366905 enheter, så extrapolationen skar det med en faktor på cirka 4327.7.

Vad Richardson-extrapolation inte är. Det extrapolerade värdet är en förfinad uppskattning, inte mittpunktsregelns resultat och inte ett garanterat exakt svar. Det förutsätter att felet verkligen beter sig som C / n², vilket gäller för glatta integrander och misslyckas för grova sådana. Om integranden har ett hörn eller en singularitet kan extrapolationen bli sämre än den ursprungliga approximationen, så kalkylatorn visar alltid det ursprungliga mittpunktsresultatet först och markerar det extrapolerade värdet separat. Richardson-extrapolation är också idén bakom Romberg-integration, som upprepar samma knep på flera nivåer för att nå mycket hög noggrannhet.

Vissa kurvor kan inte skrivas som y = f(x). En cirkel ger till exempel två y-värden för de flesta x-värden, så ingen enskild funktion beskriver den. Parametriska kurvor löser detta genom att beskriva båda koordinaterna som funktioner av en parameter t: x = x(t) och y = y(t). När t löper från t₁ till t₂ spårar punkten (x(t), y(t)) kurvan.

För en parametrisk kurva kommer det pyttelilla längdelementet från rörelsen i båda koordinaterna: ds = √(dx² + dy²). Att skriva dx = x'(t) dt och dy = y'(t) dt ger den parametriska båglängdsintegralen:

Kvantiteten under kvadratroten är hastigheten v(t) = √([x'(t)]² + [y'(t)]²), den hastighet med vilken punkten rör sig längs kurvan när t ändras. Mittpunktsregeln gäller för parameterintervallet precis som för ett x-intervall: dela [t₁, t₂] i n lika stora delar med bredden Δt = (t₂ - t₁) / n, sampla hastigheten i mittpunktsparametrarna t_i* = t₁ + (i + ½)Δt, och addera bidragen:

Genomarbetat exempel: kvartscirkeln. Ta x(t) = 5 cos(t) och y(t) = 5 sin(t) på [0, π/2], med n = 4. Derivatorna är x'(t) = -5 sin(t) och y'(t) = 5 cos(t). Hastigheten är v(t) = √(25 sin²(t) + 25 cos²(t)) = √(25) = 5, en konstant, eftersom sin²(t) + cos²(t) = 1. Parameterbredden är Δt = (π/2 - 0) / 4 = π/8 = 0.3926990817, och mittpunktsparametrarna är π/16, 3π/16, 5π/16 och 7π/16. Varje hastighetsvärde är 5, så varje bidrag är 0.3926990817 × 5 = 1.9634954085, och totalen är 4 × 1.9634954085 = 7.8539816340 enheter, vilket är exakt 5π/2. Mittpunktsregeln är exakt här för varje n, eftersom hastigheten är konstant.

Varför den parametriska formen spelar roll. Många verkliga kurvor är naturligt parametriska: cirklar, ellipser, spiraler, Bezierkurvor i designprogram och banor för rörliga objekt i fysik och robotik. Den explicita formen y = f(x) kan inte beskriva en kurva som vänder tillbaka på sig själv, men den parametriska formen kan det. Kalkylatorn hanterar båda formerna med samma mittpunktsmotor, samma tabell och samma konvergensanalys. För mer om parametriska kurvor, se den parametriska båglängdskalkylatorn.

En varning om parametern. Parametern t är inte samma sak som båglängd om inte hastigheten är 1. Lika stora steg i t ger lika stora steg längs kurvan bara när hastigheten är konstant. För en allmän parametrisk kurva fungerar mittpunktsregeln ändå korrekt, eftersom den integrerar hastigheten, men rektanglarna är lika stora i t, inte lika stora i sträcka längs kurvan. Det är okej: formeln tar hänsyn till det genom v(t).

Mittpunktsregeln är användbar överallt där en integral måste uppskattas numeriskt, och båglängd är bara ett av dess jobb. Dess enkelhet, dess andra ordningens noggrannhet och dess tolerans för ändpunktsproblem gör den till en arbetshäst inom många områden.

Numerisk analys och utbildning. Mittpunktsregeln är ofta den första metoden för numerisk integration som studenter möter, eftersom rektanglar är lätta att rita och mittpunktsidén är lätt att förstå. Den lär ut Riemannsummor, konvergens och felanalys i ett paket. Kalkylatorns tabell och diagram är byggda för precis det syftet: varje tal i svaret kan spåras till en rektangel.

Ingenjörskonst. Ingenjörer uppskattar längder av böjda delar, kablar, rör och vägböjar utifrån matematiska modeller. När kurvan är ett polynom eller en trigonometrisk funktion saknar båglängdsintegralen ofta en sluten form, så en numerisk metod är det enda alternativet. Mittpunktsregeln ger ett snabbt svar med en tydlig feluppskattning, och en höjning av n förfinar det så långt som designen kräver.

Fysik. Fysiker integrerar längs banor för att beräkna arbete, tillryggalagd sträcka och fältstorheter. En partikel som rör sig längs en böjd bana tillryggalägger en sträcka som är lika med båglängden av dess bana, och mittpunktsregeln uppskattar den sträckan från hastighetsdata. Samma metod uppskattar integraler i termodynamik, optik och kvantmekanik, överallt där en analytisk primitiv funktion är utom räckhåll.

Datorgrafik och CAD. Kurvor på en skärm ritas som många korta räta segment, och en kurvas totala längd uppskattas genom att addera de segmenten, vilket är en form av numerisk integration. Teckensnittskonturer, animationsbanor och CAD-verktygsbanor beror alla på exakt kurvlängdsmätning. Mittpunktsregeln ger ett principfast sätt att mäta de längderna i stället för att gissa från en bild.

Kurvmätning och robotik. En robotarm som flyttar sitt gripdon längs en cirkelbåge, en CNC-maskin som skär en böjd kant, eller en lantmätare som mäter en böjd gräns behöver alla längden av en kurva som bara är känd genom sin ekvation eller genom samplade punkter. Mittpunktsregeln förvandlar de samplen till en längd, med en feluppskattning som säger hur långt svaret kan litas på.

Vetenskaplig modellering och beräkningsgeometri. Modeller i biologi, finans och klimatvetenskap integrerar funktioner som bara är kända i diskreta punkter. Mittpunktsregeln är det naturliga valet när data är sampel tagna i mitten av intervall, vilket är vanligt i sensoravläsningar och simuleringsutdata. I beräkningsgeometri mäter den omkretsar av böjda former och längder av splinekurvor.

När du inte ska använda den. Mittpunktsregeln är ett dåligt val när integranden har en singularitet inne i intervallet, när funktionen är diskontinuerlig, eller när mycket hög noggrannhet behövs med få beräkningar: Simpsons regel eller en adaptiv metod gör det bättre där. Den är också fel verktyg när integralen i sig inte är frågan, till exempel när en sluten båglängd redan är känd. En cirkel med radie 5 har en kvartscirkelbåge på exakt 5π/2 = 7.8539816340 enheter, och ingen numerisk metod kan förbättra det.

De flesta felaktiga mittpunktssvar kommer från en kort lista med fel. Att kontrollera varje punkt tar sekunder och sparar de flesta misstag.

1. Att använda ändpunkter i stället för mittpunkter. Att sampla vid den vänstra eller högra kanten av varje delintervall beräknar vänster- eller höger-Riemannsumman, inte mittpunktsregeln. För parabeln y = x² på [0, 2] med n = 4 är vänstersumman 3.75 enheter och högersumman 5.75 enheter, medan mittpunktsregeln ger 4.6267234873 enheter. Bara mittpunktsvärdet är mittpunktsregelns svar.
2. Att glömma att derivera funktionen. Integranden är g(x) = √(1 + [f'(x)]²), inte √(1 + [f(x)]²). För y = x² ger användning av f(x) i stället för f'(x) integralen av √(1 + x⁴), ett helt annat tal. Derivera alltid först.
3. Att använda fel båglängdsintegrand. Integranden måste vara √(1 + [f'(x)]²). Att tappa 1:an, glömma kvadraten eller ta kvadratroten av fel kvantitet ger alla fel svar. För en rät linje med lutningen 3 är integranden konstanten √(10) = 3.1622776602, inte 3 och inte √(1 + 3).
4. Att dela intervallet felaktigt. Bredden är Δx = (b - a) / n, med hela intervallängden överst. Att dividera med n + 1, eller använda b / n i stället för (b - a) / n, flyttar varje mittpunkt och förstör svaret. För [0, 2] med n = 4 är bredden 0.5, inte 0.4 och inte 0.6666666667.
5. Att använda fel antal mittpunkter. Det finns exakt n mittpunkter, en per delintervall, inte n + 1. Ändpunkterna är inte mittpunkter. Att räkna n + 1 sampel betyder att ett delintervall räknades två gånger eller att en ändpunkt samplades av misstag.
6. Att blanda ihop den ursprungliga kurvan med integrandgrafen. Mittpunktsrektanglarna sitter under integranden g(x), inte under den ursprungliga kurvan f(x). Att rita dem under f(x) och addera deras areor ger arean under f(x), vilket inte är båglängden. De två graferna är olika kurvor och måste visas separat.
7. Att avrunda mellanliggande värden. Att avrunda mittpunkterna, derivatvärdena eller integrandvärdena innan de adderas inför små fel som samlas över många rektanglar. Behåll full precision genom hela beräkningen och avrunda bara det slutliga visade resultatet. Kalkylatorn gör detta automatiskt.
8. Att anta att varje approximation är exakt. Mittpunktsregeln ger en approximation, inte det sanna värdet, utom i specialfall som räta linjer och parametriska kurvor med konstant hastighet. För parabeln är L₄ = 4.6267234873 enheter nära det sanna 4.6467837624 enheter, men det är inte lika med det. Rapportera alltid resultatet som en approximation och ange feluppskattningen vid sidan av det.
9. Att blanda ihop uppskattat fel med garanterade felgränser. Förfiningsskillnaden |M_2n - M_n| är en indikator på hur mycket svaret fortfarande rör sig, inte en garanti för att det sanna felet är mindre än den skillnaden. Den teoretiska gränsen ((b - a)³ / (24n²)) × max|g''(x)| är en garanti, men bara när andraderivatans maximum är känt och integranden är glatt.
10. Att ignorera funktionsdefinitionsmängdens begränsningar. En funktion som √(x) är odefinierad för negativt x, och ln(x) är odefinierad för x ≤ 0. Om en mittpunkt faller utanför definitionsmängden finns inte integrandvärdet, och approximationen är meningslös. Kalkylatorn rapporterar ett fel i det fallet i stället för att hitta på ett värde.
11. Att använda fel parametriska derivator. För en parametrisk kurva behövs båda derivatorna x'(t) och y'(t), och hastigheten är √([x'(t)]² + [y'(t)]²). Att använda bara en av dem, eller derivera med avseende på fel variabel, ger en fel hastighet och en fel längd. För kvartscirkeln är x'(t) = -5 sin(t) och y'(t) = 5 cos(t), och båda behövs för att få den konstanta hastigheten 5.
12. Att behandla en icke-glatt eller singulär funktion som om standardantagandena om fel automatiskt gäller. Felbeteendet 1 / n² och den teoretiska gränsen förutsätter en glatt integrand. En funktion med ett hörn, ett hopp eller en singularitet inne i intervallet kan konvergera mycket långsammare, och Richardson-extrapolationen kan misslyckas. Kontrollera konvergenstabellen: om en dubblering av n inte skär förändringen med ungefär fyra gäller inte standardantagandena för den funktionen.

En sista kontroll som fångar de flesta av dessa. Verifiera resultatet mot kordans längd och referensvärdet. Approximationen måste vara minst lika lång som kordan, den måste närma sig referensen när n växer, summan av tabellraderna måste vara lika med huvudresultatet, och för en rät linje eller en kurva med konstant hastighet måste den vara lika med det exakta värdet. Om någon av dessa kontroller misslyckas är ett av de tolv misstagen ovan orsaken.

Dessa sex exempel täcker kalkylatorns huvudsakliga lägen och funktionstyper. Varje exempel beräknas av samma mittpunktsmotor som driver verktyget, och varje resultat har verifierats oberoende mot en adaptiv Simpson-referens med hög precision.

Exempel 1: Parabel, y = x² på [0, 2], n = 4. Derivatan är f'(x) = 2x, bredden är Δx = 0.5, och mittpunkterna är 0.25, 0.75, 1.25 och 1.75. Integrandvärdena är 1.1180339887, 1.8027756377, 2.6925824036 och 3.6400549446, deras summa är 9.2534469747, och L₄ = 0.5 × 9.2534469747 = 4.6267234873 enheter. Referensen är 4.6467837624 enheter, så det absoluta felet är 0.0200602751 enheter och det relativa felet är 0.4317%. En dubblering till n = 8 ger 4.6417306581 enheter, och Richardson-uppskattningen från de två är 4.6467330483 enheter.

Exempel 2: Rät linje, y = 3x + 1 på [0, 4], n = 4. Derivatan är f'(x) = 3, så integranden är konstanten g(x) = √(10) = 3.1622776602. Bredden är Δx = 1, mittpunkterna är 0.5, 1.5, 2.5 och 3.5, och varje integrandvärde är 3.1622776602. Summan är 4 × 3.1622776602 = 12.6491106407, och L₄ = 1 × 12.6491106407 = 12.6491106407 enheter, vilket är exakt 4√(10). Mittpunktsregeln är exakt för varje n här, eftersom integranden är konstant och varje rektangel har rätt höjd.

Exempel 3: Horisontell linje, y = 2 på [1, 5], n = 4. Derivatan är f'(x) = 0, så integranden är konstanten g(x) = √(1 + 0) = 1. Bredden är Δx = 1, mittpunkterna är 1.5, 2.5, 3.5 och 4.5, och varje bidrag är 1 × 1 = 1. Totalen är L₄ = 4.0000000000 enheter, vilket är exakt det horisontella avståndet 5 - 1 = 4. En horisontell linje har ingen vertikal förändring, så dess båglängd är helt enkelt intervallängden.

Exempel 4: Parametrisk kvartscirkel, x(t) = 5 cos(t), y(t) = 5 sin(t) på [0, π/2], n = 4. Derivatorna är x'(t) = -5 sin(t) och y'(t) = 5 cos(t), och hastigheten är konstanten v(t) = 5. Parameterbredden är Δt = π/8 = 0.3926990817, mittpunktsparametrarna är π/16, 3π/16, 5π/16 och 7π/16, och varje bidrag är 0.3926990817 × 5 = 1.9634954085. Totalen är L₄ = 7.8539816340 enheter, exakt 5π/2. Hela cirkeln skulle vara 10π = 31.4159265359 enheter, och halvcirkeln skulle vara 5π = 15.7079632679 enheter.

Exempel 5: Kub, y = x³ på [0, 1], n = 4. Derivatan är f'(x) = 3x², bredden är Δx = 0.25, och mittpunkterna är 0.125, 0.375, 0.625 och 0.875. Integrandvärdena är 1.0010980300, 1.0853471867, 1.5405489332 och 2.5051217067, deras summa är 6.1321158566, och L₄ = 0.25 × 6.1321158566 = 1.5330289642 enheter. Referensen är 1.5478656547 enheter, så det absoluta felet är 0.0148366905 enheter och det relativa felet är 0.9585%. En dubblering till n = 8 ger 1.5441590533 enheter, och Richardson-uppskattningen är 1.5478690830 enheter.

Exempel 6: Trigonometrisk, y = sin(x) på [0, π], n = 8. Derivatan är f'(x) = cos(x), bredden är Δx = π/8 = 0.3926990817, och de åtta mittpunkterna är π/16, 3π/16, 5π/16, 7π/16, 9π/16, 11π/16, 13π/16 och 15π/16. Integrandvärdena är 1.4006926023, 1.3005159423, 1.1439660326, 1.0188524102, 1.0188524102, 1.1439660326, 1.3005159423 och 1.4006926023. Deras summa är 9.7280539749, och L₈ = 0.3926990817 × 9.7280539749 = 3.8201978627 enheter. Referensen är 3.8201977890 enheter, så det absoluta felet är bara 0.0000000736 enheter, ett relativt fel på 0.00000193%. Denna integral saknar elementär primitiv funktion, så mittpunktsregeln och referensen är det praktiska sättet att få svaret.

Klicka på valfri förinställningsknapp i kalkylatorn för att ladda ett av dessa exempel och se den fullständiga steg-för-steg-lösningen, livediagrammen, beräkningstabellen och konvergensjämförelsen för den inmatningen.

Båglängd med mittpunktsregeln dyker upp överallt där en böjd sträcka måste mätas från en formel snarare än från ett måttband. Samma beräkning tjänar studenter, ingenjörer och designers.

Vägar, banor och kablar. En böjd väg, en löparbanas kurva eller en kabel som hänger mellan två stöd har alla längder som följer en kurva. När kurvan modelleras av en funktion uppskattar mittpunktsregeln dess längd från modellen, och feluppskattningen säger hur mycket svaret kan litas på. En kabel modellerad av y = x² på [0, 2] behöver 4.6467837624 enheter kabel, inte de 4.4721359550 enheterna för det räta gapet.

Tillverkning och plåt. En böjd panel, ett böjt rör eller en rullad kant börjar som plant material av känd längd. Mittpunktsregeln omvandlar den böjda designen till den plana längd som ska skäras, och konvergenstabellen bekräftar längden till den precision verkstaden behöver. En kvartscirkelböj med radie 5 behöver 7.8539816340 enheter material längs kurvan.

Arkitektur. Valv, kupoler och böjda trappor dimensioneras efter sina bågar. En halvcirkelformad båge med radie 5 behöver 5π = 15.7079632679 enheter list längs kurvan. För icke-cirkulära bågar modellerade av parabler eller sinuskurvor ger mittpunktsregeln listlängden direkt från designekvationen.

Fysik och rörelse. Den sträcka en partikel tillryggalägger längs en böjd bana är båglängden av dess bana. En pendel, en projektil eller en satellit i en elliptisk omloppsbana rör sig alla längs kurvor vars längder spelar roll för energi, tid och bränsle. Mittpunktsregeln uppskattar de sträckorna från rörelseekvationerna.

Datorgrafik, animation och spel. Objekt som rör sig längs böjda banor, kamerarörelser och animerade figurer följer alla splines och parametriska kurvor. Sträckan längs banan avgör hastighet, tid och avstånd mellan bildrutor. Mittpunktsregeln mäter de banlängderna från kurvekvationerna, och det parametriska läget hanterar kurvor som vänder tillbaka.

Medicin och biologi. Blodkärl, ryggraden och de böjda ytorna hos organ mäts från skanningar och modeller. Båglängd längs ett kärls mittlinje hjälper till att planera stentar och katetrar, och mittpunktsregeln förvandlar en modellerad mittlinje till en längd med en feluppskattning.

Relaterade kalkylatorer på den här webbplatsen. För exakta båglängder av namngivna kurvor, prova parabelns båglängdskalkylator, ellipsens båglängdskalkylator, halvcirkelns båglängdskalkylator, cirkelns båglängdskalkylator och sektorns båglängdskalkylator. För steg-för-steg exakt integration, prova båglängdskalkylatorn med steg, den parametriska båglängdskalkylatorn, den polära båglängdskalkylatorn och 3D-båglängdskalkylatorn. För andra numeriska metoder, prova trapetsregelns kalkylator och Simpsons regels kalkylator.

Dessa tabeller samlar konvergensbeteendet hos de två huvudsakliga genomarbetade exemplen på ett ställe, så att du kan se hur snabbt mittpunktsregeln stabiliseras och hur Richardson-extrapolationen förbättrar svaret. Alla värden beräknas av mittpunktsmotorn, med referensen från adaptiv Simpson-integration.

Parabel y = x² på [0, 2], referens 4.6467837624 enheter.

Sinuskurva y = sin(x) på [0, π], referens 3.8201977890 enheter.

Så läser du tabellerna. Parabelns fel faller med en faktor på ungefär fyra vid varje dubblering, från 0.0200602751 vid n = 4 till 0.0050531044 vid n = 8 till 0.0012632302 vid n = 16, vilket är det förväntade beteendet 1 / n². Sinuskurvan konvergerar ännu snabbare, eftersom dess integrand är nästan platt nära intervallets mitt: vid n = 8 är felet redan 0.0000000736 enheter, och vid n = 16 stämmer approximationen överens med referensen till alla tio visade decimaler. Richardson-raderna visar samma idé i ett steg: att kombinera två approximationer slår var och en för sig, med ett fel på 0.0000507141 enheter för parabeln och 0.0000000245 enheter för sinuskurvan. Använd dessa tabeller som riktmärke: om din egen beräkning för ett av dessa exempel inte stämmer med tabellen, kontrollera derivatan, mittpunkterna och integranden innan du litar på svaret.