C++中的Gauss-Jordan消除

Question

我在C++中使用Gauss-Jordan消除来求解线性方程组.代码工作正常.想知道为什么1,2,3 void gauss()号线不能被4号线取代(这样做后输出不正确)？

#include <iostream>
using namespace std;
class Gauss
{
    float a[50][50];
    int n;
public:
    void accept()
    {
        cout<<"Enter no. of variables: ";
        cin>>n;
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n+1;j++)
            {
                if(j==n)
                    cout<<"Constant no."<<i+1<<" = ";
                else
                    cout<<"a["<<i+1<<"]["<<j+1<<"] = ";
                cin>>a[i][j];
            }
        }
    }
    void display()
    {
        for(int i=0;i<n;i++)
        {
            cout<<"\n";
            for(int j=0;j<n+1;j++)
            {
                if(j==n)
                    cout<<" ";
                cout<<a[i][j]<<"\t";
            }
        }
    }

    void gauss()//converting augmented matrix to row echelon form
    {
        float temp;//Line 1
        for(int i=0;i<n;i++)
        {
            for(int j=i+1;j<n;j++)
            {
                temp=a[j][i]/a[i][i];//Line 2
                for(int k=i;k<n+1;k++)
                {
                      a[j][k]-=temp*a[i][k];//Line 3
                    //a[j][k]-=a[j][i]*a[i][k]/a[i][i];//Line 4
                }
            }
        }
    }

    void EnterJordan()//converting to reduced row echelon form
    {
        float temp;
        for(int i=n-1;i>=0;i--)
        {

            for(int j=i-1;j>=0;j--)
            {
                temp=a[j][i]/a[i][i];
                for(int k=n;k>=i;k--)
                {
                    a[j][k]-=temp*a[i][k];
                }
            }
        }

        float x[n];
        for(int i=0;i<n;i++)//making leading coefficients zero
            x[i]=0;
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n+1;j++)
            {
                if(x[i]==0&&j!=n)
                    x[i]=a[i][j];
                if(x[i]!=0)
                    a[i][j]/=x[i];
            }
        }
    }
    void credits()
    {
        for(int i=0;i<n;i++)
        {
            cout<<"\nx"<<i+1<<" = "<<a[i][n]<<endl;
        }
    }

};

int main()
{
    Gauss obj;
    obj.accept();
    cout<<"\n\nAugmented matrix: \n\n\n";
    obj.display();
    obj.gauss();
    cout<<"\n\nRow Echelon form: \n\n\n";
    obj.display();
    obj.EnterJordan();
    cout<<"\n\nReduced row echelon form:\n\n\n";
    obj.display();
    cout<<"\n\nSolution: \n\n\n";
    obj.credits();
    return 0;
}

注意:我的代码没有考虑到枢轴为零时的除法问题(我每次都选择对角线元素作为枢轴).然而,对于我尝试的特定示例,没有遇到这种情况.

增广矩阵是:

 2   1  -1    8 
-3  -1   2   -11    
-2   1   2   -3

输出矩阵是:

1   0   0    2  
0   1   0    3  
0   0   1    -1 

解决方案是:

x1 = 2

x2 = 3

x3 = -1

使用第4行,输出矩阵是:

1   0   0    -0.75  
0   1   -0   8  
0   0   1    -1.5   

解决方案是:

x1 = -0.75

x2 = 8

x3 = -1.5

Answer 1

您的第 4 行多次读取a[j][i]，第一次通过内部循环时，k == i将更改为a[j][i]0.0f，从而中断下一次n-i迭代。

对变量的读取和对同一位置的写入进行重新排序并不安全。