分法和牛顿迭代法求解方程的比较

《分法和牛顿迭代法求解方程的比较》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、二分法和牛顿迭代法求解方程的比较200822401018徐小良一、问题叙述求解的解；通过编写matlab程序分别用分析二分法和牛顿迭代法求解方程，通过两种方法的比较，分析二者求解方程的快慢程度。二、问题分析由matlab画图命令，容易得到此方程解的范围为（2，4）；两种迭代方法，在使用相同的误差（0.00001）的情况下，得出matlab迭代次数，通过次数的比较得出二者求解速度快慢比较。三、实验程序及注释（1）、二分法程序：clear;%清除所有内存数据；f=inline('12-3*x+2*cos(x)');f

2、ormatlong%数据显示格式设为长型；a=2;b=4;%求解区间；er=b-a;ya=f(a);k=0;er0=0.00001;%误差分析；whileer>er0x0=.5*(a+b);y0=f(x0);ifya*y0<0b=x0;%二分法求解程序；elsea=x0;ya=y0;enddisp([a,b]);er=b-a;k=k+1%显示各个区间值和求解次数；enddisp([a,b]);%显示最后一个区间值；（2）、牛顿迭代法程序：clear;%清除所有内存数据；f=inline('12-3*x+2*cos

3、(x)');formatlong%数据显示格式设为长型；b=3;a=4;k=0;%求解区间；y0=f(b);y=f(a);whileabs(b-a)>0.00001t=a-y*(a-b)/(y-y0);b=a;y0=y;%牛顿迭代法求解程序；a=t;y=f(a);k=k+1;disp([b,a]);k%显示各个区间值和求解次数；enddisp([b,a]);%显示最后一个区间值；一、实验数据结果及分析表1：二分法程序结果迭代次数区间值：a区间值：b13.000000000000003.50000000000000

4、23.250000000000003.5000000000000033.250000000000003.3750000000000043.312500000000003.37500000000000......…………143.347351074218753.34741210937500153.347381591796883.34741210937500163.347396850585943.34741210937500173.347396850585943.34740447998047183.3473968505

5、85943.34740447998047表2：牛顿迭代法程序结果迭代次数区间值：b区间值：a13.438282138662913.3199556816049223.319955681604923.3483632970400433.348363297040043.3474127204823343.347412720482333.3474028396087953.347412720482333.34740283960879二、实验结论通过表1可知，在二分法下，程序迭代了17次后和第18次的结果一致，即程序迭代了17次

6、达到要求的试验误差；通过表2可知，在牛顿迭代法下，程序迭代了4次后和第5次的结果一致，即程序迭代了4次达到要求的试验误差；二者比较明显可以看出牛顿迭代法的求解效率要远远优于二分法。多面体旋转实验200822401018徐小良一、问题叙述：编写matlab程序实现对正立方体的旋转，并用适当的方法来验证程序设计的正确性。二、问题分析：使用相对应的三个正交矩阵即可实现对三位图形进行各个方向的旋转，在此不再赘述。，，使用moviein命令，就可以对原正立方体和旋转以后的三位图形进行全方位的观测；为观测的方便，可在立方体各

7、个面上涂不同的颜色加以区分。三、试验程序以及注释：（1）、主题函数：clear;%清除以前所有数据，以防干扰[X,Y,Z]=peaks(30);m=moviein(30);%扑捉相关画面B0=[000;100;110;010;001;101;111;011];n=8;%延时相关因子view(-25,16);B=B0;cube(B);fori=1:30view(10*(i-1),10*(i-1))%设置观测点m(:,i)=getframe;%观察立方体delay(n);%延时一定时间，以便观测endQz=[cos(

8、pi/2)-sin(pi/2)0;sin(pi/2)cos(pi/2)0;001];Qy=[cos(-pi/4)0sin(-pi/4);010;-sin(-pi/4)0cos(-pi/4)];Qx=[100;0cos(pi/4)-sin(pi/4);0sin(pi/4)cos(pi/4)];B=B*Qz';cube(B);B=B*Qz';cube(B);%立方体旋转B=B