.NET用行列演算クラスライブラリ（行列の積、逆行列、転置行列、擬似逆行列など）

// 行列(matA)の設定（二次元配列で行列を設定）
var matA = new double[,] {
    { 1, 2, 3 },
    { 4, 5, 6 },
    { 7, 8, 9 }
};
matA.Print("matA =");　// コンソールへ行列の表示

// 行列(matB)の設定
var matB = new double[,] {
    { 1, 4, 8 },
    { 6, 2, 5 },
    { 9, 7, 3 }
};
matB.Print("matB =");

// 行列(matA)と行列(matA)の積
var matMult = matA.Mult(matB);
matMult.Print("matA matB =");

// 行列(matA)と行列(matA)の加算
var matAdd = matA.Add(matB);
matAdd.Print("matA + matB =");

// 行列(matA)と行列(matA)の減算
var matSub = matA.Sub(matB);
matSub.Print("matA - matB =");

// 行列(matB)の逆行列
var matInverse = matB.Inverse();
matInverse.Print("(matB)-1 =");

// 行列(matA)の転置
var matTranspose = matA.Transpose();
matTranspose.Print("(matA)T =");

// 回転行列の取得
var matRotate = Mat.RotateMat(10.0);
matRotate.Print("RotateMat =");

// 拡大行列の取得
var matScale = Mat.ScaleMat(2.0, 5.0);
matScale.Print("ScaleMat =");

// 平行移動行列の取得
var matTranslate = Mat.TranslateMat(-5.0, 12.0);
matTranslate.Print("TranslateMat =");

Bitmap bmp;
// 画像ファイルを開く
var imgMat = Mat.Imread("0.bmp", out bmp);
imgMat.Print("Image data =", false);

// 画像の表示
var window = bmp.Imshow("Image");

// 配列(画像データ)をcsvファイルに保存
imgMat.SaveCsv("matData.csv");

// 行列(matA)の設定例
var matA = new double[3, 3];
matA[0, 0] = 1; matA[0, 1] = 2; matA[0, 2] = 3;
matA[1, 0] = 4; matA[1, 1] = 5; matA[1, 2] = 6;
matA[2, 0] = 7; matA[2, 1] = 8; matA[2, 2] = 9;

// 行列(matB)の設定例
var matB = new double[,] {
    { 1, 4, 8 },
    { 6, 2, 5 },
    { 9, 7, 3 }
};

var matA = new double[,] {
    { 1, 2, 3 },
    { 4, 5, 6 },
    { 7, 8, 9 }
};
matA.Print();

matA.Print("matA =");

// 先頭行付き表示
matA.Print("matA =", false);
// 先頭行なしの場合
matA.Print(false);

var matA = new double[,] {
    { 1, 2, 3 },
    { 4, 5, 6 },
    { 7, 8, 9 }
};

var matB = new double[,] {
    { 1, 4, 8 },
    { 6, 2, 5 },
    { 9, 7, 3 }
};

// 行列の積（従来の書き方）
var mat1 = Mat.Mult(matA, matB);

// 行列の積（拡張メソッドの書き方）
var mat2 = matA.Mult(matB);

// 逆行列（従来の書き方）
var mat3 = Mat.Inverse(matB);

// 逆行列（拡張メソッドの書き方）
var mat4 = matB.Inverse();

Bitmap bmp;
//画像ファイルを開く
var matImg = Mat.Imread("0.bmp", out bmp);
// 画像データの表示
matImg.Print(false);
// 画像の表示
bmp.Imshow("画像");
// CSVファイル保存
matImg.SaveCsv("image.csv");

// 画像の表示
bmp.Imshow("画像");
bmp.Imshow("画像");
bmp.Imshow("画像");

// ウィンドウの生成
var win = new Mat.NamedWindow("画像");
// 画像の表示
bmp.Imshow(win);
bmp.Imshow(win);
bmp.Imshow(win);

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;

using ImagingSolution;

namespace ConsoleApplication1
{
    class Program
    {
        static void Main(string[] args)
        {
            // 近似データ
            var dataX = new double[] { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 };
            var dataY = new double[] { 10, 9, 6, 12, 15, 21, 34, 36, 48, 52 };

            var mat = new double[dataX.Length, 3];

            for (int i = 0; i < dataX.Length; i++)
            {
                mat[i, 0] = dataX[i] * dataX[i];
                mat[i, 1] = dataX[i];
                mat[i, 2] = 1;
            }
            mat.Print();

            // 配列から行列（N行1列）へ変換
            var matY = dataY.ArrayToMat();
            matY.Print();

            // 擬似逆行列で二次式の係数を解く
            var ans = mat.PseudoInverse().Mult(matY);

            // 解の表示
            ans.Print();

            // 近似式を表示
            Console.WriteLine($"y = {ans[0,0]}x^2 + {ans[1, 0]}x + {ans[2, 0]}");

            // キー入力待ち
            Console.WriteLine("続行するには何かキーを押してください . . .");
            Console.ReadKey();
        }
    }
}