一种圆度误差评估的计算几何方法技术

本发明专利技术涉及测试计量技术领域，且公开了一种圆度误差评估的计算几何方法，步骤一：首先利用格雷厄姆(Graham)方法求得平面点集P中位于其凸壳上的有序点集P1；步骤二：从点集P1中删除位于凸壳边的中间点，且每删除一个点，之后的点序号依次减一，所剩点集即为决定平面点集凸壳的最小有序点集，记为P2；步骤三：在点集P2中找出同时满足条件min(∠p2*p1*p3*)(p1*∈P2\{p2*,p3*})、min(∠p1*p2*p3*)(p2*∈P2\{p1*,p3*})及min(∠p1*p3*p2*)(p3*∈P2\{p1*,p2*})的三个点p1*,p2*,p3*；步骤四：Δp1*p2*p3*的最小覆盖圆，即为平面点集P的最小覆盖圆，结束。该圆度误差评估的计算几何方法，不仅能精确确定最小覆盖圆圆心位置和半径长度，而且方法的原理简单，解结构固定，计算效率高，容易编程实现，适合大规模并行计算和后续有需要严格分析的场合。的场合。

【技术实现步骤摘要】
一种圆度误差评估的计算几何方法


[0001]本专利技术涉及测试计量
，具体为一种圆度误差评估的计算几何方法。

技术介绍

[0002]在测试计量领域，最小覆盖圆方法(简称3M法)是评定圆度误差的重要方法。
[0003]对于3M法的研究已有诸多报道，但是目前有效的方法并不多见，各类方法有自己的局限性。函数逼近的3M方法得到的是最小覆盖圆圆心与半径的近似值，严格地说所求并不是真正最小的覆盖圆。最优化的方法把圆心位置作为优化参数，把外接圆的半径作为优化目标，采用各种有效算法来搜寻最小覆盖圆的圆心，这种方法耗时长、效率较低，得到的也不是准确结果。另外，基于二次规划的3M法原理深奥，给使用者造成理解和编程困难，计算效率也不高。为此，提出一种圆度误差评估的计算几何方法。

技术实现思路

