C ++矩阵乘法不同的维度

时间:2013-11-25 19:19:19

标签: c++ matrix multiplication

我想要这样的事情: 将3x3矩阵与3x1矩阵相乘。我尝试将我的1维数组转换为2维。我认为这没关系,但后来在乘法中我得到了错误。但我不确定为什么。这是我的代码。

    // 20  -20  0          150
    //-20   40  -20   x     0
    // 0   -20  30         400  
    int n = 3;
double *tab_P = new double [n];

tab_P[0] = 150;
tab_P[1] = 0;
tab_P[2] = 400;


double **tab_H;
tab_H = new double* [n];
for(int i=0; i<n; i++)
    tab_H[i] = new double[n];

tab_H[0][0] = 0.2;
tab_H[0][1] = 0.15;
tab_H[0][2] = 0.1;


tab_H[1][0] = 0.15;
tab_H[1][1] = 0.15;
tab_H[1][2] = 0.1;

tab_H[2][0] = 0.1;
tab_H[2][1] = 0.1;
tab_H[2][2] = 0.1;

double **tmp_P;
tmp_P = new double *[n];
for(int i=0; i<n; i++)
    tmp_P[i] = new double [1];

for(int i=0; i<n; i++){
    for(int j=0; j<1; j++) {
        tmp_P[i][j] = tab_P[i];
    }
}

double **result = new double * [n];
for(int i=0; i<n; i++)
    result[i] = new double [1];


int i,j,k;
int s = 0;
for(i=0; i<n; i++) {
    for(j=0; j<1; j++) {
        s = 0;
        for(k =0; k<n; k++)
            s += tab_H[i][k] * tmp_P[k][j];
        result[i][j];
    }
}

1 个答案:

答案 0 :(得分:0)

您发布的代码不包含任何内存错误(我使用valgrind进行了检查,这是一个非常有价值的工具,请查看它)。所以,不知何故,你没有将错误复制到你发布的代码中。当您在发布的内容与原始代码之间进行逐字符比较时,您应该会发现错误。

您的代码有一个逻辑错误,但是:您使用int作为变量s,用于对中间产品求和。无论何时为s赋值,浮点乘积都会舍入为整数。这将导致不正确的结果。使用double作为变量s