东南大学数值分析上机题.docxVIP

数值分析上机题刘乃东 161743舍入误差与有效数所使用软件为Dev-C++。#include iostreamusing namespace std;int main(void){ float N,S=0; cinN; for(;N1;N--) {S=S+1/(N*N-1);}coutSendl;system(pause);return 0;}#include iostreamusing namespace std;int main(void){ float N,S=0; cinN; for(double i=2;iN;i++) {S=S+1/(i*i-1);}coutSendl;system(pause);return 0;}（1） （2）（3）算法（1）0.7399490.7498520.749852算法（2）0.7400500.7499000.749999精确值0.7400500.7499000.749999算法（1）有效位数344算法（2）有效位数666（4）若干同号数相加，采用绝对值较小者先加的算法，其结果的相对误差较小，但对于双精度数则不必考虑。Newton迭代法本题及之后题目所使用的软件为Matlab R2015a。function result=Newton(a,b,c)% 使用牛顿迭代，求方程a=0的根，b为迭代初值，c为精度要求l=-log(c)/log(10);syms x;fx=a;ffx=diff(fx,x);Fx=x-fx/ffx;x1=b;x2=subs(Fx,x,x1);i=0;while(abs(x2-x1)ci100) x1=vpa(x2,l); x2=vpa(subs(Fx,x,x1),l); i=i+1;endif i=100 disp(?T?ù); result=220161743;else result=eval(vpa(x2,l)); if abs(result)c result=0; endend% result(1,2)=i;end（1） （2）①function result=Newtons(a)i=0;j=0.01;while(j0.0000001)while(i1) s1=Newton(a,i,j/2); s2=Newton(a,-i,j/2); if s1~=220161743s2~=220161743s1==0s2==0 i=i+j; else break; endendi=i-j;j=j/10;endresult=i-j;end计算结果为0unction result=delt(a)syms x;j=1;for i=-3:0.11:3 result(j,1)=i; result(j,2)=Newton(a);endend②在[-3,3]上，以0.11为步长，判断是否收敛及收敛于哪一个根，计算结果见下表。初始值根初始值根初始值根初始值根-3-1.73-1.46-1.730.0801.621.73-2.89-1.73-1.35-1.730.1901.731.73-2.78-1.73-1.24-1.730.301.841.73-2.67-1.73-1.13-1.730.4101.951.73-2.56-1.73-1.02-1.730.5202.061.73-2.45-1.73-0.911.730.6302.171.73-2.34-1.73-0.8-1.730.7402.281.73-2.23-1.73-0.6900.85-1.732.391.73-2.12-1.73-0.5800.96-1.732.51.73-2.01-1.73-0.4701.071.732.611.73-1.9-1.73-0.3601.181.732.721.73-1.79-1.73-0.2501.291.732.831.73-1.68-1.73-0.1401.41.732.941.73-1.57-1.73-0.0301.511.73由表中可以看出，在（-3，-1）上收敛于负根，在（-0.775,0.775）上收敛于零根，在（1,3）上收敛于正根，在（-1，-0.774）和（0.775,1）上不一致收敛。（3）说明在Newton迭代法中，初始值的选择很重要，选择不当则不能得到想要的结果。列主元Gauss消去法function result=Gauss(A,b)% 使用Gauss消去法求线性方程组的解，Ax=b。if size(A,1)==size(A,2) n=size(A,1); result=zeros(n,1);else result=矩阵输入错误; return;endA=[A,b];for i=1:n-1 mel=find(abs(A((i:n),i))==max(abs