利用高斯消去法计算行列式的值

#include <iostream>

#include <string>

#include <vector>

#include <iomanip>

using namespace std;

 

const int n=3;          //设置方程组的维

 

void input_data(float a[n][n])      //输入行列式的相关数据

{  

    cout<<"输入要计算的行列式a["<<n<<"]["<<n<<"]:"<<endl;

    for(int i=0;i<n;i++)

        for(int j=0;j<n;j++)

            cin>>a[i][j];      

}

 

//

//

//find_line():找出行(i,,,n)、列为i的最大主元素所在的行,并返回行值

//

//

int find_line(int i,float a[n][n])

{

    int max_a_line=i;               //max_a_line为最大元素所在的行

    float max_a=a[i][i];            //max_a为最大元素的值

    int out_i=i;                      //out_i为i的初始值

    for(i=i+1;i<n;i++)

    {

        if(a[i][out_i]<0)

            if(max_a<-a[i][out_i])

            {

                max_a=-a[i][out_i];

                max_a_line=i;

            };

        if(max_a<a[i][out_i])

        {

            max_a=a[i][out_i];

            max_a_line=i;

        };

    }

    return max_a_line;

}

//

//

//chang_line():交换两行的从line列到n列的元素

//

//

void chang_line(int line,int other_line,float a[n][n])

{

    float m;                          //m为中间变量

    for(int j=line;j<n;j++)

    {

        m=a[line][j];

        a[line][j]=a[other_line][j];

        a[other_line][j]=m;

    }

}

 

//

//

//clear_num():消去第i列、i行以下的元素

//

//

void clear_num(int i,float a[n][n])

{

    int line=i,j;

    float m;

    for(i=i+1;i<n;i++)

    {

        m=-a[i][line]/a[line][line];

        for(j=line;j<n;j++)

        {

            a[i][j]=a[i][j]+m*a[line][j];

        };

    }

}

//

//

//print_num():输出行列式的元素

//

//

void print_num(float a[n][n])

{

    cout << setiosflags(ios::fixed);   //和下一句结合起来用来设定小数位的位数

    cout << setprecision(6);

    for(int i=0;i<n;i++)

    {

        cout<<"┃ ";

        for(int j=0;j<n;j++)

            cout<<a[i][j]<<" ";

        cout<<"┃"<<endl;

    }

    cout<<endl<<endl;

}

//

//

//calculate_det():计算行列式的值

//

//

float calculate_det(float a[n][n])

{

    float det=1;

    int i=0,j;

    input_data(a);

    do{

        //print_num(a);

        j=find_line(i,a);      

        if(i!=j)

        {

            chang_line(i,j,a);

            det=-det;

        };

        if(a[i][i]==0)

            return 0;

        det=det*a[i][i];

        if(i<n-1)

            clear_num(i,a);    

        i++;

    }while(i<n);

    return det;

}

 

void main ()

{

    float a[n][n]; 

    cout<<"行列式的值为:"<<calculate_det(a)<<endl;

}