[0004](一)解决的技术问题
[0005]针对现有技术的不足，本专利技术提供了一种圆度误差评估的计算几何方法，不仅能精确确定最小覆盖圆圆心位置和半径长度，而且方法的原理简单，解结构固定，计算效率高，容易编程实现，适合大规模并行计算和后续有需要严格分析的场合，解决了
技术介绍
提出的问题。
[0006](二)技术方案
[0007]为实现上述的目的，本专利技术提供如下技术方案：一种圆度误差评估的计算几何方法，步骤一：首先利用格雷厄姆(Graham)方法求得平面点集P中位于其凸壳上的有序点集P1，具体过程如下
[0008]ⅰ
.将P中y坐标最小的点设为p1，如果这样的点有多个，将其中x坐标最小的那个点设为p1，p1也称为基点，其余的点p
i
按p
i
p1与水平方向的夹角大小(由小到大)进行排序，如果出现有多个点对应同一夹角，进一步对这些点按p
i
到p1的距离(由近到远)进行排序，得到点序列p1，p2，
…
，p
n
，其中p1和p
n
必为凸壳顶点，而p2必位于凸壳上，设p
n+1
＝p1；
[0009]ⅱ
.从p3，p4，
…
，p
n
‑1中删除不是处于凸壳上的点，方法如下：
[0010]①
k＝4
[0011]②
j＝2
[0012]③
如果p1和p
k
位于线段p
k
‑
j+1
p
k
‑
j
两侧，则删除p
k
‑1，后续点编号减1，k＝k
‑
1，j＝j
‑
1，否则p
k
‑1暂时保留为凸壳点
[0013]④
j＝j+1，转到
③
，直到j＝k
‑1[0014]⑤
k＝k+1，转到
②
，直到k＝n+1(注意，因为节点编号随着节点的删除是在不断减少的，所以在算法过程中n也是在不断减少的)
[0015]经过以上删除过程后，剩下的有序点集记为P1；
[0016]步骤二：从点集P1中删除位于凸壳边的中间点(如果存在的话)，且每删除一个点，
之后的点序号依次减一，所剩点集即为决定平面点集凸壳的最小有序点集，记为P2，如果card(P2)＝1，平面点集P无最小覆盖圆，结束；如果card(P2)＝2，以P2中两点的连线为直径的圆既为平面点集P的最小覆盖圆，结束；如果card(P2)＝3，将P2中的三点分别记为p1*，p2*，p3*，转到步骤四，否则继续执行步骤三；
[0017]步骤三：
[0018]在点集P2中找出同时满足条件min(∠p2*p1*p3*)(p1*∈P2\{p2*,p3*})、min(∠p1*p2*p3*)(p2*∈P2\{p1*,p3*})及min(∠p1*p3*p2*)(p3*∈P2\{p1*,p2*})的三个点p1*,p2*,p3*。具体可通过如下迭代过程实现，先从P2中随机选定三个点，分别记为p
10
,p
20
,p
30
，并对其依次循环检测更新，在对某个点迭代时固定其余的两个点，如果有与之构成更小角度的点就用该点替代迭代点，直到三个点都找不到可替代的点，迭代终止，此三点即为所求，分别记为p1*,p2*,p3*；
[0019]步骤四：Δp1*p2*p3*的最小覆盖圆，即为平面点集P的最小覆盖圆，结束。
[0020]进一步优选的，所述步骤1中的格雷厄姆方法的特点是除了位于相对基点(p1)最大水平夹角的凸壳边的中间点不被保留外，其它凸壳边的中间点(如果存在的话)会被保留下来。
[0021]进一步优选的，所述步骤一到步骤四中，由于点集P有n个点，先看第一步的步骤
ⅱ
中转移到
②
的次数不会超过n，每一个顶点至多删去1次，删去顶点的个数也不可能超过n，因此第一步的步骤
ⅱ
是线性时间复杂度，第一步的步骤
ⅰ
很显然涉及到角度的排序问题，要计算n
‑
1个夹角，并按夹角分类，计算每个夹角只需要常数时间，计算n
‑
1个夹角耗费线性时间，分类需要时间O(n log n)，因此，第一步格雷厄姆算法部分总的时间复杂性为O(n log n)，第二步、第三步处理的集合点数相对原始集合P来说少的多，且为线性时间开销，第四步为常数时间，故所提计算方法总的时间复杂度近似为O(n log n)。
[0022](三)有益效果
[0023]与现有技术相比，本专利技术提供了一种圆度误差评估的计算几何方法，具备以下有益效果：
[0024]该圆度误差评估的计算几何方法，基于计算几何思想给出了一种高效实用的圆度误差评估的方法，所提出的方法不仅能精确确定最小覆盖圆圆心位置和半径长度，而且方法的原理简单，解结构固定，计算效率高，容易编程实现，适合大规模并行计算和后续有需要严格分析的场合，适合于需要精确评估和大规模计算机械部件圆度误差的设备制造和在线检测场合。
附图说明
[0025]图1为本专利技术提出的第一种测试结构图；
[0026]图2为本专利技术提出的第二种测试结构图。
[0027][0028]具体实施方式
[0029]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完
整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0030]请参阅图1
‑
2，一种圆度误差评估的计算几何方法，步骤一：首先利用格雷厄姆(Graham)方法求得平面点集P中位于其凸壳上的有序点集P1，具体过程如下
[0031]ⅰ
.将P中y坐标最小的点设为p1，如果这样的点有多个，将其中x坐标最本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种圆度误差评估的计算几何方法，其特征在于，具体步骤如下：步骤一：首先利用格雷厄姆方法求得平面点集P中位于其凸壳上的有序点集P1，具体过程如下
ⅰ
.将P中y坐标最小的点设为p1，如果这样的点有多个，将其中x坐标最小的那个点设为p1，p1也称为基点，其余的点p
i
按p
i
p1与水平方向的夹角大小进行排序，如果出现有多个点对应同一夹角，进一步对这些点按p
i
到p1的距离进行排序，得到点序列p1，p2，
…
，p
n
，其中p1和p
n
必为凸壳顶点，而p2必位于凸壳上，设p
n+1
＝p1；
ⅱ
.从p3，p4，
…
，p
n
‑1中删除不是处于凸壳上的点，方法如下：
①
k＝4
②
j＝2
③
如果p1和p
k
位于线段p
k
‑
j+1
p
k
‑
j
两侧，则删除p
k
‑1，后续点编号减1，k＝k
‑
1，j＝j
‑
1，否则p
k
‑1暂时保留为凸壳点
④
j＝j+1，转到
③
，直到j＝k
‑1⑤
k＝k+1，转到
②
，直到k＝n+1经过以上删除过程后，剩下的有序点集记为P1；步骤二：从点集P1中删除位于凸壳边的中间点，且每删除一个点，之后的点序号依次减一，所剩点集即为决定平面点集凸壳的最小有序点集，记为P2，如果card(P2)＝1，平面点集P无最小覆盖圆，结束；如果card(P2)＝2，以P2中两点的连线为直径的圆既为平面点集P的最小...

【专利技术属性】
技术研发人员：李翔，
申请(专利权)人：广州工商学院，
类型：发明
国别省市：