德国数学家David Hilbert在1891年构造了一种曲线,首先把一个正方形等分成四个小正方形,依次从西北角的正方形中心出发往南到西南正方形中心,再往东到东南角的正方形中心,再往北到东北角正方形中心,这是一次迭代;如果对四个小正方形继续上述过程,往下划分,反复进行,最终就得到一条可以填满整个正方形的曲线,这就是Hibert曲线。其生成过程如图1所示。

图1  Hilbert曲线的生成

Hilbert曲线可以采用递归过程实现,在递归处理时,连接中点的方式有4种,如图2所示。

图2  连接中心点的4种方式

设正方形左上角的顶点坐标为(x1,y1),右下角顶点坐标为(x2,y2)。若将方式(3)的正方形左上角坐标置为(x2,y2),右下角坐标置为(x1,y1),则方式(3)等同于方式(1),相当于旋转180°;同理,方式(4)等同于方式(2)。因此,4种连接中心点的方式可以看成(1)和(2)两种。

两种连线方式的连线走向及下一次扩展的方式如图3所示。

图3  两种连线方式走向及扩展

其中,方式(1)的四个中心点坐标分别为:

①(x1+dx/4,y1+dy/4)           ②(x1+dx/4, y1+3*dy/4)

③ (x1+3*dx/4, y1+3*dy/4)     ④(x1+3*dx/4,y1+dy/4)   (dx=x2-1,dy=y2-y1)

方式(2)的四个中心点坐标分别为:

①(x1+dx/4,y1+dy/4)           ②(x1+3*dx/4,y1+dy/4)

③ (x1+3*dx/4, y1+3*dy/4)     ④(x1+dx/4,  y1+3*dy/4)

为此,引入一个标识变量s,s=1表示方式(1),s=-1表示方式(2),这样两种方式的中心点坐标可以统一表示为:

①(x1+dx/4,y1+dy/4)        ②(x1+(2-s)*dx/4,  y1+(2+s)*dy/4)

③(x1+3*dx/4, y1+3*dy/4)      ④(x1+(2+s)*dx/4,y1+(2-s)*dy/4)

递归扩展时,方式(1)中4个小正方形的扩展方式分别是:方式(2)、方式(1)、方式(1)和方式(4)(注意:给定两个顶点坐标顺序调整后等同于方式(2));方式(2)中4个小正方形的扩展方式分别是:方式(1)、方式(2)、方式(2)和方式(3)。

编写如下的HTML代码。

<!DOCTYPE html>

<head>

<title>Hilbert曲线</title>

</head>

<body>

<canvas id="myCanvas" width="500" height="500" style="border:3px double #996633;">

</canvas>

<script type="text/javascript">

var canvas = document.getElementById('myCanvas');

var ctx = canvas.getContext('2d');

var depth=5;

ctx.lineWidth = 2;

ctx.strokeStyle = "red";

ctx.beginPath();

ctx.moveTo(50+400/Math.pow(2,depth+1),50+400/Math.pow(2,depth+1));

drawShapes(depth,1,50,50,450,450);

ctx.stroke();

function drawShapes(n,s,x1,y1,x2,y2)

{

dx = x2 - x1,

dy = y2 - y1;

if (n>1)

{

if(s>0)

{

drawShapes(n-1,-1,x1,y1,(x1+x2)/2,(y1+y2)/2);

drawShapes(n-1,1,x1,(y1+y2)/2,(x1+x2)/2,y2);

drawShapes(n-1,1,(x1+x2)/2,(y1+y2)/2,x2,y2);

drawShapes(n-1,-1,x2,(y1+y2)/2,(x1+x2)/2,y1);

}

else

{

drawShapes(n-1,1,x1,y1,(x1+x2)/2,(y1+y2)/2);

drawShapes(n-1,-1,(x1+x2)/2,y1,x2,(y1+y2)/2);

drawShapes(n-1,-1,(x1+x2)/2,(y1+y2)/2,x2,y2);

drawShapes(n-1,1,(x1+x2)/2,y2,x1,(y1+y2)/2);

}

}

if (n==1)

{

ctx.lineTo(x1+dx/4,y1+dy/4);

ctx.lineTo(x1+(2-s)*dx/4,  y1+(2+s)*dy/4);

ctx.lineTo(x1+3*dx/4, y1+3*dy/4);

ctx.lineTo(x1+(2+s)*dx/4,y1+(2-s)*dy/4);

}

}

