列主元素消去法.doc

2012-2013（1）专业课程实践论文列主元素消去法范 宁：0818180102，R数学08-1班夏之秋：0818180110，R数学08-1班一、算法理论列主元素消去法既是选主元高斯消去法的一种，也是实际计算中常用的部分选主元消去法。 列主元素消去法则是对完全主元素消去法的又一次改进。列主元素消去法在完全主元素消去法的基础上减少了在选主元素时所要花费的一定的计算时间。设有线性方程组其中，为非奇异矩阵。方程组的增广矩阵为 首先在的第1列选取绝对值最大的元素作为主元素，即选择然后交换的第1行与第行（交换后增广矩阵为简单起见仍记为，其元素仍记为）。经过第1次消元计算得到与原方程组等价的方程组其中 上述过程可记为 重复上述计算过程，现假设已完成第步的选主元素过程，交换两行并进行消元计此时约化为其中的元素仍记为，的元素仍记为.第步选主元素（在右下角方阵的第1列内选），即确定，使 交换第行与行的元素，再进行消元计算，最后将原线性方程组化为回代可求解得 二、算法框图输出迭代失败标志交换中两行输入（增广矩阵）开始输出奇异标志结束 三、算法程序#include #include#include #include #define max_dimension 20int n;static float amax_dimensionmax_dimension;static float bmax_dimension;static float xmax_dimension; void main()int i;int j;int d;int row;float temp;float known_items;float lmax_dimensionmax_dimension;system(cls);printf(Please Input Matrix jieshu :);scanf(%d,&n);printf(n);printf(Please Input Matrix Factors : );printf(n);for (i=0; in; i+) printf(input di %d hang dezhi:,i+1); for (j=0; jn; j+) scanf(%f,&aij); printf(n);printf(Please Input Changshu xiang: );for (i=0; in; i+) scanf(%f,&bi);printf(The Augmented(zenguang) Matrix is :nn);for (i=0; in; i+) for (j=0; jn; j+) printf(%f,aij); printf(%f,bi); printf(n);printf(n); for (d=0; dn-1;d+) row=d; for (i=d+1; ifabs(arowd) row=i; if (row!=d) for (j=d;jn; j+) temp=arowj; arowj=adj; adj=temp; temp=brow; brow=bd; bd=temp; for (i=d+1; in; i+) lid=-aid/add; for (j=d;jn; j+) aij=aij+adj*lid; bi=bi+bd*lid; printf(The shangsanjiaozenguang Matrix after predigestion is:nn);for (i=0; in; i+) for (j=0; j-1; i-) known_items=0; for (j=1; jn-i; j+) known_items=known_items+aii+j*xi+j; xi=(bi-known_items)/aii;printf(The Root X is :nn);for (i=0; in; i+) printf(%.5f ,xi);printf(nn);getch();四、算法实现例1. 求解方程组：用四位浮点数进行计算，精确解舍入到四位有效数字为解：运行程序(1)显示 Please Input Matrix jieshu：输入的值为3，回车。(2)显示input di 1 hang dezhi：输入，回车。(3)显示 input di 2 hang dezhi：输入，回车。(4)显示 input di 3 hang dezhi：输入，回车。(5)显示 Please Input Changshu xiang:：输入，回车。(6)显示结果： 例2. 用列主元素消去法解该方程解：运行程序(1)显示出 Please Input Matrix jieshu：输入的值为4，回车。(2)显示input di 1 hang dezhi：输入回车。(3)显示 input di 2 hang dezhi：输入回车。(4)显示 inp