计算方法实验报告2.docxVIP

计算方法实验报告（四） 方程和方程组的迭代解法 一、实验问题 利用简单迭代法，两种加速技术，牛顿法，改进牛顿法，弦割法求解习题5-1，5-2，5-3中的一题，并尽可能准确。 选取5-3：求x3-x2-1=0在x=1.5附近的根。 二、问题的分析（描述算法的步骤等） （1）简单迭代法算法： 给定初始近似值p0,求p=φp的解。 Step 1 令i=0； Step 2 令pi+1=φpi(计算pi+1)； Step 3 如果pi+1=pi,则迭代终止，否则重复Step 2。 （2）Aitken加速法算法 Step 1 令k=0，利用简单迭代算法xk+1=φxk得到迭代序列xk； Step 2 令xk*=xk-xk-xk-12xk-2xk-1+xk-2(计算xk*得到一个新的序列xk*，其中k=0,1,2…)； Step 3 如果xk+1*=xk*,则迭代终止，否则重复Step 2。 (3)插值加速法算法 Step 1 令k=0，利用简单迭代算法xk+1=φxk得到迭代序列xk； Step 2 令xk*=xk+xk-xk-1xk-xk+1xk-1-2xk+xk+1(计算xk*得到一个新的序列xk*，其中k=1,2,3…)； Step 3 如果xk+1*=xk*,则迭代终止，否则重复Step 2。 (4)牛顿法算法 Step 1给定初始近似值x0； Step 2令xk+1=xk-fxkfxk,其中k∈N,计算得到xk的 序列； Step 3如果xk+1=xk,则迭代终止，否则重复Step 2。 （5）改进牛顿法的算法 Step 1给定初始近似值x0； Step 2令xk+1=xk-2fxkfxk+sgnfxkf2xk-2fxkfxk ,其中k∈N,迭代计算得到xk的 序列； Step 3如果xk+1=xk,则迭代终止，否则重复Step 2。 （6）弦割法算法（双点弦割法） Step 1给定初始近似值x0，x1； Step 2令xk+1=xk-fxkxk-xk-1fxk-fxk-1其中k∈N,计算得到xk的 序列； Step 3如果xk+1=xk,则迭代终止，否则重复Step 2。 三、程序设计 (1)简单迭代法 利用迭代公式x=31+x2进行迭代运算。 #include iostream.h #include math.h #includestdio.h double fun(double x) { double c=1+x*x; return pow(c,1/3.0); } void main() { double x=1.5; double y=0; double D=1; double e=0.001; while(De) { D=0; y=fun(x); if(fabs(y-x)=D) { D=fabs(y-x); } x=y; } coutxendl; } （2）牛顿法： 利用公式xk+1=xk-fxkfxk 进行迭代运算 程序设计如下： #include iostream.h #include math.h double fun(double x) { double a=2*pow(x,3.0)-pow(x,2.0)+1; double b=3*pow(x,2.0)-2*x; return a/b; } void main() { double x=1.5; double y=0; double D=1; double e=0.001; double f=0; while(De) { D=0; y=fun(x); if(fabs(y-x)=D) { D=fabs(y-x); } x=y; f++; } coutxendl; coutf=fendl; } （3）运用改进牛顿法： 迭代公式：xk+1=xk-fxkfxk -f(xk)2f3(xk)f2(xk) 程序代码如下： #include iostream.h #include math.h double fun(double x) { double a=2*pow(x,3.0)-pow(x,2.0)+1; double b=3*pow(x,2.0)-2*x; double c=pow((pow(x,3.0)-pow(x,2.0)-1),2.0); double d=(6*x-2)/12; return a/b-c*d; } void main() { double x=1.5; double y=0; double D=1; double e=0.001; double f=0;