數(shù)值求積公式及代數(shù)精度數(shù)值求導方法與截斷誤差一階常微分方程數(shù)值法局部截斷誤差與精度《數(shù)值分析》典型例題

IV典型例題與習題4插值型求積公式:求積系數(shù)求積余項等距結點插值型求積公式稱為Newton-Cotes公式,偶數(shù)階Newton-Cotes公式至少有(n+1)階代數(shù)精度2/16求積結點典型例題與習題41.梯形公式復合梯形求積公式

令h=(b-a)/n求積余項3/162.辛卜生公式典型例題與習題4求積余項:

兩點高斯型數(shù)值求積公式4/16練習:復合辛卜生公式求積余項?典型例題與習題4一階向前差商一階向后差商二階中心差商一階中心差商5/16典型例題與習題4外推算法練習：二階中心差商的外推公式？6/16典型例題與習題41.Euler方法

常微分方程初值問題2.梯形公式:

7/16典型例題與習題4預測-校正公式局部截斷誤差設

yn=y(xn),稱Rn+1=y(xn+1)-yn+1為局部截斷誤差常表示為:O(hp+1)，

p稱為單步法的精度階數(shù)又稱為修正的Euler公式

yn+1=yn+0.5h[k1+k2]k1=f(xn,yn),k2=f(xn+h,yn+hk1)8/16典型例題與習題4Ex1.推導左矩形求積公式

令F(u)=F(a)+(u-a)F’(a)+0.5(u-a)2F”()練習:9/16典型例題與習題4Ex3.求復合中矩形公式的求積誤差?Ex2.復合左矩形求積公式的求積誤差設被積函數(shù)在積分區(qū)間上的一階導數(shù)連續(xù),由連續(xù)函數(shù)介值定理10/16典型例題與習題4Ex4.利用復合梯形公式計算積分使其截斷誤差不超過

0.5×10-3,應算多少次函數(shù)值？提示：練習:

給定積分當要求誤差小于10-3時用復合梯形公式和Simpson公式計算時,需要計算多少次函數(shù)值？11/16典型例題與習題4Ex5.驗證，復合梯形公式與復合Simpson

公式之間有如下關系12/16典型例題與習題4Ex6.定積分的計算問題可化為初值問題

y’=f(t),y(a)=0試證明用Euler公式計算結果為其中,h=(b–a)/N,

tn=a+nh(n=0,1,2,···,N)Ex7.試證明4階Range-Kutta公式解[a,b]內初值問題

y’=f(x),y(a)=0結果有:其中,h=(b–a)/N,

xn=a+nh(n=0,1,2,···,N)13/16典型例題與習題4Ex8.將線性常系數(shù)非齊次高階常微分方程初值問題:y(n)+a1

y(n-1)+a2y(n-2)+·······+any=f(x,y,····,y(n-1))

y(x0)=y00,y’(x0)=y01,y”(x0)=y02,···y(n-1)(x0)=y0,n-1轉化為一階線性常微分方程組問題,并成出矩陣形式解:令

y1(x)=y(x),y2(x)=y’(x),y3(x)=y”(x),····,yn(x)=y(n-1)(x),

14/16典型例題與習題4Ex9.

初值問題有解y(x)=0.5ax2+bx

。若取

xn=nh，yn為歐拉方法得到的數(shù)值解，試證明y(xn)–yn=0.5ahxn若取

xn=nh，yn為用梯形公式計算所得的數(shù)值解，記y(xn)為初值問題的在x=xn處的解析解。試證明:

y(xn)=yn

Ex10.

初值問題15/16典型例題與習題4寫出梯形