资源描述：
《递推最小二乘法的应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、1递推最小二乘法在电厂模型辨识中的应用电厂中大多数热工对象可以用一阶或二阶有迟延和非迟延的模型来表示，对这些模型中参数的辨识，递推最小二乘法是一种较好的方法。本文以火电厂部分典型一阶模型为例子，借助于某电厂现场数据，分别对以下几种环节进行辨识。1.1一阶惯性环节火电厂中，来自锅炉的过热蒸汽，经高压调节汽门和导汽管道进入高压缸膨胀做功，高压缸的排汽回到锅炉再热器被重新加热，加热后的蒸汽经中压调节汽门进入中低压缸进一步膨胀做功，做功后的乏汽最终排入凝汽器变成凝结水，一般中压调节汽门的开度是高压调节汽门的3倍，即在机组负荷大于额定的30%或者滑压运
2、行时，汽轮机的中压调门是完全开启的。因此，在简化模型中，汽机侧调速器一级压力与机组有功功率可以简化为一阶惯性环节如下：将以上环节离散化，并写成差分方程的形式其中u为调速器一级压力，y为机组有功功率，为零均值方差为1的高斯白噪声。该论文依据递推最小二乘法原理，借助MATLAB工具编写程序，设定合适的初始值和加权因子进行参数辨识，辨识结果为，由可得到，进而得到系统的传递函数为：下面运用递推最小二乘法对所得结果进行仿真：假设，已知，采样时间为，则计算可得用M序列作为输入信号，得到输出信号，然后对参数进行辨识。程序如下：%最小二乘的递推算法%Z(k+
3、1)=0.9191*Z(k)+2.7670*u(k)+v(k)%========================================clearclc%==========400个产生M序列作为输入===============x=[010110111];%initialvaluen=403;%n为脉冲数目M=[];%存放M序列fori=1:ntemp=xor(x(4),x(9));M(i)=x(9);forj=9:-1:2x(j)=x(j-1);endx(1)=temp;end%===========产生均值为0，方差为1的高斯白
4、噪声=========v=randn(1,400);%==============产生观测序列z=================z=zeros(402,1);z(1)=-1;fori=2:401z(i)=0.9191*z(i-1)+2.7670*M(i-1)+v(i-1);end%递推求解P=10*eye(2);%协方差阵Theta=zeros(2,401);%参数的估计值，存放中间过程估值Theta(:,1)=[0.001;0.001];%参数估计值初始值a1=0.001,b1=0.001u=0.99;%遗忘因子fori=2:401h=[
5、-z(i-1);M(i-1)];K=P*h/(h'*P*h+u);%增益矩阵Theta(:,i)=Theta(:,i-1)+K*(z(i)-h'*Theta(:,i-1));P=(eye(2)-K*h')*P/u;end%==========================输出结果及作图=============================disp('参数a1b1的估计值:')Theta(:,401)a1=zeros(1,401);a1(1,:)=-0.9191;b1=zeros(1,401);b1(1,:)=2.7670;i=1:401
6、;figure(1)plot(i,Theta(1,:),'m',i,a1,'g',i,Theta(2,:),'b',i,b1,'r');e_a1=Theta(1,:)-a1;%a1的误差e_b1=Theta(2,:)-b1;%b1的误差figure(2)plot(i,e_a1,'m',i,e_b1,'b');%误差曲线图递推算法仿真曲线如图1所示：图1递推算法仿真曲线误差曲线如图2所示：图2误差曲线得到400个数据的a1b1估计值为:由估计值可得，由递推算法仿真曲线图1可知，参数估计收敛速度较快，占用内存较少。由误差曲线图2可知，a1辨识效果
7、较好，b1辨识误差较大。1.2一阶惯性迟延环节炉内燃烧与传热过程可以简化为磨煤机动态和水冷壁动态两个部分，在此我们将两者合并为一个一阶惯性迟延环节，研究表明此简化也能较好反映锅炉传热过程：式中，B为炉膛的燃料量，为已知量；D为锅炉总有效吸热量；τ为纯迟延时间；K2、T2为传递函数的系数。针对纯迟延的参数应用改进的具有最小损失函数的递推最小二乘法辨识，具体辨识算法为：B为已知量，令：对τ、K2、T2进行辨识。假设采样时间，。将以上环节离散化，并写成差分方程的形式相对于递推最小二乘法h（k）中多了一个时间常数d，其中，参数h（k）是各个参数的函数
8、，输入输出观测向量是纯迟延时间d的函数，残差e（k）取决于各参数的拟合误差，最小二乘估计是使目标函数J（k，d）=min［J（k，d）］最小，由于d为离散值，当已知
