未知初始条件下滞后过程的模型参数估计方法报告.pdfVIP

Equipment ManufacturingTechnologyNo．8，2015 未知初始条件下滞后过程的模型参数估计方法 曾雪芬 (江汉大学物理与信息工程学院，湖北 武汉 430000) 摘 要 ：提出了一种参数估计过程中消去未知初始条件的方法。在对象输入输出微分方程两边乘以时间变量，依据分部 积分原理，经过多重积分，估计方程中不再含有未知初始条件。通过优化幅值均差函数辨识对象的滞后时间，然后采用最 小二乘法辨识得到其他模型参数。仿真结果表明算法具有较高精度。 关键词：未知初始条件 ；滞后时间；参数估计；最小二乘法 中图分类号：0213 文献标识码：A 文章编号：1672—545X(2015)08—0008—02 大多数工业生产过程都会在不同程度上存在着 分部积分原理 ，经过与模型阶数相同阶次的多重积 滞后 ”【。通过参数估计方法，利用实验或现场运行数 分运算后 ，消去估计方程中的未知初始条件。根据输 据辨识这些过程的滞后模型参数，可以为后续的控 入输出数据的幅值均差函数辨识对象的滞后时间， 制器设计、控制系统性能分析提供重要依据。 然后利用过程输入输出数据构造出估计方程，通过 文献 2【】对 比分析了离散时间序列的互相关 函 最小二乘法辨识得到模型的其他参数。 数、均方差函数和幅值均差函数，求解这些函数的极 值点都可以找出滞后时间。其中，幅值均差函数的 1 辨识方法 计算量最省。在滞后时间不是采样时间的整数倍 时，可以利用抛物线曲线近似描述函数极值点附近 1．1滞后参数的辨识 的图形，求解抛物线曲线的极值点，其对应的时间点 对过程的输入输出量 u(t)、y(t)进行周期采样 ， 作为滞后时间的估计值。文献[3】分析了通过多重积 得到离散时间序列u(kT,)、y(kT!)， 为采样周期。定 分运算估计模型参数的问题。通过对输人输出微分 义幅值均差函数 ： 方程进行多重积分运算 ，将参数估计问题转化为对 ， 上 以过程模型参数为变量、以输出响应数据为参数的 JR(下)=』 k = J 1( )一y(kT!+r)1 代数方程的求解。由于只使用了输出响应的少量数 其中Ⅳ为采样点数量。与互相关函数相比，幅值 据，参数的估计结果很容易受到噪声干扰。文献[4】 均差函数不需要乘法计算，减少了计算量。 利用大量数据通过积分运算构造出估计方程，然后 输入输出之间的滞后时间6通过优化R (丁)得 采用最小二乘法估计得到模型参数，增强了参数估 到，即6=argmax(-R ())．求解滞后时间时，预先 计算法的抗干扰能力。文献 5【】针对二阶加滞后模型 估计出延迟时间丁的大致范围，即rnn／≤r≤ 一，从中 的参数估计方案 ，提出了一种消去未知初始条件的 找出使R (下)达到最大值的延迟时间即为滞后时间。 方法。在经过初次积分运算得到的估计方程两边同 当滞后时间不是采样周期的整数倍时，即6≠noT,， 时乘以一个满足一定条件的权值函数，经过再次积 (no一1) 8(no+1) ，noT!是R (丁)最大值对应 分之后 ，可以消去估计方程中出现的未知初始条件 的延迟时间，利用极值点附近的三点R ((n。一1))、 部分。增加了权值函数之后，采样点与采样间隔都 R (n。)、R ((11,。+1) )、构造出抛物线JR()=6t3 需要按特定要求选取，并且权值函数也会对信号频 +6丁+c，近似描述最大值附近的图形，然后用抛物 谱产生影响。 线顶点 (丁)=一b／2a对应的延迟时间作为滞后时 本文提出了一种在参数估计过程中消去未知初 间的估计值。R (丁)存在最大值的前提条件是输入输 始条件的方法。通过在对象的输入输出微分方程两 出数据必须包含波分或包谷 ，因此 ，实际应用时需选 端乘以与模型阶数相同阶次的时间变量 ，然后利用 择合适的输入信号。 收稿 日期：