东北大学软件学院数值分析程序实践报告

数值分析程序实践实验报告姓 名学 号班 级指 导 教 师阎家斌开 设 学 期2013-2014第二学期开 设 时 间第8周第13周实验题目用SOR法求解方程组实验日期2014.4.15-2014.5.30评 定 成 绩评定人签字评 定 日 期东北大学软件学院实验目的:用SOR法求解方程组：分别取w=0.65、1、1.2、1.45计算.要求精度为10-6；并指出迭代次数。实验代码和结果：#include#include#include / fabs();using namespace std;void SolveEquation(const int &dim) double *Arr = new doubledim*dim; /方程矩阵 double *B = new doubledim; /方程右端系数 double *X = new doubledim; /初始解 double *dX = new doubledim; double e; /解的精度 double w; /w=1时，为GS迭代 double S, S1; const int N = 5000; /最大迭代次数,防止死机 int n=0; bool hasSolved = false; cout endl 初始化方程组左端系数矩阵! endl; for (int i=0; idim; i+) for (int j=0; jdim; j+) cout Arr i j Arri*dim+j; cout 初始化方程组右端系数矩阵! endl; for (i=0; idim; i+) cout B i Bi; system(cls); cout 方程组如下: endl endl; for (i=0; idim; i+) for (int j=0; jdim; j+) cout setw(8) Arri*dim+j X j+1; if (j dim-1) cout + ; else if (j = dim-1) cout = ; cout setw(8) Bi endl; cout endl 设置解的精度要求:(比如1e-6) e; cout endl; cout 输入初始向量! endl; for (i=0; idim; i+) cout X i+1 Xi; cout endl 设置参数 ( 02 )! w; for (int k=0; kN; k+) /3-8 S = 0.0; for (i=0; idim; i+) /5-7 double sum = 0.0; for (int j=0; jdim; j+) sum += Arri*dim+j*Xj; if (fabs(Arri*dim+i) 1e-10) cout 被零除! endl Arri*dim+i fabs(S) /fabs()取绝对值函数 S = S1; Xi = Xi + S1; n+; if (fabs(S) = e) cout迭代次数为：endl; coutnendl; hasSolved = true; /system(cls); cout 线性方程组的解为： endl; for (int t=0; tdim; t+) cout X t+1 = Xt t; cout endl; break; if (!hasSolved) cout 初始向量或 不合理, 请重新设置! endl; delete Arr; delete B; delete X; del