</script>

</body>

</html>

在浏览器中打开包含这段HTML代码的html文件,可以看到在浏览器窗口中绘制出如图4所示的Hilbert曲线。

图4  递归深度maxdepth =5的Hilbert曲线

上面的程序需要推出方式(一)和方式(二)的坐标统一形式,还需注意方式(3)和方式(4)与方式(一)和方式(二)的同一性。

由于Hilbert曲线可以看成是4种方式进行组合,因此可以直接对4种方式编写递归过程。编写如下的HTML文件。

<!DOCTYPE html>

<head>

<title>Hilbert曲线</title>

</head>

<body>

<canvas id="myCanvas" width="500" height="500" style="border:3px double #996633;">

</canvas>

<script type="text/javascript">

var canvas = document.getElementById('myCanvas');

var ctx = canvas.getContext('2d');

ctx.lineWidth = 2;

ctx.strokeStyle = "red";

ctx.beginPath();

var depth=5;    //  递归深度

var h=400/Math.pow(2,depth);

var x = 50+h;

var y = 50+h;

ctx.moveTo(x,y);

One(depth);

ctx.stroke();

function One(n)   // 方式(1)的递归调用

{

if(n > 0)

{

Two(n-1);

ctx.lineTo(x, y+h); y+=h;

One(n-1);

ctx.lineTo(x+h, y); x+=h;

One(n-1);

ctx.lineTo(x, y-h); y-=h;

Four(n-1);

}

}

function Two(n)   // 方式(2)的递归调用

{

if(n > 0)

{

One(n-1);

ctx.lineTo(x+h, y); x+=h;

Two(n-1);

ctx.lineTo(x, y+h); y+=h;

Two(n-1);

ctx.lineTo(x-h, y); x-=h;

Three(n-1);

}

}

function Three(n)   // 方式(3)的递归调用

{

if(n > 0)

{

Four(n-1);

ctx.lineTo(x, y-h);  y-=h;

Three(n-1);

ctx.lineTo(x-h, y);  x-=h;

Three(n-1);

ctx.lineTo(x, y+h);  y+=h;

Two(n-1);

}

}

function Four(n)  // 方式(4)的递归调用

{

if(n > 0)

{

Three(n-1);

ctx.lineTo(x-h,y);      x-=h;

Four(n-1);

ctx.lineTo(x, y-h);     y-=h;

Four(n-1);

ctx.lineTo(x+h, y); x+=h;

One(n-1);

}

}

</script>

</body>

</html>

在浏览器中打开包含这段HTML代码的html文件,可以看到在浏览器窗口中绘制出如图5所示的Hilbert曲线。

图5  调用One(depth)时绘制的图形

将程序中的调用语句“One(depth)”改写成“Two(depth)”,则在浏览器窗口中绘制出如图6所示的Hilbert曲线。这个图形可以看成是图5向左旋转90°得到的。实际上,由图2可知,将方式(一)的图形向左旋转90°得到的就是方式(二)的图形。

图6  调用Two(depth)时绘制的图形

将程序中调用语句“One(depth)”改写成“Three(depth)”,同时修改初始坐标为

“var x = 450-h;   var y = 450-h;”,则在浏览器窗口中绘制出如图7所示的Hilbert曲线。

图7  调用THree(depth)时绘制的图形

将程序中调用语句“One(depth)”改写成“Four(depth);”,同时修改初始坐标为

“var x = 450-h;   var y = 450-h;”,则在浏览器窗口中绘制出如图8所示的Hilbert曲线。

图8  调用Four(depth)时绘制的图形

将Hilbert曲线的生成过程进行动画展示,编写如下的HTML代码。

<!DOCTYPE>

<html>

<head>

