目录
转置矩阵
-
定义
-
性质
/// <summary>
/// 转置矩阵
/// </summary>
/// <param name="arr"></param>
/// <returns></returns>
public float[,] TransposeMatrix(float[,] arr)
{
for (int i = 0; i < arr.GetLength(0); i++)
{
for (int j = i + 1; j < arr.GetLength(1); j++)
{
float temp = arr[i, j];
arr[i, j] = arr[j, i];
arr[j, i] = temp;
}
}
return arr;
}
行列式
-
余子式
-
代数余子式
代数余子式等于对应的余子式的有符号的行列式
-
n阶矩阵求行列式公式
-
意义
在2D中行列式代表着以基向量为两边的平行四边形的有符号面积.在3D环境中则代表着以基向量为三边的平行六面体有符号体积
/// <summary>
/// 求矩阵的行列式
/// </summary>
public static float GetDeterminant(float[,] arr)
{
if (arr.GetLength(0) <= 2) return GetSubArrDeterminant(arr, 0, 0);
float num = 0;
for (int i = 0; i < arr.GetLength(1); i++)
{
num += GetSubArrDeterminant(arr, 0, i) * (float)Math.Pow(-1, i);
}
return num;
}
/// <summary>
/// 得到余子式的行列式
/// </summary>
public static float GetSubArrDeterminant(float[,] arr, int x, int y)
{
if (arr.GetLength(0) == 2)
{
return arr[0, 0] * arr[1, 1] - arr[0, 1] * arr[1, 0];
}
return arr[x, y] * GetDeterminant(GetSubArr(arr, x, y));
}
/// <summary>
/// 得到余子式
/// </summary>
public static float[,] GetSubArr(float[,] arr, int x, int y)
{
if (arr.GetLength(0) <= 1) return arr;
float[,] temp = new float[arr.GetLength(0) - 1, arr.GetLength(1) - 1];
for (int i = 0; i < arr.GetLength(0); i++)
{
if (i == x) continue;
for (int j = 0; j < arr.GetLength(1); j++)
{
if (j == y) continue;
temp[i < x ? i : i - 1, j < y ? j : j - 1] = arr[i, j];
}
}
return temp;
}
伴随矩阵
- 矩阵M的代数余子式矩阵的转置矩阵(挺绕口)
/// <summary>
/// 求矩阵的伴随矩阵
/// 求每个元素的代数余子式,组成一个矩阵
/// </summary>
public static float[,] GetAdjointMatrix(float[,] arr)
{
if (arr.GetLength(0) <= 2) return arr;
float[,] result = new float[arr.GetLength(0), arr.GetLength(1)];
for (int i = 0; i < arr.GetLength(0); i++)
{
for (int j = 0; j < arr.GetLength(1); j++)
{
float[,] temp = GetSubArr(arr, i, j);
result[i, j] = (float)Math.Pow(-1, i + j) * GetDeterminant(temp);
}
}
return result;
}
逆矩阵
- 定义
假设一个方阵M,方阵M的逆,记作M^-1,方阵的逆也是一个矩阵.当M和M^-1相乘的时候,结果是单位矩阵I
- 求解方法
adj M叫做矩阵M的伴随矩阵,定义为矩阵M的代数余子式矩阵的转置矩阵(挺绕口)
public float[,] Inverse(float[,] arr)
{
float det = GetDeterminant(arr);
if (det.Equals(0))
{
Console.WriteLine("行列式为0,没有逆矩阵");
}
else
{
float[,] adjM = GetAdjointMatrix(arr);
for (int i = 0; i < arr.GetLength(0); i++)
{
for (int j = i + 1; j < arr.GetLength(1); j++)
{
arr[j, i] /= det;
}
}
}
return arr;
}
正交矩阵
两个矩阵相加减
即它们相同位置的元素相加减
两个矩阵相乘
乘法只有在第一个矩阵的列数和第二个矩阵的行数相同时才有意义
如果想要理解矩阵变换的几何意义,推荐看一下这个视频