CN101751008B - 一种数控插补方法 - 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_ksin(ω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_Ω)的定义域二个端点间插入若干个中间点并确定所述中间点的坐标函数值。插补中确定的中间点将所述定义域分成分段，每个分段定义域将对应一个线性函数，所述各个分段定义域的端点所对应的线性函数值与坐标函数值相等，整个定义域将对应一个由各个分段定义域对应的线性函数组成的分段线性函数ψ_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(t)(k＝1、2、3、……、m_ψ)定义域二个端点间的中间点，包括， 

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

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

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

(3)存储/输出运算结果。 

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

(1)将对应所述坐标ψ_k(k＝1、2、3、……、m_ψ)的各坐标函数ψ_k(t)(k＝1、2、3、……、m_ψ)定义域等分，以等分分段的交点作为定义域的中间点， 

(2)所述中间点所对应的坐标函数值增量ΔΨ_k(t)(k＝1、2、3、……、m_Ψ)按下述公式确定， 

Δψ_Λ(t_u+1)＝Δψ_Λ(t_u)cosΔT_Λ+ΔΦ_Λ(t_u)sinΔT_Λ，(1-2) 

ΔΦ_Λ(t_u+1)＝ΔΦ_Λ(t_u)cosΔT_Λ-Δψ_Λ(t_u)sinΔT_Λ，(1-3) 

式中，①ψ_Λ表示所述坐标ψ_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-4) 

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

其中，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-6) 

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

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

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

Δψ_δΛ(t_u)＝ψ_Λ(t_u+1)-ψ_Λ(t_u)，(1-7) 

其中，ψ_Λ(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-8) 

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

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

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

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

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

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

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

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

Φ_Λ(t)＝H_Λcos(ωt+α_Λ)。(1-11) 

说明： 

(1)参见图2，如果所需路径或轮廓线的位置坐标只包括了二个个坐标t和ψ_A，且坐标轴t和ψ_A构成了直角坐标系tOψ_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)定义域的端点与中间点对应着圆弧段C_A相应的端点与中间点。图中C_A是第一象限的圆弧段。 

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

参见图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_A\cos [\mathrm{\ω}(t_u+\frac{\mathrm{\Δ}{\mathrm{\τ}}_A}{2})+{\mathrm{\α}}_A]=l_A\cos (\mathrm{\ω}{t}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A),---(1-12)$

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

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

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

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

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

$=\mathrm{\Δ}{\mathrm{\Ψ}}_A\left(t_u\right)\cos \mathrm{\Δ}{T}_A+\mathrm{\Δ}{\mathrm{\Φ}}_A\left(t_u\right)\sin \mathrm{\Δ}{T}_A,---(1-15)$

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

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

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

(4)对应所述坐标ψ_Λ的坐标函数ψ_Λ(t)其定义域等分分段起点所对应的参数t等效增量ΔT_Λ的正弦函数值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^{\′}}}{d{t}^{v^{\′}}}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_S,---(1-17)$

$\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^{\′\′}}}{d{t}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′\′}}+R_C,---(1-18)$

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

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

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

v′≥0，(1-21) 

②γ′的取值范围为 

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

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

v″≥0，(1-23) 

④γ″的取值范围为 

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

按公式(1-17)、(1-18)确定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-17)、(1-18)的依据是： 

