CN101751006A - 一种插补方法 - Google Patents

Abstract

本发明公开一种插补方法，包括步骤为：针对正弦函数类所需路径或轮廓线，包括正弦曲线、椭圆曲线、圆弧曲线，提出一种通过算术运算的递推公式准确确定曲线中间点位置坐标值增量的插补；并提出设定替代曲线，以提高插补的拟合精度或减小插补运算量。所述方法运算简单，对于同样的分段段数，完成运算所需时间少，相当于提高了插补运算速度；或者，在同等的运算时间内，允许增加曲线分段段数，从而提高折线拟合曲线的精度。由于运算简单，这种方法可在单片计算机上实现，从而降低了插补器的成本。因此，这种插补方法速度快，精度高，装置成本低。

Description

一种插补方法

【技术领域】

本发明属于计算机数控领域，尤其是涉及一种控制物体运动轨迹的插补方法。

【背景技术】

1)插补的任务

数控系统广泛应用于机械运动轨迹的控制，可以控制机床、机器人、缝纫机、焊接机等的运动轨迹。所需路径或轮廓线的“插补”是数控技术的关键。

插补的任务就是在所需路径或轮廓线Q的二个已知点间插入若干个中间点，并确定所述中间点的位置坐标值。

插补所得结果将依次送给直线插补器；对应一组相邻二点的位置坐标值，直线插补器产生一组分布均匀的脉冲序列，并通过步进马达或伺服系统控制受控对象运动，产生一个直线段的运动轨迹。或者，插补所得结果直接以数值方式依次送给伺服系统，控制受控对象运动；对应一组相邻二点的位置坐标值，产生一个直线段的运动轨迹。受控对象整个运动轨迹将是一个由上述首尾相接的直线段构成的折线，折线的起点、终点、中间点与所需路径或轮廓线Q的对应起点、终点、中间点重合。换句话说，在所需路径或轮廓线Q的二个端点间插入若干个中间点，将曲线Q上相邻二点以直线段连接，这些首尾相接的直线段构成的折线就是受控对象的运动轨迹。

插补的目的就是确定受控对象的运动轨迹，使之尽可能接近所需路径或轮廓线Q；或者说，插补的目的就是确定所述的折线，也即受控对象的运动轨迹，以此拟合所需路径或轮廓线Q，且拟合误差不超过允许值。

2)现状

对于常见的正弦曲线、椭圆、圆弧等所需路径或轮廓线，准确确定所述中间点位置坐标值的运算涉及三角函数等比较复杂的计算。这样的计算，靠单片计算机很难完成，而由PC级计算机采用高级语言不难完成。在运动控制装置中，如果作为整个装置的控制核心配置PC级计算机，直接由之承担插补的实时运算将占用大量资源，从而影响整个装置工作；如果专为插补工作配置PC级计算机，将导致装置成本提高。

当前，运动控制技术已广泛应用于各领域，甚至将进入普通家庭，如家用机器人等；基于现有插补方法的产品难以适应、满足市场的需求。

可是，可以在单片计算机上快速、准确实现的插补方法，至今未见在公开资料上发表或在产品中使用。

【发明内容】

本发明的目的是针对正弦函数类所需路径或轮廓线，包括正弦曲线、椭圆曲线、圆弧曲线，提出一种通过算术运算的递推公式准确确定曲线中间点位置坐标值增量的插补方法，并且还提出一种设定替代曲线的方法，以提高拟合精度或减小插补运算量。所述方法运算简单，对于同样的分段段数，完成运算所需时间少，相当于提高了插补运算速度；或者，在同等的运算时间内，允许增加曲线分段段数，从而提高折线拟合曲线的精度。由于运算简单，这种方法可在单片计算机上实现，从而降低了插补器的成本。因此，这种插补方法速度快，精度高，装置成本低。

本发明的目的是按如下技术方案实现的：

1、本发明所述的一种插补方法，其所针对的所需路径或轮廓线Q的位置坐标Ω_P(p＝1、2、3、……、m_Ω)中包括有一个或若干个坐标Ψ_k(k＝1、2、3、……、m_Ψ)，对应所述一个或若干个坐标Ψ_k(k＝1、2、3、……、m_Ψ)的坐标函数分别可以表示为以参数t为自变量、幅值分别为H_k(k＝1、2、3、……m_Ψ)、初始相位分别为α_k(k＝1、2、3、……、m_Ψ)、周期相同为(2π/ω)的正弦函数，其表达式为

Ψ_k(t)＝H_k sin(ωt+α_k)(k＝1、2、3、……、m_Ψ)，(1-1)

所述曲线Q对应坐标Ψ_k(k＝1、2、3、……、m_Ψ)的坐标函数指的是描述所述曲线Q上的点的位置坐标Ψ_k(k＝1、2、3、……、m_Ψ)变化的以参数t为自变量的函数，

所述参数t可以是该曲线Q位置坐标Ω_P(p＝1、2、3、……、m_Ω)中的一个坐标，也可以是这些坐标之外的另一个参数，

所述的所需路径或轮廓线Q指的是二个已知点间的曲线段，所述曲线Q各坐标函数的定义域相同，定义域的二个端点分别与曲线Q二个已知点对应着相同的t值。

某一t值所对应的所需路径或轮廓线Q上的点的位置坐标值，就是同一t值所对应的所述坐标函数的函数值。针对所需路径或轮廓线Q的插补也就是针对其坐标函数Ω_P(t)(p＝1、2、3、……、m_Ω)的插补，即在其各个坐标函数Ω_P(t)(p＝1、2、3、……、m_Ω)的定义域二个端点间插入若干个中间点并确定对应所述中间点的坐标函数值。

申请人注意到，对所需路径或轮廓线Q进行插补时，为了提高拟合精度，或减小插补运算量，可以设定相应的替代曲线Q_δ；在替代曲线Q_δ的二个端点间插入若干个中间点，将替代曲线Q_δ上相邻的二个点以直线段连接，用这些首尾相接的直线段构成的折线拟合二个已知点间的所需路径或轮廓线Q。为此，需通过插补确定所述替代曲线Q_δ的中间点及其位置坐标值。

替代曲线Q_δ的坐标函数Ω_δP(t)(p＝1、2、3、……、m_Ω)称为替代坐标函数，也称为所述曲线Q的坐标函数Ω_P(t)(p＝1、2、3、……、m_Ω)的替代坐标函数。确定替代曲线Q_δ中间点及其位置坐标值，就是确定相应替代坐标函数Ω_δP(t)(p＝1、2、3、……、m_Ω)定义域的中间点及其函数值。插补中确定的中间点将所述定义域分成分段，每个分段定义域将对应一个线性函数，所述各个分段定义域的端点所对应的线性函数值与替代坐标函数值相等，整个定义域将对应一个由各个分段定义域对应的线性函数组成的分段线性函数Ψ_δLk(t)(k＝1、2、3、……、m_Ψ)，而坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)将以分段线性函数Ψ_δLk(t)(k＝1、2、3、……、m_Ψ)拟合。

申请人还注意到，所述曲线Q或Q_δ二个端点间的中间点的位置坐标值，可依据与之相邻的前一个中间点或所述曲线的起始端点其位置坐标值与位置坐标值增量之和决定，因此，插补中确定了中间点的位置坐标值增量，也就确定了相应中间点的位置坐标值。或者说，插补中确定了定义域中间点坐标函数值增量，也就确定了定义域相应中间点的坐标函数值。

针对坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)的插补步骤包括，

(1)设定对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)，

(2)确定替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)定义域二个端点间的中间点，包括，

①确定所述中间点所对应的参数t的值或其增量Δt的值，

②确定所述中间点的个数，

(3)确定所述中间点的替代坐标函数值或其增量值，

(4)存储/输出运算结果，

上述插补步骤也适用于针对坐标函数Ω_P(t)(p＝1、2、3、……、m_Ω)的插补。

本发明提出的插补方法，其特征在于：

(1)对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的替代坐标函数设定为与相应的坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)具有相同周期(2π/ω)、相同初始相位α_k(k＝1、2、3、……、m_Ψ)、相同定义域但不同幅值的正弦函数，其表达式为

Ψ_δk(t)＝H_δk sin(ωt+α_k)(k＝1、2、3、……、m_Ψ)，(1-2)

式中，H_δk＝H_k+δ_k，(k＝1、2、3、……、m_Ψ)，(1-3)

其中δ_k(k＝1、2、3、……、m_Ψ)为替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)幅值H_δk(k＝1、2、3、……、)相对相应的坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)幅值H_k(k＝1、2、3、……、)的幅值差，或说是替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)的幅值差，或说是对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的幅值差，

δ_k＝H_δk-H_k，(k＝1、2、3、……、m_Ψ)，(1-4)

δ_k≥0，(k＝1、2、3、……、m_Ψ)，(1-5)

δ_k的数值可根据第4至11点所述的方法确定。

(2)将对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的各替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)定义域等分，以等分分段的交点作为定义域的中间点，这些中间点所对应的替代坐标函数值增量ΔΨ_δk(t)(k＝1、2、3、……、m_Ψ)，按下述公式确定：

ΔΨ_δΛ(t_u+1)＝ΔΨ_δΛ(t_u)cosΔT_Λ+ΔΦ_δΛ(t_u+1)sinΔT_Λ，(1-6)

ΔΦ_δΛ(t_u+1)＝ΔΦ_δΛ(t_u)cosΔT_Λ-ΔΨ_δΛ(t_u)sinΔT_Λ，(1-7)

式中，①Ψ_Λ表示所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)中的某个坐标，Λ为序号k(k＝1、2、3、……、m_Ψ)中的某个序号，

