正三角形区域中一组完备正交基的性质与分类.docVIP

正三角形区域中一组完备正交基的性质与分类 杨志杰 (中科院软件所并行软件研究开发中心) 摘要 本文通过重心坐标法来研究正三角形区域中的一组完备正交基的性质与分类，并给出这组基底函数的可视化结果。 关键词 三角形 特征值 特征向量 正交基 1( 介绍 正交函数基底在数值分析中一直起着重要的作用。如何将一组正交基底推广到高维空间，仍然是一个尚待解决的问题。而传统的张量积方法只能适用于矩形区域，对于其他区域如三角形区域则不适用。 在一维情形下，通过求解刘微尔—斯托姆(Sturm-Liouville)问题可以得到一组完备正交基底。这种方法可以通过求解相应的偏微分问题来推广到高维。 在这篇文章中，我们考虑正三角形区域上带狄里赫勒(Dirichlet)边界条件的拉普拉斯(Laplace)算子的特征值问题: T3 ,,u,u,(x,y),,, { (1) u,0,(x,y),,, ,,{(x,y):0,y,3(x，1/3), T1 O T2 y,,3(x,1/3)} 图一 [1]1980年， Pinsky通过延拓法解出了这一特征值问题一组特征值和相应的特征函数， [3]并于1985年证明了这组特征函数的完备性。1998年， Prager又用相似的方法独立推导出 [6]了一组特征值和相应的特征函数，并证明了特征函数的完备性与正交性。本文从Prager的结果出发，首先给出正三角形区域中的一组完备正交基，然后通过重心坐标来研究这些特征函数的性质与分类，最后通过可视化来验证所得的结论。 2( Prager的研究结果 1998年，Prager通过延拓折叠法(prolongation and folding), 从矩形区域的特征函数导出了正三角区域的一组特征函数，并证明了这组特征函数构成了正三角区域的一组完备正交基， 结果如下: 定理1 函数u , k, l = 1, 2, ..., k , l mod 2, 0 k l, 和v , k = 0，1，2 ..., l = 1, k,lk,l 54 2, ... , k , l mod 2, 0 , k l, 组成上的一组完备正交体系, 其中: ,kx,u(x,y),2sinsinlyk,l3 ,,xx(k，l)/2 ,2(,1)sin(k，3l)sin(k,l)(2)223 ,x,x(k,l)/2，2(,1)sin(k,3l)sin(k，l)223 ,kx,v(x,y),2cossinlyk,l3 ,,xx(k，l)/2 ，2(,1)cos(k，3l)sin(k,l)(3)223 ,x,x(k,l)/2,2(,1)cos(3l,k)sin(k，l)223 1/4 并且，当k, l = 1, 2, ...时u和v 的模为3;v 的模为23。k,lk,l0,l 3( 基底函数的重心坐标表示以及它的几何性质 定义 平面上任意一点P可以用它对应于非退化三角形T(T, T, T)的重心坐标(, 1231 , )来表示，其中: 23 area(P,T,T)area(T,P,T)area(T,T,P)231312,,,,,,,, 123area(T,T,T)area(T,T,T)area(T,T,T)123123123 3 |,|:,,,1显然有 。 ,i1i, 由以上定义，可得到对应于第一节中的三角形(T, T, T)(见图一)的重心坐标与笛卡123 尔坐标之间的转换公式: 13y13y,,,x,,,,，x,,,,y(4) 123222222 将公式(4)分别代入公式(2)、(3)， 并化简，得到由重心坐标表示的基底如下: 55 ,kx,,,,,,,u(,,),2sin(,)sinl,123213kl3 ,kx,,,,，2sin(,，2)sinl(5)3213 ,kx，2sin,(,,,2)sinl,,1323 ,kx,,,,,,,v(,,),2cos(,)sinl,123213kl3 ,kx,,,,，2cos(,，2)sinl(6)3213 ,kx，2cos,(,,,2)sinl,,1323 由基底函数的重心坐标表示，通过简单的计算可以得到基底函数的一些几何性质: 性质一:函数u ,关于高O T反对称，函数v ,关于高O T对称。 k,l3k,l3 这由u (, , ) = -u (, , )，v (, , ) = v (, , k,l123k,l213k,l123k,l21 )可以看出。 3 性质二:若k , 0 mod 3， 则函数v是全对称函数，即函数v关于三条高都对称;k,lk,l 函数u是全反对称函数，即函数v关于三条高都反对称。 k,lk,l 证明:若k , 0 mod 3，则有 56 ,,,,,,,,kxu(,,),2sin(,)sinl,123213kl