Correction method for angle measurement error of eccentric encoder with double read-heads installed non-diametrically opposite

非对径安装双读数头圆光栅偏心测角误差修正 Correction method for angle measurement error of eccentric encoder with double read-heads installed non-diametrically opposite 1 2 first-author 王笑一（1976-），男，河南洛阳人，博士，副教授，1999年于洛阳工学院获得学士学位，2005年于河南科技大学获得硕士学位，2016年于北京工业大学获得博士学位，主要研究方向为精密测量、齿轮测量和数控技术。E-mail： wxy2.0@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850931&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850937&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850933&type= 22.01333618 32.00402451 王笑一 （1976-），男，河南洛阳人，博士，副教授，1999年于洛阳工学院获得学士学位，2005年于河南科技大学获得硕士学位，2016年于北京工业大学获得博士学位，主要研究方向为精密测量、齿轮测量和数控技术。E-mail： wxy2.0@163.com wxy2.0@163.com 1 1 corresp E-mail： ly_lxq@163.com E-mail： ly_lxq@163.com ly_lxq@163.com 1 2 河南科技大学 河南省机械设计及传动系统重点实验室， 河南 洛阳 471003 河南科技大学 河南省机械设计及传动系统重点实验室， 河南 洛阳 471003 Henan Key Laboratory of Mechanical Design and Transmission System of Henan University of Science and Technology， Luoyang 471003， China Henan Key Laboratory of Mechanical Design and Transmission System of Henan University of Science and Technology， Luoyang 471003， China 贵阳新天光电科技有限公司， 贵州 贵阳 550018 贵阳新天光电科技有限公司， 贵州 贵阳 550018 Guiyang Xintian Photoelectric Technology Co.， Ltd， Guiyang 550018， China Guiyang Xintian Photoelectric Technology Co.， Ltd， Guiyang 550018， China 圆光栅安装偏心引起的测角误差是影响圆光栅角度测量精度的关键因素，而对圆光栅安装偏心引起的测角误差进行修正是提高圆光栅测角精度的有效方法。本文分析了圆光栅安装偏心与测角误差之间的关系，建立了基于双读数头的圆光栅偏心测角误差理论模型，提出了一种基于非对径安装双读数头的圆光栅偏心测角误差修正方法。仿真及实测试验表明，在双读数头非对径安装时，新方法对读数头安装误差和圆光栅偏心误差造成的测角误差的修正效果明显优于双读数头均值法。试验中，当两读数头对径安装误差约为4°时，新方法修正后的测角误差为均值法修正后测角误差的二分之一以下，其中仿真试验中由均值法修正后测角误差为1.785″， 新方法修正后的测角误差为0.720″。 The angle measurement error caused by the eccentric installation of the circular grating is the key factor affecting the angle measurement accuracy of the circular grating. It is an effective method for correcting the angle measurement error caused by the eccentric installation of an encoder. In this paper， the relation between the eccentricity of the circular grating installation and the angle measurement error is analyzed， the development of a theoretical model of the eccentric angle measurement error of the circular grating based on a dual read head is presented， and a new method for correcting the eccentric angle measurement error of the circular grating based on a non-diametrically installed dual read head is proposed. The simulation and experimental results showed that the new method is better than the conventional mean method in correcting the angle measurement error caused by the installation error of the dual read head and by the eccentric error of the circular grating. In the test， when the installation error of the diametrical alignment of the dual read head was approximately 4°， the angular measurement error corrected by the new method was less than that corrected by the mean method. In the simulation， the corrected angular measurement error of the mean method was 1.785 inches， whereas that of the new method was 0.720 inch. 圆光栅 双读数头 测角误差 误差修正 circular grating double reading head angle measurement error error correction 1 引　言 1 实时精确地获取仪器、机床、机器人等设备的旋转运动角位移信息是实现精密测量和加工的关键之一。圆光栅是最常用的高精度角位移传感器，被广泛应用于多种精密仪器设备中 ［ 1 1 - 7 7 1-7 ］ 。在测量精度要求较高时，由圆光栅安装偏心引起的测角误差是不能忽略的，即使是很小的安装误差都会造成很大影响 ［ 8 8 ］ 。因此，对圆光栅安装偏心引起的测角误差进行修正是非常有必要的。 目前，解决由于圆光栅安装偏心引起的测角误差的方法主要有两种：一是用更高精度的仪器检测出测角误差并通过软件方法进行修正补偿，二是采用多读数头的方法修正测角误差等 ［ 9 9 - 10 10 9-10 ］ 。 折文集等 ［ 11 11 ］ 使用多面体及自准直仪作为基准对圆光栅测角误差进行标定，并在测角误差标定的基础上建立测角误差补偿模型，对测角误差进行补偿。张文颖等 ［ 12 12 ］ 提出了一种自准直仪的在位校准方法，基于圆周封闭原则和傅里叶级数的性质，建立了测量值与理想值之间的函数关系并搭建校准测角系统。王福全等 ［ 13 13 ］ 利用激光干涉仪与校准回转装置对转台进行标定，提出了一种基于稀疏分解的角分度误差补偿方法，对转台的角分度误差进行补偿。以上方法主要是采用更高精度的仪器作为基准检测出测角误差后并通过软件方法进行修正补偿，能够有效的减小测角误差；但是当测角精度要求较高时，利用此方法对测角误差进行修正补偿成本会大大的提高，且补偿效果也不理想。 艾晨光等 ［ 14 14 ］ 研究了圆光栅安装偏心引起的测角误差，并采用双读数头对测角误差进行修正。郭阳宽等 ［ 15 15 ］ 研究了运动偏心对圆光栅测量影响的分析方法，推导了由于偏心引起的测角误差补偿公式。唐松等 ［ 16 16 ］ 分析了安装偏心对圆光栅测角精度的影响，提出了一种两个读数头对径安装读数头的补偿方法。王文等 ［ 17 17 ］ 分析了圆光栅分度盘的安装偏心对角度测量的影响，提出了一种基于模拟退火算法角度传感器偏心参数辨识方法并将其用于测量关节圆光栅12个偏心参数的辨识和修正。王园等 ［ 18 18 ］ 建立了圆光栅偏心误差模型，提出了一种双读数头平均误差补偿方法，对读数头误差进行修正。以上方法主要是采用多读数头的方法对圆光栅安装偏心引起的测角误差进行修正补偿，能够有效对测角误差进行补偿修正，但是对读数头安装精度要求较高，当安装角误差 > 1°时，无法有效消除安装偏心对测角精度的不利影响。 本文提出了一种新的误差补偿方法，可使用非精确对径安装的双读数头对圆光栅的偏心误差进行修正，当读数头安装误差>1°时（试验中最大为9°），相对于传统的均值法，使用本文方法能够大幅提高测角精度。本文方法一方面可以减小对径安装误差对测量结果的不利影响；另一方面突破了读数头必须对径安装的限制，提高了仪器设备在结构设计时的自由度。 2 误差分析与均值法数学模型 1 2.1　圆光栅测角误差分析 2 圆光栅测角误差主要可分为系统误差和随机误差两部分。其中，圆光栅测角装置的系统误差来源主要包括圆光栅盘的刻划误差、圆光栅盘的安装偏心、圆光栅盘端面和回转轴线的垂直度、轴系轴承的径向跳动中的系统误差成分等。一般来说，误差补偿主要是针对其中的系统误差成分进行的。焦杨等研究过类似的数学模型，其中没有考虑随机误差的影响 ［ 19 19 ］ 。但是在精密测量条件下，随机误差和系统误差的幅值经常是接近的；在系统误差得到正确的修正之后，随机误差的幅值有可能超过系统误差的幅值。为使仿真试验的结果更加接近于实测试验，在圆光栅测角误差模型中加入了随机误差的成分。综合系统误差和随机误差成分，在多读数头圆光栅测角系统中，第 i 个读数头的测角误差如 式（1） 式（1） 所示 ［ 19 19 - 25 25 19-25 ］ ： <math display="block" id="M1"><mtable><mtr><mtd><msub><mrow><mi>θ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><mi>θ</mi><mo>=</mo><mi>κ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>+</mo><mi>ε</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>+</mo><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mo stretchy="false">[</mo><mi>x</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">]</mo><mo>/</mo><mi>R</mi><mo>+</mo><mi>ψ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>+</mo><mi>δ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850944&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850940&type= 73.15200043 10.24466610 ， 其中： <math id="M2"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851160&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851157&type= 1.77800000 2.79399991 为无误差的理论旋转角度， <math id="M3"><msub><mrow><mi>θ</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850957&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850955&type= 2.53999996 3.72533321 为读数头 i 的读数值， <math id="M4"><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850965&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850961&type= 2.37066650 3.72533321 为读数头 i 所在的安装位置，即读数头 i 的安装位置与起始位置的夹角（在双读数头系统中，无安装误差时，读数头1对应 <math id="M5"><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850970&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850968&type= 10.66800022 3.72533321 ，读数头2对应 <math id="M6"><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">180</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850976&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850973&type= 14.39333248 3.72533321 ）； <math id="M7"><msub><mrow><mi>κ</mi></mrow><mrow><mi mathvariant="normal">i</mi></mrow></msub><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850981&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850979&type= 16.25600052 4.48733330 为圆光栅安装偏心引起的读数头i的测角误差； <math id="M8"><mi>ε</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850986&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850984&type= 6.85799980 3.55599999 为圆光栅的刻划误差，不同位置的读数头测得的刻划误差存在一个相位偏移； <math id="M9"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850993&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850989&type= 2.70933342 2.79399991 为光栅盘的半径， <math id="M10"><mi>x</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851001&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16850997&type= 7.45066643 3.55599999 ， <math id="M11"><mi>y</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851009&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851004&type= 7.28133297 3.80999994 分别是回转轴系沿 X ， Y 方向的径向运动误差， <math id="M12"><mo stretchy="false">[</mo><mi>x</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><mi>y</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">]</mo><mo>/</mo><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851017&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851013&type= 43.60333252 4.82600021 表示径向运动误差对读数头测角的影响； <math id="M13"><mi>ψ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851036&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851031&type= 7.61999989 3.72533321 为读数头内部的不一致性引起的测角误差，该项误差频率较高、幅值较小，但规律不明确； <math id="M14"><mi>δ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851028&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851026&type= 7.36600018 3.55599999 为由于其他不确定因素引起的随机误差。 测角误差是衡量圆光栅测量精度的标准，在实际应用中刻划误差与随机误差虽无法避免，但是相较于圆光栅安装偏心引起的测角误差其幅值较小，通常情况下圆光栅安装偏心及轴系径向运动是引起单读数头圆光栅测角误差的关键因素。读数头不一致性 <math id="M15"><mi>ψ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851036&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851031&type= 7.61999989 3.72533321 属于高频误差，在进行整周大范围角度测量或者测角间隔为信号周期整数倍时通常可以忽略 ［ 19 19 ］ 。 2.2　均值法数学模型 2 对于双读数头圆光栅测角系统，采用均值法是常规的误差修正方法。均值法误差修正的理论模型为： <math display="block" id="M16"><mtable><mtr><mtd><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover><mo>=</mo><mo stretchy="false">(</mo><mi>κ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><mi>ε</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>+</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mo stretchy="false">[</mo><mi>x</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo></mtd></mtr><mtr><mtd><mi>y</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">]</mo><mo>/</mo><mi>R</mi><mo>+</mo><mi>ψ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>δ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>κ</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo></mtd></mtr><mtr><mtd><mi>ε</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>+</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mo stretchy="false">[</mo><mi>x</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><mi>y</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">]</mo><mo>/</mo><mi>R</mi><mo>+</mo></mtd></mtr><mtr><mtd><mi>ψ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>δ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo>/</mo><mn mathvariant="normal">2</mn><mo>+</mo><mi>θ</mi></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851043&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851040&type= 72.38999939 22.35199928 . （2） 均值法计算简单，但对两读数头的安装精度要求较高。当读数头没有安装误差时，式中 <math id="M17"><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">0</mn><mo>°</mo><mo>,</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mn mathvariant="normal">180</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851052&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851047&type= 26.83933067 4.14866638 。此时，公式（2）可变形为： <math display="block" id="M18"><mtable columnalign="left"><mtr><mtd><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover><mo>=</mo><mo stretchy="false">(</mo><mi>ε</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>+</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>ψ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>δ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><mi>ε</mi><mo stretchy="false">(</mo><mi>θ</mi><mo>+</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>ψ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mi>δ</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo>/</mo><mn mathvariant="normal">2</mn><mo>+</mo><mi>θ</mi></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851059&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851056&type= 51.73133087 10.83733368 . 式（3） 式（3） 中，由于圆光栅安装偏心及轴系径向跳动引起的测角误差都被消除掉了。但是，当读数头存在安装误差的时候，均值法是不能完全消除这两项误差的影响的。 通过以上对圆光栅测角系统中的误差来源和影响机理的分析可知，双读数头均值法应用在有较大的读数头安装误差的环境下存在一定的局限性。为此，提出了基于非对径安装双读数头的圆光栅测角误差修正方法，进行了数字仿真对比试验，并搭建试验台开展了实测试验。 3 新方法数学模型 1 如 图1 图1 所示， <math id="M19"><msub><mrow><mi>O</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851140&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851136&type= 4.06400013 3.72533321 为轴系回转中心， <math id="M20"><mi>O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851074&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851070&type= 2.96333337 2.87866688 为圆光栅盘的几何中心， <math id="M21"><mi>C</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851147&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851143&type= 2.62466669 2.79399991 ， <math id="M22"><mi>D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851154&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851151&type= 2.79399991 2.79399991 为两读数头的安装位置。 <math id="M23"><mi>O</mi><msub><mrow><mi>O</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851095&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851091&type= 6.94266701 3.72533321 为圆光栅盘的偏心距，用 <math id="M24"><mi>e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851102&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851099&type= 1.52400005 2.79399991 表示； <math id="M25"><mi mathvariant="normal">∠</mi><msub><mrow><mi>O</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>O</mi><mi>D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851110&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851106&type= 13.54666710 4.06400013 为圆光栅盘的偏心角，用 <math id="M26"><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851117&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851115&type= 2.53999996 2.79399991 表示。 <math id="M27"><mi mathvariant="normal">∠</mi><msub><mrow><mi>O</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>O</mi><mi>C</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851124&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851121&type= 13.29266548 4.06400013 用 <math id="M28"><mi>η</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851131&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851128&type= 1.86266661 3.64066648 表示。为方便表达， 图1 图1 中假设轴和圆光栅盘固定，而读数头绕回转中心 <math id="M29"><msub><mrow><mi>O</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851140&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851136&type= 4.06400013 3.72533321 逆时针旋转，从 <math id="M30"><mi>C</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851147&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851143&type= 2.62466669 2.79399991 ， <math id="M31"><mi>D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851154&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851151&type= 2.79399991 2.79399991 点分别旋转 <math id="M32"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851160&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851157&type= 1.77800000 2.79399991 角到达 <math id="M33"><mi>A</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851166&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851163&type= 2.87866688 2.79399991 ， <math id="M34"><mi>B</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851171&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851169&type= 2.53999996 2.79399991 点。 <math id="M35"><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851276&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851274&type= 7.78933382 4.14866638 为两个读数头在这个旋转过程中测得的旋转角度。两读数头的安装角 <math id="M36"><mi mathvariant="normal">∠</mi><mi mathvariant="normal">C</mi><msub><mrow><mi>O</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi mathvariant="normal">D</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851182&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851179&type= 13.20800018 4.06400013 用 <math id="M37"><mi>φ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851188&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851186&type= 2.37066650 3.64066648 表示。当双读数头是严格的对径安装，即没有读数头安装位置误差时， <math id="M38"><mi>φ</mi><mo>=</mo><mn mathvariant="normal">180</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851194&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851191&type= 13.63133430 3.64066648 。 由三角形外角等于不相邻两内角和，可得读数头1的测角误差： <math display="block" id="M39"><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><mi>θ</mi><mo>=</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>ω</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851204&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851202&type= 74.59133911 8.89000034 . 同理可得读数头2的测角误差： <math display="block" id="M40"><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><mi>θ</mi><mo>=</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>η</mi></mrow></mfenced><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>η</mi><mo>+</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851208&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851206&type= 70.78133392 8.89000034 . 从 图1 图1 可知： <math display="block" id="M41"><mtable><mtr><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>η</mi><mo stretchy="false">[</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>α</mi><mo>+</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>+</mo><mi>γ</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mo>+</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>η</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>+</mo><mi>γ</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mn mathvariant="normal">0</mn></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851212&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851210&type= 67.22533417 10.24466610 . 令 <math id="M42"><mi>m</mi><mo>=</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>α</mi><mo>+</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>+</mo><mi>γ</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851216&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851214&type= 64.09266663 4.65666676 ， <math id="M43"><mi>n</mi><mo>=</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>+</mo><mi>γ</mi><mo>+</mo><mi>φ</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851218&type= 50.79999924 4.65666676 则： <math display="block" id="M44"><mi>m</mi><mo>⋅</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>η</mi><mo>+</mo><mi>n</mi><mo>⋅</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>η</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851225&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851223&type= 34.71333313 3.72533321 . 根据辅助角公式可得： <math display="block" id="M45"><mroot><mrow><msup><mrow><mi>m</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot><mo>⋅</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>η</mi><mo>+</mo><mi>σ</mi><mo stretchy="false">)</mo><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851229&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851226&type= 41.31733322 4.82600021 ， <math display="block" id="M46"><mtable><mtr><mtd><mi>η</mi><mo>=</mo><mn mathvariant="normal">180</mn><mo>-</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfrac><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfenced open="[" close="]" separators="|"><mrow><mi>ω</mi><mo>+</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo></mrow></mfenced><mo>+</mo><mi>φ</mi></mrow></mfenced></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>ω</mi><mfenced open="[" close="]" separators="|"><mrow><mo>+</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo></mrow></mfenced><mo>+</mo><mi>φ</mi></mrow></mfenced><mo>+</mo><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac></mrow></mfrac></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851234&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851231&type= 65.95532990 25.82333374 ， 由 式（4） 式（4） ， 式（5） 式（5） ， 式（9） 式（9） 可得两读数头测量值与偏心参数之间的关系式： <math display="block" id="M47"><mtable><mtr><mtd><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>ω</mi></mrow></mfenced><mo>+</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>η</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>η</mi></mrow></mfenced><mo>=</mo><mn mathvariant="normal">0</mn></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851241&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851238&type= 73.82933807 19.64266586 ， 式（10） 式（10） 中只包含两个读数头的读数和圆光栅的安装偏心参数。其中两个读数头读数 <math id="M48"><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851276&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851274&type= 7.78933382 4.14866638 可由数据采集系统进行采集。将轴系旋转一周即可获得一系列的 <math id="M49"><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851276&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851274&type= 7.78933382 4.14866638 ，根据 式（10） 式（10） 使用最小二乘法即可拟合出圆光栅安装偏心参数（最小二乘法为常用算法，在此不再赘述），从而实现圆光栅安装偏心参数的标定。 由 式（4） 式（4） ， 式（5） 式（5） ， 式（9） 式（9） 可得测角误差补偿模型： <math display="block" id="M50"><mtable columnalign="left"><mtr><mtd><mover accent="true"><mi>θ</mi><mo>̂</mo></mover><mo>=</mo><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo stretchy="false">(</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>ω</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>ω</mi></mrow></mfenced><mo>-</mo></mtd></mtr><mtr><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>η</mi><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced><mo>+</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>e</mi></mrow><mrow><mi>R</mi></mrow></mfrac><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mi>η</mi><mo stretchy="false">)</mo><mo>-</mo><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851258&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851255&type= 72.98266602 19.64266586 ， 其中： <math id="M51"><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851264&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851262&type= 1.77800000 3.64066648 表示经过新方法补偿后的转动角度测量值，补偿后的测角误差可表示为 <math id="M52"><mover accent="true"><mi>θ</mi><mo>̂</mo></mover><mo>-</mo><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851270&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851267&type= 8.46666718 3.64066648 。实际应用中，首先让被测轴系带动圆光栅旋转一周，以获取圆光栅安装偏心参数和两个读数头的安装位置误差参数，然后将两读数头的读数 <math id="M53"><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851276&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851274&type= 7.78933382 4.14866638 代入 式（11） 式（11） ，即可得到补偿后的测角结果，其中已经消除了圆光栅安装偏心和读数头对径安装误差的不利影响。 4 试　验 1 为了对比新方法和均值法在基于非对径安装双读数头的圆光栅偏心测角误差修正中的应用效果，开展了数字仿真试验和实测试验。 相关研究中通常使用多棱镜和自准直仪来验证误差修正效果，但这种方案只能在少数几个特定角度位置进行验证，不能对任意角位置的修正效果进行验证。如果要实现对任意角位置的修正效果进行验证，需要有能够连续测量的角位移基准仪器，且仪器的测量精度要足够高。本研究的实物试验中使用的圆光栅系统的分辨率为0.158″，精度（按产品规格书）为±4″，校准难度大，实验室尚无能够满足连续角度测量精度要求的校准仪器。 为此，提出两个方案解决试验难题。第一是使用数字仿真试验，在仿真条件下对新方法和均值法的误差修正效果进行对比；第二是在实测试验中提出了同时安装4个读数头，基于试验系统自身测量数据进行“自标定”、“自校准”的技术方案，在不使用更高等级仪器的条件下就可以对新方法和均值法的误差修正效果进行比对。 仿真试验中，首先不引入圆光栅刻划误差及随机误差，用新方法和均值法两种方法对圆光栅安装偏心和读数头位置误差引起的测角误差进行修正并对比修正效果；其次，为了更接近实际，引入了圆光栅刻划误差及随机误差，再用两种方法对测角误差进行修正并对比修正效果。 实测试验中，共进行两次数据采集。第1次采集数据时，4个读数头安装在相隔90°均布的4个位置上，尽可能减小读数头的安装位置误差。第2次采集数据时，保持其中3个读数头的位置不变，改变另1个读数头的位置，使其存在明显的安装误差。4个读数头组成两组双读数头测量系统，其中一组读数头在两次数据采集中都没有对径安装位置误差，其数据用于验证两次采集的试验数据的重复精度；另一组读数头在第1次数据采集中没有对径安装位置误差，在第2次数据采集中有较大的对径安装位置误差，其前后两次数据用于对比新方法和均值法这两种误差修正方法的效果。 4.1　无刻划误差及随机误差仿真试验 2 设偏心距 <math id="M54"><mi>e</mi><mo>=</mo><mn mathvariant="normal">10</mn><mtext> </mtext><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851352&type= 15.66333294 3.64066648 ，圆光栅盘半径 <math id="M55"><mi>R</mi><mo>=</mo><mn mathvariant="normal">26</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851360&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851357&type= 17.61066628 2.79399991 ，偏心角为 <math id="M56"><mi>ω</mi><mo>=</mo><mn mathvariant="normal">60</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851366&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851363&type= 11.93799973 2.79399991 ，两读数头安装误差为2°，将以上预设参数代入 式（10） 式（10） 中，用最小二乘法对仿真数据进行拟合，得到偏心距 <math id="M57"><mi>e</mi><mo>=</mo><mn mathvariant="normal">10.070</mn><mtext> </mtext><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851299&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851296&type= 22.26733208 3.64066648 ，偏心角 <math id="M58"><mi>ω</mi><mo>=</mo><mn mathvariant="normal">59.465</mn><mtext> </mtext><mn mathvariant="normal">867</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851305&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851302&type= 24.80733490 2.79399991 ，两读数头安装误差为 <math id="M59"><mn mathvariant="normal">2.000</mn><mtext> </mtext><mn mathvariant="normal">030</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851312&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851309&type= 15.83266735 2.79399991 ，与输入参数基本一致，证明了偏心参数拟合算法的正确性。拟合结果与输入参数的差别可能是由于运算中圆周采样点数仅为360点，数据密度不高造成的。将以上求出的偏心参数代入 式（4） 式（4） ， 式（5） 式（5） ， 式（11） 式（11） 中，可得到两个读数头的测角误差、均值法补偿后的测角误差和新方法补偿后的测角误差，如 图2 图2 所示。从 图2 图2 可以看出，利用新方法能够有效消除由偏心及读数头安装误差引起的测角误差。 4.2　有刻划误差及随机误差的仿真试验 2 为了更接近实际，在仿真试验中引入了圆光栅刻划误差及随机误差。 图3 图3 所示为仿真中使用的圆光栅刻划误差，图中标出了读数头1，2分别用到的刻线误差区段。 图4 图4 ~ 图5 图5 为叠加在读数头1，2测量结果上的随机误差。 设偏心距 <math id="M60"><mi>e</mi><mo>=</mo><mn mathvariant="normal">10</mn><mtext> </mtext><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851352&type= 15.66333294 3.64066648 ，圆光栅盘半径 <math id="M61"><mi>R</mi><mo>=</mo><mn mathvariant="normal">26</mn><mtext> </mtext><mi mathvariant="normal">m</mi><mi mathvariant="normal">m</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851360&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851357&type= 17.61066628 2.79399991 ，偏心角 <math id="M62"><mi>ω</mi><mo>=</mo><mn mathvariant="normal">60</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851366&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851363&type= 11.93799973 2.79399991 ，读数头2的安装误差从0°到9°变化，间隔1°为一组，共做10组实验。在这个试验条件下，当读数头安装误差为9°以上时，补偿后的测角误差将 > 1.5″，认为不能满足精密测角的需要，故不再继续试验。 试验结果显示，两读数头安装误差逐渐变大的情况下，补偿后测角误差曲线形状基本一致，但幅值逐渐增大。 图6 图6 所示为读数头安装误差为0°（即无安装误差）时的仿真实验结果。 图7 图7 所示为读数头安装误差为4°时的仿真实验结果。 图8 图8 是两读数头安装误差为1°～9°时，补偿后测角误差相对于安装误差为0°时补偿后测角误差的变化量。 从以上实验中可知，当安装误差在9°以内时，利用新方法补偿后的测角误差均≤1.5″；而双读数头均值法在安装误差≥2°时补偿后测角误差已>1.5″，在一些精密测量条件下已不能满足要求。其中，当安装误差为4°时，利用新方法补偿后的测角误差为0.720″；而双读数头均值法补偿后测角误差为1.785″。 仿真试验结果表明，当读数头安装误差较大时，无论是否引入圆光栅刻划误差及随机误差，新方法的补偿效果都比均值法要好。 4.3　实测试验 2 为了对比均值法和新方法的误差修正效果，通常要采用更高精度的仪器作为基准。但本试验中使用的圆光栅测角分辨率为0.158″，精度（按产品规格书）约为±4″，实验室尚无能够满足要求的校准仪器。为此，搭建了有4个读数头的光栅测角试验装置，读数头的安装方式如 图9 图9 所示。其中读数头1~3为一对，读数头2~4为一对，分别组成双读数头测角系统。圆光栅选用雷尼绍RESM20，直径为52 mm；读数头选用雷尼绍VIONiC20系列，分辨率为20 nm。 试验中共进行两次数据采集。第1次采集数据时，4个读数头按图示实线位置安装在相隔90°均布的4个位置上，尽可能减小读数头的安装位置误差。第2次采集数据时，保持读数头1，3，4位置不变，将读数头2移动到图中虚线位置处，试验中测得其安装误差（即偏离实线位置的角位移）名义值为4.154 138°（该名义值在试验中由圆光栅实测得到）。 在每次数据采集后，分别利用新方法及双读数头均值法对两对读数头1~3，2~4的测角值进行误差补偿。补偿后测量结果的记号规定如下：新方法的补偿结果记为 <math id="M63"><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851407&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851403&type= 1.77800000 3.64066648 ，均值法的补偿结果记为 <math id="M64"><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851412&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851409&type= 1.77800000 3.38666677 ；读数头1~3的结果下标为1，读数头2~4的结果下标为2；第1次数据采集的下标为F（意思是First），第2次数据采集的下标为S（意思是Second）。按以上约定，试验数据可记为 <math id="M65"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851418&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851415&type= 4.23333359 4.57200003 ， <math id="M66"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851423&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851421&type= 4.23333359 4.57200003 ， <math id="M67"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851605&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851603&type= 4.14866638 4.57200003 ， <math id="M68"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851434&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851432&type= 4.14866638 4.57200003 ， <math id="M69"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851438&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851436&type= 4.57200003 4.40266657 ， <math id="M70"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851442&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851440&type= 4.57200003 4.40266657 ， <math id="M71"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851445&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851443&type= 4.48733330 4.40266657 ，和 <math id="M72"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851449&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851447&type= 4.48733330 4.40266657 ，详细的符号说明见 表1 表1 。 方法 符号 符号说明 备注 4 均值法 <math id="M73"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851516&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851514&type= 4.14866638 3.80999994 第1次数据采集中，读数头1~3的数据用均值法的补偿结果 在以下图中作为横坐标数据； 试验中 <math id="M74"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851456&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851455&type= 3.80999994 3.64066648 与 <math id="M75"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851460&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851458&type= 3.89466691 3.64066648 差值曲线的幅值很小，说明两次试验的重复性较好 <math id="M76"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851492&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851490&type= 4.14866638 3.80999994 第1次数据采集中，读数头2~4的数据用均值法的补偿结果 由于4个读数头安装位置不同， <math id="M77"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851468&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851466&type= 3.89466691 3.64066648 与 <math id="M78"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851516&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851514&type= 4.14866638 3.80999994 的差值曲线可认为是刻划误差曲线的函数 <math id="M79"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851476&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851474&type= 4.06400013 3.80999994 第2次数据采集中，读数头1~3的数据用均值法的补偿结果 参见 <math id="M80"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851480&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851478&type= 3.89466691 3.64066648 备注 <math id="M81"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851488&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851486&type= 4.14866638 3.80999994 第2次数据采集中，读数头2~4的数据用均值法的补偿结果 试验中 <math id="M82"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851488&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851486&type= 4.14866638 3.80999994 与 <math id="M83"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851492&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851490&type= 4.14866638 3.80999994 曲线差异大，说明均值法不适用于读数头安装误差大的场合 4 本文方法 <math id="M84"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851497&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851494&type= 3.89466691 3.89466691 第1次数据采集中，读数头1~3的数据用新方法的补偿结果 试验中 <math id="M85"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851501&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851499&type= 3.80999994 3.89466691 与 <math id="M86"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851504&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851502&type= 3.89466691 3.89466691 差值的曲线幅值很小，说明两次试验的重复性较好 <math id="M87"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851507&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851505&type= 3.89466691 3.89466691 第1次数据采集中，读数头2~4的数据用新方法的补偿结果 由于4个读数头安装位置不同， <math id="M88"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851532&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851531&type= 3.89466691 3.89466691 与 <math id="M89"><msub><mrow><mover><mrow><mi>θ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851516&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851514&type= 4.14866638 3.80999994 的差值曲线可认为是刻划误差的函数 <math id="M90"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851520&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851518&type= 3.80999994 3.89466691 第2次数据采集中，读数头1~3的数据用新方法的补偿结果 参见 <math id="M91"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851523&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851522&type= 3.89466691 3.89466691 备注 <math id="M92"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851536&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851534&type= 3.80999994 3.89466691 第2次数据采集中，读数头2~4的数据用新方法的补偿结果 试验中 <math id="M93"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851532&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851531&type= 3.89466691 3.89466691 与 <math id="M94"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn>2</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851536&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851534&type= 3.80999994 3.89466691 曲线差异小，说明本文方法可用于读数头安装误差较大的场合 为了比对新方法和均值法的误差修正效果，按照 表1 表1 中定义的数据符号及其说明，可绘制出曲线。如 图10 图10 ~ 图11 图11 所示分别为 <math id="M95"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851605&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851603&type= 4.14866638 4.57200003 与 <math id="M96"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851547&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851545&type= 4.23333359 4.57200003 差值的曲线和 <math id="M97"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851551&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851549&type= 4.14866638 4.23333359 与 <math id="M98"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851601&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851599&type= 4.23333359 4.23333359 差值的曲线。由于两次测量时读数头1~3位置没有改变，所以当两次试验的重复性很好，即未引入其他误差时，两条差值曲线的幅值应当很小。从图中可以看出两次测量值曲线差值的幅值非常小（≤0.000 018°），说明两次试验的重复性较好。 图10 图10 中曲线的幅值大于 图11 图11 中曲线的幅值，这一点与预期不符，其原因是：为了和 图12 图12 保持一致， 图10 图10 中的曲线选用 <math id="M99"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851601&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851599&type= 4.23333359 4.23333359 作为横坐标，需要对 <math id="M100"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851605&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851603&type= 4.14866638 4.57200003 曲线进行插值才能得到 图10 图10 中曲线上对应点的数据，因此 图10 图10 中曲线的幅值显得大于 图11 图11 中曲线的幅值了。 如 图12 图12 所示为 <math id="M109"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851631&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851629&type= 13.20800018 4.23333359 和 <math id="M110"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851642&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851641&type= 13.20800018 4.57200003 的差值曲线。虽然读数头都没有安装误差，但是读数头1~3和2~4的安装位置存在90°相位差，这两条插值曲线的起伏主要来源于刻划误差的作用，可近似看作刻划误差的函数。 图12 图12 中， <math id="M111"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851639&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851637&type= 13.20800018 4.23333359 ， <math id="M112"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851642&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851641&type= 13.20800018 4.57200003 差值曲线的形状基本一致，由此可知当读数头2~4没有安装误差时，新方法和均值法具有几乎相同的误差修正效果。 如 图13 图13 所示分别为 <math id="M117"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851645&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851643&type= 13.20800018 4.23333359 ， <math id="M118"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851649&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851647&type= 13.20800018 4.57200003 ， <math id="M119"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851654&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851652&type= 13.03866673 4.57200003 ， <math id="M120"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851658&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851656&type= 13.03866673 4.23333359 的差值曲线。相比第1次测量，第2次测量时仅有读数头2的位置发生了约4.15°的变化，两对读数头1~3和2~4的安装位置仍存在约90°的相位差。如果均值法和本文方法都能完全消除测角误差，那么 图13 图13 中4条曲线的形状应当是基本一致的。从 图13 图13 中可以看出，有3条曲线，即 <math id="M121"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851661&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851660&type= 13.20800018 4.23333359 ， <math id="M122"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851665&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851663&type= 13.20800018 4.57200003 ， <math id="M123"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851669&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851667&type= 13.03866673 4.57200003 差值曲线基本重合，而 <math id="M124"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851674&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851672&type= 13.03866673 4.23333359 差值曲线与其他3条曲线存在较大差异。如 图14 图14 所示分别为差值曲线（ <math id="M125"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851678&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851676&type= 13.03866673 4.57200003 ） - （ <math id="M126"><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>̂</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851681&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851679&type= 13.20800018 4.57200003 ）和（ <math id="M127"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">S</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851685&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851683&type= 13.03866673 4.23333359 ） - （ <math id="M128"><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">F</mi></mrow></msub><mo>-</mo><msub><mrow><mover accent="true"><mi>θ</mi><mo>¯</mo></mover></mrow><mrow><mn mathvariant="normal">1</mn><mi mathvariant="normal">F</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851690&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16851687&type= 13.20800018 4.23333359 ），其中新方法曲线（红线）的变化幅度比均值法曲线（黑线）的变化幅度小45%。 图13 图13 和 图14 图14 说明当读数头存在安装误差时，新方法能够更有效地修正测角误差，效果明显优于均值法。 实测试验数据表明，对于安装误差较大的一对读数头，新方法前后两次修正后的数据一致性较好，而均值法前后两次修正后的数据一致性较差，说明当读数头安装误差较大时，新方法的补偿效果优于均值法。 5 结　论 1 为了提高安装双读数头的圆光栅测角精度，提出了基于非对径安装双读数头的圆光栅偏心测角误差修正方法。 首先通过分析圆光栅安装偏心与测角误差之间的关系，建立了两读数头读数值与偏心系数之间的关系式，并利用最小二乘法拟合出偏心系数，实现了圆光栅安装偏心参数的标定；然后建立了非对径安装双读数头偏心测角误差补偿模型，实现了对测角误差的修正。试验中，为了克服缺少高精度角度测量基准这个困难，开展了数字仿真试验，设计了基于试验系统自身条件的“自校准”方案，在不使用更高等级仪器的条件下就可以对新方法和均值法的误差修正效果进行比对。 仿真及实测试验表明新方法能够对非对径安装双读数头偏心测角误差进行有效地修正，且修正效果明显优于双读数头均值法。 新方法也存在一些不足：（1）应用新方法时，开机后首先需要圆光栅旋转一周以上，获得圆光栅安装的偏心参数和读数头的安装误差数据之后，才能使用进行精度较好的误差修正，而传统的均值法没有这个限制。（2）新方法应用于实践时，由于运算量较大，暂时只能用于实时性要求不高的场合，比如用于高精度分度、对准等仪器设备中；而不能用于对实时性要求较高的场合，如参与多轴联动加工的机床精密转台。 对于以上不足，可以采取一些措施进行弥补。在开机后圆光栅旋转一周之前，可以用均值法作为临时的误差修正结果，而在完成一周旋转之后再切换到新方法；对于运算量较大的问题，未来将尝试使用FPGA实现修正算法，借助FPGA的快速运算能力提高修正计算的实时性。 参考文献 XI-JUN C ， ZHEN-HUAN W ， QING-SHUANG Z . Angle measurement error and compensation for decentration rotation of circular gratings ［J］. Journal of Harbin Institute of Technology ， 2010 ， 17 （ 4 ）： 536 - 539 . Angle measurement error and compensation for decentration rotation of circular gratings 江炜 ， 周维虎 ， 劳达宝 ， 等 . 多光栅测角系统精密转台轴系设计 ［J］. 仪表技术与传感器 ， 2018 ，（ 9 ）： 24 - 28 . 多光栅测角系统精密转台轴系设计 JIANG W ， ZHOU W H ， LAO D B ， et al . Design of precision turntable shaft for multi-gratings angle measuring system ［J］. Instrument Technique and Sensor ， 2018 ，（ 9 ）： 24 - 28 . （in Chinese） Design of precision turntable shaft for multi-gratings angle measuring system 杨婧 ， 冯其波 . 数控机床空间几何误差测量研究进展 ［J］. 仪器仪表学报 ， 2017 ， 38 （ 8 ）： 1901 - 1911 . 数控机床空间几何误差测量研究进展 YANG J ， FENG Q B . Research progress on volumetric geometric error measurement of numerical control （NC） machine tools ［J］. China Journal of Scientific Instrument ， 2017 ， 38 （ 8 ）： 1901 - 1911 . （in Chinese） Research progress on volumetric geometric error measurement of numerical control （NC） machine tools 于连栋 ， 鲍文慧 ， 赵会宁 ， 等 . 新型圆光栅测角误差补偿方法及其应用 ［J］. 光学 精密工程 ， 2019 ， 27 （ 8 ）： 1719 - 1726 . 新型圆光栅测角误差补偿方法及其应用 YU L D ， BAO W H ， ZHAO H N ， et al . Application and novel angle measurement error compensation method of circular gratings ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 8 ）： 1719 - 1726 . （in Chinese） Application and novel angle measurement error compensation method of circular gratings 娄志峰 ， 郝秀朋 ， 刘力 ， 等 . 圆光栅配合自准直仪测量主轴径向运动误差 ［J］. 光学 精密工程 ， 2019 ， 27 （ 9 ）： 2053 - 2061 . 圆光栅配合自准直仪测量主轴径向运动误差 LOU ZH F ， HAO X P ， LIU L ， et al . Spindle radial motion error measurement using a circular grating and a autocollimator ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 9 ）： 2053 - 2061 . （in Chinese） Spindle radial motion error measurement using a circular grating and a autocollimator 石照耀 ， 尚岳 ， 王笑一 ， 等 . 长光栅传感器动态误差检测系统设计 ［J］. 北京工业大学学报 ， 2016 ， 42 （ 3 ）： 321 - 328 . 长光栅传感器动态误差检测系统设计 SHI ZH Y ， SHANG Y ， WANG X Y ， et al . Design of Dynamic Error Measurement System for Liner-encoders ［J］. Journal of Beijing University of Technology ， 2016 ， 42 （ 3 ）： 321 - 328 . （in Chinese） Design of Dynamic Error Measurement System for Liner-encoders 陈希军 ， 王振桓 ， 曾庆双 . Angle measurement error and compensation for decentration rotation of circular gratings ［J］. Journal of Harbin Institute of Technology ， 2010 ， 17 （ 4 ）： 536 - 539 . （in Chinese） Angle measurement error and compensation for decentration rotation of circular gratings CHEN X J ， WANG ZH H ， ZENG Q SH . Angle measurement error and compensation for decentration rotation of circular gratings ［J］. Journal of Harbin Institute of Technology ， 2010 ， 17 （ 4 ）： 536 - 539 . （in Chinese） Angle measurement error and compensation for decentration rotation of circular gratings 高贯斌 . 关节臂式坐标测量机自标定方法与误差补偿研究 ［D］. 杭州 ： 浙江大学 ， 2010 . 关节臂式坐标测量机自标定方法与误差补偿研究 GAO G B . Study on Self-calibration Method and Error Compensation of Articulated Arm Coordinate Measuring Machine ［D］. Hangzhou ： Zhejiang university ， 2010 . （in Chinese） Study on Self-calibration Method and Error Compensation of Articulated Arm Coordinate Measuring Machine 费斌 ， 徐文雄 ， 蒋庄德 ， 等 . 计量圆光栅动态精度的检测技术 ［J］. 仪器仪表学报 ， 1996 ， （ 5 ）： 78 - 82 . 计量圆光栅动态精度的检测技术 FEI B ， XU W ， JIANG ZH D ， et al . Measurement technology of dynamic precision of metering circular grating ［J］. China Journal of Scientific Instrument ， 1996 ， （ 5 ）： 78 - 82 . （in Chinese） Measurement technology of dynamic precision of metering circular grating 钱建强 ， 王东生 . 一种闭环控制圆光栅角度分度测量装置 ［J］. 仪器仪表学报 ， 2006 ， 27 （ S2 ）： 1515 - 1516 . 一种闭环控制圆光栅角度分度测量装置 QIAN J Q ， WANG D SH . Circular grating angle-division measuring equipment based on closed loop controlling system ［J］. China Journal of Scientific Instrument ， 2006 ， 27 （ S2 ）： 1515 - 1516 . （in Chinese） Circular grating angle-division measuring equipment based on closed loop controlling system 折文集 ， 高立民 ， 杨永清 ， 等 . 激光跟踪仪测角误差标定及补偿 ［J］. 仪器仪表学报 ， 2015 ， 36 （ Z1 ）： 52 - 56 . 激光跟踪仪测角误差标定及补偿 SHE W J ， GAO L M ， YANG Y Q ， et al . Calibration and compensation of angular error in laser tracker ［J］. Chinese Journal of Scientific Instrument . 2015 ， 36 （ Z1 ）： 52 - 56 . （in Chinese） Calibration and compensation of angular error in laser tracker 张文颖 ， 朱浩然 ， 李美萱 ， 等 . 基于自准直仪的测角传感器实时在位校准方法 ［J］. 中国激光 ， 2019 ， 46 （ 4 ）： 166 - 172 . 基于自准直仪的测角传感器实时在位校准方法 ZHANG W Y ， ZHU H R ， LI M X ， et al . Real-time In-Situ Calibration for Angle Measuring Sensor Based on Autocollimator ［J］. Chinese Journal of Lasers ， 2019 ， 46 （ 4 ）： 166 - 172 . （in Chinese） Real-time In-Situ Calibration for Angle Measuring Sensor Based on Autocollimator 王福全 ， 王珏 ， 谢志江 ， 等 . 精密转台角分度误差补偿 ［J］. 光学 精密工程 ， 2017 ， 25 （ 8 ）： 2165 - 2172 . 精密转台角分度误差补偿 WANG F Q ， WANG J ， XIE ZH J ， et al . Compensation of angular indexing error for precision turntable ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 8 ）： 2165 - 2172 . （in Chinese） Compensation of angular indexing error for precision turntable 艾晨光 ， 褚明 ， 孙汉旭 ， 等 . 基准圆光栅偏心检测及测角误差补偿 ［J］. 光学 精密工程 ， 2012 ， 20 （ 11 ）： 2479 - 2484 . 基准圆光栅偏心检测及测角误差补偿 AI CH G ， CHU M ， SUN H X ， et al . Eccentricity testing of benchmark circular grating and compensation of angular error ［J］. Opt. Precision Eng. ， 2012 ， 20 （ 11 ）： 2479 - 2484 . （in Chinese） Eccentricity testing of benchmark circular grating and compensation of angular error 郭阳宽 ， 李玉和 ， 李庆祥 ， 等 . 运动偏心对圆光栅副测量的影响及误差补偿 ［J］. 清华大学学报（自然科学版） ， 2005 ， （ 2 ）： 178 - 181 . 运动偏心对圆光栅副测量的影响及误差补偿 GUO Y K ， LI Y H ， LI Q X ， et al . Error compensation for eccentric motions of circular gratings ［J］. Journal of Tsinghua University（Science and Technology） ， 2005 ， （ 2 ）： 178 - 181 . （in Chinese） Error compensation for eccentric motions of circular gratings 唐松 ， 方明 ， 姚磊 . 安装偏心对圆光栅测角精度的影响及误差补偿技术研究 ［J］. 计量与测试技术 ， 2013 ， 40 （ 6 ）： 11 - 12 . 安装偏心对圆光栅测角精度的影响及误差补偿技术研究 TANG S ， FANG M ， YAO L . Influence of installation eccentricity on angular measurement accuracy of circular gratings and error compensation technology ［J］. Metrology & Measurement Technique ， 2013 ， 40 （ 6 ）： 11 - 12 . （in Chinese） Influence of installation eccentricity on angular measurement accuracy of circular gratings and error compensation technology 王文 ， 林铿 ， 高贯斌 ， 等 . 关节臂式坐标测量机角度传感器偏心参数辨识 ［J］. 光学 精密工程 ， 2010 ， 18 （ 1 ）： 135 - 141 . 关节臂式坐标测量机角度传感器偏心参数辨识 WANG W ， LIN K ， GAO G B ， et al . Eccentricity parameter identification of angle sensors for articulated arm CMMs ［J］. Opt. Precision Eng. ， 2010 ， 18 （ 1 ）： 135 - 141 . （in Chinese） Eccentricity parameter identification of angle sensors for articulated arm CMMs 王园 . 双读数头圆光栅偏心参数标定修正算法 ［J］. 工具技术 ， 2018 ， 52 （ 11 ）： 142 - 146 . 双读数头圆光栅偏心参数标定修正算法 WANG Y . Correction algorithm for double reading head eccentricity parameters calibration ［J］. Tool Eengineering ， 2018 ， 52 （ 11 ）： 142 - 146 . （in Chinese） Correction algorithm for double reading head eccentricity parameters calibration 焦扬 ， 黄明 ， 刘品宽 ， 等 . 双角度编码器超精密转台测角误差校准 ［J］. 光学 精密工程 ， 2019 ， 27 （ 10 ）： 2180 - 2191 . 双角度编码器超精密转台测角误差校准 JIAO Y ， HUANG M ， LIU P K ， et al . Calibration of angular errors of high-precision rotary table with two-angle encoders ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 10 ）： 2180 - 2191 . （in Chinese） Calibration of angular errors of high-precision rotary table with two-angle encoders RALF D G ， ANDREAS J . A shearing-based method for the simultaneous calibration of angle measuring devices ［J］. Measurement Science and Technology ， 2014 ， 25 （ 10 ）. A shearing-based method for the simultaneous calibration of angle measuring devices RALF D G ， ALFRED L ， MICHAEL K ， et al . Capabilities and limitations of the self-calibration of angle encoders ［J］. Measurement Science and Technology ， 2014 ， 25 （ 5 ）. Capabilities and limitations of the self-calibration of angle encoders YANG J ， YE D ， ZEGUANG D ， et al . Optimal-arrangement-based four-scanning-heads error separation technique for self-calibration of angle encoders ［J］. Measurement Science and Technology . 2018 ， 29 （ 8 ）. Optimal-arrangement-based four-scanning-heads error separation technique for self-calibration of angle encoders 王显军 . 光电轴角编码器细分信号误差及精度分析 ［J］. 光学 精密工程 ， 2012 ， 20 （ 2 ）： 379 - 386 . 光电轴角编码器细分信号误差及精度分析 WANG X J . Errors and precision analysis of subdivision signals for photoelectric angle encoders ［J］. Opt. Precision Eng. ， 2012 ， 20 （ 2 ）： 379 - 386 . （in Chinese） Errors and precision analysis of subdivision signals for photoelectric angle encoders 乔栋 ， 续志军 ， 吴宏圣 ， 等 . 绝对式光栅尺细分误差补偿方法 ［J］. 光学学报 ， 2015 ， 35 （ 1 ）： 220 - 225 . 绝对式光栅尺细分误差补偿方法 QIAO D ， XU ZH J ， WU H SH ， et al . A Method for Compensating Interpolation Error of Absolute Linear Encoder ［J］. Acta Optica Sinica ， 2015 ， 35 （ 1 ）： 220 - 225 . （in Chinese） A Method for Compensating Interpolation Error of Absolute Linear Encoder 冯英翘 ， 万秋华 . 小型光电编码器细分误差校正方法 ［J］. 仪器仪表学报 ， 2013 ， 34 （ 6 ）： 1374 - 1379 . 小型光电编码器细分误差校正方法 FENG Y Q ， WAN Q H . Interpolation error calibration method of small photoelectric encoders ［J］. Chinese Journal of Scientific Instrument ， 2013 ， 34 （ 6 ）： 1374 - 1379 . （in Chinese） Interpolation error calibration method of small photoelectric encoders 10.37188/OPE.20212905.1103 TH712 A 1004-924X（2021）05-1103-12 王笑一,王永军,雷贤卿等.非对径安装双读数头圆光栅偏心测角误差修正[J].光学精密工程,2021,29(05):1103-1114. 王笑一,王永军,雷贤卿等.非对径安装双读数头圆光栅偏心测角误差修正[J].光学精密工程,2021,29(05):1103-1114. WANG Xiao-yi,WANG Yong-jun,LEI Xian-qing,et al.Correction method for angle measurement error of eccentric encoder with double read-heads installed non-diametrically opposite[J].Optics and Precision Engineering,2021,29(05):1103-1114. WANG Xiao-yi,WANG Yong-jun,LEI Xian-qing,et al.Correction method for angle measurement error of eccentric encoder with double read-heads installed non-diametrically opposite[J].Optics and Precision Engineering,2021,29(05):1103-1114. 国家自然科学基金资助项目（No.51775172）；中国博士后科学基金资助项目（No.188521） 2020-10-14 2021-01-08 2021-05-15 2021-06-24 光学精密工程 5 29