(②u+1为替代坐标函数Ψ_δΛ(t)定义域二个端点间的某个中间点的序号，u就是与之相邻的前一个中间点或定义域的起点的序号，

以n_Λ表示替代坐标函数Ψ_δΛ(t)定义域等分分段的段数，以i_Λ(i_Λ＝1、2、3、……、n_Λ)作为分段的序号，以i_Λ(i_Λ＝1、2、3、……、n_Λ，n_Λ+1)作为包括定义域的二个端点及二个端点间的中间点在内的点的序号，二个端点分别对应着序号1及n_Λ+1，二个端点间的中间点分别对应着序号i_Λ(i_Λ＝2、3、……、n_Λ)，u+1就是序号i_Λ(i_Λ＝2、3、……、n_Λ)中的某一个序号，u就是序号i_Λ(i_Λ＝1，2、3、……、n_Λ)中与u+1相邻的前一个中间点或定义域起点的序号，u+2就是序号i_Λ(i_Λ＝1，2、3、……、n_Λ，n_Λ+1)中与u+1相邻的后一个中间点或定义域终点的序号，

所述定义域起点指的是与所述替代曲线Q_δ起点对应的定义域的端点，所述定义域的终点指的是与所述替代曲线Q_δ终点对应的定义域的端点，

③t_u+1为替代坐标函数Ψ_δΛ(t)定义域二个端点间序号为u+1的中间点所对应的参数t的值，

t_u为与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的参数t的值，

t_u+2为与序号为u+1的中间点相邻的后一个中间点或定义域终点所对应的参数t的值，

t₁为替代坐标函数Ψ_δΛ(t)定义域起点所对应的参数t的值，

t(n_Λ+1)为替代坐标函数Ψ_δΛ(t)定义域终点所对应的参数t的值，

④ΔT_Λ为替代坐标函数Ψ_δΛ(t)定义域各等分分段起点所对应的参数t的等效增量，或说是替代坐标函数Ψ_δΛ(t)定义域中序号为i_Λ(i_Λ＝1、2、3、……、n_Λ)的点所对应的参数t的等效增量，或说是替代坐标函数Ψ_δΛ(t)定义域各等分分段所对应的参数t的等效增量，其数值为所述参数t的增量Δτ_Λ的ω倍，是常数，

ΔT_Λ＝ωΔτ_Λ＝常数，(1-8)

(i_Λ＝1、2、3、……、n_Λ)，(1-9)

其中，t(i_Λ+1)(i_Λ＝1、2、3、……、n_Λ)为序号为(i_Λ+1)(i_Λ＝1、2、3、……、n_Λ)的点对应的参数t的值，

t(i_Λ)(i_Λ＝1、2、3、……、n_Λ)为序号为i_Λ(i_Λ＝1、2、3、……、n_Λ)的点对应的参数t的值，

⑤ΔΨ_δΛ(t_u+1)为替代坐标函数Ψ_δΛ(t)定义域中序号为u+1的中间点所对应的替代坐标函数值的增量，

ΔΨ_δΛ(t_u+1)＝Ψ_δΛ(t_u+2)-Ψ_δΛ(t_u+1)，(1-10)

其中，Ψ_δΛ(t_u+2)为与序号为u+1的中间点相邻的后一个中间点或定义域终点所对应的替代坐标函数值，

Ψ_δΛ(t_u+1)为序号为u+1的中间点所对应的替代坐标函数值，

⑥ΔΨ_δΛ(t_u)为替代坐标函数Ψ_δΛ(t)定义域中与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的替代坐标函数值增量，

ΔΨ_δΛ(t_u)＝Ψ_δΛ(t_u+1)-Ψ_δΛ(t_u)，(1-11)

其中，Ψ_δΛ(t_u+1)为序号为u+1的中间点所对应的替代坐标函数值，

Ψ_δΛ(t_u)为与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的替代坐标函数值，

⑦ΔΦ_δΛ(t_u+1)为替代坐标函数Ψ_δΛ(t)定义域中序号为u+1的中间点所对应的虚拟替代坐标函数值的增量，

ΔΦ_δΛ(t_u+1)＝Φ_δΛ(t_u+2)-Φ_δΛ(t_u+1)，(1-12)

其中，Φ_δΛ(t_u+2)为与序号为u+1的中间点相邻的后一个中间点或定义域终点所对应的虚拟替代坐标函数值，

Φ_δΛ(t_u+1)为序号为u+1的中间点所对应的虚拟替代坐标函数值，

⑧ΔΦ_δΛ(t_u)为替代坐标函数Ψ_δΛ(t)定义域中与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的虚拟替代坐标函数值增量，

ΔΦ_δΛ(t_u)＝Φ_δΛ(t_u+1)-Φ_δΛ(t_u)，(1-13)

其中，Φ_δΛ(t_u+1)为序号为u+1的中间点所对应的虚拟替代坐标函数值，

Φ_δΛ(t_u)为与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的虚拟替代坐标函数值，

⑨sinΔT_Λ为对应ΔT_Λ的正弦函数值，cosΔT_Λ为对应ΔT_Λ的余弦函数值，

上述⑦、⑧中所述的虚拟替代坐标函数Φ_δΛ(t)是为简化公式(1-6)、(1-7)的表述设定的函数，Ψ_δΛ(t)与Φ_δΛ(t)是一组具有相同幅值H_δΛ、相同周期(2π/ω)、相同初始相位α_Λ及相同定义域的正弦函数与余弦函数，其表达式为

Ψ_δΛ(t)＝H_δΛsin(ωt+α_Λ)，(1-14)

Φ_δΛ(t)＝H_δΛcos(ωt+α_Λ)。，，(1-15)

Φ_δΛ(t)或说是对应坐标函数Ψ_Λ(t)设定的虚拟坐标函数Φ_Λ(t)的替代坐标函数，Ψ_Λ(t)与Φ_Λ(t)是一组具有相同幅值H_Λ、相同周期(2π/ω)、相同初始相位α_Λ及相同定义域的正弦函数与余弦函数，其表达式为

Ψ_Λ(t)＝H_Λsin(ωt+α_Λ)，(1-16)

Φ_Λ(t)＝H_Λcos(ωt+α_Λ)，(1-17)

其中H_δΛ＝H_Λ+δ_Λ， (1-18)

说明：

(1)参见图2。如果所需路径或轮廓线的位置坐标只包括了二个个坐标t和Ψ_A，且坐标轴t和Ψ_A构成了直角坐标系tOΨ_A；则所需路径或轮廓线就是平面tOΨ_A上幅值为H_A、初始相位为α_A、周期为(2π/ω)的的正弦曲线段Z_A，其替代曲线就是平面tOΨ_A上幅值为H_δA、初始相位为α_A、周期为(2π/ω)的替代正弦曲线段Z_δA。

(2)参见图5。以序号为A的坐标函数Ψ_A(t)为例，坐标轴Ψ_A和Φ_A构成了虚拟直角坐标系Φ_AOΨ_A，坐标函数Ψ_A(t)与虚拟坐标函数Φ_A(t)在Φ_AOΨ_A虚拟平面上的图形即为圆心在坐标轴原点O、半径为H_A的圆弧段C_A。圆弧段C_A上的点与圆心O的连线相对轴Φ_A的夹角即为(ωt+α_A)。替代坐标函数Ψ_δA(t)与虚拟替代坐标函数Φ_δA(t)在Φ_AOΨ_A虚拟平面上的图形即为圆心在坐标轴原点O、半径为H_δA的替代圆弧段C_δA。圆弧段C_δA上的点与圆心O的连线相对轴Φ_A的夹角即为(ωt+α_A)。替代坐标函数Ψ_δA(t)定义域的端点与中间点对应着圆弧段C_δA相应的端点与中间点。图中C_A、C_δA分别是第一象限的圆弧段。

(3)公式(1-6)、(1-7)的依据是：

参见图6。以序号为A的坐标函数Ψ_A(t)为例，替代坐标函数Ψ_δA(t)定义域中序号为i(i＝2、3、……、n)的中间点对应着圆弧段C_δA上相应的中间点。这些中间点将圆弧段C_δA等分成序号为i(i＝1、2、3、……、n)n个分段。圆弧段C_δA的等分分段对应的圆心角即为ΔT_A。以l_δA表示圆弧段C_δA一个分段的内接弦的长度，对应于t_u有

${\mathrm{\Δ\Ψ}}_{\mathrm{\δA}}\left(t_u\right)=l_{\mathrm{\δA}}\cos [\mathrm{\ω}(t_u+\frac{\mathrm{\Δ}{\mathrm{\τ}}_A}{2})+{\mathrm{\α}}_A]=l_{\mathrm{\δA}}\cos ({\mathrm{\ωt}}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A),---(1-19)$

${\mathrm{\Δ\Φ}}_{\mathrm{\δA}}\left(t_u\right)=-l_{\mathrm{\δA}}\sin [\mathrm{\ω}(t_u+\frac{\mathrm{\Δ}{\mathrm{\τ}}_A}{2})+{\mathrm{\α}}_A]=-l_{\mathrm{\δA}}\sin ({\mathrm{\ωt}}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A),---(1-20)$

式中，ΔT_A＝ωΔτ_A。(1-21)

以u+1表示所述的某个中间点的序号，则

${\mathrm{\Δ\Ψ}}_{\mathrm{\δA}}\left(t_{u+1}\right)=l_{\mathrm{\δA}}\cos (\mathrm{\ω}{t}_u+\frac{{\mathrm{\ΔT}}_A}{2}+{\mathrm{\ΔT}}_A+{\mathrm{\α}}_A)=$

$=l_{\mathrm{\δA}}\cos (\mathrm{\ω}{t}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A)\cos \mathrm{\Δ}{T}_A-l_{\mathrm{\δA}}\sin (\mathrm{\ω}{t}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A)\sin \mathrm{\Δ}{T}_A,$

$=\mathrm{\Δ}{\mathrm{\Ψ}}_{\mathrm{\δA}}\left(t_u\right)\cos \mathrm{\Δ}{T}_A+\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{\δA}}\left(t_u\right)\sin \mathrm{\Δ}{T}_A.---(1-22)$

${\mathrm{\Δ}\mathrm{\Φ}}_{\mathrm{\δA}}\left(t_{u+1}\right)=-l_{\mathrm{\δA}}\sin (\mathrm{\ω}{t}_u+\frac{\mathrm{\Δ}{T}_A}{2}+\mathrm{\Δ}{T}_A+{\mathrm{\α}}_A)=$

$=-l_{\mathrm{\δA}}\sin (\mathrm{\ω}{t}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A)\cos \mathrm{\Δ}{T}_A-l_{\mathrm{\δA}}\cos (\mathrm{\ω}{t}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A)\sin \mathrm{\Δ}{T}_A,$

$=\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{\δA}}\left(t_u\right)\cos \mathrm{\Δ}{T}_A-\mathrm{\Δ}{\mathrm{\Ψ}}_{\mathrm{\δA}}\left(t_u\right)\sin \mathrm{\Δ}{T}_A.---(1-23)$

(4)公式(1-6)、(1-7)中所述正弦函数值sinΔT_Λ与余弦函数值cosΔT_Λ，可按下述公式确定，

$\sin \mathrm{\Δ}{T}_{\mathrm{\Λ}}=\mathrm{\Δ}{T}_{\mathrm{\Λ}}-\frac{1}{6}\×{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+......+\frac{1}{v^{\′}!}{\left[\frac{d^{v^{\′}}}{{\mathrm{dt}}^{v^{\′}}}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_S,---(1-24)$

$\cos \mathrm{\Δ}{T}_{\mathrm{\Λ}}=1-\frac{1}{2}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+\frac{1}{24}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^4+......+\frac{1}{v^{\′\′}!}{\left[\frac{d^{v^{\′\′}}}{{\mathrm{dt}}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′\′}}+R_C,---(1-25)$

其中， $R_S=\frac{1}{(v^{\′}+1)!}{\left[\frac{d^{v^{\′}+1}}{{\mathrm{dt}}^{v^{\′}+1}}\sin t\right]}_{t={\mathrm{\γ}}^{\′}}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}+1},---(1-26)$

$R_C=\frac{1}{(v^{\′\′}+1)!}{\left[\frac{d^{v^{\′\′}+1}}{{\mathrm{dt}}^{v^{\′\′}+1}}\cos t\right]}_{t={\mathrm{\γ}}^{\′\′}}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′\′}+1},---(1-27)$

式中，①v′+2为正整数，数值等于等式(1-24)右边部分的项数，

v′≥0，(1-28)

②γ′的取值范围为

0≤γ′≤ΔT_Λ，(1-29)

③v″+2为正整数，数值等于等式(1-25)右边部分的项数

v″≥0，(1-30)

④γ″的取值范围为

0≤γ″≤ΔT_Λ。(1-31)

按公式(1-24)、(1-25)确定sinΔT_Λ、cosΔT_Λ时，可以将R_S、R_C取值为0。sinΔT_Λ、cosΔT_Λ数值精度由舍去的R_S、R_C数值决定，而R_S、R_C数值由ΔT_Λ及v′、v″决定。v′、v″的选择可参见实施例。

公式(1-24)、(1-25)的依据是：

$\sin \mathrm{\Δ}{T}_{\mathrm{\Λ}}=\sin (0+\mathrm{\Δ}{T}_{\mathrm{\Λ}})=$

