C语言用高斯消元法求行列式

这篇具有很好参考价值的文章主要介绍了C语言用高斯消元法求行列式。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

目录

数学原理

选择主元

程序设计

整体流程与代码

测试函数

测试结果


数学原理

高斯消元法求行列式:利用初等行变换,化为上三角行列式,求其主对角线的乘积

行列式的初等行变换:

1)换行变换:交换两行(行列式需变号)

2)倍法变换:将行列式的某一行(列)的所有元素同乘以数k(行列式需乘K倍)

3)消法变换:把行列式的某一行(列)的所有元素乘以一个数k并加到另一行(列)的对应元素上(行列式不变)

上述三种变化中,本章将会用到换行变换消法变换。 

例如:

行列式A为:

C语言用高斯消元法求行列式

化为上三角行列式(选主元):

C语言用高斯消元法求行列式

A经过选主元与高斯消去后,化为上三角行列式(选主元见下文)

行列式是值为:

det(A)=1 * 1 * 2 * 6 = 12

选择主元

主元就是在矩阵消去过程中,每列的要保留的非零元素,用它可以把该列其他消去。在阶梯型矩阵中,主元就是每个非零行第一个非零元素就是主元。

高斯消去法在消元过程中可能出现零主元,即,这时消元过程将无法进行;也可能主元绝对值非常小,用它做除法将会导致舍入误差的扩散,使数值解不可靠。解决该问题的办法是避免使用绝对值过小的元素作主元。

选择主元的方法:

1)找到主对角线以下每列最大的数Max所在的行数k

2)利用初等行变换——换行变换,将k行与当前主元行互换(记录总共换行次数n)

3)以当前主元行为基,利用初等行变换——消法变换,将主对角线下方消0

4)行列式每次换行需变号,行列式最后的符号为

5)每次进行高斯消去前都必须选择主元,计算n维的行列式,则需要进行n-1次主元选择

程序设计

整体流程与代码

1)判断传入指针是否为空

2)判断矩阵维数,判断是否为方阵

3)为临时矩阵开辟空间

4)将原矩阵数据拷贝到临时矩阵中(保护原矩阵)

5)选择主元:利用初等行变换,找出每列绝对值最大的数,与主元行互换(1.提高一定的精度 2.避免原函数主对角线有0)

6)利用初等行变换进行消0

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<string.h>

double Det(double** src)
{
	//step 1
	//判断指针是否为空
	if (src == NULL)exit(-1);
	int i, j, k, row, col;
	double sum, tmp,** res;
	int count = 0,flag;
	//判断矩阵维数
	row = (double)_msize(src) / (double)sizeof(double*);
	col = (double)_msize(*src) / (double)sizeof(double);
	if (row != col)exit(-1);
	//step 2
	res = (double**)malloc(sizeof(double*) * row);
	for (i = 0; i < row; i++)
	{
		res[i] = (double*)malloc(sizeof(double) * row);
		memset(res[i], 0, sizeof(res[0][0]) * row);//初始化
	}
	//step 3
	//进行数据拷贝
	for (i = 0; i < row; i++)
	{
		memcpy(res[i], src[i], sizeof(res[0][0]) * row);
	}
	//step 4
	//找主元,绝对值最大的那一行,与主元行互换
	for (j = 0; j < col; j++)
	{
		flag = j;
		double Max = fabs(res[flag][j]);//用绝对值比较
		//默认当前主元行的数最大,从主对角线下方选择主元
		for (i = j; i < row; i++)
		{
			if (fabs(res[i][j]) > Max)
			{
				flag = i;
				Max = fabs(res[i][j]);
			}
		}
		if (i == j && i != flag)
		{
			count++;//记录互换次数
			for (k = 0; k < col; k++)
			{
				tmp = res[flag][k];
				res[flag][k] = res[i][k];
				res[i][k] = tmp;
			}
		}
		//将主对角下方元素消成0
		for (i = j + 1; i < row; i++)
		{
			double b = res[i][j] / res[j][j];
			for (k = 0; k < col; k++)
			{
				res[i][k] += b * res[j][k] * (-1);//初等行变换
			}
		}
	}
	//计算主对角线元素乘积
	sum = 1;
	for (i = 0; i < row; i++)
	{
		for (j = 0; j < col; j++)
		{
			if (i == j)
				sum *= res[i][j];
		}
	}
	//必须释放res内存!
	free(res);
	return pow(-1,count)*sum;
}

