发明内容
本发明是为避免上述现有技术所存在的不足之处,提供一种应用简便、建模容易、稳定性高、比传统的ADL模型具有更高的精度的数控机床热误差补偿高次多阶自回归分布滞后建模方法。
本发明为解决技术问题采用如下技术方案:
本发明数控机床热误差补偿高次多阶自回归分布滞后建模方法的特点是按如下步骤进行:
步骤1:定义高次多阶自回归分布滞后模型的表达式如式(1):
式(1)中xj,t-k为参加建模的第j个温度传感器第t-k次温度测量值;j=1,2,3L u,u为建模使用的温度传感器的个数;t=1,2,3L l,l为数控机床热误差测量次数;k为相对于t滞后k次,K=1,2,3L n;yt为数控机床热误差第t次测量值;yt-i是数控机床热误差第t-i次测量值,表示相对于t滞后i次的数控机床热误差测量值,i=1,2,3L m;m和n分别为yt和xj的最大滞后期,取值范围分别为1、2、3或4,具体取值根据第4步赤池信息准则判断;α0,αi,1,L αi,w和βj,k,1,βj,k,2Lβj,k,w均为待求系数;w为xj和yt的最高次数,w取值为2;
对于式(1),记:Km,n=[α0],
Cm=[α1,1,α1,2,L,α1,w,L,αm,1,αm,2,L,αm,w]T,
Dn=[β1,0,1,β1,0,2,L,β1,0,w,L,βu,n,1,βu,n,2,L,βu,n,w]T
将式(1)简代表达为式(2):
yt=Km,n+Am×Cm+Bn×Dn (2)
步骤2:通过最小二乘法算计待求系数Km,n、Cm、Dn,以确定式(2):
取数控机床热误差滞后期m和数控机床实测温度滞后期n分别为1、2、3和4,由式(2)根据最小二乘法计算得出:K1,1、C1、D1,K1,2、C1、D2,L,K4,4、C4、D4;
步骤3:求滞后期分别为m和n时yt估计值y′m,n,t:
将步骤2中得到的Am、Bn、Km,n、Cm、Dn代入公式y′m,n,t=Km,n+Am×Cm+Bn×Dn;计算数控机床热误差估计值序列y′m,n,l,y′m,n,l-1,L y′m,n,1,求出数控机床热误差yt滞后期m和数控机床实测温度xj滞后期n取值分别为1,、2、3和4时数控机床热误差估计值序列:
y′1,1,l,y′1,1,l-1,Ly′1,1,1
y′1,2,l,y′1,2,l-1,Ly′1,2,1
N
y′4,4,l,y′4,4,l-1,Ly′4,4,1;
步骤4:通过确定最佳滞后阶数m和n建立最终模型:
由数控机床热误差实测序列yl,yl-1,Ly1和由步骤3得到的数控机床热误差估计值序列y′m,n,l,y′m,n,l-1,Ly′m,n,1计算残差平方和 将残差平方和RSSm,n代入赤池信息量准则 计算出数控机床热误差yt滞后期m和机床温度xj滞后期n取值分别为1,、2、3和4时的赤池信息量AIC1,1、AIC1,2LAIC4,4;通过赤池信息量判断最佳滞后阶数,当赤池信息量AICm,n取值为最小时,对应的m和n即为数控机床热误差的高次多阶自回归分布滞后模型最佳滞后期;取得最佳滞后期时的m和n对应的Am、Bn、Km,n、Cm、Dn代入式(2),即为数控机床热误差补偿高次多阶自回归分布滞后模型。
与已有技术相比,本发明有益效果体现在:
1、因为本发明模型为高次自回归分布滞后模型,综合了传统多元回归模型和自回归分布滞后模型各自的优点,在应用中具有比传统多元回归模型和自回归分布模型更高的精度,所以本发明模型可以实现高精度的数控机床热误差补偿。
2、因为本发明数控机床热误差高次多阶自回归分布模型为高次,而传统的自回归分布滞后模型只为1次,所以本发明方法极大地丰富了数控机床热误差建模技术,为数控机床热误差建模提供了一种新的建模方法。
3、本发明建模技术使用简便、稳定性高、可靠性强,可通过软件编程实现在线数控机床X轴方向上和Z轴方向上的数控机床热误差补偿。
4、本发明提出了一种新型高次多阶自回归分布滞后模型,完善并丰富了传统的自回归分布滞后模型(ADL模型),为数学建模技术提供了一种新的参考方案。
具体实施方式
实施例中数控机床热误差yt(以X轴向热误差为例)和数控机床实测温度xj,t(j取值为1、2、3)记录如下:
表1数控机床热误差实测值和温度实测值
步骤1:定义高次多阶自回归分布滞后模型的表达式如式(1):
式(1)中xj,t-k为参加建模的第j个温度传感器第t-k次温度测量值;j=1,2,3Lu,u为建模使用的温度传感器的个数;t=1,2,3Ll,l为数控机床热误差测量次数;k为相对于t滞后k次,K=1,2,3Ln;yt为数控机床热误差第t次测量值;yt-i是数控机床热误差第t-i次测量值,表示相对于t滞后i次的数控机床热误差测量值,i=1,2,3Lm;m和n分别为yt和xj的最大滞后期,取值范围分别为1、2、3或4,具体取值根据第4步赤池信息准则判断;α0,αi,1,Lαi,w,βj,k,1,βj,k,2Lβj,k,w为待求系数;w为xj和yt的最高次数,w取值为2;
对于式(1),记:Km,n=[α0],
Cm=[α1,1,α1,2,L,α1,w,L,αm,1,αm,2,L,αm,w]T,
Dn=[β1,0,1,β1,0,2,L,β1,0,w,L,βu,n,1,βu,n,2,L,βu,n,w]T
则式(1)简代表达为如下式(2):
yt=Km,n+Am×Cm+Bn×Dn (2)
步骤2:通过最小二乘法算计待求系数Km,n、Cm、Dn,以确定式(2):
取数控机床热误差滞后期m和数控机床实测温度滞后期n分别为1、2、3和4,由式(2)根据最小二乘法计算得出:K1,1=[-1.1562],C1=[-1.6937,L,5.3559],D1=[0.0209,L,0.3030],L,D4=[0.0086,L,-1.0891]。
步骤3:求滞后期分别为m和n时yt估计值y′m,n,t:
将步骤2中得到的Am、Bn、Km,n、Cm、Dn代入公式y′m,n,t=Km,n+Am×Cm+Bn×Dn;计算数控机床热误差估计值序列y′m,n,l,y′m,n,l-1,L y′m,n,1,求出数控机床热误差yt滞后期m和数控机床实测温度xj滞后期n取值分别为1,、2、3和4时数控机床热误差估计值序列:
S′1,1=[0.9607,L,9.7104],L,S′4,4=[0.2400,L,9.6413];
步骤4:通过确定最佳滞后阶数m和n建立最终模型:
由数控机床热误差实测序列yl,yl-1,Ly1和由步骤3得到的数控机床热误差估计值序列y′m,n,l,y′m,n,l-1,L y′m,n,1计算残差平方和 将残差平方和RSSm,n代入赤池信息量准则 计算出数控机床热误差yt滞后期m和机床温度xj滞后期n取值分别为1、2、3和4时的赤池信息量AIC1,1=-1.7993,L,AIC4,4=-2.7316;通过赤池信息量判断最佳滞后阶数,当m和n取值为3时AICm,n取最小值,所以可以确定本实例中数控机床热误差与数控机床各实测温度之间的高次自回归分布滞后模型为:
yt=K3,3+A3×C3+B3×D3
式中:
K3,3=[35.4036],
C3=[0.0450,L,0.3533]T,
D3=[-0.7287,L,-4.7592]T。
该模型与其它算法拟合函数曲线如图1所示,图1中,1为试验测量获得的原始数据,2为多元线性回归模型拟合曲线,3为多元二次回归模型拟合曲线,4为自回归分布滞后模型拟合曲线,5为本发明二次自回归分布滞后模型拟合曲线。由图1可见自回归分布滞后模型精度较高,二次自回归分布滞后模型数据几乎与原始数据重合,精度最高。各模型计算的标准差如表2所示,其中,二次自回归分布滞后模型精度最高,比传统的ADL精度提高近一倍,比多元线性回归模型精度提高一个数量级。
表2不同模型标准差比较