$={\left[\sin t\right]}_{t=0}+{\left[\frac{d}{\mathrm{dt}}\sin t\right]}_{t=0}\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)+\frac{1}{2}{\left[\frac{d^2}{{\mathrm{dt}}^2}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+\frac{1}{3!}{\left[\frac{d^3}{{\mathrm{dt}}^3}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+$

$+......+\frac{1}{v^{\′}!}{\left[\frac{d^{v^{\′}}}{{\mathrm{dt}}^{v^{\′}}}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_S=$

$=\sin 0+\left(\cos 0\right)\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)-\frac{1}{2}\left(\sin 0\right){\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2\frac{1}{3!}\left(\cos 0\right){\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+......+$

$+\frac{1}{v^{\′}!}{\left[\frac{d^{v^{\′}}}{{\mathrm{dt}}^{v^{\′}}}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_S=$

$=0+\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)-0-\frac{1}{6}\×{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+0+......+\frac{1}{v^{\′}!}{\left[\frac{d^{v^{\′}}}{{\mathrm{dt}}^{v^{\′}}}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_S,---(1-32)$

即 $\sin \mathrm{\Δ}{T}_{\mathrm{\Λ}}=\mathrm{\Δ}{T}_{\mathrm{\Λ}}-\frac{1}{6}\×{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+......+\frac{1}{v^{\′}!}{\left[\frac{d^{v^{\′}}}{{\mathrm{dt}}^{v^{\′}}}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_S,---(1-24)$

其中， $R_S=\frac{1}{(v^{\′}+1)!}{\left[\frac{d^{v^{\′}+1}}{{\mathrm{dt}}^{v^{\′}+1}}\sin t\right]}_{t={\mathrm{\γ}}^{\′}}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}+1},---(1-26)$

$\cos \mathrm{\Δ}{T}_{\mathrm{\Λ}}=\cos (0+\mathrm{\Δ}{T}_{\mathrm{\Λ}})=$

$={\left[\cos t\right]}_{t=0}+{\left[\frac{d}{\mathrm{dt}}\cos t\right]}_{t=0}\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)+\frac{1}{2}{\left[\frac{d^2}{{\mathrm{dt}}^2}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+\frac{1}{3!}{\left[\frac{d^3}{{\mathrm{dt}}^3}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+$

$+......+\frac{1}{v^{\′\′}!}{\left[\frac{d^{v^{\′\′}}}{{\mathrm{dt}}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_C=$

$=\cos 0-\left(\sin 0\right)\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)-\frac{1}{2}\left(\cos 0\right){\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+\frac{1}{3!}{\left(\sin 0\right)\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+......+$

$+\frac{1}{v^{\′\′}!}{\left[\frac{d^{v^{\′\′}}}{{\mathrm{dt}}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_C$

$=1-0-\frac{1}{2}\×{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+0+\frac{1}{4!}\×1\≥{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^4+......+\frac{1}{v^{\′\′}!}{\left[\frac{d^{v^{\′\′}}}{{\mathrm{dt}}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′\′}}+R_C,---(1-33)$

即 $\cos \mathrm{\Δ}{T}_J=1-\frac{1}{2}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+\frac{1}{24}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^4+......+\frac{1}{v^{\′\′}!}{\left[\frac{d^{v^{\′\′}}}{{\mathrm{dt}}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′\′}}+R_C,---(1-25)$

其中， $R_C=\frac{1}{(v^{\′\′}+1)!}{\left[\frac{d^{v^{\′\′}+1}}{{\mathrm{dt}}^{v^{\′\′}+1}}\cos t\right]}_{t={\mathrm{\γ}}^{\′\′}}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′\′}+1}.---(1-27)$

(5)公式(1-6)、(1-7)中对应所述坐标Ψ_Λ的替代坐标函数Ψ_δΛ(t)定义域起点的函数值增量ΔΨ_δΛ(t₁)、ΔΦ_δΛ(t₁)，可按下述公式确定，

ΔΨ_δΛ(t₁)＝Ψ_δΛ(t₁)cosΔT_Λ+Φ_δΛ(t₁)sinΔT_Λ-Ψ_δΛ(t₁)，(1-34)

ΔΦ_δΛ(t₁)＝Φ_δΛ(t₁)cosΔT_Λ-Ψ_δΛ(t₁)sinΔT_Λ-Φ_δΛ(t₁)。(1-35)

其依据是：

ΔΨ_δΛ(t₁)＝Ψ_δΛ(t₁+Δt₁)-Ψ_δΛ(t₁)＝

            ＝H_δΛsin(ω(t₁+Δt₁)+α_Λ)-Ψ_δΛ(t₁)＝

            ＝H_δΛsin(ωt₁+α_Λ+ΔT_Λ)-Ψ_δΛ(t₁)＝

            ＝H_δΛ[sin(ωt₁+α_Λ)cosΔT_Λ+cos(ωt₁+α_Λ)sinΔT_Λ]-Ψ_δΛ(t₁)＝

            ＝Ψ_δΛ(t₁)cosΔT_Λ+Φ_δΛ(t₁)sinΔT_Λ-Ψ_δΛ(t₁)，(1-36)

ΔΦ_δΛ(t₁)＝Φ_δA(t₁+Δt₁)-Φ_δΛ(t₁)＝

            ＝H_δΛcos(ω(t₁+Δt₁)+α_Λ)-Φ_δΛ(t₁)＝

            ＝H_δΛcos(ωt₁+α_Λ+ΔT_Λ)-Φ_δΛ(t₁)＝

            ＝H_δΛ[cos(t₁+α_Λ)cosΔT_Λ-sin(t₁+α_Λ)sinΔT_Λ]-Φ_δΛ(t₁)＝

            ＝Φ_δΛ(t₁)cosΔT_Λ-Ψ_δΛ(t₁)sinΔT_Λ-Φ_δΛ(t₁)。(1-37)

2、如第1点所述的插补方法，其特点为：

(1)所述坐标Ψ_k(k＝1、……、m_Ψ)共有二个坐标Ψ_k(k＝1、2)或Ψ₁、Ψ₂，对应Ψ₁、Ψ₂的坐标函数Ψ₁(t)、Ψ₂(t)可以表示为一组幅值分别为H₁与H₂、周期相同为(2π/ω)、初始相位相同为α₀的余弦函数与正弦函数，其表达式为

${\mathrm{\Ψ}}_1\left(t\right)=H_1\sin (\mathrm{\ωt}+{\mathrm{\α}}_1)=H_1\sin (\mathrm{\ωt}+{\mathrm{\α}}_0+\frac{\mathrm{\π}}{2})=H_1\cos (\mathrm{\ωt}+{\mathrm{\α}}_0),---2-1)$

      Ψ₂(t)＝H₂sin(ωt+α₂)＝H₂sin(ωt+α₀)，(2-2)

式中， ${\mathrm{\α}}_1={\mathrm{\α}}_0+\frac{\mathrm{\π}}{2},---(2-3)$

      α₂＝α₀，(2-4)

(2)对应所述坐标Ψ₁及Ψ₂的替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)为

      Ψ_δ1(t)＝H_δ1cos(ωt+α₀)，(2-5)

      Ψ_δ2(t)＝H_δ2sin(ωt+α₀)，(2-6)

式中，H_δ1＝H₁+δ₁，(2-7)

      H_δ2＝H₂+δ₂。(2-8)

坐标函数Ψ₁(t)、Ψ₂(t)、Ψ_δ1(t)、Ψ_δ2(t)具有相同的定义域。

参见图7。如果所需路径或轮廓线的位置坐标只包括了二个个坐标Ψ₁和Ψ₂，且坐标轴Ψ₁和Ψ₂构成了直角坐标系Ψ₁OΨ₂，则所需路径或轮廓线即为Ψ₁OΨ₂平面上中心在坐标轴原点O的椭圆弧段Z_E。图中Z_E是一个第一象限的椭圆弧段。其在直角坐标系Ψ₁OΨ₂下的曲线方程为

$\frac{{{\mathrm{\Ψ}}_1}^2}{{H_1}^2}+\frac{{{\mathrm{\Ψ}}_2}^2}{{H_2}^2}=1,---(2-9)$

式中，H₁和H₂为椭圆的半轴。

椭圆弧段Z_E上的点与坐标轴原点O的连线相对Ψ₁轴的夹角即为(ωt+α₀)。

当H₁＝H₂＝R₀，(2-10)

所述的椭圆弧段Z_E即成为半径为R₀的圆弧段Z_C。

3、如第2点所述的插补方法，其特点为：

(1)对应所述坐标Ψ₁、Ψ₂的坐标函数Ψ₁(t)、Ψ₂(t)具有相同幅值R₀，函数Ψ₁(t)、Ψ₂(t)的表达式为

${\mathrm{\Ψ}}_1\left(t\right)=H_1\sin (\mathrm{\ωt}+{\mathrm{\α}}_1)=R_0\sin (\mathrm{\ωt}+{\mathrm{\α}}_0+\frac{\mathrm{\π}}{2})=R_0\cos (\mathrm{\ωt}+{\mathrm{\α}}_0),---(3-1)$

      Ψ₂(t)＝H₂sin(ωt+α₂)＝R₀sin(ωt+α₀)，(3-2)

式中，H₁＝H₂＝R₀，(3-3)

${\mathrm{\α}}_1={\mathrm{\α}}_0+\frac{\mathrm{\π}}{2},---(3-4)$

      α₂＝α₀，(3-5)

(2)对应所述坐标Ψ_X及Ψ_Y的替代坐标函数为Ψ_δ1(t)、Ψ_δ2(t)具有相同的幅值R_δ0，函数Ψ_δ1(t)、Ψ_δ2(t)的表达式为

      Ψ_δ1(t)＝H_δ1cos(ωt+α₀)＝R_δ0cos(ωt+α₀)，(3-6)

      Ψ_δ2(t)＝H_δ2sin(ωt+α₀)＝R_δ0sin(ωt+α₀)，(3-7)

式中，H_δ1＝H_δ2＝R_δ0，(3-8)

      H_δ1＝H₁+δ₁＝H_δ2＝H₂+δ₂＝R₀+δ_R，(3-9)

其中，δ₁＝δ₂＝δ_R，(3-10)

(3)对应所述坐标Ψ₁、Ψ₂的替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域等分为相同的n_C个分段，或者说，替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域各分段起点对应相同的参数t的增量Δτ_C或等效增量ΔT_C，

n₁＝n₂＝n_C，(3-11)

Δτ₁＝Δτ₂＝Δτ_C，(3-12)

ΔT₁＝ΔT₂＝ΔT_C。(3-13)

由所述特点可知：

①在按第1点所述的插补方法确定替代坐标函数值增量时，对应于坐标函数Ψ₂(t)的虚拟坐标函数Φ₂(t)就是坐标函数Ψ₁(t)，

Φ₂(t)＝R₀cos(ωt+α₀)＝Ψ₁(t)，(3-10)

虚拟替代坐标函数Φ_δ2(t)就是替代坐标函数Ψ_δ1(t)，

Φ_δ2(t)＝Ψ_δ1(t)。(3-11)

因此，在确定ΔΨ_δ1(t)、ΔΨ_δ2(t)过程中，无需为Ψ_δ1(t)、Ψ_δ2(t)另设虚拟替代坐标函数。

②如果所需路径或轮廓线的位置坐标只包括了二个个坐标Ψ₁和Ψ₂，且坐标轴Ψ₁和Ψ₂构成了直角坐标系Ψ₁OΨ₂，则所需路径或轮廓线即为Ψ₁OΨ₂平面上圆心在坐标轴原点O、半径为R₀的圆弧段Z_C。圆弧段Z_C上的点与圆心O的连线相对轴Ψ₁的夹角即为(ωt+α₀)。每个等分圆弧段所对应的圆心角为

ΔT_C＝ωΔτ_C。(3-12)

曲线在直角坐标系Ψ₁OΨ₂下的方程为

Ψ₁ ²+Ψ₂ ²＝R₀ ²。(3-13)

替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)的图形就是半径为R_δ0的与Z_C同心的替代圆弧段Z_δC。其曲线方程为

Ψ₁ ²+Ψ₂ ²＝R_δ0 ²(3-14)

所述图形可参见图5，Ψ₁、Ψ₂、R₀、R_δ0、α₀、t、Z_C、Z_δC分别对应图中的Φ_A、Ψ_A、H_A、H_δA、α_A、t、C_A、C_δA。图中C_A、C_δA分别是第一象限的圆弧段。

4、如上述第1至3点中任何一点所述的插补方法，其特点为：对应所述坐标其幅值差的取值、其替代坐标函数定义域等分分段起点所对应的参数t的等效增量的取值及替代坐标函数定义域等分分段段数的取值满足下述公式，

0≤δ_A|sin(ωt+α_A)|_MAX≤ε_A，(4-1)

$\left|\mathrm{\Δ}{T}_A\right|\≤\sqrt{8\×\frac{{\mathrm{\ϵ}}_A+{\mathrm{\δ}}_A\×{\left|\sin (\mathrm{\ωt}+{\mathrm{\α}}_A)\right|}_{\mathrm{MAX}}}{H_A\×{\left|\sin (\mathrm{\ωt}+{\mathrm{\α}}_A)\right|}_{\mathrm{MAX}}}},---(4-2)$

$n\≥\mathrm{\ω}|t_{n+1}-t_1|\×\sqrt{\frac{1}{8}\×\frac{H_A\×{\left|\sin (\mathrm{\ωt}+{\mathrm{\α}}_A)\right|}_{\mathrm{MAX}}}{{{\mathrm{\ϵ}}_A+\mathrm{\δ}}_A\×{\left|\sin (\mathrm{\ωt}+{\mathrm{\α}}_A)\right|}_{\mathrm{MAX}}}},---(4-3)$

式中，①Ψ_A表示所述坐标中的某一个坐标，A是序号k(k＝1、2、3、……、m_Ψ)中的某个序号，对应坐标Ψ_A的坐标函数与替代坐标函数为

Ψ_A(t)＝H_Asin(ωt+α_A)，(4-4)

Ψ_δA(t)＝H_δAsin(ωt+α_A)，(4-5)

②δ_A为对应坐标Ψ_A的幅值差，

δ_A＝H_δA-H_A，(4-6)

③ΔT_A为替代坐标函数Ψ_δA(t)定义域等分分段起点所对应的参数t的等效增量，其数值为所述参数t的增量Δτ_A的ω倍，

ΔT_A＝ωΔτ_A，(4-7)

④n表示替代坐标函数Ψ_δA(t)定义域等分分段的段数，

⑤t₁为替代坐标函数Ψ_δA(t)定义域起点所对应的参数t的值，

  t_n+1为替代坐标函数Ψ_δA(t)定义域终点所对应的参数t的值，

⑥|sin(ωt+α_A)|_MAX为对应替代坐标函数Ψ_δA(t)定义域|sin(ωt+α_A)|的最大值，

⑦ε_A为以分段线性函数Ψ_δLA(t)拟合坐标函数Ψ_A(t)的误差绝对值|δΨ_δL-A(t)|的允许值，

  δΨ_δL＝A(t)＝Ψ_δLA(t)-Ψ_A(t)，(4-8)

对应替代坐标函数Ψ_δA(t)的第i(i＝1、2、3、……、n)个分段定义域，所述的分段线性函数Ψ_δLA(t)为

${\mathrm{\Ψ}}_{\mathrm{\δLA}}\left(t\right)={\mathrm{\Ψ}}_{\mathrm{\δA}}\left(t_i\right)+\frac{\mathrm{\Δ}{\mathrm{\Ψ}}_{\mathrm{\δA}}\left(t_i\right)}{\mathrm{\Δ}{T}_A}(t-t_i),,(i=\mathrm{1,2,3},......,n),---(4-9)$

式中，(a)i(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域n个等分分段的序号，等分分段的交点就是定义域的中间点，以i(i＝1、2、3、……、n、n+1)作为包括定义域的二个端点及二个端点间的中间点在内的点的序号，二个端点分别对应着序号1及n+1，二个端点间的中间点分别对应着序号i(i＝2、3、……、n)，

(b)t_i(i＝1、2、3、……、n、n+1)为替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n、n+1)的点所对应的参数t的值，

(c)ΔΨ_δA(t_i)(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n)的点所对应的替代坐标函数值的增量，

ΔΨ_δA(t_i)＝Ψ_δA(t_i+1)-Ψ_δA(t_i )，(i＝1、2、3、……、n)，(4-10)

其中，Ψ_δA(t_i+1)(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域中与序号为i(i＝1、2、3、……、n)的点相邻的后一个点所对应的替代坐标函数值，

Ψ_δA(t_i)(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n)的点所对应的替代坐标函数值。

满足公式(4-1)、(4-2)、(4-3)的δ_A、|ΔT_A|、n取值，将使得所述的拟合误差|δΨ_δL-A(t)|不超过允许值ε_A。

参见图2。如果所需路径或轮廓线的位置坐标只包括了二个个坐标t和Ψ_A，且坐标轴t和Ψ_A构成了直角坐标系tOΨ_A；则所需路径或轮廓线即为平面tOΨ_A上幅值为H_A、初始相位为α_A的正弦曲线段Z_A。其替代曲线段即为幅值为H_δA、初始相位为α_A的正弦曲线段Z_δA。参见图3，分段线形函数Ψ_δLA(t)的图形就是由替代正弦曲线段Z_δA的n个首尾相连的等长内接弦构成的折线Z_δLA。以Ψ_δLA(t)拟合Ψ_A(t)，即相当于以折线Z_δLA拟合Z_A。正弦曲线段Z_δA分段的交点即为需要设定的正弦曲线段Z_δA的中间点，i(i＝1、2、3、……、n、n+1)就是包括正弦曲线段Z_δA的二个端点及二个端点间的中间点在内的点的序号。

根据公式(4-1)、(4-2)、(4-3)确定δ_A、|ΔT_A|、n取值的依据是：

(1)以替代圆弧段C_δA的内接弦拟合圆弧段C_A的径向误差

①坐标函数Ψ_A(t)与对应的虚拟坐标函数Φ_A(t)是一组具有相同幅值H_Λ、相同周期(2π/ω)、相同初始相位α_A及相同定义域的正弦函数与余弦函数，

Ψ_A(t)＝H_Asin(ωt+α_A)，(4-11)

Φ_A(t)＝H_Acos(ωt+α_A)。(4-12)

替代坐标函数Ψ_δA(t)与对应的虚拟替代坐标函数Φ_δA(t)是一组具有相同幅值H_δA、相同周期(2π/ω)、相同初始相位α_A及相同定义域的正弦函数与余弦函数，

Ψ_δA(t)＝H_δAsin(ωt+α_A)，(4-13)

Φ_δA(t)＝H_δAcos(ωt+α_A)，(4-14)

其中，H_δA＝H_A+δ_A，(4-15)

δ_A为替代坐标函数的幅值差。

②参见图5，坐标轴Ψ_A和Φ_A构成了虚拟直角坐标系Φ_AOΨ_A，坐标函数Ψ_A(t)与虚拟坐标函数Φ_A(t)在Φ_AOΨ_A虚拟平面上的图形即为圆心在坐标轴原点O、半径为H_A的圆弧段C_A，圆弧段C_A上的点与圆心O的连线相对轴Φ_A的夹角即为(ωt+α_A)。替代坐标函数Ψ_δA(t)与虚拟替代坐标函数Φ_δA(t)在Φ_AOΨ_A虚拟平面上的图形即为圆心在坐标轴原点O、半径为H_δA的替代圆弧段C_δA，圆弧段C_δA上的点与圆心O的连线相对轴Φ_A的夹角即为(ωt+α_A)。替代坐标函数Ψ_δA(t)定义域的端点与中间点对应着圆弧段C_δA的端点与中间点。图中C_A、C_δA分别是第一象限的圆弧段。

③以C_δA的内接弦拟合C_δA的径向误差

参见图6，将替代圆弧段C_δA分成序号为i(i＝1、2、3、……、n)的n个等分分段，每个分段对应的参数t的等效增量为ΔT_A，ΔT_A也是每个等分圆弧段对应的圆心角，

ΔT_A＝ωΔτ_A，(4-16)

式中，Δτ_A为每个分段对应的参数t的增量。

以圆弧段C_δA的n个等分分段的内接弦构成的折线拟合C_δA。拟合产生的最大径向误差绝对值|δH_δA|发生在内接弦的中点处，且

$\left|{\mathrm{\δH}}_{\mathrm{\δA}}\right|\≈\left|\mathrm{\δ}{H}_A\right|=\frac{1}{8}{H}_A{\left({\mathrm{\ΔT}}_A\right)}^2,---(4-17)$

式中，|δH_A|为以圆弧段C_A的n个等分圆弧段的内接弦构成的折线拟合圆弧段C_A的最大径向误差的绝对值。

④以C_δA的内接弦拟合C_A的误差

以圆弧段C_δA的内接弦构成的折线拟合圆弧段C_A，其最大径向误差绝对值|δH_δ-A|发生在内接弦的中点或内接弦的二个端点处。

(a)在中点处圆弧段C_δA各内接弦拟合圆弧段C_A的径向误差绝对值|δH_δ-A，M|

i)当C_δA内接弦与圆弧段C_A相交，C_δA内接弦就成为C_A的割线，

$\left|{\mathrm{\δH}}_{\mathrm{\δ}-A,M}\right|=\left|{\mathrm{\δH}}_{\mathrm{\δA}}\right|-{\mathrm{\δ}}_A\≈\left|\mathrm{\δ}{H}_A\right|{-\mathrm{\δ}}_A=\frac{1}{8}{H}_A{\left(\mathrm{\Δ}{T}_A\right)}^2-{\mathrm{\δ}}_A,---(4-18)$

此时，

$\left|{\mathrm{\δH}}_{\mathrm{\δ}-A,M}\right|=\frac{1}{8}{H}_A{\left(\mathrm{\Δ}{T}_A\right)}^2-{\mathrm{\δ}}_A>0,---(4-19)$

或 $\left|\mathrm{\Δ}{T}_A\right|>\sqrt{\frac{8{\mathrm{\δ}}_A}{H_A}}.---(4-20)$

ii)当C_δA内接弦与圆弧段C_A相切，

$\left|{\mathrm{\δH}}_{\mathrm{\δ}-A,M}\right|=\left|{\mathrm{\δH}}_{\mathrm{\δA}}\right|-{\mathrm{\δ}}_A=\frac{1}{8}{H}_A{\left(\mathrm{\Δ}{T}_A\right)}^2-{\mathrm{\δ}}_A=0,---(4-21)$

或 $\left|\mathrm{\Δ}{T}_A\right|=\sqrt{\frac{8{\mathrm{\δ}}_A}{H_A}}.---(4-22)$

iii)当内接弦不与圆弧段C_A相交，在圆弧段外侧，

$\left|{\mathrm{\δH}}_{\mathrm{\δ}-A,M}\right|={\mathrm{\δ}}_A-\left|{\mathrm{\δH}}_{\mathrm{\δA}}\right|\≈{\mathrm{\δ}}_A-\left|\mathrm{\δ}{H}_A\right|{=\mathrm{\δ}}_A-\frac{1}{8}{H}_A{\left(\mathrm{\Δ}{T}_A\right)}^2,---(4-23)$

此时 $\left|\mathrm{\&Delta;}{T}_A\right|<\sqrt{\frac{8{\mathrm{\&delta;}}_A}{H_A}},---(4-24)$

且|δ_Hδ-A，M|≤δ_A。(4-25)

(b)在端点处圆弧段C_δA各内接弦拟合圆弧段C_A的径向误差绝对值|δH_δ-A，D|恒为δ_A，

|δH_δ-A，D|＝δ_A。(4-26)

(2)δ_A、|ΔT_A|、n的取值

①δ_A的取值

由于C_δA内接弦端点处径向误差绝对值|δH_δ-A，D|恒为δ_A，为使所述|δH_δ-A，D|不超过允许值ε_Φ其值不能超过允许值ε_Φ，δ_A取值应满足下述公式，

0≤δ_A≤ε_Φ。(4-27)

②|ΔT_A|的取值

(a)在端点处，以C_δA内接弦拟合C_A的径向误差绝对值|δH_δ-A，D|恒为δ_A，不受|ΔT_A|取值的影响。

(b)在中点处，以C_δA内接弦拟合C_A的径向误差绝对值|δH_δ-A，M|应不超过允许值ε_Φ，

|δH_δ-A，M|≤ε_Φ。(4-28)

i)当圆弧段C_δA内接弦与圆弧段C_A相交，

$\left|\mathrm{\&Delta;}{T}_A\right|>\sqrt{\frac{8{\mathrm{\&delta;}}_A}{H_A}}.---(4-20)$

由式(4-18)可知，为使所述的拟合误差|δH_δ-A，M|不超过允许ε_Φ，|ΔT_A|的取值应满足下述公式，

$\left|{\mathrm{\&Delta;T}}_A\right|\&le;\sqrt{\frac{8\&times;({\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}+{\mathrm{\&delta;}}_A)}{H_A}}.---(4-29)$

ii)当圆弧段C_δA内接弦与圆弧段C_A相切，

$\left|\mathrm{\&Delta;}{T}_A\right|=\sqrt{\frac{8{\mathrm{\&delta;}}_A}{H_A}}.---(4-22)$

此时， $\left|{\mathrm{\&delta;H}}_{\mathrm{\&delta;}-A,M}\right|=\left|\mathrm{\&delta;}{H}_A\right|-{\mathrm{\&delta;}}_A=\frac{1}{8}{H}_A{\left({\mathrm{\&Delta;T}}_A\right)}^2-{\mathrm{\&delta;}}_A=0.---(4-21)$

iii)当内接弦不与圆弧段C_A相交，在圆弧段外侧，

$\left|\mathrm{\&Delta;}{T}_A\right|<\sqrt{\frac{8{\mathrm{\&delta;}}_A}{H_A}},---(4-24)$

且|δH_δ-A，M|≤δ_A。(4-25)

由于式(4-27)的限定，拟合误差|δH_δ-A，M|肯定不超过允许值ε_Φ，

|δH_δ-A，M|≤δ_A≤ε_Φ。(4-30)

VI)综合i)、ii)、iii)所述可知，只要ΔT_A的取值满足公式(4-29)，就可满足公式(4-28)的要求，使所述的拟合误差|δH_δ-A，M|不超过允许值ε_Φ。

③综合①、②所述可知，为使以C_δA内接弦拟合C_A的径向误差绝对值|δH_δ-A|不超过允许值ε_Φ，δ_A的取值应满足下述公式，

0≤δ_A≤ε_Φ，(4-31)

ΔT_A的取值应满足下述公式，

$\left|{\mathrm{\&Delta;T}}_A\right|\&le;\sqrt{\frac{8\&times;({\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}+{\mathrm{\&delta;}}_A)}{H_A}},---(4-32)$

或|Δτ_A|的取值应满足下述公式，

$\left|{\mathrm{\&Delta;\&tau;}}_A\right|\&le;\frac{1}{\mathrm{\&omega;}}\&times;\sqrt{\frac{8\&times;({\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}+{\mathrm{\&delta;}}_A)}{H_A}}.---(4-33)$

相应地，圆弧段C_δA的分段段数由于

$n\&GreaterEqual;\frac{|t_{n+1}-t_1|}{\left|{\mathrm{\&Delta;\&tau;}}_A\right|},---(4-34)$

其取值应满足

$n\&GreaterEqual;\mathrm{\&omega;}|t_{n+1}-t_1|\&times;\sqrt{\frac{H_4}{8\&times;({\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}+{\mathrm{\&delta;}}_A)}}.---(4-35)$

④|ΔT_A|的最大允许取值

(a)为加大|ΔT_A|取值，以减少分段数；δ_R的取值应尽量大，取

δ_A＝ε_Φ。(4-36)

相应地，ΔT_A的取值应满足下述公式，

$\left|\mathrm{\&Delta;}{T}_A\right|\&le;4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}}{H_A}},,---(4-37)$

或 $\left|\mathrm{\&Delta;}{\mathrm{\&tau;}}_A\right|\&le;4\&times;\frac{1}{\mathrm{\&omega;}}\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}}{H_A}}.---(4-38)$

(b)|ΔT_A|或|Δτ_A|的最大允许取值即为

$\left|{\mathrm{\&Delta;T}}_A\right|{=\left|{\mathrm{\&Delta;T}}_A\right|}_{\mathrm{MAX}}=4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}}{H_A}},---(4-39)$

或 $\left|{\mathrm{\&Delta;\&tau;}}_A\right|{=\left|{\mathrm{\&Delta;\&tau;}}_A\right|}_{\mathrm{MAX}}=4\&times;\frac{1}{\mathrm{\&omega;}}\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}}{H_A}}.---(4-40)$

