最新东南大学数值分析上机题matlab(前三章) .pdfVIP

精品好文档，推荐学习交流

数值分析上机题

第一章（17题）

(1)从2依次累加到N的程序

functionsn=sum1(n)

sn=0;sn=single(sn);

fori=2:n

ai=1/(i^2-1);

sn=sn+ai;

end

end

（2）从N依次累加到2的程序

functionsn=sum2(n)

sn=0;sn=single(sn);

fori=n:-1:2

ai=1/(i^2-1);

sn=sn+ai;

end

end

（3）编制求精确值的求和函数sum0

functionsn=sum0(n)

sn=0;sn=single(sn);

sn=1/2*(3/2-1/n-1/(n+1));

end

按第一种顺序得到的值及有效位数如下：

N=100时

sn0=sum0(100);

sn=sum1(100)

n=fix(-log10(2*abs(sn-sn0)))

得到：sn=0.7400495n=7

N=10e4时

sn0=sum0(10e4);

sn=sum1(10e4)

n=fix(-log10(2*abs(sn-sn0)))

得到：sn=0.7498521n=3

N=10e6时

sn0=sum0(10e6);

sn=sum1(10e6)

n=fix(-log10(2*abs(sn-sn0)))

仅供学习与交流，如有侵权请联系网站删除谢谢1

精品好文档，推荐学习交流

得到：sn=0.7498521n=3

按第二种顺序得到的值及有效位数如下：

N=100时

sn0=sum0(100);

sn=sum2(100)

n=fix(-log10(2*abs(sn-sn0)))

得到：sn=0.7400495n=7

N=10e4时

sn0=sum0(10e4);

sn=sum2(10e4)

n=fix(-log10(2*abs(sn-sn0)))

得到：sn=0.7499900n=7

N=10e6时

sn0=sum0(10e6);

sn=sum2(10e6)

n=fix(-log10(2*abs(sn-sn0)))

得到：sn=0.7499999n=7

（4）通过这道上机题，我明白了，应用计算机进行求和运算时，求和的顺序不

同对结果的精度是有影响的。具体而言，小数先求和，大数后求和算得的结果更

精确；反之，则会出现“大数吃小数”的现象，从而精度会降低。

第二章（20题）

（1）用牛顿法解f（x）=0的根的通用程序如下：

functionx=newton(f,df,x0,ep)

%用牛顿法解f（x）=0的根

%f为待求解函数，df为f的导函数，x0为迭代初值，ep为允许误差

epi=100;

%为了保证能进入下列循环，epi取一个大于ep的初值，取epi=100

forepiep

x=x0-feval(f,x0)/feval(df,x0);

epi=abs(x-x0);

x0=x;

end

end

（2）

1）编程算法如下：

显然x0=0时，迭代序列收敛于根x2。设置步长b=0.00001，从x0=0开始每

次加一个步长b,然后判断迭代序列是否收敛于根x2,若收敛于根x2，则继续加

仅供学习与交流，如有侵权请联系网站删除谢谢2

精品好文档，推荐学习交流

步长b；若不收敛与根x2，则退出循环，此时的x0即为题中所求的尽可能大的

δ。

程序如下：

正向找：

f=@(x)x^3/3-x;

df=@(x)x^2-1;

x0=0;x=0;

whilex==0

x=newton(f,df,x0,1e-10);

x0=x0+0.00001;

end

x0

得到：x0=0.77461

负向找

f=@(x)x^3/3-x;

df=@(x)x^2-1;

x0=0;x=0;

whilex==0

x=newton(f,df,x0,1e-10);

x0=x