一种基于卡尔曼滤波的灰色模型热误差数据处理方法技术

本发明专利技术涉及一种基于卡尔曼滤波的灰色模型热误差数据处理方法。首先采用卡尔曼滤波对实验所测得的热误差数据进行滤波处理，然后建立灰色热误差模型，并根据滤波后的热误差数据利用灰色热误差模型进行预测。实验结果表明，与直接对热误差实验数据进行预测的现有方法相比，本发明专利技术所述方法获得的热误差预测值与实测值吻合得很好，可消除或减弱实验所获得的热误差数据中含有的系统及测量噪声等因素的影响，从而获得更加准确的热误差预测值，对精密卧式加工中心高速电主轴热误差建模和热误差补偿具有重要意义。

【技术实现步骤摘要】
【专利摘要】本专利技术涉及。首先采用卡尔曼滤波对实验所测得的热误差数据进行滤波处理，然后建立灰色热误差模型，并根据滤波后的热误差数据利用灰色热误差模型进行预测。实验结果表明，与直接对热误差实验数据进行预测的现有方法相比，本专利技术所述方法获得的热误差预测值与实测值吻合得很好，可消除或减弱实验所获得的热误差数据中含有的系统及测量噪声等因素的影响，从而获得更加准确的热误差预测值，对精密卧式加工中心高速电主轴热误差建模和热误差补偿具有重要意义。【专利说明】
本专利技术涉及数控机床误差补偿技术，具体涉及一种基于卡尔曼滤波的精密卧式加工中心高速电主轴的灰色模型热误差数据处理方法。
技术介绍
热误差是指由机床温度升高引起机床零部件热变形，并且导致工件与刀具间相对位置变化而产生的加工误差。有关热误差的研究内容主要有热误差的机理分析、热误差补偿中的测量、温度关键点优化、热误差建模和热误差补偿实施五大部分。热误差是数控机床最大的误差源，也是加工精度重要的影响因素。因此，必须合理有效地控制因温度升高导致的热误差，对热误差补偿中所测得的热误差数据的测量噪声进行相应的处理，以此建立精度较高的热误差模型，是实现热误差补偿和提高机床加工精度重要的关键技术。随着数控机床的不断发展，实现机床的高速精密加工，这就需要较高的机床加工精度，而影响机床加工精度的误差项不再是简单的导轨的直线度和滚珠丝杠的误差，主要是高速电主轴以及其他各部件的发热导致的热变形误差。近年来国内外学者对于高速电主轴热误差进行了大量的研究。其中，在对所使用的热误差数据的测量方面，主要使用了诸如位移传感器和温度传感器、球杆仪、坐标测量机、激光干涉仪等精密测量仪器作为机床加工精度的检测工具。关于热误差数据的滤波技术有很多。2013年6月刘志峰等发表在《高技术通讯》上的论文“基于灰色线性回归组合模型的机床热误差建模方法”，公开了一种热误差数据处理技术，该技术采用灰色线性回归组合模型对所测得的热误差原始数据进行预测，并利用BP神经网络对该模型的热误差残差进行了修正，取得了一定的预测效果。但是，该技术只是直接采用灰色线性回归组合模型对实验所测得的热误差数据进行建模和预测，并没有考虑热误差数据的前期处理，这样会减弱热误差数据的预测效果。因为这些实验所测得的热误差数据有测量噪声、环境等因素的影响，需要对所测得的热误差数据进行前期处理，来获得更好的热误差的预测值。
技术实现思路
针对现有技术中存在的因对实验所测量的热误差数据缺少前期处理而造成热误差预测值误差较大的问题，本专利技术提供一种基于卡尔曼滤波的高速电主轴灰色模型的热误差数据处理方法，通过卡尔曼滤波方法获得比较准确的热误差测量值，然后采取灰色热误差模型对滤波后的热误差数据做进一步处理，获得更加精确的热误差预测值。为了实现上述目的，本专利技术采用以下技术方案:采用卡尔曼滤波对实验所测得的热误差数据进行滤波处理，然后再根据滤波后的热误差数据利用灰色热误差模型进行预测。灰色热误差模型是一种通过少量样本热误差数据建立的预测模型。其基本思想是通过热误差数据组成原始序列，经过累加生成法生成另一种热误差数据序列，这样可以减弱热误差原始数据的随机性。通过生成序列的热误差数据建立微分方程模型，从而获得具有较强规律性的热误差预测值。基于卡尔曼滤波的灰色模型热误差数据处理方法，包括下列步骤:步骤一，对实验测量所获得的热误差数据进行卡尔曼滤波处理。(I)由热误差的差分方程和测量方程获得热误差测量值。热误差的差分方程和测量方程分别为:【权利要求】1.，其特征在于，首先对热误差测量值进行卡尔曼滤波，然后采取灰色热误差模型对滤波后的热误差数据做进一步处理，获得更加精确的热误差预测值；所述方法包括以下步骤: 步骤一，对实验测量所获得的热误差数据进行卡尔曼滤波处理； (1)由热误差的差分方程和测量方程获得热误差测量值； 热误差的差分方程和测量方程分别为: 【文档编号】G05B19/404GK103885386SQ201410139067【公开日】2014年6月25日 申请日期:2014年4月8日 优先权日:2014年4月8日 【专利技术者】刘志峰, 潘明辉, 张爱平, 方翠, 赵永胜, 蔡力钢 申请人:北京工业大学本文档来自技高网...

