数值计算方法实验报告---穆明国.docx

数值计算方法实验报告 学院：信息学院 班级：计科高职12级2班 学号：201203014080 姓名：穆明国 基本情况 ? 课程名称：计算方法III ? 使用教材：《数值计算方法》 ? 实验环境：装有C语言或装有其它编程语言系统的系统的微型计算机 ? 实验要求： 在上机实验时完成相应实验的算法的程序编制，并上机运行，学会应用这些算法于实际问题，以便对算法有更进一步的认识和理解。考察和体会数值计算中出现的一些问题和现象： 误差的估计,算法的稳定性、收敛性、收敛速度以及迭代初值对收敛的影响等。 实验内容 实验一：解线性方程组的迭代法 1、目的与要求： （1）理解向量范数、矩阵范数、谱半径等的有关概念 （2）掌握常用的几种迭代格式 （3）掌握迭代收敛性及误差估计 2、实验程序： （1）Jacobi迭代方法： #include stdio.h #include math.h int main() { float a[3][3]={{10,-1,-2},{-1,10,-2},{-1,-1,5}},b[3]={7.2,8.3,4.2}; float x[3]={0,0,0},sum; int i,j,k,n=3; printf(\t\t X[1]\t\t X[2]\t\t X[3]\n); for(k=0;k8;k++) { for(i=0;i3;i++) { sum=0; for(j=0;jn;j++) { if(i==j)continue; sum=sum+a[i][j]*x[j]; } x[i]=(b[i]-sum)/a[i][i]; } printf(第%d次迭代：\t,k+1); for(i=0;in;i++) { printf(%f\t,x[i]); } printf(\n); } } （2）Seidel迭代方法： #include iostream #include iomanip using namespace std; int main() { int i,j,k=0,m,n; double t1,t2,e1,e2=0.0; cout请输入精度e：; cine1; cout请输入系数矩阵行数：; cinm; cout请输入系数矩阵列数：; cinn; coutendl; double (**a)=new double *[m]; for(i=0;i=m;i++) { a[i]=new double[n]; } double (*b)=new double [m]; double (*x)=new double [n]; cout请输入系数矩阵：endl; cout-------------------------------------------------------------endl; for(int num1=0;num1m;num1++) { for(int num2=0;num2n;num2++) { cina[num1][num2]; } coutendl; } cout输入的系数矩阵为：endl; for (int num3=0;num3m;num3++) { for(int num4=0;num4n;num4++) { couta[num3][num4] ; } coutendl; } cout请输入矩阵b：endl; for(int num5=0;num5m;num5++) { cinb[num5]; } cout输入的矩阵b为：endl; for(int num6=0;num6m;num6++) { coutb[num6] ; } coutendl; for(int num7=0;num7n;num7++) { x[num7]=0.0000; } do { cout第k次迭代值：; e2=0.0; for(i=0;im;i++) { double sum=0.0; for(j=0;jn;j++) { if(j!=i) sum+=a[i][j]*x[j]; } t1=x[i]; t2=e2; x[i]=(b[i]-sum)/a[i][i]; e2=(x[i])-t1=0?(x[i])-t1:t1-(x[i]); e2=(e2=t2?e2:t2); coutsetprecision(8)x[i] ; } coutendl; k++; }while(e2=e1k30); cout共迭代了k次; delete[]a; delete[]b; delete[]x; return 0 ;} 3、运行结果： （1）Jacobi迭代方法： （2）Seidel迭代方法： 实验二：解线性方程组的直接法 1