系统辨识第五章作业分析.docx

1、系统辨识第五章作业分析摘 要系统辨识是描述各种各样系统运动规律的一种方法论，是研究系统的一种有效工具。利用这个工具可以对我们要研究的系统进行定量描述。随着现代控制理论的迅速发展，系统辨识得到迅速而蓬勃发展，并已经成功运用与多种工程应用领域。但针对有色噪声干扰系统，传统的辨识方法不能得到良好的参数估计，而工程上大多系统都为有色噪声干扰系统。有色噪声干扰系统的一类系统为广义输出误差模型，本文对白噪声和有色噪声两种噪声干扰下运用最小二乘法，递推最小二乘法分别比较噪声的不同，以及进一步得出最小二乘法的适用性。进一步研究广义最小二乘法对有色噪声系统辨识的改进。关键词：系统辨识；白噪声；有色噪声；最小二乘

2、；递推最小二乘；广义最小二乘ABSTRACT the law of motion of the system identification is to describe the various system and an effective tool for the system. It also can be used to study quantitative description of the system. With the rapid development of modern control theory,the development of system identificat

3、ion has been rapid and vigorous , and has been successfully applied to a variety of engineering applications. However, the traditional identification methods cant get good parameter estimation in engineering systems for colored noise . Colored noise jamming system is for generalized output error mod

4、el.under two kinds of noise of white noise and color noise ,by using the least square method and the recursive least square method ,respectively to compare the different noise, and to further concluded that the applicability of the least-square method. To further studied the improvement of generaliz

5、ed least squares identification of colored noise system.KEY WORDS: system identification,white noise, colored noise.,least squares;,recursive least squares,the generalized least squares目录摘 要 1第一章 系统辨识 21.1辨识的定义 21.2辨识的内容和步骤 31.3辨识目的 31.4先验知识 31.5数据预处理 41.6模型结构辨识 41.7模型参数辨识 41.8模型检验 4第二章 题目详解 52.1.1最

6、小二乘法 52.1.2递推法 82.2.1最小二乘法 122.2.2递推法 132.3广义最小二乘法 15附录 15第一章 系统辨识1.1辨识的定义Zadeh对辨识的定义（1962年）辨识就是在输入和输出数据的基础上，从一组给定的模型类中，确定一个与所测系统等价的模型。Liung 对辨识的的定义（1978年）系统辩识有三个要素数据、模型类和（等价）准则。系统辩识是按照一个准则，在模型类中选择一个与数据拟合得最好的模型。Liung 认为，实际系统的复杂性很难找到一个适用的模型与之等价。因此，系统辩识的任务只是要求从输入输出数据出发，找到一个与实际系统相逼近的模型。该定义体现了逼近的观点。系统辨识

7、是在输入和输出观测的基础上，从指定的一类模型中确定一个与系统等价的模型。1.2辨识的内容和步骤简单地说，辨识就是一种从观测到的含有噪声的输入输出数据中提取数学模型的方法。根据现场的情况，辨识可以离线进行，也可以在线进行；辨识的内容主要包括四个方面：实验设计，模型结构辨识，模型参数辨识，模型检验。1.3辨识目的明确模型应用的最终目的很重要，因为它将决定模型的类型，精度要求和采用哪种辨识方法等问题。1.4先验知识对于给定的过程进行辨识之前，要通过一些手段对过程有一些了解粗略的掌握一些过程和数据，这些先验知识对实验设计将其指导作用。1.5数据预处理输入输出数据通常含有直流或交流成分，任何方法都无法消

8、除对系统精度的影响。此外，高频成分也是不利的。因此需对输入数据进行零均值化和剔除高频成分预处理。1.6模型结构辨识模型结构辨识包括模型验前结构和模型结构参数两部分。1.7模型参数辨识当模型结构确定之后，就需要进行模型参数辨识。方法很多，最基本，应用最广的是最小二乘法。但是最小二乘法也有一些重大缺陷，比如过程辨识受到有色噪声污染时，他几乎比能适应。1.8模型检验模型检验是辨识过程不可或缺的步骤之一。但是，它没有一般方法可循。它和模型结构问题密切相关。如果模型结构不合适，模型检验是不能通过的。第二章 题目详解2.1.1最小二乘法 一个单输入单输出线性定常系统的差分方程 式中为输入信号，为理论上的输

