Aufgabe 12

• Schreiben Sie eine Funktion
  • [M, C] = createMatrices(chain),
• die zu einer Schwingerkette die Massenmatrix M und die Steifigkeitsmatrix C berechnet.
• Zur einfachen Lösung der Differentialgleichung der Schwingerkette erzeugen Sie eine Funktion
  • [t, x] = solveVibrationODE(M, C, x0, tEnd, tStep),
• mit folgenden Parametern

  • M	Massenmatrix
    C	Steifigkeitsmatrix
    x0	Anfangsauslenkungen
    tEnd	Endzeit der Lösung
    tStep	Ausgabe-Schrittweite

  • (Anfangszeit sei t₀ = 0, Anfangsgeschwindigkeiten v_i(0) = 0)
• und Ergebniswerten

  • t	Vektor mit Zeitwerten
    x	Array mit Werten der Auslenkung zur Zeit t für jede Masse

• Verwenden Sie zum Testen die Demo-Schwingerkette und plotten Sie die Ergebnisse.