Fourierkalkylator

Introduktion till Fourierserier

Att dela upp en periodisk signal i sina enklaste beståndsdelar är en grundläggande metod inom signalbehandling, fysik och tillämpad matematik. Genom Fourieranalys kan en godtycklig periodisk eller styckvis kontinuerlig funktion dekomponeras till en oändlig summa av sinus- och cosinussvängningar. Denna matematiska representation kallas för en Fourierserie.

I praktiken arbetar man ofta med en trunkerad version av serien, en så kallad partiell summa, där man begränsar antalet harmoniska svängningar till ett ändligt antal N. Ju fler harmoniska komponenter som inkluderas i summan, desto närmare kommer den rekonstruerade funktionen att efterlikna ursprungssignalen. Fourierkalkylatorn är utformad för att beräkna dessa koefficienter numeriskt, visualisera approximationen och verifiera beräkningarnas stabilitet direkt i webbläsaren.

---

Matematisk konvention för Fourierkoefficienter

Det finns olika matematiska konventioner för hur en Fourierserie och dess koefficienter definieras i läroböcker. Denna Fourierkalkylator använder följande etablerade konvention för att uttrycka funktionen $f(x)$:

f(x) ≈ a₀ / 2 + Σₙ₌₁^N [ aₙ cos(nω₀ x) + bₙ sin(nω₀ x) ]

Här är vinkelfrekvensen definierad som ω₀ = 2π/T, där T är funktionens periodlängd beräknad som T = b - a.

En konsekvens av denna konvention är att samtliga koefficienter – a₀, aₙ och bₙ – utnyttjar samma integrationsfaktor, 2/T. Formlerna som tillämpas för beräkningen är:

• Konstanttermens koefficient: a₀ = 2 / T∫ₐᵇ f(x) dx
• Cosinuskoefficienter: aₙ = 2 / T∫ₐᵇ f(x) cos(nω₀ x) dx
• Sinuskoefficienter: bₙ = 2 / T∫ₐᵇ f(x) sin(nω₀ x) dx

Eftersom konstanttermen i serien skrivs som a₀/2, motsvarar denna term funktionens rena medelvärde över perioden.

---

Styckvis definierade funktioner i praktiken

Många tekniska och fysikaliska vågformer kan inte beskrivas med ett enda matematiskt uttryck över en hel period. Istället definieras de styckvis. Exempel på sådana vågformer är fyrkantsvåg, triangelvåg, sågtandsvåg och likriktad sinus.

Kalkylatorn tillåter definition av antingen en sammanhängande funktion över en period eller en funktion uppdelad i upp till fem separata delar. För varje del definieras unika gränser med fälten "från" och "till". När koefficienterna beräknas för en styckvis funktion integreras varje del separat över sitt eget intervall, varpå resultaten adderas för att erhålla de slutgiltiga värdena för aₙ och bₙ.

Godkända matematiska element

När du definierar uttrycket $f(x)$ kan följande matematiska operatorer, variabler och funktioner användas:

• Variabler och konstanter: x, pi, e
• Aritmetik: +, −, *, /, ^ (multiplikation kan skrivas explicit som 2*x eller implicit som 2x)
• Standardfunktioner: sin, cos, tan, exp, ln, log, sqrt, abs
• Extremvärden: min, max

---

Numerisk integration och stabilitetskontroll

Eftersom analytisk integration av komplexa eller styckvisa uttryck kan vara utomordentligt svår att lösa algebraiskt, använder verktyget numerisk integration via Simpsons sammansatta regel.

Simpsons delintervall och precision

Användaren kan ställa in parametern "Simpsons delintervall" till ett jämnt tal mellan 200 och 4000. Detta värde avgör hur fint integrationsintervallet delas upp under beräkningen. Fler delintervall ökar precisionen men kräver mer beräkningskapacitet.

Automatisk stabilitetskontroll

För att säkerställa att den numeriska approximationen är tillförlitlig utförs en automatisk stabilitetskontroll. Kalkylatorn upprepar hela integrationsprocessen med en fördubblad upplösning (dubbelt så många delintervall). Den maximala förändringen bland koefficienterna mäts och rapporteras som ‹delta›.

Beroende på resultatet av denna kontroll visas olika statusmeddelanden:

• Stabil integration: Om förändringen är försumbar visas texten "Beräknad. Koefficientkontrollen med dubblerad upplösning är stabil vid den visade precisionen."
• Instabil integration: Om skillnaden är märkbar visas varningen "Beräknad, men koefficienterna ändrades med ‹delta› när integrationsupplösningen dubblerades. Öka antalet delintervall eller undersök om funktionen har skarpa detaljer."

---

Gibbs fenomen vid diskontinuiteter

