实验3：4种解方程组的方法.docVIP

《数值分析》实验报告 实验序号：实验三 题目名称: 四种方法解方程 学号: 姓名: 任课教师: 马季骕 专业班级:计算机科学与技术（非师范） 1、实验目的： 用对分区间法、迭代法、Newton法、弦位法求X^3-X^2-1=0在X0=1.5附近的根，准确到4位有效数字。 2、算法分析： 对分区间法： 算法思想： (a)找出 f(x)=0的根的存在区间（a,b），并计算出端点处的函数值f(a),f(b)； (b)计算f(x)在区间中点的值f((a+b)/2)； (c)判断：若f((a+b)/2)≈0 ，则停止计算。否则，若f((a+b)/2)与发f（a）异号，则根位于（a，(a+b)/2）中，以(a+b)/2代替b；若f((a+b)/2)与发f（b）异号，则根位于（(a+b)/2，b）中，以(a+b)/2代替a。 (d)重复步骤(b)(c)，直到区间缩小到容许的误差范围之内。此时，区间中点即可作为所要求的根。 实验代码： while(1) { n=a*a*a-3*a-1; m=b*b*b-3*b-1; k=(a+b)*0.5; ?????????? t=k*k*k-3*k-1; if((t=jt0)||(pow(0.5,i))=j||t==0)? { i++; cout该方程的解为：X=kendl; break; } else????????????? { if((t0n0m0)||(t0n0m0)) b=k; if((t0n0m0)||(t0n0m0)) a=k; i++; } } 迭代法： 算法思想： (1)选取初值 ，并确定方程f（x）=0 的等价形式x=φ（x）。 (2)计算x1=φ（x0） 。 (3)判别：如果|x1-x0|ε 则停止计算;否则用x1代替x0，重复步骤(2)和(3)（可以预设最多重复次数）。 实验代码： while(1) { a=pow(3.0*b+1,1/3.0);?? if(abs(a-b)j||abs(a-b)==0)? { i++; cout该方程的解为：X=aendl; break; } else????? { b=a; i++; } } Newton法： 算法思想： (a) 选定X0，计算f(X0) 。 (b) 对于n=0,1,2,…,Nmax,用Xn+1=Xn-f(Xn)/f’(Xn)求Xn+1，并计算f(Xn+1).若f(Xn+1)≈0或|Xn+1-Xn|ε（容许误差），则停止计算。 实验代码： while(1) { a=x0*x0*x0-3*x0-1;?????? b=3*x0*x0-3;??????????? x1=x0-(a/b); t=x1*x1*x1-3*x1-1; if(t0tj||abs(x1-x0)j||t==0)? { i++; cout该方程的解为：X=x1endl; break; } else { i++;?? //迭代 x0=x1; } } 弦位法： 算法思想： (1) 选定初始值X0，X1，并计算f(X0)和f(X1) ； (2) 用迭代公式计算X2，再求f(X2)； (3) 判别：如果f(X2) ≈则迭代停止；否则，用(X2, f(X2))和(X1, f(X1))分别代替(X1, f(X1))和(X0, f(X0))。重复计算(2)(3)，直至相邻两次迭代值之差在容许范围之内。 实验代码： while(1) { a=x0*x0*x0-3*x0-1; b=x1*x1*x1-3*x1-1; x2=x1-(b/(b-a))*(x1-x0);?? t=x2*x2*x2-3*x2-1; if(t0tj||t==0||abs(x1-x0)j)? { i++; cout该方程的解为：X=x2endl; break; } else??????? { x0=x1; x1=x2; i++; } } 3、程序代码 #includeiostream #includecmath using namespace std; int i=0; double j; void Duifen() { double a,b,m,n,k,t; cout请输入一个区间：endl; cinab; cout请输入要求的精度:endl; cinj; cout所求解的区间为:(a,b)endl; while(1) { n=a*a*a-3*a-1; m=b*b*b-3*b-1; k=(a+b)*0.5; t=k*k*k-3*k-1; coutdxi=b-aendl; if((t=jt0)||(pow(0.5,i))=j||t==0) { i++; cout对分后的区间为:(a,b)endl; cout第i步结果是： kendl; cout该方程的解为：X=k