系统辨识实验报告

-.z.一、相关分析法(1)实验原理图1实验原理图本实验的原理图如图1。过程传递函数中；输入变量，输出变量，噪声服从，为过程的脉冲响应理论值，为过程脉冲响应估计值，为过程脉冲响应估计误差。过程输入采用M序列，其输出数据加白噪声得到输出数据。利用相关分析法估计出过程的脉冲响应值，并与过程脉冲响应理论值比拟，得到过程脉冲响应估计误差值。M序列阶次选择说明：首先粗略估计系统的过渡过程时间TS(通过简单阶跃响应)、截止频率fM(给系统施加不同周期的正弦信号或方波信号，观察输出)。本次为验证试验，系统模型，经计算，。根据式及式，则取值为1，此时，由于与N选择时要求完全覆盖，则选择六阶M移位存放器，即N=63。编程说明人机对话人机对话噪声标准差：sigma；生成数据周期数：r生成数据生成数据生成M序列；生成白噪声序列过程仿真过程仿真得到理论输出数据计算脉冲响应估计值计算脉冲响应估计值计算互相关函数，得到脉冲响应估计值计算脉冲响应估计误差计算脉冲响应估计误差计算脉冲响应理论值，得到脉冲响应估计误差图2程序流程图分步说明=1\*GB3①生成M序列：M序列的循环周期，时钟节拍，幅度，移位存放器中第5、6位的内容按"模二相加〞，反应到第一位作为输入。其中初始数据设为{1,0,1,0,0,0}。程序如下：=2\*GB3②生成白噪声序列：程序如下：=3\*GB3③过程仿真得到输出数据：如图2所示的过程传递函数串联，可以写成形如，其中。图2过程仿真方框图程序如下：=4\*GB3④计算脉冲响应估计值：互相关函数采用公式，互相关函数所用的数据是从第二个周期开场的，其中为周期数，取1-3之间。则脉冲响应估计值为：，。补偿量。程序如下：=5\*GB3⑤计算脉冲响应估计值：脉冲响应的理论值由式可计算得到。这时可得到过程脉冲相应估计误差。脉冲响应估计误差为：程序如下：数据记录=1\*GB3①当噪声标准差sigma=0.1，生成数据周期r为2时：脉冲响应估计误差为0.0121。脉冲响应估计曲线为图3所示。=2\*GB3②当噪声标准差sigma=0.5，生成数据周期r为1时：脉冲响应估计误差为0.0347。脉冲响应估计曲线为图4所示。=3\*GB3③当噪声标准差sigma=0.5，生成数据周期r为3时：脉冲响应估计误差为0.0258。脉冲响应估计曲线为图5所示。=4\*GB3④当噪声标准差sigma=1，生成数据周期r为3时：脉冲响应估计误差为0.0279。脉冲响应估计曲线为图6所示。图3sigma=0.1，r=2时脉冲响应估计曲线图4sigma=0.5，r=1时脉冲响应估计曲线图5sigma=0.5，r=3时脉冲响应估计曲线图6sigma=1，r=3时脉冲响应估计曲线结果分析实验中可以看到脉冲响应估计的曲线与理论曲线的重合度还是比拟高的，脉冲响应估计误差也比拟小，实验证明相关分析法的估计效果还是不错的。同时，经过实验可以得出结论：固定数据周期r，给定不同的噪声标准差sigma可以发现，噪声的方差越大，也就是信噪比越大，估计的效果越不好；固定噪声标准差sigma，选择不同的数据生成周期r可以发现，数据周期越大，估计的周期越多，估计的效果越好。二、最小二乘法根本最小二乘〔离线辨识〕残差为：最小二乘目标：残差平方和最小〔一阶导为0，二阶导>0〕。从上式看出，逆存在才有解，满足条件的u(k)：(1)伪随机；(2)白噪声；(3)有色随机信号。程序如下：结果如下：result1=[-0.8287;0.1275;-0.0024;1.9884;-1.2723]递推最小二乘〔在线辨识〕RLS为了启动RLS，需给初值：。计算框图见书P66。程序如下：结果如下：result2=[-0.8284;0.1274;-0.0024;1.9883;-1.2717]图7递推最小二乘法参数过渡过程数据饱和：〔1〕原因：，不再起修正作用，引起误差变大。〔2〕为了克制数据饱和现象，可以用降低老的数据影响的方法：=1\*GB3①渐消记忆法〔遗忘因子法〕=2\*GB3②限定记忆法〔固定窗法〕当为不相关序列，最小二乘有一致性与无偏性，但往往为相关序列，为克制最小二乘有偏估计的缺点，引入辅助变量法和广义最小二乘法，增广最小二乘法等。渐消记忆法〔遗忘因子法〕一般程序如下：结果如下：result3=[-0.6862;0.1114;0.0640;2.0356;-0.9429]图8渐消记忆法参数过渡过程限定记忆法〔固定窗法〕程序如下：结果如下：〔较前三种方法偏差较大〕result4=[-1.0190;0.1725;-0.4531;2.1092;-1.3724]辅助变量法〔IV〕(1).辅助变量Z.计算步骤：=1\*GB3①先根据实测数据最小二乘求粗略〔为有偏估计〕=2\*GB3②=3\*GB3③=4\*GB3④.递推：RIV初始条件：缺点:P0的选择非常敏感，一个改良方法是，用递推最小二乘辨识算法作为启动方法，然后转换到辅助变量法。程序如下：结果如下：result5=[-0.9369;0.1207;-0.0254;1.9781;-1.5331]图9辅助变量法参数过渡过程广义最小二乘法〔GLS〕(1)广义最小二乘法的根本思想：由于在n+k个采样周期的时差范围内具有自相关性，从而使的最小二乘估计为有偏的，所以引入一个所谓成形滤波器〔白化滤波器〕，把相关噪声转化成白噪声。如果知道有色噪声序列的相关性：令，有广义最小二乘法(GLS)是建立在最小二乘法(LS)的根底上的。根本最小二乘法只是广义最小二乘法在时的特例。(2)广义最小二乘法计算步骤：广义最小二乘法的关键问题是如何用比拟简单的方法找到成形滤波器的系数。其计算是逐次逼近法。=1\*GB3①应用输入输出数据按最初模型求出的最小二乘估计。这个估计值是不准确的，它只是被估参数的一次近似。=2\*GB3②计算残差e(k)，并拟合成形滤波器的模型：得到其中=3\*GB3③应用所得的成形滤波器，对输入输出数据滤波：其中，m为噪声模型的阶，一般事先不知道，实际经历说明指定m为2或3可以得到比拟满意的输出。=4\*GB3④按新的输入、输出模型求出参数的第二次估计值。结果如下：Result6=[-0.6538;0.2926;0.0454;1.8776;-1.6042]广义递推最小二乘法〔GLS〕广义最小二乘法的递推计算过程可分成两个局部：按递推最小二乘法(RLS)，随着N的增大，不断计算(逐步接近于无偏)和(逐步使噪声白化)；在递推过程中，和是时变的，则过滤信号及残差是由时变系统产生，要不断计算。因而，递推广义最小二乘法由两组普通的递推最小二乘法组成，它们是通过滤波算法联系起来的：结果如下：result7=[-1.0164;0.1754;-0.0159;2.0056;-1.6438]噪声传递系数的估计结果：[-0.0307;0.0900]图10广义递推最小二乘参数过渡过程增广矩阵法(ELS/RELS)(增广最小二乘法)增广矩阵法是把观测矩阵适当增大，使得有偏估计的程度得到一定改善。这