21xrx.com
2025-03-27 20:54:33 Thursday
文章检索 我的文章 写文章
用C++编写矩阵的转置
2023-06-28 12:35:08 深夜i     8     0
C++ 矩阵 转置

在计算机科学中，矩阵是用于进行计算的一种重要数据结构。矩阵由数行数列组成，通常被使用于线性代数、图像处理、机器学习等领域。在许多情况下，对于某些算法和应用程序，需要对矩阵进行转置操作。

矩阵转置可以将矩阵的行变为列，列变为行，即将一个矩阵的主对角线翻折过去后得到的新矩阵。矩阵转置在数学运算和计算机程序中经常使用，它可以使计算机的矩阵操作更加高效、便捷。下面我们使用C++编写一个矩阵转置的程序。

为了表示矩阵和运算，我们先定义一个Matrix类。该类中包含矩阵及其转置的各种操作和运算，并提供构造和析构函数等操作：

class Matrix {
private:
  int n, m;
  double** data;
public:
  // 构造函数
  Matrix(int n, int m);
  // 析构函数
  ~Matrix();
  // 矩阵转置
  Matrix transpose() const;
  // 矩阵元素访问
  double& operator()(int i, int j);
  double operator()(int i, int j) const;
  // 矩阵运算
  Matrix operator+(const Matrix& other) const;
  Matrix operator-(const Matrix& other) const;
  Matrix operator*(const Matrix& other) const;
  // 矩阵输入输出
  friend ostream& operator<<(ostream& os, const Matrix& matrix);
  friend istream& operator>>(istream& is, Matrix& matrix);
};

在上述代码中，我们首先定义矩阵的数据成员，包括行数$n$、列数$m$和一个$n\times m$的二维指针$data$，用于存储矩阵元素。接下来，我们定义了构造函数和析构函数，分别用于创建和销毁矩阵对象。而且，我们还定义了矩阵的转置操作，它将原矩阵的行列互换。我们还定义了一些其他矩阵运算，如加、减、乘等，同时也提供了输入输出运算符。

下面是Matrix类的具体实现：

//构造函数
Matrix::Matrix(int n, int m) : n(n), m(m) {
  data = new double* [n];
  for (int i = 0; i < n; i++) {
    data[i] = new double[m];
  }
}
//析构函数
Matrix::~Matrix() {
  for (int i = 0; i < n; i++) {
    delete[] data[i];
  }
  delete[] data;
}
//矩阵转置
Matrix Matrix::transpose() const {
  Matrix trans(m, n);
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < m; j++) {
      trans.data[j][i] = data[i][j];
    }
  }
  return trans;
}
//矩阵元素访问
double& Matrix::operator()(int i, int j) {
  return data[i][j];
}
double Matrix::operator()(int i, int j) const {
  return data[i][j];
}
//矩阵运算
Matrix Matrix::operator+(const Matrix& other) const {
  Matrix result(n, m);
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < m; j++) {
      result(i, j) = (*this)(i, j) + other(i, j);
    }
  }
  return result;
}
Matrix Matrix::operator-(const Matrix& other) const {
  Matrix result(n, m);
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < m; j++) {
      result(i, j) = (*this)(i, j) - other(i, j);
    }
  }
  return result;
}
Matrix Matrix::operator*(const Matrix& other) const {
  Matrix result(n, other.m);
  for (int i = 0; i < n; i++) {
    for (int j = 0; j < other.m; j++) {
      for (int k = 0; k < m; k++) {
        result(i, j) += (*this)(i, k) * other(k, j);
      }
    }
  }
  return result;
}
// 矩阵输入输出
ostream& operator<<(ostream& os, const Matrix& matrix) {
  for (int i = 0; i < matrix.n; i++) {
    for (int j = 0; j < matrix.m; j++) {
      os << matrix(i, j) << " ";
    }
    os << endl;
  }
  return os;
}
istream& operator>>(istream& is, Matrix& matrix) {
  for (int i = 0; i < matrix.n; i++) {
    for (int j = 0; j < matrix.m; j++) {
      is >> matrix(i, j);
    }
  }
  return is;
}

在上述代码中，我们首先定义了Matrix类的各种函数的具体实现，比如构造函数、析构函数、矩阵转置操作，还有运算符重载函数、输入和输出流运算符。在实现矩阵转置操作时，我们利用了新的矩阵对象tran，它的行数为m，列数为n，首先遍历原矩阵的每一行和列，将其元素拷贝到新矩阵的每一列和行，然后返回新矩阵。

现在，我们可编写主函数对矩阵转置进行测试：

int main() {
  Matrix mat(3, 4);
  cout << "请输入" << 3 * 4 << "个数据:" << endl;
  cin >> mat;
  cout << "原矩阵:" << endl;
  cout << mat;
  cout << "转置后的矩阵:" << endl;
  Matrix trans = mat.transpose();
  cout << trans;
  return 0;
}

在上述代码中，我们首先定义一个3×4的矩阵mat，为其输入元素。接下来，我们输出原矩阵，然后利用transpose函数求出其转置矩阵，最后输出转置后的矩阵。运行上述代码，我们将得到如下结果：

请输入12个数据:
1 2 3 4
5 6 7 8
9 10 11 12
原矩阵:
1 2 3 4
5 6 7 8
9 10 11 12
转置后的矩阵:
1 5 9
2 6 10
3 7 11
4 8 12

我们发现矩阵转置函数得到的输出结果符合预期，证明我们的程序实现了矩阵转置操作。此外，上述实现方式只是一个简单的C++实现，实际上，还有许多高效的实现方式，可以具有更好的性能表现和代码复用性。

  
  

评论区