特征值问题的预变换方法Ⅱ任意三角形域Laplace特征值的计算分析.pdfVIP

2012年2月 计算数学 第34卷第1期 F曲．2012 MATHEMATICA NUMERICASINICA Vol34．Nol 特征值问题的预变换方法(II)： 孙家昶 (中国科学院软件研究所并行计算实验室，北京100190) 摘 要 系的构造，提出任意三角形区域上LapIace特征值的近似公式与算法给出任意三角形域上所 其中s是该三角形(^1Sk茎b)的面积，可作为数值PDE中三角剖分质量的一种新标准 q(7n：=蒜里垡二：；；互±蛆结合数值计算与符号计算，将这三类三角形的基底综合形成统一的 新基底，以反映几何(三条边)对于特征问题的影响，从而提高任意三角形域的求解精度 关键词：特征值问题的预变换方法，Laplace特征值问题，任意三角形域 MR(2000)主题分类：35J05，35P99，65N25，65N30 1．引 言 L印lace算子的特征值问题的估计与分析是经典数学的一个课题，很多大数学家，包括R． 征函数系，从而给出相应的L印lace特征值的精确初等解或近似表示、是计算数学与科学工程 年来，再次出现了关于等边三角形域特征值的极小性质的文献flo．131． 我们在文【7]中指出任意三角形域上的函数在重心坐标下可分解为各向同性、偏反对称及 偏对称三部分并采用三次B一样条计算过L印lace特征值稍后作者在专著㈣第二章系统研 究了正三角形域零法向复值特征函数系与零边界复值特征函数系的构造以及相应高频函数与 近发表的工作口…强调基底选取结合采用预变换方法对于高效求解特征值问题的重要性． 基于我们前面的研究成果，求解任意三角形域特征问题时，本文考虑以下三个方案一是 采用传统的单一基底，如等边三角形正交特征函数系，二是将以上等边三角形、等腰直角三角 形及(300，600，90。)三角形这三类三角形的基底综合形成统一的新基底，代替只对一类三角形 精确的基底，以更好反映任意三角形中三条边长对于特征值的影响，第三则是通过这三组精确 特征值通过插值给出任意三角形区域上Laplace特征值的近似解析解 本文节2叙述三种特殊三角形区域上L印lace特征值问题的精确解，节3提出任意三角 +2011年9月15日收到 1)基金项目：国家自然基金项目资助(N61170075) 万方数据 。 二—————————j坐塑兰 婴! 征值的一类近似解析解，最后，节5给出部分计算试验结果 2·三种特殊三角形区域上Laplace特征值问题的精确解 ^}=(可z∥+(z2一z3)2，^；=‰一∥+(蜘一z，)t，壕=(Ⅳ， 掰+(。。一。2)。 引人相应的重心坐标(“，”，w)为 “：去睢斟。：去|}锄一去帷珥㈦。， 11 23 掣3【 l z3 耵3l fl 。 可| 这时7-2((“，”，”)：o茎Ⅱ，，w；“+u+Ⅲ=1)记s为三角形7-的标量面积， 肛嘉m，n‰=器m，叫，氩=蒜m，乩 易知L印】ace算予在(％口)坐标下的表示 △，(“，u)2j刍{hi^n+h；，t，。+(^；～^：一^!)^．。)(221 一般，不妨设^lS^2 1b，引入双参数 “1 2“2订i r％=^2厩，(o!El，旬1)f23) 由三角不等式o s乜一^l^2对应的双参数(sl，旬J取值范围为(见示意图1)