C#计算伪逆矩阵（PseudoInverse）

最新推荐文章于 2024-09-17 20:34:02 发布

汪汪星河

最新推荐文章于 2024-09-17 20:34:02 发布

阅读量949

点赞数 1

分类专栏： C# 文章标签： c# 数学矩阵

本文链接：https://blog.csdn.net/jemi_a/article/details/118933556

版权

C# 专栏收录该内容

14 篇文章 0 订阅

订阅专栏

这篇博客介绍了如何利用MathNet.Numerics库中的Matrix.PseudoInverse()方法来计算矩阵的伪逆，与Python的numpy.linalg.pinv和Matlab的pinv函数等效。示例中给出了具体的矩阵xEX，并展示了计算得到的伪逆矩阵pia的详细数值，验证了计算的正确性。

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

我们用的是MathNet.Numerics命名空间里面的Matrix.PseudoInverse()方法计算伪逆，跟Python的numpy.linalg.pinv或者Matlab的pinv是一样的。

            double[,] xEX = {
                { 1.0, 0,0 ,0,0,6.0},
                { 0, -1,0 ,0,0,17.0},
                { 0, 0,1 ,-6,17,0},

                { 1.0, 0,0 ,0,0,-6.0},
                { 0, -1,0 ,0,0,17.0},
                { 0, 0,1 ,6,17,0},

                { 1.0, 0,0 ,0,0,6.0},
                { 0, -1,0 ,0,0,-17.0},
                { 0, 0,1 ,-6,-17,0},

                { 1.0, 0,0 ,0,0,-6.0},
                { 0, -1,0 ,0,0,-17.0},
                { 0, 0,1 ,6,-17,0},
            };
            var matrix6 = DenseMatrix.OfArray(xEX);
            var pi = matrix6.PseudoInverse();
            var pia = pi.ToArray();
            double[,] xEX_PseudoInverse ={
                { 0.2500, 0,0 ,0.2500,0,0,0.2500,0,0,0.2500,0,0},
                { 0, -0.2500,0 ,0,-0.2500,0,0,-0.2500,0,0,-0.2500,0},
                { 0, 0,0.2500 ,0,0,0.2500,0,0,0.2500,0,0,0.2500},
                { 0,0,-0.0417,0,0,0.0417,0,0,-0.0417,0,0,0.0417},
                { 0,0,0.0147,0,0,0.0147,0,0,-0.0147,0,0,-0.0147},
                { 0.0046,0.0131,0,-0.0046,0.0131,0,0.0046,-0.0131,0,-0.0046,-0.0131,0},
            };
            //可以通过resultStr检查变换结果
            int row = pia.GetLength(0);  //可以理解为：第一维的长度（即行数)
            int col = pia.GetLength(1);  //可以理解为：第二维的长度（即列数）
            string resultStr = "";
            for (int i = 0; i < row; i++)
            {
                for (int j = 0; j < col; j++)
                {
                    resultStr += (pia[i, j].ToString("f4"));
                    resultStr += "\t";
                }
                resultStr += "\n";
            }