C语言 n*n矩阵求值及求逆矩阵[通俗易懂]

求出矩阵的值以及输出逆矩阵,英语不好,略拗口。上代码:#include#include#includeintmap1[110][110],i,j,k,l,m,n,p,vis[110][110],num;intjust[110][110];intjust2[110][110];voidprint(intmap[][110],intn)//打印矩阵{ intj;

大家好,又见面了,我是你们的朋友全栈君。

C语言 n*n矩阵求值及求逆矩阵[通俗易懂]

求出矩阵的值以及输出逆矩阵,英语不好,略拗口。

上代码:

#include<stdio.h>
#include<string.h>
#include<math.h>
int InitialMatrix[110][110];//初始矩阵,即输入的矩阵
int CurrentMatrix[110][110];//当前 矩阵 多用来表示当前余子式

//打印矩阵matrix
void print(int matrix[][110], int n)//打印矩阵 
{
    int j;
    for (int i = 1; i <= n; i++)
    {
        for (j = 1; j < n; j++)
            printf("%d ", matrix[i][j]);
        printf("%d\n", matrix[i][j]);
    }
}

//得到矩阵matrix第numi行第numj列的余子式 并存入到CurrentMatrix中
void GetCurrentMatrix(int numi, int numj, int matrix[][110], int MatrixSize)
{
    //tempMatrix为临时矩阵
    int tempMatrix[110][110];
    //先全部置为1
    for (int i = 0; i < MatrixSize + 10; i++)
    {
        for (int j = 0; j < MatrixSize + 10; j++)
            tempMatrix[i][j] = 1;
    }

    for (int i = 1; i <= MatrixSize - 1; i++)//求余子式矩阵 
    {
        for (int j = 1; j <= MatrixSize - 1; j++)
        {
            if (i >= numi&&j < numj)
                tempMatrix[i][j] = matrix[i + 1][j];
            else if (i >= numi&&j >= numj)
                tempMatrix[i][j] = matrix[i + 1][j + 1];
            else if (i < numi&&j >= numj)
                tempMatrix[i][j] = matrix[i][j + 1];
            else if (i < numi&&j < numj)
                tempMatrix[i][j] = matrix[i][j];
        }
    }
        
    //临时矩阵tempMatrix赋值给 当前矩阵CurrentMatrix
    for (int i = 0; i <= MatrixSize - 1; i++)
    {
        for (int j = 0; j <= MatrixSize - 1; j++)
            CurrentMatrix[i][j] = tempMatrix[i][j];
    }
}

//求当前行列式为CurrentMatrix时,CurrentMatrix的值。
int GetMatrixValue(int MatrixSize)
{
    int matrix[110][110];
    //赋值给matrix  交给matrix去计算 防止破坏CurrentMatrix
    for (int i = 1; i <= MatrixSize; i++)
    {
        for (int j = 1; j <= MatrixSize; j++)
            matrix[i][j] = CurrentMatrix[i][j];
    }

    //一阶二阶 行列式求值
    if (MatrixSize == 2)
        return matrix[1][1] * matrix[2][2] - matrix[1][2] * matrix[2][1];
    else if (MatrixSize == 1)
        return matrix[1][1];
    else//高阶行列式求值
    {
        int Matrixvalue = 0;//CurrentMatrix的值
        for (int i = 1; i <= MatrixSize; i++)
        {
            GetCurrentMatrix(1, i, matrix, MatrixSize);//得到第一行 第i列 的 当前余子式 CurrentMatrix 

            if (i % 2 == 1)//因为 是 行项为1 所以 列项i为奇数时  逆序数为偶数
                Matrixvalue += matrix[1][i] * GetMatrixValue(MatrixSize - 1);
            else
                Matrixvalue -= matrix[1][i] * GetMatrixValue(MatrixSize - 1);
        }
        return Matrixvalue;
    }
}
int gcd(int m, int n)//求逆矩阵时约分 
{
    if (m < n)
        gcd(n, m);
    if (n == 0)
        return m;
    else
        return gcd(n, m%n);
}