$\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}{d{t}^2}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+\frac{1}{3!}{\left[\frac{d^3}{d{t}^3}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+$

$+\·\·\·\·\·\·+\frac{1}{v^{\′}!}{\left[\frac{d^{v^{\′}}}{d{t}^{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^{\′}}}{d{t}^{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^{\′}}}{d{t}^{v^{\′}}}\sin t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′}}+R_S,$

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

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

$\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}{d{t}^2}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^2+\frac{1}{3!}{\left[\frac{d^3}{d{t}^3}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^3+$

$+\·\·\·\·\·\·+\frac{1}{v^{\′\′}!}{\left[\frac{d^{v^{\′\′}}}{d{t}^{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^{\′\′}}}{d{t}^{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^{\′\′}}}{d{t}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left({\mathrm{\ΔT}}_{\mathrm{\Λ}}\right)}^{v^{\′\′}}+R_C,$

$=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^{\′\′}}}{d{t}^{v^{\′\′}}}\cos t\right]}_{t=0}{\left(\mathrm{\Δ}{T}_{\mathrm{\Λ}}\right)}^{v^{\′\′}}+R_C,---(1-26)$

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

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

Δψ_Λ(t₁)＝ψ_Λ(t₁)cosΔT_Λ+Φ_Λ(t₁)sinΔT_Λ-ψ_Λ(t₁)，(1-27) 

ΔΦ_Λ(t₁)＝Φ_Λ(t₁)cosΔT_Λ-ψ_Λ(t₁)sinΔT_Λ-Φ_Λ(t₁)。(1-28) 

其依据是： 

Δψ_Λ(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-29) 

ΔΦ_Λ(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-30) 

2、如第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) 

坐标函数ψ₁(t)、ψ₂(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-5)$

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

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

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

所述的椭圆弧段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)对应所述坐标Ψ₁、Ψ₂的坐标函数Ψ₁(t)、Ψ₂(t)定义域等分为相同的n_C个分段，或者说，坐标函数Ψ₁(t)、Ψ₂(t)定义域各分段起点对应相同的参数t的增量Δτ_C或等效增量ΔT_C， 

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

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

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

由所述特点可知， 

①在按第1点所述的公式(1-2)、(1-3)确定坐标函数值增量时，对应于坐标函数ψ₂(t)的虚拟坐标函数Φ₂(t)就是坐标函数ψ₁(t)， 

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

因此，在确定Δψ₁(t)、Δψ₂(t)过程中，无需为ψ₁(t)、ψ₂(t)另设虚拟坐标函数。 

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

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

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

ψ₁ ²+ψ₂ ²＝R₀ ²。(3-11) 

所述圆弧段Z_C可参见图5，ψ₁、ψ₂、R₀、α₀、t、Z_C分别对应图中的Φ_A、ψ_A、H_A、α_A、t、C_A；图中C_A是一个第一象限的圆弧段。 

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

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

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

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

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

②ΔT_A为坐标函数ψ_A(t)定义域等分分段起点所对应的参数t的等效增量，其数值为所述参数t的增量Δτ_A的ω倍， 

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

③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-5) 

对应坐标函数ψ_A(t)第i个分段定义域，所述的分段线性函数ψ_LA(t)表达式为 

${\mathrm{\Ψ}}_{\mathrm{LA}}\left(t\right)={\mathrm{\Ψ}}_A\left(t_i\right)+\frac{\mathrm{\Δ}{\mathrm{\Ψ}}_A\left(t_i\right)}{\mathrm{\Δ}{\mathrm{\τ}}_A}(t-t_i),,(i=\mathrm{1,2,3},\·\·\·\·\·\·,n),---(4-6)$

式中，(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-7) 

其中，ψ_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)的|ΔT_A|、n取值，将使得所述的拟合误差|δψ_L-A(t)|不超过允许值ε_A， 

参见图2。如果所需路径或轮廓线的位置坐标只包括了二个个坐标t和Ψ_A，且坐标轴t和Ψ_A构成了直角坐标系tOΨ_A；则所需路径或轮廓线即为平面tOΨ_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)确定|ΔT_A|及n的取值的依据是： 

(A)以ψ_LA(t)拟合ψ_A(t)的误差 

①以ψ_LA(t)拟合ψ_A(t)，其拟合误差记以δψ_L＝A(t)， 

δΨ_L＝A(t)＝Ψ_LA(t)-Ψ_A(t)。                         (4-8) 

②参见图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的增量Δτ_A或等效增量ΔT_A。以Z_A的内接弦构成的折线拟合Z_A，其拟合误差为δΨ_L＝A(t)；以C_A的内接弦构成的折线拟合C_A，在中点处其对应坐标轴Ψ_A拟合误差等于内接弦径向误差在坐标轴Ψ_A上的投影值δH_A，Ψ(t_M)。图6中，以t_M，u表示一个分段中点处的参数t的值，对应t_M，u的所述径向误差绝对值在坐标轴上的投影值即为|δH_A，Ψ(t_M，u)|。 

