不使用导入的Java中的矩阵乘法

Question

我试图做一个方法，可以从2个文件中读取输入，每个文件都包含一个矩阵。 目的是检查它们是否可以相乘（如果一个行的长度与另一个行的长度相同），然后创建两个输入矩阵的乘积矩阵。 这是我到目前为止所拥有的。 它不起作用，我没想到它会起作用。 我只想在该方法的逻辑部分获得帮助。

public class MatrixOps {

    public static double[][] multiply(double[][] matrix1, double[][] matrix2) {
        int matrix1Cols = matrix1.length;
        int matrix1Rows = matrix1[0].length;
        int matrix2Cols = matrix2.length;
        int matrix2Rows = matrix2[0].length;
        double[][] productMatrix = new double[0][0];

        if (matrix1Rows == matrix2Cols) {
            productMatrix = new double[matrix1Cols][matrix2Rows];
            for (int i = 0; i < matrix1Cols; i++) {
                for (int j = 0; j < matrix1Rows; j++) {
                    productMatrix[i][j] += matrix1[i][j] * matrix2[j][i];
                    for (int k = 0; k < matrix2Rows; k++) {
                        productMatrix[i][j] += matrix1[i][j] * matrix2[j][k];
                    }
                }
            }
        }

        if (matrix1Cols == matrix2Rows) {
            productMatrix = new double[matrix2Cols][matrix1Rows];
            for (int i = 0; i < matrix2Rows; i++) {
                for (int j = 0; j < matrix1Cols; j++) {
                    productMatrix[i][j] += matrix1[i][j] * matrix2[j][i];
                    for (int k = 0; k < matrix2Rows; k++) {
                        productMatrix[i][j] += matrix1[i][j] * matrix2[j][k];
                    }
                }
            }
        }
        return productMatrix;
    }

}

关于如何使其正常工作的任何想法？

Answer 1

首先，不必在两个循环上都做乘积矩阵：

 for (int i = 0; i < matrix1Cols; i++) {
     for (int j = 0; j < matrix1Rows; j++) {
         productMatrix[i][j] += matrix1[i][j] * matrix2[j][i];
     for (int k = 0; k < matrix2Rows; k++) {
         productMatrix[i][j] += matrix1[i][j] * matrix2[j][k];
     }
 }

这足够了：

for (int i = 0; i < matrix1Cols; i++) {
    for (int j = 0; j < matrix1Rows; j++) {
        for (int k = 0; k < matrix2Rows; k++) {
            productMatrix[i][j] += matrix1[i][j] * matrix2[j][k];
       }
   }
}

同样，如果您要更改矩阵乘法的顺序，那么它必须反映在您的操作中（似乎您只是在进行复制粘贴）

所以在第二部分：

if (matrix1Cols == matrix2Rows) {

矩阵乘法必须更改为：

productMatrix[i][j] += matrix1[i][j] * matrix2[j][k];

至：

productMatrix[i][j] += matrix2[i][j] * matrix1[j][k];