在 C++ 中使用向量进行二维或多维矩阵乘法

Question

我需要声明一个包含 3 x 3 矩阵的向量。 并且需要将该向量传递给 function 并且需要使用相同的向量执行乘法运算。 我知道，我走的路不对。 但我无法弄清楚问题所在。 我希望有关 2d 矩阵向量乘法的一些方向。

#include <iostream>
#include <vector>
using namespace std;
void display(vector<vector<int>> &vec)
{
    for (int i = 0; vec.size(); i++)
    {
        for (int j = 0; j < vec.size(); j++)
        {
            cout << vec[i][j] << " ";
        }
        cout << endl;
    }
}
vector<vector<int>> matMultiply(vector<vector<int>> &vec)
{
    for (int i = 0; vec.size(); i++)
    {
        for (int j = 0; j < vec.size(); j++)
            vec[i][j] = 0;
        {
            for (int k = 0; k < vec.size(); i++)
            {
                vec[i][j] *= vec[i][k] * vec[k][j];
                cout << vec[i][j] << " ";
            }
        }
        cout << endl;
        return vec;
    }
}
int main()
{
    vector<vector<int>> matrix = {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}};
    display(matrix);
    matMultiply(matrix);
    return 0;
}

Answer 1

因为矩阵乘法需要多个矩阵值来计算单个索引处的值，所以您需要某种方式来存储中间结果。 可能在另一个vector中。 您需要在matMultiply function 中创建一个新的std::vector<std::vector<int>> 。 它应该像这样：

std::vector<std::vector<int>> matMultiply(std::vector<std::vector<int>> &vec){
    std::vector<std::vector<int>> result(3,vector<int>(3));
    for(int i = 0; result.size(); i++){
        for (int j = 0; j < result.size(); j++) {
            result[i][j] = 0;
            
            for (int k = 0; k < result.size(); i++){
                result[i][j] *= vec[i][k] * vec[k][j];
                std::cout << result[i][j] << " ";
            }
        }   
        std::cout << std::endl;
        return result;    
    }
}

注意：对于以后的问题，请在此处发布之前正确缩进您的代码。

Answer 2

@Taimoor Zaeem。 再次感谢你。 我按照你的建议做了代码。 但结果是错误的。 在显示器旁边 function 不工作，可能是我做错了.. 为了查看结果，我添加了std:: cout<< mat[i][j]<<" "; 在matMultiply结果来了

1 9 30 
2 12 36 
3 15 42 
4 24 66 
8 33 81 
12 42 96 
7 39 102 
14 54 126
21 69 150

结果将是

30 36 42 
66 81 96
102 126 150

#include<iostream>
#include<vector>


void display (std::vector<std::vector<int>> mat){
     for(int i = 0; mat.size(); i++){
        for (int j = 0; j < mat.size(); j++) {
            std :: cout<< mat[i][j]<<" ";
        }
        std :: cout<< std ::endl;   
    }
}
std::vector<std::vector<int>> matMultiply(const std::vector<std::vector<int>> &vec){
    std::vector<std::vector<int>> result(3, std::vector<int>(3));
    for(int i = 0; result.size(); i++){
        for (int j = 0; j < result.size(); j++) {
            for (int k = 0; k < result.size(); k++){
                result[i][j] += vec[i][k] * vec[k][j];
                std :: cout<< result[i][j] <<" ";
            }
            std :: cout<< std:: endl;
        }   
    }
    return result;
}

int main(){
    
    std :: vector < std :: vector< int > > matrix = {{1,2,3},{4,5,6},{7,8,9}};
    //display(matrix);
    matMultiply(matrix);
    display(matMultiply(matrix));
    return 0;
}