C语言实现高斯消去法和列主元高斯消去法

本篇主要实现高斯消去法和列主元高斯消去法

高斯消去法和列主元高斯消去法都是为了解线性方程组的有效方法，但列主元高斯消去法是高斯消去法的一个优化版本，强烈建议后面许多地方用到解方程组时，都用列主元高斯消去法。

高斯消去法：

我个人觉得，例子比数学公式更好让人理解本质。上述是线性代数的高斯消元法，而我们只需要用程序将这个代码实现出来。上述方程组，第二行减去第一行，第三行减去两倍第一行。那么可以推断，若第三行要消去第一个数，则让第一行每个数除上第一行第一个数并乘上第三行第一个数，再让第三行减去第一行，即可。

代码实现：

#include<stdio.h>
# define N 100
int main()
{int k,n,i,j;double a[N][N],b[N],m[N][N],x[N],c[N],s=0.0;printf("请输入未知数的个数："); //该方程组的未知数scanf("%d",&n);printf("请输入数据："); //输入增广矩阵for(i=0;i<n;i++){for(j=0;j<n+1;j++){scanf("%lf",&a[i][j]);}}for(i=0;i<n;i++){b[i]=a[i][n];//将增广矩阵最后一列作为y的值}for(k=0;k<n-1;k++){for(i=k+1;i<n;i++){   m[i][k]=1.0*a[i][k]/a[k][k];//这一步减少了计算机运算量，因为后面要用到b[i]=b[i]-1.0*m[i][k]*b[k];//将y与每一行的变化保持一致for(j=k;j<n;j++){a[i][j]=a[i][j]-1.0*m[i][k]*a[k][j];//这一步实现消元，可自己用具体数值推导}}}//可以在这里加两个循环，查看矩阵中每一个位置的数是多少。就是上面输入增广矩阵下面那个双循环//方程回代过程 x[n-1]=1.0*b[n-1]/a[n-1][n-1];//先计算最后未知数，从下上算for(i=n-2;i>=0;i--){    s=0;                   //这一步非常重要，每一步都要更新为0，不然s保留参与下面计算出错for(j=i+1;j<n;j++){c[j]=1.0*a[i][j]*x[j];//计算过程s=c[j]+s;}x[i]=1.0*(b[i]-s)/a[i][i];//将每一行的方程除未知数其他全部放到等式右边}for(i=0;i<n;i++){printf("x[%d]=%lf\n",i+1,x[i]); //将未知数按顺序打印出来}return 0;
}

列主元高斯消元法，只是将每一列的最大数先放到前面，因为消元的过程中，如第三行第一个要消去第一个数字，则要第一行每一个数除上第一个数再乘上第三行第一个数，而第一行第一个数太小则会引起误差，若为0，则不能做分母，因此有必要引入列主元高斯消元法。

列主元高斯消元法：

代码实现：

#include<stdio.h>
# define N 100
int main()
{//我的坏毛病，参数过多，后面不一定有用上int k,n,i,j,r,w,e,g;double a[N][N+1],b[N],m[N][N],x[N],c[N],s=0,f,max,uu[N];printf("请输入未知数的个数：");scanf("%d",&n);printf("请输入数据：");for(i=0;i<n;i++){for(j=0;j<n+1;j++){scanf("%lf",&a[i][j]);}}for(i=0;i<n;i++){b[i]=a[i][n];}/*验证代码是否正确*//可将这段代码放到任何一个位置进行正确性检验for(i=0;i<n;i++){   for(j=0;j<n;j++){printf("%lf ",a[i][j]);}printf("%lf\n",b[i]);}*/ for(r=0;r<n-1;r++)//这一步开始寻找每一列最大的数{   k=r;max=a[r][r];for(w=r;w<n-1;w++)//寻找，若找到该列下面有比起始大的数，则进行标记{if(a[w+1][r]>=a[w][r]){max=a[w+1][r];k=w+1;          //若找到下列大的数，标记}}if(k>=r)//此时表示找到{ for(e=r;e<n;e++)//将下列最大的数的那一行与这列起始数的那一行所有数进行行更换{f=a[r][e];a[r][e]=a[k][e];a[k][e]=f;}f=b[r];//让y保持与行同变化b[r]=b[k];b[k]=f;    }   for(i=r+1;i<n;i++)//与上面相同，进行高斯消去法，因为列已经变化完毕{   m[i][r]=1.0*a[i][r]/a[r][r];b[i]=b[i]-1.0*m[i][r]*b[r];for(j=r;j<n;j++){a[i][j]=a[i][j]-1.0*m[i][r]*a[r][j];}}}//回代过程x[n-1]=1.0*b[n-1]/a[n-1][n-1];for(i=n-2;i>=0;i--){    s=0;for(j=i+1;j<n;j++){c[j]=1.0*a[i][j]*x[j];s=c[j]+s;}x[i]=1.0*(b[i]-s)/a[i][i]; }for(i=0;i<n;i++){printf("x[%d]=%lf\n",i+1,x[i]); }return 0;
}

