首先，当声明multiply_matrix函数时，你的参数不能是mat_1[N][R1]和mat_2[N][R2]，因为N，R1和R2都没有定义（这些数组的大小是动态的，这是c++不允许的，因此你需要创建一个指向数组的指针，如下所示：（一个月一个月）
假设您要对N进行硬编码，则应该如下声明temp：

std::vector<std::vector<float> > temp(N, std::vector<float>(C2))

相反，您应该使用std：：vector，它是一个动态调整大小的数组，因此您的函数看起来如下所示：

#include <vector>
    //std::vector<std::<vector<float> > is just a vector of vector or a 2 dimensional array of floats
    std::vector<std::vector<float> > multiply_matrix(std::vector<std::vector<float> > mat_1, std::vector<std::vector<float> > mat_2){
    //your code here
    }

你的函数也返回一个浮点数（意味着一个浮点数，而不是一个数组），所以返回类型应该是vector〈vector〉（一个动态的浮点数组）
最后，尝试返回temp[N][C2]，它不是数组，而是浮点数（第N-1行，第C2- 1列），stress_12是数组（也应该是向量）。
希望这对你有帮助。

问题中的代码有多个问题。我将假设矩阵维数N、R1、R2、C2在编译时已知。如果不是这样，则应在问题中说明。我将替换表单中的所有问题

float matrix[M][N]

由

Matrix<M, N> matrix

其中Matrix被定义为

#include <array>

template <std::size_t M, std::size_t N>
using Matrix = std::array<std::array<float, N>, M>;

这样做的好处是，你可以复制你的矩阵，将它作为一个函数的值返回，等等。换句话说，它的工作方式就像任何其他值一样，无论是结构体还是int。
函数multiply_matrix中存在一些矩阵维数不兼容的问题。我将把该函数作为模板函数重新实现，其中矩阵维数是编译时参数。请注意，我更正了mat_1和mat_2的维数，使mat_1的列数与mat_2的行数相同，结果矩阵的行数与mat_1的行数相同和mat_2一样多的列。任何传递维数不兼容的矩阵的尝试都将在编译时被捕获：

template <std::size_t M, std::size_t K, std::size_t N>
Matrix<M, N> multiply_matrix(
  const Matrix<M, K>& mat_1, const Matrix<K, N>& mat_2) {
    Matrix<M, N> temp;
    for (int i = 0; i < M; i++){
        for (int j = 0; j < N; j++) {
            temp[i][j] = 0;
            for (int a = 0; a < K; a++) {
                temp[i][j] += mat_1[i][a]*mat_2[a][j];
            }
        }     
    }
    return temp;
}

另请注意，最后一行是return temp;，而不是return temp[N][C2];。
那我就可以这样用了......

Matrix<3, 3> T /* = ... */;
Matrix<3, 2> stress_12 /* = ... */;
stress_12 = multiply_matrix(T, stress_12);

如果你打算做很多矩阵数学，你可能要考虑使用一个好的矩阵库，比如Eigen，而不是重新发明你自己的东西。