③以C_A的内接弦构成的折线拟合C_A，其径向误差绝对值的最大值|δH_A|发生在内接弦的中点t_M处，且 

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

其在坐标轴ψ_A上的投影为 

|δH_A，Ψ(t_M)|＝|δH_A|×|sin(ωt_M+α_A)|。                         (4-10) 

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

|δΨ_L-A(t_M.)|＝|δH_A，Ψ(t_M)|＝|δH_A|×|sin(ωt_M+α_A)|。         (4-11) 

⑤在定义域范围内，|δψ_L-A(t_M.)|的最大值|δψ_L-A(t_M.)|_MAX为 

|δψ_L-A(t_M.)|_MAX＝|δH_A|×|sin(ωt_M+α_A)|_MAX，(4-12) 

式中，|sin(ωt_M+α_A)|_MAX为对应坐标函数ψ_A(t)定义域|sin(ωt_M+α_A)|的最大值， 

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

|sin(ωt_M+α_A)|_MAX＝|sin(ωt+α_A)|_MAX，(4-13) 

式中，|sin(ωt+α_A)|_MAX为对应坐标函数ψ_A(t)定义域|sin(ωt+α_A)|的最大值， 

因此，式(4-12)可以改写为 

|δψ_L-A(t_M.)|_MAX＝|δH_A|×|sin(ωt+α_A)|_MAX。(4-14) 

(B)ΔT_A、n的取值 

ΔT_A、n的取值应使得Ψ_LA(t)拟合Ψ_A(t)的误差绝对值|δΨ_L-A(t)|的最大值|δΨ_L-A(t_M.)|_MAX不超过允许值ε_A， 

|δΨ_L-A(t_M.)|_MAX≤ε_A。                    (4-15) 

①由公式(4-9)、(4-14)、(4-15)，坐标函数ψ_A(t)定义域等分分段对应的参数t的等效增量|ΔT_A|的取值应满足下述公式， 

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

或，所述参数t的增量Δτ_A应满足 

$\left|\mathrm{\Δ}{\mathrm{\τ}}_A\right|＜\frac{1}{\mathrm{\ω}}\×\sqrt{\frac{8{\mathrm{\ϵ}}_A}{H_A\×{\left|\sin (\mathrm{\ωt}+{\mathrm{\α}}_A)\right|}_{\mathrm{MAX}}}},---(4-17)$

|ΔT_A|或|Δτ_A|的最大允许取值为 

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

或 $\left|\mathrm{\Δ}{\mathrm{\τ}}_A\right|={\left|\mathrm{\Δ}{\mathrm{\τ}}_A\right|}_{\mathrm{MAX}}=\frac{1}{\mathrm{\ω}}\×\sqrt{\frac{{8\mathrm{\ϵ}}_A}{H_A\×{\left|\sin (\mathrm{\ωt}+{\mathrm{\α}}_A)\right|}_{\mathrm{MAX}}}}.---(4-19)$

②坐标函数ψ_A(t)定义域等分分段段数n应满足下述公式， 

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

③如果对坐标函数ψ_A(t)的任何可能的定义域，上述公式中的|sin(ωt+α_A)|_MAX都以1替换， 

则 $\left|\mathrm{\Δ}{T}_A\right|＜\sqrt{\frac{8{\mathrm{\ϵ}}_A}{H_A}},---(4-21)$

|ΔT_A|或|Δτ_A|的最大允许取值为 

$\left|\mathrm{\Δ}{T}_A\right|={\left|\mathrm{\Δ}{T}_A\right|}_{\mathrm{MAX}}=\sqrt{\frac{8{\mathrm{\ϵ}}_A}{H_A}},---(4-22)$

