数值分析实验四(金).docVIP

实验名称： 实验四 方程求根 指导教师： 数值分析实验组 实验时数： 2 实验设备：安装了Matlab、C++、VF软件的计算机 实验日期：2014年 月 日 实验地点： 第五教学楼北802或902 实验目的： 1. 掌握非线性方程数值解法的基本思想和基本步骤。 2. 理解各类数值方法的优缺点，并能自行编程求解。 3. 初步了解非线性方程的简单迭代法及其收敛性，体会迭代函数对收敛性的影响，体会不同初值对同一迭代函数的影响。 实验准备： 1. 在开始本实验之前，请回顾教科书的相关内容； 2. 需要一台准备安装Windows XP Professional操作系统和装有数学软件的计算机。 实验内容及要求 A题的根，准确到10。 （1）设计一种线性收敛的迭代法求方程的根，输出迭代初值、各次迭代值及迭代次数。 （2）用牛顿迭代法求方程的根，输出迭代初值、各次迭代值及迭代次数，并与（1）的结果比较。 （3）用MATLAB内部函数solve直接求出方程的所有根，并与（1）、（2）的结果进行比较。 在=1.5附近的根。 （1）对牛顿迭代公式： ，编写程序进行实验，分别取，迭代10次，观察比较其计算值，并分析原因。 （2）用MATLAB内部函数solve直接求出方程的所有根，并与（1）的结果进行比较。 的一个根， 当时颇为轰动，但无人知道他是用什么方法得到的。现在，请你试试用二分法和Newton迭代法求解上述方程能否得到这个结果。 D题 用简单迭代法求方程 的根。 方案一： 化为等价方程 方案二： 化为等价方程 （1）分别对方案一、方案二取初值，迭代10次，观察其计算值，并加以分析。 （2）用MATLAB内部函数solve直接求出方程的所有根，并与（1）的结果进行比较。 二分法： function[x,k]=demimethod(a,b,f,eps) %a,b表示求解区间[a,b]的端点 %f 表示所求方程函数名 %eps是精度指标 %x表示所求近似解 %k表示循环次数 fa=feval(f,a); fab=feval(f,(a+b)/2); k=0; while abs(b-a)eps if fab==0 x=(a+b)/2; return; elseif fa*fab0 b=(a+b)/2; else a=(a+b)/2; end fa=feval(f,a); fab=feval(f,(a+b)/2); k=k+1; end x=(a+b)/2; 在MATLAB命令窗口中输入下列命令： f=@(x)x^3+2*x^2+10*x-20; [x0,k]=demimethod(1,2,f,10^-9) 求得结果如下： x0 = 1.3688 k = 30 Newton迭代法： function?root=NewtonRoot(f,a,b,eps)format?longif(nargin==3)????eps=1.0e-5;endf1=subs(sym(f),findsym(sym(f)),a);f2=subs(sym(f),findsym(sym(f)),b);if(f1==0)????root=a;endif(f2==0)????root=b;endtol=1;fun=diff(sym(f)); %求导数fa=subs(sym(f),findsym(sym(f)),a);fb=subs(sym(f),findsym(sym(f)),b);dfa=subs(sym(fun),findsym(sym(fun)),a);dfb=subs(sym(fun),findsym(sym(fun)),b);if(dfadfb) %初始值取两端点导数较大者????root=a-fa/dfa;else????root=b-fb/dfb;endwhile(toleps)????r1=root;????fx=subs(sym(f),findsym(sym(f)),r1);????dfx=subs(sym(fun),findsym(sym(fun)),r1); %求该点的导数值root=r1-fx/dfx; %迭代的核心公式tol=abs(root-r1);end 在MATLAB命令窗口下输入下列命令： r=NewtonRoot(‘x^3+2*x^2+10*x-20’,1,2) 得到结果如下： r =