【技术保护点】
一种基于卡尔曼滤波的灰色模型热误差数据处理方法，其特征在于，首先对热误差测量值进行卡尔曼滤波，然后采取灰色热误差模型对滤波后的热误差数据做进一步处理，获得更加精确的热误差预测值；所述方法包括以下步骤：步骤一，对实验测量所获得的热误差数据进行卡尔曼滤波处理；（1）由热误差的差分方程和测量方程获得热误差测量值；热误差的差分方程和测量方程分别为：Ay(k)＝Bu(k)+w(k)z(k)＝Cy(k)+v(k)其中，A＝1+a1x‑1+a2x‑2+…+anx‑n，B＝1+b1x‑1+b2x‑2+…+bnx‑n，a1,a2,…,an,b1,b2,…,bn为系数，x为多项式中的未知标量元，n为多项式的次数，且为不小于1的正整数，y(k)是热误差估计值，u(k)是传感器带来的测量噪声，w(k)是系统噪声，z(k)是热误差测量值，C是与热误差估计值相对应的单元矩阵，v(k)是测量噪声，k＝1,2,…,N，表示数据元素的次序，N为k的最大值；将热误差测量值写成如下形式：φ＝[a1,a2,…,an,b1,b2,…,bn]T其中，为测量噪声序列，φ为参数矩阵；（2）利用最小二乘法估计误差测量方程中的参数；根据热误差数据z(k)及测量噪声序列构造热误差数据向量矩阵Z及含有测量噪声的向量矩阵Ψ分别为：Ζ＝[z(1),z(2),…,z(k)]TΨ=-z(0),-z(-1),...,-z(1-n),u(0),u(1),...,u(1-n)-z(1),-z(0),...,-z(2-n),u(1),u(2),...,u(2-n)......-z(k-1),-z(k-2),...,-z(k-n),u(k-1),u(k-2),...,u(k-n)]]>其中，n＝1,2,…,N；利用最小二乘法估计求解参数矩阵φ的估计值：φ^=(ΨTΨ)-1ΨTZ]]>（3）求热误差测量值的卡尔曼滤波值；卡尔曼滤波的递推公式为：x^(k|k-1)=Ax^(k-1|k-1)+Bu(k)]]>P(k|k‑1)＝AP(k‑1)AT+BQBTK(k)＝P(k|k‑1)CT[CP(k|k‑1)CT+R]‑1x^(k|k)=x^(k|k-1)+K(k)[z(k)-Cx^(k|k-1)]]]>P(k)＝[I‑K(k)C]P(k|k‑1)其中，Q、R分别表示系统噪声和测量噪声的协方差即Q＝Var(w)，R＝Var(v)，表示状态一步预测估计值，表示状态预测估计值，P(k|k‑1)表示预测的误差方差阵，K(k)为滤波增益；当测量噪声v(k)较小时，利用最小二乘算法求解热误差测量方程参数可以获得较理想的结果；随着测量噪声v(k)的增加，滤波效果将会变差；先滤掉热误差数据测量值z(k)中的噪声，然后再用于热误差数据处理，将会取得较好的效果；可通过均值法滤掉z(k)中的噪声，即用y(k)的估计值代替测量值z(k)，其中，E表示求数学期望，将y(k)的估计值代入差分方程和测量方程得：Ay(k)~=Bu(k)+Av(k)~+w(k)]]>其中，因此，热误差数据预测精度会得到提高；步骤二，建立灰色热误差模型，对滤波后的热误差数据作进一步处理，获得热误差的预测值；（1）求紧邻均值生成序列；设X(0)＝(x(0)(1),x(0)(2),…,x(0)(n))是一个卡尔曼滤波后热误差数据序列，X(1)为X(0)的累加序列，X(1)＝(x(1)(1),x(1)(2),…,x(1)(n))，其中，令z(1)(k)=12(x(1)(k)+x(1)(k-1)),k=2,3,...,n,]]>则紧邻均值生成序列为：z(1)＝(z(1)(2),z(1)(3),…,z(1)(k))（2）建立灰色热误差模型，由模型求热误差预测值；灰色热误差模型的微分方程为：dx(1)dt+ax(1)=b]]>式中，a为发展系数，b为灰色作用量；假设为参数列，利用灰色模型方程x(0)(k)+az(1)(k)＝b，通过最小二乘估计算法得到参数列a^=(BTB)-1BTY]]>其中，Y=x(0)(2)x(0)(3)···x(0)(n),B=-z(1)(2)1-z(1)(3)1······-z(1)(n)1]]>由灰色热误差模型的微分方程求得时间响应序列为：x^(1)(k+1)=(x(0)(1)-ba)e-ak+ba,k=1,2,...,n]]>通过时间响应序列作累减运算得到还原值：x^(0)(k+1)=x^(1)(k+1)-x^(1)(k)]]>从而得到灰色热误差模型的预测方程：x(0)(k+1)=(1-ea)(x(0)(1)-ba)e-ak=...

【技术特征摘要】

【专利技术属性】
技术研发人员：刘志峰，潘明辉，张爱平，方翠，赵永胜，蔡力钢，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11