或 $\left|\mathrm{\Δ}{\mathrm{\τ}}_A\right|={\left|\mathrm{\Δ}{\mathrm{\τ}}_A\right|}_{\mathrm{MAX}}=\frac{1}{\mathrm{\ω}}\×\sqrt{\frac{8{\mathrm{\ϵ}}_A}{H_A}}.---(4-23)$

此时，对任何可能的定义域，都可做到ψ_LA(t)拟合ψ_A(t)的误差的绝对值|δψ_L-A(t)|不超过允许值ε_A。 

5、如第4点所述的一种插补方法，其特点为：对应所述坐标ψ_A确定ΔT_A、n所应满足的公式中的正弦函数绝对值的最大值|sin(ωt+α_A)|_MAX以1替换。 

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

$\left|\mathrm{\Δ}{T}_A\right|＜\sqrt{\frac{8{\mathrm{\ϵ}}_A}{H_A}},---(5-1)$

$n\≥\mathrm{\ω}|t_{n+1}-t_1|\×\sqrt{\frac{H_A}{{8\mathrm{\ϵ}}_A}}.---(5-2)$

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

6、如第1至3点或第5至9点中任何一点所述的一种插补方法，其特点为：对应所述坐标的坐标函数中 

ω＝1。(6-1) 

也就是说，所述坐标对应的正弦坐标函数的周期为2π。此时，所述坐标对应的坐标函数定义域各等分分段起点所对应的参数t的等效增量，等于所述参数t的增量；以ψ_A表示各点所述坐标中的某一个坐标，则有 

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

7、如第4点所述的一种插补方法，其特点为：对应所述坐标的坐标函数中 

ω＝1。(7-1) 

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

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

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

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

【附图说明】

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

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

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

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

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

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

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

【具体实施方式】

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

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

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

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

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

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

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

2→i_Λ，(J1-3) 

(4)确定序号为i_Λ的中间点所对应的坐标函数值增量Δψ_Λ(t)的值， 

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

(6)判定 

i_Λ＝n_Λ？(J1-4) 

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

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

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

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

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

(8)插补完成。 

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

ψ_A(t)＝H_Asin(ωt+α_k)。(J2-1) 

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

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

ψ_A(t)＝H_Asin(ωt+α_k)，(J3-1) 

n是坐标函数ψ_A(t)定义域的等分分段数。正弦曲线段Z_A的n个等分分段的内接弦构成折线Z_LA，相应的分段线性函数ψ_LA(t)为 

${\mathrm{\Ψ}}_{\mathrm{LA}}\left(t\right)={\mathrm{\Ψ}}_A\left(t_i\right)+\frac{\mathrm{\Δ}{\mathrm{\Ψ}}_A\left(t_i\right)}{\mathrm{\Δ}{\mathrm{\τ}}_A}(t-t_i),,(t_i＜t＜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)的定义域， 

③Δτ_A是坐标函数Ψ_A(t)定义域各等分分段起点所对应的参数t的增量，Δτ_A的ω倍等于所述参数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-4) 

其中，ψ_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_LA之间的关系。图中以u表示所示分段的序号，该分段对应的参数t的区间为[t_u，t_u+1]，[t_u，t_u+1]的中点以t_M，u表示，ψ_A(t_M，u)、ψ_LA(t_M，u)分别是对应于t_M，u函数ψ_A(t)、ψ_LA(t)的值，δψ_L-A(t_M，u)是对应于点t_M，u以ψ_LA(t)拟合ψ_A(t)的误差δψ_L＝A(t)， 

δψ_L-A(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)。 

坐标函数ψ_A(t)、Φ_A(t)定义域为[t₁，t_n+1]，该定义域对应的圆弧曲线段记以C_A；图中C_A是一个第一象限的圆弧段。 

六、参见图6，这是圆弧曲线局部示意图。示意图表示圆弧段C_A在第一象限的一个分段与其内接弦之间的关系。 

①表明以圆弧段C_A内接弦构成的折线拟合圆弧段C_A的关系。 

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

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

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

