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

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

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

3、3;a=4;k=0;y0=f(b);y=f(a);while abs(b-a)>0.00001 t=a-y*(a-b)/(y-y0); b=a;y0=y;%求解区间;、%牛顿迭代法求解程序;四、实验数据结果及分析表1 :二分法程序结果迭代次数区间值：a区间值：b13.000000000000003.5000000000000023.250000000000003.5000000000000033.250000000000003.3750000000000043.312500000000003.37500000000000143.347351074218753.34741210937500

4、153.347381591796883.34741210937500163.347396850585943.34741210937500173.347396850585943.34740447998047183.347396850585943.34740447998047表2 :牛顿迭代法程序结果迭代次数区间值：b区间值：a13.438282138662913.3199556816049223.319955681604923.3483632970400433.348363297040043.3474127204823343.347412720482333.3474028396087953.34

5、7412720482333.34740283960879五、实验结论通过表1可知，在二分法下，程序迭代了 17次后和第18次的结果一致，即程序迭代了 17次达到要求的试验误差；通过表 2可知，在牛顿迭代法下，程序迭代了 4次后和第5次 的结果一致，即程序迭代了 4次达到要求的试验误差；二者比较明显可以看出牛顿迭代法的求解效率要远远优于二分法。a=t;y=f(a); k=k+1;dis p(b,a);k%显示各个区间值和求解次数;enddis p(b,a);%显示最后一个区间值;多面体旋转实验200822401018 徐小良、问题叙述：编写matlab程序实现对正立方体的旋转，并用适当的方法来验

6、证程序设计的正确性。二、问题分析:使用相对应的三个正交矩阵即可实现对三位图形进行各个方向的旋转，在此不再赘述。10Qx (a) = 0 cosa b sinot 使用moviein命令,cosY -sin?，Qz(Y)= sinVI 0cosp 0 sin P-sinot f Qy(P)=0co翘L_sinp 0就可以对原正立方体和旋转以后的三位图形进行全方位的观测；为观1 0cosp JcosY0001测的方便，可在立方体各个面上涂不同的颜色加以区分。三、试验程序以及注释：(1)、主题函数：clear;清除以前所有数据，以防干扰X,Y,Z=peaks(30);m=moviei n(30);B

7、0=0 0 0;1 0 0;1 1 0;0 1 0;0 0 1;1 0 1;1 1 1;0 1 1;n=8;%view(-25,16);B=B0;cube(B);for i=1:30扑捉相关画面延时相关因子endview(10*(i-1),10*(i-1) % m(:,i)=getframe;设置观测点观察立方体delay( n);延时一定时间，以便观测Qz=cos( pi/2) -si n(pi/2) O;si n(pi/2) cos( pi/2) 0;0 0 1;Qy=cos( -p i/4) 0 sin(-pi /4);0 1 0;-si n(-pi/4) 0 cos(- pi/4);Q

8、x=1 0 0;0 cos( pi/4) -si n(p i/4);0 si n(pi/4) cos( pi/4);B=B*Qz：cube(B);B=B*Qz：cube(B);B=B*Qz'cube(B);B=B0*Qy：B=B*Qx：B(:,3)=B(:,3)+1; cube(B);view(-25,16);for i=1:30view(3-10*(i-1),10) m(:,i)=getframe; delay (n);end(2)、Cube 函数:立方体旋转在不同面涂不同颜色延时函数toe延时时间计算fun cti on cube(B)fac=1 2 3 4;1 2 6 5;1 4

9、 8 5;7 8 5 6;7 3 2 6;7 3 4 8;patch('faces',fac(1,:),'vertices',B,'faceColor','m'); % patch('faces',fac(2,:),'vertices',B,'faceColor','b'); patch('faces',fac(3,:),'vertices',B,'faceColor',T); patch('faces'

10、,fac(4,:),'vertices',B,'faceColor','c'); p atch('faces',fac(5,:),'vertices',B,'faceColor','g'); patch('faces',fac(6,:),'vertices',B,'faceColor','y');(3)、delay 函数:fun cti on delay (n)ticM=10000000;N=M* n;for k=0:NM=M*1-1;end四、试验结果：通过程序仿真可观测到，立方体不同侧面的颜色，以此检测图色的正确性。 El 口oe0.60 40.2当旋转后，经过程序仿真观测旋转后的图形各个侧面的颜色，以此观测旋转的正确性。£.11* Edit ICif* Insert&#