而且，|ΔT_A|_MAX对应着|δH_δA|或|δH_A|的最大值|δH_δA|_MZX或|δH_A|_MAX，

${\left|{\mathrm{\&delta;H}}_{\mathrm{\&delta;A}}\right|}_{\mathrm{MZX}}\&ap;{\left|\mathrm{\&delta;}{H}_A\right|}_{\mathrm{MAX}}=\frac{1}{8}{{\left|{\mathrm{\&Delta;T}}_A\right|}_{\mathrm{MAX}}}^2{H}_A=\frac{1}{8}\&times;16\&times;\frac{{\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}}{H_A}\&times;H_A=2{\mathrm{\&epsiv;}}_{\mathrm{\&Phi;}}=2{\mathrm{\&delta;}}_A.---(4-41)$

此时，C_δA内接弦中点与端点所对应的径向误差绝对值相等，

|δH_δ-A，M|＝|δH_A|_MAX-δ_A＝2δ_A-δ_A＝δ_A＝|δH_δ-A，D|。(4-42)

(3)以分段线形函数Ψ_δLA(t)拟合坐标函数Ψ_A(t)的误差δΨ_δL＝A(t)

①以Ψ_δLA(t)拟合Ψ_A(t)的误差记以δΨ_δL＝A(t)，

δΨ_δL＝A(t)＝Ψ_δLA(t)-Ψ_A(t)。(4-43)

②参见图2至图6。在tOΨ_A平面上，将Ψ_δA(t)所对应的正弦曲线段Z_δA分成序号为i(i＝1、2、3、……、n)的n个等分分段；在虚拟Φ_AOΨ_A平面上，将Ψ_δA(t)与Φ_δA(t)所对应的圆弧段C_δA分成相同的n个等分分段。序号相同的Z_δA分段端点与C_δA分段端点所对应的t值相同，都是t_i(i＝1、2、3、……、n、n+1)。Z_δA分段与C_δA分段对应的参数t的增量或等效增量相同。以Z_δA的内接弦构成的折线拟合Z_A，其拟合误差为δΨ_δL＝A(t)；以C_δA的内接弦构成的折线拟合C_A，其对应坐标轴Ψ_A拟合误差等于内接弦径向误差在坐标轴上的投影值δH_δ＝A，Ψ(t)。图6中，以t_M，u表示一个分段中点处的参数t的值，对应t_M，u的所述径向误差绝对值在坐标轴上的投影值即为|δH_δ-A，Ψ(t_M，u)|。

③以C_δA的内接弦构成的折线拟合C_A，其径向误差绝对值的最大值|δH_δ-A|发生在C_δA内接弦的中点或端点t_M/D处，且|δH_δ-A|在坐标轴Ψ_A上的投影为

|δH_δ-A，Ψ(t_M/D)|＝|δH_δ-A|×|sin(ωt_M/D+α_A)|。(4-44)

④相应地认为，一个分段内在中点或端点t_M/D处，以正弦曲线段Z_A的内接弦拟合正弦曲线段Z_δA的误差绝对值|δΨ_δ-LA(t_M./D)|也是最大；且等于对应t_M/D处|δH_δ-A|在坐标轴Ψ_A上的投影，

|δΨ_δ-LA(t_M./D)|＝|δH_δ-A|×|sin(ωt_M/D+α_A)|。(4-45)

⑤在定义域范围内，|δΨ_δ-LA(t_M./D)|的最大值|δΨ_δ-LA(t_M./D)|_MAX为|δΨ_δ-LA(t_M./D)|_MAX＝|δH_δ-A|×|sin(ωt_M/D+α_A)|_MAX，(4-46)

式中，|sin(ωt_M/D+α_A)|_MAX为对应替代坐标函数Ψ_δA(t)定义域|sin(ωt_M/D+α_A)|的最大值。

由于|ΔT_A|或说|Δτ_A|很小，可以认为

|sin(ωt_M/D+α_A)|_MAX＝|sin(ωt+α_A)|_MAX，(4-47)

式中，|sin(ωt+α_A)|_MAX为对应替代坐标函数Ψ_δA(t)定义域|sin(ωt+α_A)|的最大值，式(4-46)可以改写为

|δΨ_δ-LA(t_M./D)|_MAX＝|δH_δ-A|×|sin(ωt+α_A)|_MAX。(4-48)

