参数辨识示例报告

1、参数辨识参数辨识的步骤飞行器气动参数辨识是一个系统工程，包括四部分：试验设计，使试验能为辨识提供含有足够信息量且信息分布均匀的试验数据；气动模型结果确定，即从候选模型集中，根据一定的准则和经验，选出最优的气动模型构式；气动参数辨识，根据辨识准则和数据求取模型中待定参数，这是气动辨识定量研究的核心阶段；模型检验，确认所得气动模型是否确实反映了飞行器动力学系统中气动力的本质属性。这四个部分环环相扣，缺一不可，要反复进行，直到对所得气动模型满意为止。参数辨识的方法参数辨识方法主要有最小二乘算法、极大似然法、集员辨识法、贝叶斯法、岭估计法、超椭球法和鲁棒辨识法等多种辨识方法。虽然目前参数辨识的领域己经

2、发展了多种算法，但是用于气动参数估计的算法主要有：极大似然法(ML),广义Kalman滤波(EKF)法，模型估计法(EBM)、分割及多分割算法(PIA及MPIA)、最小二乘法，微分动态规划法等。因为最小二乘法和极大似然法是两种经典的算法，目前己经发展得相当成熟。最小二乘法适于线性模型的参数辨识，可以用于飞行器系统辨识中很多的线性模型，如惯性仪表误差系数的辨识，线性时变离散系统初始状态的辨识及多项式曲线拟合等。目前最小二乘法已经广泛应用于工程实际中。而极大似然算法因其具有渐进一致性、估计的无偏性、良好的收敛特性等特点而被广泛应用于飞行器参数辨识领域。最小二乘法大约是1975年高斯在其著名的星体运

3、动轨道预报研究工作中提出来的。后来，最小二乘法就成了估计理论的奠基石。由于最小二乘法原理简单，编程容易，所以它颇受人们重视，应用相当广泛。极大似然估计算法在实践中不断地被加以改进，这种改进主要表现在三个方面，即算法的优化、是否考虑过程噪声以及灵敏度的计算等。极大似然法具有集中处理数据的特点，即对一段时间历程上的数据集中分析，直接得到所需的气动参数，此方法也具有较高的精度，其估计具有一致性和无偏性。但是，它也有缺点，主要是计算所需时间较长，易于进行集中处理而不适合在线实时处理。递推最小二乘法设时不变动态过程的数学模型为A(z-1)z(k)=B(zt)u(k),v(k)式中u(k)和z(k)为系统

4、的输入量和输出量；v(k)是噪声；A(z-1)和3(z-1)是多项式，定义为：A(z-1)az-2+.+az-n1nB(z-1)bz-1+bz-2+.+bz-m12m现在问题是如何利用系统的u(k),z(k);k=1,.,L序列，估计多项式A(z-1)和B(zi)中的未知系数。定义:0二a,a，,a,b,b，,bT12n12mh(k)=z(k1),.,z(kn),u(k1),.,u(km)t则可将数学模型化为一个标准的最小二乘格式z(k)=hT(k)0,v(k)上式可写成一个线性方程组，形式如下:其中Z=z(1),z(2),.,z(LLh41)-z(0)一z(1一n)u(0)u(1-m)h2(

5、2)=-z-z(2-n)u(1)u(2-m)_hT(L)_-z(L-1)-票L禺)u(L一1)u(L-m)最小二乘的准则函数J(0)可取为:V二v(1),v(2),.,v(L)FLH=LJ()二(ZH)t(ZH)LlLl极小化J(),可求出系统系数的估计值。下面给出最小二乘估计一般估值计算公式：=(HtH)，1H1zlllL对应的方法叫做最小二乘法。随着待估计的参数和观测数据的增多，最小二乘一次完成算法的计算量增加计算困难不利于在线辨识。现在把它化成递推算法，其格式为：新的估计值(k)=老的估计值(k，1)+修正项这样的优点是：每一步的计算量比较小，因而能够用比较少的计算量完成较大的任务。具有

6、跟踪时变参数的能力，能辨识含有时变参数的系数模型。在自适应控制系统和基本诊断目的的故障检测中，需要递推算法。在参数估计达到给定的精确度时算法可以给出收敛中止的判据。示例：一模型形式为y=au+au+e1122取得3组实验数据为y(k):1,1,3;x(k):1,2,1;x(k):0,-1,212求出参数a,a的最小二乘估计。0，1，所以212该模型中，1=(a,a)T,y=(1,1,3)T,O=2121解正则方程组tQ，Ty0fa1可得参数的最小二乘估计为：a，1,a，1。12MATLAB算例通过上例我们可初步了解应用最小二乘法进行了解，下面通过MATLAB编程实现递推最小二乘法求解系统参数辨

7、识，这里选用单输入单输出系统y，au2+bu+c作为待辨识模型，要求对a,b,c进行参数辨识，首先需要取得若干组输入输出数据，故我们事先选定实际系统为y，0.2345u2-4.5678u+5.1234，通过对系统进行测量我们得到以下初始数据。u，1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,1&19,20y，0.7901,-3.0742,-6.4695,-9.3958,-11.8531,-13.8414,-15.3607,-16.4110,-16.9923,-17.1046,-16.7479,-15.9222,-14.6275,-12.8638,-10.6

