东南大学 数值分析上机题作业 MATLAB版.docx

2015 2015.1.9 上机作业题报告  上机作业题报告 1.Chapter 1 1.1题目 设SN=j=2N1 （1）编制按从大到小的顺序，计算SN的通用程序。 （2）编制按从小到大的顺序，计算SN的通用程序。 （3）按两种顺序分别计算,并指出有效位数。（编制程序时用单精度） （4）通过本次上机题，你明白了什么？ clear;N=input( clear; N=input(请输入N值:); Ac=single((3/2-1/N-1/(N+1))/2); Snl2s=single(0); Sns2l=single(0); for i=2:N Snl2s=Snl2s+1/(i*i-1); end for i=N:-1:2 Sns2l=Sns2l+1/(i*i-1); end fprintf(精确值为： %f\n,Ac); fprintf(从大到小的顺序累加得SN=%f\n,Snl2s); fprintf(从小到大的顺序累加得SN=%f\n,Sns2l); disp(========================================================); P20T17请输入N P20T17 请输入N值:10^2 精确值为： 0.740049 从大到小的顺序累加得SN=0.740049 从小到大的顺序累加得SN=0.740050 ============================================================ P20T17 请输入N值:10^4 精确值为： 0.749900 从大到小的顺序累加得SN=0.749852 从小到大的顺序累加得SN=0.749900============================================================ P20T17 从小到大的顺序累加得SN=0.749900 ============================================================ P20T17 请输入N值:10^6 精确值为： 0.749999 从大到小的顺序累加得SN=0.749852 从小到大的顺序累加得SN=0.749999 ============================================================ 按从大到小的顺序，有效位数分别为：6，4，3。 按从小到大的顺序，有效位数分别为：5，6，6。 可以看出，不同的算法造成的误差限是不同的，好的算法可以让结果更加精确。当采用从大到小的顺序累加的算法时，误差限随着N的增大而增大，可见在累加的过程中，误差在放大，造成结果的误差较大。因此，采取从小到大的顺序累加得到的结果更加精确。 2.Chapter 2 2.1题目 （1）给定初值及容许误差，编制牛顿法解方程f(x)=0的通用程序。 （2）给定方程,易知其有三个根 eq \o\ac(○,1)由牛顿方法的局部收敛性可知存在当时，Newton迭代序列收敛于根x2*。试确定尽可能大的δ。 eq \o\ac(○,2)试取若干初始值，观察当时Newton序列的收敛性以及收敛于哪一个根。 （3）通过本上机题，你明白了什么？ 2.2程序 f(x)函数m文件：fu.m function Fu=fu(x) Fu=x^3/3-x; end f(x)函数m function Fu=dfu(x) Fu=x^2-1; end 用Newton法求根的通用程序Newton.m clear; x0=input(请输入初值x0:); ep=input(请输入容许误差:); flag=1; while flag==1 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0)ep flag=0; end x0=x1; end fprintf(方程的一个近似解为：%f\n,x0); 寻找最大δ值的程序：Find.m clear eps=input(请输入搜索精度：); ep=input(请输入容许误差：); flag=1; k=0; x0=0; while flag==1 sigma=k*eps; x0=sigma; k=k+1; m=0; flag1=1; while flag1==1 m=10^3 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0)ep flag1=0; end m=m+1; x0=x1; end if flag1==1||a