C语言实现高斯消去法和列主元高斯消去法相关推荐

【计算方法】#01 高斯消去法和列主元高斯消去法的原理简介及C++实现
[计算方法]#01 高斯消去法和列主元高斯消去法的原理简介及C++实现 1. 高斯消去法 1.1 算法的适用条件 1.2 算法步骤和公式 1.3 算法复杂度分析 1.4 算法的C++实现 2 列主元高 ...
【数值分析】顺序高斯消去法和列主元高斯消去法的三个主要不同点
概要求解线性方程组 A x = b Ax=b Ax=b 可以使用[顺序高斯消去法]和[列主元高斯消去法],本文试列举二者的三个主要不同点. 不同点 1. 使用条件
列主元高斯消去法(c语言)(可以实现所有阶的)(超级详细)
其实列主元高斯消去法无非就是比之前的高斯消去法多了一个判断主元这个步骤,但是里面还是有一些小细节的,比如:你要求一个3*4的增广矩阵,你的主元只需要选两次,第一次是在第一列的0.1.2里面选,第二次就 ...
列主元高斯消去法解线性方程组——C语言实现
原理高斯消去法的基本原理就是用初等变换将用行的,逐次消去未知数的方法,把原来的方程组,化为与其等价的上三角方程组. 设有线性方程Ax=B\boldsymbol{A}\boldsymbol{x}=\ ...
列主元高斯消去法数学原理及超级纯手工Python实现
一.引言从高斯消去法,我们看到还是有缺陷,高斯消去法中,当对角元素=0时,消去无法进行,当对角元素很小的时候,导致其他它元素数量级严重增长和舍入误差扩散1,使结果不可靠.因此引出了列主元高斯消去法. ...
列主元高斯消去法Matlab实现
列主元高斯消去法解线性方程组在求解线性方程组中,可以使用列主元高斯消去法. 步骤如下: 现有方程组:Ax = b (1)提取出方程组的增广矩阵 A = [A b]: (2)进行n-1次消元,分别对应 ...
列主元高斯消去法的C++实现
下述所有内容都是建立在线性方程组有唯一解的情况高斯消去法主要用来求解线性方程组比如求解下图中的四维线性方程组该方程组写成行列式形式如下图所示首先将行列式变为阶梯行列式,以第一行为例: 第一行同时 ...
关于列主元高斯消去法的matlab实现
function x = nagauss1(a,b ,flag) %用途:列主元高斯消去法解线性方程组ax=b % A:系数矩阵,b:右端列变量 % flag:若为0,则显示中间过程,否则不显示,默 ...
Matlab 列主元高斯消去法
列主元高斯消去法是高斯消去法的优化版本,通过找出列中的最大值,避免了除以很小的数时容易引起的数值偏差. function x=Elim_Gauss(Matrix,n,b) %列主元高斯消去法 %输入- ...

C语言实现高斯消去法和列主元高斯消去法

C语言实现高斯消去法和列主元高斯消去法相关推荐

最新文章

热门文章