matlab实验报告方程求根实验,数值分析实验报告——方程求根

《数值分析实验报告——方程求根》由会员分享，可在线阅读，更多相关《数值分析实验报告——方程求根(6页珍藏版)》请在人人文库网上搜索。

1、数值分析实验报告班级 信科1501 学号 姓名 梁恩昊 日期 2017.10.3 学 院数学科学学院专 业信息与计算科学课程名称数值分析成 绩实验一 方程求根 一、 实验目的: 掌握二分法、Newton法、不动点迭代法、弦截法求方程的根的各种计算方法、并实施程序调试和运行，学习应用这些算法于实际问题。二、 实验内容: 二分法、Newton法、不动点迭代法、弦截法求方程的根、程序的调试和运行，给出实例的计算结果。观察初值对收敛性的影响。三、 实验步骤： 、二分法：定义：对于区间a，b上连续不断且f(a)f(b)0 a=p; else b=p; end;end;if i=n0&(fp=0|(abs。

2、(b-a)/2)=1e-6)&(n BDD(1)n =21ans =1.3652 、弦截法：定义：弦截法是求非线性方程近似根的一种线性近似方法。它是以与曲线弧AB对应的弦AB与x轴的交点横坐标作为曲线弧AB与x轴的交点横坐标的近似值来求出方程的近似解。该方法一般通过计算机编程来实现。弦截法的原理是以直代曲即用弦(直线)代替曲线求方程的近似解，也就是利用对应的弦 与 轴的交点横坐标来作为曲线弧 与 轴的交点横坐标 的近似值。实现方法：我们与二分法一样，先设一方程400*(x4)-300*(x3)+200*(x2)-10*x-1=0，并求其在区间0.1,1上的根。PS：本方法应用的软件为matla。

3、b。函数文件：functionp1,err,k,y=secant(f,p0,p1,delta,max1)% f是非线性函数 %p0,p1是初始值 %delta是给定允许误差 %max1是迭代次数的上限 %p1是所求得的方程的近似解 %err是p1-p0的误差估计 %k是所需要的迭代次数 %y=f(p1) K=0,p0,p1,feval(f,p0),feval(f,p1) for k=1:max1p2=p1-feval(f,p1)*(p1-p0)/(feval(f,p1)-feval(f,p0);err=abs(p2-p1); p0=p1; p1=p2; k,p1,err,y=feval(f,p。

4、1) if(errdelta)|(y=0), break,endendendM文件：secant(400*(x4)-300*(x3)+200*(x2)-10*x-1=0,0,1,1.0e-4,100);disp(用弦截法求得的方程的根为：x=)disp(p1)运行结果：用弦截法求得的方程的根为：x=0.1108四、 实验中遇到的问题及解决方法：MATLAB运用不熟练，于是翻开以前的课本与ppt，并结合网上查阅的资料，才又掌握了一些基础的运用方法。一开始忘了怎么创建函数文件，于是实现二分法的程序时变得十分麻烦。后来复习过后，运用在牛顿法上，就简单了许多。二分法与牛顿法是比较容易掌握并实现的，而后两个方法则是查阅了很多资料，却仍然似懂非懂，程序也是大部分参考了资料。我将在此之后继续钻研它们，直到能熟练运用为止。