[发明专利]系统辨识装置和方法以及记录有系统辨识程序的存储介质

权利要求书

1.一种系统辨识装置，所述系统辨识装置根据输入输出数据进行系统的高速实时辨识，所述系统辨识装置的特征在于，

包括滤波器，所述滤波器通过作为评价基准将从干扰向滤波器误差的最大能量增益确定为抑制成比预先提供的上限值γ_f²小而相对于干扰稳健，

所述滤波器对于由以下公式(11)～(13)表示的状态空间模型，由以下公式(49)、(50)、(55)、(58)、(60)～(62)提供，

推定所述系统的状态x_k的状态推定值x＾_k|k，

其中，C_k和W由以下公式定义，

是C_k＝[c_k，...，c_k-N+1]的第一列向量，

H_k＝[u_k、…、u_k-N+1]^T

其中，

x_k：状态向量或仅为状态，未知，其成为推定的对象；

x₀：初始状态，未知；

w_k：系统噪声，未知；

V_k：观测噪声，未知；

y_k：观测信号，成为滤波器的输入，已知；

z_k：输出信号，未知；

G_k：驱动矩阵，执行时已知；

H_k：观测矩阵，已知；

x＾_k|k：使用观测信号y₀～y_k的时刻k的状态x_k(未知系统的脉冲响应)的推定值，由滤波器方程式提供；

x＾_k+1|k：使用观测信号y₀～y_k的时刻k的状态x_k+1的一步预测值，由滤波器方程式提供；

x＾_-1|-1：状态的初始推定值，原本为未知，但是为了便于说明而使用0；

K_s、k：滤波器增益，由矩阵Σ＾_k|k-1得到；

R_e、k：辅助变量，2×2矩阵；

N：状态向量的维度(抽头数)，事先提供；

μ～_k：K^U_k的第N+1个1×2行向量，由K^U_k得到；

W：权重矩阵，由γ_f确定；

e_k：事后前方预测误差，二维列向量；

e_k：事后后方预测误差，二维列向量；

γ_f：衰减水平，设计时提供；

μ：步长(松弛因子)，μ>0的实数，事先提供。

2.根据权利要求1所述的系统辨识装置，其特征在于，

使用公式(60)，对于公式(58)、公式(55)、公式(49)和公式(50)，使用将公式(58)、公式(55)、公式(49)和公式(50)变形后的以下公式，

其中，μ是步长，是μ>0的实数。