2016-04-14 206 views
0

与输入矢量u,输出矢量y和干扰E中的系统的状态空间模型以在连续时间以下形式:如何从ssest matlab函数中提取状态空间的残差(e(t))?

d(X(t))/ dt的=斧(T)+柯(吨)
Y(T)= CX(T)+ E(t)的

在MATLAB函数ssest我能够提取A,C和K个参数,但是我在真正感兴趣提取e(t)参数以建立残基的协方差。

data = iddata(X(:,1:100,10)',[],1); 
    opt = ssestOptions('EstCovar',true); 
    sys = ssest(data,[1:10],opt); 

    A = sys.A;      % state transition matrix 
    C = sys.C;      % observation matrix 
    K = sys.K;      % inovation form 

有人可以告诉我,如果使用此功能,我能够这样做吗?

回答

0

e(t)是一种干扰,它不是模型的一部分(它是模型的输入),所以“提取”它是没有意义的,你不能这样做。这是你想要的。

查看documentation page on ssest了解更多详情。

+0

谢谢你的回复!但@ am304,关于*状态空间估计*,它估计了使残基最小化的参数,对吗?所以如果这是真的,我应该期待最终获得残留信息,作为最小化迭代的结果。 – tTimoteof

+0

查看['idssdata'](http://uk.mathworks.com/help/ident/ref/idssdata.html)和[各种模型识别方法比较](http://uk.mathworks.com /help/ident/examples/comparison-of-various-model-identification-methods.html) – am304

+0

关于第一个建议,它返回参数的uncertenty。但是,我怎样才能将这些信息转换为模型的不确定性?这是因为考虑到在这个问题的介绍中提到的等式,矩阵维数必须匹配,并且这些等式也清楚地表明了模型的不确定性。 我在这里看[在Matlab状态空间模型](http://www.mathworks.com/help/ident/state-space-models.html)但是没有运气。 – tTimoteof