Wird ein gedämpfter harmonischer Oszillator ohne Anfangsamplitude schwingen, wenn es Hintergrund-"Rauschen" gibt?

Angenommen, ich habe einen gedämpften harmonischen Oszillator, der in Ruhe ist, bequem ohne Anfangsamplitude sitzt und der Gleichung gehorcht

x ¨ + 1 Q x ˙ + x = 0

wobei x die vertikale Amplitude und Q der Qualitätsfaktor ist. Bei t = 0 , x = 0 .

Nehmen wir nun an, ich modelliere mein System so, dass es eine Art kleiner Störung, wie z. B. Wärme, enthält. Wir können dies als zufällige, Gaußsche Schwingungen modellieren: zum Beispiel weißes Rauschen. Die Gleichung wird:

x ¨ + 1 Q x ˙ + x = N ( t )

wo N ( t ) ist eine Zufallszahlenfunktion, die für eine Gaußsche Verteilung geformt ist.

Wird dieses Rauschen den Oszillator stören und ihm eine kleine Amplitude verleihen, und können wir für einen solchen Fall ein Diagramm wie das unten stehende erwarten?

Geben Sie hier die Bildbeschreibung ein

Dies ist eine Simulation, die ich durchgeführt habe, und ich frage mich, ob diese kleinen zufälligen Störungen den Oszillator auslösen und das typische, aber willkürliche, sinusförmige Verhalten verursachen.

Antworten (3)

Die Position der Masse als Funktion der Zeit ist einfach eine gefilterte Version des "Eingangs"-Signals des zufälligen Rauschens. Um dies im Frequenzbereich zu sehen, nehmen Sie die (Größe der) Fourier-Transformation beider Seiten und ordnen Sie sie neu an:

| X ( ω ) | = 1 ( 1 ω 2 ) 2 + 1 Q 2 ω 2 | N ( ω ) |

Zum ω = 1 , wir haben

| X ( 1 ) | = Q | N ( 1 ) |

Also für groß Q (stark unterdämpft), weist die Positionsfunktion eine starke sinusförmige Komponente auf ω = 1 .

Für kleinere Q , sollte die Positionsfunktion eine "geglättete" Version der Eingangsrauschfunktion sein, da die Frequenzen deutlich darüber liegen ω = 1 fallen mit einer Rate von etwa 40 dB/Dekade ab.

Ein Diagramm der (Betrags-) Übertragungsfunktion X ( ω ) N ( ω ) für verschiedene Werte von Q = 1 2 ζ sieht aus wie:

Geben Sie hier die Bildbeschreibung ein

Vielleicht ergänzen X ( ω ) bedeutet die Fourier-Transformation von x ( t ) , so dass ein imaginärer Wert für X ( ω ) bedeutet eine Phasenverschiebung von 90° gegenüber dem Eingang ( N ( t ) ). Und da PPG von Rauschen als Eingabe spricht, denke ich, dass ein Bode-Plot oder eine Übertragungsfunktion ebenfalls nützlich sein könnten.
Eine kleine Randnotiz, ich halte es für schlechte Angewohnheit, die Einheiten einer physikalischen Größe nicht zu erwähnen, also ω = 1 sollte sein ω = 1   r a d / s .
@fibonatic, in diesem Fall, ω normalisiert ist, was meiner Meinung nach angesichts des Kontexts offensichtlich ist. Andernfalls erfordert die gegebene Differentialgleichung Faktoren von ω 0 .

Ich glaube nicht, dass du wirklich eine Antwort brauchst. Die Antwort ist ja, und darüber hinaus haben Sie ein ziemlich solides Modell der Wirkung von Rauschen auf den gedämpften Oszillator erstellt. Ich gehe davon aus, dass Sie Frequenzen normalisiert haben, damit die Eigenfrequenz des Oszillators ω n ist eine Einheit.

Der einzige Faktor, den Sie nicht erwähnt haben und den Sie anscheinend übersehen haben, ist die spektrale Dichte des Rauschens N ( t ) . Sie können entweder davon ausgehen, dass das Rauschen im Vergleich zur Bandbreite Ihres Oszillators sehr breitbandig ist, in diesem Fall N ( t ) eine unkorrelierte Folge von Zahlen mit einer durch gegebenen Varianz sein σ 2 = η ich n f t j | H ( ω ) | 2 d ω wo η ist die spektrale Leistungsdichte des Rauschens und H ( ω ) die Fourier-Transformation der Einheitsimpulsantwort Ihres Oszillators. Ansonsten solltest du dir vielleicht Gedanken über die Art der mechanischen/sonstigen Effekte machen, die das Rauschen erzeugen und somit machen N ( t ) eine Folge von Variablen, die korreliert werden, indem Gaußsches weißes (unkorreliertes) Rauschen durch Modelle des Rauscherzeugungsmechanismus geleitet wird, die selbst wahrscheinlich Systeme sind, die durch lineare ODEs genau wie Ihres beschrieben werden. In beiden Fällen werden Sie also Gaußsches, unkorreliertes Rauschen entweder direkt in Ihr System oder durch ein modifiziertes System einspeisen.

Die Antwort muss ja sein, da eine Kraft ausgeübt wird, wird es eine Bewegung geben (wenn es keine Bewegung gäbe, würden die Newtonschen Gesetze natürlich verletzt).

Die Anwendung wird höchstwahrscheinlich winzig sein und der Dämpfungsterm wird dafür sorgen, dass sie winzig bleibt. Wenn es keinen Dämpfungsterm gäbe, könnten Sie starke Schwingungen bekommen.

Die Form des Diagramms wird höchstwahrscheinlich dem ähneln, was Sie basierend auf Schwingungen bei der Eigenfrequenz des Systems vorschlagen - aber ich fürchte, ich verstehe nicht genau, was Sie meinen N ( t ) und das wird natürlich einen starken Einfluss auf die Bewegung haben.

Das erinnert mich an das Problem mit einer Brücke in London , die keine Dämpfung hatte, wo die Eigenfrequenz nahe an der Frequenz von Menschen lag, die gingen. Sie haben die Brücke repariert, indem sie Dämpfung eingebaut haben.

Das war auch meine Intuition. Nun, in meinem Code habe ich implementiert N ( t ) einfach als (Zufalls-)Zahl, die mit der Zeit variiert. Ich bringe es zur LHS und löse die Gleichung: Am Ende habe ich zwei Gleichungen und benutze RK4 zum Simulieren. Außerdem, warum zum Teufel versuchen die Leute in diesem Video absichtlich, die Brücke zum Schwingen zu bringen!?
@PPG - die Leute fanden es schwierig, auf der Brücke zu gehen, während sie sich bewegte, und fanden, dass der einfachste Weg, sich zu bewegen, darin bestand, mit ihr im Takt zu bleiben ... das verschlimmerte das Problem natürlich. Es gibt eine ganze Menge im Internet darüber - eine interessante Geschichte darüber, was beim Design schief gelaufen ist ...