实验2解非线性方程实验.docVIP

实验 2 解非线性方程 姓名 报告日期 3月29日. 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：进一步熟练掌握解线性方程，计算 0的根 实验原理 二分法算法 牛顿迭代法 弦截法算法 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 实验总结 实验报告：根据实验情况和结果撰写并递交实验报告。 参考程序 1.二分法算法解方程 ①先建立二分法算法程序文件bisect.m，内容如下： function [c,err,yc] bisect f,a,b,delta ya feval f,a ; yb feval f,b ; if ya*yb 0,return,end max1 1+round log b-a -log delta /log 2 ; for k 1:max1 c a+b /2; yc feval f,c ; if yc 0 a c; b c; elseif yb*yc 0 b c; yb yc; else a c; ya yc; end if b-a delta, break,end end c a+b /2; err abs b-a ; yc feval f,c ; ②建立的函数文件f.m，内容如下： function f f x ; f x.*x-5; ③在command window按如下操作，得结果： [c,err,yc] bisect f,2,3,0.00001 c 2.2361 err 7.6294e-006 yc -1.3714e-005 [c,err,yc] bisect f,-2,-3,0.00001 c -2.2500 err 0.5000 yc 0.0625 即求得两根近似值。 2. 弦截法算法 ①先建立弦截法算法程序文件，内容如下： function [p1,err,k,y] secant f,p0,p1,delta,epsilon,max1 for k 1:max1 p2 p1-feval f,p1 * p1-p0 / feval f,p1 -feval f,p0 ; err abs p2-p1 ; relerr 2*err/ abs p2 +delta ; p0 p1; p1 p2; y feval f,p1 ; if err delta | relerr delta | abs y epsilon ,break,end end ②建立的函数文件f.m，内容如下： function f f x ; f x.*x-5; ③在command window按如下操作，得结果： [p1,err,k,y] secant f,2,3,0.00001,0.00001,10 p1 2.2361 err 4.3185e-005 k 4 y -2.2882e-007