上述代码中:

  • malloc函数在动态内存规划一文中有详细讲解
  • 判断矩阵维数在C语言判断矩阵维数中有详细讲解
  • 行列式必须是方阵,因此row和col是相等的,代码中row对应行操作,col对应列操作
  • 因为高斯消元法是化为上三角行列式,所以每次消元时,起始的行数i=j+1,上三角部分不用消0
  • 最后只需要返回三角行列式主对角元素的乘积即可,在函数末尾需要释放内存

测试函数

为了方便测试,创建三个测试函数

创建矩阵函数:

double** MakeMat(int n)
{
	int i = 0;
	if (n <= 0)exit(-1);
	double** res = (double**)malloc(sizeof(double*) * n);
	if (res == NULL)exit(-1);
	for (i = 0; i < n; i++)
	{
		res[i] = (double*)malloc(sizeof(double) * n);
	}
	return res;
}

初始化函数:

void InitMat(double** src)
{
	if (src == NULL)exit(-1);
	int i, j, n;
	n = (double)_msize(src) / (double)sizeof(double*);
	for (i = 0; i < n; i++)
	{
		for (j = 0; j < n; j++)
		{
			src[i][j] = pow(i, j);
		}
	}
}

初始化为i的j次方 

打印函数:

void print(double** src)
{
	if (src == NULL)exit(-1);
	putchar('\n');
	int i, j, row, col;
	row = (double)_msize(src) / (double)sizeof(double*);
	col = (double)_msize(*src) / (double)sizeof(double);
	for (i = 0; i < row; i++)
	{
		for (j = 0; j < col; j++)
		{
			printf("%9.4lf", src[i][j]);
		}
		putchar('\n');
	}
}

测试结果

int main()
{
	int n = 4;
	double** arr = MakeMat(n);
	InitMat(arr);
	printf("原行列式:>");
	print(arr);
	printf("上三角行列式:>");
	double res = Det(arr);
	printf("计算结果:>");
	printf("%lf\n", res);
	return 0;
}

这里没有返回上三角行列式,只是在函数最后加了一个打印,对其进行观察

C语言用高斯消元法求行列式文章来源地址https://www.toymoban.com/news/detail-476459.html

