#include <stdio.h>
#include <stdlib.h>
#define MAXN 100 //最大变量数量
int arr[MAXN][MAXN]; //保存增广矩阵
int result[MAXN]; //保存方程的解
int unuse_result[MAXN];//判断是否是不确定的变元
int unuse_num;
void swap(int *a,int *b) //交换两数
{
int t;
t=*a;
*a=*b;
*b=t;
}
int gcd(int a,int b) //返回最大公约数
{
int t;
while(b!= 0)
{
t=b;
b=a%b;
a=t;
}
return a;
}
int lcm(int a,int b) //返回最小公倍数
{
return a*b/gcd(a,b);
}
void debug(int equ,int var)
{
int i,j;
for(i=0;i<equ;i++)
{
for(j=0;j<var+1;j++)
printf("%d ",arr[i][j]);
printf("\n");
}
printf("\n");
}
int Gauss(int equ,int var)
{
int i,j,k,col;
int max_r,ta,tb,lcm1;
int temp,unuse_x_num,unuse_index;
col=0; //设当前处理列的值为0,表示从第1列开始处理
for(k=0;k<equ && col<var;k++,col++) //循环处理矩阵中的行
{
max_r=k; //绝对值最大行
for(i=k+1;i<equ;i++)
if(abs(arr[i][col])>abs(arr[max_r][col]))
max_r=i; //保存绝对值最大的行号
if(max_r!=k) //最大行不是当前行,则与第k行交换
for(j=k;j<var+1;j++)
swap(&arr[k][j], &arr[max_r][j]); //交换矩阵右上角数据
if(arr[k][col]==0) //说明col列第k行以下全是0了,则处理当前行的下一列
{
k--;
continue;
}
for(i=k+1;i<equ;i++) //查找要删除的行
{
if(arr[i][col]!=0) //左列不为0,进行消元运算
{
lcm1=lcm(abs(arr[i][col]),abs(arr[k][col])); //求最小公倍数
ta=lcm1/abs(arr[i][col]);
tb=lcm1/abs(arr[k][col]);
if(arr[i][col]*arr[k][col]<0) //相乘为负,表示两数符号不同
tb=-tb; //异号的情况是两个数相加
for(j=col;j<var+1;j++)
arr[i][j]=arr[i][j]*ta-arr[k][j]*tb;
}
}
}
for(i=k;i<equ;i++)//判断最后一行最后一列,若不为0,表示无解
if(arr[i][col]!=0)
return -1; //返回无解
if(k<var)//自由变元有var-k个,即不确定的变元至少有var-k个.
{
for(i=k-1;i>=0;i--)
{
unuse_x_num=0; //判断该行中不确定变量数量,若超过1个,则无法求解
for(j=0;j<var;j++)
{
if(arr[i][j]!=0 && unuse_result[j])
{
unuse_x_num++;
unuse_index=j;
}
}
if(unuse_x_num>1)
continue; // 无法求解出确定的解
temp=arr[i][var];
for(j=0;j<var;j++)
{
if(arr[i][j]!=0 && j!=unuse_index)
temp-=arr[i][j]*result[j];
}
result[unuse_index]=temp/arr[i][unuse_index]; // 求出该变元.
unuse_result[unuse_index]=0; //该变元是确定的
}
return var-k; //自由变元有var-k个
}
for(i=var-1;i>=0;i--) //回代求解
{
temp=arr[i][var];
for(j=i+1;j<var;j++)
{
if(arr[i][j]!=0)
temp-=arr[i][j]*result[j];
}
if(temp % arr[i][i]!=0) //若不能整除
return -2; //返回有浮点数解,但无整数解
result[i]=temp/arr[i][i];
}
return 0;
}
int main()
{
int i,j;
int equ, var;
printf("方程数:");
scanf("%d",&equ); //输入方程数量
printf("变量数:");
scanf("%d",&var); //输入变量数量
for(i=0;i<equ;i++) //循环输入各方程的系数
{
printf("第%d个方程的系数:",i+1);
for(j=0;j<var+1;j++) //循环输入一个方程的系数
{
scanf("%d", &arr[i][j]);
}
}
unuse_num=Gauss(equ,var); //调用高斯函数
if(unuse_num==-1) //无解
printf("无解!\n");
else if(unuse_num==-2) //只有浮点数解
printf("有浮点数解,无整数解!\n");
else if(unuse_num>0) //无穷多解
{
printf("无穷多解! 自由变量数量为%d\n",unuse_num);
for(i=0;i<var;i++)
{
if(unuse_result[i])
printf("x%d 是不确定的\n",i+1);
else
printf("x%d: %d\n",i+1,result[i]);
}
}
else
{
for(i=0;i<var;i++) //输出解
{
printf("x%d=%d\n",i+1,result[i]);
}
}
printf("\n");
return 0;
}
分享到:
相关推荐
数值方法中用高斯消元法求解线性方程组,c++源码。已调试成功。
在数值计算方法中用高斯消元法和列选主元来求线性方程组的解的C++源代码
高斯消去法是求解线性方程组的一种实用方法,绝大多数以matlab程序居多,虽然有部分C/C++程序,但是他们都有一些不足之处。本程序采用C/C++语言,实现高斯消去法,输入输出结果均可以分数和浮点数的形式展现,满足更...
#include #include #include using namespace std; double A[1010][1010]; double B[1010]; double l;... printf("Please enter the order of the matrix.... printf("Please enter the order of the matrix....
用c++实现了高斯消元法求解齐次方程,文件里有整个工程,解决方案在debug文件夹,用到的编译软件是VS2019
1. 用高斯消元法解方程组: 21.0x1+67.0x2+88.0x3+73.0x4 =141.0 76.0x1+63.0x2 + 7.0x3+20.0x4 =109.0 85.0x2+56.0x3+54.0x4 =218.0 19.3x1+43.0x2+30.2x3+29.4x4 =93.7
高斯消元法解线性方程 数值方法练习 是用C++编写的代码
包含完整实验报告和源码 要求: 使用高斯消元法求解n元一次线性方程组。
高斯消元法 求线性方程组 方程组采用随机数的方式进行赋值 注释较多 方便理解
消元法是指将多个方程式组成的方程组中的若干个变量通过有限次地变换,消去方程式中的变量,通过简化方程式,从而获取结果的一种解题方法。 消元法主要有代入消元法、加减消元法、整体消元法、换元消元法、构造消元...
文档内容为数值分析算法的C++实现。 算法包括:非线性方程求解、高斯消元发、高斯列主消元法、牛顿迭代法、割线法
利用高斯列主元消元法求解方程组的C++代码,用VC++6.0实现,通过更改输入参数可求一般线性方程组的解。
在VC++环境下打开高斯消元法解n阶线性方程组.dsw 可以用高斯消元求解线性方程组
单线程解线性方程组的代码,提供两种算法,一种方法为标准的,精度高但速度慢;另一种做了简化,精度稍低但速度快。另外还有一个用显卡(CUDA)解线性方程组的程序,在大维数(大于500维)的情况下速度更快。
列主元高斯消元法,可以求解线性方程组,比较好用
自选线性方程组,分别编制高斯列主元消元法和LU分解法的功能模块函数和完整程序实现求解。(以下方程组供选择) 2.自选线性方程组(系数矩阵为实对称正定矩阵),编制LDLT分解法的程序完成求解。(以下线性方程组...
在初值时,我们首次接触求解二元一次方程组,当时我只是知道求解的过程那样来写是理所当然的,就想知道这个方法是谁想出来的。...级,而高斯消元法只有 N 的平方级,远小于公式法,所以选择了高斯消元法。
采用一般的正规矩阵解法,由于数据量不是很大,所以采用动态内存的方法。同时也是联系一下二维数组的应用,即...先把数据读入一个矩阵A中,然后计算出ATA,然后读取向量b,计算AT*b,然后利用高斯消元法,得到解的值。
6.8.1 线性方程求解——高斯消元法 195 6.8.2 非线性方程求解——二分法 200 6.8.3 非线性方程求解——牛顿迭代法 202 6.9 小结 205 第7章 复杂的数值计算算法 206 7.1 拉格朗日插值 206 7.1.1 拉格朗日插值...
该代码设计是基于C++6.0平台,制作的一系列求解方程组的数值计算方法,几乎涵盖了数值计算中所有方程组的求解方法,相信会有很大的实用价值!