不动点非线性方程求解.docVIP

《MATLAB程序设技实践》课程考核 1编程实现以下科学计算算法，并举一例应用之，（参考书籍《精通MATLAB科学计算》），王正林等著，电子工业出版社，2009年）。“不动点非线性方程求解” ①算法说明： 在MATLAB中编程实现不动点迭代法的函数是：StablePoint。 功能：用不动点迭代法求函数的一个零点。 调用格式：[root，n]=StablePoint（f,x0,eps）。 其中，f为函数名； x0为初始迭代向量； eps为根的精度； root为求出的函数零点； n为迭代步数。 ②流程图： N Y ③源程序代码： Function [root,n]=StablePoint(f,x0,eps) %用不动点迭代法求函数的一个零点 %初始迭代量：x0 %根的精度：eps %求出的函数零点：root %迭代步数：n if(nargin==2) %输入参数个数为2 eps=1.0e-4; end tol=1; %绝对误差限 root=x0; n=0; while(toleps) n=n+1; r1=root; Root=subs(sym(f),findsym(sym(f)),r1)+r1; %迭代的核心公式 tol=abs(root-r1); end ④举例：采