第13章 实验十一非线性方程求根.doc

第13章 实验十一非线性方程求根 实验目的：学会用两分法、牛顿迭代法求解非线性方程的根 两分法： 在数学运算中，经常遇到解方程的问题。例如要求以下方程的根： （13.1） 设是方程的根，也称为函数的零点，即。 二分法定理：设在区间[a,b]内连续，与的符号相反，表示两分法生成的中点序列，则有： 其中 这样的序列收敛到零点，即。 例13.1. 用两分法求方程的近似根，保留小数点后四位有效数字。 解：（1）先通过Matlab画出函数的图形。 画图程序为： 所得图形见图13.1。 图13.1 函数的曲线图 可知方程的根位于1.2～1.4之间。 再用两分法求出近似根。 （2）写出两分法的源程序代码 .m文件： function bisec_n(f_name,a,c) fprintf(Bsenction Scheme\n\n); tolerance=0.00001;it_limit=20; fprintf(It. a b c fa=f(a) fc=f(c) abs(a-c)\n); it=0; y_a=feval(f_name,a); y_c=feval(f_name,c); if(y_a*y_c0); fprintf(\n\n Stopped becaude f(a)f(c)0 \n); else while 1 it=it+1; b=(a+c)/2;y_b=feval(f_name,b); fprintf(%3.0f,%10.6f,%10.6f,it,a,b); fprintf(%10.6f, %10.6f,%10.6f,c,y_a,y_c); fprintf(%12.3e\n,abs(y_c-y_a)); if(abs(c-a)/2= tolerance ); fprintf(Tolerance is satisfied. \n);break fprintf(\n Change a or b and run again.\n); end if(itit_limit); fprintf(Iteration limit exceeded.\n);break end if(y_a*y_b 0) c=b; y_c=y_b; else a=b; y_a=y_b; end end fprintf(Finial result:Root=%12.6f \n,b); end （3）写出函数f= sqrt(9+x.^2)-tan(x)的.m文件： function f=eqn_w1(x) f= sqrt(9+x.^2)-tan(x); （4）写出执行函数命令： bisec_n(eqn_w1,1.2,1.4); 结果显示如下 Bsenction Scheme It. a b c fa=f(a) fc=f(c) abs(a-c) 1, 1.200000, 1.300000 1.400000, 0.658947, -2.487295 2.000e-001 2, 1.200000, 1.250000 1.300000, 0.658947, -0.332546 1.000e-001 3, 1.250000, 1.275000 1.300000, 0.240430, -0.332546 5.000e-002 4, 1.250000, 1.262500 1.275000, 0.240430, -0.021829 2.500e-002 5, 1.262500, 1.268750 1.275000, 0.114618, -0.021829 1.250e-002 6, 1.268750, 1.271875 1.275000, 0.047807, -0.021829 6.250e-003 7, 1.271875, 1.273438 1.275000, 0.013353, -0.021829 3.125e-003 8, 1.271875, 1.272656 1.273438, 0.013353, -0.004145 1.563e-003 9, 1.272656, 1.273047 1.273438, 0.004627, -0.004145 7.813e-004 10, 1.273047, 1.273242