曲线的圆弧拟合

曲线的圆弧拟合

—数学应用于实践之一

一、问题的提出

在实践中常出现需要将曲线拟合成圆弧的场合，例如数控机床通常只能作直线、圆弧或圆柱螺旋线的运动，因此必须把不同曲线轨迹转化成机床运动能够接受的形式。我们可以把直线看着为半径值非常大的圆弧，而圆柱螺旋线在圆柱底面的投影就是一段圆弧，因此下面着重由简到繁地介绍曲线拟合成圆弧的几种方法。

二、椭圆曲线的拟合

椭圆曲线是一种简单常见的曲线。现以椭圆曲线长轴为对称轴，取曲线的一半。这部分曲线可以用3圆弧法或5圆弧法拟合。这部分曲线拟合后，另部分曲线以长轴为对称，其拟合结果也容易得到了。

⑴ 3圆弧法

如图1示，3圆弧法用3段相切圆弧拟合椭圆曲线段

。

设椭圆长半轴为a ，短半轴为b 。,则各圆弧半径计算公式如下：

R 1=a

b a b a b a 2)(2

222+--+

R 2=

b b

b a b a b a ++-+-2)(2

222 R 3=R 1

各圆心坐标为：)0,();,0();0,(132211a R O R b O R a O ---

用3圆弧法拟合椭圆曲线，计算方法简单，拟合圆弧段少，但对于长、短轴长度相差较大的椭圆曲线，拟合精度降低，如采用5圆弧法拟合，可以取得比较好的效果。

⑵ 5圆弧法

如图2示，5圆弧法用5段相切圆弧拟合椭圆曲线段。

同样设椭圆的长半轴为a ,短半轴为b 。 各圆弧半径计算公式如下： R 1=a b /2; R 2=ab

R 3=b a /2; R 4= R 2 R 5= R 1

各圆心坐标为：

)

0,();,();,0();,();0,(152243322211a R O y x O R b O y x O R a O ----

从图2知：23213123321221)()(;;R b R a O O R R O O R R O O -+-=-=-=

∴υ=arc cos(32312

2

12322312O O O O O O O O O O ⋅-+)

ω=arc tan(

b

R R a --31

) )sin()(232υω--=R R x b R R R y +---=3232)cos()(υω

与3圆弧法相比，5圆弧拟合比3圆弧更接近理论曲线，因此5圆弧法有较高的拟合精度。

三、复杂曲线的拟合

在这里复杂曲线是指非圆函数曲线和列表曲线。非圆函数曲线通过计算可转化成列表曲线。列表曲线由一系列有序点列P 1（x 1，y 1），P 2（x 2，y 2），…，P i （x i ，y i ），…，P n （x n ，y n ）组成，可列成表格形式。列表曲线的数据也可以通过检测工具对实物逐点测量后获得。

下面介绍两种较常见的列表曲线拟合方法：

⑴双圆弧法

双圆弧处理方法依次取列表曲线的4个点，如图3所示。

要确定P2P3之间的圆弧段P2T和TP3，必须给出P1，P2，P3和P4共4个点。

P1～

P4点的分布有多种情况，但经过适当处理后可归纳成图4所示的两种情况。

情况①是P1，P4点在P2，P3点连线的同侧，所拟合的圆弧也在连线的同侧，称为同侧圆弧拟合（见图5）。情况②的P1，P4点在P2，P3点连线的异侧，所拟合的两条圆弧彼此外切，称为异侧圆弧拟合。

同侧圆弧拟合的计算公式推导如下：

如图5示，过P 2点作∠P 1P 2P 3的角平分线P 2A 2的垂线P 2M ，过P 3点作∠P 2P 3P 4的角平分线P 3A 3的垂线P 3M ，过P 2点再作∠P 3P 2M 的角平分线P 2N ，过P 3点作∠P 2P 3M 的角平分线P 3N ，过P 2N 与P 3N 的交点N 作P 2P 3的垂线NP ，NP 与P 2A 2和P 3A 3分别交于O 2、O 3点，分别以O 2、O 3为圆心，O 2P 2、O 3P 3为半径作圆弧C 2、C 3，2圆弧相切于N 点，且P 2M 和P 3M 也分别在P 2点和P 3点与C 2、C 3相切，则弧P 2N 和弧NP 3构成了P 2、P 3给定点之间的逼近线。

设C 2、C 3圆弧的半径分别为r 2、r 3，则可推导出：

h v u v d v u v d r /)]sin (cos sin )cos (sin cos [2--+= h v q v d v q v d r /)]sin (cos sin )cos (sin cos [3--+=

式中，)sin )(cos cos (sin )cos )(sin sin (cos v u v q v u v q h -+-+-=；d 是P 2点和P 3点之间的距离；q 、u 和v 分别是矢量22A P 、33A P 及32P P 与x 轴正向之间的夹角。

因此，由给出的P 1～P 4点座标值，可以确定d 、q 、u 、v 的值，这样就可以确定r 2、r 3。根据已知数据再确定圆心O 2、O 3及N 点的座标值。

异侧圆弧计算公式的推导类似同侧圆弧，在此不再赘述。

双圆弧拟合法可保证圆弧逼近线过给出点，并在各点一阶导数连续。它的计算方法简单，适应性强，有较好的拟合精度。缺点是拟合后的圆弧段数较多，拟合精度不能调控。

⑵最小二乘法

与双圆弧拟合法不同，它拟合成的近似圆弧不一定要通过各给出点。如图6示，近似圆弧R 与每个给出点之间绝对误差的平方和为最小。

圆弧与点准确值的误差可以人为地用“容差”ta 来限定，如果所有被拟合的点均在R ±ta 的“容差”带范围内，则认为R 圆弧拟合成功。

“容差”是根据实际需要而人为设定的，一般选用很小的数值，如0.01mm 。如果“容差”越小，则逼近精度就越高，每条圆弧段能满足精度的点减少，因此总的圆弧段数就会增加。反之，则逼近精度降低，每段圆弧能满足精度的点增加，