9、出值。通过观测得到，在观测过程中往往附加有随机干扰。观测值用表示： （k）=v（k）+则；如果 u（k）也有测量误差，则在（k）中应包含这一测量误差。 则 原题设为二阶系统，且的真值已知，输入已知，求，并且由输入输出序列得出参数估值，其程序已在matlab中实现。实验结果如下图：参数a1a2b1b2真实值1.6420.7150.390.35估计值1.60080.64830.41430.3309实验结果分析：由于所用的输出观测值有白噪声成分，所以辨识结果有误差。2.1.2递推法 图1.2 最小二乘递推算法辨识的Malab6.0程序流程图Matlab仿真结果：c = Columns 1 throu

10、gh 12 1.6420 0 1.6420 0 0 0 0 0 0 0 0 0 0.7150 0 0.7150 0 0 0 0 0 0 0 0 0 0.3900 0 0.3689 0 0 0 0 0 0 0 0 0 0.3500 0 0.2298 0 0 0 0 0 0 0 0 0 Columns 13 through 24 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Columns 25 through 30 0 0 0 0 0 0 0 0 0

11、0 0 0 0 0 0 0 0 0 0 0 0 0 0 0e = Columns 1 through 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0540 0 0 0 0 0 0 0 0 0 0 0 -0.3433 0 0 0 0 0 0 0 0 0 Columns 13 through 24 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Columns 25 through 3

12、0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0结果分析：仿真结果表明，大约递推到第五步时，参数辨识的结果基本达到稳定状态。此时，参数的相对变化量E0.000000005。从整个辨识过程来看，精度的要求直接影响辨识的速度。虽然最终的精度可以达到很小，但开始阶段的相对误差会很大。2.2.1最小二乘法题设将白噪声改为有色噪声，；有色噪声可以看做是白噪声经过成型滤波器后得到的。Matlab仿真结果：参数a1a2b1b2真实值1.6420.7150.390.35估计值1.28750.37020.34620.2773结果分析：由于所用的输出观测值有有色噪

13、声成分，所以辨识结果有误差，且用最小二乘法计算得到的参数估计误差较大。2.2.2递推法结果：结果分析：第一问和第二问的区别只是噪声不同，但是从仿真结果图像来看，都是从第五步开始达到稳定状态的，但有色噪声参数辨识显然误差更大更明显。2.3广义最小二乘法 广义最小二乘法辨识的计算步骤如下：（1）应用已得到的输入和输出数据 u(k) 和y(k) (2)计算残差e(k),并用其代替（3）计算和和 (4)应用得到的和，再用最小二乘 法重新估计，得的第2次估值 。然后按步骤计算残差 e(k) 。重新估计，得到估值 。再按步骤计算和，按步骤求的第3次估计。重复上循环步骤，直到 的估值 收敛为止。附录程序%1

14、.1 二阶系统的最小二乘一次完成算法辨识程序，在光盘中的文件名：FLch3LSeg1.mu=1.147,0.201,-0.787,-1.589,-1.052,0.866,1.152,1.573,0.626,0.433,-0.958,0.81,-0.044,0.947,-1.474,-0.719,-0.086,-1.099,1.45,1.151,0.485,1.633,0.043,1.326,1.706,-0.340,0.89,1.444,1.177; %系统辨识的输入信号为一个周期的M序列y=zeros(1,30); %定义输出观测值的长度 w=0.1*randn(1,30) for k=3:

15、30 y(k)=-1.642*y(k-1)-0.715*y(k-2)+0.39*u(k-1)+0.35*u(k-2)+w(k); %用理想输出值作为观测值endsubplot(3,1,1) %画三行一列图形窗口中的第一个图形stem(u) %画输入信号u的径线图形subplot(3,1,2) %画三行一列图形窗口中的第二个图形i=1:1:30; %横坐标范围是1到16，步长为1plot(i,y) %图形的横坐标是采样时刻i, 纵坐标是输出观测值z, 图形格式为连续曲线subplot(3,1,3) %画三行一列图形窗口中的第三个图形stem(y),grid on %画出输出观测值z的径线图形，并

16、显示坐标网格u,y %显示输入信号和输出观测信号%L=30%数据长度HL=-y(2) -y(1) u(2) u(1);-y(3) -y(2) u(3) u(2);-y(4) -y(3) u(4) u(3);-y(5) -y(4) u(5) u(4);-y(6) -y(5) u(6) u(5);-y(7) -y(6) u(7) u(6);-y(8) -y(7) u(8) u(7);-y(9) -y(8) u(9) u(8);-y(10) -y(9) u(10) u(9);-y(11) -y(10) u(11) u(10);-y(12) -y(11) u(12) u(11);-y(13) -y(1

17、2) u(13) u(12);-y(14) -y(13) u(14) u(13);-y(15) -y(14) u(15) u(14);-y(16) -y(15) u(16) u(15);-y(17) -y(16) u(17) u(16);-y(18) -y(17) u(18) u(17);-y(19) -y(18) u(19) u(18);-y(20) -y(19) u(20) u(19);-y(21) -y(20) u(21) u(20);-y(22) -y(21) u(22) u(21);-y(23) -y(22) u(23) u(22);-y(24) -y(23) u(24) u(23)

