C++ 求一个矩阵的行列式的值（星星笔记）

最新推荐文章于 2024-10-12 21:25:47 发布

GrowthDiary007

最新推荐文章于 2024-10-12 21:25:47 发布

阅读量2.2w

点赞数 6

分类专栏：算法 C++ 文章标签： C++ 算法矩阵运算矩阵行列式

本文链接：https://blog.csdn.net/XX_123_1_RJ/article/details/39267665

版权

算法同时被 2 个专栏收录

189 篇文章 3 订阅

订阅专栏

C++

15 篇文章 0 订阅

订阅专栏

本文介绍了一个C++程序，用于将输入的方阵通过初等行变换转化为上三角矩阵，并计算该矩阵的行列式值。程序首先接收用户输入的矩阵元素，然后通过一系列的行变换操作达到目标形式，并最终输出行列式的值及变换后的矩阵。

摘要由CSDN通过智能技术生成

输入矩阵要求为方阵（n×n）

原理：将一个矩阵经过初等行变换化为上三角矩阵。

#include<iostream>  
#include <iomanip>  
using namespace std;  
int const n=3;  //确定矩阵的节数  
/* 
作者 星星笔记 
 
*/  
int main()  
{   
    double fun(double array[n][n]);  
    void temp(double aa[],double bb[],int n);  
    double a[n][n];  
    double det;  
    int i,j,kk=0;  
    cout<<"请输入一个"<<n<<"节方阵"<<endl;  
    for(i=0; i<n; i++)    //一个矩阵的输入  
        for(j=0; j<n; j++)  
            cin>>a[i][j];  
      det = fun(a);       //获取行列式的值  
      cout<<"矩阵的行列式的值为："<<det<<endl;  
      cout<<"转换的上三角矩阵为："<<endl;  
      for(i=0; i<n; i++)  
          for(j=0; j<n; j++)  
          {  
              cout<<setw(7)<<a[i][j];  
              kk = kk + 1;  
              if( kk % n == 0 )   
                  cout<<endl;  
          }  
     cout<<endl;  
          return 0;    
}  
//交换数组指定的两行，即进行行变换（具体为行交换）  
void temp(double aa[],double bb[],int n)  
{     
    int i;  
    double temp1;  
    for(i=0 ; i<n ; i++)  
    {  
        temp1 = aa[i];  
        aa[i] = bb[i];  
        bb[i] = temp1;    
    }    
}  
  
double fun(double array[n][n])  
{     
    int ii,jj,k,u;  
	int iter = 0;  //记录行变换的次数（交换）
    double det1=1,yin;  

    for(ii=0 ; ii<n; ii++)  
    {     
        if(array[ii][ii] == 0) 
			
        for(jj=ii; jj<n; jj++)  
        {    
            if(array[jj][ii] != 0)  
			{
				temp(array[ii],array[jj],n);//交换两行  
				iter ++;
			}
        }  

        for(k=ii+1; k<n; k++)  
		{    
          yin = -1 * array[k][ii] / array[ii][ii] ;  

          for(u=0; u<n; u++)  
          {   
              array[k][u] = array[k][u] + array[ii][u] * yin;  
          }  
		}  
   }  
  for(ii=0; ii<n; ii++)  //求对角线的积 即 行列式的值
      det1 = det1 * array[ii][ii];  
  //行变换偶数次符号不变
  if(iter%2 == 1)
	  det1= -det1;
   return (det1);  
}

运行结果如下：