$\mathrm{\Δ}{\mathrm{\Φ}}_A\left(t_u\right)=-l_A\sin [\mathrm{\ω}(t_u+\frac{\mathrm{\Δ}{\mathrm{\τ}}_A}{2})+{\mathrm{\α}}_A]=-l_A\sin (\mathrm{\ω}{t}_u+\frac{\mathrm{\Δ}{T}_A}{2}+{\mathrm{\α}}_A),---(J6-2)$

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

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

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

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

曲线的坐标函数为 

${\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),---(J7-2)$

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

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

α₂＝α₀，(J7-5) 

椭圆曲线Q_E上的点与坐标轴原点O的连线相对ψ_X轴的夹角即为(ωt+α₀)。坐标函数ψ₁(t)、ψ₂(t)定义域所对应的椭圆曲线段记以Z_E，图中Z_E是一个第一象限椭圆弧段。 

八、实施例 

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

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

曲线的坐标函数为 

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

ψ_Y(t)＝R₀sint，(L-3) 

R₀＝50000。(L-5) 

圆弧段Z_C上的点与圆心O的连线相对轴ψ_X的夹角即为参数t。 

要求以折线拟合圆弧段Z_C的最大径向误差绝对值不超过ε_C， 

ε_C＝0.5。(L-6) 

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

以圆弧段Z_C等分分段的内接弦拟合圆弧段Z_C。其插补过程如下： 

1、确定内接弦对应的圆心角Δ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-7) 

如果要求Z_C内接弦拟合Z_C的径向误差绝对值的最大值|δR₀|不超过0.5，则 

$\left|\mathrm{\δ}{R}_0\right|=\frac{1}{8}{R}_0{\left(\mathrm{\Δ}{T}_C\right)}^2＜{\mathrm{\ϵ}}_C=0.5,---(L-8)$

相应地， 

$\mathrm{\Δ}{T}_C＜\sqrt{\frac{8\×{\mathrm{\ϵ}}_C}{R_0}}=\sqrt{\frac{8\×0.5}{50000}}=0.00894,---(L-9)$

取ΔT_C＝0.00894。(L-10) 

(2)确定圆弧段Z_C分段 

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

$n_{\mathrm{CL}}=\frac{\mathrm{\π}}{2\×\mathrm{\Δ}{T}_C}=175.62.---(L-11)$

分段段数只能是整数，取分段段数的实际值为 

n_C＝176。(L-12) 

相应的ΔT_C取值为 

$\mathrm{\Δ}{T}_C=\frac{\mathrm{\π}}{2\×n_C}=0.00892,---(L-13)$

此时，拟合的径向误差误差绝对值的最大值 

$\left|\mathrm{\δ}{R}_0\right|=\frac{1}{8}{R}_0{\left(\mathrm{\Δ}{T}_C\right)}^2=0.498＜{\mathrm{\ϵ}}_C=0.5.---(L-14)$

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

以i_C(i_C＝1、2、3、……、176、177)表示构成拟合圆弧Z_C的折线的176个内接弦的起点及最后一个内接弦的终点在内的各点的序号，也就是坐标函数ψ_X(t)、ψ_Y(t)定义域的二个端点及二个端点间的中间点在内的点的序号，二个端点分别对应着序号1及177，二个端点间的中间点分别对应着序号 

i_C＝2、3、……、176。(L-15) 

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

ψ_X，1＝R₀＝50000，(L-16) 

ψ_Y，1＝0。(L-17) 

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

Δψ_X(t_V+1)＝Δψ_X(t_V)cosΔT_C-Δψ_Y(t_V)sinΔT_C，(L-18) 

Δψ_Y(t_V+1)＝Δψ_Y(t_V)cosΔT_C+Δψ_X(t_V)sinΔT_C。(L-19) 

相应的各中间点的位置坐标值ψ_X(t_h+1)、ψ_Y(t_h+1)由下式确定 

ψ_X(t_V+1)＝ψ_X(t_V)+Δψ_X(t_V)，(L-20) 

ψ_Y(t_V+1)＝ψ_Y(t_V)+Δψ_Y(t_V)。(L-21) 

式(L-18)、(L-19)中， 

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

