发明名称 一种激光自准直测角零位基准误差测量方法
摘要 本发明属于测量技术,涉及对激光自准直测角技术方法的改进。本发明的实施基于:一个带有反射镜的三自由度旋转台,激光测角传感器,光电接收器件,数据采集卡和计算机。根据建立的激光自准直测角系统基准零位的畸变模型,按照测量步骤完成系统零位基准误差的测量。本发明方法对提高激光自准直测角系统的测量精度有着十分重要的作用,具有工程实用价值。
申请公布号 CN1862222A 申请公布日期 2006.11.15
申请号 CN200510069432.9 申请日期 2005.05.10
申请人 北京航空航天大学 发明人 张广军;尚鸿雁;魏振忠
分类号 G01B11/26(2006.01);G01B21/22(2006.01) 主分类号 G01B11/26(2006.01)
代理机构 中国航空专利中心 代理人 梁瑞林
主权项 1、一种激光自准直测角系统零位基准误差测量方法,其特征在于,1.1、本测量方法的实施基于以下装置:一个精度在角秒级的三自由度旋转台,在旋转台的旋转中心放置着平面反射镜;激光测角传感器,即测角系统的光学部分;光电接收器件;数据采集卡和计算机;1.2、建立激光自准直测角系统基准零位的畸变模型,首先,建立xyz坐标系,右手系。1.2.1、当光电接收器件的安装位置垂直于光轴且平面镜绕光轴旋转一周时,折射光线绕光轴旋转形成的正锥体的方程为:          [z+M]<sup>2</sup>=N(x<sup>2</sup>+y<sup>2</sup>)                 (1)式中:<maths num="001"><![CDATA[ <math><mrow><mi>M</mi><mo>=</mo><mfrac><mrow><mi>f</mi><mo>&CenterDot;</mo><mi>L</mi><mo>&CenterDot;</mo><mi>tan</mi><mi>&alpha;</mi></mrow><mrow><mo>|</mo><mrow><mo>(</mo><mi>f</mi><mo>-</mo><mi>L</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&alpha;</mi><mo>|</mo></mrow></mfrac><mo>+</mo><mi>f</mi><mo>,</mo><mi>N</mi><mo>=</mo><mfrac><msup><mi>f</mi><mn>2</mn></msup><mrow><mo>|</mo><mrow><mo>(</mo><mi>f</mi><mo>-</mo><mi>L</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&alpha;</mi><msup><mo>|</mo><mn>2</mn></msup></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>f为光学系统的焦距,L为测量距离,α为反射光线与光轴之间的夹角;1.2.2、当光电接收器件的安装位置发生倾斜时,光电接收器件倾斜平面方程为:       Ax+By+Cz+D=0                     (2)当光电接收器件倾斜平面位置确定时,式中A、B、C、D为确定的常数。1.2.3、联立公式(1)和公式(2),则像点在光电接收器件上的运动轨迹方程为:(A<sup>2</sup>-C<sup>2</sup>N)x<sup>2</sup>+2ABxy+(B<sup>2</sup>-C<sup>2</sup>N)y<sup>2</sup>+2A(D-CM)x+2B(D-CM)y+(D-CM)<sup>2</sup>=0   (3)式(3)表示的是一椭圆曲线,该椭圆曲线的中心像坐标(a,b,c)由下式给出:<maths num="002"><![CDATA[ <math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>a</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mi>A</mi><mrow><mo>(</mo><mi>D</mi><mo>-</mo><mi>CM</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>A</mi><mn>2</mn></msup><mo>+</mo><msup><mi>B</mi><mn>2</mn></msup><mo>-</mo><msup><mi>C</mi><mn>2</mn></msup><mi>N</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mi>b</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mi>B</mi><mrow><mo>(</mo><mi>D</mi><mo>-</mo><mi>CM</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>A</mi><mn>2</mn></msup><mo>+</mo><msup><mi>B</mi><mn>2</mn></msup><mo>-</mo><msup><mi>C</mi><mn>2</mn></msup><mi>N</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mi>c</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mi>Aa</mi><mo>-</mo><mi>Bb</mi><mo>-</mo><mi>D</mi></mrow><mi>C</mi></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>1.2.4、激光自准直测角系统实际零位基准(x<sub>o</sub>,y<sub>o</sub>,z<sub>o</sub>)是光轴与光电接收器件平面的交点,即<maths num="003"><![CDATA[ <math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>x</mi><mi>o</mi></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>o</mi></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>z</mi><mi>o</mi></msub><mo>=</mo><mo>-</mo><mfrac><mi>D</mi><mi>C</mi></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>1.2.5、系统基准零位的畸变模型d,<maths num="004"><![CDATA[ <math><mrow><mi>d</mi><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mi>a</mi><mo>-</mo><msub><mi>x</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>b</mi><mo>-</mo><msub><mi>y</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>c</mi><mo>-</mo><msub><mi>z</mi><mi>o</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths><maths num="005"><![CDATA[ <math><mrow><mo>=</mo><msqrt><mfrac><mrow><msup><mrow><mo>(</mo><mi>D</mi><mo>-</mo><mi>CM</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>&CenterDot;</mo><mrow><mo>(</mo><msup><mi>A</mi><mn>4</mn></msup><mo>+</mo><msup><mi>A</mi><mn>2</mn></msup><msup><mi>C</mi><mn>2</mn></msup><mo>+</mo><msup><mi>B</mi><mn>4</mn></msup><mo>+</mo><msup><mi>B</mi><mn>2</mn></msup><msup><mi>C</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><mi>A</mi><mn>2</mn></msup><msup><mi>B</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mrow><msup><mi>C</mi><mn>2</mn></msup><msup><mrow><mo>(</mo><msup><mi>A</mi><mn>2</mn></msup><mo>+</mo><msup><mi>B</mi><mn>2</mn></msup><mo>-</mo><msup><mi>C</mi><mn>2</mn></msup><mi>N</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>d为式(4)所代表的点与式(5)所代表的点两个点之间的距离;1.3、本发明方法的测量步骤如下:1.3.1、调整平面镜到初始零位状态,打开激光器,使得自准直测量系统的出射光束照射在平面镜上;1.3.2、转动平面镜,绕转台内框轴以角速度W旋转一周,1°/s<W<150°/s,通过数据采集与显示部分,每隔一定时间间隔t,记录光电接收器件探测到的像点坐标,设一共记录像点坐标p对,P的取值大于10,t=360/(W×P);1.3.3、对上述p对像点坐标进行圆拟合,求解得到圆心坐标(x<sub>1</sub>,y<sub>1</sub>)和圆半径R;1.3.4、根据已知的测量系统的光学系统的焦距f,在1.3.3的基础上,通过下列公式求解得到反射光线与光轴之间的夹角α:<maths num="006"><![CDATA[ <math><mrow><mi>&alpha;</mi><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mfrac><mi>R</mi><mi>f</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>1.3.5、对1.3.2中获得的p对像点坐标进行椭圆拟合,求解得到椭圆的圆心坐标(x<sub>2</sub>,y<sub>2</sub>)和下列方程的各个系数:     a<sub>1</sub>x<sup>2</sup>+2a<sub>2</sub>xy+a<sub>3</sub>y<sup>2</sup>+2a<sub>4</sub>x+2a<sub>5</sub>y+a<sub>6</sub>=0              (8)求取圆心坐标和椭圆圆心坐标之间的距离偏差,如果距离偏差小于光电接收器件的位置分辨率,就认为反射光线与光轴之间夹角的真实值近似等于夹角α,测量步骤跳转到第1.3.8步;否则进行第1.3.6步;1.3.6、在俯仰和偏航两个方向转动平面镜,使得反射光线与光轴之间的夹角减去α值,得到平面镜新的空间位置,令转动的俯仰角θ和偏航角满足下式:<img file="A2005100694320003C5.GIF" wi="1419" he="74" />1.3.7、在平面镜新的空间位置的基础上,绕转台内框轴转动一周,重复1.3.2~1.3.5的过程,直到反射光线与光轴之间夹角的真实值近似等于夹角α,根据式<maths num="007"><![CDATA[ <math><mrow><mi>N</mi><mo>=</mo><mfrac><msup><mi>f</mi><mn>2</mn></msup><msup><mrow><mo>|</mo><mrow><mo>(</mo><mi>f</mi><mo>-</mo><mi>L</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&alpha;</mi><mo>|</mo></mrow><mn>2</mn></msup></mfrac></mrow></math>]]></maths>求解得到N的值;1.3.8、根据1.3.5的拟和结果和夹角α的值,结合公式(3),得到光电接收器件平面方程系数间的关系:<maths num="008"><![CDATA[ <math><mrow><msup><mi>A</mi><mn>2</mn></msup><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mn>4</mn></msub><msub><mi>a</mi><mn>2</mn></msub></mrow><msub><mi>a</mi><mn>5</mn></msub></mfrac><mo>,</mo><msup><mi>B</mi><mn>2</mn></msup><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mn>5</mn></msub><msub><mi>a</mi><mn>2</mn></msub></mrow><msub><mi>a</mi><mn>4</mn></msub></mfrac><mo>,</mo><msup><mrow><mo>(</mo><mi>D</mi><mo>-</mo><mi>CM</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><msub><mi>a</mi><mn>6</mn></msub><mo>,</mo><msup><mi>C</mi><mn>2</mn></msup><mo>=</mo><mfrac><mrow><msup><mi>A</mi><mn>2</mn></msup><mo>-</mo><msub><mi>a</mi><mn>1</mn></msub></mrow><mi>N</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>1.3.9、把公式(10)的结果带入到公式(6),计算得到激光自准直测角系统基准零位的测量误差。
地址 100083北京市海淀区学院路37号