到了这里,关于C语言用高斯消元法求行列式的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 线性代数——(期末突击)行列式(上)-行列式计算、行列式的性质

    目录 行列式 行列式计算 逆序数  行列式的性质 转置 两行(列)互换 两行(列)对应相等 提公因子 两行(列)对应成比例 某行(列)为零 行列式分裂 行列式变换及三角行列式 行列式:(i是行标,j是列标)  计算方法(以二阶行列式为例):主对角线(ad)减去次对角线

    2024年02月03日
    浏览(46)
  • 0105行列式按行(列)展开-行列式-线性代数

    在n阶行列式中,把 ( i , j ) 元 a i j (i,j)元a_{ij} ( i , j ) 元 a ij ​ 所在的第 i 行和第 j i行和第j i 行和第 j 列划去后,留下来的 n − 1 n-1 n − 1 阶行列式叫做 ( i , j ) 元 a i j (i,j)元a_{ij} ( i , j ) 元 a ij ​ 的余子式,记作 M i j M_{ij} M ij ​ ;记 A i j = ( − 1 ) i + j M i j , A i j 叫做 (

    2024年03月21日
    浏览(50)
  • 线性代数笔记1-二阶行列式和三阶行列式

    本笔记记录自B站《线性代数》高清教学视频 “惊叹号”系列 宋浩老师第一课 有2行2列,4个元素 ∣ a 11 a 12 a 21 a 22 ∣ begin{vmatrix} a_{11} a_{12}\\\\ a_{21} a_{22} end{vmatrix} ∣ ∣ ​ a 11 ​ a 21 ​ ​ a 12 ​ a 22 ​ ​ ∣ ∣ ​ a i j a_{ij} a ij ​ : i是行标,j是列标 ∣ a 11 a 12 a 21 a 22 ∣

    2023年04月09日
    浏览(45)
  • 【线性代数】P4 行列式相乘+范德蒙德行列式+克莱姆法则 cramer

    行列式相乘的原则,就是将第一个行列式中依次将每行的每个元素分别与第二个行列式每列的每个元素进行相加再相乘。 其实这样理解:已知两个行列式,如上,相乘有新行列式,新行列式左上角第一个值为: a 11 *b 11 +a 12 *b 21 +a 13 *b 31 实例2: 当然,三阶行列式无法与四阶

    2024年02月02日
    浏览(49)
  • 摆(行列式、杜教筛)

    有一个 n × n ntimes n n × n 的矩阵 A A A ,满足: A i , j = { 1 i = j 0 i ≠ j ∧ i ∣ j C otherwise A_{i,j}=begin{cases} 1 i=j\\\\ 0 inot=jland imid j\\\\ C text{otherwise} end{cases} A i , j ​ = ⎩ ⎨ ⎧ ​ 1 0 C ​ i = j i  = j ∧ i ∣ j otherwise ​ 求 det ⁡ ( A ) det(A) det ( A ) 。答案模 998244353 998244353 998244

    2024年02月19日
    浏览(37)
  • LA@行列式性质

    设行列式 ∣ A ∣ = d e t ( a i j ) |A|=mathrm{det}(a_{ij}) ∣ A ∣ = det ( a ij ​ ) ,行列式性质主要有5条 转置不变性质 行列式与它的转置行列式相等 或说经过转置,行列式的值不变(方阵 A A A 转置前后取行列式的值相等) ∣ A T ∣ = ∣ A ∣ |A^T|=|A| ∣ A T ∣ = ∣ A ∣ n n n 阶方阵 B = A T B=

    2024年02月14日
    浏览(45)
  • 线性代数——行列式

    一、行列式的性质 性质1 行列互换,其值不变,即 |A|=|A^{T}| 性质2 若行列式中某行(列)元素全为 0, 则行列式为 0 性质3 若行列式中某行(列)元素有公因子 k(kneq0) ,则 k 可提到行列式外面( 倍乘性质 ) $$ begin{vmatrix}a_{11}a_{12}cdotsa_{1n}\\\\vdotsvdotsvdots\\\\ka_{i1}ka_{i2}cdotska_{in}\\\\

    2024年04月26日
    浏览(35)
  • 讨论行列式的几何意义

    先给结论:行列式是线性变化的伸缩因子  给定一个二阶行列式: 由行列式的计算方法可得:  矩阵可表示为二维空间上的两个向量和,这两个向量可以组成一个平行四边形。 下面证明 (平行四边形) 由: 三角形的面积公式:                 得:  平行四边形的面积

    2024年02月01日
    浏览(73)
  • 矩阵——对称行列式快解

    1、先化成爪型行列式 2、再化成上三角或下三角 第一步:把第1行的1倍分别加至第2、3、4行,化为爪型行列式 第二步:把第2、3、4列的(-1)倍都加到第1列,化为上三角 第三步:得出结果

    2024年02月16日
    浏览(49)
  • 线性代数复习:行列式

    求行列式就是求这个行列式的值 二,三阶行列式:可以用:对角线法则和沙路法做 对角线法则: 主对角线和的值减去 副对角线积的和值。 a b c d : 值就是ad-bc 注意:n阶:n行n列. 1.下三角法则(主对角线以上都为0): 把行列式化为下三角行列式值等于主对角线的元素的值的

    2024年02月07日
    浏览(41)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包