摘 要 常规数控系统的伺服控制中，当速度矢量变化较大时会产生大的轨迹跟踪误差. 本文应用M序列信号、最小二乘法和ARX模型结构辨识出整个闭环系统的传递函数. 在轨迹控制器的前面加上带预见控制的前馈补偿器，以使数控轨迹控制器在不同的频率段具有相同的增益， 消除由延迟带来的相位差，从而可减小轨迹跟踪误差，实验表明其性能较为理想.
关键词 动态特性；辨识；传递函数；前馈；预见补偿
中图资料分类号 TP 273
FEEDFORWARD TRACK-PREVIEWING COMPENSATOR baseD
ON THE IDENTIFICATION OF NC MILL DYNAMICS
Li Weiguang Chen Tongjian Xu Jianhua Yao Xifan
(Dept. of Mechatronic Eng., South China Univ. of Tech., Guangzhou 510640)
Abstract A great tracking error will occur when the velocity vectors sharply charged in a general servo system of a NC machine tool. The transfer function of closed-loop of the NC mill is identified by using the M serial signal, least square and ARX model. Considering the identification,a feedforward preview compensator is designed to keep the same magnitude at different frequencies,.and to eliminate the phase delay. The experimental results is showthat the tracking performance is improved by the feedforward preview compensator.
Key words dynamics; identification; transfer function; feedforward; preview compensation
数控铣床是机械行业加工异型轮廓零件的一种重要设备，其控制器的主要作用是根据零件的数控加工程序发出运动指令，跟踪给定的轨迹. 当今数控加工的潮流已向着高速、高精度方向发展，以取得高的生产效率和加工精度. 采用常规控制算法的带反馈控制的轨迹控制器已不能适应这种要求，这种控制器在较低的轨迹矢量速度下能达到一定的控制精度， 随着控制速度矢量的变化率加大，延迟现象明显加大，轨迹误差随之增大. 在反馈控制器的基础上加上轨迹前馈误差补偿器的目的是为了消除延迟，使控制系统的闭环传递函数保持为定值. 即使得控制系统在不同的频率段保持相同的增益，从而达到更高的轨迹跟踪性能［1～4］. 数控铣床的一个特点是所跟踪的轨迹预先已知，因此可在前馈控制器中加进预见控制环节，以进一步提高性能. 前馈补偿器的补偿效果取决于对整个控制系统(含被控对象)传递函数的获取，当未确定被控对象的动态特性时， 前馈补偿器中缺乏用于校正误差的补偿项参数，因此，即使采用预见控制也不能提高轨迹控制器的控制性能［4］. 本文采用伪随机M序列信号对一台数控铣床的伺服系统进行动态特性辨识.在得到数控铣床闭环传递函数的基础上，针对其零-极点的分布，设计并计算出前馈预见补偿器的形式和参数，通过实验验证了此前馈预见补偿器的有效性.
1 数控铣床伺服系统及实验装置
所研究的对象是一台经技改后配用美国著名的AUTOCON公司于1993年生产的DYNAPATH Delta20型数控系统的三坐标铣床. 铣床数控系统，伺服系统及为进行辨识实验而附加的信号发生、采集装置如图1所示，为简化仅画出X轴伺服部分.
图1 铣床数控系统， 伺服系统及信号发生、 采集装置
Fig.1 Structure on servo system of NC mill and identification device
辨识过程中，由微机发出M序列二值伪随机信号，该信号经PCL812板的D/A口以合适的电平送至数控系统的速度控制端，叠加在系统正常工作的信号上，在测速电机的输出端产生的响应信号经A/D转换后与M信号一起存放在微机数据文件中.
2 M序列信号与相关分析法辨识
2.1 M序列信号及功率谱
辨识过程中采用最大长度双电平伪随机M序列信号［5］，M序列信号统计特性较好，生成方法简便. 本研究中，由微机预先产生一个周期的M序列信号，然后循环输出. 经初步计算及实验后确定M序列信号周期长度为Np＝27－1＝127，基本电平时间(钟周期)大于4ms时可取得较理想的辨识效果. 图2中给出了钟周期＝11ms时M序列信号的功率谱及数控铣床在此信号激励下的输出功率谱. 此M序列信号近似于白噪声，它能激励出铣床在工作通频带范围内的动态行为.
图2 M序列信号功率谱及数控铣床的输出功率谱
Fig.2 Power spectrum of M serial signal and output power spectrum of NC mill
2.2 相关分析－最小二乘法辨识系统模型
相关分析法辨识的理论基础是维纳－何浦(Wiener_Hopt)方程［5］：
(1)
式中，Ruu(λ-τ)为输入量的自相关函数；Ruy(λ)为输入/输出量的互相关函数；为系统冲激响应的最佳估计.
设有一SISO线性离散时间系统可表示为：
y(t)＝-a1y(t-1)-a2y(t-2)-…-any(t-n)+b1u(t-1)+b2u(t-2)+…+bnu(t-n)+e(t)=φTtθ+e(t) (2)
式中：{u(t)}和{y(t)}分别为实际测量的输入和输出序列，e(t)为误差项.
φTt＝［-y(t-1),…,-y(t-n), u(t-1),…,u(t-n)
θT=［a1,…an, b1,…,bn］
向量矩阵形式为
Y=φTθ+E
残差为：
应用最小二乘法的辨识就是确定系统参数向量θ的估计值， 使残差平方和为最小.
将相关分析法和最小二乘法结合起来可同时辨识出非参数模型及参数模型.
图3是辨识X轴时输入输出量的相关函数曲线(利用R9211B仪器获得).从图中可看到：M序列输入信号的自相关性、伺服运动输出速度的自相关性以及输入、输出的互相关性都很显著. 这说明，用于辨识的输入 输出信号可信度较高.
图3 X轴的相关函数曲线
Fig.3 Correlation functions on X axis
3 数控铣床轨迹控制模型的辨识
当数控铣床的机械性能较好时可略去其非线性因素. 因而， 参考位置输入量和实际位置输出量间采样序列u(t)和y(t)可表示为常系数线性差分方程：
y(t)+a1y(t-1)+…+any(t-n)=b0u(t-d)+
b1u(t-1-d)+…+bmu(t-m-d) (3)
式中d为延迟量.
对式(3)进行Z变换后可得到：
(4)
式(4)中的即为系统的传递函数.
辨识实验中系统的I/O信号波形如图4所示（M序列信号的钟周期为50ms，采样周期为7.81ms)，采用ARX模型结构，对采样的前半部分信号进行辨识计算，得到X轴的传递函数，其中延迟步数为4，主导极点为7，零点数等于6，见表1, 计算出差分方程式(3)的系数见表2.
图4 采样数据波形
Fig.4 Waves of signal measured
表1 X轴主导零-极点数据
Table 1 Major zero-pole on X axis
用辨识出的模型画出仿真曲线，和后半部分采样信号进行比较，其拟合较理想，见图5.
图5 采样数据曲线与模型仿真曲线
Fig.5 Measured curve and simulated curve
表2 模型系数
Table 2 Coefficient of model
4 前馈预见补偿器设计
4.1 前馈预见补偿器的结构
在不带前馈补偿的数控铣床反馈控制器中，由数控程序发出的位置指令yd(t)就是参考位置输入量u(t)，因此当t＞0时，伺服系统动态实际位置y(t)不等于给定的位置yd(t)(伺服系统中，A(Z-1)不等于B(Z-1))，采用前馈补偿器的目的就是要使实际位置尽可能紧密跟踪给定的位置yd(t).为此要根据不同情况设计前馈补偿器的形式. 图6中表示的是前馈补偿器与闭环反馈系统间的关系.
数控铣床轨迹控制的一个显著特点是未来的轨迹点可预先知道，因此在本前馈补偿器的设计中也考虑利用此未来值. 即采用预见控制方式，以进一步减小轨迹控制误差.
在数控机床系统中，零点在单位圆外的情况(称为非最小相位系统)极为普遍［4］.由辨识结果可知：本研究所采用的数控铣床也是一非最小相位系统.
将图6中的零点多项式B(Z-1)分解为：
B(Z-1)=Bin(Z-1).Bout(Z-1).B1(Z-1) (5)
图6 前馈补偿器与闭环反馈系统间的位置关系
Fig.6 Location about feedforward preview compensator
and closed-feedback system
式中Bin(Z-1)为单位圆内的零点,B1(Z-1)为单位圆上的零点,Bout(Z-1)为单位圆外的零点.
这样，铣床闭环系统的传递函数可表示为：
(6)
因此， 前馈补偿器中分别对相应的几部分进行补偿：
{F(Z-1), F(Z)}=F1(Z).F2(Z-1).F3(Z) (7)
其中, F1(Z)＝Zd, 用于消除系统延迟d步而产生相位差的影响; F2(Z-1)=A(Z-1)/Bin(Z-1)，用于删除落在单位圆内和圆上的极-零点；F3(Z)用于消除系统中单位圆外零点的影响；对于F3(Z)，如取其逆系统1/Bout(Z-1)，就会导致系统产生不稳定的极点，因此采用与其等价的复数形式进行补偿［2］.
设系统共有k个单位圆外的零点，则
其复数的指数形式展开为幂级数后：
(8)
由于数控轨迹控制中未来s步的目标值可预先算出，则代入式(8)后就实现了与 Bout(Z-1)的逆系统的等价. 由于这些零点在单位圆外，所以级数收敛，因此取式(8)的前有限项就可获得近似结果. 考虑到取有限项后引起的增益改变，将式(8)除以系数 {1-(rleiθl)-nsZns}得到Bout(Z-1)的近似逆系统为：
(9)
4.2 前馈预见补偿器的参数计算
对于辨识获得的数控铣床X轴伺服系统的传递函数， 其前馈补偿器参数为：
所以
(10)
式(10)中,r1=r2=1.0442，θ1= -θ2＝2.04726rad，h1=h2=0.5
单位圆外一对共轭复数零点的前馈预见补偿器为：Fout(Z)=Fout(Z)1+Fout(Z)2
其中：
确定了预见步数ns后，Fout(Z)1和Fout(Z)2中的各虚数项可一一消去，Fout(Z)就成为简单的实系数补偿器.
5 实验结果及分析
根据以上辨识得出的数控铣床模型，加上前馈预见补偿器(预见步数取3)，构造出一新的数控伺服控制器，进行实验. 补偿前后X轴的Bode图如图7所示.
图7 补偿前后X轴的Bode图
Fig.7 Bode diagrams compensated and un-compensated on X axis
从图7中可以看到, 补偿前的幅频增益特性和相频特性随着频率增大而明显下降, 频率在4 Hz时增益约下降18dB, 相位滞后约180°；频率在10 Hz时增益约下降31dB, 相位滞后约270°；频率在40Hz时增益约下降39dB, 相位滞后约700°. 加前馈预见补偿后幅频增益在25 Hz内时基本保持定值,相频特性在频率达100Hz时，只有约100°的相位差.
经比较后可知，采用前馈预见补偿器后铣床伺服控制器的频率特性比不加补偿器时有明显改善，尤其在低频段效果更显著.
6 结 论
采用合适参数的M序列伪随机信号和最小二乘法辨识出了数控铣床的动态传递函数，相关分析及谱图表明辨识过程中的I/O信号可信度较高. 根据数控轨迹可知的特点，分别对闭环系统中延迟步数、极点以及零点进行相应的补偿，构成了一个前馈预见补偿器，实验结果表明补偿后的频率特性有很大提高.
