生きていくことって結構難しいなと感じることがあるかもしれません。
たくさんやらなければいけないことがあって、プレッシャーを感じながら、
ついつい焦ってしまい、失敗して、結局やり直し。
「忙しい時こそ、焦らない。早くやってもいいけど、焦ってやらない」
ほっとお茶でも飲んで、落ち着いて、やっていきましょう^^
こんにちは、データサイエンティストひろき🐶です。
伝達関数の状態空間表現
さて、
前回こちらのブログで紹介した下図の1自由度振動モデル、

この伝達関数
G=xy=cs+kms2+cs+k
を状態空間表現で記述してみたいと思います。
まず、状態空間表現の基本形は下記2式(状態方程式と観測方程式)でした。
˙x=Ax+Bu ・・・状態方程式
y=Cx+Du ・・・観測方程式
ちなみに、
x=(x˙x)
u:入力 とします。
つまり今回は、
- 伝達関数Gから状態空間表現A,B,C,Dを求める方法
を紹介します。
まず状態方程式
˙x=Ax+Bu
をラプラス変換します。
下記のようにラプラスの演算子sを用いて書くことができます。
sx=Ax+Bu
sx−Ax=Bu
(sI−A)x=Bu
x=(sI−A)−1Bu
これを観測方程式へ代入すると、
y={C(sI−A)−1B+D}u
uが入力、yが出力なので、伝達関数Gは
G=C(sI−A)−1B+D
となります。
AとBはこちらの記事で求めた下記を使えることが、
こちらの記事でChatGPTが答えてくれています。
A=(01−k/m−c/m)
B=(01/m)
また、伝達関数がプロパー(分子の次数が分母の次数より小さいこと)の時はD=0となることがわかっている。
残るはC。元々の数理モデルから求めた伝達関数G
G=xy=cs+kms2+cs+k
と比較して、Cを求めていきます。
つまり、
C=(c1c2)
と置いて、上記A~Dを下記Gへ代入して
G=C(sI−A)−1B+D
伝達関数
G=xy=cs+kms2+cs+k
と係数を比較します。実際に紙と鉛筆で手計算してみてください。
すると、 C=(kc)
となります。
そしてMatlabには、
伝達関数から状態空間表現の係数(A,B,C,D)を求める便利な関数 “tf2ss”
たるものがあります。
まじ便利です。
詳しい使い方はこちらを参照ください。
もし詳しい式変形が知りたければ、
ChatGPTに「matlabのtf2ssのアルゴリズムは?」と聞いてみてください。
きっと答えがかえってくるはず。
それでは今日はこの辺で
最後まで読んでくださりありがとうございました。
Love&Respect♡
Hiroki🐶