弦截法非线性方程求解.docxVIP

《MATLAB 程序设计实践》课程考核 一、,编程实现以下科学计算算法， 并举一例应用之。（参考书籍《精通 MATLAB 科学计算》，王正林等著，电子工业出版社， 2009 年） “弦截法非线性方程求解” 1、算法说明 : 弦截法的算法过程如下： （1） 过两点 (a,f(a)),(b,f(b)) 作一直线，它与 x 轴有一个交点，记为 x1.. （2） 如果 f(a)f(x1)0, 过两点 (a,f(a)),(x1,f(x1)) 作一直线，它与 x 轴的交点记为 x2，否则过两点 (b,f(b)),(x1,f(x1)) 作一直线，它与 x 轴的交点 记为 x2 ； （3） 如此下去，直到 |xn-xn-1| ，就可以认为 xn 为 f(x)=0 在区间 [a,b] 上的一个根。 xk xk 1 a 1 ) （4） Xk 的递推公式为： f ( xk xk 1 xk b 1 ) f ( xk  a 0 f (a), f (a) f ( xk 1 ) f ( a) 且 b 0 f (b), f (a) f ( xk 1 ) f (b) x1 b a a f (a) f (b) f (a) 在 MATLAB 中编程实现的弦截法的函数为： Secant. 功能：用弦截法求函数在某个区间的一个零点。 调用格式： root=Secant(f,a,b,eps). 其中， f