程轶平副教授的状态空间辨识算法MATLAB代码

这是我在《北京交通大学学报》2014年的论文《基于多元线性回归的线性系统状态空间辨识算法》配套MATLAB代码，用来生成论文中的图的。 读者有不明白的，可以联系我。

function [Abar B K] = abkest（X u y）

if （nargin ~= 3）
    error（‘ABKEST: Incorrect number of arguments‘）;
end

nx=size（X1）;
nu=size（u1）;
ny=size（y1）;

LX=size（X2）-1;
if （LX ~= size（u2））
    error（‘ABKEST: the lengths of x and u must match!‘）;
end
if （LX ~= size（y2））
    error（‘ABKEST: the lengths of x and y must match!‘）;
end

Xp=X（:1:LX）;
wp = [Xp; u; y];
Xf=X（:2:LX+1）;

[USV] = svd（wp）;
theta = Xf*V*S‘*inv（S*S‘）*U‘;

Abar=theta（:1:nx）;
B=theta（:nx+1:nx+nu）;
K=theta（:nx+nu+1:nx+nu+ny）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        535  2011-12-15 01:01  SSARX\abkest.m

     文件       1049  2013-08-12 16:43  SSARX\backup.m

     文件        507  2013-06-18 16:49  SSARX\block_hank.m

     文件        288  2011-12-15 00:27  SSARX\cest.m

     文件       1086  2013-08-12 16:50  SSARX\chengsubs.m

     文件       2220  2012-01-12 15:04  SSARX\chengsubs_10runs.m

     文件      18500  2013-08-27 16:56  SSARX\figure1.emf

     文件      17492  2013-08-27 16:56  SSARX\figure2.emf

     文件      19472  2013-08-27 16:56  SSARX\figure3.emf

     文件      19860  2013-08-27 16:56  SSARX\figure4.emf

     文件      61424  2013-10-18 15:45  SSARX\newfig.emf

     文件      19396  2013-10-18 14:12  SSARX\newfig1.emf

     文件     670614  2013-10-18 16:12  SSARX\newfig2.bmp

     文件      18388  2013-10-18 14:12  SSARX\newfig2.emf

     文件      22955  2013-10-18 16:06  SSARX\newfig2.eps

     文件      11779  2013-10-18 16:07  SSARX\newfig2.png

     文件      20368  2013-10-18 14:13  SSARX\newfig3.emf

     文件      20756  2013-10-18 14:16  SSARX\newfig4.emf

     文件        529  2013-08-23 00:55  SSARX\plotdemo.m

     文件        415  2011-12-14 19:55  SSARX\predest.m

     文件       1758  2013-06-18 16:49  SSARX\ssarx.m

     文件       6215  2013-10-18 14:09  SSARX\SSARX_MLR2.m

     文件       6613  2013-10-18 16:16  SSARX\SSARX_MLR3.m

     文件        582  2013-06-18 16:49  SSARX\stackdata.m

     文件      47640  2013-10-18 16:24  SSARX\untitled.fig

     目录          0  2013-10-18 16:24  SSARX

----------- ---------  ---------- -----  ----

               990441                    26