数值分析上机题课后作业全部-东南大学.docxVIP

1.Chapter11.1题目设，其精确值为。（1）编制按从大到小的顺序，计算SN的通用程序。（2）编制按从小到大的顺序，计算SN的通用程序。（3）按两种顺序分别计算,并指出有效位数。（编制程序时用单精度）（4）通过本次上机题，你明白了什么？1.2程序1.3运行结果1.4结果分析按从大到小的顺序，有效位数分别为：6，4，3。按从小到大的顺序，有效位数分别为：5，6，6。可以看出，不同的算法造成的误差限是不同的，好的算法可以让结果更加精确。当采用从大到小的顺序累加的算法时，误差限随着N的增大而增大，可见在累加的过程中，误差在放大，造成结果的误差较大。因此，采取从小到大的顺序累加得到的结果更加精确。2.Chapter22.1题目（1）给定初值及容许误差，编制牛顿法解方程f(x)=0的通用程序。（2）给定方程,易知其有三个根由牛顿方法的局部收敛性可知存在当时，Newton迭代序列收敛于根x2*。试确定尽可能大的。试取若干初始值，观察当时Newton序列的收敛性以及收敛于哪一个根。（3）通过本上机题，你明白了什么？2.2程序函数m文件：fu.mfunction Fu=fu(x)Fu=x^3/3-x;end函数m文件：dfu.mfunction Fu=dfu(x)Fu=x^2-1;end用Newton法求根的通用程序Newton.mclear;x0=input(请输入初值x0:);ep=input(请输入容许误差:);flag=1;while flag==1 x1=x0-fu(x0)/dfu(x0);if abs(x1-x0)epflag=0;endx0=x1;endfprintf(方程的一个近似解为：%f\n,x0);寻找最大δ值的程序：Find.mcleareps=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;endif flag1==1||abs(x0)=ep flag=0; endendfprintf(最大的sigma值为：%f\n,sigma);2.3运行结果（1）寻找最大的值。算法为：将初值x0在从0开始不断累加搜索精度eps，带入Newton迭代公式，直到求得的根不再收敛于0为止，此时的x0值即为最大的sigma值。运行Find.m，得到在不同的搜索精度下的最大sigma值。 Find请输入搜索精度：10^-6请输入容许误差：10^-6最大的sigma值为：0.774597 Find请输入搜索精度：10^-4请输入容许误差：10^-6最大的sigma值为：0.774600 Find请输入搜索精度：10^-2请输入容许误差：10^-6最大的sigma值为：0.780000（2）运行Newton.m在内取初值，运行结果如下：X0Xk-1000-1.732051-500-1.732051-100-1.732051-10-1.732051-5-1.732051-2.5-1.732051-1.5-1.732051可见，在区间内取初值，Newton序列收敛，且收敛于根。在内取初值，运行结果如下：X0Xk-0.951.732051-0.851.732051-0.81.732051-0.7745981.732051可见，在内取初值，Newton序列收敛，且收敛于根。在内内取初值，运行结果如下：X0Xk-0.7745960.000000-0.550.000000-0.350.000000-0.150.0000000.050.0000000.250.0000000.450.0000000.650.0000000.7745960.000000可见，在内取初值，Newton序列收敛，且收敛于根0。在内取初值，运行结果如下：X0Xk0.774598-1.7320510.8-1.7320510.85-1.7320510.95-1.732051可见，在内取初值，Newton序列收敛，且收敛于根在内取初值，运行结果如下：X0Xk1.51.7320512.51732051101.7320511001.7320515001.73205110001.732051可见，在内取初值，Newton序列收敛，且收敛于根3.Chapter33.1题目对于某电路的分析，归结为求解线性方程组RI=V，其