由于所述的|δΨ_δ-LA(t_M./D)|的最大值|δΨ_δ-LA(t_M./D)|_MAX也就是所述的拟合误差δΨ_δL＝A(t)绝对值的最大值|δΨ_δ-LA(t)|_MAX，因此，

|δΨ_δ-LA(t)|_MAX＝|δH_δ-A|×|sin(ωt+α_A)|_MAX。(4-49)

(4)δ_A、|ΔT_A|、n的取值

①δ_A的取值

C_δA内接弦端点处径向误差绝对值|δH_δ-A，D|恒为δ_A，相应地，端点处以Ψ_δLA(t)拟合Ψ_A(t)的误差绝对值的最大值|δΨ_δ-LA(t)|_MAX为δ_A|sin(t+α_A)|_MAX，其值不能超过允许值ε_A；因此，δ_A取值应满足

  0≤δ_A|sin(ωt+α_A)|_MAX≤ε_A，(4-50)

或0≤δ_A≤ε_A/|sin(ωt+α_A)|_MAX。(4-51)

②|ΔT_A|的取值

(a)在端点处，以Ψ_δLA(t)拟合Ψ_A(t)的误差绝对值|δΨ_δ-LA(t)|恒为δ_A|sin(ωt+α_A)|_MAX，不受|ΔT_A|取值的影响，只要δ_A取值满足式(4-50)，|δΨ_δ-LA(t)|就不会超过允许值ε_A。

(b)在中点处以Ψ_δLA(t)拟合Ψ_A(t)的误差绝对值的最大值|δΨ_δ-LA(t)|_MAX应不超过允许值ε_A，

|δΨ_δ-LA(t)|_MAX＝|δH_δ-A，M||sin(ωt+α_A)|_MAX≤ε_A，(4-52)

i)当 $\left|\mathrm{\&Delta;}{T}_A\right|>\sqrt{\frac{8{\mathrm{\&delta;}}_A}{H_A}},---(4-20)$

C_δA内接弦与圆弧段C_A相交，C_δA内接弦就成为C_A的割线；此时，为满足式(4-52)的要求，由式(4-18)、(4-52)可知，应有

$\left|\mathrm{\&Delta;}{T}_A\right|\&le;\sqrt{8\&times;\frac{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}},---(4-53)$

或 $\left|\mathrm{\&Delta;}{\mathrm{\&tau;}}_A\right|\&le;\frac{1}{\mathrm{\&omega;}}\&times;\sqrt{8\&times;\frac{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}}.---(4-54)$

ii)当 $\left|\mathrm{\&Delta;}{T}_A\right|=\sqrt{\frac{8{\mathrm{\&delta;}}_A}{H_A}},---(4-22)$

C_δA内接弦与圆弧段C_A相切，中点处|δH_δ-A，M|将为0，

$\left|{\mathrm{\&delta;H}}_{\mathrm{\&delta;}-A,M}\right|=\left|\mathrm{\&delta;}{H}_A\right|-{\mathrm{\&delta;}}_A=\frac{1}{8}{H}_A{\left({\mathrm{\&Delta;T}}_A\right)}^2-{\mathrm{\&delta;}}_A=0.---(4-21)$

因而，相应的Ψ_δLA(t)拟合Ψ_A(t)的误差|δΨ_δ-LA(t)|为0，

|δΨ_δ-LA(t)|＝|δH_δ-A，M||sin(ωt+α_A)|_MAX＝0＜ε_A。(4-55)

iii)当 $\left|\mathrm{\&Delta;}{T}_A\right|>\sqrt{\frac{8{\mathrm{\&delta;}}_A}{H_A}},---(4-24)$

C_δA内接弦不与圆弧段C_A相交，在圆弧段C_A外侧，此时总有

|δH_δ-A，M|＜δ_A，(4-25)

只要δ_A取值满足式(4-50)，|δΨ_δ-LA(t)|_MAX都不会超过允许值ε_A，

|δΨ_δ-LA(t)|_MAX＝|δH_δ-A，M||sin(ωt+α_A)|_MAX＜δ_A|sin(ωt+α_A)|_MAX≤ε_A。(4-56)

VI)综合i)、ii)、iii)所述可知，只要ΔT_A的取值满足公式(4-53)，就可满足公式(4-52)的要求，使所述的拟合误差|δΨ_δ-LA(t)|_MAX不超过允许值ε_A。

③综合①、②所述可知，为使所述的拟合误差|δΨδ_L-A(t)|不超过允许值ε_A，δ_A的取值应满足下述公式，

  0≤δ_A|sin(ωt+α_A)|_MAX≤ε_A， (4-57)

或0≤δ_A≤ε_A/|sin(ωt+α_A)|_MAX。(4-58)

ΔT_A的取值应满足下述公式，

$\left|\mathrm{\&Delta;}{T}_A\right|\&le;\sqrt{8\&times;\frac{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}},---(4-59)$

或|Δτ_A|的取值应满足下述公式

$\left|\mathrm{\&Delta;}{\mathrm{\&tau;}}_A\right|\&le;\frac{1}{\mathrm{\&omega;}}\&times;\sqrt{8\&times;\frac{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}}.---(4-60)$

相应地，替代坐标函数定义域的等分分段段数n的取值应满足

$n\&GreaterEqual;\mathrm{\&omega;}|t_{n+1}-t_1|\&times;\sqrt{\frac{1}{8}\&times;\frac{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}{{{\mathrm{\&epsiv;}}_A+\mathrm{\&delta;}}_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}.}---(4-61)$

④|ΔT_A|的最大允许取值值

(a)为加大|ΔT_A|取值，以减少分段数；δ_A的取值尽量大，取

δ_A|sin(ωt+α_A)|_MAX＝ε_A。(4-62)

相应地，ΔT_A的取值应满足下述公式，

$\left|\mathrm{\&Delta;}{T}_A\right|\&le;4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}},---(4-63)$

或 $\left|\mathrm{\&Delta;}{\mathrm{\&tau;}}_A\right|\&le;\frac{1}{\mathrm{\&omega;}}\&times;4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}}.---(4-64)$

(b)|ΔT_A|或|Δτ_A|的最大允许取值即为

$\left|\mathrm{\&Delta;}{T}_A\right|={\left|\mathrm{\&Delta;}{T}_A\right|}_{\mathrm{MAX}}=4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}},---(4-65)$

或 $\left|\mathrm{\&Delta;}{\mathrm{\&tau;}}_A\right|={\left|\mathrm{\&Delta;}{\mathrm{\&tau;}}_A\right|}_{\mathrm{MAX}}=4\&times;\frac{1}{\mathrm{\&omega;}}\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}}.--(4-66)$

而且，|ΔT_A|_MAX对应着|δH_δA|或|δH_A|的最大值

${\left|{\mathrm{\&delta;H}}_A\right|}_{\mathrm{MAX}}=\frac{1}{8}{{\left|\mathrm{\&Delta;}{T}_A\right|}_{\mathrm{MAX}}}^2{H}_A=\frac{1}{8}\&times;16\&times;\frac{{\mathrm{\&epsiv;}}_A}{H_A}\&times;H_A=2{\mathrm{\&epsiv;}}_A/{\left|\sin (t+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}=2{\mathrm{\&delta;}}_A,---(4-67)$

此时，C_δA内接弦中点与端点所对应的径向误差绝对值相等，

|δH_δ-A，M|＝|δH_A|_MAX-δ_A＝2δ_A-δ_A＝δ_A＝|δH_δ-A，D|。(4-68)

5、如第4点所述的插补方法，其特点为：对应坐标Ψ_A其幅值差δ_A取值为

δ_A＝ε_A/|sin(ωt+α_A)|_MAX。(5-1)

相应地，对应所述坐标Ψ_A其替代坐标函数定义域等分分段起点所对应的参数t的等效增量ΔT_A的取值及其替代坐标函数定义域等分分段段数n的取值满足下述公式，

$\left|\mathrm{\&Delta;}{T}_A\right|\&le;4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}},---(5-2)$

$n\&GreaterEqual;\frac{\mathrm{\&omega;}|t_{n+1}-t_1|}{4}\&times;\sqrt{\frac{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&epsiv;}}_A)\right|}_{\mathrm{MAX}}}{{\mathrm{\&epsiv;}}_A}}.---(5-3)$

|ΔT_A|的最大允许取值即为

$\left|\mathrm{\&Delta;}{T}_A\right|={\left|\mathrm{\&Delta;}{T}_A\right|}_{\mathrm{MAX}}=4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}}.---(5-4)$

其依据见第4点中对确定δ_A、|ΔT_A|、n取值的依据的分析。

6、如第4点所述的一种插补方法，其特点为：所述|sin(ωt+α_A)|_MAX以1替换。

相应地，所述δ_A、ΔT_A及n的取值满足下述公式，

0≤δ_A≤ε_A，(6-1)

$\left|\mathrm{\&Delta;}{T}_A\right|\&le;\sqrt{8\&times;\frac{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A}{H_A}},---(6-2)$

$n\&GreaterEqual;\mathrm{\&omega;}|t_{n+1}-t_1|\&times;\sqrt{\frac{1}{8}\&times;\frac{H_A}{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A}}.---(6-3)$

满足公式(6-1)、(6-2)、(6-3)的δ_A、|ΔT_A|、n的取值，可以使得对任何可能的定义域，都可做到以分段线性函数Ψ_δLA(t)拟合坐标函数Ψ_A(t)的误差绝对值|δΨ_δL-A(t)|不超过允许值ε_A。

7、如第6点所述的插补方法，其特点为：对应坐标Ψ_A其幅值差δ_A取值为ε_A，

δ_A＝ε_A。(7-1)

相应地，对应所述坐标Ψ_A其替代坐标函数定义域等分分段起点所对应的参数t的等效增量ΔT_A的取值及其替代坐标函数定义域等分分段段数n的取值满足下述公式，

$\left|\mathrm{\&Delta;}{T}_A\right|\&le;4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A}},---(7-2)$

$n\&GreaterEqual;\frac{\mathrm{\&omega;}|t_{n+1}-t_1|}{4}\&times;\sqrt{\frac{H_A}{{\mathrm{\&epsiv;}}_A}}.---(7-3)$

|ΔT_A|的最大允许取值即为

$\left|{\mathrm{\&Delta;T}}_A\right|={{|\mathrm{\&Delta;T}}_A|}_{\mathrm{MAX}}=4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_A}{H_A}}.---(7-4)$

依据式(7-4)确定的|ΔT_A|_MAX是对任何可能的的定义域所述拟合误差|δΨ_δL-A(t)|都不超过允许值ε_A的|ΔT_A|的最大允许取值

8、如第3点所述的插补方法，其特点为：对应所述坐标Ψ₁、Ψ₂其幅值差δ_R的取值、其替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域等分分段起点所对应的参数t的等效增量|ΔT_C|的取值及其替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域等分分段段数n_C的取值满足下述公式，

0＜δ_R≤ε_C，(8-1)

$\left|{\mathrm{\&Delta;T}}_C\right|\&le;\sqrt{\frac{8\&times;({\mathrm{\&epsiv;}}_C+{\mathrm{\&delta;}}_R)}{R_0}},---(8-2)$

$n_C\&GreaterEqual;\frac{\mathrm{\&omega;}|t(n_C+1)-t_1|}{4}\&times;\sqrt{\frac{R_0}{8\&times;({\mathrm{\&epsiv;}}_C+{\mathrm{\&delta;}}_R)}},---(8-3)$

式中，①t₁为替代坐标函数Ψ_δ1(t)或Ψ_δ2(t)定义域起点对应的参数t的值，

t(n_C+1)为替代坐标函数Ψ_δ1(t)或Ψ_δ2(t)定义域终点对应的参数t的值，

②ε_C为以圆弧段Z_δC的等分圆弧分段内接弦拟合圆弧段Z_C的径向误差的绝对值的允许值，圆弧段Z_δC为替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)在坐标轴Ψ₁和Ψ₂构成的直角坐标系Ψ₁OΨ₂的坐标平面上的图形，圆弧段Z_C为坐标函数Ψ₁(t)、Ψ₂(t)在坐标平面Ψ₁OΨ₂上的图形。所述圆弧段Z_C、Z_δC可参见图5，Ψ₁、Ψ₂、R₀、R_δ0、α₀、t、Z_C、Z_δC分别是图中的Φ_A、Ψ_A、H_A、H_δA、α_A、t、C_A、C_δA。

满足公式(8-1)、(8-2)、(8-3)的δ_R、|ΔT_C|或n_C取值，将使得以圆弧段Z_δC的等分圆弧分段内接弦拟合圆弧段Z_C的径向误差的绝对值不超过允许值ε_C。

公式(8-1)、(8-2)、(8-3)的依据见第4点对δ_A、|ΔT_A|、n取值依据的分析。

9、如第8点所述的插补方法，其特点为：所述δ_R取值为ε_C，

δ_R＝ε_C。(9-1)

相应地，所述|ΔT_C|及n_C的取值满足下述公式，

$\left|\mathrm{\&Delta;}{T}_C\right|\&le;4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_C}{R_0}},---(9-2)$

$n_C\&GreaterEqual;\frac{\mathrm{\&omega;}|t(n_C+1)-t_1|}{4}\&times;\sqrt{\frac{R_0}{{\mathrm{\&epsiv;}}_C}}.---(9-3)$

|ΔT_A|的最大允许取值即为

$\left|{\mathrm{\&Delta;T}}_A\right|={|\mathrm{\&Delta;T}}_A{|}_{\mathrm{MAX}}=4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_C}{R_0}}.---(9-4)$

其依据见第4点中对δ_A、|ΔT_A|、n取值依据的分析。

10、如第1至3点或第5至9点中的任何一点所述的插补方法，其特点为：

ω＝1。(10-1)

也就是说，所述坐标对应的正弦坐标函数的周期为2π。此时，所述坐标对应的替代坐标函数定义域各等分分段起点所对应的参数t的等效增量，等于所述参数t的增量；以A表示所述坐标其序号k(k＝1、2、3、……、m_Ψ)中的某个序号，则有

