数值分析上机实践报告.docVIP

数值分析上机实践报告 班级：计算机1002姓名：陈斯琪 学号 课题三 实验题目：线性方程组的迭代法 实验要求 应用迭代法求解线性方程组，并与直接法作比较； 分别对不同精度要求，如，利用所需迭代次数体会该迭代法的收敛快慢； 对方程组（2），（3）使用SOR方法时，选取松弛因子=0.8,0.9,1,1.1,1.2等，试观察对算法收敛性的影响，并找出你所选用松弛因子的最佳值； 编制出各种迭代法的程序并给出计算结果。 目的和意义 通过上机了解迭代法求解线性方程组的特点；掌握求解线性方程组的各类迭代法； 体会上机计算时，终止准则‖X^(k+1)-X^k ‖∞ε，对控制迭代精度的有效性； 体会初始值和松弛因子的选择，对迭代收敛速度的影响 实验方程组 （1）线性方程组 = 精确解. (2) 对称正定线性方程组 精确解. （3）三对角线性方程组 精确解. 实验程序代码及截图 应用Jacobi迭代法求解方程组 代码如下： #includeiostream.h #includemath.h #define N 10 //十阶矩阵 static double A[N][N]={4,2,-3,-1,2,1,0,0,0,0,8,6,-5,-3,6,5,0,1,0,0,4,2,-2,-1,3,2,-1,0,3,1,0,-2,1,5,-1,3,-1,1,9,4,-4,2,6,-1,6,7,-3,3,2,3,8,6,-8,5,7,17,2,6,-3,5,0,2,-1,3,-4,2,5,3,0,1,16,10,-11,-9,17,34,2,-1,2,2,4,6,2,-7,13,9,2,0,12,4,0,0,-1,8,-3,-24,-8,6,3,-1};//方程组左侧系数矩阵 static double B[N]={5,12,3,2,3,46,13,38,19,-21}; //右侧值 static double Y[N]; //输出比较项static double Y[N]; static double X[N]; //输出项 static double G[N]; //X = BX + G的G矩阵 int i,j,k; //计数器 double eps; int M=100; bool distance() { //求两输出项的差的范数是否满足精度要求 double temp=0; for (i=0;iN;i++) { temp=temp+fabs(X[i]-Y[i]); } if (tempeps) return false; else return true; //满足精度要求则结束程序 } void main() { cout最大迭代次数为100次endl; cout你希望的精度是多少？endl; couteps=; cineps; //形成迭代矩阵B，存放到A中 for (i=0;iN;i++) { if (fabs(A[i][i])eps) { cout 打印失败endl; return; } double T=A[i][i]; for (j=0;jN;j++) { A[i][j]=-A[i][j]/T; } A[i][i] = 0; G[i]=B[i]/T; } int counter=0; while (counterM) { //迭代 for (i=0;iN;i++) { double temp=0; for (j=0;jN;j++) { temp=temp+A[i][j]*Y[j]; } X[i]=G[i]+temp; } if (distance()==true) break; else { //交换X，Y向量； for(i=0;iN;i++) { Y[i]=X[i]; } } counter++; } //打印X cout 迭代次数为：counter次。该线性方程组的解（X1，X2，X3......X10）为：endl; coutX1=X[1]endl; coutX2=X[2]endl; coutX3=X[3]endl; coutX4=X[4]endl; coutX5=X[5]endl; coutX6=X[6]endl; coutX7=X[7]endl; coutX8=X[8]endl; coutX9=X[9]endl;