请求调试Numpy矩阵分类函数：one-to-one/onto等类型判断错误

原代码核心问题 #

1. 非方阵处理错误：非方阵也可能是one-to-one（单射）或onto（满射），原代码直接返回neither完全不符合线性代数规则
2. 浮点精度隐患：直接用det != 0判断行列式是否为0，会因为numpy浮点计算的微小误差（如实际为0却算出1e-16）导致误判
3. 方阵逻辑冗余：方阵的秩等于列数时，必然是可逆矩阵，原代码后续的分支判断完全无效

修正后的代码 #

import numpy as np

def classifier(Matrix):
    rows, cols = Matrix.shape
    rank = np.linalg.matrix_rank(Matrix)
    
    # 可逆矩阵：仅方阵且秩等于阶数（等价于行列式非零）
    if rows == cols:
        if rank == rows:
            return "invertible"
        else:
            # 方阵不满秩时，既不是单射也不是满射
            return "neither"
    else:
        # 非方阵分类逻辑
        if rank == cols:
            # 列满秩 → 单射（one-to-one），此时行数≥列数
            return "one-to-one"
        elif rank == rows:
            # 行满秩 → 满射（onto），此时列数≥行数
            return "onto"
        else:
            # 既不列满秩也不行满秩
            return "neither"

关键逻辑说明 #

• 可逆矩阵：必须是方阵，且秩等于矩阵的阶数（此时既是单射也是满射）
• 非方阵单射：列数等于矩阵的秩，且行数大于列数（比如3×2列满秩矩阵）
• 非方阵满射：行数等于矩阵的秩，且列数大于行数（比如2×3行满秩矩阵）
• 其他情况：包括方阵不满秩、非方阵既不列满秩也不行满秩，返回neither

内容的提问来源于stack exchange，提问作者future ds