求行列式的值(递归)

本文介绍了一个用于计算矩阵行列式的C语言程序实现。该程序通过递归的方法计算任意大小的方阵行列式,并提供了输入矩阵、计算子矩阵及输出结果等功能。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

/*
#include <math.h>
#include <stdio.h>
#define N 100
#define LIM -100000000
float det(float a[N][N],int n){
    if(n==1)
        return a[0][0];
    if(n==2)
        return a[0][0]*a[1][1]-a[0][1]*a[1][0];// the base situation
    else{
        int j,i,flag=1;
        float ret=0;
        for(j=0;j<n;j++){//j is the column
            int st_c=0,st_r=0;
            float a0j=a[0][j];
            flag=-flag;
            a[0][j]=LIM;
            for(i=1;i<n;i++){//compose the trans_a
                a[i][j]=LIM;
            }//end reseting a
            float trans_a[N][N];
            for(i=1;i<n;i++){
                for(j=0;j<n;j++){
                    if(a[i][j]!=LIM){
                        if(st_c==n-1){
                            st_r++;
                            st_c=0;
                        }
                        trans_a[st_r][st_c]=a[i][j];
                        st_c++;
                    }
                }//end composing
            }
            ret+=det(trans_a,n-1)*a0j*flag;
            printf("%f\n",ret);
        }

    return det(a,n-1);
    }

}
int main(){
    printf("Please enter matrix size n(1<=n<20):");
    int n;
    scanf("%d",&n);
    printf("Please input matrix line by line:\n");
    float a[N][N];
    int i,j;
    for(i=0;i<n;i++)
        for(j=0;j<n;j++)
            scanf("%f",&a[i][j]);
    printf("matrix a:\n");
    for(i=0;i<n;i++){
        for(j=0;j<n;j++)
            printf("%6.1f",a[i][j]);//the elements
        printf("\n");
    }
    printf("result = %f\n",det(a,n));//the res
    return 0;
}*/

#include <math.h>
#include <stdio.h>
 
#define  CONST 1e-6
#define  SIZE 20
 
void InputMatrix (double a[][SIZE], int n);
double DeterminantValue(double a[][SIZE], int n);
void SubMatrix(double a[][SIZE], double b[][SIZE], int n, int row, int col);
void PrintMatrix(double a[][SIZE], int n);
int main(void)
{
    double  a[SIZE][SIZE];
    int     n;
    double  result;
    printf("Please enter matrix size n(1<=n<%d):", SIZE);
    scanf("%d", &n);
    printf("Please input matrix line by line:\n");
    InputMatrix(a, n);
    printf("matrix a:\n");
    PrintMatrix(a, n);
    printf("\n");
    result = DeterminantValue(a, n);
    printf("result = %f\n", result);
    return 0;
}
//  函数功能: 输入一个n×n矩阵的元素
void InputMatrix (double a[][SIZE], int n)
{
    int i, j;
    for (i = 0; i < n; i++)
    {
        for (j = 0; j < n; j++)
        {
            scanf("%lf", &a[i][j]);
        }
    }
}
//  函数功能: 计算n×n矩阵的行列式的值
double DeterminantValue(double a[][SIZE], int n)
{
    int    i = 0, j = 0;
    double temp, result, b[SIZE][SIZE];
    if (n == 1)
    {
        result = a[0][0];
    }
    else if (n == 2)
    {
        result = a[0][0] * a[1][1] - a[0][1] * a[1][0];
    }
    else
    {
        result = 0.0;
        for (j = 0; j < n; j++)//遍历列
        {
            SubMatrix(a, b, n, i, j);//the key
            printf("Submatrix:\n");
            PrintMatrix(b, n - 1);
            temp = DeterminantValue(b, n - 1);        //递归
            result += pow(-1, i + j) * a[0][j] * temp;
            printf("DValue of the Submatrix is %6.1f\n", temp);
        }
    }
    return result;
}
//  函数功能: 计算n×n矩阵a中第row行col列元素的(n-1)×(n-1)子矩阵b
void SubMatrix(double a[][SIZE], double b[][SIZE], int n, int row,
               int col)
{
    int i, j, ii = 0, jj = 0;
    for (i = 0; i < n; i++)
    {
        jj = 0;//回归0
        for (j = 0; j < n; j++)//复制第j列的数
        {
            if (i != row && j != col)//过滤行列为i j的项
            {
                b[ii][jj] = a[i][j];
                jj++;
            }
        }
        if (i != row && j != col)//复制完后向下推进一行
        {
            ii++;////////////////////////////////////////////////////////
        }        ///////两次的i!=row&&...直接不考虑row,col上的数/////////
    }            ////////////////////////////////////////////////////////
}
//  函数功能: 输出一个n×n矩阵的元素
void PrintMatrix(double a[][SIZE], int n)
{
    int i, j;
    for (i = 0; i < n; i++)
    {
        for (j = 0; j < n; j++)
        {
            printf("%6.1f\t", a[i][j]);
        }
        printf("\n");
    }
}

 

转载于:https://www.cnblogs.com/xzenith/p/3702688.html

资源下载链接为: https://pan.quark.cn/s/1bfadf00ae14 “STC单片机电压测量”是一个以STC系列单片机为基础的电压检测应用案例,它涵盖了硬件电路设计、软件编程以及数据处理等核心知识点。STC单片机凭借其低功耗、高性价比和丰富的I/O接口,在电子工程领域得到了广泛应用。 STC是Specialized Technology Corporation的缩写,该公司的单片机基于8051内核,具备内部振荡器、高速运算能力、ISP(在系统编程)和IAP(在应用编程)功能,非常适合用于各种嵌入式控制系统。 在源代码方面,“浅雪”风格的代码通常简洁易懂,非常适合初学者学习。其中,“main.c”文件是程序的入口,包含了电压测量的核心逻辑;“STARTUP.A51”是启动代码,负责初始化单片机的硬件环境;“电压测量_uvopt.bak”和“电压测量_uvproj.bak”可能是Keil编译器的配置文件备份,用于设置编译选项和项目配置。 对于3S锂电池电压测量,3S锂电池由三节锂离子电池串联而成,标称电压为11.1V。测量时需要考虑电池的串联特性,通过分压电路将高电压转换为单片机可接受的范围,并实时监控,防止过充或过放,以确保电池的安全和寿命。 在电压测量电路设计中,“电压测量.lnp”文件可能包含电路布局信息,而“.hex”文件是编译后的机器码,用于烧录到单片机中。电路中通常会使用ADC(模拟数字转换器)将模拟电压信号转换为数字信号供单片机处理。 在软件编程方面,“StringData.h”文件可能包含程序中使用的字符串常量和数据结构定义。处理电压数据时,可能涉及浮点数运算,需要了解STC单片机对浮点数的支持情况,以及如何高效地存储和显示电压。 用户界面方面,“电压测量.uvgui.kidd”可能是用户界面的配置文件,用于显示测量结果。在嵌入式系统中,用
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值