「シミュレーション」の版間の差分

359行目: 359行目:
<math>\frac{d}{dt} \left(\begin{array}{c}x_1(t)  \\x_2(t)\end{array}\right) =\left(\begin{array}{cc}0 & 1  \\-\frac{K}{M}  &  -\frac{D}{M}\end{array}\right)\left(\begin{array}{c}x_1(t)  \\x_2(t)\end{array}\right)+ \left(\begin{array}{c}0  \\\frac{1}{M}\end{array}\right)u(t)</math>
<math>\frac{d}{dt} \left(\begin{array}{c}x_1(t)  \\x_2(t)\end{array}\right) =\left(\begin{array}{cc}0 & 1  \\-\frac{K}{M}  &  -\frac{D}{M}\end{array}\right)\left(\begin{array}{c}x_1(t)  \\x_2(t)\end{array}\right)+ \left(\begin{array}{c}0  \\\frac{1}{M}\end{array}\right)u(t)</math>


状態変数を組み合わせたものを状態ベクトル&math(\boldsymbol{x}(t));とすれば,
状態変数を組み合わせたものを状態ベクトル<math>\boldsymbol{x}(t)</math>とすれば,


<math>\frac{d\boldsymbol{x}(t)}{dt}  =\left(\begin{array}{cc}0 & 1  \\-\frac{K}{M}  &  -\frac{D}{M}\end{array}\right)\boldsymbol{x}(t)+ \left(\begin{array}{c}0  \\\frac{1}{M}\end{array}\right)u(t)</math>
<math>\frac{d\boldsymbol{x}(t)}{dt}  =\left(\begin{array}{cc}0 & 1  \\-\frac{K}{M}  &  -\frac{D}{M}\end{array}\right)\boldsymbol{x}(t)+ \left(\begin{array}{c}0  \\\frac{1}{M}\end{array}\right)u(t)</math>


となり,これが状態方程式になります.出力は,状態量&math(x_1(t));そのものですから,
となり,これが状態方程式になります.出力は,状態量<math>x_1(t)</math>そのものですから,


<math>y(t) = \left(\begin{array}{cc}1 & 0\end{array}\right)\boldsymbol{x}(t)</math>
<math>y(t) = \left(\begin{array}{cc}1 & 0\end{array}\right)\boldsymbol{x}(t)</math>
372行目: 372行目:


<math>\boldsymbol{A} = \left(\begin{array}{cc}0 & 1  \\-\frac{K}{M}  &  -\frac{D}{M}\end{array}\right), \,\boldsymbol{B} = \left(\begin{array}{c}0  \\\frac{1}{M}\end{array}\right), \, \boldsymbol{C} = \left(\begin{array}{cc}1 & 0\end{array}\right), \, \boldsymbol{D} = 0</math>
<math>\boldsymbol{A} = \left(\begin{array}{cc}0 & 1  \\-\frac{K}{M}  &  -\frac{D}{M}\end{array}\right), \,\boldsymbol{B} = \left(\begin{array}{c}0  \\\frac{1}{M}\end{array}\right), \, \boldsymbol{C} = \left(\begin{array}{cc}1 & 0\end{array}\right), \, \boldsymbol{D} = 0</math>


=== オイラー法によるシミュレーション ===
=== オイラー法によるシミュレーション ===
95

回編集