系统辨识实验报告.doc

《系统辨识实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.专业整理.一、相关分析法(1)实验原理图1实验原理图本实验的原理图如图1。过程传递函数中；输入变量，输出变量，噪声服从，为过程的脉冲响应理论值，为过程脉冲响应估计值，为过程脉冲响应估计误差。过程输入采用M序列，其输出数据加白噪声得到输出数据。利用相关分析法估计出过程的脉冲响应值，并与过程脉冲响应理论值比较，得到过程脉冲响应估计误差值。M序列阶次选择说明：首先粗略估计系统的过渡过程时间TS(通过简单阶跃响应)、截止频率fM(给系统施加不同周期的正弦信号或方波信号，观察输出)。本次为验证试验，已知系统

2、模型，经计算，。根据式及式，则取值为1，此时，由于与N选择时要求完全覆盖，则选择六阶M移位寄存器，即N=63。.学习帮手..专业整理.(2)编程说明人机对话噪声标准差：sigma；生成数据周期数：r生成数据生成M序列；生成白噪声序列过程仿真得到理论输出数据计算脉冲响应估计值计算互相关函数，得到脉冲响应估计值计算脉冲响应估计误差计算脉冲响应理论值，得到脉冲响应估计误差图2程序流程图(3)分步说明①生成M序列：M序列的循环周期，时钟节拍，幅度，移位寄存器中第5、6位的内容按“模二相加”，反馈到第一位作为

3、输入。其中初始数据设为{1,0,1,0,0,0}。程序如下：.学习帮手..专业整理.②生成白噪声序列：程序如下：③过程仿真得到输出数据：如图2所示的过程传递函数串联，可以写成形如，其中。图2过程仿真方框图程序如下：④计算脉冲响应估计值：.学习帮手..专业整理.互相关函数采用公式，互相关函数所用的数据是从第二个周期开始的，其中为周期数，取1-3之间。则脉冲响应估计值为：，。补偿量。程序如下：⑤计算脉冲响应估计值：脉冲响应的理论值由式可计算得到。这时可得到过程脉冲相应估计误差。脉冲响应估计误差为：程序如

4、下：.学习帮手..专业整理.(2)数据记录①当噪声标准差sigma=0.1，生成数据周期r为2时：脉冲响应估计误差为0.0121。脉冲响应估计曲线为图3所示。②当噪声标准差sigma=0.5，生成数据周期r为1时：脉冲响应估计误差为0.0347。脉冲响应估计曲线为图4所示。③当噪声标准差sigma=0.5，生成数据周期r为3时：脉冲响应估计误差为0.0258。脉冲响应估计曲线为图5所示。④当噪声标准差sigma=1，生成数据周期r为3时：脉冲响应估计误差为0.0279。脉冲响应估计曲线为图6所示。图

5、3sigma=0.1，r=2时脉冲响应估计曲线图4sigma=0.5，r=1时脉冲响应估计曲线.学习帮手..专业整理.图5sigma=0.5，r=3时脉冲响应估计曲线图6sigma=1，r=3时脉冲响应估计曲线(2)结果分析实验中可以看到脉冲响应估计的曲线与理论曲线的重合度还是比较高的，脉冲响应估计误差也比较小，实验证明相关分析法的估计效果还是不错的。同时，经过实验可以得出结论：固定数据周期r，给定不同的噪声标准差sigma可以发现，噪声的方差越大，也就是信噪比越大，估计的效果越不好；固定噪声标准差

6、sigma，选择不同的数据生成周期r可以发现，数据周期越大，估计的周期越多，估计的效果越好。二、最小二乘法1.基本最小二乘（离线辨识）残差为：最小二乘目标：残差平方和最小（一阶导为0，二阶导>0）。从上式看出，逆存在才有解，满足条件的u(k)：(1)伪随机；(2).学习帮手..专业整理.白噪声；(3)有色随机信号。程序如下：结果如下：result1=[-0.8287;0.1275;-0.0024;1.9884;-1.2723]1.递推最小二乘（在线辨识）RLS为了启动RLS，需给初值：。计算框图见书

7、P66。程序如下：.学习帮手..专业整理.结果如下：result2=[-0.8284;0.1274;-0.0024;1.9883;-1.2717]图7递推最小二乘法参数过渡过程数据饱和：（1）原因：，不再起修正作用，引起误差变大。（2）为了克服数据饱和现象，可以用降低老的数据影响的方法：①渐消记忆法（遗忘因子法）②限定记忆法（固定窗法）当为不相关序列，最小二乘有一致性与无偏性，但.学习帮手..专业整理.往往为相关序列，为克服最小二乘有偏估计的缺点，引入辅助变量法和广义最小二乘法，增广最小二乘法等。1

8、.渐消记忆法（遗忘因子法）一般程序如下：结果如下：result3=[-0.6862;0.1114;0.0640;2.0356;-0.9429]图8渐消记忆法参数过渡过程2.限定记忆法（固定窗法）.学习帮手..专业整理.程序如下：结果如下：（较前三种方法偏差较大）result4=[-1.0190;0.1725;-0.4531;2.1092;-1.3724]1.辅助变量法（IV）(1).辅助变量Z(2).计算步骤：①先根据实测数据最小二乘求粗略（为有偏估计）②③④(3).