在一些需要用单片机、DSP、ARM等系统中用C语言实现高维滤波器的场合,例如卡尔曼滤波器,我们会经常遇到用C语言求解方阵的逆的情况,当阶数为2阶或者3阶的时候我们可以用公式法直接求解,但是当阶数一旦达到4阶或者4阶以上时利用公式求解将会非常的麻烦,而且极易出错,这时我们需要寻求一种可以求解任意阶数的算法。由线性代数的知识可以知道求任意阶数矩阵逆矩阵有2种算法,一种是初等行(列)变换,一种是伴随矩阵法。其中初等行(列)变换由于需要进行行(列)的加减或交换,以及乘上(除以)非零的数这3种步骤,且顺序不固定,因此不适合用C程序进行实现。而伴随矩阵的方法求解步骤固定因此适合用C程序实现。
算法描述
根据线性代数的知识,我们可以知道对于N阶方阵,它的伴随矩阵为