浙江工业大学数值计算方法实验2.docxVIP

实验报告一、实验目的用下列方法解线性方程并比较计算结果精度(1)顺序消元法(2)列主元消元法二、实验方法（要求用自己的语言描述算法）顺序消元法本代码中是将题中矩阵转化为上三角矩阵，具体流程如下：通过双重for循环【外层for循环控制作为基准的行数，内层循环计算与基准行的倍数，并且进行行相减操作(行相减的for循环不考虑在内)】，最后通过回代求解出所有的解。列主元消元法本代码中在将矩阵转化为上三角矩阵的同时，事先需进行，预处理操作如下(化为上三角矩阵过程同顺序消元法)：在每次操作时先找出列中绝对值最大的元素，并与基准行先进行交换。三、实验代码#includeiostream#includeiomanip#includefstreamusingnamespace std;void display(double **a, double *b, int N);int main() {//ifstream cin(input.txt);//ofstream cout(output.txt);int N;//方程的阶double max, p;int i, j, k, l;cout 请输入方程的阶（输入0退出）：;while (cin N N) {double **a = newdouble*[N], *b = newdouble[N];double **c = newdouble*[N], *d = newdouble[N];double **m = newdouble*[N], *x = newdouble[N];for (i = 0; i N; i++) {a[i] = newdouble[N];c[i] = newdouble[N];m[i] = newdouble[N];}cout 请输入数组A： endl;for(i = 0; i N; i++)for (j = 0; j N; j++) {cin a[i][j];c[i][j] = a[i][j];}cout 请输入数组B： endl;for (i = 0; i N; i++) {cin b[i];d[i] = b[i];}//消元计算，将方程组约化为三角形方程组for (k = 0; k = N; k++)for (i = k + 1; i N; i++) {m[i][k] = a[i][k] / a[k][k];cout m[i][k] endl;for (j = k; j N; j++)a[i][j] -= m[i][k] * a[k][j];b[i] -= m[i][k] * b[k];display(a, b, N);}//回代求解double t = 0.0;x[N - 1] = b[N - 1] / a[N - 1][N - 1];for (i = N - 2; i = 0; i--) {x[i] = b[i];for (j = i + 1;j N;j++)x[i] -= a[i][j] * x[j];x[i] /= a[i][i];}//输出方程组的解cout 方程的解为： endl;for (i = 0; i N; i++)cout x[ i + 1 ]= x[i] endl;cout endl;//------------------------------分割线---------------------------------cout endl;for (k = 0; k N - 1; k++) {max = c[k][k];l = k;for (i = k + 1; i N; i++)//选取主元if (fabs(c[i][k]) fabs(max)) {max = c[i][k];l = i;}if (fabs(max) * 100000 1)cout det A = 0 endl;else {if (l != k) {//交换行p = d[l];d[l] = d[k];d[k] = p;for (j = 0;j N;j++) {p = c[l][j];c[l][j] = c[k][j];c[k][j] = p;}display(c, d, N);}for (i = k + 1; i N; i++) {//消元，化为三角方程组m[i][k] = c[i][k] / c[k][k];for (j = k;j N;j++)c[i][j] -= m[i][k] * c[k][j];d[i] -= m[i][k] * d[k];}display(c, d, N);}}if (fabs(c[N - 1][N - 1]) * 100000 1)cout det A=0 endl;else {//回代求解x[N -