<title>Hilbert曲线</title>

</head>

<body>

<canvas id="myCanvas" width="500" height="500" style="border:3px double #996633;"></canvas>

<script type="text/javascript">

var canvas = document.getElementById('myCanvas');

var ctx = canvas.getContext('2d');

var depth=1;

function drawShapes(n,s,x1,y1,x2,y2)

{

dx = x2 - x1,

dy = y2 - y1;

if (n>1)

{

if(s>0)

{

drawShapes(n-1,-1,x1,y1,(x1+x2)/2,(y1+y2)/2);

drawShapes(n-1,1,x1,(y1+y2)/2,(x1+x2)/2,y2);

drawShapes(n-1,1,(x1+x2)/2,(y1+y2)/2,x2,y2);

drawShapes(n-1,-1,x2,(y1+y2)/2,(x1+x2)/2,y1);

}

else

{

drawShapes(n-1,1,x1,y1,(x1+x2)/2,(y1+y2)/2);

drawShapes(n-1,-1,(x1+x2)/2,y1,x2,(y1+y2)/2);

drawShapes(n-1,-1,(x1+x2)/2,(y1+y2)/2,x2,y2);

drawShapes(n-1,1,(x1+x2)/2,y2,x1,(y1+y2)/2);

}

}

if (n==1)

{

ctx.lineTo(x1+dx/4,y1+dy/4);

ctx.lineTo(x1+(2-s)*dx/4,  y1+(2+s)*dy/4);

ctx.lineTo(x1+3*dx/4, y1+3*dy/4);

ctx.lineTo(x1+(2+s)*dx/4,y1+(2-s)*dy/4);

}

}

function go()

{

ctx.clearRect(0,0,canvas.width,canvas.height);

ctx.lineWidth = 2;

ctx.strokeStyle = "red";

ctx.beginPath();

ctx.moveTo(50+400/Math.pow(2,depth+1),50+400/Math.pow(2,depth+1));

drawShapes(depth,1,50,50,450,450);

ctx.stroke();

depth++;

if (depth>6)

{

depth=1;

}

}

window.setInterval('go()', 1000);

</script>

</body>

</html>

在浏览器中打开包含这段HTML代码的html文件,可以看到在浏览器窗口中呈现出如图9所示的Hilbert曲线动态生成效果。

图9  Hilbert曲线动态生成

