状態空間モデルで行列内の値を固定するには？

1 次查看（过去 30 天）
Haruto Yamazaki
Haruto Yamazaki 2019-12-27
回答： Hiroumi Mita 2019-12-27
状態空間モデルのパラメータ推定を行っています．
X'=AX+Bu
Y=CX+Du
において，
A=[0 1 0 0;NaN 0 NaN 0;0 0 0 1;NaN 0 NaN 0];
B=[0;NaN;0;0];
C=[1 0 0 0;0 0 1 0];
D=[0;0];
設定したモデルは上記の通りです．
NaN以外の行列の値を固定して，ssestなどを用いて値を推定したいのですが，その固定方法等がよくわかりません．
おそらくidssのstructureを使う様なのですが，例えば
m.Structure.A.Free(2,1) = false;
m.Structure.A.Value(2,2)
m.Structure.A.Free(2,3) = false;
m.Structure.A.Value(2,4)
のようにしても固定が出来ません.
どなたかご教授をお願いします．

请先登录，再回答此问题。

采纳的回答

Hiroumi Mita
Hiroumi Mita 2019-12-27
状態空間内部の任意の箇所だけ推定するには
idgreyで設定し、pem等のシステム同定の関数を使います。
idgreyとそのサンプルをhelpで調べてみてください。
  1 个评论
Haruto Yamazaki
Haruto Yamazaki 2019-12-27
编辑：Haruto Yamazaki 2019-12-27
ご回答ありがとうございます。
https://jp.mathworks.com/help/ident/ref/idgrey.html#btd7ziz-1-odefun
idgreyでは上記のサイトのConfigure Estimable Parameter of Grey-Box Model
を参考にしたのですが，
＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
エラー: idgrey (line 372)
指定されたパラメーター、サンプル時間およびオプション引数のセットを使用して ODE 関数 "LinearPendulum" を適切に評価で
きませんでした。評価中に次のエラー メッセージが生成されました:
入力引数が多すぎます。
＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
といったエラーが出てきてしまいます．
回答を頂きながらおこがましいのですが，
このエラーの原因を教えていただけないでしょうか。
以下，プログラムです．
＿＿＿＿＿＿＿＿＿＿＿＿＿ODE_file＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
function [A,B,C,D] = LinearPendulum(a,b,c,d,e)
A=[0 1 0 0;a 0 b 0;0 0 0 1;c 0 d 0];
B=[0;e;0;0];
C=[1 0 0 0;0 0 1 0];
D=[0;0];
end
＿＿＿＿＿＿＿＿＿＿＿＿＿command window＿＿＿＿＿＿＿＿＿＿＿＿
odefun = 'LinearPendulum' ;
a=1;
b=1;
c=1;
d=1;
e=1;
parameter={'data1',a;'data2',b;'data3',c;'data4',d;'data5',e};
fcn_type = 'c' ;
sys=idgrey(odefun,parameter,fcn_type);

请先登录，再进行评论。

更多回答（1 个）

Hiroumi Mita
Hiroumi Mita 2019-12-27
例題を見ると、
サンプル時間Tsが関数LinearPendulum
の引数に必要なようです。
function [A,B,C,D] = LinearPendulum(a,b,c,d,e,Ts) %<-ココ
A=[0 1 0 0;a 0 b 0;0 0 0 1;c 0 d 0 ];
B=[0;e;0;0 ];
C=[1 0 0 0;0 0 1 0 ];
D=[0;0 ];
end

类别

在 Help Center 和 File Exchange 中查找有关 状態空間モデル 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!