matlab - 質量ばね系の線形方程式を解く

Answer 1

cとkだけ？mはどうですか？それも要らない？

これは同次 ODE です。配列は時間 0 での初期変位を示しますが、速度はどうでしょうか? 2 つの初期条件が必要です。開始時も速度がゼロであると想定するつもりですか? それとも、あなたが持っている 2 つの点に基づいて最初の差を計算することを選択しますか?

値をプロットし、それらがどのように見えるかを確認し、観察内容に基づいて定数について知識に基づいた推測を試みることをお勧めします。この ODE の閉じた形式の解を知っています。与えられたデータに基づいて、値がどうあるべきかを判断できる場合があります。

それを行った後、いくつかの点で計算された差に基づいて、未知数の行列解を作成することをお勧めします。これにより、問題が定数の最小二乗解になります。定数の 1 つのセットがすべての時点で方程式を満たさない可能性があります。

もう 1 つのアプローチは、データの FFT を取得し、それを周波数ドメインのソリューションと比較することです。