東北大學(xué)軟件學(xué)院數(shù)值分析程序?qū)嵺`報告

數(shù)值分析程序?qū)嵺`實驗報告姓 名學(xué) 號班 級指 導(dǎo) 教 師閻家斌開 設(shè) 學(xué) 期2013-2014第二學(xué)期開 設(shè) 時 間第8周第13周實驗題目用SOR法求解方程組實驗日期2014.4.15-2014.5.30評 定 成 績評定人簽字評 定 日 期東北大學(xué)軟件學(xué)院實驗?zāi)康?用SOR法求解方程組：分別取w=0.65、1、1.2、1.45計算.要求精度為10-6；并指出迭代次數(shù)。實驗代碼和結(jié)果：#include#include#include / fabs();using namespace std;void SolveEquation(const int &dim) double *Arr = new doubledim*dim; /方程矩陣 double *B = new doubledim; /方程右端系數(shù) double *X = new doubledim; /初始解 double *dX = new doubledim; double e; /解的精度 double w; /w=1時，為GS迭代 double S, S1; const int N = 5000; /最大迭代次數(shù),防止死機 int n=0; bool hasSolved = false; cout endl 初始化方程組左端系數(shù)矩陣! endl; for (int i=0; idim; i+) for (int j=0; jdim; j+) cout Arr i j Arri*dim+j; cout 初始化方程組右端系數(shù)矩陣! 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 設(shè)置解的精度要求:(比如1e-6) e; cout endl; cout 輸入初始向量! endl; for (i=0; idim; i+) cout X i+1 Xi; cout endl 設(shè)置參數(shù) ( 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()取絕對值函數(shù) S = S1; Xi = Xi + S1; n+; if (fabs(S) = e) cout迭代次數(shù)為：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 初始向量或 不合理, 請重新設(shè)置! endl; delete Arr; delete B; delete X; del