JavaScript图形实例:Hilbert曲线的更多相关文章

  1. JavaScript图形实例:线段构图

    在“JavaScript图形实例:四瓣花型图案”和“JavaScript图形实例:蝴蝶结图案”中,我们绘制图形时,主要采用的方法是先根据给定的曲线参数方程计算出两点坐标,然后将两点用线段连接起来,线段 ...

  2. JavaScript动画实例:曲线的绘制

    在“JavaScript图形实例:曲线方程”一文中,我们给出了15个曲线方程绘制图形的实例.这些曲线都是根据其曲线方程,在[0,2π]区间取一系列角度值,根据给定角度值计算对应的各点坐标,然后在计算出 ...

  3. JavaScript图形实例:再谈IFS生成图形

    在“JavaScript图形实例:迭代函数系统生成图形”一文中,我们介绍了采用迭代函数系统(Iterated Function System,IFS)创建分形图案的一些实例.在该文中,仿射变换函数W的 ...

  4. JavaScript图形实例:随机SierPinski三角形

    在“JavaScript图形实例:SierPinski三角形”中,我们介绍了SierPinski三角形的基本绘制方法,在“JavaScript图形实例:迭代函数系统生成图形”一文中,介绍了采用IFS方 ...

  5. JavaScript图形实例:Koch曲线

    Koch曲线的构造过程是:取一条长度为L0的直线段,将其三等分,保留两端的线段,将中间的一段改换成夹角为60度的两个等长直线:再将长度为L0/3的4个直线段分别进行三等分,并将它们中间的一段均改换成夹 ...

  6. JavaScript图形实例:四瓣花型图案

    设有坐标计算公式如下: X=L*(1+SIN(4α))*COS(α) Y=L*(1+SIN(4α))*SIN(α) 用循环依次取α值为0~2π,计算出X和Y,在canvas画布中对坐标位置(X,Y)描 ...

  7. JavaScript图形实例:图形的旋转变换

    旋转变换:图形上的各点绕一固定点沿圆周路径作转动称为旋转变换.可用旋转角表示旋转量的大小. 旋转变换通常约定以逆时针方向为正方向.最简单的旋转变换是以坐标原点(0,0)为旋转中心,这时,平面上一点P( ...

  8. JavaScript图形实例:圆内螺线

    数学中有各式各样富含诗意的曲线,螺旋线就是其中比较特别的一类.螺旋线这个名词来源于希腊文,它的原意是“旋卷”或“缠卷”.例如,平面螺旋线便是以一个固定点开始向外逐圈旋绕而形成的曲线. 阿基米德螺线和黄 ...

  9. JavaScript图形实例:正多边形

    圆心位于坐标原点,半径为R的圆的参数方程为 X=R*COS(θ) Y=R*SIN(θ) 在圆上取N个等分点,将这N个点首尾连接N条边,可以得到一个正N边形. 1.正多边形阵列 构造一个8行8列的正N( ...

随机推荐

  1. Linux文件处理命令 ls 详解

    Linux系统的应用场景最多的就是用作服务器的系统了,简洁,安全,高效,一般我们服务器端不会安装Linux的图形化界面,虽然现在一些Linux发行版的图形界面也很漂亮,但是,服务器最主要的是高效.所以 ...

  2. 减少if...的使用

    最近维护一批代码,其中包括一堆if...的使用,多的情况嵌套8.9层,痛苦不堪,所以搜寻一些可以降低if...else的方法来改善一下代码,写个简单总结. 第一种: 优化前 if (measuredV ...

  3. Git Gui工具从远程克隆代码总是提示路径已经存在。问题完美解决!

    最近使用Git Gui工具从远程克隆代码总是提示路径已经存在,提示如下图, 原因是使用Git Guii克隆代码时的路径不能是选择一个已经建好的目录(D盘的project下的test文件夹是我提前建好的 ...

  4. 总结梳理:webpack中如何使用vue

    1. 安装vue的包 cnpm i vue -S  2. 由于在webpack中,推荐使用 .vue这个组件模板文件定义的组件,所以,需要安装,   能解析这个文件的loader: cnpm i vu ...

  5. 源码分析(1)-HashMap(JDK1.8)

    UML类图 java.util.Map<K, V>接口,有4个实现类:HashMap.Hashtable.LinkedHashMap和TreeMap. 1.说明 (1)HashMap除允许 ...

  6. 河青的持久层框架hqbatis

    谈到对数据库的操作,powerbuilder 的嵌入式SQL还是最方便的,增.删.改.查都无比的方便,可惜它落败于BS架构的盛起.java 以mvc的框架,实现对数据库的操作,写起来是相当麻烦,jav ...

  7. 【Java Spring Cloud 实战之路】添加一个SpringBootAdmin监控

    0. 前言 在之前的几章中,我们先搭建了一个项目骨架,又搭建了一个使用nacos的gateway网关项目,网关项目中并没有配置太多的东西.现在我们就接着搭建在Spring Cloud 微服务中另一个重 ...

  8. javascript 面向对象学习(三)——this,bind、apply 和 call

    this 是 js 里绕不开的话题,也是非常容易混淆的概念,今天试着把它理一理. this 在非严格模式下,总是指向一个对象,在严格模式下可以是任意值,本文仅考虑非严格模式.记住它总是指向一个对象对于 ...

  9. [TopCoder]Seatfriends

    题目   点这里看题目. 分析   可以想到用 DP 解决.   由于把空位放到状态里面太麻烦了,因此我们单独将 " 组 " 提出来进行 DP .   \(f(i,j)\):前\( ...

  10. GCDLCM 【米勒_拉宾素数检验 (判断大素数)】

    GCDLCM 题目链接(点击) 题目描述 In FZU ACM team, BroterJ and Silchen are good friends, and they often play some ...