在Java中确定Cofactor矩阵

时间:2014-01-16 05:12:40

标签: java matrix linear-algebra

我正在尝试确定辅助因子矩阵。我的代码正确地生成了所有的辅助因子;但是,在某些情况下,将得到的矩阵旋转90度(好吧,切换cols /行)。

例如,矩阵:{{8, 5, 1}, {3, 6, 7}, {5, 6, 6}}产生了正确的结果。

输出>

a
8 3 5 
5 6 6 
1 7 6 

a
-6 17 -12 
-24 43 -23 
29 -53 33 

但是,矩阵:{{1, 0, 5}, {9, 3, 0}, {0, 9, 3}}会切换行和列。

输出>

b
1 0 5 
9 3 0 
0 9 3 

b
9 45 -15 
-27 3 45 
81 -9 3 

正确的结果是:

9 -27 81
45 3 -9
-15 45 3

这样存储矩阵:

Matrix:
    int matrix[][]
    int rows
    int cols

行和列是非常必要的,但每次我尝试确定我正在使用的值时,它比使用matrix.length更好。

这是生成这些矩阵的代码:

public Matrix cofactor() {
    Matrix result = new Matrix(this.rows, this.cols);
    for (int i = 0; i < result.rows; i++) {
        for (int j = 0; j < result.cols; j++) {
            result.matrix[j][i] = (int)(Math.pow(-1, i + j) * removeRowCol(i, j).determinant());
        }
    }

    return result;
}

public Matrix removeRowCol(int row, int col) {
    Matrix result = new Matrix(this.rows - 1, this.cols - 1);

    int k = 0, l = 0;
    for (int i = 0; i < this.rows; i++) {
        if (i == row) continue;
        for (int j = 0; j < this.cols; j++) {
            if (j == col) continue;
            result.matrix[l][k] = this.matrix[i][j];

            k = (k + 1) % (this.rows - 1);
            if (k == 0) l++;
        }
    }

    return result;
}

现在,行列式部分有点像黑客,但它适用于3x3和2x2矩阵。

public int determinant() {
    if (this.rows == 2) return this.matrix[0][0] * this.matrix[1][1] - this.matrix[0][1] * this.matrix[1][0];

    int determinant1 = 0, determinant2 = 0;
    for (int i = 0; i < this.rows; i++) {
        int temp = 1, temp2 = 1;
        for (int j = 0; j < this.cols; j++) {
            temp *= this.matrix[(i + j) % this.cols][j];
            temp2 *= this.matrix[(i + j) % this.cols][this.rows - 1 - j];
        }

        determinant1 += temp; 
        determinant2 += temp2;
    }

    return determinant1 - determinant2;
}

无论如何,我一直试图找出为什么只有某些矩阵被“旋转”。

1 个答案:

答案 0 :(得分:1)

坦率地说,我有点白痴。

我用于测试的矩阵是{{8, 3, 5}, {5, 6, 6}, {1, 7, 6}},但我在Wolfram上检查的值是{{8, 5, 1}, {3, 6, 7}, {5, 6, 6}} ...