Question

 for(count = 0; count < max; count ++)
 {
    for (row=(count+1); row < max; row++)
    {
        for(column = 0; column < max; column ++)
        {
            double t = matrix[row][count]/matrix[count][count];
            matrix[row][column] = (matrix[row][column] - (t*matrix[count][column]));
        }
    }
  }

这是我的高斯消除过程的代码，所有变量都是先前声明的。原始矩阵是：

1.4 2.1 2.1 7.4 9.6

1.6 1.5 1.1 0.7 5.0

3.8 8.0 9.6 5.4 8.8

4.6 8.2 8.4 0.4 8.0

2.6 2.9 0.1 9.6 7.7

我得到的输出是：

我执行行交换，因此矩阵现在看起来像：

4.6 8.2 8.4 0.4 8

3.8 8 9.6 5.4 8.8

1.4 2.1 2.1 7.4 9.6

2.6 2.9 0.1 9.6 7.7

1.6 1.5 1.1 0.7 5

我得到的输出是：

4.6 8.2 8.4 0.4 8

0 8 9.6 5.4 8.8

0 0 2.1 7.4 9.6

0 0 0 9.6 7.7

0 0 0 0 5

可以看出，该过程部分工作，消除了每行中的主要术语。但是，它并没有减去剩余的条款。例如矩阵[1] [1]在原始矩阵中，在行交换= 8之后，矩阵[1] [1]在高斯消除后仍为8。

我想知道是否有人会知道我的问题是什么，因为我找不到它。

Answer 1

这一行：

double t = matrix[row][count]/matrix[count][count];

应该在行之前：

for(column = 0; column < max; column ++)

问题是在column循环的第一次迭代之后，您设置了matrix[1][0] = 0;。但是对于列循环的第二次迭代，您使用t计算matrix[1][0]/matrix[0][0]，但现在是0，因为您只需将matrix[1][0]设置为0。

NB。为了帮助我调试added a printf output，以确切了解正在进行的计算。与使用调试器相比，这是原始的，但它是一个有用的工具。

高斯消元不能消除正确的术语

1 个答案: