一类求积公式的构造方法.docVIP

一类求积公式的构造方法 一类求积公式的构造方法 第16卷第6期 2000年12月 工科数学 JOURNALOFMATHEMATICSFORTECHNOLOGY V0I.16.?6 Dec.2000 @ I一l 一 类求积公式的构造方法 车翔玖,刘明才,jL (1北华大学师范理学院数学系,吉林13201a四平师范学院计算机科学系,四乎 [摘要]奉文利用一求和公式构造了一类求积公式,称为修正复合梯形公式.它和复合 梯形公式的求积节点及计算量是一样的,但收敛阶有很大的提高,特别适合于计算带有各种类型小放的数值 积分. [关t调]求积公式小波分析收敛阶--._.-..?_,?--_____一 [中圈分娄号].4[文献标识码]A l引言 数值积分在实际应用及科学计算中经常遇到,如复合梯形公式和复合公式,方法衙单又 具有递推性,但收敛阶较低.型公式虽然代数精确度高,但没有递推性,特别是在计算带有小波函 数的积分时,有时很困难,因为有些小波函数只能求出在二分点,^,召)上的值,求任意点的值比 较困难,如小渡等.文献[1]利用多项式构造出一种求积公式,获得了较高的精 度.但该公式的构造相当复杂,对不同的小波,要构造不同的公式,公式通用性受到 限制.本文利用一 一 求和公式构造了一类等距节点求积公式,称为修复合梯形公式,它和复合梯形公式 的 求积节点及计算量是一样的,但收敛阶有很大的提高.特别适合于计算带有各种类 型小渡的数值积分. 构造修正复合梯形公式需要用到一求和公式,现介绍如下: 设,)?c舯..阳,,…,^=气,?2五,Bz为rn.ulIi数,即日一吉,占=一,B一 ,s一 盍...?记 s()=(一1)(2)., 则有下面的.求和公式: 』:,()d=++E?(1.1) 其中 = 告[,(4)+,(6)+2?,+)],(1.2) P茬赘^[,l】(n卜.)],(1_3) E一n[,1)(n卜I)(.一(詈),(1_4) 脯式 州I?…一妻南l+a【…,(州I?…2南,(1_5) [收稿日期]-- 二h 拓 6 正m 名 积 44I科数学第16卷 其中Mzso ~20If()I,/(.是Riemann级数? 定义1.1若对一种复合求积公式J,当^一0时,有—c(c为不等于零的常数),则称J 的收 敛阶是O(h). 当E知L+P的收敛阶是O(h),即(1.1)式可写为 r,()d一十+()(^+).(1.6) 我们的基本思想就是通过已知点的函数值给出P中导数值的近似计算,并保证收 敛阶仍为O(h). 2求积公式的构造 对,()?CEa,6],由Taylor公式有 f(x+ih)--fc耋c, 其中i=1,2,…,2k,?b)?+m.记 c:=一!?:兰,E—csc,(,…,sc, 一 ,)^,y一(】,Y2,…,蹒),一,(x+ih)--f(x), F一(,l,f2,…,^),A一(口).i,一1,2,…,2h,4一 则(2.1)式可写为 AY=F+^E.(2.2) 命囊1.1(2.2)式中的矩阵A是非奇异的. 证将A的第i行(=1,2,…,2)除以f,将第列(』=1,2,…,2)乘以』,即得Vandermode行 列 式所对应的矩阵.从而A的逆存在. 由(2.2)式得 r—A一,十^A一E(2.3) 记A一(),则(2.3)式为 ,()一?6Ef(x+jh)一f(x)Je-h?6,,(),=12,2h.(2-4) 取一4,有 2^ (口)?bEf(a+jh)一f(a)Jq-h?b^(4)r1=1 (2.5) 取x=b,h为一^,(此时b--ih?gab)?6),有 (6)(--;?6,厂(6一ja)一,(6)]一h?6.,,().(2.6) 将(2.5)和(2.6)代入(1,3)中的,有 ^ 妾{善%Ef(卅)卅m一)_,((EJ(+1]}.(2.7) 记 又记 ,Lf(口+)一,(口)q-f(b—ih)一,(6)]) 耋{f[E,](2.8) (2.9) 6? 一?. ? ^ 第6期车翔玖等:一类求积公式的构造方法45 7=7+户.(2.10) 则(1.1)可写为 lf(x)dx=一邑+E,(2.11) 其中丁就是我们所要构造的修正复合梯形公式.而就是对(1.1)式中的导数值用已 知点的函数 值代替所产生的误差.下面给出(2.II)中的邑一的估计式.由(14)和(2.9),可得 ; ?[一毒l川 +^Ma2Or)一.(6一n)721,(2.12) 其中M2蚪l—maxI()I. 由此可得 定理l_2设,b)EC肚*],则复台修正梯形公式的收敛阶是()(^坩).即,()d 了