如何实现支持加减乘运算的C#泛型Matrix类？

一、核心问题分析 #

1. LineMultiplyCol方法运算错误：C#泛型无法直接对未约束的T使用+、*运算符，编译器无法推断T是否支持这些操作。
2. 矩阵乘法运算符未实现：需要按照矩阵乘法规则，完成维度校验、元素点积计算的逻辑。

二、分步解决方案 #

1. 完善泛型类型约束（.NET 6+推荐方案） #

利用.NET 6引入的运算符接口，明确约束T必须支持加法、乘法、减法运算，以及拥有加法单位元（替代default(T)，保证运算初始值的正确性）：

using System.Numerics;

class Matrix<T> 
    where T : IAdditionOperators<T, T, T>, 
              IMultiplyOperators<T, T, T>, 
              ISubtractionOperators<T, T, T>,
              IAdditiveIdentity<T, T>
{
    // 原有属性、构造函数、索引器代码保持不变
}

2. 修复LineMultiplyCol方法 #

使用加法单位元替代default(T)，确保初始值符合加法运算的要求（比如数值类型的0，自定义类型的加法单位元素）：

private static T LineMultiplyCol(List<T> l, List<T> c)
{
    T result = T.AdditiveIdentity; // 用加法单位元替代default(T)，更准确
    for (int i = 0; i < l.Count; i++)
    {
        result += l[i] * c[i];
    }
    return result;
}

3. 实现矩阵乘法运算符* #

按照矩阵乘法规则：

• 校验维度合法性：第一个矩阵的列数必须等于第二个矩阵的行数
• 初始化结果矩阵（行数为第一个矩阵的行数，列数为第二个矩阵的列数）
• 遍历每个元素，计算对应行与列的点积

代码实现：

// 新增获取矩阵列的辅助方法
private List<T> GetColumn(int columnIndex)
{
    List<T> column = new List<T>();
    foreach (var line in Lines)
    {
        column.Add(line[columnIndex]);
    }
    return column;
}

public static Matrix<T> operator *(Matrix<T> m1, Matrix<T> m2)
{
    // 维度合法性检查
    if (m1.Lines.Count == 0 || m2.Lines.Count == 0 || m1[0].Count != m2.Lines.Count)
    {
        throw new InvalidOperationException("矩阵维度不匹配，无法执行乘法运算");
    }

    int resultRows = m1.Lines.Count;
    int resultCols = m2[0].Count;
    Matrix<T> result = new Matrix<T>();

    // 初始化结果矩阵的行
    for (int i = 0; i < resultRows; i++)
    {
        result.Lines.Add(new List<T>(new T[resultCols]));
    }

    // 计算每个元素的值
    for (int i = 0; i < resultRows; i++)
    {
        for (int j = 0; j < resultCols; j++)
        {
            result[i][j] = LineMultiplyCol(m1[i], m2.GetColumn(j));
        }
    }

    return result;
}

4. 兼容旧版本.NET的备选方案（动态类型） #

如果使用.NET 6以下版本，无法使用运算符接口，可以通过dynamic绕过编译检查（注意：此方式失去编译时类型安全，运行时若T不支持运算会报错）：

// 修改LineMultiplyCol方法
private static T LineMultiplyCol(List<T> l, List<T> c)
{
    dynamic result = default(T);
    for (int i = 0; i < l.Count; i++)
    {
        result += (dynamic)l[i] * (dynamic)c[i];
    }
    return result;
}

// 泛型约束可简化为（仅作提示，无实际编译约束）
class Matrix<T> where T : struct // 可选，限制为值类型
{
    // 原有代码保持不变
}

三、使用示例 #

// 整数矩阵乘法示例
var m1 = new Matrix<int>();
m1.Lines.Add(new List<int> { 1, 2 });
m1.Lines.Add(new List<int> { 3, 4 });

var m2 = new Matrix<int>();
m2.Lines.Add(new List<int> { 5, 6 });
m2.Lines.Add(new List<int> { 7, 8 });

Matrix<int> product = m1 * m2;
// 结果矩阵为：
// [19, 22]
// [43, 50]

内容的提问来源于stack exchange，提问作者Mahdi Jafari