$\sin \mathrm{\Δ}{T}_C=\mathrm{\Δ}{T}_C-\frac{1}{6}\×{\left(\mathrm{\Δ}{T}_C\right)}^3,---(L-22)$

$\cos \mathrm{\Δ}{T}_C=1-\frac{1}{2}{\left(\mathrm{\Δ}{T}_C\right)}^2+\frac{1}{24}{\left(\mathrm{\Δ}{T}_C\right)}^4,---(L-23)$

②根据公式(1-27)、(1-28)，圆弧段Z_C起点的位置坐标值增量Δψ_X(t₁)、Δψ_Y(t₁)可按下式确定， 

Δψ_X(t₁)＝＝ψ_X，1cosΔT_C+ψ_Y，1sinΔT_C-ψ_X，1，(L-24) 

Δψ_Y(t₁)＝＝ψ_Y，1cosΔT_C-ψ_X，1sinΔT_C-ψ_Y，1。(L-25) 

(3)当i_C＝176，即对应最后一个中间点，其相应的位置坐标值增量Δψ_X，176、Δψ_Y，176无需计 算，因为最后一个内接弦的终点即所需圆弧段Z_C的终点，其位置坐标值已给定， 

ψ_X，177＝0，(L-26) 

ψ_Y，177＝50000。(L-27) 

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

Claims (7)

1.一种数控插补方法，其所针对的所需路径或轮廓线Q上的点的m_Ω个位置坐标Ω_P(p＝1、2、3、……、m_Ω)中包括有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_Ψ)， 

所需路径或轮廓线Q对应坐标Ψ_k(k＝1、2、3、……、m_Ψ)的坐标函数指的是描述所需路径或轮廓线Q上的点的位置坐标Ψ_k(k＝1、2、3、……、m_Ψ)变化的以参数t为自变量的函数，参数t或说是所需路径或轮廓线Q上的点的位置坐标Ω_P(p＝1、2、3、……、m_Ω)中的一个坐标，或说是这些坐标之外的另一个参数， 

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

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

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

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

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

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

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

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

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

插补中确定的中间点将所述定义域分成分段，每个分段定义域将对应一个线性函数，所述各个分段定义域的端点所对应的线性函数值与坐标函数值相等，整个定义域将对应一个由各个分段定义域对应的线性函数组成的分段线性函数Ψ_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(t)(k＝1、2、3、……、m_Ψ)按下述公式确定， 

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

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

式中，①Ψ_Λ表示所述坐标Ψ_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_Λ＝ωΔτ_Λ＝常数， 

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

其中，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)， 

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

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

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

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

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

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

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

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

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

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

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

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

其中，Φ_Λ(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+α_Λ)， 

Φ_Λ(t)＝H_Λcos(ωt+α_Λ)。                 

2.如权利要求1所述的一种插补方法，其特征在于：所述坐标Ψ_k(k＝1、……、m_Ψ)共有二个坐标Ψ_k(k＝1、2)或Ψ₁、Ψ₂，对应Ψ₁、Ψ₂的坐标函数Ψ₁(t)、Ψ₂(t)可以表示为一组幅值分别为H₁、H₂、周期相同为2π/ω、初始相位相同为α₀的余弦函数与正弦函数，其表达式为 

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

式中，

α₂＝α₀。 

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

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

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

式中，H₁＝H₂＝R₀， 

α₂＝α₀， 

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

n₁＝n₂＝n_C， 

Δτ₁＝Δτ₂＝Δτ_C， 

ΔT₁＝ΔT₂＝ΔT_C。 

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

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

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

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

ΔT_A＝ωΔτ_A， 

③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)， 

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

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

式中，(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)， 

其中，Ψ_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(t)定义域|sin(ωt+α_A)|的最大值|sin(ωt+α_A)|_MAX以1替换。 

6.如权利要求1至3或权利要求5中任何一项权利要求所述的一种插补方法，其特征在于：对应所述坐标的坐标函数中 

ω＝1。 

7.如权利要求4所述的一种插补方法，其特征在于：对应所述坐标的坐标函数中 

ω＝1。 

Images