基于最小二乘法的平面抛物线轮廓度的误差评定

本资料包含pdf文件1个，下载需要1积分

在工程应用当中，广泛涉及到平面线轮廓如渐开线、椭圆、抛物线、摆线和萁舌线的轮廓度测量、识别和误差评 ”。如在模式识别和计算机视觉中，图形 (图像)数据的模型拟合和匹配是- 项基本工作。CAD中在理论和实际应用中经常遇到二次曲线的拟合问题l 。在径流式叶轮的设计和制造过程中，抛物线型叶片以其良好的气动性能、强度及工艺性被普遍采用口。在汽车制造业的深孔加工中抛物线型麻花钻以其排屑流畅、刚性好、使用寿命长代替了传统的麻花 。在数控系统加工零件过程中，抛物线轮廓的插补拟合算法的研究对提高拥有抛物线轮廓零件的精度有着很高的实用价值 。因此研究抛物线轮廓度误差评价方法对保证抛物线型叶片及其他抛物线型零件的加工质量和精度有着重要的意义。

关于抛物线的轮廓度误差的评定，国家标准尚未给出明确的定义和特定的评定算法求解抛物线轮廓度误差，近年来国内外学者专门对抛物线的研究较少，比较有代表性的成果分为几何距离拟合和代数距离拟合；几何距离拟合由H.Spith在1996提出的用正交距离最小二乘法拟合抛物线 和Sung Joon Ahn在H.Spath的基础上对椭圆、抛物线等二次曲线用正交距离的最小二乘法进行拟合 ”；在代数距离拟合中，刘海香对二次曲线的最小二乘法拟合进行了阐述和比较，指出在二次曲线拟合中抛物型的曲线较难拟合，抛物型曲线的曲率比较大，拟合出来的误差曲线理想程度较差[2。这些拟合方法对于抛物线误差的评定都有-定的效果和借鉴作用。

本文通过平面任意位置的抛物线方程，抛物线法线方程和抛物线本身的性质，利用最小二乘原理和约束条件实现对平面任意位置抛物线轮廓度误差的评定。

1 最dx-乘抛物线拟合平面任意位置抛物线的表达式F(x， )用平面- 般二次曲线方程表示较为合适。设平面抛物线方程为：F(x，J，)Ax BxyC Dx置 F (1)设 ( ， )(f1，23.忉 为抛物线上的N个测量点，根据最小二乘原理拟合的目标函数为：ⅣF(A，B，C，D，E，F)∑(ax Bxy 观 ，F) (2)il为使得F为最小，使： 篆o IaC aD (3) ----I- 1aA 8B aE 8F 、由此可得矩阵方程：， H W w H M∑ Zx， M ∑ Yt ∑ ∑ ∑H N N N∑ ∑ ∑ ∑ ∑ ∑∑ YJ ∑xlYt Ey， ∑ ∑ ∑Ex， ∑ Yt ∑ ∑ ∑ M ∑∑ ∑ Zy， ∑ Zy， Zy，I1 lml I-I I-I l t-I∑ ∑ ∑ ∑ Zy， N：。 (4)收稿日翔：2013-01-26基盒项目：国家自然科学基金 (50875076)；河南戍础与前沿技术研究计划项目 (122300413209，122300410114)作者简介：王海洋 (1987-)，男，硕士研究生，研究方向为精密测试技术。

第35卷 第5期 2013-05(上) 23 解方程 (4)会得出0解，为了得 出A、B、C、D、E和F的值，需要对方程加限制条件。根据- 般二次曲线的平面解析几何知识9二元二次方程为抛物线时， 令HAC，则H都是在坐标轴的平移和旋转变换之下的变量 ，在这里我们令ACI带入到方程 (4)中，便可以解出A、B、C、D、E和F的值，方程(1)中的六个系数与抛物线的位置参数(0， ，Y )及形状参数p存在以下关系 ：IAsin 0lB--sinOcos0 c： c。sz1D- sin sinOcosO-pcosO (5)lE sin0cos0- cos 0-pcosOF(xsinO-Ycos ) 2p( ∞s sin )[每/2 CD/2 BE/4 CE/2 BD/2 (- A CC ] l仰- 肚-fI。l1C -，D ： ---(-A-E-/-2 -： ：B：D：/：4：)f c)JA /4e ( B/2) (6) -，由方程组 (6)可以得到用最小二乘法拟合出的抛物线的四个参数。

2 计算各测量点到最dx-乘抛物线的法向距离图1 平面抛物线及测量点与最小二乘抛物线关系图2-1计算法向距离如图1所示，我们定义点Mi(Xi， )为过抛物线外-点 (x ，Y )的抛物线法线与拟合抛物线的交点，由于交点gi(Xi， )既在法线上又在最小二乘表1 测量数据 (mm)1'241 第35卷 第5期 2013-05(上)务l 訇 化抛物线隅 上，则点Mi(Xi， )满足方程组 (7)： -yt)t2似 B 七D)- t-xt) 2cr，E) 01 丑 c 居 F0 (7)求解非线性方程组 (7)求出各个满足条件的( ， )(i-1，2，3Ⅳ)的值并带入式 (8)计算各个测量点到抛物线的法向距离d(f)。

2.2判断测量点 位于最小二乘抛物线的内外侧通过将测量点到抛物线的代数距离的正负来判断测量点位于拟合曲线的内测还是外侧。

将 测 量 点 (x ，Y )的值 带 入 到 拟 合 曲 线F(x， )Ax Bxyc D 点 F中。

当F(x， )>0，则测量点 (x ，Y )在抛物线的外侧 ；F(x，Y)<0，则测量点 (xi，Y )在抛物线的内侧。

抛物线轮廓度误差是指实际被测轮廓线对其理 想 轮 廓 线 的 变 动 量 。 则 每 - 个 测 量 点( ，Y )(j1，2，3Ⅳ)到拟合最小二乘抛物线的法相距离为：d(i)±√( - ) ( - ) f81当测量点位于最小二乘抛物线外侧时， (f)取正值；当测量点位于最小二乘抛物线内侧时，d(O取负值。

3 抛物线的误差抛物线轮廓度公差带是与理想抛物线等距的两个抛物线等距线之间的区域。也即理想抛物线外侧测量点的最小法向距离中的最大值与抛物线内侧测量点的最小法向距离的最小值之间的距离：Am p(f)-mind(i)4 实例验证为了检验算法的正确性和可靠性，用计算机模拟发生具有不同几何位置的抛物线数据，用本算法对其进行形状误差评定，本文采用matlab软件编程计算，由于没有已有的参考数据进行比较，所 以 在 - 组 标 准 抛 物 线 中 加 入 测 量 误 差 ±0.05m的数据进行仿真验证，测量数据处理结果如表2所示。

表2 数据处理结果误差评定方法 误差 A/mm理想误差本文方法0.10.1O85 结论本文结合最小二乘原理，通过利用抛物线方程系数限制条件得到平面任意位置抛物线的拟合方程 ，并通过抛物线方程及抛物线法线方程特点，找到测量点沿法线方向到最小二乘抛物线的距离，实现对平面任意位置抛物线的最小二乘误差评定；本文方法适用于平面二次曲线并且无需进行坐标变换。