用高斯列主元消元法解线性方程组.doc

《用高斯列主元消元法解线性方程组.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、用C语言编写软件完成以下任务：请用高斯列主元消元法解下列线性方程组：方法说明（以4阶为例）：第1步消元——在增广矩阵（A，b）第一列中找到绝对值最大的元素，将其所在行与第一行交换，再对（A，b）做初等行变换使原方程组转化为如下形式：第2步消元——在增广矩阵（A，b）中的第二列中（从第二行开始）找到绝对值最大的元素，将其所在行与第二行交换，再对（A，b）做初等行变换使原方程组转化为：第3步消元——在增广矩阵（A，b）中的第三列中（从第三行开始）找到绝对值最大的元素，将其所在行与第二行交换，再对（A，b）做初等行变换使原方程组转化为：按x4àx3à

2、x2àx1的顺序回代求解出方程组的解总体流程图（一）定义数组x[NUMBER]，变量r，k，i，j，celect判断celect是否等于esc是否i=1i<=n输入系数和向量j=1j<=n+1将输入的数据存入A[i][j]中退j++i++k=1k<=n-1出调用mymax函数比较系数大小判断ark是否等于零是否判断f是否不等于k是否此方程组不合法调用change不执行函数i=k+1i<=nj=k+1j<=n+1j++k++x[n]=A[n][n+1]/A[n][n]k=n-1k>=1j=k+1j<=nme=me+A[k][j]*x[j]j++k

3、--输出结果调用message函数A[i][j]=A[i][j]–A[k][j]*A[i][k]/A[k][k]出j++k++x[u]=A[n][n+1]/A[n][n]求解最后一个方程组的解k=n-1k>=1j=k+1j<=np1=p1+A[k][j]*x[j]j++k--输出线性方程组的结果调用message函数