本文格式為Word版，下載可任意編輯——XX工業(yè)大學數(shù)值分析試驗要求試驗一舍入誤差與數(shù)值穩(wěn)定性

（2學時）

1.試驗名稱2.試驗?zāi)康?.算法描述4.源程序5.運行結(jié)果6.對算法的理解與分析（包括改進與建議）

程序與實例

1例1對n=0,1,2,…,20計算定積分yn=

?x?5dx0xn

算法1利用遞推公式

1yn=

1n-5

yn?1n=1,2,…,20

取

y0??x?5dx01?ln6-ln5?0.182322

15n120算法2利用遞推公式y(tǒng)n?1注意到

1126?11201?1?15y1nn=20,19,…,1

?x60dx??x?50x20dx??x50dx?105

取

y20?)?0.00873020235126(?111

試驗二拉格朗日插值與牛頓插值

（4學時）

一、目的與要求：

?熟悉拉格朗日插值多項式和牛頓插值多項式，注意其不同特點；

二、試驗內(nèi)容：

?通過拉格朗日插值和牛頓插值多項式的兩個實例的計算，了解兩種求解方法，分析他們的優(yōu)缺點。

三、程序與實例

算法1．2．l

輸入xi,yi(i=0,1,2,?,n),令L(xn)=0;對=0,1,2,?,n計算

li(x)=?i?0j?inx?xjxi?xj

Ln?Ln+li(x)yi程序與實例例1已知函數(shù)表

xiyi

用三次拉格朗日多項式求x=0.5635的函數(shù)近似值。?牛頓插值多項式算法

1.輸入n,xi,yi(i=0,1,2?,n);

0.561600.827410.562800.826590.564010.825770.565210.824952.對k=1,2,3?,n,i=1,2,?,k計算各階差商f(x0,x1?,xk);3.計算函數(shù)值

Nn(x)=f(x0)+f[x0,x1](x-x0)+?+f[x0,x1,?,xn](x-x0)(x-x1)?(x-xn?1)程序與實例

例2已知函數(shù)表

xiyi

用牛頓插值多項式求Nn(0.596)和Nn(0.895)。

0.40.410750.550.578150.650.888110.80.91.026521.02652

試驗三復化辛卜生法，龍貝格法

（4學時）

一、目的與要求：

?通過實際計算體會各種方法的確切度；

?會編寫用復化辛卜生、龍貝格算法求定積分的程序。

二、試驗內(nèi)容：

?通過實際計算體會各種方法的確切度并且會編寫用復化辛卜

生、龍貝格算法求定積分的程序

三、程序與實例

復化辛卜生公式

算法：復化辛卜生公式為Sn=h/6??[f(xk)?4f(xkk?0n?1?h/2)?f(xk?1)],計

算過程為：1．令h?(b?a)/n,2．對k

?1,2,?,n?1

s1?f(a?h/2),

s2?0;

計算

。

s2?s2?f(a?kh)

s1?s1?f(a?kh?h/2),

3．s?h/6(f(a)?4s1?2s2程序與實例

?f(b))例用復化辛卜生法計算積分

運行結(jié)果為s(2)=0.785392s(4)=0.785398s(8)=0.785398

I??101/(1?x)dx2

說明：本例運行了三次，當n?2?8時，就與n?2?4時有6位數(shù)字一致，若用復化梯形法計算，當n=512時有此結(jié)果。

32?龍貝格算法計算I??1/(1?x012)dx

??5e?6

算法

n?1T2n?(Tn?hn?f(xk?1/2)k?0

hn?(b?a)/n;xk?1/2?(k?1/2)hn

Sn?T2n?1/3(T2n?Tn)Cn?S2n?1/15(S2n_Sn)Rn?C2n?1/63(C2n?Cn)

R2n?Rn|?5e-6。

用事后估計法控制精度

試驗四改進歐拉法，二分法，牛頓法

?（4學時）

一、目的與要求：

?熟悉求解常微分方程初值問題的有關(guān)方法和理論，主要是改進歐拉法

?會編制上述方法的計算程序

?

針對實習題編制程序，并上機計算其所需要的結(jié)果；

二、試驗內(nèi)容：

?熟悉求解常微分方程初值問題的有關(guān)方法和理論，主要是改進歐拉法,體會其解法的功

能。程序與實例

?改進歐拉方法算法概要

解一階常微分方程初值問題

?y??f(x,y)?y(x)?y00?

a?x?b

b?an將區(qū)間[a,b]作n等分，取步長h?歐拉公式為

yi?1?yi