估计状态空间模型与数据:问题的自由/固定参数模型

2视图(30天)

显示旧的评论

foboo 2016年9月7日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/302206-state-space-model-estimation-with-data-issue-with-free-fixed-parameters-of-model

评论道: 马塞尔•亨德里克斯2020年3月27日

你好,我想估计给定模型的几个参数(噪音)。我想(想暂时修复它):

                         = (1 0 1 0;0 1 0 1;0 0 1 0;0 0 0 1);
                        
                         B = [];
                        
                         C = [(2), 0 (2, 2)];
                        
                         D = [];
                        
                         K = [0 (2, 2);的(2,2)];

为了使用党卫军后来,我也设置以下值:

                         x0 = 0 (4,1);
                        
                         Ts = 1;
                        
                         K (K = = 1) =南;%值标记为未知w.r.t.初始值
                        
                         init_sys = idss (A, B, C, D, K, x0, Ts);

我想估计参数 K ,我固定其他如下:

                         init_sys.Structure.K。自由= K ~ = 0;
                        
                         init_sys.Structure.C。自由= false(大小(C));
                        
                         init_sys.Structure.A。自由= false(大小(A));

此外,我设置:

                         init_sys。InitialState =“估计”;
                        
                         选择= ssestOptions;
                        
                         opt.SearchOption。麦克斯特= 100;

然后运行

系统= ss(数据、init_sys选择)

它给了我

                         K > >系统
                        
                         系统=
                        
                         离散时间确定:状态空间模型
                        
                         x (t + Ts) = x (t) + K e (t)
                        
                         y (t) = C x (t) + e (t)
                        
                         一个=
                        
                         x1x2x3x4
                        
                         x10.9014-0.03647-0.034790.005559
                        
                         x20.1841.0730.23870.1506
                        
                         x30.0393-0.13730.8788-0.3308
                        
                         x4-0.06804-0.25450.27840.6195
                        
                         C =
                        
                         x1x2x3x4
                        
                         日元301年18.766.241-2.569
                        
                         y2-34.7317.92.1361.577
                        
                         K =
                        
                         日元y2
                        
                         x10.005925-0.00643
                        
                         x20.0021050.06914
                        
                         x30.028190.08107
                        
                         x40.17480.4147
                        
                         样本时间:1秒
                        
                         参数化:
                        
                         结构化形式(一些固定系数A、B、C)。
                        
                         数量免费的系数:4
                        
                         使用“idssdata”,“getpvec”,“getcov”为参数和他们的不确定性。
                        
                         状态:
                        
                         估计使用党卫军时域数据“数据”。
                        
                         适合估算数据:[93.3692.87);%(预测聚焦)
                        
                         MSE消防工程:2.099:8.258

我为什么党卫军修改/估计新的参数矩阵一个 , C 吗?和我怎么能改变这/我做错了什么?

编辑:当我正确理解它, 党卫军转换的输入系统的另一个代表和执行评估噪音和 K 使用它。有办法防止这种情况(或将再次到原始的形式)?

相反是写在文档党卫军修复一些规范的形式表示失败:

                         > >系统= ss(数据、init_sys“形式”,“模态”选择)
                        
                         错误使用ss(第132行)
                        
                         没有房地产类的“形式”“采用”匹配的标识符。使用属性来获取这个类的属性列表。