三维空间中二维网格的迭代结果不是环面.docxVIP

三维空间中二维网格的迭代结果不是环面 ?　　昨天，张文琦(228254913) 同学的一篇 HYPERLINK /9Qxbf4 \t _blank 日志又把椭圆迭代问题推向了高潮。这篇日志把平面中的多边形推广到了三维空间中的二维网格。作者得到的结论是，二维网格的迭代结果是环面（torus）。同时作者注意到，有时候环面上会打一个“结”，这个结需要大量的迭代才能解开。 　　我对这个“结”感到非常好奇，于是继续用上一篇日志中的傅里叶分析方法对二维网格的迭代进行了分析。结果发现：这样的“结”其实有两个；它们并不只是“有时”会出现，而是必然会出现；而且，无论迭代多少次，这两个“结”都不会消失。二维网格的迭代结果并不是环面，而是一种叫作elliptic translational celestial的曲面，它与环面最显著的区别就是有两个“结”。 一、三维空间中二维网格迭代问题的描述 　　与上一篇日志一样，我用小括号表示迭代次数，中括号表示点的下标。 　　首先定义二维网格。在三维空间中取n*m个点P[i][j] = (x[i][j], y[i][j], z[i][j])，其中0=in，0=jm。另设P[n][j] = P[0][j]，P[i][m] = P[i][0]。将所有一个下标相同、另一个下标相差1的点对相连，就组成了一个网格。这个网格的每一个“孔”是一个（空间）四边形，四边形[i][j]的四个顶点为P[i][j]、P[i+1][j]、P[i+1][j+1]、P[i][j+1]。 　　迭代过程如下：首先取定初始网格的各个顶点P(0)[i][j]，并平移使得各个顶点的平均值位于原点。每次迭代时，取各个四边形的中心为新网格的顶点，即 P(k)[i][j] = ( P[i][j] + P[i+1][j] + P[i+1][j+1] + P[i][j+1] ) / 4，然后缩放使得各点的x、y、z坐标的平方和均为1。 　　问：二维网格经过上述迭代过程，收敛成什么形状？ 二、对二维网格形状的朴素思考 　　读到这里，你脑中想象出来的二维网格也许是分布在一个球面上的。这很正常，因为日常生活中最典型的网格就是地球上的经纬线。但不幸的是，地球的经纬线并不能组成上面定义的二维网格。因为按上面的定义，在网格上，往上下左右四个方向都能够无限地前进，而地球却有南北两极。 　　“吃豆子”游戏（Pacman）的地图符合上面的定义，因为从地图的左（右）边出去可以从右（左）边回来，从上（下）边出去可以从下（上）边回来。你可曾想过，这样的地图怎样才能实现？答案是环面： 　　张文琦同学用欧拉示性数证明了满足上面定义的二维网格所“撑”起来的曲面，必定同胚于环面。通过计算机绘图观察，他认为二维网格经过迭代也收敛成环面。 三、用傅里叶变换分析迭代结果的形状 　　我用Matlab绘制得到的一次迭代结果如下图所示。我们看到，图形打“结”了，而且这个“结”经过了上千次迭代都没有消失。 　　怎么回事呢？我决定再次拿出傅里叶变换这个大杀器。与多边形的情形一样，二维网格的迭代关于x、y、z三个坐标也是独立的。所以，我们可以只分析x坐标。 　　在分析多边形时，我们把各个点的x坐标看成了一个离散的、周期性重复的一维序列。在这里，我们可以把二维网格上各点的x坐标看成一个离散的、周期性重复的二维序列。而迭代操作，则是把这个二维序列与二维滤波器[1/4, 1/4; 1/4, 1/4]进行二维循环卷积。 　　这个滤波器的二维傅里叶变换的模如下： 除去直流分量（左上角），模最大的两个值位于第2行第1列（以及左下角）和第1行第2列（以及右上角）。迭代过程中能留下来的也正是这两个分量，它们在原域中分别是sin(2*pi*i/n + P1)以及sin(2*pi*j/m + P2)（P1、P2为任意相位）。于是我们可以得到迭代足够多次后的网格表达式： 　　x(k)[i][j] ≈ Cx1(k) * sin(2*pi*i/n + Px1(k)) + Cx2(k) * sin(2*pi*j/m + Px2(k))　　y(k)[i][j] ≈ Cy1(k) * sin(2*pi*i/n + Py1(k)) + Cy2(k) * sin(2*pi*j/m + Py2(k))　　z(k)[i][j] ≈ Cz1(k) * sin(2*pi*i/n + Pz1(k)) + Cz2(k) * sin(2*pi*j/m + Pz2(k)) 其中Cx1、Cy1、Cz1、Cx2、Cy2、Cz2、Px1、Py1、Pz1、Px2、Py2、Pz2由初始网格决定，在迭代过程中按滤波器的特性演变。 　　注意到上面三个式子中，如果仅保留右边的第一项，则可以得到一个三维空间中的椭圆，其中心位于原点；仅保留右边的第二项，也可以得到一个中