数值分析上机题目详解.docxVIP

第一章一、题目设，其精确值为。编制按从大到小的顺序，计算SN的通用程序。编制按从小到大的顺序，计算SN的通用程序。按两种顺序分别计算,并指出有效位数。（编制程序时用单精度）通过本次上机题，你明白了什么？二、通用程序N=input(Please Input an N (N1):);AccurateValue=single((0-1/(N+1)-1/N+3/2)/2);Sn1=single(0);for a=2:N;Sn1=Sn1+1/(a^2-1);endSn2=single(0);for a=2:N; Sn2=Sn2+1/((N-a+2)^2-1);endfprintf(The value of Sn (N=%d)\n,N);fprintf(Accurate Calculation %f\n,AccurateValue);fprintf(Caculate from large to small %f\n,Sn1);fprintf(Caculate from small to large %f\n,Sn2);disp(____________________________________________________)三、结果从结果可以看出有效位数是6位。感想：可以得出，算法对误差的传播有一定的影响，在计算时选一种好的算法可以使结果更为精确。从以上的结果可以看到从大到小的顺序导致大数吃小数的现象，容易产生较大的误差，求和运算从小数到大数所得到的结果才比较准确。第二章一、题目（1）给定初值及容许误差，编制Newton法解方程f(x)=0的通用程序。（2）给定方程,易知其有三个根由Newton方法的局部收敛性可知存在当时，Newton迭代序列收敛于根x2*。试确定尽可能大的。b)试取若干初始值，观察当时Newton序列的收敛性以及收敛于哪一个根。（3）通过本上机题，你明白了什么？二、通用程序1、定义函数和导函数%% 定义函数f(x)function Fx=fx(x) Fx=x^3/3-x;----------------------------------%% 定义导函数df(x)function Fx=dfx(x) Fx=x^2-1;2、寻找最大的clearflag=1;k=1;x0=0;while flag==1 delta=k*10^-6; x0=delta; k=k+1; m=0; flag1=1; while flag1==1 m=10^3 x1=x0-fx(x0)/dfx(x0); if abs(x1-x0)10^-6 flag1=0; end m=m+1; x0=x1; end if flag1==1||abs(x0)=10^-6 flag=0; endendfprintf(The maximun delta is %f\n,delta);3、Newton法求方程的根clearef=10^-6; %% 给定容许误差10^-6k=0;x0=input(Please input initial value Xo:);disp(k Xk);fprintf(0 %f\n,x0); flag=1;while flag==1 k=10^3 x1=x0-fx(x0)/dfx(x0); if abs(x1-x0)ef flag=0; end k=k+1; x0=x1;fprintf(%d %f\n,k,x0); end三、结果寻找最大的δ值：在题目给出的不同区间内进行牛顿迭代： 综上所述：首先得到最大的δ值为0.774597。在各区间任意取初值进行迭代得到(-∞,-1)区间收敛于-1.73205，(-1,-δ)区间局部收敛于0， (-δ,δ)区间收敛于0,(δ,1)区间类似于(-1,δ)区间，收敛于0，(1,∞)收敛于1.73205。感想：通过本上机题，了解到用Newton法求多根方程的根时，迭代序列收敛于某一个根有一定的区间限制。如果不清楚这个限制随意取值的话，会出现在一个区间上局部收敛于不同的根的情况。如下面的迭代：第三章一、题目列主元Gauss消去法对于某电路的分析，归结为求解线性方程组。其中编制解n阶线性方程组的列主元高斯消去法的通用程序；用所编程序线性方程组，并打印出解向量，保留5位有效数；二、通用程序%% 列主元Gauss消去法求解线性方程组%%%%参数输入n=input(Please input the order of matrix A: n=); %输入线性方程组阶数nb=zeros(1,n);A=input(Input matrix A (such as