算法详情请参考书籍《系统辨识及其MATLAB仿真》侯媛彬.，这里不做详细介绍。（如果有需要的话，请在评论区留言，如果留言数足够多，可以更新一期）
好的，直接来看编程实例：
题目：利用最小二乘一次完成算法对如下系统进行参数辨识：
假设v（k）服从正态分布的白噪声，输入信号u（k）采用幅值为1的M序列，利用LS算法进行参数估计。
分析：
式中，na=2，nb=2 利用式（3.7）（3.8）可写出 h(k)=[-z(k-1)，-z(k-2)，u(k-1)，u(k-2)]T θ=[a1，a2，b1，b2]T
数据长度取L=500； 加权矩阵为单位阵I（即无加权）
观测ZL和HL的构造，可参考公式（3.16）（3.17）：
