C语言程序设计整理(3)求n阶行列式的值(n<100)

【内容概述】

   本文主要讲述行列式应该怎么样用C语言来求值。

   行列式可以按行展开,利用代数余子式来求值,其中运用到递归的算法

【代码呈现】

/*利用递归的算法来求行列式的值*/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
void readin ( int n , int det[100][100] );
int calculate ( int n , int det[100][100] );
int minor ( int det[100][100] , int i , int n );
int main (void)
{
    int n , det[100][100] = {0};
    printf ("请输入行列式的阶数:(少于100)");
    scanf ("%d", &n);
    readin ( n , det );
    printf("行列式的值为 %d",calculate ( n , det ));
    system ( "pause" );
    return 0;
}
/*输入行列式的每一个元素*/
void readin ( int n , int det[100][100] )
{
    printf ("请输入行列式每行每列的元素,同行元素用空格隔开");
    int i , j;
    for ( i = 0; i < n; i++ )
    {
        for ( j = 0; j < n; j++ )
        {
            scanf("%d", &det[i][j]);
        }
    }
    return;
}
/*用递归的算法来求出行列式的值*/
int calculate ( int n , int det[100][100] )
{
    int i , m , sum = 0;
    if ( n == 1 )
    {
        return det[0][0];
    }
    else 
    {
        for ( i = 0; i < n; i++ )      //按第一行展开
        {
            m = minor ( det , i , n );
            sum += pow (-1 , i+2) * m * det[0][i];
        }
        return sum;
    }
}
int minor ( int det[100][100] , int i , int n )
{
	int  j , k;
	int det2[100][100];
    for (j = 0; j < n - 1; j++)        //构造出余子式
	{
		for (k = 0; k < n - 1; k++)
		{
			if (k < i)
            	det2[j][k] = det[j + 1][k];
			else if (k >= i)
				det2[j][k] = det[j + 1][k + 1];
		}
	}
	return calculate ( n-1 , det2 );   //回到calculate函数

  • 2
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
可以使用递归的方法来解n方阵的行列式。具体步骤如下: 1. 如果矩阵的数为1,则行列式就是该矩阵中唯一元素的。 2. 如果矩阵的数大于1,则可以选取第一行或第一列作为展开元素,假设选取了第一行。然后对于每个元素a[1][j],构造一个(n-1)子矩阵A',其中A'是将第一行和第j列删除后的子矩阵。 3. 对于每个元素a[1][j],计算其代数余子式A[1][j],即A[1][j] = (-1)^(1+j) * det(A'),其中det(A')表示A'的行列式。 4. 根据展开元素所在的行列的奇偶性,将所有代数余子式的相加,得到矩阵的行列式。 下面是用C语言实现解矩阵行列式代码: ``` #include <stdio.h> #define MAXN 100 // 矩阵的行列式 double det(double a[][MAXN], int n) { double ans = 0; if (n == 1) { ans = a[0][0]; } else { double temp[MAXN][MAXN]; int sign = 1; for (int k = 0; k < n; k++) { int i, j; for (i = 1; i < n; i++) { for (j = 0; j < k; j++) { temp[i - 1][j] = a[i][j]; } for (j = k + 1; j < n; j++) { temp[i - 1][j - 1] = a[i][j]; } } double subdet = det(temp, n - 1); ans += sign * a[0][k] * subdet; sign = -sign; } } return ans; } int main() { double a[MAXN][MAXN]; int n; printf("请输入矩阵的数n:"); scanf("%d", &n); printf("请输入矩阵的元素:\n"); for (int i = 0; i < n; i++) { for (int j = 0; j < n; j++) { scanf("%lf", &a[i][j]); } } printf("矩阵的行列式为:%lf\n", det(a, n)); return 0; } ``` 这里使用了一个名为det的函数来计算矩阵的行列式,该函数使用了递归的方法实现。在主函数中,用户需要输入矩阵的数和元素,程序会输出矩阵的行列式

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值