//打印当前两个值相除得到的最简分数
void final(int n, int m) 
{
    if (n*m < 0)
    {
        printf("-");
        final(fabs(n), fabs(m));
        return;
    }
    //printf("%d %d /%d\n",n,m,gcd(n,m));
    if (m == 1)
        printf("%d\t", n);
    else if (n%m == 0)
        printf("%d\t", n / m);
    else
        printf("%d/%d\t", n / gcd(n, m), m / gcd(n, m));
}
int main()
{
    while (1)
    {
        int MatrixSize = 0;//矩阵的size
        printf("please input the determinant`s size:\n");//输入矩阵的规格 
        scanf("%d", &MatrixSize);
        //初始数组 全部 置为 1
        for (int i = 0; i < MatrixSize + 10; i++)
        {
            for (int j = 0; j < MatrixSize + 10; j++)
                InitialMatrix[i][j] = 1;
        }
            
        printf("please input the determinant:\n");//输入矩阵 
        for (int i = 1; i <= MatrixSize; i++)
        {
            for (int j = 1; j <= MatrixSize; j++)
            {
                scanf("%d", &InitialMatrix[i][j]);
                CurrentMatrix[i][j] = InitialMatrix[i][j];//CurrentMatrix即 为当前矩阵 
            }
        }


        printf("the value of the determinant is:\n");
        int MatrixValue = GetMatrixValue(MatrixSize);//矩阵的值 


        printf("ans is %d\n", MatrixValue);
        int TransposeMatrix[110][110];//转置行列式

        for (int i = 1; i <= MatrixSize; i++)//求转置行列式 
        {
            for (int j = 1; j <= MatrixSize; j++)
            {
                TransposeMatrix[i][j] = InitialMatrix[j][i];
                CurrentMatrix[i][j] = TransposeMatrix[i][j];
            }
        }
            
        printf("Transpose the determinant is:\n"); //转置矩阵 
        printf("the value of the determinant is:\n"); //转置矩阵的值 
        print(TransposeMatrix, MatrixSize);

        int TransposeMatrixValue = GetMatrixValue(MatrixSize);
        printf("%d\n", TransposeMatrixValue);
        if (TransposeMatrixValue == 0)
        {
            printf("this matrix dont have inverse matrix!!!\n");
            continue;//矩阵值为0,无逆矩阵 
        }
        printf("***************\n");
        printf("inverse matrix is:\n");//输出逆矩阵 
        for (int i = 1; i <= MatrixSize; i++)
        {
            for (int j = 1; j <= MatrixSize; j++)
            {
                GetCurrentMatrix(i, j, TransposeMatrix, MatrixSize);//得到转置矩阵第i行第j列的 余子式 存入到CurrentMatrix里边
                if ((i + j) % 2 == 0)
                    final(GetMatrixValue(MatrixSize - 1), MatrixValue);
                else
                    final(-1 * GetMatrixValue(MatrixSize - 1), MatrixValue);
            }
            printf("\n");
        }
    }
}

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/129010.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • java执行curl命令带json(执行shell脚本命令)

    需要注意的是,命令需要以数组的形式传递参数,就是把正常的命令以空格切分成数组就行了。而且参数前后不能有空格,不然会报错。具体的测试代码如下:importjava.io.BufferedReader;importjava.io.IOException;importjava.io.InputStreamReader;publicclassTestCurl{p…

  • 在鱼眼和全向视图图像的深度学习方法

    点击上方“3D视觉工坊”,选择“星标”干货第一时间送达作者:黄浴https://zhuanlan.zhihu.com/p/88675419本文仅做学术分享,如有侵权,请联系删除。‍‍普通…

  • 搭建SpringCloud项目,并实现自动化部署[通俗易懂]

    搭建SpringCloud项目,并实现自动化部署[通俗易懂]前一阵子一直在忙着做项目,最近有空,把项目中用到的SpringBoot微服务相关架构再梳理一遍。项目资源:https://github.com/sunroyi/SpringCloud主要分为以下几步:(1)搭建SpringBootService,这里是各个微服务的业务逻辑。(这里搭建了2个Service,用来测试熔断)(2)搭建SpringBootEureka,用来发现服务。(…

  • JavaScript[容易忽视的错误]:当续行遇到换行,换行符丢失

    JavaScript[容易忽视的错误]:当续行遇到换行,换行符丢失1、引言在JavaScript中,当定义一个大的字符串,特别是有换行情况时,为了看起来整齐、易读,一般使用续行符,例如: varscript=”varchart=anychart.pieChart([\[‘Chocolatepaste’,5],\[‘Whitehoney’,2],\

  • cocos2dx 3.0 研究(4)渲染分析[通俗易懂]

    cocos2dx 3.0 研究(4)渲染分析

  • 学生网上选课管理系统_选课管理系统

    学生网上选课管理系统_选课管理系统**数据库系统原理课程设计报告**学生选课管理系统(上)设计内容与要求:1、系统用户由三类组成:教师、学生和管理员。2、管理员负责的主要功能:①用户管理(老师、学生及管理员的增、删、改);②课程管理(添加、删除和修改);③选课管理(实现选课功能开放和禁止、老师成绩输入开放和禁止)。3、学生通过登录,可以查询课程的基本信息、实现选课、退课和成绩查询;4、老师通过登录,可以查看选…

    2022年10月16日

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号