ΔT_A＝Δτ_A。(10-2)

11、如第4点所述的插补方法，其特点为：

ω＝1。(11-1)

也就是说，所述坐标Ψ_A对应的正弦坐标函数Ψ_A(t)的周期为2π。此时，所述替代坐标函数Ψ_δA(t)定义域各等分分段起点所对应的参数t的等效增量ΔT_A，等于所述参数t的增量Δτ_A，

ΔT_A＝Δτ_A。(11-2)

值得注意的是，如果一个坐标函数是由若干个函数之和构成，则该坐标函数的插补结果等于构成该坐标函数的各个函数各自的插补结果之和。特别是构成该坐标函数的某个函数如果是常量，则该常量不会影响该函数的增量值。此外，所需路径或轮廓线其所在的直角坐标系坐标轴平移时，也不影响其位置坐标值增量的数值。

本发明的有益效果为：本发明针对正弦函数类所需路径或轮廓线，包括正弦曲线、椭圆曲线、圆弧曲线，提出一种通过算术运算的递推公式准确确定曲线中间点位置坐标值增量的插补方法，并且还提出一种设定替代曲线的方法，以提高拟合精度或减小插补运算量。由于运算简单，对于同样的分段段数，完成运算所需时间少，相当于提高了插补运算速度；或者，在同等的运算时间内，允许增加曲线分段段数，从而提高折线拟合曲线的精度。由于运算简单，这种方法可在单片计算机上实现，从而降低了插补器的成本。因此，这种插补方法速度快，精度高，装置成本低。

【附图说明】

图1是本发明的坐标函数Ψ_Λ(t)插补的程序流程图，

图2是本发明的正弦曲线示意图，

图3是本发明的分段线性函数示意图，

图4是本发明的正弦曲线段局部示意图，

图5是本发明的圆弧曲线示意图，

图6是本发明的圆弧曲线局部示意图，

图7是本发明的椭圆曲线示意图。

【具体实施方式】

一、参见图1，这是坐标函数Ψ_Λ(t)插补的程序流程图。所述的坐标函数为

Ψ_Λ(t)＝H_Λsin(ωt+α_Λ)，(J1-1)

程序由步骤(1)至步骤(9)组成：

(1)设定Ψ_Λ(t)的替代坐标函数Ψ_δΛ(t)

Ψ_δΛ(t)＝H_δΛsin(ωt+α_Λ)，(J1-2)

H_δΛ＝H_Λ+δ_Λ，(J1-3)

确定替代坐标函数Ψ_δΛ(t)的幅值差δ_Λ的数值；

(2)确定替代坐标函数Ψ_δΛ(t)定义域起点及中间点所对应的参数t的等效增量ΔT_Λ，其数值为所述参数t的增量Δτ_Λ的ω倍，是常数，

ΔT_Λ＝ωΔτ_Λ＝常数，(J1-4)

(3)确定替代坐标函数Ψ_δΛ(t)定义域中间点的个数n_Λ-1，以i_Λ(i_Λ＝1、2、3、……、n_Λ、n_Λ+1)表示包括定义域端点与中间点在内的点的序号，中间点的序号为i_Λ(i_Λ＝2、3、……、n_Λ)；

(4)初始化i_Λ值，即将i_Λ设置为2，

2→i_Λ；(J1-5)

(5)确定序号为i_Λ的中间点所对应的替代坐标函数函数值增量ΔΨ_δΛ(t)的值；

(6)存储/输出运算结果；

(7)判定

i_Λ＝n_Λ？(J1-6)

若i_Λ≠n_Λ，表明插补尚未完成，则转至步骤(8)，

若i_Λ＝n_Λ，则转至步骤(9)；

(8)将i_Λ加1，并保存，

i_Λ+1→i_Λ，(J1-7)

再转至步骤(5)，继续进行插补；

(9)插补完成。

二、参见图2，这是正弦曲线示意图。在直角坐标系tOΨ_A下的一个幅值为H_A、初始相位为α_A的正弦曲线Q_A，其坐标函数为

Ψ_A(t)＝H_Asin(ωt+α_A)，(J2-1)

曲线Q_A的替代曲线为幅值为H_δA、初始相位为α_A的正弦曲线Q_δA，相应地，替代坐标函数为

Ψ_δA(t)＝H_δAsin(ωt+α_A)，(J2-2)

式中，H_δA＝H_A+δ_A，(J2-3)

δ_A是替代坐标函数Ψ_δΛ(t)的幅值差。

对应于所需路径或轮廓线，坐标函数Ψ_A(t)、Ψ_δA(t)定义域为[t₁，t_n+1]，该定义域对应的正弦曲线段记以Z_A，对应的替代正弦曲线段记以Z_δA。

三、参见图3，这是等分分段线性函数示意图。在直角坐标系tOΨ_A下的替代正弦曲线段Z_δA，其替代坐标函数为

Ψ_δA(t)＝H_δAsin(ωt+α_A)，(J3-1)

n是替代坐标函数Ψ_δA(t)定义域的等分分段数。正弦曲线段Z_δA的n个等分分段的内接弦构成折线Z_δLA，其相应的分段线性函数Ψ_δLA(t)为

${\mathrm{\&Psi;}}_{\mathrm{\&delta;LA}}\left(t\right)={\mathrm{\&Psi;}}_{\mathrm{\&delta;A}}\left(t_i\right)+\frac{\mathrm{\&Delta;}{\mathrm{\&Psi;}}_{\mathrm{\&delta;A}}\left(t_i\right)}{\mathrm{\&Delta;}{T}_A}(t-t_i),,(t_i<t\&le;t_{i+1}),$

(i＝1、2、3、……、n)，(J3-2)

式中，①i(i＝1、2、3、……、n)是替代坐标函数Ψ_δA(t)定义域n个等分分段的序号，等分分段的交点就是定义域的中间点，以i(i＝1、2、3、……、n、n+1)作为包括定义域的二个端点及二个端点间的中间点在内的点的序号，二个端点分别对应着序号1及n+1，二个端点间的中间点分别对应着序号i(i＝2、3、……、n)，

②t_i(i＝1、2、3、……、n、n+1)是替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n、n+1)的点所对应的参数t的值，

t₁是替代坐标函数Ψ_δA(t)定义域起点所对应的参数t的值，

t_n+1是替代坐标函数Ψ_δA(t)定义域终点所对应的参数t的值，

[t₁、t_u+1]是替代坐标函数Ψ_δA(t)的定义域，

③ΔT_A是替代坐标函数Ψ_δA(t)定义域各等分分段起点所对应的参数t的等效增量的值，或说是替代坐标函数Ψ_δA(t)定义域各等分分段所对应的参数t的等效增量，是所述参数t的增量Δτ_A的ω倍，是常数，

ΔT_A＝ωΔτ_A，(J3-3)

④ΔΨ_δA(t_i)(i＝1、2、3、……、n)是替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n)的点所对应的替代坐标函数值的增量，

ΔΨ_δA(t_i)＝Ψ_δA(t_i+1)-Ψ_δA(t_i)，(i＝1、2、3、……、n)，(J3-5)

其中，Ψ_δA(t_i+1)(i＝1、2、3、……、n)是替代坐标函数Ψ_δA(t)定义域中与序号为i(i＝1、2、3、……、n)的点相邻的后一个点所对应的替代坐标函数值，

Ψ_δA(t_i)(i＝1、2、3、……、n)是替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n)的点所对应的替代坐标函数值。

四、参见图4，这是正弦曲线段局部示意图。示意图表示相应于定义域[t₁，t_n+1]的一个分段，正弦曲线段Z_A、替代正弦曲线段Z_δA及由替代正弦曲线段Z_δA内接弦构成的折线Z_δLA之间的关系。图中以u表示所示分段的序号，该分段对应的参数t的区间为[t_u，t_u+1]，[t_u，t_u+1]的中点以t_M，u表示。Ψ_A(t_M，u)、Ψ_δA(t_M，u)、Ψ_δLA(t_M，u)分别是对应于t_M，u函数Ψ_A(t)、Ψ_δA(t)、Ψ_δLA(t)的值，δΨ_δ-LA(t_M，u)是对应于点t_M，u以Ψ_δLA(t)拟合Ψ_A(t)的误差δΨ_δ-LA(t)，

δΨ_δ-LA(t_M，u)＝Ψ_δLA(t_M，u)-Ψ_A(t_M，u)。(J4-1)

五、参见图5，这是圆弧曲线示意图。在直角坐标系Φ_AOΨ_A下的一个圆心在坐标轴原点O、半径为H_A的圆弧曲线Q_C，其曲线方程为

Φ_A ²+Ψ_A ²＝H_A ²。(J5-1)

曲线的坐标函数为

Ψ_A(t)＝H_Asin(ωt+α_A)，(J5-2)

Φ_A(t)＝H_Acos(ωt+α_A)。(J5-3)

圆弧Q_C上的点与圆心O的连线相对轴Φ_A的夹角即为(ωt+α_A)。

替代圆弧曲线Q_Cδ为圆心在坐标轴原点O、半径为H_δA的圆弧，其曲线方程为

Φ_A ²+Ψ_A ²＝H_δΛ ²。(J5-4)

曲线Q_Cδ的坐标函数为

Ψ_δA(t)＝H_δAsin(ωt+α_A)，(J5-5)

Φ_δA(t)＝H_δAcos(ωt+α_A)，(J5-6)

式中H_δA＝H_A+δ_A，(J5-7)

δ_A是Ψ_δA(t)、Φ_δA(t)的幅值差。

圆弧Q_Cδ上的点与圆心O的连线相对轴Φ_A的夹角即为(ωt+α_A)。

坐标函数Ψ_A(t)、Φ_A(t)、Ψ_δA(t)、Φ_δA(t)定义域为[t₁，t_n+1]，该定义域对应的圆弧曲线段记以C_A，对应的替代圆弧曲线段记以C_δA；图中C_A、C_δA分别是第一象限的圆弧段。

六、参见图6，这是圆弧曲线局部示意图。示意图表示相应于圆弧曲线段C_A在第一象限的一个分段其圆弧段C_A、替代圆弧段C_δA及由替代圆弧曲线段C_δA内接弦之间的关系。

①图中圆弧段C_δA的内接弦就是圆弧段C_A的割线。以圆弧段C_δA内接弦构成的折线拟合圆弧段C_A，实际上就是由圆弧段C_A的割线构成的折线拟合圆弧段C_A，从而减小拟合的误差。

②以t_M，u表示一个分段中点处的参数t的值。对应t_M，u，所述C_δA内接弦拟合C_A的径向误差绝对值在坐标轴上的投影值即为|δH_δ-A，Ψ(t_M，u)|。

③圆弧段C_δA的等分分段对应的圆心角为ΔT_A。以l_δA表示圆弧段C_δA一个分段的内接弦的长度，对应于t_u有

$\mathrm{\&Delta;}{\mathrm{\&Psi;}}_{\mathrm{\&delta;A}}\left(t_u\right)=l_{\mathrm{\&delta;A}}\cos [\mathrm{\&omega;}(t_u+\frac{\mathrm{\&Delta;}{\mathrm{\&tau;}}_A}{2})+{\mathrm{\&alpha;}}_A]=l_{\mathrm{\&delta;A}}\cos ({\mathrm{\&omega;t}}_u+\frac{\mathrm{\&Delta;}{T}_A}{2}+{\mathrm{\&alpha;}}_A),---(J6-1)$

$\mathrm{\&Delta;}{\mathrm{\&Phi;}}_{\mathrm{\&delta;A}}\left(t_u\right)={-l}_{\mathrm{\&delta;A}}\sin [\mathrm{\&omega;}(t_u+\frac{\mathrm{\&Delta;}{\mathrm{\&tau;}}_A}{2})+{\mathrm{\&alpha;}}_A]={-l}_{\mathrm{\&delta;A}}\sin ({\mathrm{\&omega;t}}_u+\frac{\mathrm{\&Delta;}{T}_A}{2}+{\mathrm{\&alpha;}}_A),---(J6-2)$

式中，ΔT_A＝ωΔτ_A。(J6-3)

七、参见图7，这是椭圆曲线示意图。在直角坐标系Ψ₁OΨ₂下的中心在坐标轴原点O的椭圆曲线Q_E，其曲线方程为

$\frac{{{\mathrm{\&Psi;}}_1}^2}{{H_1}^2}+\frac{{{\mathrm{\&Psi;}}_2}^2}{{H_2}^2}=1,---(J7-1)$

式中：H₁和H₂为椭圆的半轴。

曲线的坐标函数为

${\mathrm{\&Psi;}}_1\left(t\right)=H_1\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_1)=H_1\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_0+\frac{\mathrm{\&pi;}}{2})=H_1\cos (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_0),---(J7-2)$

Ψ₂(t)＝H₂sin(ωt+α₂)＝H₂sin(ωt+α₀)，(J7-3)

式中， ${\mathrm{\&alpha;}}_1={\mathrm{\&alpha;}}_0+\frac{\mathrm{\&pi;}}{2},---(J7-4)$

α₂＝α₀，(J7-5)

椭圆曲线Q_E上的点与坐标轴原点O的连线相对Ψ₁轴的夹角即为(ωt+α₀)。

坐标函数Ψ₁(t)、Ψ₂(t)定义域所对应的椭圆曲线段记以Z_E；图中Z_E是一个第一象限的椭圆弧段。

八、实施例

所需路径或轮廓线是在直角坐标平面Ψ_XOΨ_Y上的圆心在坐标轴原点O、半径为R₀的第1象限圆弧段Z_C，其曲线方程为

Ψ_X ²+Ψ_Y ²＝R₀ ²。(L-1)

曲线的坐标函数为

Ψ_X(t)＝R₀cost，(L-2)

Ψ_Y(t)＝R₀sint。(L-3)

R₀＝50000，

圆弧段Z_C上的点与圆心O的连线相对轴Ψ_X的夹角即为参数t。要求以折线拟合圆弧段Z_C的最大径向误差绝对值不超过ε_C，