18、;-y(25) -y(24) u(25) u(24);-y(26) -y(25) u(26) u(25);-y(27) -y(26) u(27) u(26);-y(28) -y(27) u(28) u(27);-y(29) -y(28) u(29) u(28) %给样本矩阵HL赋值Y=y(3);y(4);y(5);y(6);y(7);y(8);y(9);y(10);y(11);y(12);y(13);y(14);y(15);y(16);y(17);y(18);y(19);y(20);y(21);y(22);y(23);y(24);y(25);y(26);y(27);y(28);y(29);y(

19、30) % 给样本矩阵z L赋值%Calculating Parametersc1=HL*HL; c2=inv(c1); c3=HL*Y; c=c2*c3 %计算并显示a1=c(1), a2=c(2), b1=c(3),b2=c(4) %从 中分离出并显示a1 、a2、 b1、 b2%1.2 RLS递推最小二乘辨识c0=1.642 0.715 0.39 0.35; %被辨识参数的初始值p0=106*eye(4,4); %直接给出初始状态P0，即一个充分大的实数单位矩阵E=0.000000005; %取相对误差E=0.000000005c=c0,zeros(4,29); %被辨识参数矩阵的初始值

20、及大小e=zeros(4,30); %相对误差的初始值及大小u=1.147,0.201,-0.787,-1.589,-1.052,0.866,1.152,1.573,0.626,0.433,-0.958,0.81,-0.044,0.947,-1.474,-0.719,-0.086,-1.099,1.45,1.151,0.485,1.633,0.043,1.326,1.706,-0.340,0.89,1.444,1.177; %系统辨识的输入信号为一个周期的M序列y=zeros(1,30); %定义输出观测值的长度 w=0.5*randn(1,30); for k=3:30 y(k)=-1.64

21、2*y(k-1)-0.715*y(k-2)+0.39*u(k-1)+0.35*u(k-2)+w(k); %用理想输出值作为观测值endfor k=3:30; %开始求K h1=-y(k-1),-y(k-2),u(k-1),u(k-2); x=h1*p0*h1+1; x1=inv(x); %开始求K(k) k1=p0*h1*x1;%求出K的值 d1=y(k)-h1*c0; c1=c0+k1*d1; %求被辨识参数c e1=c1-c0; %求参数当前值与上一次的值的差值 e2=e1./c0; %求参数的相对变化 e(:,k)=e2; %把当前相对变化的列向量加入误差矩阵的最后一列 c0=c1; %

22、新获得的参数作为下一次递推的旧参数 c(:,k)=c1; %把辨识参数c 列向量加入辨识参数矩阵的最后一列 p1=p0-k1*k1*(h1*p0*h1+1); %求出 p(k)的值 p0=p1; %给下次用 if e2=E break; %如果参数收敛情况满足要求，终止计算 end %小循环结束end %大循环结束c,e %显示被辨识参数及其误差(收敛)情况%分离参数a1=c(1,:); a2=c(2,:); b1=c(3,:); b2=c(4,:); ea1=e(1,:); ea2=e(2,:); eb1=e(3,:); eb2=e(4,:); figure(2); %第二个图形i=1:30

23、; %横坐标从1到30plot(i,a1,r,i,a2,:,i,b1,g,i,b2,:) %画出a1，a2，b1，b2的各次辨识结果title(Parameter Identification with Recursive Least Squares Method) %图形标题figure(3); %第三个图形i=1:30; %横坐标从1到30plot(i,ea1,r,i,ea2,g,i,eb1,b,i,eb2,r:) %画出a1，a2，b1，b2的各次辨识结果的收敛情况title(Identification Precision) %图形标题%2.1 二阶系统的最小二乘一次完成算法辨识程序，

