chengxu.zip_高斯消去法

＃i nclude <stdio.h> #define IMAX 3 #define JMAX 4 float matrix[IMAX][JMAX] = { {-1.0, 1.0, 2.0, 2.0 }, {3.0, -1.0, 1.0, 6.0 }, {-1.0, 3.0, 4.0, 4.0 } 　　　}; float solvec[IMAX] = { 0.0, 0.0, 0.0 }; main() { void wrt_output(void); void pivot(void); void triang(void); void back(void); void wrt_vector(void); (void)printf("The original matrix %d by %d :\n", IMAX, JMAX); (void)wrt_output(); (void)pivot(); (void)printf("The matrix after pivoting:\n"); (void)wrt_output(); (void)triang(); (void)printf("The matrix after lower decomposition:\n"); (void)wrt_output(); (void)back(); (void)printf("The solution vector after back substitution:\n"); (void)wrt_vector(); } /*********************************************************** * 决定参数矩阵中第一列的最大元素并移动第一列含有最大值的行 * * 到第一行. 然后重复对其他的行和列做这种处理, 对于每一次叠 * * 代, 列的位置和行的位置增加一(即, 第1行-第1列, 然后第2行- * * 第2列, 然后第3行-第3列, 等* ************************************************************/ void pivot() { int i, j, k; float maxval, tempval; for(j = 1; j < IMAX; j++) { 　maxval = matrix[j][j]; 　for ( i = (j + 1); i < IMAX; i++) { 　　if ( maxval < matrix[i][j] ) { 　　　maxval = matrix[i][j]; 　　　for( k = 0; k <= IMAX; k++) { 　　　　tempval = matrix[i][k]; 　　　　matrix[i][k] = matrix[j][k]; 　　　　matrix[j][k] = tempval; 　　　} 　　　} 　　} 　} } /*********************************** * 完成一个输入矩阵的低级分解的函数 * ************************************/ void triang(void) { int i, j, k; float pivot, pcelem; for ( j = 0; j < IMAX; j++) { pivot = matrix[j][j]; for ( k = ( j + 1 ); k <= IMAX; K++) { matrix[j][k] = matrix[j][k] / pivot; } for ( i = ( j + 1 ); i < IMAX; i++) { pcelem = matrix[i][j]; for ( k = ( j + 1 ); k <= IMAX; k++) { matrix[i][k] = matrix[i][k] - ( pcelem * matrix[j][k] ); } } } } /********************************************************* * 子过程从一个已经经历了低级分解的参数矩阵计算一个解向量 * *********************************************************/ void back(void) { int i, j; float sum; solvec[IMAX - 1] = matrix[IMAX - 1][JMAX -1]; for ( i = (IMAX -1); i > -1; i--) { 　　sum = 0.0; 　　for ( j = (i + 1); j < IMAX; j++) { 　　　　sum = sum + matrix[i][j] * solvec[j]; 　　} 　　solvec[i] = matrix[i][IMAX] - sum; } } void wrt_output(void) { 　int i, j; 　(void)printf("**************************************\n"); for ( i = 0; i < IMAX; i++) { 　　for ( j = 0; j < (JMAX - 1); j++) { 　　　　(void)printf("%f", matrix[i][j]); 　　} 　　(void)printf("%f\n", matrix[i][JMAX - 1]); } (void)printf("****************************************\n"); } void wrt_vector(void); { (void)printf("*************************************\n"); 　　(void)printf("%f", solvec[0]); 　　(void)printf(" %f", solvec[1]); 　　(void)printg(" %f\n", solvec[2]); (void)printf("****************************************\n"); } /*********************************************************** * 测试函数bisec.f的程序, bisec.f 决定一个方程(f中)的解 　　* * 可是这个函数确实假设函数-f由两个值支撑. 即在用户给定的终 * * 点之间的解不超过一个　　　　　　　　　　　　　　　　　　 * ************************************************************/ ＃i nclude <stdio.h> ＃i nclude <math.h> main() { 　void bisec(float init_left_val, float init_right_val); 　float f(float value): 　　　　float xleft, xright; 　　　　char line[100]; 　　　　(void)printf("Please enter an initial left and right value:"); 　　　　(void)fgets(line, sizeof(line), stdin); 　　　　(void)sscanf(line, "%f %f", &xleft, &xright ); 　　　　(void)bisec(xleft, xright); 　　　　return(0); }