C++模板矩阵类operator()赋值与operator+重载实现答疑

两个疑惑点的底层原理 #

• m2(0, 0) = 13 赋值实现逻辑
  这个语法的核心是重载函数调用运算符operator()，且必须返回元素的左值引用，才能承接赋值操作。
  你需要在Matrix类中声明两个版本的调用运算符：

  • 非const版本：用于读写元素，返回T&，支持赋值
  • const版本：用于const对象读元素，返回const T&，避免权限放大
    实现时先做行列下标越界校验，再返回一维存储数组中对应位置的元素即可。

• cout << m2 + m3 运算实现逻辑
  这个语句分两步执行，需要分别重载两个运算符：

  1. 先执行m2 + m3：重载成员函数operator+，先校验两个矩阵的行列数完全一致，再逐元素相加，返回一个存了相加结果的新Matrix对象（注意不能返回局部对象的引用）。
  2. 再执行cout << 相加结果：重载全局友元函数operator<<，因为左操作数是ostream类型无法作为Matrix的成员函数实现，函数中按行遍历矩阵元素输出到流，最后返回流对象支持链式输出。

完整代码实现要点 #

Matrix模板类核心结构 #

• 私有成员：
  • size_t _rows, _cols：存储矩阵行列数
  • T* _data：用一维数组存储矩阵元素，比二维动态数组更易管理内存，下标映射规则为第i行j列元素对应下标 i * _cols + j
• 必备基础函数：
  • 构造函数：支持传入行列数初始化零矩阵、支持拷贝构造
  • 析构函数：释放_data指向的堆内存
  • 拷贝赋值运算符：处理自赋值、释放旧内存、拷贝新数据，避免浅拷贝导致的 double free 问题
• 关键运算符实现参考：

// 非const版本的下标访问
template<typename T>
T& Matrix<T>::operator()(size_t row, size_t col) {
    if (row >= _rows || col >= _cols) throw std::out_of_range("Matrix index invalid");
    return _data[row * _cols + col];
}

// const版本的下标访问
template<typename T>
const T& Matrix<T>::operator()(size_t row, size_t col) const {
    if (row >= _rows || col >= _cols) throw std::out_of_range("Matrix index invalid");
    return _data[row * _cols + col];
}

// 矩阵加法
template<typename T>
Matrix<T> Matrix<T>::operator+(const Matrix<T>& other) const {
    if (_rows != other._rows || _cols != other._cols) 
        throw std::invalid_argument("Matrix dimension not match");
    Matrix<T> res(_rows, _cols);
    for (size_t i = 0; i < _rows; ++i) {
        for (size_t j = 0; j < _cols; ++j) {
            res(i,j) = (*this)(i,j) + other(i,j);
        }
    }
    return res;
}

// 流输出运算符，声明为类的友元才能访问私有成员
template<typename T>
std::ostream& operator<<(std::ostream& os, const Matrix<T>& mat) {
    for (size_t i = 0; i < mat._rows; ++i) {
        for (size_t j = 0; j < mat._cols; ++j) {
            os << mat(i,j) << " ";
        }
        os << "\n";
    }
    return os;
}

SymmetricMatrix对称矩阵类实现要点 #

对称矩阵是方阵，且满足a(i,j) = a(j,i)，无需存储全量元素，可以继承Matrix类做空间优化：

• 构造时仅需传入阶数n，实际只需要分配n*(n+1)/2个元素的存储空间（仅存下三角/上三角元素）
• 重写operator()：访问元素时如果row < col，交换row和col的位置再映射到存储下标，赋值时自动同步对称位置的取值，无需手动维护对称特性。
• 加法、流输出逻辑可以直接复用基类的实现，无需重复编写。

内容的提问来源于stack exchange，提问作者Van Wilder