加权最小二乘法系统辨识

案例：对下列数学模型进行系统辨识

式中，a1=1.5，a2=0.7，b1=1.0，b2=0.5，V(K)服从N(0,1)的随机噪声；输入u(k)采用4阶M序列，其幅值为5，加权矩阵Wm=I

解：由于输入信号为4阶M序列，所以M序列的循环长度为L=2^4-1=15，因此，设输入信号的取值为从k=1到k=16的M序列，于是可得：

根据上式可一次计算出的结果。

辨识流程如下：

开始➡️产生输入信号M序列➡️产生输出信号z(k)➡️给出样本矩阵Hm和Zm➡️估计参数➡️分离估计参数a1，a2，b1，b2➡️画图：输入、输出信号和估计参数➡️结束

程序代码看文末

仿真结果如下：

输入信号为幅值为5的4阶M序列：

在输入信号作用下的输出信号

辨识结果如下：

参数	a1	a2	b1	b2
真值	1.5	0.7	1.0	0.5
估计值	1.496	0.697	0.966	0.482