基于矩阵运算的平面Hilbert曲线向自由曲面上的映射.docVIP

基于矩阵运算的平面Hilbert曲线向自由曲面上的映射 　　摘要：基于矩阵变化描述生成平面Hilbert曲线的方法简单，易于编程实现，将其进行了正确的空间矢量化，使它能够映射在自由曲面上。分析了矩阵变化描述曲线的思想，结合自由曲面的形成过程，提出将平面Hilbert曲线通过与曲面参数域相对应的方法映射到曲面上的思想，建立了二者之间的映射关系，构建了映射算法，最后用一个具体实例验证了构建算法的可行性。该映射算法为Hilbert曲线的应用扩展到空间区域提供了方法支持。 　　Abstract: Plane Hilbert curve is described by matrix operation with the advantage of generating method simply and programming easily, in order to map the curve to free surface, the curve is made vector in space. The idea of generating curve is analyzed, combined with the formation of free surface, plane Hilbert space-filling curve was mapped to the surface through correspondence with the surface parameters domain, the mapping relation is established between curve and surface, mapping algorithm is constructed, finally, an example has been given to testify the feasibility of the algorithm. The algorithm is helpful to extend the application of Hilbert curve to space domain. 　　关键词：平面Hilbert曲线；矩阵运算；映射；自由曲面 　　Key words: plane Hilbert curve；matrix operation；map；free surface 　　中图分类号：TP391.73 文献标识码：A文章编号：1006-4311（2011）18-0155-02 　　0引言 　　Hilbert曲线是由意大利数学家Peano于1890年提出，德国数学家Hilbert于1891年首先给出了构造这种填充曲线的几何过程[1]。Hilbert曲线已有广泛的应用，例如在图像存储和检索，空间数据库索引等领域得到了成功的应用。近年来，Hilbert曲线应用于高质量数控加工的刀具路径规划方面被很多学者提出和研究[2]。因此研究Hilbert曲线有重要的理论意义和应用价值。 　　对于Hilbert曲线的生成，经典算法大部分是按照Hilbert本人初始的思想不断的四分一个正方形区域，从始点到终点递归地计算画线的位置的过程。这些方法比较抽象，编程实现繁琐[3]。针对这一缺陷，后来有学者提出了用矩阵变化的方法[4]来描述生成该曲线，该方法简单易于编程。本文针对矩阵法描述生成的平面Hilbert曲线特点设计了映射算法，将它在设计曲面上进行了空间矢量化，从而实现了不同阶次Hilbert曲线在自由曲面上的映射图。 　　不管用什么方法生成的平面Hilbert曲线，如何将它进行有效的映射，使它的应用可以扩展到空间区域也是非常有意义的研究，本文针对矩阵法描述生成的平面Hilbert曲线特点设计了映射算法，将它在设计曲面上进行了空间矢量化，从而实现了不同阶次Hilbert曲线在自由曲面上的映射图。 　　1平面Hilbert曲线的矩阵法描述 　　Hilbert曲线的定义是通过正方形逐次分割的标号次序给出的，例如，一单位正方形，可将其逐级分割为4个相等的小正方形，图1（a）和图1（b）分别为一、二级分割所得到的结果，依次类推，每个自然数k，对k-1次分割再做细分割，即得到4k个k级正方形，每个正方形的边长为1/2k，可记这些正方形为Qi（k），（i=1，2，…，4k），且它们的标号次序满足以下条件： 　　①每级Qi（k）包含点O（0，0）；②k级正方形Q4i-3（k），Q4i-2（k），Q4i-1（k），Q4i（k）包含在k-1级的正方形Q i（k-1）中，k1，i=1，2，…，4k-1；③任何下标为连续整数的正方形Qi（k）与Qi+1（k）具有公共边。将以上标号按从小到大的顺序，用折线连接每个小正方形的中