C# 矩阵求逆

未来超低端科技研究所

已于 2022-03-30 18:43:38 修改

阅读量2.4k

点赞数 3

分类专栏：九点标定的矩阵运算文章标签：算法矩阵

于 2021-01-18 22:46:21 首次发布

本文链接：https://blog.csdn.net/qq_35044766/article/details/112795413

版权

九点标定的矩阵运算专栏收录该内容

3 篇文章 2 订阅

订阅专栏

一、矩阵求逆算法

本次实现矩阵求逆的算法为伴随矩阵法

代数余子式求逆矩阵：如果矩阵A可逆，则
$A^{-1} = \frac{A^{*}}{|A|}$
(|A| != 0) ,|A|为该矩阵对应的行列式的值

设三阶矩阵A
$A=\left[ \begin{matrix} a & b & c \\ d & e & f \\ g & h & i \end{matrix} \right]$

1、首先求出伴随矩阵
$A^{*}=\left[ \begin{matrix} ei-hf & hc-bi & bf-ce \\ fg-id & ia-cg & cd-af \\ dh-ge & gb-ah & ae-bd \end{matrix} \right]$
将矩阵 A 中的数值带入即可得出伴随矩阵

2、求出其行列式的值
$∣ A ∣ = a e i + b f g + c d h - c e g - b d i - a f h$
3、将结果带入式
$A^{-1} = \frac{A^{*}}{|A|}$
中即可得出矩阵的逆

二、矩阵求逆的代码实现

1、按照上述方式，首先求出矩阵的伴随矩阵

—首先定义一个3X3的二维数组表示一个矩阵—

 double[,] input = new double[3, 3];   //初始矩阵
 double[,] output = new double[3, 3];  //矩阵的逆

———求出伴随矩阵——

 output[0, 0] = input[2, 2] * input[1, 1] - input[2, 1] * input[1, 2];
 output[0, 1] = input[2, 1] * input[0, 2] - input[0, 1] * input[2, 2];
 output[0, 2] = input[0, 1] * input[1, 2] - input[0, 2] * input[1, 1];

 output[1, 0] = input[1, 2] * input[2, 0] - input[2, 2] * input[1, 0];
 output[1, 1] = input[2, 2] * input[0, 0] - input[0, 2] * input[2, 0];
 output[1, 2] = input[0, 2] * input[1, 0] - input[0, 0] * input[1, 2];

 output[2, 0] = input[1, 0] * input[2, 1] - input[2, 0] * input[1, 1];
 output[2, 1] = input[2, 0] * input[0, 1] - input[0, 0] * input[2, 1];
 output[2, 2] = input[0, 0] * input[1, 1] - input[1, 0] * input[0, 1];

———————求出行列式的值———————

 double Avalue   = input[0, 0] * input[1, 1] * input[2, 2] 
 				+ input[0, 1] * input[1, 2] * input[2, 0] 
				+ input[0, 2] * input[1, 0] * input[2, 1]
				- input[0, 2] * input[1, 1] * input[2, 0] 
				- input[0, 1] * input[1, 0] * input[2, 2] 
			    - input[0, 0] * input[1, 2] * input[2, 1];

———————求出逆矩阵 ———————

            for (int i = 0; i < 3; i++)
            {
                for (int j = 0; j < 3; j++)
                {
                    output[i, j] = output[i, j] / Avalue;
                }
            }

得出 output[ i , j ] 即为**矩阵的逆**

以下完整代码

public static double[,] inv3(double[,] input)
        {

            double[,] output = new double[3, 3];

            output[0, 0] = input[2, 2] * input[1, 1] - input[2, 1] * input[1, 2];
            output[0, 1] = input[2, 1] * input[0, 2] - input[0, 1] * input[2, 2];
            output[0, 2] = input[0, 1] * input[1, 2] - input[0, 2] * input[1, 1];

            output[1, 0] = input[1, 2] * input[2, 0] - input[2, 2] * input[1, 0];
            output[1, 1] = input[2, 2] * input[0, 0] - input[0, 2] * input[2, 0];
            output[1, 2] = input[0, 2] * input[1, 0] - input[0, 0] * input[1, 2];

            output[2, 0] = input[1, 0] * input[2, 1] - input[2, 0] * input[1, 1];
            output[2, 1] = input[2, 0] * input[0, 1] - input[0, 0] * input[2, 1];
            output[2, 2] = input[0, 0] * input[1, 1] - input[1, 0] * input[0, 1];

            double Avalue = input[0, 0] * input[1, 1] * input[2, 2]
                + input[0, 1] * input[1, 2] * input[2, 0]
               + input[0, 2] * input[1, 0] * input[2, 1]
               - input[0, 2] * input[1, 1] * input[2, 0]
               - input[0, 1] * input[1, 0] * input[2, 2]
               - input[0, 0] * input[1, 2] * input[2, 1];

            for (int i = 0; i < 3; i++)
            {
                for (int j = 0; j < 3; j++)
                {
                    output[i, j] = output[i, j] / Avalue;
                }
            }
            return output;

        }