分段错误 C++ 矩阵

Question

我想让任意矩阵[n][n+1]变成上三角矩阵[n][n+1]。

我做了这段代码，但这会导致分段错误。

void diagonalizarMatriz(float** Matriz, int n){
for(int i = 0; i < n-1; i++)
    for(int k = 0; k < n; k ++)
        for(int j = n; j >= i; j++)
            Matriz[k][j] = Matriz[k][j] - ((Matriz[k][i] * Matriz[i][j]) / Matriz[i][i]);
}

int main(){
float** Matriz = new float* [3];
for(int i = 0; i < 3 ; i++)
    Matriz[i] = new float [4];

//test matrix
Matriz[0][0] = 1;
Matriz[0][1] = 4;
Matriz[0][2] = 52;
Matriz[0][3] = 57;
Matriz[1][0] = -27;
Matriz[1][1] = -110;
Matriz[1][2] = -3;
Matriz[1][3] = -134;
Matriz[2][0] = 22;
Matriz[2][1] = 2;
Matriz[2][2] = 14;
Matriz[2][3] = 38;

diagonalizarMatriz(Matriz, 3);

Answer 1

这里

for(int j = n; j >= i; j++)

您从数组维度的上边界的 n 开始并向上计数，
很快你就可以访问你的阵列之外的东西，如果你幸运的话，这会给你带来一个段错误。

随意猜测

for(int j = n; j >= i; j--)

倒计时。

Answer 2

更新

假设您的问题是关于矩阵对角化（正如 Eugene 在评论部分中声称的那样）：

除了Yunnosch指出的，对角化矩阵意味着矩阵必须是正方形n x n。但是，在 main 中，您将其初始化为 3 x 4 矩阵。

原码

float** Matriz = new float* [3];
for(int i = 0; i < 3 ; i++)
  Matriz[i] = new float [4];

为了摆脱段错误，在main中更改以下部分（将矩阵设置为3 x 3）：

float** Matriz = new float* [3];
for(int i = 0; i < 3 ; i++)
  Matriz[i] = new float [3];

//test matrix
Matriz[0][0] = 1;
Matriz[0][1] = 4;
Matriz[0][2] = 52;
Matriz[1][0] = -27;
Matriz[1][1] = -110;
Matriz[1][2] = -3;
Matriz[2][0] = 22;
Matriz[2][1] = 2;
Matriz[2][2] = 14;

最后得到如下矩阵（下三角）：

1 0 0
4 5 0
7 8 9

从第三个嵌套循环中删除等号：

for(int j = n; j > i; j--)

我假设您可以从这里开始工作，使其成为upper-triangular matrix。

旧答案

试试这个：

int matrix[3][3] = {{1,2,3},{4,5,6},{7,8,9}};
int row = 3;
int col = 3;
int i, j; 
for (i = 0; i < row; i++) { 
  for (j = 0; j < col; j++) {
    if (i > j) { 
        cout << "0" << " "; 
    } else
        cout << matrix[i][j] << " "; 
  } 
  cout << endl; 
}

会给你这个矩阵

1 2 3
0 5 6
0 0 9