第六章数值分析.docxVIP

第六章方程求根3.试编写MATLAB程序，分别用二分法和试位法求下述方程的根，并给出各自达到精度要求所需计算数值的次数，这里设;（1）解得（2）又因为，所以确定为单根区间。f=@(x)x^3-cos(x);[c,err,k]=bisec(f,0,0.5*pi,0.5e-2)yb =3.8758c = 0.8682err = 0.0031k =8可得用在根区间二分法使用8次得到符合精度要求的根0.8682（3）可确定方程有单根，且在[1,3]之间。f=@(x)x-sin(x)-1;[c,err,k]=bisec(f,1,3,0.5e-2)yb =1.8589c = 1.9336err = 0.0039k = 8所以在[1,3]使用8次二分法得到符合精度要求的根1.93364.用迭代方法求解下述方程在给定附近的根，要求误差不超过；（3）解：迭代对应的迭代公式为，由，得，于是迭代在3附近收敛。由此迭代公式可得迭代结果见表1-1.k012345632.68282.64082.63482.63392.63382.6338所以为符合精度要求的附近的根。10分别用（1）Newton法，取；（2）割线法，取，求方程解：这里Newton迭代公式为，取，计算结果如下k01234-2.0-2.2222-2.1962-2.1958-2.1958割线法的迭代公式，取，计算结果如下k012345-2.0-2.1-2.2081-2.1951-2.1958-2.1958比较Newton法和割线法可以看出，Newton法收敛的率快一些，两者都具有较快的收敛速度。附Matlab二分法函数bisec函数：function [x,X,err,k,y]=newton(f,df,x0,epsilon,maxit)%Newtons method for finding the root of the equation f(x)=0;if nargin 3 error(more input are needed)endif nargin 4 epsilon =1e-12;endif nargin 5 maxit =100;endX=x0;flag=0;for k=1:maxit x=x0-f(x0)/df(x0);X=[X,x];err=abs(x-x0); x0=x;y=f(x0);if(err epsilon)||(abs(y)epsilon) flag=1; break;endendif flag==1 disp([The root x is num2str(x)]); disp([The error of x is num2str(err)]); disp([The number of iteration is num2str(k)]); disp([The function value at x is num2str(y)]); else disp([Newton method fails to find the root within... num2str(maxit) iterations]);endreturn