实验一非线性方程的数值解法（一）信息与计算科学金融崔振威201002034031一、实验目的：熟悉二分法和简单迭代法的算法实现。
二、实验内容：教材P40 2.1.5三、实验要求1 根据实验内容编写二分法和简单迭代法的算法实现2 简单比较分析两种算法的误差3 试构造不同的迭代格式，分析比较其收敛性(一)、二分法程序：function ef=bisect(fx,xa,xb,n,delta)% fx是由方程转化的关于x的函数，有fx=0。
% xa 解区间上限% xb 解区间下限% n 最多循环步数，防止死循环。
%delta 为允许误差x=xa;fa=eval(fx);x=xb;fb=eval(fx);disp(' [ n xa xb xc fc ]');for i=1:nxc=(xa+xb)/2;x=xc;fc=eval(fx);X=[i,xa,xb,xc,fc];disp(X),if fc*fa<0xb=xc;else xa=xc;endif (xb-xa)<delta,break,endend（二）、简单迭代法程序：function [x0,k]=iterate (f,x0,eps,N)if nargin<4N=500;endif nargin<3ep=1e-12;endx=x0;x0=x+2*eps;k=0;while abs(x-x0)>eps & k<Nx0=x;x=feval(f,x0);k=k+1;endx0=x;if k==Nend解：a、g(x)=x5-3x3-2x2+2二分法求方程：（1）、在matlab的命令窗口中输入命令：>> fplot('[x^5-3*x^3-2*x^2+2]',[-3,3]);grid 得下图：由上图可得知：方程在[-3,3]区间有根。
（2）、二分法输出结果>> f='x^5-3*x^3-2*x^2+2'f =x^5-3*x^3-2*x^2+2>> bisect(f,-3,3,20,10^(-12))2.0000 -3.0000 0 -1.5000 0.03133.0000 -3.0000 -1.5000 -2.2500 -31.61824.0000 -2.2500 -1.5000 -1.8750 -8.43015.0000 -1.8750 -1.5000 -1.6875 -2.96326.0000 -1.6875 -1.5000 -1.5938 -1.21817.0000 -1.5938 -1.5000 -1.5469 -0.53828.0000 -1.5469 -1.5000 -1.5234 -0.24059.0000 -1.5234 -1.5000 -1.5117 -0.101510.0000 -1.5117 -1.5000 -1.5059 -0.034311.0000 -1.5059 -1.5000 -1.5029 -0.001412.0000 -1.5029 -1.5000 -1.5015 0.015013.0000 -1.5029 -1.5015 -1.5022 0.006814.0000 -1.5029 -1.5022 -1.5026 0.002715.0000 -1.5029 -1.5026 -1.5027 0.000716.0000 -1.5029 -1.5027 -1.5028 -0.000317.0000 -1.5028 -1.5027 -1.5028 0.000218.0000 -1.5028 -1.5028 -1.5028 -0.000119.0000 -1.5028 -1.5028 -1.5028 0.000120.0000 -1.5028 -1.5028 -1.5028 -0.00002、迭代法求方程：迭代法输出结果：>> f=inline('x^5-3*x^3-2*x^2+2');>> [x0,k]=iterate(fun1,2)x0 =2k =1>> [x0,k]=iterate(fun1,1.5)x0 =NaNk =6>> [x0,k]=iterate(fun1,2.5)x0 =NaNk =5（3）、误差分析：由二分法和迭代法输出结果可知，通过定点迭代法得出方程的解误差比二分法大，而利用二分法求出的结果中，可以清楚看出方程等于零时的解，其误差比迭代法小。
b、g(x)=cos(sin(x))二分法求方程：（1）、在matlab的命令窗口中输入命令：>> fplot('[cos(sin(x))]',[-4,4]);grid得下图：由上图可得知：方程在[-4,4]区间无根。
（2）、二分法输出结果>>f='cos(sin(x))'f =cos(sin(x))>> bisect(f,-4,4,20,10^(-12))2.0000 0 4.0000 2.0000 0.61433.0000 2.00004.0000 3.0000 0.99014.0000 3.0000 4.0000 3.5000 0.93915.0000 3.5000 4.0000 3.7500 0.84116.0000 3.7500 4.0000 3.8750 0.78427.0000 3.8750 4.0000 3.9375 0.75548.0000 3.9375 4.0000 3.9688 0.74129.0000 3.9688 4.0000 3.9844 0.734110.0000 3.9844 4.0000 3.9922 0.730511.0000 3.9922 4.0000 3.9961 0.728812.0000 3.9961 4.0000 3.9980 0.727913.0000 3.9980 4.0000 3.9990 0.727514.0000 3.9990 4.0000 3.9995 0.727315.0000 3.9995 4.0000 3.9998 0.727116.0000 3.9998 4.0000 3.9999 0.727117.0000 3.9999 4.0000 3.9999 0.727118.0000 3.9999 4.0000 4.0000 0.727019.0000 4.0000 4.0000 4.0000 0.727020.0000 4.0000 4.0000 4.0000 0.7270 2、迭代法求方程：迭代法输出结果：>> f=inline('cos(sin(x))');>> [x0,k]=iterate(f,0.5)x0 =0.7682k =15>> [x0,k]=iterate(f,1)x0 =0.7682k =15>> [x0,k]=iterate(f,1.5)x0 =0.7682k =16>> [x0,k]=iterate(f,2)x0 =0.7682k =15>> [x0,k]=iterate(f,2.5)x0 =0.7682k =14（3）、由于该方程无解，所以无法比较误差。
c、g(x)=x2-sin(x+0.15)二分法求方程：（1）、在matlab的命令窗口中输入命令：>> fplot('[x^2-sin(x+0.15)]',[-10,10]);grid 得下图：由上图可得知：方程在[-3,3]区间有根。
（2）、二分法输出结果>> f='x^2-sin(x+0.15)'f =x^2-sin(x+0.15)>> bisect(f,-3,3,30,10^(-12))1.0000 -3.0000 3.0000 0 -0.14942.0000 -3.0000 0 -1.5000 3.22573.0000 -1.5000 0 -0.7500 1.12714.0000 -0.7500 0 -0.3750 0.36375.0000 -0.3750 0 -0.1875 0.07266.0000 -0.1875 0 -0.0938 -0.04747.0000 -0.1875 -0.0938 -0.1406 0.01048.0000 -0.1406 -0.0938 -0.1172 -0.01919.0000 -0.1406 -0.1172 -0.1289 -0.004510.0000 -0.1406 -0.1289 -0.1348 0.002911.0000 -0.1348 -0.1289 -0.1318 -0.000812.0000 -0.1348 -0.1318 -0.1333 0.001113.0000 -0.1333 -0.1318 -0.1326 0.000114.0000 -0.1326 -0.1318 -0.1322 -0.000315.0000 -0.1326 -0.1322 -0.1324 -0.000116.0000 -0.1326 -0.1324 -0.1325 0.000017.0000 -0.1325 -0.1324 -0.1324 -0.000018.0000 -0.1325 -0.1324 -0.1325 -0.000019.0000 -0.1325 -0.1325 -0.1325 0.000020.0000 -0.1325 -0.1325 -0.1325 0.000021.0000 -0.1325 -0.1325 -0.1325 0.000022.0000 -0.1325 -0.1325 -0.1325 0.000023.0000 -0.1325 -0.1325 -0.1325 -0.000024.0000 -0.1325 -0.1325 -0.1325 0.000025.0000 -0.1325 -0.1325 -0.1325 -0.000026.0000 -0.1325 -0.1325 -0.1325 0.000027.0000 -0.1325 -0.1325 -0.1325 0.000028.0000 -0.1325 -0.1325 -0.1325 0.000029.0000 -0.1325 -0.1325 -0.1325 0.000030.0000 -0.1325 -0.1325 -0.1325 -0.00002、迭代法求方程：迭代法输出结果：>> f=inline('x^2-sin(x+0.15)');>> [x0,k]=iterate(f,1.96)x0 =NaNk =12>> [x0,k]=iterate(f,0,2)x0 =-0.1494k =1>> [x0,k]=iterate(f,0.2)x0 =0.3234k =500>> [x0,k]=iterate(f,0.3)x0 =0.3234k =500>> [x0,k]=iterate(f,0.001)x0 =0.3234k =500（3）、误差分析：由二分法和迭代法输出结果可知，利用二分法求出的结果中，可以清楚看出方程等于零时的解，其误差比迭代法小。