ε_C＝0.5。(L-4)

参见图5，Ψ_X、Ψ_Y、R₀、t、Z_C分别是图中的Φ_A、Ψ_A、H_A、ωt、C_A。

以替代圆弧段Z_δC的内接弦构成的折线拟合圆弧段Z_C，其插补过程如下：

1、设定替代曲线确定幅值差

(1)设定替代曲线

以与圆弧段Z_C同心的第1象限圆弧段Z_δC设定为圆弧段Z_C的替代曲线，以替代圆弧段Z_δC的内接弦拟合圆弧段Z_C。Z_δC的曲线方程为

Ψ_X ²+Ψ_Y ²＝R_δ0 ²。(L-5)

曲线的坐标函数为

Ψ_δX(t)＝R_δ0cost，(L-6)

Ψ_δY(t)＝R_δ0sint，(L-7)

R_δ0＝R₀+δ_R。(L-8)

(2)确定Ψ_δX(t)、Ψ_δY(t)的幅值差δ_R

根据公式(9-1)，取

δ_R＝ε_C＝0.5，(L-9)

则R_δ0＝R₀+δ_R＝50000.5。(L-10)

参见图5，Ψ_X、Ψ_Y、R_δ0、t、Z_δC分别是图中的Φ_A、Ψ_A、H_δA、t、C_δA。

2、确定内接弦对应的圆心角ΔT_C与圆弧段Z_δC的分段

(1)确定ΔT_C

将圆弧段Z_δC等分，以ΔT_C表示圆弧段Z_δC各等分分段起点所对应的参数t的等效增量，ΔT_C也就是圆弧段Z_δC内接弦对应的圆心角。对于本例，ΔT_C等于所述参数t的增量Δτ_C，

ΔT_C＝Δτ_C。(L-11)

如果要求Z_δC内接弦拟合Z_C在内接弦中点处的径向误差绝对值的最大值|δR_δ-0，M|不超过0.5，则根据公式(9-2)，

$\left|\mathrm{\&Delta;}{T}_C\right|\&le;4\&times;\sqrt{\frac{{\mathrm{\&epsiv;}}_C}{R_0}}=4\&times;\sqrt{\frac{0.5}{5000}}=0.012649,---(L-12)$

取|ΔT_C|＝0.012649，(L-13)

(2)确定圆弧段Z_δC分段段数

依据式(L-13)确定的ΔT_C值，对圆弧段Z_δC进行分段，其分段段数理论值应为

$n_{\mathrm{CL}}=\frac{\mathrm{\&pi;}}{2\&times;0.0126}=124.18,---(L-14)$

实际分段段数n_C只能是整数，根据公式(9-3)，取

n_C＝125，(L-15)

此时，|ΔT_C|实际值为

$\left|\mathrm{\&Delta;}{T}_C\right|=\frac{\mathrm{\&pi;}}{2\&times;n_C}=\frac{\mathrm{\&pi;}}{2\&times;125}=0.01257,---(L-16)$

根据公式(4-18)，Z_δC内接弦拟合Z_C中点处径向误差为

|δR_δ-0，M|＝|δR_δ0|-δ_R＝0.487＜ε_C＝0.5，(L-17)

根据公式(4-26)，Z_δC内接弦拟合Z_C端点处的径向误差绝对值为

|δR_δ-0，D|＝δ_R＝ε_C＝0.5。(L-18)

3、确定替代坐标函数定义域中间点的位置坐标值

以i_C(i_C＝1、2、3、……、125、126)表示构成逼近圆弧Z_C的折线的125个内接弦的起点及最后一个内接弦的终点在内的各点的序号，也就是坐标函数Ψ_δX(t)、Ψ_δY(t)定义域的二个端点及二个端点间的中间点在内的点的序号，二个端点分别对应着序号1及126，二个端点间的中间点分别对应着序号

i_C＝2、3、……、125。(L-19)

(1)当i_C＝1，即圆弧段Z_C起点，其位置坐标值为已知

Ψ_δX(t₁)＝R_δ0＝50000.5，(L-20)

Ψ_δY(t₁)＝0，(L-21)

式中t₁＝0。

(2)当i_C＝2～175，以v+1表示i_C(i_C＝2、3、……、n_C-1)中的某个中间点的序号，v代表与之相邻的前一个点的序号.，根据第3点的说明可知，Ψ_X(t)可视为Ψ_Y(t)的虚拟坐标函数，因此，可根据公式(1-6)、(1-7)确定圆弧段各中间点的坐标函数值增量，

ΔΨ_δX(t_v+1)＝ΔΨ_δX(t_v)cosΔT_C-ΔΨ_δY(t_v)sinΔT_C，(L-22)

ΔΨ_δY(t_v+1)＝ΔΨ_δY(t_v)cosΔT_C+ΔΨ_δX(t_v)sinΔT_C。(L-23)

相应的各中间点的位置坐标值Ψ_δX(t_v+1)、Ψ_δY(t_v+1)由下式确定

Ψ_δX(t_v+1)＝Ψ_δX(t_v)+ΔΨ_δX(t_v)，(L-24)

Ψ_δY(t_v+1)＝Ψ_δY(t_v)+ΔΨ_δY(t_v)，(L-25)

式(L-22)、(L-23)中，

①根据公式(1-24)、(1-25)，sinΔT_C、cosΔT_C可按下式确定，

$\sin \mathrm{\&Delta;}{T}_C=\mathrm{\&Delta;}{T}_C-\frac{1}{6}\&times;{\left(\mathrm{\&Delta;}{T}_C\right)}^3,---(L-26)$

$\cos \mathrm{\&Delta;}{T}_C=1-\frac{1}{2}{\left(\mathrm{\&Delta;}{T}_C\right)}^2+\frac{1}{24}{\left(\mathrm{\&Delta;}{T}_C\right)}^4.---(L-27)$

②根据公式(1-34)、(1-35)，对应圆弧段Z_δC起点的位置坐标值增量ΔΨ_δX(t₁)、ΔΨ_δY(t₁)可按下式确定，

ΔΨ_δX(t₁)＝Ψ_δX(t₁)cosΔT_C-Ψ_δY(t₁)sinΔT_C-Ψ_δX(t₁)，(L-28)

ΔΨ_δY(t₁)＝Ψ_δY(t₁)cosΔT_C+Ψ_δX(t₁)sinΔT_C-Ψ_δY(t₁)。(L-29)

式中，t₁＝0。(L-30)

(3)当i_C＝125，即对应最后一个中间点，其相应的位置坐标值增量ΔΨ_δX(t₁₂₅)、ΔΨ_δY(t₁₂₅)无需计算，因为最后一个内接弦的终点即圆弧段Z_δC的终点，其位置坐标值是已知的，

Ψ_δX(t₁₂₆)＝0，(L-31)

Ψ_δY(t₁₂₆)＝50000.5。(L-32)

如果不设替代曲线，直接以圆弧段Z_C的内接弦拟合圆弧段Z_C，为了保证内接弦相对圆弧段Z_C的径向误差绝对值不超过ε_C，对本例ε_C为0.5，则应有

$\left|\mathrm{\&Delta;}{T}_C\right|\&le;\sqrt{\frac{8{\mathrm{\&epsiv;}}_C}{R_0}}.---(L-33)$

若取 $\left|\mathrm{\&Delta;}{T}_C\right|=\sqrt{\frac{8{\mathrm{\&epsiv;}}_C}{R_0}},---(L-34)$

则其取值是式(L-13)取值的

倍。由此可见，设定替代曲线将减小分段段数，从而减小插补运算量。

由于第一象限的1/4圆弧相对角XOY的平分线是对称的；因此，插补可以只针对圆心角为0～(π/4)部分的圆弧进行就可以了。圆心角为(π/4)～(π/2)部分的圆弧的位置坐标值或其增量可利用所述的对称性获得，计算量可以节省50％。

Claims (11)

1.一种插补方法，其所针对的所需路径或轮廓线Q上的点的位置坐标Ω_P(p＝1、2、3、……、m_Ω)中包括有一个或若干个坐标Ψ_k(k＝1、2、3、……、m_Ψ)，对应所述一个或若干个坐标Ψ_k(k＝1、2、3、……、m_Ψ)的坐标函数分别可以表示为以参数t为自变量、幅值分别为H_k(k＝1、2、3、……m_Ψ)、初始相位分别为α_k(k＝1、2、3、……、m_Ψ)、周期相同为(2π/ω)的正弦函数，其表达式为

Ψ_k(t)＝H_k sin(ωt+α_k)，(k＝1、2、3、……、m_Ψ)，(1-1)

所述曲线Q对应坐标Ψ_k(k＝1、2、3、……、m_Ψ)的坐标函数指的是描述所述曲线Q上的点的位置坐标Ψ_k(k＝1、2、3、……、m_Ψ)变化的以参数t为自变量的函数，

所述参数t可以是该曲线Q上的点的位置坐标Ω_P(p＝1、2、3、……、m_Ω)中的一个坐标，也可以是这些坐标之外的另一个参数，

所述曲线Q各坐标函数的定义域相同，定义域的二个端点分别与曲线Q二个已知点对应着相同的t值，

针对所需路径或轮廓线Q的插补就是针对其坐标函数Ω_P(t)(p＝1、2、3、……、m_Ω)的插补，

针对坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)的插补步骤包括，

(1)设定对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)，

(2)确定替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)定义域二个端点间的中间点，包括，

①确定所述中间点所对应的参数t的值或其增量Δt的值，

②确定所述中间点的个数，

(3)确定所述中间点的替代坐标函数值或其增量值，

(4)存储/输出运算结果，

所述的替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)是替代所需路径或轮廓线Q的替代曲线Q_δ对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的坐标函数，或说是坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)的替代坐标函数，

插补中确定的中间点将所述定义域分成分段，每个分段定义域将对应一个线性函数，所述各个分段定义域的端点所对应的线性函数值与替代坐标函数值相等，整个定义域将对应一个由各个分段定义域对应的线性函数组成的分段线性函数Ψ_δLk(t)(k＝1、2、3、……、m_Ψ)，而坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)将以分段线性函数Ψ_δLk(t)(k＝1、2、3、……、m_Ψ)拟合，

其特征在于：

(1)对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的替代坐标函数设定为与相应的坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)具有相同周期(2π/ω)、相同初始相位α_k(k＝1、2、3、……、m_Ψ)、相同定义域但不同幅值的正弦函数，其表达式为

Ψ_δk(t)＝H_δk sin(ωt+α_k)，(k＝1、2、3、……、m_Ψ)，(1-2)

式中，H_δk＝H_k+δ_k，         (k＝1、2、3、……、m_Ψ)，(1-3)

其中δ_k(k＝1、2、3、……、m_Ψ)为替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)幅值H_δk(k＝1、2、3、……、m_Ψ)相对相应的坐标函数Ψ_k(t)(k＝1、2、3、……、m_Ψ)幅值H_k(k＝1、2、3、……、m_Ψ)的幅值差，或说是替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)的幅值差，或说是对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的幅值差，

δ_k＝H_δk-H_k，     (k＝1、2、3、……、m_Ψ)，(1-4)

δ_k≥0，           (k＝1、2、3、……、m_Ψ)，(1-5)

(2)将对应所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)的各替代坐标函数Ψ_δk(t)(k＝1、2、3、……、m_Ψ)定义域等分，以等分分段的交点作为定义域的中间点，这些中间点所对应的替代坐标函数值增量ΔΨ_δk(t)(k＝1、2、3、……、m_Ψ)，按下述公式确定，

ΔΨ_δΛ(t_u+1)＝ΔΨ_δΛ(t_u)cosΔT_Λ+ΔΦ_δΛ(t_u+1)sinΔT_Λ，(1-6)

ΔΦ_δΛ(t_u+1)＝ΔΦ_δΛ(t_u)cosΔT_Λ-ΔΨ_δΛ(t_u)sinΔT_Λ，  (1-7)

式中，①Ψ_Λ表示所述坐标Ψ_k(k＝1、2、3、……、m_Ψ)中的某个坐标，Λ为序号k(k＝1、2、3、……、m_Ψ)中的某个序号，

②u+1为替代坐标函数Ψ_δΛ(t)定义域二个端点间的某个中间点的序号，u是与之相邻的前一个中间点或定义域的起点的序号，

以n_Λ表示替代坐标函数Ψ_δΛ(t)定义域等分分段的段数，以i_Λ(i_Λ＝1、2、3、……、n_Λ)作为分段的序号，以i_Λ(i_Λ＝1、2、3、……、n_Λ，n_Λ+1)作为包括定义域的二个端点及二个端点间的中间点在内的点的序号，二个端点分别对应着序号1及n_Λ+1，二个端点间的中间点分别对应着序号i_Λ(i_Λ＝2、3、……、n_Λ)，u+1就是序号i_Λ(i_Λ＝2、3、……、n_Λ)中的某一个序号，u就是序号i_Λ(i_Λ＝1，2、3、……、n_Λ)中与u+1相邻的前一个中间点或定义域起点的序号，u+2就是序号i_Λ(i_Λ＝1，2、3、……、n_Λ，n_Λ+1)中与u+1相邻的后一个中间点或定义域终点的序号，

所述定义域起点指的是与所述替代曲线Q_δ起点对应的定义域的端点，所述定义域的终点指的是与所述替代曲线Q_δ终点对应的定义域的端点，

③t_u+1为替代坐标函数Ψ_δΛ(t)定义域二个端点间序号为u+1的中间点所对应的参数t的值，

t_u为与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的参数t的值，

t_u+2为与序号为u+1的中间点相邻的后一个中间点或定义域终点所对应的参数t的值，

t₁为替代坐标函数Ψ_δΛ(t)定义域起点所对应的参数t的值，

t(n_Λ+1)为替代坐标函数Ψ_δΛ(t)定义域终点所对应的参数t的值，

