矩阵求逆及行列式求值

最新推荐文章于 2021-05-23 07:40:05 发布

寒潭落叶

最新推荐文章于 2021-05-23 07:40:05 发布

阅读量3.5k

点赞数

本文链接：https://blog.csdn.net/wangyx1995/article/details/22813041

版权

本文介绍了一个C语言程序，用于计算矩阵的行列式和逆矩阵。程序中包含了递归调用的行列式计算函数，以及余子式矩阵计算函数。通过输入函数获取矩阵数据，然后计算并显示结果。程序还注意到了内存管理的安全性问题，避免了使用局部变量返回数组指针。此外，程序还包括了矩阵显示和逆矩阵计算功能，如果行列式为0，则矩阵不可逆。

摘要由CSDN通过智能技术生成

正在研究GTK+，争取尽快写出带有框体的程序，这个行列式的程序算是先写着练练手，感受一下递归调用函数吧，应该算是DFS吧，写起来实在是方便。

2014.3.3修改

之前犯了很严重的错误，b[MAXNUM][MAXNUM]数组放在complemet_calculate函数内部，为自动变量，而返回的时候返回的b数组的指针，返回函数体内部声明定义的数组地址是很不安全的，之前能够正确运行的原因是函数结束时释放该自动变量，但该地址所拥有的数据暂时没有改动，系统知识标记了该部分数据可以修改，并不是进行数据的清空，所以指向该处的指针仍然可以正确读取数据。而一旦系统修改了这部分数据，程序就会出现不可预知的错误。

安全的做法应当是在不需要同时存在的情况下，把b数组声明为外部数组，就不会产生上述问题了。

2014.1.3 修改作
1.因为有了行列式计算的函数，于是就写了矩阵求逆的函数加了进去，现在已经可以计算行列式大小，输出伴随矩阵，输出余子式，代数余子式。矩阵的迹很好算，就先没有写，未来再加入计算A^2，计算矩阵的秩，计算矩阵特征值等等。（开始时因为只是想算行列式，输入函数写的时候只能输入为方阵，为了让函数有广泛适用性，添加着几个功能的时候应该还要修改一下输入函数。）
2.发现代数余子式的正负号似乎是错的，改了一下

#include "stdio.h"
#include"math.h"

#include "conio.h"
#define MAXNUM100 //行列式最大大小

float b[MAXNUM][MAXNUM];//complemet_calculate函数的临时数组，声明为外部变量防止数据丢失

float* complement_calculate(float a[][MAXNUM],int n,int x,int y)//余子式矩阵计算函数，返回一个n-1阶二维数组指针
{
    inti,j,m,p;
   memset(b,0,n-1);
   for(i=0,m=0;i<=n-1;m++,i++)
    {
       if(