8、311,-7.9294,-4.7587,-1.1190,2.9897,7.5674加入随机误差，我们得到将要使用的初始数据。u=1.012.033.024.0156.027.038.049.0310.0211.0212.0213.0214.0215.0216.0317.0218.0119.0220.04;y=0.791-3.072-6.465-9.398-11.83-13.414-15.607-16.10-16.23-17.046-16.479-15.222-14.275-12.638-10.311-7.294-4.87-1.902.977.74对该模型参数辨识采用递推最小二乘法，即RLS(r

9、ecurisiveleastsquare)，它是一种能够对模型参数进行在线实时估计的辨识方法。其基本思想可以概括为：新的估计值(k)=旧的估计值0(k-1)+修正项下面将批处理最小二乘法改写为递推形式即递推最小二乘参数估计的计算方法。批处理最小二乘估计为，(t)-1tY，设k时刻的批处理最小二乘估计为：，(T)-1tYkkkk令P(k)，(T)-1，P-1(k-1)+(k)T(k-1)-1kkK时刻的最小二乘估计可以表示为(k)，P(k)tY，P(k)TY+(k)y(k)=(k1)+K(k)y(k)-t(k)(k1)；kkk-1k-1式中K(k)，P(k)(k),因为要推导出P(k)和K(k)

10、的递推方程，因此这里使用矩阵求逆引理：设A,(A+BC),I+CA-1B均为非奇异方阵，则(A+BC)-1，A-1-A-1B(I+CA-1B)-1CA-1通过运用矩阵求逆引理，把复杂的矩阵求逆转化为标量求倒数，大大减小了计算量。P与P间的递推关系。最终得到递N1N推最小二乘参数递推估计公式如下：(k)，(k-1)+K(k)y(k)申t(k)(k-1)P(k)，I-K(k)Qt(k)P(k-1)P(k-1)申(k)1+申t(k)P(k-1)申(k)程序流程图Step1：设置初值(0)和P(0)，输入初始数据；Step2:采样当前输出y(k)、和输入u(k)Step3:利用递推公式计算K(k)、(

11、k)和P(k)；Step4:k=k+1，返回step2，继续循环。程序代码：X=1.012.033.024.0156.027.038.049.0310.0211.0212.0213.0214.0215.0216.0317.0218.0119.0220.04;Y=0.791-3.072-6.465-9.398-11.83-13.414-15.607-16.10-16.23-17.046-16.479-15.222-14.275-12.638-10.311-7.294-4.87-1.902.977.74;%实验输入数据、实验输出数据symsabc%定义待辨识参数theta=a;b;c;%theta

12、包含待辨识参数a，b，ctheta1=zeros(3,1);%对象参数初始化P=10A6*eye(3);%构造初始P阵fork=1:20%仿真步长范围1到20phi=X(k)*X(k);X(k);1;%y=aX*X+bX+c=phi*theta%theta=a；b；c；phi=X(k)*X(k);X(k);1K=P*phi/(1+phi*P*phi);%递推最小二乘法K阵的递推公式%theta的递推公式theta=theta1+K*(Y(k)-phi*theta1);P=(eye(3)-K*phi)*P;%递推最小二乘法P阵的递推公式theta1=theta;%theta的最终估计向量thet

13、a2(:,k)=theta%theta估计向量矩阵化，目的是为了%下面的plot仿真图像输出endtheta1%输出参数估计值plot(1:20,theta2)%输出参数逐步递推估计的轨迹图像xlabel(k);%设置横坐标为步长kylabel(参数估计a，b，c);%纵坐标为估计参数a，b，clegend(a,b,c);%标示相应曲线对应的参数得到结果:th&tal=0.2287-4.45754.9319参数估计变化曲线:6c,b，即估数参可以计算得到各辨识参数误差分别为：2.47%,0.438%,3.74%。考虑在初始参数中引入了随机误差，故模型参数的最小二乘估计在误差允许的范围内基本符合

14、要求。10上图为辨识出的系统与真实系统曲线比对，可观察到曲线拟合程度较好，说明此方法对系统的辨识结果是可以接受的。算例扩展通过上例我们对应用递推最小二乘法进行系统参数辨识有了一定了解，下面对此算例进行扩展，选用多输入单输出系统yau2+bu+cu+d作为待辨识模112型，要求对a,b,c,d进行参数辨识，首先需要取得若干组输入输出数据，故我们事先选定实际系统为yu2，2u+3u-4，通过对系统进行测量并加入随机误差112我们得到以下初始数据。u=1234567891011121314151617181920;1u=10101010101010101010;2y=-2.01-4.062.073.

15、9114.0219.9134.0144.0462.0876.0498.92115.943142.98164194.1220.91253.88284.1322.91356对该模型参数辨识采用递推最小二乘法。得到结果theta1=0.9984-1.93913.2436-4.3414参数估计变化曲线：,b，舟估数参可以计算得到各辨识参数误差分别为：0.16%，3.045%，8.12%，8.535%。考虑在初始参数中引入了随机误差，故模型参数的最小二乘估计在误差允许的范围内基本符合要求。由此可见，递推最小二乘法在单输入单输出和多输入单输出系统的参数辨识中，均可达到较好的辨识效果。本周总结本周对系统辨识方法进行了更深入的学习，使用了递推最小二乘法进行了系统参数辨识的学习应用，并编写了MATLAB程序对实例系统进行了辨识，得到了较为满意的结果。不过也有些许不足：参数变化曲线不够平滑，分析原因可能是因为我引入的是随机误差，而不是高斯白噪声，对结果有了一定影响，也有可能是采样数据个数的限制，