曲线的圆弧拟合.docVIP

曲线的圆弧拟合 —数学应用于实践之一 一、问题的提出 在实践中常出现需要将曲线拟合成圆弧的场合，例如数控机床通常只能作直线、圆弧或圆柱螺旋线的运动，因此必须把不同曲线轨迹转化成机床运动能够接受的形式。我们可以把直线看着为半径值非常大的圆弧，而圆柱螺旋线在圆柱底面的投影就是一段圆弧，因此下面着重由简到繁地介绍曲线拟合成圆弧的几种方法。 二、椭圆曲线的拟合 椭圆曲线是一种简单常见的曲线。现以椭圆曲线长轴为对称轴，取曲线的一半。这部分曲线可以用3圆弧法或5圆弧法拟合。这部分曲线拟合后，另部分曲线以长轴为对称，其拟合结果也容易得到了。 ⑴ 3圆弧法 如图1示，3圆弧法用3段相切圆弧拟合椭圆曲线段。 设椭圆长半轴为，短半轴为b。,则各圆弧半径计算公式如下： R1= R2= R3=R1 各圆心坐标为： 用3圆弧法拟合椭圆曲线，计算方法简单，拟合圆弧段少，但对于长、短轴长度相差较大的椭圆曲线，拟合精度降低，如采用5圆弧法拟合，可以取得比较好的效果。 ⑵ 5圆弧法 如图2示，5圆弧法用5段相切圆弧拟合椭圆曲线段。 同样设椭圆的长半轴为,短半轴为b。 各圆弧半径计算公式如下： R1=; R2= R3=; R4= R2 R5= R1 各圆心坐标为： 从图2知： ∴υ=arc cos() ω=arc tan() 与3圆弧法相比，5圆弧拟合比3圆弧更接近理论曲线，因此5圆弧法有较高的拟合精度。 三、复杂曲线的拟合 在这里复杂曲线是指非圆函数曲线和列表曲线。非圆函数曲线通过计算可转化成列表曲线。列表曲线由一系列有序点列P1（x1，y1），P2（x2，y2），…，Pi（xi，yi），…，Pn（xn，yn）组成，可列成表格形式。列表曲线的数据也可以通过检测工具对实物逐点测量后获得。 下面介绍两种较常见的列表曲线拟合方法： ⑴双圆弧法 双圆弧处理方法依次取列表曲线的4个点，如图3所示。 要确定P2P3之间的圆弧段P2T和TP3，必须给出P1，P2，P3和P4共4个点。P1～P4点的分布有多种情况，但经过适当处理后可归纳成图4所示的两种情况。 情况①是P1，P4点在P2，P3点连线的同侧，所拟合的圆弧也在连线的同侧，称为同侧圆弧拟合（见图5）。情况②的P1，P4点在P2，P3点连线的异侧，所拟合的两条圆弧彼此外切，称为异侧圆弧拟合。 同侧圆弧拟合的计算公式推导如下： 如图5示，过P2点作∠P1P2P3的角平分线P2A2的垂线P2M，过P3点作∠P2P3P4的角平分线P3A3的垂线P3M，过P2点再作∠P3P2M的角平分线P2N，过P3点作∠P2P3M的角平分线P3N，过P2N与P3N的交点N作P2P3的垂线NP，NP与P2A2和P3A3分别交于O2、O3点，分别以O2、O3为圆心，O2P2、O3P3为半径作圆弧C2、C3，2圆弧相切于N点，且P2M和P3M也分别在P2点和P3点与C2、C3相切，则弧P2N和弧NP3构成了P2、P3给定点之间的逼近线。 设C2、C3圆弧的半径分别为r2、r3，则可推导出： 式中，；d是P2点和P3点之间的距离；q、u和v分别是矢量、及与x轴正向之间的夹角。 因此，由给出的P1～P4点座标值，可以确定d、q、u、v的值，这样就可以确定r2、r3。根据已知数据再确定圆心O2、O3及N点的座标值。 异侧圆弧计算公式的推导类似同侧圆弧，在此不再赘述。 双圆弧拟合法可保证圆弧逼近线过给出点，并在各点一阶导数连续。它的计算方法简单，适应性强，有较好的拟合精度。缺点是拟合后的圆弧段数较多，拟合精度不能调控。 ⑵最小二乘法 与双圆弧拟合法不同，它拟合成的近似圆弧不一定要通过各给出点。如图6示，近似圆弧R与每个给出点之间绝对误差的平方和为最小。 圆弧与点准确值的误差可以人为地用“容差”ta来限定，如果所有被拟合的点均在R±ta的“容差”带范围内，则认为R圆弧拟合成功。 “容差”是根据实际需要而人为设定的，一般选用很小的数值，如0.01mm。如果“容差”越小，则逼近精度就越高，每条圆弧段能满足精度的点减少，因此总的圆弧段数就会增加。反之，则逼近精度降低，每段圆弧能满足精度的点增加，因此总的圆弧段数就会减少。 最小二乘法的计算方法如下： 如图6示，近似圆弧的圆心坐标为，半径为R，为曲线上任一点的坐标，则可求得： ； R= (推导过程略，请参阅最小二乘法有关资料) 通过一组列表曲线的拟合可得到若干条拟合圆弧。前面曾提