24、在光盘中的文件名：FLch3LSeg1.mu=1.147,0.201,-0.787,-1.589,-1.052,0.866,1.152,1.573,0.626,0.433,-0.958,0.81,-0.044,0.947,-1.474,-0.719,-0.086,-1.099,1.45,1.151,0.485,1.633,0.043,1.326,1.706,-0.340,0.89,1.444,1.177; %系统辨识的输入信号为一个周期的M序列y=zeros(1,30); %定义输出观测值的长度 w=0.1*randn(1,30); for k=3:30 W(k)=w(k)+1.642*w(k

25、-1)+0.715*w(k-2) y(k)=-1.642*y(k-1)-0.715*y(k-2)+0.39*u(k-1)+0.35*u(k-2)+W(k); %用理想输出值作为观测值endsubplot(3,1,1) %画三行一列图形窗口中的第一个图形stem(u) %画输入信号u的径线图形subplot(3,1,2) %画三行一列图形窗口中的第二个图形i=1:1:30; %横坐标范围是1到16，步长为1plot(i,y) %图形的横坐标是采样时刻i, 纵坐标是输出观测值z, 图形格式为连续曲线subplot(3,1,3) %画三行一列图形窗口中的第三个图形stem(y),grid on %画

26、出输出观测值z的径线图形，并显示坐标网格u,y %显示输入信号和输出观测信号%L=30%数据长度HL=-y(2) -y(1) u(2) u(1);-y(3) -y(2) u(3) u(2);-y(4) -y(3) u(4) u(3);-y(5) -y(4) u(5) u(4);-y(6) -y(5) u(6) u(5);-y(7) -y(6) u(7) u(6);-y(8) -y(7) u(8) u(7);-y(9) -y(8) u(9) u(8);-y(10) -y(9) u(10) u(9);-y(11) -y(10) u(11) u(10);-y(12) -y(11) u(12) u(1

27、1);-y(13) -y(12) u(13) u(12);-y(14) -y(13) u(14) u(13);-y(15) -y(14) u(15) u(14);-y(16) -y(15) u(16) u(15);-y(17) -y(16) u(17) u(16);-y(18) -y(17) u(18) u(17);-y(19) -y(18) u(19) u(18);-y(20) -y(19) u(20) u(19);-y(21) -y(20) u(21) u(20);-y(22) -y(21) u(22) u(21);-y(23) -y(22) u(23) u(22);-y(24) -y(2

28、3) u(24) u(23);-y(25) -y(24) u(25) u(24);-y(26) -y(25) u(26) u(25);-y(27) -y(26) u(27) u(26);-y(28) -y(27) u(28) u(27);-y(29) -y(28) u(29) u(28) %给样本矩阵HL赋值Y=y(3);y(4);y(5);y(6);y(7);y(8);y(9);y(10);y(11);y(12);y(13);y(14);y(15);y(16);y(17);y(18);y(19);y(20);y(21);y(22);y(23);y(24);y(25);y(26);y(27);

29、y(28);y(29);y(30) % 给样本矩阵z L赋值%Calculating Parametersc1=HL*HL; c2=inv(c1); c3=HL*Y; c=c2*c3 %计算并显示a1=c(1), a2=c(2), b1=c(3),b2=c(4) %从 中分离出并显示a1 、a2、 b1、 b2%2.2 RLS递推最小二乘辨识c0=1.642 0.715 0.39 0.35; %被辨识参数的初始值p0=106*eye(4,4); %直接给出初始状态P0，即一个充分大的实数单位矩阵E=0.000000005; %取相对误差E=0.000000005c=c0,zeros(4,29)

30、; %被辨识参数矩阵的初始值及大小e=zeros(4,30); %相对误差的初始值及大小u=1.147,0.201,-0.787,-1.589,-1.052,0.866,1.152,1.573,0.626,0.433,-0.958,0.81,-0.044,0.947,-1.474,-0.719,-0.086,-1.099,1.45,1.151,0.485,1.633,0.043,1.326,1.706,-0.340,0.89,1.444,1.177; %系统辨识的输入信号为一个周期的M序列y=zeros(1,30); %定义输出观测值的长度 w=0.5*randn(1,30); for k=3:30 W(k)=w(k)+1.642*w(k-1)+0.715*w(k-2) y(k)=-1.642*y(k-1)-0.715*y(k-2)+0.39*u(k-1)+0.35*u(k-2)+W(k); %用理想输出值作为观测值endfor k=3:30; %开始求K h1=-y(k-1),-y(k-2),u(k-1),u(k-2); x=h1*p0*h1+1; x1=inv(x); %开始求K(k) k1=p0*h1*x1;%求出K的值 d1=y(k)-h1*c0; c1=c0+k1*d