上篇文章求解了矩阵的行列式,链接为:点击打开链接。
有了行列式,可以求解矩阵的逆矩阵,逆矩阵可以表示为 A^(-1) = A^* / |A|,其中 A^* 为伴随矩阵,|A| 为行列式。
/**计算逆矩阵,求逆矩阵必须是n*n,因此放在行列式中求解
* @param A 要求解的矩阵
* @return 逆矩阵
*/
public double[][] inverse_Matrix(double A[][]){
double result = get_Result(A);
if(result == 0){
System.out.println("行列式为0,不符合计算要求");
return null;
}
double B[][] = new double[A.length][A.length];
// 求伴随矩阵
B = adjoint_Matrix(A.length);
//求逆矩阵
for(int i = 0; i < B.length; i++){
for(int j = 0; j < B.length; j++){
B[i][j] = B[i][j]/result;
}
}
return B;
}
在这里没有判断矩阵的行列是否相等,因为已经在行列式求解,也就是函数get_Result(A)中进行了判断。
现在需要求解伴随矩阵,本例通过 adjoint_Matrix()求解,传入参数为A数组的长度,伴随矩阵的求解方法可以参照该链接的内容:点击打开链接,相关实现代码如下:
/**求伴随矩阵
* @param length 矩阵的行数
* @return 伴随矩阵
*/
private double[][] adjoint_Matrix(int length){
double B[][] = new double[length][length];
//如果是二阶
if(length == 2){
B[0][0] = A[1][1];
B[0][1] = - A[1][0];
B[1][0] = - A[0][1];
B[1][1] = A[0][0];
return B;
}
double sub_A[][] = new double[length - 1][length - 1];
for(int i = 0; i < length; i++){
for(int j = 0; j < length; j++){
//获取当前行列式对应元素的子矩阵
sub_A = get_SubMatrix(A, i, j, A.length);
B[j][i] = Math.pow(-1, (i + j)) * get_Result(sub_A);
}
}
return B;
}
这里用到了get_SubMatrix(),该函数用于求解Aij对应的子矩阵,代码如下:
/**获取当前矩阵某个元素的子矩阵
* @param A 当前矩阵
* @param m 第m行
* @param n 第n列
* @param length 矩阵的行
* @return 子矩阵
*/
public double[][] get_SubMatrix(double A[][], int m, int n, int length){
double B[][] = new double[length - 1][length - 1];
int i_b = 0;
int j_b = 0;
for(int i = 0; i < length; i++){
//去掉第m行元素
if(i == m)continue;
for(int j = 0; j < length; j++){
//去掉第n列元素
if(j == n)continue;
B[i_b][j_b] = A[i][j];
j_b++;
}
i_b++;
j_b = 0;
}
return B;
}