求逆矩阵以及两矩阵相乘的算法实现

最新推荐文章于 2025-07-14 15:34:56 发布

pathuang68

最新推荐文章于 2025-07-14 15:34:56 发布

阅读量1.1w

点赞数

CC 4.0 BY-SA版权

分类专栏： C/C++ 文章标签：算法 matrix delete c 存储测试

本文链接：https://blog.youkuaiyun.com/pathuang68/article/details/4202331

本文介绍了如何使用C++实现矩阵的逆矩阵计算和矩阵乘法。程序包括计算行列式、Gauss-Jordan消去法求逆矩阵以及矩阵相乘的函数，并通过测试代码展示了实际运算过程。

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

下面的程序借鉴了多位高贤的代码，玄机逸士加以整理和修改而成。

// matrixcomputation.h

//计算行列式，参数一为存储行列式的数组，参数二为阶数
double calculateDeterminant(double *p,int n);

// 使用Gauss-Jordan消去法求n阶实矩阵的逆矩阵
// 返回结果存放在a中，n是矩阵的阶数
int inverseMatrix(double *a, int n);

// 矩阵相乘
// 矩阵a[m*n]乘以矩阵b[n*k]，得到矩阵c[m*k]
// 矩阵c[m*k]为返回结果
void multiplyMatrix(double a[], double b[], int m, int n, int k, double c[]);

// matrixcomputation.cpp

#include "matrixcomputation.h"
#include <math.h>

//计算行列式，参数一为存储行列式的数组，参数二为阶数
double calculateDeterminant(double *p,int n) 
{
   
   
    if(n<1)
    {
   
   
        return -1;
    }
    //若行列式阶数为1，则返回其内的唯一元素值
    if(n==1) 
    {
   
   
        return *p;
    }
    //定义临时指针变量指向当前行列式余子式数组
    double *curr = new double[(n-1)*(n-1)]; 
    double result = 0;
    int i, r, c, cc;
 
    //求当前第一行第i列余子式(按第一行展开)
    for(i = 1; i <= n; ++i) 
    {
   
   
        for(r = 2; r <= n; ++r)
        {
   
   
            cc = -1;
            for(c = 1; c <= n; ++c)
            {
   
   
                if(c != i)
                {
   
   
                    *(curr + (r-2) * (n-1) + (++cc)) = *(p + (r - 1) * n + (c-1));
                }
            }
  &nb