④ΔT_Λ为替代坐标函数Ψ_δΛ(t)定义域各等分分段起点所对应的参数t的等效增量，或说是替代坐标函数Ψ_δΛ(t)定义域中序号为i_Λ(i_Λ＝1、2、3、……、n_Λ)的点所对应的参数t的等效增量，或说是替代坐标函数Ψ_δΛ(t)定义域各等分分段所对应的参数t的等效增量，其数值为所述参数t的增量Δτ_Λ的ω倍，是常数，

ΔT_Λ＝ωΔτ_Λ＝常数，            (1-8)

(i_Λ＝1、2、3、……、n_Λ)，        (1-9)

其中，t(i_Λ+1)(i_Λ＝1、2、3、……、n_Λ)为序号为(i_Λ+1)(i_Λ＝1、2、3、……、n_Λ)的点对应的参数t的值，

t(i_Λ)(i_Λ＝1、2、3、……、n_Λ)为序号为i_Λ(i_Λ＝1、2、3、……、n_Λ)的点对应的参数t的值，

⑤ΔΨ_δΛ(t_u+1)为替代坐标函数Ψ_δΛ(t)定义域中序号为u+1的中间点所对应的替代坐标函数值的增量，

ΔΨ_δΛ(t_u+1)＝Ψ_δΛ(t_u+2)-Ψ_δΛ(t_u+1)，         (1-10)

其中，Ψ_δΛ(t_u+2)为与序号为u+1的中间点相邻的后一个中间点或定义域终点所对应的替代坐标函数值，

Ψ_δΛ(t_u+1)为序号为u+1的中间点所对应的替代坐标函数值，

⑥ΔΨ_δΛ(t_u)为替代坐标函数Ψ_δΛ(t)定义域中与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的替代坐标函数值增量，

ΔΨ_δΛ(t_u)＝Ψ_δΛ(t_u+1)-Ψ_δΛ(t_u)，             (1-11)

其中，Ψ_δΛ(t_u+1)为序号为u+1的中间点所对应的替代坐标函数值，

Ψ_δΛ(t_u)为与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的替代坐标函数值，

⑦ΔΦ_δΛ(t_u+1)为替代坐标函数Ψ_δΛ(t)定义域中序号为u+1的中间点所对应的虚拟替代坐标函数值的增量，

ΔΦ_δΛ(t_u+1)＝Φ_δΛ(t_u+2)-Φ_δΛ(t_u+1)，        (1-12)

其中，Φ_δΛ(t_u+2)为与序号为u+1的中间点相邻的后一个中间点或定义域终点所对应的虚拟替代坐标函数值，

Φ_δΛ(t_u+1)为序号为u+1的中间点所对应的虚拟替代坐标函数值，

⑧ΔΦ_δΛ(t_u)为替代坐标函数Ψ_δΛ(t)定义域中与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的虚拟替代坐标函数值增量，

ΔΦ_δΛ(t_u)＝Φ_δΛ(t_u+1)-Φ_δΛ(t_u)，            (1-13)

其中，Φ_δΛ(t_u+1)为序号为u+1的中间点所对应的虚拟替代坐标函数值，

Φ_δΛ(t_u)为与序号为u+1的中间点相邻的前一个中间点或定义域起点所对应的虚拟替代坐标函数值，

⑨sinΔT_Λ为对应ΔT_Λ的正弦函数值，cosΔT_Λ为对应ΔT_Λ的余弦函数值，

上述⑦、⑧中所述的虚拟替代坐标函数Φ_δΛ(t)是一个与替代坐标函数Ψ_δΛ(t)对应的函数，Ψ_δΛ(t)与Φ_δΛ(t)是一组具有相同幅值H_δΛ、相同周期(2π/ω)、相同初始相位α_Λ及相同定义域的正弦函数与余弦函数，其表达式为

Ψ_δΛ(t)＝H_δΛ sin(ωt+α_Λ)，                   (1-14)

Φ_δΛ(t)＝H_δΛ cos(ωt+α_Λ)。，，               (1-15)。

2.如权利要求1所述的一种插补方法，其特征在于：

(1)所述坐标Ψ_k(k＝1、……、m_Ψ)共有二个坐标Ψ_k(k＝1、2)或Ψ₁、Ψ₂，对应Ψ₁、Ψ₂的坐标函数Ψ₁(t)、Ψ₂(t)可以表示为一组幅值分别为H₁与H₂、周期相同为(2π/ω)、初始相位相同为α₀的余弦函数与正弦函数，其表达式为

${\mathrm{\&Psi;}}_1\left(t\right)=H_1\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_1)=H_1\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_0+\frac{\mathrm{\&pi;}}{2})=H_1\cos (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_0),---(2-1)$

Ψ₂(t)＝H₂ sin(ωt+α₂)＝H₂ sin(ωt+α₀)，         (2-2)

式中， ${\mathrm{\&alpha;}}_1={\mathrm{\&alpha;}}_0+\frac{\mathrm{\&pi;}}{2},---(2-3)$

α₂＝α₀，                                         (2-4)

(2)对应所述坐标Ψ₁及Ψ₂的替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)为

Ψ_δ1(t)＝H_δ1 cos(ωt+α₀)，                              (2-5)

Ψ_δ2(t)＝H_δ2 sin(ωt+α₀)，                              (2-6)

式中，H_δ1＝H₁+δ₁，                                       (2-7)

H_δ2＝H₂+δ₂。                                             (2-8)。

3.如权利要求2所述的一种插补方法，其特征在于：

(1)对应所述坐标Ψ₁、Ψ₂的坐标函数Ψ₁(t)、Ψ₂(t)具有相同幅值R₀，函数Ψ₁(t)、Ψ₂(t)的表达式为

${\mathrm{\&Psi;}}_1\left(t\right)=H_1\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_1)=R_0\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_0+\frac{\mathrm{\&pi;}}{2})=R_0\cos (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_0),---(3-1)$

Ψ₂(t)＝H₂ sin(ωt+α₂)＝R₀ sin(ωt+α₀)，                 (3-2)

式中，H₁＝H₂＝R₀，                                         (3-3)

${\mathrm{\&alpha;}}_1={\mathrm{\&alpha;}}_0+\frac{\mathrm{\&pi;}}{2},---(3-4)$

α₂＝α₀，                                                 (3-5)

(2)对应所述坐标Ψ_X及Ψ_Y的替代坐标函数为Ψ_δ1(t)、Ψ_δ2(t)具有相同的幅值R_δ0，函数Ψ_δ1(t)、Ψ_δ2(t)的表达式为

Ψ_δ1(t)＝H_δ1 cos(ωt+α₀)＝R_δ0 cos(ωt+α₀)，           (3-6)

Ψ_δ2(t)＝H_δ2 sin(ωt+α₀)＝R_δ0 sin(ωt+α₀)，           (3-7)

式中，H_δ1＝H_δ2＝R_δ0，                                   (3-8)

H_δ1＝H₁+δ₁＝H_δ2＝H₂+δ₂＝R₀+δ_R，                       (3-9)

其中，δ₁＝δ₂＝δ_R，                                      (3-10)

(3)对应所述坐标Ψ₁、Ψ₂的替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域等分为相同的n_C个分段，或者说，替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域各分段起点对应相同的参数t的增量Δτ_C或等效增量ΔT_C，

n₁＝n₂＝n_C，                                               (3-11)

Δτ₁＝Δτ₂＝Δτ_C，                                      (3-12)

ΔT₁＝ΔT₂＝ΔT_C。                                         (3-13)。

4.如权利要求1至3中任何一项权利要求所述的一种插补方法，其特征在于：对应所述坐标其幅值差的取值、其替代坐标函数定义域等分分段起点所对应的参数t的等效增量的取值及替代坐标函数定义域等分分段段数的取值满足下述公式，

0≤δ_A|sin(ωt+α_A)|_MAX≤ε_A，                             (4-1)

$\left|{\mathrm{\&Delta;T}}_A\right|\&le;\sqrt{8\&times;\frac{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}},---(4-2)$

$n\&GreaterEqual;\mathrm{\&omega;}|t_{n+1}-t_1|\&times;\sqrt{\frac{1}{8}\&times;\frac{H_A\&times;{\left|\sin (\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\right|}_{\mathrm{MAX}}}{{\mathrm{\&epsiv;}}_A+{\mathrm{\&delta;}}_A\&times;{\left|\sin \right|(\mathrm{\&omega;t}+{\mathrm{\&alpha;}}_A)\left|\right|}_{\mathrm{MAX}}}},---(4-3)$

式中，①Ψ_A表示所述坐标中的某一个坐标，A是序号k(k＝1、2、3、……、m_Ψ)中的某个序号，对应坐标Ψ_A的坐标函数与替代坐标函数为

Ψ_A(t)＝H_A sin(ωt+α_A)，                                    (4-4)

Ψ_δA(t)＝H_δA sin(ωt+α_A)，                                (4-5)

②δ_A为对应坐标Ψ_A的幅值差，

δ_A＝H_δA-H_A，                                               (4-6)

③ΔT_A为替代坐标函数Ψ_δA(t)定义域等分分段起点所对应的参数t的等效增量，其数值为所述参数t的增量Δτ_A的ω倍，

ΔT_A＝ωΔτ_A，                                              (4-7)

④n表示替代坐标函数Ψ_δA(t)定义域等分分段的段数，

⑤t₁为替代坐标函数Ψ_δA(t)定义域起点所对应的参数t的值，

t_n+1为替代坐标函数Ψ_δA(t)定义域终点所对应的参数t的值，

⑥|sin(ωt+α_A)|_MAX为对应替代坐标函数Ψ_δA(t)定义域|sin(ωt+α_A)|的最大值，

⑦ε_A为以分段线性函数Ψ_δLA(t)拟合坐标函数Ψ_A(t)的误差绝对值|δΨ_δL-A(t)|的允许值，

δΨ_δL＝A(t)＝Ψ_δLA(t)-Ψ_A(t)，                            (4-8)

对应替代坐标函数Ψ_δA(t)的第i(i＝1、2、3、……、n)个分段定义域，所述的分段线性函数Ψ_δLA(t)表达式为

${\mathrm{\&Psi;}}_{\mathrm{\&delta;LA}}\left(t\right)={\mathrm{\&Psi;}}_{\mathrm{\&delta;A}}\left(t_i\right)+\frac{{\mathrm{\&Delta;\&Psi;}}_{\mathrm{\&delta;A}}\left(t_i\right)}{{\mathrm{\&Delta;T}}_A}(t-t_i),,$ (i＝1、2、3、……、n)，(4-9)

式中，(a)i(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域n个等分分段的序号，等分分段的交点就是定义域的中间点，以i(i＝1、2、3、……、n、n+1)作为包括定义域的二个端点及二个端点间的中间点在内的点的序号，二个端点分别对应着序号1及n+1，二个端点间的中间点分别对应着序号i(i＝2、3、……、n)，

(b)t_i(i＝1、2、3、……、n、n+1)为替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n、n+1)的点所对应的参数t的值，

(c)ΔΨ_δA(t_i)(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n)的点所对应的替代坐标函数值的增量，

ΔΨ_δA(t_i)＝Ψ_δA(t_i+1)-Ψ_δA(t_i)，(i＝1、2、3、……、n)， (4-10)

其中，Ψ_δA(t_i+1)(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域中与序号为i(i＝1、2、3、……、n)的点相邻的后一个点所对应的替代坐标函数值，

Ψ_δA(t_i)(i＝1、2、3、……、n)为替代坐标函数Ψ_δA(t)定义域序号为i(i＝1、2、3、……、n)的点所对应的替代坐标函数值。

5.如权利要求4所述的一种插补方法，其特征在于：对应所述坐标Ψ_A其幅值差δ_A取值为

δ_A＝ε_A/|sin(t+α_A)|_MAX。                                  (5-1)

6.如权利要求4所述的一种插补方法，其特征在于：所述|sin(t+α_A)|_MAX以1替换。

7.如权利要求6所述的一种插补方法，其特征在于：对应所述坐标Ψ_A其幅值差δ_A取值为ε_A，

δ_A＝ε_A。                                 (7-1)。

8.如权利要求3所述的一种插补方法，其特征在于：对应所述坐标Ψ₁、Ψ₂其幅值差δ_R的取值、其替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域等分分段起点所对应的参数t的等效增量|ΔT_C|的取值及其替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)定义域等分分段段数n_C的取值满足下述公式，

0＜δ_R≤ε_C，                              (8-1)

$\left|{\mathrm{\&Delta;T}}_C\right|\&le;\sqrt{\frac{8\&times;({\mathrm{\&epsiv;}}_C+{\mathrm{\&delta;}}_R)}{R_0}},---(8-2)$

$n_C\&GreaterEqual;\frac{\mathrm{\&omega;}|t(n_C+1)-t_1|}{4}\&times;\sqrt{\frac{R_0}{8\&times;({\mathrm{\&epsiv;}}_C+{\mathrm{\&delta;}}_R)}},---(8-3)$

式中，①t₁为替代坐标函数Ψ_δ1(t)或Ψ_δ2(t)定义域起点对应的参数t的值，t(n_C+1)为替代坐标函数Ψ_δ1(t)或Ψ_δ2(t)定义域终点对应的参数t的值，

②ε_C为以圆弧段Z_δC的等分圆弧分段内接弦拟合圆弧段Z_C的径向误差的绝对值的允许值，圆弧段Z_δC为替代坐标函数Ψ_δ1(t)、Ψ_δ2(t)在坐标轴Ψ₁和Ψ₂构成的直角坐标系Ψ₁OΨ₂的坐标平面上的图形，圆弧段Z_C为坐标函数Ψ₁(t)、Ψ₂(t)在坐标平面Ψ₁OΨ₂上的图形。

9.如权利要求8所述的一种插补方法，其特征在于：所述δ_R取值为ε_C，

δ_R＝ε_C。                                 (9-1)。

10.如权利要求1至3或权利要求5至9中任何一项权利要求所述的一种插补方法，其特征在于：

ω＝1。                                    (10-1)。

11.如权利要求4所述的一种插补方法，其特征在于：

ω＝1。                                    (11-11)。

Images