数值分析实验三(解线性方程组).doc

《数值分析》实验报告 实验编号：实验三 课题名称：解线性方程组 一、算法介绍 1、 定义四个函数分别为：DieDai()，Newton()，XianWei()，DuiFen()。在主函数中输入要选用方法所对应的序号后，用switch语句对函数进行调用。 2、 迭代法的主要思想为：保留一个变量在等号的左边，其他都移到等号右边。 3、 Newton法的主要算法为：把f(x)在x0附近展开成Taylor级数，取其线性部分，选取一点x0，该点所对应的值为f(x0)，对于n=1,2,…,Nmax，按Xn+1=Xn-f(Xn)/f’(Xn)求出Xn+1，并计算f(Xn+1)，若|Xn+1-Xn|小于容许误差，则停止计算。 4、 弦位法：选定初始值x0,x1，并计算f(x0)，f(x1)；按迭代公式xn+1=xn-f(xn)(xn-xn-1)/(f(xn)-f(xn-1))计算x2，再求f(x2)；如果相邻两次迭代值之差在容许的误差范围内则迭代停止，否则用(x2,f(x2)),(x1,f(x1))分别代替(x1,f(x1)),(x0,f(x0))，重复前两个步骤，直至相邻两次迭代值之差在容许的误差范围内。 5、 对分区间法： 选取方程的根所在的区间(a,b)，取其中点c代入方程中得其值为f(c)，如果f(c)与f(a)异号说明方程的根在(a,c)区间中，则令b=c，否则令a=c，如果f(c)的绝对值小于0.0001，则停止运算，否则继续计算，直至f(c)的绝对值小于0.0001。 二、程序代码 #include iostream #include iomanip #include cmath using namespace std; double f(double x) { x=x*x*x-3*x-1; return x; } void DieDai() { cout迭代序列：\n; double x0=2,x1=0; int i=1; while(fabs(x0-x1)=0.0001) { x1=x0; x0=pow((3*x0+1),1.0/3); coutxi=setiosflags(ios::fixed)setprecision(4)x0 ,f(xi)=f(x0),误差为：x0*x0*x0-3*x0-1endl; i++; } cout方程近似解x*=x0endl; cout共进行i-1次迭代\n ; } void Newton() { cout迭代序列：\n; double x0=2,x1=0; int i=1; while(fabs(x0-x1)=0.0001) { x1=x0; x0=x0-f(x0)/(3*x0*x0-3); coutxi=setiosflags(ios::fixed)setprecision(4)x0 ,f(xi)=f(x0),误差为：x0*x0*x0-3*x0-1endl; i++; } cout方程近似解x*=x0endl; cout共进行i-1次迭代\n ; } void XianWei() { double x0=1,x1=3,x2=2; cout选定曲线y=f(x)上的两个点P0(x0,f(x0))和P1(x1,f(x1))\n; int i=2; while(fabs(f(x2))=0.0001) { x2=x1-f(x1)*(x1-x0)/(f(x1)-f(x0)); coutsetiosflags(ios::fixed)setprecision(4)当前区间(; coutx0,x1),与x轴交点(x2,f(x2)),误差为：; if(f(x2)*f(x0)0) { coutx2*x2*x2-3*x2-1endl; x1=x2; } else { coutx2*x2*x2-3*x2-1endl; x0=x2; } i++; } cout方程近似解x*=x2endl; c