增广矩阵法实现最小二乘的系统辨识.docVIP

北京工商大学 《系统建模与辨识》课程 上机实验报告 （2014年秋季学期） 专业名称 ： 控制工程 上机题目 ： 增广矩阵法实现最小二乘的系统辨识 专业班级 ： 2015 年 1 月 一、实验目的 通过仿真实验掌握增广矩阵法进行参数估计的原理和方法。 二、实验原理 考虑CARMA模型 其中： 还可表示为 ： 其中： 以上两向量均由2n维扩展为（2n+r）维。 在式（2）中若是已知量，则可直接用最小二乘法估计出，但是｛｝是不可测的未知量。一个简单可行的方法是用计算的代替，｛｝由下式递推得出： ? （计算残差） 其中： ? ? 三、实验内容 按附表9-11、表9-12给出的二阶线性离散系统的输入输出数据，用增广矩阵法进行参数估计。对任务进行方案设计，画出实验流程图，用MATLAB编程实现。撰写实验报告。 四、流程图 五、代码实现 clc clear all; close all; u=[1.147,0.201 ,-0.787 ,-1.584 -1.052, 0.866 ,1.152 ,1.573 ,0.626, 0.433 ... -0.958,0.810, -0.044, 0.947, -1.474, -0.719, -0.086, 1.099, 1.450, 1.151... 0.485,1.633, 0.043 ,1.326, 1.706, -0.340, 0.890,0.433, -1.177, -0.390... -0.982,1.435 ,-0.119 ,-0.769, -0.899, 0.882,-1.008, -0.844, 0.628,-0.679 ... 1.541,1.375, -0.984, -0.582, 1.609,0.090, -0.813, -0.428 ,-0.848,-0.410... 0.048,-1.099 ,-1.108 ,0.259,-1.627, -0.528, 0.203, 1.204,1.691, -1.235... -1.228,-1.267, 0.309,0.043, 0.043, 1.461, 1.585, 0.552,-0.601, -0.319... 0.744 0.829,-1.626, -0.127, -1.578, -0.822, 1.469, -0.379,-0.212, 0.178... 0.493 -0.056, -0.1294, 1.228, -1.606, -0.382, -0.229, 0.313,-0.161, -0.810... -0.277 0.983, -0.288, 0.846, 1.325, 0.723, 0.713, 0.643 0.463,0.786... 1.161,0.850, -1.349, -0.596, 1.512, 0.795, -0.713, 0.453,-1.604, 0.889... -0.938, 0.056, 0.829,-0.981, -1.232, 1.327, -0.681,0.114,-1.135, 1.284... -1.201 0.758, 0.590, -1.007, 0.390, 0.836,-1.52, -1.053, -0.083, 0.619... 0.840 -1.258, -0.354, 0.629, -0.242,1.680, -1.236, -0.803,0.537, -1.100... 1.417,-1.024, 0.671, 0.688,-0.123, -0.952, 0.232, -0.793,-1.138, 1.154... 0.206,1.196, 1.013,1.518, -0.553, -0.987, 0.167, -1.445,0.630, 1.255... 0.311,-1.726,0.975, 1.718, 1.360, 1.667, 1.111, 1.018, 0.078,-1.66