C++对角优势矩阵?


如果对每个矩阵行,行对角项的幅值大于或等于该行每一各非对角项的幅值之和,则称该矩阵为对角优势矩阵。

让我们首先定义一个常量 int 变量 N 值为 3,表示我们的矩阵维度。

const int N = 3;

isDDM(int mat[N][N], int n) 是一个布尔函数,它获取我们矩阵的副本和我们矩阵的大小。在里面,我们使用嵌套 for 循环迭代我们矩阵的行和列。然后,我们计算出每一行每一列的和,并将其添加到我们的 sum 变量中。

bool isDDM(int mat[N][N], int n){
for (int i = 0; i < n; i++){
   int sum = 0;
   for (int j = 0; j < n; j++)
      sum += abs(mat[i][j]);

接下来,我们从当前 sum 中移除对角元素的和。

sum -= abs(mat[i][i]);

接下来,我们检查是否有任何对角元素小于 sum。如果有任何一个元素小于 sum,那么我们就返回 false,退出我们的循环和函数,否则,在完成循环后,因为没有元素小于 sum,我们返回 true。

bool isDDM(int mat[N][N], int n){
for (int i = 0; i < n; i++){
   int sum = 0;
   for (int j = 0; j < n; j++)
      sum += abs(mat[i][j]);
      sum -= abs(mat[i][i]);
   if (abs(mat[i][i]) < sum)
      return false;
   }
   return true;
}

最后,根据返回值,我们在我们的 main 函数中显示它是否是对角优势矩阵。

if(isDDM(mat,matSize)){
   cout << "yes,its a diagonally dominant matrix";
} else {
   cout << "NO, its not a diagonally dominant matrix";
}

示例

我们来看一下以下实现,以检查矩阵是否是或者不是对角优势矩阵。

 在线演示

#include <iostream>
const int N = 3;
using namespace std;
bool isDDM(int mat[N][N], int n){
   for (int i = 0; i < n; i++){
      int sum = 0;
      for (int j = 0; j < n; j++)
         sum += abs(mat[i][j]);
         sum -= abs(mat[i][i]);
         if (abs(mat[i][i]) < sum)
            return false;
   }
   return true;
}
int main(){
   int matSize = 3;
   int mat[N][N] = {
      { 3, -2, 1 },
      { 1, -3, 2 },
      { -1, 2, 4 }
   };
   if(isDDM(mat,matSize)){
      cout << "yes,its a diagonally dominant matrix";
   } else {
      cout << "NO, its not a diagonally dominant matrix";
   }
   return 0;
}

输出

以上代码将生成以下输出 −

yes,its a diagonally dominant matrix

更新于: 16-Jan-2021

315 次浏览

开启你职业生涯

通过完成课程获得认证

入门
广告