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)


相关推荐

  • 根据两点的经纬度计算距离_经纬度两点距离

    根据两点的经纬度计算距离_经纬度两点距离问题提出目前手头的一个项目要用到GPS地理定位信息,很自然的就需要知道两个地点之间的距离,于是上网找了一下。背景知识这些经纬线是怎样定出来的呢?地球是在不停地绕地轴旋转(地轴是一根通过地球南北两极和地球中心的假想线),在地球中腰画一个与地轴垂直的大圆圈,使圈上的每一点都和南北两极的距离相等,这个圆圈就叫作“赤道”。在赤道的南北两边,画出许多和赤道平行的圆圈,就是“纬圈”;构成这些圆圈的线段,叫做纬线

  • IIS中WEB服务器的日志存放到SQL Server 2005中

    IIS中WEB服务器的日志存放到SQL Server 2005中

  • codetemp-2: CObject::AssertValid()

    codetemp-2: CObject::AssertValid()原型:[virtual]voidCObject::AssertValid()const//说明:该成员函数只能对成员变量做”引用型”操作,不能做”加工型”操作{ASSERT(this!=NULL);}MSDN:Whenyouwriteyourownclass,youshouldoverridetheAssertValidfunc…

  • 2020阿里笔试编程题[通俗易懂]

    2020阿里笔试编程题[通俗易懂]选择题很难做,阿里的秋招貌似非常难,大部分岗位都留给了实习生,但是两道编程题不怎么难。第一题有一个n*n的地图,一只兔子想要穿过这个地图,给出的地图是一个二维数组map[i][j],数值表示该位置的毒雾持续时间,当兔子在(x,y)位置时,它可以跳到(x+2,y)或者(x,y+2)位置,跳的时候需要对应等待map[x+1][y]或者map[x][y+1]的时间,兔子开始跳的位置从map[1][1…

  • tomcat最大并发数连接数_lvs最大并发数

    tomcat最大并发数连接数_lvs最大并发数Tomcat默认配置的最大请求数是150,也就是说同时支持150个并发如何设置提高并发数修改catalina.sh:rem以下配置为JVM参数调优setJAVA_OPTS=-serverrem以服务器模式启动,启动速度慢,但更稳定,性能更好-Xms8192Mrem由于本机内存为16G,这里就设置成8G(实际并未达到最大内存的80%)-Xm…

    2022年10月31日
  • 关于sstream的灵活使用

    关于sstream的灵活使用问题有10000个队伍参加。经过工作人员认真负责的统计,本来已经统计好了这一万个队伍的分数和排名,并按照排名从高到低依次进行了编号(从1到10000)但是由于一个非常偶然的因素,导致其中三个编号的数据丢失,而且剩余编号的顺序也全被打乱了。你需要编写一个程序,根据还保留的统计数据,来判断哪些编号的数据丢失了,并将这些编号按照从小到大的顺序重新拼接为一个新数字,然后计算这个新数字除以11的余数。如…

发表回复

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

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