東南大學(xué)數(shù)值分析上機(jī)題答案說(shuō)課講解-2023修改整理

千里之行，始于足下讓知識(shí)帶有溫度。第第2頁(yè)/共2頁(yè)精品文檔推薦東南大學(xué)數(shù)值分析上機(jī)題答案說(shuō)課講解東南高校數(shù)值分析上

機(jī)題答案

數(shù)值分析上機(jī)題

第一章

17.（上機(jī)題）舍入誤差與有效數(shù)設(shè)∑=-=

N

jNjS2

2

11

，其精確值為）111-23（21+-NN。（1）編制按從大到小的挨次1

-1

···1-311-212

22NSN+++=，計(jì)算NS的通用程序；

（2）編制按從小到大的挨次1

21

···1)1(111

222

-++--+-=NNSN，計(jì)

算NS的通用程序；

（3）按兩種挨次分離計(jì)算210S,410S,610S,并指出有效位數(shù)（編制程序時(shí)用單精度）；

（4）通過(guò)本上機(jī)題，你明了了什么？解：程序：

（1）從大到小的挨次計(jì)算1

-1

···1-311-212

22NSN+++=

：functionsn1=fromlarge(n)%從大到小計(jì)算sn1formatlong;

sn1=single(0);form=2:1:n

sn1=sn1+1/(m^2-1);endend

（2）從小到大計(jì)算1

21

···1)1(111

222-++--+

-=

N

NSN

functionsn2=fromsmall(n)%從小到大計(jì)算sn2

formatlong;

sn2=single(0);

form=n:-1:2

sn2=sn2+1/(m^2-1);

end

end

（3）

總的編程程序?yàn)椋?/p>

functionp203()

clearall

formatlong;

n=input('pleaseenteranumberasthen:')sn=1/2*(3/2-1/n-1/(n+1));%精確值為sn

fprintf('精確值為%f\n',sn);

sn1=fromlarge(n);

fprintf('從大到小計(jì)算的值為%f\n',sn1);

sn2=fromsmall(n);

fprintf('從小到大計(jì)算的值為%f\n',sn2);

functionsn1=fromlarge(n)%從大到小計(jì)算sn1

formatlong;

sn1=single(0);

form=2:1:n

sn1=sn1+1/(m^2-1);

end

end

functionsn2=fromsmall(n)%從小到大計(jì)算sn2

formatlong;

sn2=single(0);

form=n:-1:2

sn2=sn2+1/(m^2-1);

end

end

end

運(yùn)行結(jié)果：

從而可以得到

N值真值挨次值有效位數(shù)2

100.740050從大到小0.7400495

從小到大0.74005064

100.749900從大到小0.7498523

從小到大0.74990066

100.749999從大到小0.7498523

從小到大0.7499996（4）感想：

通過(guò)本上機(jī)題，我明了了，從小到大計(jì)算數(shù)值的精確位數(shù)比較高而且與真值較為臨近，而從大到小計(jì)算數(shù)值的精確位數(shù)比較低。機(jī)器數(shù)在舉行加法運(yùn)算時(shí)，用從大到小的挨次簡(jiǎn)單浮現(xiàn)大數(shù)吃小數(shù)的狀況，簡(jiǎn)單產(chǎn)生較大的誤差，是由于對(duì)于相加的兩個(gè)數(shù)值，計(jì)算機(jī)首先提供與大數(shù)相全都的位數(shù)，此時(shí)將小數(shù)的尾數(shù)向右移位，并舉行四舍五入，之后對(duì)尾數(shù)舉行依次相加。從大到小時(shí)，越往后計(jì)算，相加的數(shù)越小，從而浮現(xiàn)大數(shù)吃小數(shù)的狀況。相比之下。從小到大計(jì)算時(shí)，每次小數(shù)與大數(shù)相加，都會(huì)增強(qiáng)位數(shù)，從而精確度比較高。

其次章

20.（上機(jī)題）Newton迭代法

（1）給定初值0x及容許誤差ε，編制Newton法解方程)(xf=0根的通用程序。

（2）給定方程xxxf-=3/)(3，易知其有三個(gè)根*1x=-3，*2x=0,*

3x=3。

①有Newton辦法的局部收斂性可知存在δ>0，當(dāng)0x∈（-δ,δ）時(shí)Newton迭代

序列收斂于根*2x，試確定盡可能大的δ；

②試取若干初始值，觀看當(dāng)0x∈（-∞，-1），（-1，-δ），（-δ，δ），（δ，1），（1，+∞）時(shí)Newton序列是否收斂以及收斂于哪一個(gè)根。（3）通過(guò)本上機(jī)題，你明了了什么？解：（1）程序

先編寫(xiě)函數(shù)function文件：文件fx.m

%定義函數(shù)f(x)

functionFx=fx(x)

Fx=x^3/3-x;

文件dfx.m

%定義導(dǎo)函數(shù)df(x)%

functionfx=dfx(x)

fx=x^2-1;

接下來(lái)是詳細(xì)步驟

文件newton1.m求盡可能大的delta值

%%課本56頁(yè)計(jì)算最大delta值

clear

flag=1;

k=1;

x0=0;

whileflag==1

delta=k*10^-6;

x0=delta;

k=k+1;

m=0;

flag1=1;

whileflag1==1

ifabs(x1-x0)=10^-6

flag=0;

end

end

fprintf('%f\n',delta);

文件newton2.m求方程的根

%%課本56頁(yè)newton法求方程的根，確定收斂于哪個(gè)根formatlong;

ef=1e-6;

k=0;

x0=input('pleaseentertheinitialnumberasthex0:');whilekk

A([p,k],:)=A([k,p],:);b([p,k],:)=b([k,p],:);end

m=A(k+