関数 > 求解と最適化 > 微分方程式の解法 > 例: 状態空間ソルバー
  
例: 状態空間ソルバー
ソルバー statespaceを使用して 1 次の微分方程式 (ODE) の解を状態空間に表示したものを求めます。
特に、調和振動方程式の右辺が 0 である場合の外部から力を加えない調和振動子の解を求めます。
クリックしてこの式をコピー
* 
過減衰、臨界減衰、不足減衰の 3 つの場合の解があります。
過減衰の解
1. 過減衰の解を求める数式を次のように記述します。
クリックしてこの式をコピー
2. 初期条件 (オブジェクトの質量、弾性定数、積分区間の始点と終点、点の数) を定義します。
クリックしてこの式をコピー
クリックしてこの式をコピー
クリックしてこの式をコピー
クリックしてこの式をコピー
クリックしてこの式をコピー
クリックしてこの式をコピー
3. 系の自然 (共振) 周波数を設定します。
クリックしてこの式をコピー
4. 過減衰条件が存在することを検証します。
クリックしてこの式をコピー
5. ODE を行列形式で記述します。
クリックしてこの式をコピー
6. statespace 関数を次のように呼び出します。
クリックしてこの式をコピー
7. 解をプロットします。
クリックしてこの式をコピー
臨界減衰の解
1. 系の自然 (共振) 周波数を設定します。
クリックしてこの式をコピー
2. 臨界減衰条件が存在することを検証します。
クリックしてこの式をコピー
3. ODE を行列形式で記述します。
クリックしてこの式をコピー
4. statespace 関数を次のように呼び出します。
クリックしてこの式をコピー
5. 解をプロットします。
クリックしてこの式をコピー
不足減衰の解
1. 系の自然 (共振) 周波数を設定します。
クリックしてこの式をコピー
2. 不足減衰条件が存在することを検証します。
クリックしてこの式をコピー
3. ODE を行列形式で記述します。
クリックしてこの式をコピー
4. statespace 関数を次のように呼び出します。
クリックしてこの式をコピー
5. 解をプロットします。
クリックしてこの式をコピー