|
|
u=[wfm_change A8_change c_etaf c_mf c_pif c_etac c_mc c_pic sigma_B];
# b9 }/ }9 g* E; s: yy=[T41C Pit];( ~% q/ R* C/ X8 O. b
x=iddata(y,u,0.025);
3 U, @8 u0 e4 w% D+ F( Q0 l5 {1 _& c* c
% ident' C( V! u& l; i" ^ B8 r
% 辨识工具箱! v6 m7 J1 c8 o) I. a5 B' T
% import data object->workapce value (输入名字在此代码中为x)->estimate->阶数2->得到ss1结构体
% i4 X3 D/ |% Y( g0 W. YA=ss1.A;
, k1 g% X% }1 w( d( ]B=ss1.B;* w" z& \4 l2 B$ }6 {/ t
C=ss1.C;
" D# O6 L7 V: `/ z/ H" nD=ss1.D;
, t2 L$ I/ }# Y( G5 |, l具体操作是这样的,所有的输入变量作为整个向量(上面代码中的u),输出变量(y)为整个向量,然后通过iddata命令重整成x结构体(0.025是我的采样周期),打开辨识工具箱后,导入x就好了。可以试一下 |
|