求解非线性方程实验报告
非线性方程求解实验
第二部分：提高实验（例题讲解） 第二部分：提高实验（例题讲解） 1、编写二分法法的Matlab程序。 2、利用所编写的程序求解非线性方程的近似解。 例如：求方程 f ( x ) = x 3 − x − 1 = 0在区间[1,1.5]内 的根，要求误差不超过0.005. 首先定义函数erfenfa.m 然后在Matlab窗口下输入： >> erfenfa 输入函数f(x)='x^3-x-1' 输入区间=[1,1.5] 请输入误差=0.005
第二部分：提高实验（学生操作） 第二部分：提高实验（学生操作） 1、编写二分法法的Matlab程序。 2、利用所编写的程序求解非线性方程的近似解。 利用二分法求方程 x 3 − 2 x 2 − 4 x − 7 = 0在区间 [3,4]内的根，要求误差不超过0.0005.
第三部分：创新实验（学生操作） 第三部分：创新实验（学生操作） 1、编写Newton迭代法的Matlab程序。 2、利用源程序求解非线性方程的近似解。 利用牛顿法求方程 x 2 − 2 x − 3 = 0 的近似根，取初 始值为4，要求误差不超过0.0001.
实验二： 实验二：非线性方程典型解法的数值实验 数学实验室
本次实验的目的： 本次实验的目的： 1、掌握基本的绘图命令Plot,ezplot。 掌握基本的绘图命令Plot,ezplot。 Plot,ezplot 2、编写Matlab程序，利用二分法求非线性方程的 编写Matlab程序， Matlab程序 近似解。 近似解。 3、在Matlab软件中，用迭代法求解非线性方程的 Matlab软件中， 软件中 近似解。 近似解。 4、编写Matlab程序，利用牛顿法求非线性方程的近 编写Matlab程序， Matlab程序 似解。 似解。
解非线性方程实验
实验 2 解非线性方程实验实验目的：进一步熟练掌握解非线性方程的二分法算法、Steffensen算法、牛顿迭代法，提高编程能力和解算非线性方程问题的实践技能。
实验内容：用二分法算法(取[a,b]=[1,2])、Steffensen算法(取)、牛顿迭代法(取)解算非线性方程的根实验原理二分法算法Steffensen算法牛顿迭代法实验步骤1 要求上机实验前先编写出程序代码2 编辑录入程序3 调试程序并记录调试过程中出现的问题及修改程序的过程4 经反复调试后，运行程序并验证程序运行是否正确。
5 记录运行时的输入和输出。
实验思考试验改变Steffensen算法的迭代初值进行计算，试验改变牛顿迭代算法的迭代初值进行计算，分别对比计算结果并作出分析总结数学建模与应用：请建立数学模型并应用你的程序解决一个物理计算问题--木质球体浸入水中深度问题。
假若球体是由密度为ρ=0.638且半径为r=10cm 的长叶松木质球体构成的，问球体浸入水中深度是多少？又假若球体是由密度为ρ=0.638且三个长半轴分别为a=10cm,b=15cm,c=20cm的长叶松木质椭圆球体构成的，问球体浸入水中深度是多少？（建模可参考课本32页）实验报告书写要求：根据实验情况和结果撰写并递交实验报告。
实验报告应当有算法原理简介，算法流程图，程序代码，运行调试记录，计算结果阐释；对于数学建模的应用问题要有分析建模的过程、求解计算的过程、实际问题解答或阐释的内容。
实验总结(学会了......; 掌握了......; 训练了......; 发现了......; 今后学习中......有待提高。
）参考程序1.二分法1、实验原理：二分法2、实验步骤1建立函数function y=f3(x)y=x.^3-x-12建立函数function[c,err,yc]=bisect(f,a,b,delta)ya=feval(f,a);yb=feval(f,b);if ya*yb>0returnendmax1=1+round((log(b-a)-log(delta))/log(2)); for k=1:max1c=(a+b)/2;yc=feval(f,c)if yc==0a=c;b=c;elseif yb*yc>0b=c;yb=yc;elsea=c;ya=yc;endif b-a<deltabreakendendc=(a+b)/2;err=abs(b-a);yc=feval(f,c);3调用函数>> [c,err,yc]=bisect('f3',1,2,0.0001)y =-1y =5y =0.8750 yc =0.8750 y =-0.2969 yc =-0.2969 y =0.2246 yc =0.2246 y =-0.0515 yc =-0.0515 y =0.0826 yc =0.0826y =0.0146yc =0.0146y =-0.0187yc =-0.0187y =-0.0021yc =-0.0021y =0.0062yc =0.0062y =0.0020yc =0.0020y =-4.6595e-005yc =-4.6595e-005y =9.9479e-004yc =9.9479e-004y =4.7404e-004yc =4.7404e-004y =2.1371e-004yc =2.1371e-004y =8.3552e-005c =1.3247err =6.1035e-005yc =8.3552e-005 2、牛顿迭代1、实验原理：牛顿迭代法2、实验步骤1.建立函数function y=f(x)y=x.^3-x-12.建立函数function y=f2(x)y=3*x.^2-13.牛顿迭代法函数的程序function [p0,err,k,y]=Newton(f,df,p0,delta,epsilon,max1)for k=1,max1p1=p0-feval(f,p0)/feval(df,p0);err=abs(p1-p0);relerr=2*err/(abs(p1)+delta);p0=p1;y=feval(f,p0);if(err<delta)|(relerr<delta)|(abs(y)<epsilon),break,end end4.函数的调用>> [p0,err,k,y]=Newton('f','f2',1.5,0.001,0.0001,1000)y =0.8750y =5.7500y =0.1007y =0.1007y =4.4499y =0.0021y =0.0021y =4.2685y =9.2438e-07 p0 =1.3247err =4.8222e-04 k =3y =9.2438e-07三、Steffensen算法1、实验原理：Steffensen算法2、实验步骤①function y=f1(x)y=x.^3-x-1②function y=f2(x)y=3*x.^2-1③function [p,Q]=steffensen(f,df,p0,delta,epsilon,max1) R=zeros(max1,3);R(1,1)=p0;for k=1:max1for j=2:3nrdenom=feval(df,R(k,j-1));if nrdenom==0'division by zero in Newton-Raphson method'breakelesR(k,j)=R(k,j-1)-feval(f,R(k,j-1))/nrdenom;endaadenom=R(k,3)-2*R(k,2)+R(k,1);if aadenom==0'division by zero in Aitken’s Acceleration'breakelesR(k+1,1)=R(k,1)-(R(k,2)-R(k,1))^2/aadenom;endendif(nrdenom==0)|(aadenom==0)breakenderr=abs(R(k,1)-R(k+1,1));relerr=err/(abs(R(k+1,1))+delta);y=feval(f,R(k+1,1));if (err<delta)|(relerr<delta)|(y<epsilon)p=R(k+1,1)Q=R(1:k+1,:);breakendend④>> [p,Q]=steffensen('f1','f2',1.5,0.0001,0.001,100)y =5.7500y =-1y =-1p =p =Q =1.5000 0 00 0 0实验总结:学会了用二分法，牛顿迭代法，Steffensen算法来求解线性方程的解; 掌握了一定的解题方法，和建立函数的方法，在调用函数中出现的错误也能及时判断理解并加以改正; 训练了应用MATLAB的技巧和方法，在一定程度上得以提升; 发现了应用此方法的简易方便之处，也认识到自己的很多缺点要继续努力学习掌握好知识技能; 今后学习中不仅要掌握好理论知识理解并记忆，更要多多在实际中操作实验，灵活应用所学知识，来提高自己。
第二次实验报告(非线性方程求根)
数学实验报告实验五学院：数学与统计学院班级：数学与统计学院(2)班姓名：石紫雲学号：0234非线性方程的数值解法实验1实验目的1）进一步熟练掌握求解非线性方程的牛顿迭代法和弦截法。
2）根据牛顿迭代法和弦截法的原理，编写程序求解非线性方程，提高编程解决问题的能力。
2 实验内容（1）用牛顿法和割线法求下列方程的根x^2-e^x=0;x*e^x-1=0;（23实验原理（1）牛顿迭代公式：1()/'()k k k k x x f x f x +=- 双点弦法公式：111()()()()k k k k k k k f x x x x x f x f x +--=--- （2）令2()f x x A =-，再用牛顿法求根。
4实验步骤1)根据牛顿迭代法，双点弦法的算法编写相应的求根函数；2)用牛顿迭代法和双点弦法分别对方程进行求解；5 程序设计牛顿迭代法x0=;N=100;k=0;eps=5e-6;delta=1e-6;while(1)x1=x0-fc1(x0)/fc2(x0);k=k+1;if k>Ndisp('Newton method failed')breakendif(abs(x1-x0)<delta || abs(fc1(x1))<delta) break;endx0=x1;endfprintf('%f',x0)fprintf('%f',abs(fc1(x1)))双点弦法function cutline(x0,x1)N=100;k=0;delta=5e-8;while(1)(abs(x1-x0)>=delta)c=x1;x1=cutnext(x0,x1);x0=c;k=k+1;if k>Ndisp('Cutline method failed')break;endif(abs(x1-x0)<delta || abs(fc1(x1))<delta) break;endendfprintf('%10f\n',x1);function y=cutnext(a,b)y=b-fc(b)/(fc(b)-fc(a))*(b-a);1）原函数function fc1=fc1(x)fc1=x^2-exp(x);end导函数function fc2=fc2(x)fc2=2*x-exp(x);end2）原函数function fc1=fc1(x) fc1=x*exp(x)-1;end导函数function fc2=fc2(x)fc2=1*exp(x)+x*exp(x);end3）原函数function fc1=fc1(x)fc1=x^2-5;end导函数function fc2=fc2(x)fc2=2*x;end6实验结果及分析注：牛顿迭代法由于设置delta=1e-6，所以算出的误差e<*10^-6;割线法由于设置delta=5e-8，所以误差e<*10^-8.7总结。
数值分析实验报告——非线性方程求根
数值分析实验报告——非线性方程求根一、实验目的：1.掌握求解非线性方程的常用方法；2.了解非线性方程求根问题的数值解法；3.熟悉使用数值分析软件进行非线性方程求根的实现。
二、实验原理：非线性方程指的是形如f(x)=0的方程，其中f(x)是一个非线性函数。
非线性方程求根的常用方法包括二分法、割线法和牛顿法等。
其中，二分法是通过不断缩小区间范围来逼近方程的解；割线法是通过使用割线来逼近方程的解；牛顿法则是通过使用切线来逼近方程的解。
对于给定的非线性方程，可以根据实际情况选择合适的方法进行求根。
三、实验内容：1.编写求解非线性方程的函数，包括二分法、割线法和牛顿法；2.使用编写的函数求解给定的非线性方程，比较各个方法的收敛速度和精确程度；3.根据实际情况分析和选择合适的方法进行求根。
四、实验步骤：1.针对给定的非线性方程，编写二分法的函数实现：（1）首先确定方程的解存在的区间；（2）根据方程的解存在的区间，使用二分法逐步缩小区间范围；（3）根据设定的精度要求，不断循环迭代，直至满足要求或达到迭代次数限制；2.针对给定的非线性方程，编写割线法的函数实现：（1）首先需要确定方程的解存在的初始点；（2）根据方程的解存在的初始点，根据割线的定义进行迭代；（3）设定迭代的精度要求和限制次数，结束迭代；3.针对给定的非线性方程，编写牛顿法的函数实现：（1）首先需要确定方程的解存在的初始点；（2）根据方程的解存在的初始点，根据牛顿法的定义进行迭代；（3）设定迭代的精度要求和限制次数，结束迭代；4.根据给定的非线性方程，分别使用二分法、割线法和牛顿法进行求解，并比较各个方法的收敛速度和精确程度；5.分析实际情况，选择合适的方法进行求解。
五、实验结果：4.通过比较，发现割线法和牛顿法的收敛速度较快，精确程度较高，因此选择割线法进行求解。
六、实验总结：通过本次实验，我掌握了求解非线性方程的常用方法，并使用数值分析软件实现了二分法、割线法和牛顿法。
高数实验报告非线性方程组求解
高数实验报告一．实验问题.1.利用弦截法编程求求解x^5+x+-1=0的近似解；2.方程f(x)=x^2+x-4=0在（0,4）内有唯一的实根，现构造以下三种迭代函数（1）g1（x）=4-x^2,迭代初值为x0=4;(2)g2(x)=4/(1+x),迭代初值为x0=4;(3)g3(x)=x-(x^2+x-4)/(2x+1), 迭代初值为x0=4;3.分别取不同的参数值r，做迭代序数（11 -6），观察分叉和馄饨现象二．问题分析这三个题目中前两个题目都是利用弦截法和迭代法处理近似解，亦即利用无限替代与收敛判定达到求近似解的目的，所以按老师所讲授的方法类比操作即可。
第三个题目要求完成不同的参数r下迭代的实验步骤。
通过多次迭代变换以及收敛判定，易得实验结果及图像如下（四实验结果）中所示。
三．程序设计程序源代码1.弦截法f=inline('x^5+x-1');a=0;b=1;dlt=1.0e-5;x=1;k=1;while abs(f(x))>a*dltx=b-(b-a)*f(b)/(f(b)-f(a))a=x;k=k+1vpa(x,7)end2.迭代1f=inline('x^2+x-4');fzero(f,[0,4])a=0;b=4;dlt=1.0e-5;x=1;k=0;while k<4;x=4-a^2a=x;endans =1.5616k =1x =4k =2x =-12k =3x =-140k =4x =-19596迭代2f=inline('x^2+x-4'); fzero(f,[0,4])a=0;b=4;dlt=1.0e-5; x=1;k=0;while k<4;k=k+1x=4/(1+x)a=x;endans =1.5616k =1x =2k =2x =1.3333k =3x =1.7143k =4x =1.4737迭代3f=inline('x^2+x-4'); fzero(f,[0,4])a=0;b=4;dlt=1.0e-5; x=1;k=0;k=k+1x=x-(x^2+x-4)/(2*x+1)a=x;endans =1.5616k =1x =1.6667k =2x =1.5641k =3x =1.5616k =4x =1.56163.步骤1clear;for r=0:0.3:3.9;x=0.9;for i=1:150x=r*x*(1-x);x1(i)=i;y(i)=x;endfor k=1:50xx(k)=r; yy(k)=y(100+k); endhold on; grid on;plot(xx,yy,'ko')end步骤3clear;for j=270:390;x=0.9; r=j/100;for i=1:100x=r*x*(1-x);x1(i)=i;y(i)=x;endfor k=1:50xx(k)=r; yy(k)=y(50+k);endhold on;grid on;plot(xx,yy,'ko')end四．实验结果实验输出图像1k =2ans =.5000000k =3ans =.6595745k =4ans =.7199516k =5ans =.7422761k =6ans =.7503618k =7ans =.7532636k =8ans =.7543013k =9ans =.7546719k =10ans =.7548042k =11ans =.7548515k =12ans =.7548683k =13ans =.7548743k =14ans =.7548765由实验结果可知，弦截法比牛顿收敛法的运算速度慢，效率较低。
实验一非线性方程组求解实验报告
计算方法实验报告专业班级：姓名：学号：实验成绩：1．【实验题目】非线性方程组求解2．【实验目的】（1）．掌握二分法、迭代法、牛顿迭代法求方程近似根的基本思想与原理。
（2）．掌握常用迭代算法的程序实现。
3．【实验内容】迭代法是求解非线性方程的基本方法，其构造方法可以有多种多样，但关键是怎样才能使迭代收敛且有较快的收敛速度。
考虑一个简单的代数方程，针对该方程，可以构造多种迭代法，如：取初始值，取，分别用以上迭代格式作实验，记录各算法的迭代过程4. 【实验要求】（1）取定某个初始值，按方案1～3对非线性方程求根，它们的收敛性如何？重复选取不同的初始值，反复实验。
请读者自行设计一种比较形象的记录方式（如利用Matlab的图形功能），分析三种迭代法的收敛性与初值选取的关系。
（2）对三个迭代格式的某一种，分别取不同的初始值进行迭代，结果如何？试分析迭代法对不同的初值是否有差异？（3）对代数方程，分别用方案1 用二分法求解；方案2 用牛顿法求解；5. 【算法描述】二分法算法步骤1）计算有根区间的端点a,b及预先给定的精度e。
2）计算中点(a+b)/2。
3）若f(x)f(a)<0,则x b,转向4)；否则，x a，转向4).⇒⇒4）若b-a<e，则输出满足精度的根x,结束；否则转向2)。
牛顿法迭代法的计算步骤x01）给出初始近根及精度e。
2）计算。
x x x x f f 1000)(')(⇒-3）若|-|<e ，则转向4)；否则转向2)。
x 1x 0x x 01⇒4）输出满足精度的根，结束。
x 16. 【源程序（带注释）】二分法#include<stdio.h>#include<math.h>#include<conio.h>#include<windows.h>float f(float x){float a;a=x*x*x-x-1;return a;} /*求函数值，如果求其它函数，只需改成其它函数即可*/ main(){float a,b,e,x; /* a,b 分别表示有根区间的左、右端点, e 是精度要求,x 区间中点值*/system("CLS");//清屏printf("对代数方程x^3-x-1=0，分别用\n 方案1 用二分法求解\n");printf(" \n please input data a =");scanf("%f",&a);printf(" \n please input data b=");scanf("%f",&b);if(f(a)*f(b)<0){while(f(x)!=0){x=(a+b)/2;if(f(x)*f(a)<0){b=x;if(fabs(b-a)<0.000001)break;elsecontinue;}else{a=x;if(fabs(b-a)<0.000001)break;else continue;}}printf("\n");x=(b+a)/2;printf("the root of f(x)=0 is x=%f\n",x);}elseprintf("\ not root! afresh input\n"); /*表示[a,b] 区间无根，重新选择有根区间*/getch();return(x);}牛顿法#include<stdio.h>#include<math.h>#include<conio.h>#include<windows.h>#define maxrept 1000 /*最大迭代次数*/float f(float x) {float a;a=x*x*x-x-1;return a; /*函数f(x) */}float df(float x) {return(1+exp(-x)); /* 函数f(x)的导数) (x f ′*/ }float iterate(float x) {float x1;x1=x-f(x)/df(x); /* 牛顿迭代函数iterate(x)=x-f(x) / ) (x f ′*/return(x1);}main() {float x0,x1,d;int k=0;//clrscr();system("CLS");printf("对代数方程x^3-x-1=0，分别用\n方案2 用牛顿法求解\n");printf("\n please input x0="); /* 输入迭代初值x0 */scanf("%f",&x0);printf("\n k xk\n");printf("\ %d %f\n",k,x0);do {k++;x1=iterate(x0);printf(" %d %f\n",k,x1);d=fabs(x1-x0);x0=x1;}while((d>=0.000001)&(k<maxrept));if(k<maxrept)printf("the root of f(x)=0 is x=%f, k=%d\n",x1,k);elseprintf("\n the iteration is failed!\n");getch();}7．【实验结果与分析总结（含运行结果截图）】。
科学计算与数学建模实验报告牛顿法求解非线性方程
科学计算与数学建模实验报告牛顿法求解非线性方程一、引言非线性方程是数学中的一个重要研究内容，其求解方法有很多，其中之一就是牛顿法。
牛顿法是一种迭代方法，通过不断逼近函数的零点来求解非线性方程。
在本实验中，我们将使用牛顿法来求解给定的非线性方程，并验证其有效性。
二、实验方法1.确定问题：给定非线性方程f(x)=0，需要求解方程的根。
2.初始化：选择一个初始解x_0，并给定停止准则，如迭代次数、函数误差等。
3.迭代计算：a)计算函数f(x)在x_i处的导数f'(x_i)。
b)利用牛顿迭代公式进行迭代计算：x_{i+1}=x_i-f(x_i)/f'(x_i)。
c)检查迭代终止条件，若满足条件则停止迭代，否则返回步骤a)继续迭代。
4.输出结果：输出迭代过程中的迭代次数和解x。
三、实验结果我们选择一个较为简单的非线性方程f(x)=x^2-2来进行牛顿法求解。
初始解选取为x_0=1，停止准则为函数误差小于等于0.0001根据上述计算方法，我们进行迭代计算，并记录迭代次数和解x的变化情况。
具体结果如下表所示：迭代次数解x-----------------11.521.416731.414241.4142(收敛)从表中可以看出，当迭代4次时，解x已经收敛于1.4142，符合停止准则，因此我们可以认为此时已经找到了方程的根。
四、实验讨论通过上述实验可以发现，牛顿法是一种有效的求解非线性方程的方法。
它利用了函数在特定点处的导数的信息来逼近函数的零点，从而实现了迭代计算。
同时，牛顿法的收敛速度比较快，迭代次数较少，可以在较短的时间内找到方程的根。
然而，牛顿法也存在一些不足之处。
首先，它对初始解的选择较为敏感，不同的初始解可能导致迭代结果的差异。
其次，牛顿法可能出现发散现象，即迭代过程无法收敛到方程的根。
因此，对于一些复杂的非线性方程，我们需要选择合适的方法来求解。
五、总结通过本次实验，我们了解了牛顿法求解非线性方程的基本过程，并验证了其有效性。
非线性方程求解数值分析上机实验报告
实验报告一题目：非线性方程求解摘要：非线性方程的解析解通常很难给出，因此线性方程的数值解法就尤为重要。
本实验采用两种常见的求解方法二分法和Newton法及改进的Newton法。
前言：（目的和意义）掌握二分法与Newton法的基本原理和应用。
数学原理：对于一个非线性方程的数值解法很多。
在此介绍两种最常见的方法：二分法和Newton 法。
对于二分法，其数学实质就是说对于给定的待求解的方程f(x)，其在[a,b]上连续，f(a)f(b)<0，且f(x)在[a,b]内仅有一个实根x*，取区间中点c，若，则c恰为其根，否则根据f(a)f(c)<0是否成立判断根在区间[a,c]和[c,b]中的哪一个，从而得出新区间，仍称为[a,b]。
重复运行计算，直至满足精度为止。
这就是二分法的计算思想。
Newton法通常预先要给出一个猜测初值x0，然后根据其迭代公式产生逼近解x*的迭代数列{x k}，这就是Newton法的思想。
当x0接近x*时收敛很快，但是当x0选择不好时，可能会发散，因此初值的选取很重要。
另外，若将该迭代公式改进为其中r为要求的方程的根的重数，这就是改进的Newton法，当求解已知重数的方程的根时，在同种条件下其收敛速度要比Newton法快的多。
程序设计：本实验采用Matlab的M文件编写。
其中待求解的方程写成function的方式，如下function y=f(x);y=-x*x-sin(x);写成如上形式即可，下面给出主程序。
二分法源程序：clear%%%给定求解区间b=1.5;a=0;%%%误差R=1;k=0;%迭代次数初值while (R>5e-6) ;c=(a+b)/2;if f12(a)*f12(c)>0;a=c;elseb=c;endR=b-a;%求出误差k=k+1;endx=c%给出解Newton法及改进的Newton法源程序：clear%%%% 输入函数f=input('请输入需要求解函数>>','s')%%%求解f(x)的导数df=diff(f);%%%改进常数或重根数miu=2;%%%初始值x0x0=input('input initial value x0>>');k=0;%迭代次数max=100;%最大迭代次数R=eval(subs(f,'x0','x'));%求解f(x0)，以确定初值x0时否就是解while (abs(R)>1e-8)x1=x0-miu*eval(subs(f,'x0','x'))/eval(subs(df,'x0','x'));R=x1-x0;x0=x1;k=k+1;if (eval(subs(f,'x0','x'))<1e-10);breakendif k>max;%如果迭代次数大于给定值，认为迭代不收敛，重新输入初值ss=input('maybe result is error,choose a new x0,y/n?>>','s');if strcmp(ss,'y')x0=input('input initial value x0>>');k=0;elsebreakendendendk;%给出迭代次数x=x0；%给出解结果分析和讨论：1.用二分法计算方程在[1，2]内的根。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一.实验目的：
通过本节实验课的学习，要求我们理解并掌握二分法、不动点迭代、牛顿切线法及弦截法解非线性方程求根的原理，掌握相应的算法原理，通过计算机解决实验问题
二.实验内容：
1、用对分区间法方程
1-x-sinx=0
在区间[0,1]上的误差小于10^(-4)的一个根，并记录对方区间的次数。
2、用不动点迭代法求解方程下
x-log(x)=2(x>1)
要求相对误差容限e=10^(-8)。
3、用Newton法求方程
x^3-x-1=0
在区间[-3,3]上的误差不大于10^(-5)的根，分别取初值x0=1.5, x0=0, x0=-1进行计算，比较他们的迭代次数。
三. 实验方案（程序设计说明）[包括算法设计思路，必要的流程图，界面设计
说明、使用模块及变量的说明等。
]
1、二分法是对区间收索法的一种改进，具体做法为：先求一区间的中点，并计算其
函数值，若恰好有函数值为0，就是方程的根，若不为0，在判断此点的函数值与两端的函数值乘积的情况，取小于0的那个端点在进行上述对分，直到满足要求为止。
2、迭代法分为两种，一种是从任何可取的初值出发都能保证收敛，称之为大范围收
敛的方法。
另一类称之为局部收敛法，即为了保证收敛必须选取初值充分接近于所要求的解。
迭代法的基本思想是一种逐渐逼近的方法，首先给定一个粗造的初值，然后用一个迭代公式，反复矫正这个初值，直到满足预先给出的精确要求为止。
3、双点弦接法与Newton法不同，两者有本质的区别，它分为两步，不属于不动点
迭代法。
四. 实验步骤或程序（经调试后正确的源程序）
（填写主要步骤与程序代码等，不够可附页）
1、f=inline('x+sin(x)-1');
a=0;
b=1;
dlt=1.0e-4;
k=1;
while abs(b-a)>dlt
c=(a+b)/2;
if f(c)==0
break;
elseif f(c)*f(b)<0
a=c;
else b=c;
end
fprintf('k=%d,x=%.5f\n',k,c); k=k+1;
end
2、
eps=10^(-8);dx=1;x0=3.5;k=0;
while(dx>eps)
k=k+1;
x=log(x0)+2;
dx=abs(x-x0)/(1+abs(x));
x0=x;
end
k
x
3、
f=inline('x^3-x-1');
df=inline('3*x^2-1');
d2f=inline('6*x');
a=-3;
b=3;
dlt=1.0e-5;
if f(a)*d2f(a)>0
x0=a;
else
x0=b;
end
m=min(abs(df(a)),abs(df(b)));
k=0;
while abs(f(x0))>m*dlt
k=k+1;
x1=x0-f(x0)/df(x0);
x0=x1;
fprintf('k=%d x=%.5f\n',k,x0); end
for x0=1.5
fprintf('k=%d x=%.5f\n',k,x0); end
for x0=0
fprintf('k=%d x=%.5f\n',k,x0); end
for x0=-1
fprintf('k=%d x=%.5f\n',k,x0); end
五．程序运行结果：
1、k=1,x=0.50000
k=2,x=0.75000
k=3,x=0.62500
k=4,x=0.56250
k=5,x=0.53125
k=6,x=0.51563
k=7,x=0.50781
k=8,x=0.51172
k=9,x=0.50977
k=10,x=0.51074
k=11,x=0.51123
k=12,x=0.51099
k=13,x=0.51086
k=14,x=0.51093
2、k =
15
x =
3.1462
3、k=1 x=-2.03846
k=2 x=-1.39028
k=3 x=-0.91161
k=4 x=-0.34503
k=5 x=-1.42775
k=6 x=-0.94242
k=7 x=-0.40495
k=8 x=-1.70690
k=9 x=-1.15576
k=10 x=-0.69419
k=11 x=0.74249
k=12 x=2.78130
k=13 x=1.98273
k=14 x=1.53693
k=15 x=1.35726
k=16 x=1.32566
k=17 x=1.32472
当x0=1.5时：
k=17 x=1.50000
当x0=0时：
k=17 x=0.00000
当x0=-1时：
k=17 x=-1.00000
2、六．实验总结：通过实验学会理解并掌握二分法、不动点迭代、牛顿切线法及弦截法解非线性方程求根的原理，掌握相应的算法原理，通过计算机解决实验问题并通过反复的上机实验操作，解决了在实验过程中遇到的实验问题，并了解了一些函数的特殊用法，学会了用这三种基本方法解决实际遇到的问题，并了解了二分法、不动点迭代、牛顿切线法及弦截法的各种变形算法。