實验三四种解方程组的方法.docVIP

《数值分析》实验报告 实验序号：实验三 题目名称: 四种方法解方程 1、实验目的： 用对分区间法、迭代法、Newton法、弦位法求在附近的根，准确到4位有效数字。 算法分析 对分区间法： 算法思想： (a)找出 的根的存在区间（ab），并计算出端点处的函数值f(a),f(b)(b)计算在区间中点的值 (c)判断：若 ，则停止计算。否则，，。 (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++; } } 迭代法： 算法思想： ()选取初值 ，并确定方程 的等价形式。(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) 选定，计算 。 (b) 对于…,Nmax,用’(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; } } 弦位法： 算法思想： () 选定初始值并计算和 ； () 用迭代公式计算，再求； () 判别：如果则迭代停止；否则，用和分别代替和。重复计算()(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=kendl; break; } else { if((t0n0m0)||(t0n0m0)) b=k; if((t0n0m0)||(t0n0m0)) a=k; i++; cout对分后的区间为:(a,b)endl; cout第i步结果是： kendl; }} i=0; } void Diedai() { double a,b=2.0; cout请输入要求的精度:endl; cinj; while(1) { a=pow(3.0*b+1,1/3.0); coutdxi=a-bendl; if(abs(a-b)j||abs(a-b)==0) { i++; cout第i步结果是:aendl; cout该方程的解为：X=aendl; break; } else { b=a; i++; cout第i步结果是：aendl; } } i=0; } void Newton() { double a,b,t,x0=2,