投影型插值的新的特殊性质及其在高次有限元中的应用.docxVIP

湖南师范大学2005届硕士学位论文 ·3 从(1 1)来看，似乎插值多项式的次数越高，逼近精度都会越好． 但事实上并非如此。这是因为随着插值节点的加密(即≈的增Jon)，插值多项式在某些非节点处的振荡可能加大即Ju(¨，iff)I之值有可能以极高的速度增加，从雨hu(x)来必能一致地逼近函数“(z)，这种现象称为高次插值Runge现象．所以在实际数值计算中很少使用高次 Lagcrange型元。 本文研究的这种新型投影型插值与Lagrange插值相比有以下特 点： Iku(a)=iku(a)=¨(n)，Iku(b)=iku(b)=u(6j{因而便于作复合摇 值(分段☆次插值)． 2它ff．j都是区阊a，67上的七次多项式插值。 3对Lagrange插值只能有误差估计(1．1)，而对于投影型插值，本 文将证得如下逐点误差估计： №)一歃“(z)l≤、／A(“。；毒‰t“一x‰6J ≤~／(6一o)凡(。)Xfi一，1¨Ⅱ：，岛‰hi)l“一xhhq 其中月k(z)为2k次非负多项式，满足 x—a≥A1(￡)≥A2(x)2 A2(x)三···≥Ak(x)一0，比∈[a，b】 也就是说，只要u∈Hl([a，b1)，就可以保证iku(x)在陋，6】上一致收敛于u(z)，不存在Runge现象． 4高次Lagrange插值多项式^u的基函数复杂不容易进行计算， 而投影型插值iku(x)的基函数咄十分简单，易于计算． 5从理论上讲，投影型插值可以广泛地应用于高次有限元的计 算，获得高效计算结果，而Lagraa、gc插值贝Ⅱ不能． 本文在论证中使用通常的Sobolev空间范数符号： W”9(n)={：D。u∈L’(Q)，Ⅱ≤m) D“”表示”的a次导数，Sobolev空间彤一(固的范数¨}Im,p,n、半 范数I．I。艄可简记为”JJ。、j-1。，当P=2时，范数、半范数记 湖南师范大学2005届硕士学位论文 为川I。，n、l·l，n，n也可简记为川I。、}·I。．用最(e)表示定义在e上 的所有％次多项式． 湖南师范大学2005届硕士学位论文 r5 第二章Legendre多项式和Lobatto多项式 本文讨论的Lobatto展开及投影型插值是基于Legendre多项式和Lobatto多项式这两类正交多项式系，本节将具体地介绍。 假定一是区间b嘲上的一个拟一致剖分，％是该剖分的k次 有限元空间，设一维单元e=k—he，z。+^。】毫In，6】，{fJ(z))‰为工2(e) 上规范完备的Legendre正交多项式系． 20(z)=、／{K“2 f。(z)=~／弦3肛(z—z。) f。(z)={、／i^i副2【3扛一$。)2一瑚 b(z)=ji V1j7，k-772【5(z一￡。)3—3^：(。一∞。)】 f4睁)=；、／2，ai9舟135@一z。)4—30hi(x一。。)2+3^4] f5(z)=}√竽^i“72[63(。一z。)5—70^：扛一z。)3+15，瑶(∞一。。)】 f6(∞)=击、／萼^i”胆[231(z—z。)6—315h：(x—z。)4+105h：(x—z。)2—5^!] f7(z)：击√警K15∞[429(z一)7—693(∞一)5瑾+315(x一。。)3^4。一35(霉一。。)磋 f，(。)=以(鑫)’【A(z)r，(i≥1)， 以=、／堕≯扣“叫2=O(hc-。V2) (2．1) 其中A(z)={【扛一z。)2一姻．又考虑铲(e)上的规范正交多项式系：令 苷嚣!／Y翟意=瑚 ：三舞裂娑≯ =；、／i酊’[@一叫3一h：如一ze)】 ：!： 塑壹坚耋查童!!塑鱼塑圭兰竺篁圭 u4=j1．v／。7h。-1／ [5(z—xo)3—3h：(。一z。)]如 =ilvj7‰--5f5(z一％)4—6h；(x一％)2十醒j u5=。1 v／。9h。-i／ [35(x—z。)4—30hm(=一。。)2+3h：]dx =i1 V i9，。c--5【7(茹一z。)5—10，瑶(z一。。)3+3h：(x一‰)] 旷{饵巧警／。。[63(x-x,)s一70h：(x飞)s+15蠼(z-x,)]dx 一去√警hii[21(z一。。)6—35h；(z—z。)4+15^：(。一。。)2一磋] u，=Z二。，去、／萼‘警∞l扛一㈦6～315扛一z。)t醒2+1。5扛一嘲z毽一5^6】出 =去~／竽h2警【33扛一。。)7—63(x～孔)5蜒2+35(x一岛)3嘭+5扛一)蠼】 us=Z二k去、／萼【429扣一引7_693@一孔)5醒+315@一。c)3磋～35@一叫确如 =南~／萼【429扛一％)8～924(x一％)6硅+630(x一％)4蠼4—140(x一。。)2磋+5h：1 “t+l=／，li(z)ax=以(杀)。1陋(z)h(i≥1) d2c—n