Zustand
Beim mechanischen Oszillator genügt der Zustandsvektor aus Auslenkung und Geschwindigkeit:
x⃗(t) = [ x(t), v(t) ]ᵀ
M01 · K01
Schwingung als
Schwingung als
dynamisches System
Eine Schwingung ist kein gezeichneter Sinus, sondern die zeitliche Entwicklung eines Zustands. Bereits ein Masse-Feder-System bildet einen geschlossenen Rückkopplungskreis: Aus der Lage entsteht Kraft, aus Kraft Beschleunigung, daraus Geschwindigkeit und wieder eine neue Lage.
01 · Systemidee
Ein System „erinnert“ sich über seinen Zustand
Für die zukünftige Entwicklung reichen Eingangsgrößen allein nicht. Entscheidend ist, in welchem Zustand sich das System gerade befindet.
Rückkopplung
Die Feder erzeugt eine Kraft entgegen der Auslenkung. Der aktuelle Zustand verändert damit seine eigene weitere Entwicklung.
FR = −k · x
Dynamik
Newton koppelt Kraft an Beschleunigung. Damit entsteht ein System erster Ordnung in zwei Zustandsgrößen.
dx/dt = v
dv/dt = −(k/m) · x
02 · Rückkopplungsstruktur
Warum das System überhaupt schwingt
Die Rückstellkraft dreht die Bewegungsrichtung nicht sofort um. Zuerst ändert sie die Geschwindigkeit. Genau diese gespeicherte Bewegung erzeugt das fortgesetzte Überschwingen.
1Auslenkung xgeometrischer Zustand
2Rückstellkraft −kxnegative Rückkopplung
3Beschleunigung aF / m
4Geschwindigkeit vintegrierte Beschleunigung
Die beiden Integrationen a → v und v → x liefern dem System sein Gedächtnis. Ohne Anfangszustand gibt es daher keine vollständige Vorhersage.
03 · Interactive Lab
State-Feedback Oscillator
Verändere Masse, Federsteifigkeit und Anfangszustand. Die drei Darstellungen zeigen dieselbe Simulation gleichzeitig als Mechanik, Zeitverlauf und Zustandsbahn.
0.55 kg
580 N/m
−1.2+1.2 m
−3+3 m/s
Mechanisches SystemMasse · Feder · Zustand
läuft
x0.750m
v0.000m/s
a−12.000m/s²
E6.750J
x(t)v(t)
Die Achsen sind mit A bzw. Aω₀ normiert. Im ungedämpften idealen Modell bleibt die Gesamtenergie konstant; die Zustandsbahn schließt sich deshalb.
04 · Modell lesen
Was die Simulation tatsächlich berechnet
ẋ = v
v̇ = −(k/m)x
Zu jedem Zeitpunkt werden nur der aktuelle Zustand und die Systemparameter benötigt. Die nächste Zustandsänderung folgt daraus deterministisch.
Ekin = ½mv²
Epot = ½kx²
Die Bewegung entsteht durch fortlaufenden Austausch zwischen kinetischer und elastischer Energie. Im idealen K01-Modell wird keine Energie dissipiert.
- Punktmasse
- lineare Feder
- keine Reibung
- keine äußere Anregung
- eine räumliche Freiheitsrichtung
Diese bewusste Reduktion isoliert die Dynamik. Dämpfung und externe Kräfte werden später als zusätzliche Systemkopplungen eingeführt.
05 · Interpretation
Vier Beobachtungen, die wir mitnehmen
Eine Position allein genügt nicht: gleiche Auslenkung bei anderer Geschwindigkeit bedeutet eine andere Zukunft.
Die aktuelle Lage verändert über die Kraft die eigene zeitliche Entwicklung.
Beschleunigung verändert Geschwindigkeit; Geschwindigkeit verändert Position. Der Zustand springt im idealen Modell nicht.
Reale Systeme verlieren Energie und werden angeregt. Diese zusätzlichen Terme erweitern später dieselbe Grundstruktur.
M01 · K02
Harmonischer
Harmonischer
Oszillator
Das dynamische System aus K01 besitzt im idealen linearen Fall eine geschlossene Lösung. Aus Masse und Federsteifigkeit folgt eine charakteristische Eigenkreisfrequenz; der Anfangszustand legt fest, welche konkrete Bahn auf dieser Systemdynamik realisiert wird.
01 · Bewegungsgleichung
Aus zwei Zustandsableitungen wird eine Gleichung zweiter Ordnung
Das Zustandsmodell aus K01 lässt sich zu einer einzigen Differentialgleichung zusammenfassen. Diese Form macht die Eigenstruktur des Systems unmittelbar sichtbar.
Newton + Hooke
m ẍ = −kx
m ẍ + kx = 0
Die lineare Feder liefert eine Kraft proportional zur Auslenkung. Ohne Dämpfung und äußere Kraft ist dies das vollständige Modell.
Normierte Form
ẍ + ω0²x = 0
ω0 = √(k/m)
Die gesamte Zeitstruktur des idealen Systems steckt in einer einzigen Systemgröße: der Eigenkreisfrequenz ω₀.
Charakteristische Lösung
r² + ω0² = 0
r = ±iω0
Die rein imaginären Eigenwerte bedeuten: keine exponentielle Verstärkung und kein Abklingen — der Zustand rotiert periodisch im Phasenraum.
02 · Geschlossene Lösung
Der Anfangszustand bestimmt die konkrete Schwingung
Masse und Feder legen die Systemfrequenz fest. Auslenkung und Geschwindigkeit bei t = 0 bestimmen, wie diese Eigenbewegung gestartet wird.
x(t) = x0 cos(ω0t) + v0ω0 sin(ω0t)
Diese Form ist besonders nützlich, weil sie direkt mit dem Zustand aus K01 arbeitet. Für jeden zulässigen Anfangszustand liefert sie die vollständige zukünftige Bewegung.
ω0 = √(k/m)
f0 = ω0 / 2π
k ↑ macht das System schneller, m ↑ macht es träger.
x(t) = A cos(ω0t + φ)
Amplitude und Phase sind eine alternative Kodierung desselben Anfangszustands. K03 untersucht diese Größen systematisch.
03 · Interactive Lab
Analytical Harmonic Oscillator
Verändere die Systemparameter und den Anfangszustand. Die analytische Lösung wird direkt berechnet; zusätzlich erzeugt derselbe Zustand eine numerische RK4-Referenz wie in K01. Beide müssen im idealen Modell übereinanderliegen.
analytische Lösung
RK4 aus K01
04 · Systemparameter lesen
Eigenfrequenz gehört zum System — nicht zur Anregung
Bei gleicher Auslenkung wächst die Rückstellkraft. Der Zustand wird schneller zur Gleichgewichtslage beschleunigt.
ω0 ∝ √k
Dieselbe Kraft erzeugt weniger Beschleunigung. Die Zustandsentwicklung wird langsamer.
ω0 ∝ 1/√m
Anfangsbedingungen verändern Amplitude und Phasenlage, aber nicht die Eigenkreisfrequenz des linearen Systems.
ω0 = √(k/m)
05 · Modellgrenze
Warum die ideale Sinusschwingung in der Realität selten ewig bleibt
Hookes Gesetz wird als exakt linear angenommen. Reale Federn können bei großen Auslenkungen nichtlinear werden.
Ohne Reibung bleibt die Energie konstant. Dämpfung wird in K05 als zusätzlicher Term eingeführt.
Das System schwingt nur aus seinem Anfangszustand heraus. Erzwungene Bewegung folgt in K06.
Ein Freiheitsgrad besitzt hier genau eine Eigenfrequenz. Gekoppelte Systeme entwickeln mehrere Moden — Thema K09.
M01 · K03
Amplitude, Periode,
Amplitude, Periode,
Frequenz & Phase
K02 liefert die geschlossene Bewegung. K03 macht daraus präzise Messgrößen: Wie groß ist die Schwingung, wie oft wiederholt sie sich und an welcher Stelle ihres Zyklus befindet sie sich gerade?
01 · Kanonische Beschreibung
Eine harmonische Schwingung lässt sich mit vier Größen lesen
Für eine gegebene Gleichgewichtslage beschreibt die Cosinusform die vollständige Bewegung kompakt.
x(t) = A cos(ωt + φ)
A Amplitude
T Periodendauer
f Frequenz
φ Anfangsphase
Amplitude
Maximaler Betrag der Auslenkung relativ zur Gleichgewichtslage. Sie legt die räumliche Skala der Bewegung fest.
A = max |x(t)|
Einheit hier: m
Periodendauer
Kleinste positive Zeitspanne, nach der sich der gesamte Zustand wiederholt.
x(t + T) = x(t)
Einheit: s
Frequenz
Anzahl vollständiger Zyklen pro Sekunde. Sie ist der Kehrwert der Periodendauer.
f = 1 / T
Einheit: Hz = s⁻¹
Phase
Winkelkoordinate des Zustands innerhalb des periodischen Zyklus. Phasen unterscheiden sich nur modulo 2π.
θ(t) = ωt + φ
Einheit: rad oder °
02 · Zeit und Winkel
Eine Periode entspricht genau einer vollständigen Phasendrehung
f = 1/T
Verdoppelt sich die Frequenz, halbiert sich die Zeit für einen Zyklus.
ω = 2πf
Die Kreisfrequenz misst Phasenfortschritt in Radiant pro Sekunde.
θ = ωt + φ
Nach einer Periodendauer gilt Δθ = 2π: geometrisch ist der Zyklus geschlossen.
03 · Interactive Lab
Phase & Period Lab
Eine Zeitfunktion und ein rotierender Phasenzeiger zeigen dieselbe Bewegung in zwei Koordinatensystemen. Amplitude, Frequenz und Anfangsphase bleiben in beiden Darstellungen geometrisch gekoppelt.
Die Breite des Fensters entspricht drei Perioden. Dadurch bleibt die geometrische Bedeutung von T unabhängig von der eingestellten Frequenz sichtbar.
θ = ωt + φ
0.524 rad
04 · Anfangszustand ↔ Phase
Phase ist keine zusätzliche Dynamik — sie kodiert den Startpunkt
Bei festem A und ω legt φ fest, wo auf der geschlossenen Zustandsbahn die Beobachtung bei t = 0 beginnt.
x0 = A cos φ
v0 = −Aω sin φ
Damit ist die Verbindung zu K02 vollständig: Anfangsauslenkung und Anfangsgeschwindigkeit lassen sich eindeutig in Amplitude und Phase umrechnen, solange ω bekannt ist.
A = √(x0² + (v0/ω)²)
φ ≡ φ + 2πn
Winkel, die sich um ganze Umdrehungen unterscheiden, beschreiben denselben Punkt im Zyklus.
05 · Präzisionsbegriffe
Was leicht verwechselt wird
f zählt Zyklen pro Sekunde, ω zählt Radiant pro Sekunde. Numerisch unterscheiden sie sich um den Faktor 2π.
Phase ist dimensionslos bzw. ein Winkel. Zeit wird erst durch ω in Phasenfortschritt übersetzt.
A ist die maximale Auslenkung; x(t) liegt innerhalb des Intervalls −A … +A.
Der Nullpunkt der Phase kann gewählt werden. Physikalisch besonders relevant werden später Phasendifferenzen zwischen Signalen.
M01 · K04
Energy
Energy
exchange
Im idealen harmonischen Oszillator verschwindet Energie nicht und wird nicht von außen zugeführt. Sie wandert periodisch zwischen zwei Speichern: kinetischer Energie der Masse und potentieller Energie der Feder.
01 · Energiespeicher
Feder und Masse speichern Energie in unterschiedlichen Zustandsgrößen
Die potentielle Energie hängt ausschließlich von der Auslenkung ab, die kinetische ausschließlich von der Geschwindigkeit.
Potentielle Federenergie
Epot = ½kx²
Maximal an den Umkehrpunkten x = ±A. Dort steht die Masse momentan still.
Kinetische Energie
Ekin = ½mv²
Maximal beim Durchgang durch die Gleichgewichtslage x = 0. Dort ist der Betrag der Geschwindigkeit maximal.
Gesamtenergie
E = Ekin + Epot
Für den ungedämpften, nicht angeregten Oszillator bleibt E konstant.
02 · Erhaltung
Die Bewegungsgleichung enthält die Energieerhaltung bereits implizit
Bewegungsgleichung
m ẍ + kx = 0
mit v = ẋ multiplizieren
m ẍẋ + kxẋ = 0
als Zeitableitung lesen
d/dt (½mv² + ½kx²) = 0
Interpretation:
Die Summe der beiden Energiespeicher ändert sich nicht. Energie wird nur intern zwischen Feder und Masse umverteilt.
03 · Interactive Lab
Energy Exchange Lab
Verändere Masse, Federsteifigkeit und Anfangsamplitude. Das Lab zeigt Bewegung und Energieverteilung synchron — einschließlich des charakteristischen Energieaustauschs mit doppelter Frequenz.
Mechanischer Zustandx(t), v(t)
läuft
x0.800m
v0.000m/s
a−12.800m/s²
EnergiespeicherVerteilung von E
100% pot.
Epot
Ekin
Die beiden Anteile ergänzen sich jederzeit exakt zur Gesamtenergie.
Eges
Epot
Ekin
04 · Energiephase
Die Energieverteilung oszilliert doppelt so schnell wie die Auslenkung
Epot = ½kA² cos²(ωt)
Ekin = ½kA² sin²(ωt)
cos²θ = ½(1 + cos 2θ)
sin²θ = ½(1 − cos 2θ)
Innerhalb einer mechanischen Periode T besitzt jede Energiekomponente zwei Maxima.
05 · Systemische Lesart
Der Oszillator ist ein rückgekoppeltes Speichernetzwerk
Zustandsgröße v, Energie ½mv².
↔
periodischer
Energiefluss
Energiefluss
Zustandsgröße x, Energie ½kx².
Diese Sicht ist später direkt übertragbar: auf elektrische LC-Schwingkreise, akustische Resonatoren und elektromagnetische Feldspeicher.
M01 · K05
Damping
Reale Oszillatoren verlieren mechanische Energie. Ein linearer Dämpfer ergänzt die Rückstellkraft um eine geschwindigkeitsabhängige Gegenkraft. Dadurch ändert sich nicht nur die Amplitude, sondern die gesamte Geometrie der Zustandsbahn.
01 · Erweiterte Bewegungsgleichung
Der Dämpfer koppelt die Gegenkraft an die Geschwindigkeit
Für lineare viskose Dämpfung wirkt die Dämpferkraft stets entgegen der momentanen Bewegungsrichtung.
Kräfte
FFeder = −kx
FDämpfer = −cv
Der Dämpfer reagiert nicht auf die Position, sondern auf die Änderungsrate des Zustands.
Bewegungsgleichung
m ẍ + c ẋ + kx = 0
Mit c = 0 fällt das Modell exakt auf den idealen Oszillator aus K01–K04 zurück.
Dimensionslose Kennzahl
ζ = c / (2√(km))
ckrit = 2√(km)
Das Dämpfungsmaß ζ entscheidet über das qualitative Bewegungsregime.
02 · Drei Regime
Eine Gleichung — drei qualitativ verschiedene Antworten
Unterdämpfung
ωd = ω0√(1−ζ²)
x(t) ∝ e−ζω₀t cos(ωdt + φ)
Das System schwingt weiter, aber seine Hüllkurve fällt exponentiell ab.
Kritische Dämpfung
x(t) = (C₁ + C₂t)e−ω₀t
Doppelter reeller Eigenwert. Für den Referenzfall „Auslenkung aus Ruhe“ ergibt sich die schnellste aperiodische Rückkehr ohne Überschwingen.
Überdämpfung
r₁,₂ = −ω₀(ζ ∓ √(ζ²−1))
Zwei reelle Zerfallsmoden; die Rückkehr ist langsamer und nicht schwingend.
03 · Energiefluss
Dämpfung entfernt mechanische Energie irreversibel aus dem Oszillator
E = ½mv² + ½kx²
dE/dt = mv v̇ + kx ẋ
dE/dt = v(mẍ + kx)
dE/dt = −cv² ≤ 0
Die Energie nimmt nur dann momentan nicht ab, wenn v = 0. Über einen vollständigen Bewegungsabschnitt ist der Nettofluss bei c > 0 stets aus dem mechanischen System heraus gerichtet.
04 · Interactive Lab
Damping Regime Lab
Verändere m, k, c und den Anfangszustand. Das Lab klassifiziert das Regime automatisch und zeigt dieselbe Dynamik synchron im Zeitbereich, Phasenraum und Energieverlauf.
Exakte Hüllkurve des unterdämpften Falls.
Spiralbahn zur Ruhelage
05 · Zeitkonstanten & Hüllkurve
Unterdämpfung trennt schnelle Schwingung und langsamen Amplitudenzerfall
α = c/(2m) = ζω₀
Die Hüllkurve fällt wie e−αt.
τ = 1/α = 2m/c
Nach τ ist die Amplitudenhüllkurve auf 1/e ihres Anfangswerts gefallen.
ωd = ω₀√(1−ζ²)
Nur für ζ < 1 existiert eine oszillatorische Eigenbewegung.
06 · Cross-domain
Dämpfung ist eine universelle Dissipationsstruktur
Stoßdämpfer kontrollieren Aufbau- und Radschwingungen.
Viskose, thermische und Strahlungsverluste begrenzen die Güte.
Ohmscher Widerstand spielt mathematisch dieselbe dissipative Rolle.
Reibung und turbulente Prozesse entziehen geordneten Bewegungen Energie.
M01 · K06
Forced
Forced
oscillation
Bis K05 entwickelte sich der Oszillator ausschließlich aus seinem Anfangszustand. Jetzt koppeln wir eine äußere periodische Quelle an das System. Die Bewegung besteht dann aus einem abklingenden Einschwinganteil und einer dauerhaft erzwungenen stationären Antwort.
01 · Anregung
Eine periodische Kraft öffnet das dissipative System für Energiezufuhr
Der Dämpfer entzieht weiterhin Energie. Die äußere Quelle kann diese Verluste fortlaufend ersetzen.
äußere Kraft
F(t) = F₀ cos(Ωt)
Amplitude F₀ und Kreisfrequenz Ω sind Eigenschaften der Quelle — nicht des Oszillators.
erzwungene Bewegung
m ẍ + c ẋ + kx = F₀ cos(Ωt)
Die rechte Seite macht aus dem autonomen System ein angeregtes System.
zwei Frequenzen unterscheiden
ω₀ = √(k/m)
Ω = Frequenz der Quelle
ω₀ gehört zum System. Ω wird von außen aufgeprägt.
02 · Zerlegung der Antwort
Transienter Anteil plus stationäre erzwungene Bewegung
x(t) = xtr(t) + xss(t)
Löst die homogene Gleichung aus K05 und verschwindet bei c > 0 mit der Zeit.
Bleibt nach dem Einschwingen erhalten und oszilliert exakt mit der Quellenfrequenz Ω.
03 · Stationäre Lösung
Amplitude und Phasenlage entstehen aus der komplexen Systemantwort
Für eine harmonische Quelle ist die stationäre Antwort wieder harmonisch — jedoch im Allgemeinen mit anderer Amplitude und zeitlicher Verzögerung.
xss(t) = X cos(Ωt − δ)
Die Quelle bestimmt Ω. Das System bestimmt daraus X und δ.
X = F₀ / √((k − mΩ²)² + (cΩ)²)
δ = atan2(cΩ, k − mΩ²)
Mit dieser Form liegt δ konsistent zwischen 0 und π.
04 · Interactive Lab
Forced Response Lab
Die numerische RK4-Antwort startet aus Ruhe. Gleichzeitig wird die analytische stationäre Lösung berechnet. So wird sichtbar, wie der Einschwinganteil verschwindet und die Bewegung an die Quelle gekoppelt bleibt.
RK4: x(t)
analytisch: xss(t)
Im Phasorbild ist die Kraft Referenz bei 0°. Die stationäre Auslenkung folgt um δ verzögert.
05 · Energie- und Leistungsbilanz
Im stationären Zustand ersetzt die Quelle im Mittel genau die dissipierte Energie
E = ½mv² + ½kx²
dE/dt = F(t)v − cv²
⟨Fv⟩ = ⟨cv²⟩
Über einen vollständigen stationären Zyklus bleibt die mittlere gespeicherte mechanische Energie konstant.
06 · Systemische Lesart
Forced oscillation ist ein Quelle–Übertragungsstrecke–Antwort-System
Amplitude und Frequenz der externen Quelle.
→
Speicher, Rückkopplung und Dissipation formen die Übertragung.
→
Stationäre Antwortamplitude und Phasenverzug.
M01 · K07
Resonance
In K06 betrachteten wir einzelne Anregungsfrequenzen. Jetzt wird Ω systematisch variiert. Dadurch erscheint die Frequenzantwort als Eigenschaft des gesamten Quelle–System–Antwort-Netzes: Dämpfung verschiebt und begrenzt das Maximum der Verschiebungsamplitude, während die maximale mittlere Leistungsaufnahme an einer anderen, präzise definierten Stelle liegt.
01 · Frequenzsweep
Die stationäre Antwort wird zur Übertragungsfunktion
Die absolute Antwort aus K06 lässt sich mit der statischen Auslenkung Xstat = F₀/k normieren. Übrig bleibt eine dimensionslose Vergrößerungsfunktion, die nur noch vom Frequenzverhältnis r und der Dämpfungszahl ζ abhängt.
Frequenzverhältnis
r = Ω / ω0
r = 1 bedeutet: Die Quellenfrequenz entspricht der ungedämpften Eigenkreisfrequenz des Systems.
statische Referenz
Xstat = F₀ / k
Diese Auslenkung würde dieselbe Kraft als konstante Last erzeugen. Sie liefert die natürliche Amplitudenskala.
dynamische Vergrößerung
M(r) = X / Xstat = 1 / √((1 − r²)² + (2ζr)²)
Die Resonanzkurve ist damit unabhängig von F₀; die Kraft skaliert die absolute Antwort linear.
02 · Lage des Amplitudenmaximums
Bei Dämpfung liegt der Verschiebungspeak unterhalb von ω₀
Das oft verwendete „Resonanz bei Ω = ω₀“ ist für die Verschiebungsamplitude nur im schwach gedämpften Grenzfall exakt genug.
rpeak = √(1 − 2ζ²)
Damit gilt Ωpeak = ω₀√(1 − 2ζ²). Für kleine ζ nähert sich der Peak r = 1.
ζ < 1 / √2
Nur dann besitzt die Verschiebungsantwort ein inneres Maximum bei Ω > 0. Oberhalb dieser Grenze fällt M(r) von r = 0 an monoton.
Mpeak = 1 / (2ζ√(1 − ζ²))
Im idealisierten Grenzfall ζ → 0 wächst der Peak ohne Begrenzung. Reale Verluste verhindern diese mathematische Singularität.
03 · Interactive Lab
Resonance Sweep Lab
Verändere Systemparameter, Dämpfung und aktuelle Anregungsfrequenz. Die obere Kurve zeigt die Verschiebungsvergrößerung, die untere die normierte mittlere Leistungsaufnahme. So werden Amplitudenresonanz und maximale Energieaufnahme direkt vergleichbar.
VerschiebungsantwortM(r) = X / Xstat
Die gestrichelte Referenz markiert r = 1; der Peakmarker zeigt das tatsächliche Maximum der Verschiebungsamplitude.
Verschiebungsvergrößerung M(r)
aktuelle Anregung
Energieaufnahmenormierte mittlere Leistung
Für viskose Dämpfung liegt das Leistungsmaximum exakt bei r = 1 — unabhängig von der Verschiebungs-Peakverschiebung.
normierte mittlere Leistung
aktuelle Anregung
Aktueller PunktQuelle → System → Antwort
Ω3.000rad/s
Xstat0.333m
P̄1.000W
P̄/P̄max0.600–
04 · Energieperspektive
Amplitudenresonanz und maximale Leistungsaufnahme sind verschiedene Fragen
P̄ = ½ c Ω² X²
P̄ / P̄max = (2ζr)² / ((1 − r²)² + (2ζr)²)
Im stationären Zustand entspricht dieser Mittelwert der mittleren Leistung, welche die äußere Quelle zuführt.
Ω = ω0
Bei viskoser Dämpfung ist die Energieübertragung von der Quelle in das dissipative System genau bei r = 1 maximal.
Die Verschiebung maximiert eine kinematische Observable. Die Leistung enthält zusätzlich den Faktor Ω² und hängt über v = ΩX direkt von der Geschwindigkeit ab. Daher verschieben sich die Optima gegeneinander.
05 · Systemische Lesart
Resonanz ist selektive Energie- und Zustandsverstärkung eines dynamischen Systems
Masse und Feder tauschen kinetische und potenzielle Energie. Diese Speicherstruktur erzeugt die charakteristische Eigenfrequenz.
Dämpfung begrenzt die Akkumulation, senkt den Peak und kann einen ausgeprägten Verschiebungs-Resonanzpeak vollständig verschwinden lassen.
Die Quelle liefert Energie. Entscheidend ist nicht nur ihre Stärke, sondern wie ihre Zeitstruktur zur Systemdynamik passt.
06 · Transfer in andere Domänen
Die gleiche Resonanzstruktur erscheint weit über den Masse-Feder-Oszillator hinaus
Unwucht, Straßenanregung und elastische Strukturen erzeugen frequenzselektive Schwingungsantworten.
Luftvolumen, Halsmasse und Verluste bilden mechanisch analoge Speicher- und Dissipationsstrukturen.
→ M05K07 · Resonators and pipesInduktivität und Kapazität speichern Energie, Widerstand dissipiert sie; die mathematische Frequenzantwort ist strukturell äquivalent.
M01 · K08
Quality factor and bandwidth
K07 zeigte, wo Resonanzmaxima liegen. K08 misst nun, wie selektiv das System ist. Güte Q, Dämpfungszahl ζ, Polabklingrate σ und Bandbreite Δω beschreiben dieselbe zweite Ordnung aus vier Blickwinkeln: Frequenzselektion, Verlust, Zeitabklingen und Systempole.
01 · Güte als Systemparameter
Q ist die dimensionslose Form der Dämpfung
Für den linearen viskos gedämpften Oszillator ist Q direkt im normierten Nenner der zweiten Ordnung enthalten. In dieser Systemdefinition ist die Beziehung zu ζ exakt; erst die Gleichsetzung mit beliebigen gemessenen −3-dB-Breiten benötigt zusätzliche Bedingungen.
Güteparameter
Q = 1 / (2ζ)
Mit ζ = c/(2√(km)) folgt ebenso Q = √(km)/c = mω0/c.
Polabklingrate
σ = ζω0 = c/(2m)
Im unterkritischen freien System trägt die Amplitudenhülle den Faktor e−σt. Kleine σ bedeuten langes Nachschwingen.
gleiche Struktur
Q = ω0 / (2σ)
Q vergleicht die charakteristische Schwingungsrate mit der Verlust- bzw. Abklingrate.
02 · Halbwertsbreite der Leistungsaufnahme
Für die dissipierte Leistung lässt sich die Bandbreite exakt bestimmen
Wir verwenden dieselbe normierte mittlere Leistungsaufnahme wie in K07. Die −3-dB-Punkte liegen dort, wo nur noch die Hälfte der maximalen mittleren Leistung aufgenommen wird.
r− = √(1 + ζ²) − ζ
r+ = √(1 + ζ²) + ζ
Das geometrische Mittel bleibt exakt √(r−r+) = 1. Die Eigenfrequenz liegt also logarithmisch in der Mitte.
Δr = r+ − r− = 2ζ
Daraus folgt Δω = 2ζω0 = c/m.
Q = ω0 / Δω = 1 / Δr
Für diese explizit definierte Leistungs-Halbwertsbreite ist die Beziehung exakt, nicht nur eine Klein-Dämpfungs-Näherung.
03 · Interactive Lab
Q & Bandwidth Lab
Ändere Masse, Steifigkeit und Dämpfungszahl. Frequenz- und Zeitbereich reagieren gemeinsam: eine schmalere −3-dB-Breite geht unmittelbar mit langsamerem exponentiellem Abklingen einher.
FrequenzbereichLeistungsaufnahme und −3-dB-Bandbreite
r− und r+ markieren exakt P̄/P̄max = 0,5. Die horizontale Linie entspricht −3,010 dB.
normierte Leistung
Halbwertspunkte
aktueller Arbeitspunkt
Zeitbereichfreie Amplitudenhülle e−σt
Die x-Achse zählt ungedämpfte Eigenperioden T₀. Der 1/e-Punkt liegt bei N = Q/π.
04 · Was −3 dB tatsächlich bedeutet
Bandbreite ist ohne Angabe der Observable nicht eindeutig
10 log₁₀(P/Pmax) = −3.010 dB
Das entspricht P/Pmax = 1/2. Für die in K07 definierte dissipierte Leistung liefert dies die exakten Punkte r−, r+.
20 log₁₀(A/Amax) = −3.010 dB
Hier ist A/Amax = 1/√2. Bei der Verschiebungsantwort liegen diese Punkte im Allgemeinen nicht exakt an denselben Frequenzen wie die Leistungs-Halbwertspunkte.
Q ≈ ωpeak / Δωx,−3dB
Für ζ ≪ 1 verschmelzen ωpeak, ω₀ und die verschiedenen Bandbreitendefinitionen praktisch. Bei stärkerer Dämpfung muss die Observable genannt werden.
05 · Zeit–Frequenz-Dualität
Breite im Frequenzraum und Lebensdauer im Zeitraum sind dieselbe Verlustinformation
Langsames Ring-down, schmale Leistungsbandbreite und starke Frequenzselektion.
Schnelles Abklingen, breite Leistungsbandbreite und geringe Frequenzselektion.
Für die Leistungs-Halbwertsbreite gilt damit zugleich Q = ω₀/(2σ) = ω₀/Δω.
06 · Transfer in andere Domänen
Q misst überall das Verhältnis von Speicherung zu Verlust
Hohe Q-Werte erzeugen schmale akustische Resonanzen und langes Ausschwingen; Verluste an Wänden und Öffnungen verbreitern sie.
Widerstände bestimmen Verluste, L und C die Speicherstruktur. Q erscheint direkt als Selektivitätsparameter zweiter Ordnung.
Modalgüte und Dämpfungsmaß verbinden Ring-down-Messungen mit Peaks im Frequenzgang.
M01 · K09
Gekoppelte
Gekoppelte
Oszillatoren
Zwei Freiheitsgrade verändern die Logik des Systems grundlegend: Die physikalischen Koordinaten sind gekoppelt, doch geeignete Linearkombinationen entkoppeln die Dynamik wieder. Diese Normalmoden sind die diskrete Vorstufe von Eigenmoden, stehenden Wellen und schließlich kontinuierlicher Wellenausbreitung.
01 · Zwei Freiheitsgrade
Eine Kopplungsfeder macht aus zwei Oszillatoren ein gemeinsames System
Beide Massen besitzen dieselbe Masse m und dieselbe äußere Federsteifigkeit k. Die mittlere Feder mit Steifigkeit κ koppelt die Auslenkungen x₁ und x₂.
Gekoppelte Bewegungsgleichungen
m ẍ₁ + (k + κ)x₁ − κx₂ = 0
m ẍ₂ + (k + κ)x₂ − κx₁ = 0
Die Beschleunigung einer Masse hängt nun auch von der Position der anderen ab. x₁ und x₂ sind keine unabhängigen Ein-Freiheitsgrad-Systeme mehr.
Matrixform
M ẍ + Kx = 0
M = m I
K = [[k+κ, −κ], [−κ, k+κ]]
Die Kopplung erscheint als Nebendiagonale der Steifigkeitsmatrix. Damit wird aus der Skalargleichung aus K02 ein Eigenwertproblem.
Eigenproblem
x(t) = φ cos(ωt + α)
(K − ω²M)φ = 0
det(K − ω²M) = 0
Nicht jede relative Bewegung behält ihre Form. Nur Eigenvektoren φ reproduzieren sich während der Schwingung mit einer einzigen Eigenkreisfrequenz ω.
02 · Normalmoden
Die richtigen Koordinaten entkoppeln das gekoppelte System
Wegen der Symmetrie ergeben sich zwei orthogonale Bewegungsformen. Ihre modalen Koordinaten qₛ und qₐ entwickeln sich unabhängig wie zwei harmonische Oszillatoren.
Gleichphasige Bewegung
φₛ = (1/√2) [1, 1]ᵀ
ωₛ = √(k/m)
Beide Massen bewegen sich gleich. Die Kopplungsfeder ändert ihre Länge nicht und speichert deshalb in dieser Mode keine Energie.
Gegenphasige Bewegung
φₐ = (1/√2) [1, −1]ᵀ
ωₐ = √((k + 2κ)/m)
Die Massen laufen gegeneinander. Die Kopplungsfeder wird zusätzlich gedehnt und erhöht dadurch die effektive Steifigkeit dieser Mode.
qₛ = (x₁ + x₂)/√2
qₐ = (x₁ − x₂)/√2
Die Transformation ist eine orthogonale Drehung des zweidimensionalen Konfigurationsraums; für die Geschwindigkeiten gilt dieselbe Transformation. Information geht nicht verloren.
q̈ₛ + ωₛ²qₛ = 0
q̈ₐ + ωₐ²qₐ = 0
In den modalen Koordinaten verschwinden die Kreuzterme vollständig. Jede Normalmode besitzt ihre eigene, konstante modale Energie.
Δω = ωₐ − ωₛ
κ ↑ ⇒ Δω ↑
Ohne Kopplung sind beide Frequenzen entartet. Wachsende Kopplung hebt diese Entartung auf und trennt die Eigenfrequenzen.
03 · Interactive Lab
Two-Oscillator Mode Lab
Verändere Masse, äußere Steifigkeit, Kopplung und Anfangsauslenkungen. Die Mechanik zeigt x₁ und x₂; die Diagramme zeigen gleichzeitig physikalische und modale Koordinaten.
0.55 kg
580 N/m
060 N/m
−1+1 m
−1+1 m
Presets
Mechanisches Systemm — k — κ — k — m
läuft
x₁0.800m
x₂0.000m
qₛ0.566m
qₐ0.566m
v₁0.000m/s
v₂0.000m/s
Physikalische Koordinatenx₁(t), x₂(t)
Normiert auf eine gemeinsame Referenzamplitude A₀; das Zeitfenster passt sich der Modenaufspaltung an.
x₁x₂
Modale Koordinatenqₛ(t), qₐ(t)
Die beiden modalen Amplituden bleiben unabhängig; nur ihre Phasen laufen mit verschiedenen Eigenfrequenzen.
qₛqₐ
Interpretation
Modenmischung und Schwebung
Beide Normalmoden sind angeregt. Die physikalischen Koordinaten zeigen ihre Überlagerung.
04 · Energie und Schwebung
Die Moden tauschen keine Energie aus — die physikalischen Massen können trotzdem „Energie wandern“ lassen
Die scheinbare Wanderung entsteht aus der Überlagerung zweier unabhängiger Eigenbewegungen. Das ist ein zentraler Unterschied zwischen physikalischen und modalen Koordinaten.
E = ½m(v₁² + v₂²)
+ ½k(x₁² + x₂²) + ½κ(x₂ − x₁)²
Ohne Dämpfung und Anregung bleibt E konstant; die Kopplungsfeder ist lediglich ein zusätzlicher Energiespeicher.
Eₛ = ½m q̇ₛ² + ½mωₛ²qₛ²
Eₐ = ½m q̇ₐ² + ½mωₐ²qₐ²
Die Eigenbasis diagonalisiert nicht nur die Gleichungen, sondern auch die Energie des linearen symmetrischen Systems.
x₁ = A cos(ω̄t) cos(Δωt/2)
x₂ = A sin(ω̄t) sin(Δωt/2)
Mit ω̄ = (ωₛ + ωₐ)/2 wird die langsame Hülle durch Δω bestimmt. Je schwächer die Kopplung, desto langsamer die Schwebung.
05 · Modellgrenzen und Übergang
Zwei Massen sind bereits ein räumliches System — aber noch keine Welle
Gleiche Massen, gleiche Außenfedern, lineare Federn, keine Dämpfung und keine äußere Kraft halten das Eigenproblem analytisch transparent.
Die Normalmode gehört nicht einer einzelnen Masse. Sie ist eine Eigenschaft des gesamten gekoppelten Netzwerks aus Trägheiten und Steifigkeiten.
Mit jeder weiteren gekoppelten Masse entstehen weitere Moden. K10 macht daraus eine Kette und führt den Langwellen-Grenzübergang zur Wellengleichung aus; M02K01 setzt dort mit der räumlichen Wellenausbreitung an.
M01 · K10
Vom Oszillator zum
Vom Oszillator zum
verteilten System
Aus zwei gekoppelten Freiheitsgraden wird eine Kette aus N lokalen Zuständen. Die Dynamik bleibt linear und deterministisch, erhält aber eine räumliche Struktur: Jeder Ort koppelt an seine Nachbarn, und aus einem einzelnen Eigenwert entsteht ein ganzes Modenspektrum.
01 · Diskrete räumliche Dynamik
Jede Masse ist lokal — die Dynamik gehört dem gesamten Netzwerk
Wir betrachten N identische Massen m im Abstand a, gekoppelt durch identische lineare Federn mit Steifigkeit κ. Die Endpunkte sind fest: x₀ = xN+1 = 0.
Lokale Bewegungsgleichung
m ẍₙ = κ(xₙ₊₁ − 2xₙ + xₙ₋₁)
n = 1, …, N
Die diskrete zweite Differenz misst die lokale Krümmung der Auslenkung. Eine Masse beschleunigt nur dann, wenn ihre beiden Nachbarfedern keine ausgeglichene Kraft liefern.
Systemdimension
x = [x₁, x₂, …, xN]ᵀ
z = [x, ẋ]ᵀ ∈ ℝ2N
K01 begann mit zwei Zustandsgrößen. Die Kette besitzt 2N Zustände; im Kontinuumsgrenzfall wird daraus kein längerer Vektor, sondern ein Feld ψ(x,t).
Tridiagonale Kopplung
M ẍ + Kx = 0
K ∝ tridiag(−1, 2, −1)
Die Topologie ist lokal: K koppelt nur direkte Nachbarn. Trotzdem sind die Eigenvektoren global und erstrecken sich über die gesamte Kette.
02 · Modenspektrum und Gitterdispersion
N Freiheitsgrade erzeugen N Normalmoden
Für feste Enden sind die Eigenformen diskrete Sinusprofile. Der Modenindex p bestimmt zugleich die räumliche Phase zwischen benachbarten Gitterpunkten.
φₙ(p) = sin(qp n a)
qp = pπ / [(N+1)a]
p = 1,…,N. q ist hier die räumliche Wellenzahl der diskreten Mode; θ = qa ist der dimensionslose Phasenschritt pro Gitterabstand.
ω(q) = 2√(κ/m) · sin(qa/2)
0 < qa < π
Die Frequenz wächst nicht unbegrenzt proportional zu q. Nahe der Gittergrenze qa → π krümmt die Dispersionsrelation ab; die diskrete Mikrostruktur wird dynamisch sichtbar.
sin(qa/2) ≈ qa/2
ω ≈ c q
c = a√(κ/m)
Für qa ≪ 1 sieht die Kette lokal glatt aus. Dann wird die Dispersionskurve nahezu linear und die mikroskopische Gitterstruktur verschwindet aus der führenden Dynamik.
Wichtig: Eine einzelne reelle Normalmode mit festen Enden ist eine stehende Bewegung. Erst M02 entwickelt aus der räumlich gekoppelten Dynamik systematisch laufende Wellen, Phasenfortschritt und Transport.
03 · Interactive Lab
Discrete Chain & Continuum Lab
Wähle Kettenlänge und Modenindex. Die Animation zeigt die exakte Normalmode; die Diagramme koppeln räumliche Eigenform und diskrete Dispersionsrelation an den Langwellen-Grenzfall.
630
116
0.84 kg
660 N/m
0.251.00 m
Diskrete Kette · normierte GeometrieExakte longitudinale Normalmode
läuft
Die Massen bewegen sich entlang der Federachse. Die Auslenkung ist zur Lesbarkeit skaliert; der gezeichnete Kettenabstand ist normiert, während a im Modell die physikalische Gitterkonstante bleibt.
Räumliche Eigenformφn(p) = sin(qpna)
Linie: kontinuierliche Sinus-Hilfe. Punkte: tatsächlich vorhandene Freiheitsgrade der Kette.
DispersionsrelationDiskretes Gitter vs. Kontinuum
Normiert mit Ω = ω / [2√(κ/m)]. Die lineare Näherung gilt nur für θ = qa ≪ 1.
sin(θ/2)θ/2
Interpretation
Langwellenregime
Viele Gitterpunkte liegen innerhalb einer Wellenlänge; die kontinuierliche Näherung beschreibt die Eigenfrequenz nahezu exakt.
04 · Kontinuumsgrenzfall
Die diskrete zweite Differenz wird zur räumlichen Ableitung
Setze xₙ(t) = ψ(na,t) und entwickle die Nachbarwerte um x = na. Die Mikrodynamik liefert direkt die makroskopische Feldgleichung.
ψ(x±a) = ψ ± aψx + ½a²ψxx ± …
ψ(x+a) − 2ψ + ψ(x−a)
= a²ψxx + (a⁴/12)ψxxxx + …
Die ungeraden Ableitungen heben sich wegen der symmetrischen Nachbarschaft exakt auf.
ψtt = c² ψxx
c = a√(κ/m)
Die 1D-Wellengleichung ist kein neues Postulat: Sie ist die führende Langwellenbeschreibung derselben lokalen Kopplungsdynamik.
ψtt = c²ψxx
+ (c²a²/12)ψxxxx + …
Der nächste Term erinnert daran, dass das Kontinuum nur eine Näherung ist. Er erzeugt genau die Krümmung der diskreten Dispersionsrelation bei kurzen Wellenlängen.
05 · Systemische Synthese von M01
Vom lokalen Zustand zur räumlich gekoppelten Dynamik
Ein einzelner Freiheitsgrad etabliert Rückkopplung, Eigenfrequenz und konservativen Energieaustausch.
Dämpfung und äußere Kräfte machen aus der Eigenbewegung ein offenes dynamisches System mit Frequenzselektion.
Mehrere Freiheitsgrade erzeugen globale Eigenformen und schließlich eine Feldbeschreibung im Langwellenlimit.
Übergabe an M02 · Wave Fundamentals
ψtt = c²ψxx ist jetzt hergeleitet — M02 untersucht, was diese Gleichung räumlich transportiert.
Als Nächstes folgen laufende Wellen, Wellenlänge, Phase und Phasengeschwindigkeit, Wellenzahl, Energiefluss, Gruppenlaufzeit, Dispersion und Dämpfung. M01 liefert dafür die dynamische Mikrobasis.