发明名称 一种基于先验误差分解定权的工业测量拟合方法
摘要 本发明涉及一种基于先验误差分解定权的工业测量拟合方法,包括以下步骤:1)输入测量点集数据,按照各测点等权,初步拟合空间几何体,然后执行步骤2);2)确定各测站‑测点连线以及几何体的空间姿态,然后执行步骤3);3)针对不同几何体形状,根据误差分解和定权算法对各测点的先验误差进行分解,并以此为依据为各测点加权,然后执行步骤4);4)依据新定的权,对测量点集进行拟合,获得加权的空间几何体拟合参数。现有技术相比,本发明具有削弱仪器测量误差对拟合结果的影响,提高工件拟合精度等优点。
申请公布号 CN104034291B 申请公布日期 2017.02.22
申请号 CN201410253961.3 申请日期 2014.06.09
申请人 同济大学 发明人 潘国荣;郭巍
分类号 G01B21/00(2006.01)I 主分类号 G01B21/00(2006.01)I
代理机构 上海科盛知识产权代理有限公司 31225 代理人 赵志远
主权项 一种基于先验误差分解定权的工业测量拟合方法,其特征在于,包括以下步骤:1)输入测量点集数据,按照各测点等权,初步拟合空间几何体,然后执行步骤2);2)确定各测站‑测点连线以及几何体的空间姿态,然后执行步骤3);3)针对不同几何体形状,根据误差分解和定权算法对各测点的先验误差进行分解,并以此为依据为各测点加权,然后执行步骤4);4)依据新定的权,对测量点集进行拟合,获得加权的空间几何体拟合参数;所述的根据误差分解和定权算法对各测点的先验误差进行分解具体为:11)误差分量定义:111)径向误差:圆、球形工件在测点处指向圆心或球心方向的误差分向量,记为E<sub>R</sub>;112)切向误差:与圆、球形工件在测点处切线或切面方向的误差分量,记为E<sub>T</sub>;113)垂向误差:在测点处垂直于直线、平面或圆面的误差分向量,记为E<sub>V</sub>;114)平向误差:在测点处平行于直线、平面或圆面的误差分向量,记为E<sub>H</sub>;12)令测量仪器的测角误差为ζ,测距常数误差为a,测距系数误差为k,记仪器中心位置为P<sub>0</sub>(X<sub>0</sub>,Y<sub>0</sub>,Z<sub>0</sub>),记某测点测量结果为P<sub>i</sub>(X<sub>i</sub>,Y<sub>i</sub>,Z<sub>i</sub>),i∈N,不同的几何形状分解方式具体如下:121)空间直线与平面测点的误差分解,考虑测站位置和工件姿态,分析某测距、测角精度下线、面上不同位置测点的垂向误差和平向误差,并依此作为定权依据,得其函数关系为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>E</mi><msub><mi>V</mi><mi>i</mi></msub></msub><mo>=</mo><msub><mi>&epsiv;</mi><msub><mi>&alpha;</mi><msub><mi>V</mi><mi>i</mi></msub></msub></msub><mo>+</mo><msub><mi>&epsiv;</mi><msub><mi>d</mi><msub><mi>V</mi><mi>i</mi></msub></msub></msub><mo>=</mo><msub><mi>&epsiv;</mi><msub><mi>&alpha;</mi><mi>i</mi></msub></msub><mo>&CenterDot;</mo><msub><mi>cos&alpha;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><msub><mi>d</mi><mn>1</mn></msub></msub><msub><mi>sin&alpha;</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0001141122590000011.GIF" wi="638" he="61" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>E</mi><msub><mi>H</mi><mi>i</mi></msub></msub><mo>=</mo><msub><mi>&epsiv;</mi><msub><mi>&alpha;</mi><msub><mi>H</mi><mi>i</mi></msub></msub></msub><mo>+</mo><msub><mi>&epsiv;</mi><msub><mi>d</mi><msub><mi>H</mi><mi>i</mi></msub></msub></msub><mo>=</mo><msub><mi>&epsiv;</mi><msub><mi>&alpha;</mi><mi>i</mi></msub></msub><mo>&CenterDot;</mo><msub><mi>sin&alpha;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>&epsiv;</mi><msub><mi>d</mi><mi>i</mi></msub></msub><msub><mi>cos&alpha;</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141122590000012.GIF" wi="1388" he="70" /></maths>其中,<img file="FDA0001141122590000013.GIF" wi="341" he="130" />Dist<sub>i</sub>为测点P<sub>i</sub>到测站P<sub>0</sub>的距离,a为测距常数误差,k为测距系数误差;对于直线和平面,只有垂向误差分量影响其拟合形状和姿态,在此设单位权中误差为1,则其测点拟合定权为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>E</mi><msub><mi>V</mi><mi>i</mi></msub></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141122590000021.GIF" wi="1053" he="152" /></maths>122)空间圆与球面测点的误差分解空间圆拟合可看为空间平面拟合和平面圆拟合的组成,将误差先依据圆平面分解成垂向误差分量和平向误差分量,方法与式(1)相同,然后再按式(3)将平向误差分解为径向误差分量和切向误差分量<maths num="0004"><math><![CDATA[<mrow><msub><mi>E</mi><msub><mi>R</mi><mi>i</mi></msub></msub><mo>=</mo><msub><mi>E</mi><msub><mi>H</mi><mi>i</mi></msub></msub><mo>&CenterDot;</mo><msub><mi>sin&beta;</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0001141122590000022.GIF" wi="286" he="63" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>E</mi><msub><mi>T</mi><mi>i</mi></msub></msub><mo>=</mo><msub><mi>E</mi><msub><mi>H</mi><mi>i</mi></msub></msub><mo>&CenterDot;</mo><msub><mi>cos&beta;</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141122590000023.GIF" wi="1044" he="54" /></maths>其中β<sub>i</sub>为平向误差分量与圆上测点切线的夹角;对于球面,在测点处切平面上的垂向误差分量和平向误差分量即为对该球面的径向误差分量和切向误差分量,因此求球面测点的切面是球面测点误差分解的第1步,随后的分解方法与式(1)相同;对空间圆测点的权值计算方法为:<maths num="0006"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>E</mi><msub><mi>V</mi><mi>i</mi></msub></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mo>(</mo><msqrt><mrow><msubsup><mi>E</mi><msub><mi>T</mi><mi>i</mi></msub><mn>2</mn></msubsup><mo>+</mo><msup><mi>R</mi><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mi>R</mi><mo>)</mo></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>E</mi><msub><mi>R</mi><mi>i</mi></msub></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001141122590000024.GIF" wi="1814" he="174" /></maths>对球面测点的权值计算方法为:<maths num="0007"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mo>(</mo><msqrt><mrow><msubsup><mi>E</mi><msub><mi>T</mi><mi>i</mi></msub><mn>2</mn></msubsup><mo>+</mo><msup><mi>R</mi><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mi>R</mi><mo>)</mo></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><msub><mi>E</mi><msub><mi>R</mi><mi>i</mi></msub></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0001141122590000025.GIF" wi="1806" he="167" /></maths>
地址 200092 上海市杨浦区四平路1239号