När man approximerar funktioner som innehåller skarpa språng – till exempel en fyrkantsvåg – uppstår ett karakteristiskt fenomen vid diskontinuiteterna. Detta kallas för Gibbs fenomen.

   f(x)
    ▲
1.0 ┼      ┌──────┐             --- Ursprunglig funktion f(x)
    │     *│      │*            *** Partiell Fouriersumma S_N(x)
    │   *  │      │  *              (med krusningar/översvängar)
0.5 ┼  *   │      │   *
    │ *    │      │    *
0.0 ┼──────┴──────┴──────► x

När en Fourierserie rekonstruerar ett abrupt språng kommer den partiella summan alltid att skjuta över målvärdet och uppvisa oscillerande krusningar precis före och efter språnget. Att lägga till fler harmoniska svängningar gör det påverkade området smalare men tar inte bort Gibbs karakteristiska översväng.

Om dina definierade funktionsdelar inte möts i ändpunkterna identifierar kalkylatorn detta och visar meddelandet: "Delarna har ett språng. Krusningar nära gränsen kan bero på Gibbs fenomen snarare än ett integrationsfel." Vid språnget närmar sig serien medelvärdet av vänster- och högergränsvärdena under vanliga konvergensvillkor.

Gränssnittets in- och utmatningar

För att genomföra en beräkning ställer du in parametrarna i kalkylatorns gränssnitt. Nedan beskrivs de centrala fälten och de resultat som genereras.

Inmatningsparametrar

• Funktion över en period / Funktionsdelar: Välj om du vill mata in en enskild funktion eller dela upp den i upp till fem segment med "Lägg till del" och "Ta bort del ‹n›".
• Periodlängd: Definieras via "Startpunkt a" och "Slut b". Trigonometriska argument förutsätter radianer, och gränserna använder samma enhet som x.
• Harmoniska N: Ett heltal från 1 till 30 som anger hur många termer som ska beräknas till den partiella summan.
• Visade decimaler: Anger hur många decimaler som ska redovisas i tabeller och formler.
• Beräkna Sₙ(x) vid: Ett valfritt inmatningsfält för att utvärdera den partiella summan i en specifik punkt x.

Utmatade resultat och tabeller

När du klickar på "Beräkna serie" genereras följande sektioner:

1. Fourier-rekonstruktion: Visar den resulterande formeln under rubriken "Partiell summa" med en tillhörande knapp för att "Kopiera serie". Här visas även grafen "Originalfunktion och partiell Fouriersumma" med teckenförklaringar för "f(x)" och "Sₙ(x)".
2. Koefficienter: En tabell med kolumnerna "aₙ", "bₙ" och "Amplitud" för varje harmonisk komponent upp till N.
3. Formel och insättning: Detaljerad redovisning av beräkningsstegen:
   • Periodens egenskaper: "T = ‹period›; ω₀ = 2π/T = ‹omega› radianer per x-enhet."
   • Konstanttermen: "a₀ = (2/T)∫f(x)dx = ‹a0›; konstanttermen är a₀/2 = ‹constant›."
   • Integrationsmetod: "För n = 1…‹terms›, integrera varje del separat och addera resultaten för aₙ och bₙ."
   • Precisionsrapport: "Simpsons sammansatta regel: ‹base› delintervall, kontrollerat igen med ‹refined›; maximal koefficientförändring ‹delta›."
   • Om en specifik punkt utvärderats visas även "Sₙ(‹x›) = ‹value›".

---

Vanliga frågor

Är koefficienterna exakta?

Nej. De är numeriska integraler beräknade med Simpsons sammansatta regel. Kalkylatorn upprepar beräkningen med dubbelt så hög upplösning och rapporterar den största förändringen bland koefficienterna. En liten förändring ger stöd för de visade siffrorna men är inte en formell felgräns.

Vilken konvention för Fourierkoefficienter använder denna kalkylator?

Den skriver f(x) ≈ a₀/2 + Σ[aₙ cos(nω₀x) + bₙ sin(nω₀x)], med ω₀ = 2π/T. Följaktligen använder a₀, aₙ och bₙ alla faktorn 2/T. Vissa böcker lägger istället medelvärdet direkt i a₀; jämför den visade formeln innan du matchar koefficienter.

Varför uppstår det krusningar i grafen nära ett språng?

En ändlig Fouriersumma oscillerar nära en diskontinuitet. Att lägga till fler harmoniska svängningar gör det påverkade området smalare men tar inte bort Gibbs karakteristiska översväng. Vid språnget närmar sig serien medelvärdet av vänster- och högergränsvärdena under vanliga konvergensvillkor.

---

Integritet och lokal körning

Dina matematiska uttryck, parametrar och beräkningsresultat är helt privata. Alla beräkningar och visualiseringar sker lokalt i din egen webbläsare. Inga data eller formler laddas någonsin upp till någon extern server.