积分域封闭的二元GAUSS数值积分方法研究.pdfVIP

第 28卷 第 3期 西华 师范 大 学学报 (自然科 学版 ) 2007年9月 Vol.28 刊住3 」皿顶目ofChinawestNO们malUni7e几it了(Naturalsciences) Sep.2007 文章绷号:1673巧0，2(2007)03刁2刁9拼 积分域封闭的二元GAUSS数值积分方法研究 谭代伦，，张世禄’，体 卫2 《1西华师范大学致学与信息学院 四川南充 637训2:2.内江师范学院孜学系.四J”内江 641.心刀) 摘 要:从构建二元肠，朋，播值多项式出发，给出了积分城封闭的具有较高代数精确度的二元cAUss数值积 分公式，其代数精确度的定义和证明完善，实例计算结果良好 关键词:二元L叱Za叱c插值多项式;二元GAuss数值积分;代数精确度 中图分类号:0241，4 文献标识码:A 引 言 一元数值积分的理论和方法在计算方法中已得到较为彻底的解决，二元及以上的多元数值积分的理论 和方法则还在不断发展和完善中，诸如乘积方法、单项式方法、抽样法、点阵法等川.归结起来，通用性较高 的二元及以上数值积分算法在理论和方法上仍较少，而其实际应用意义比一元数值积分大，因此给出二元及 以上数值积分公式就具有较高实用的价值. 在一元数值积分中，和Ne们on一cot，积分一样，GAuss积分属于给定求积节点的插值型求积方法， 且都是用Lag姗ge插值多项式代替被积函数的一类数值算法，其区别在于Ne吻n一Cote日积分用等距节点的 代数插值多项式，CAUss积分的节点选用 。+1阶正交多项式零点.基于稳定性、收敛性等原因.很少采 用高阶Newlon一cote。积分，而是采用GAuss型求积公式或复合求积公式以进一步提高数值积分的代数精 确度阁.因此下面从构建二元Lagrang。插值多项式出发.然后讨论二元cAuss数值积分公式及其代数精 确度. 2 二元GAuss数值积分公式及代数精确度 2.1构建二元L旧侣range插值多项式 对于连续函数:二f(二，力.二。[a二，二‘〕，y。。〔，，b，〕，取节点(二，为‘》，二‘。.1，…，。，1二。，1，·，n，当节 点互异时，有 L一(，，·)=弧‘1.“，()‘“.，(，)八，，若)· (1) 其中1.“，(，)二只我，二‘”，，‘一 ”‘，)‘二众狱 j=。，1，·，枯.是关于x为胜阶、y为几阶 的二元插值多项式. 证明 任给点(二，，九)，p=0，1，·，。，，二0，1二，n，由孟护，(二)和叮”(了)有 L一、〔，，。)=毓 :司‘(、)心，·(y，)f(!。，:卜“，·(、)“，·(、)八、，y，)=八!，.、)· 即在节点处插值多项式的值是等于函数的值. 类似地，可定义二元以上的擂值多项式，因上述插值多项式和一元La脚n，插值多项式类似，故称为二 收摘 日期:2007一供 一叮 作者简介:谭代伦(l97一)，男 四川南充人，西华师范大学数学与信息学院讲师，主要从事应用数学与计算机应用方面 的教学与研究工作. 万方数据 250 西华师范大学学报 〔自然科学版 》 2007年 元Lagran郎插值多项式. 2.2给出二元GAuss数值积分公式 对函数:二f(二，了)，二。a〔.，b二」，李。a〔，，b，〕，取节点二j为二十1阶肠gran驴插值多项式的第(‘二‘。， 1，·，司个零点击为.+1阶肠乎an邵插值多项式的第j(j二。.1，…，n)个零点，则有 二“以L一，y·‘，d琦二铸圈*j(气，yj，气铸:一，f(x，yj)· (2) 上式称为二元GAuss数值积分公式，式中9尸和司.)分别为。阶和。阶GAuss系数的第(‘二‘。，1，·，。) 个和第j(j二。，1，二的个系数，9汁，为二元GAuss积分系数， 当二‘[a.，，‘]、y。a〔，，，‘]时.利用坐标变换可将(2)式变成 S二毕 导滋‘卿争 毕·‘，导芳+毕，· (3) 同样，也可给出二元以上的C人USs数值积分公式. 2.3 代数精度的