模糊逻辑隶属函数：三角与梯形函数除零错误解决方案问询

问题根源 #

你的梯形隶属函数在输入[0,0,10,20]这类参数时触发除零错误，是因为当a == b时，代码中(x - a)/(b - a)会出现除数为0的情况。这种参数组合对应的是左边界垂直的梯形隶属函数（x在[a,c]区间内隶属度为1，[c,d]区间内从1线性下降到0），不属于标准三角函数的形状，但可以通过特殊处理解决，也能在特定条件下复用三角函数逻辑。

解决方案1：修改梯形函数，直接处理除零情况 #

最直接的方式是在梯形函数中加入参数相等的判断，避免除零运算：

def trap(x, a, b, c, d):
    # 处理左边界重合（a == b）的情况：x >= a时直接进入平台段，左项为1
    if a == b:
        left_segment = 1.0 if x >= a else 0.0
    else:
        left_segment = (x - a) / (b - a) if x >= a else 0.0
    
    # 处理右边界重合（c == d）的情况：x <= d时右项为1
    if c == d:
        right_segment = 1.0 if x <= d else 0.0
    else:
        right_segment = (d - x) / (d - c) if x <= d else 0.0
    
    # 计算最终隶属度
    return max(min(left_segment, 1, right_segment), 0)

修改后的函数可以正常处理[0,0,10,20]这类参数：

• 当x < 0时，返回0
• 当0 <= x <=10时，返回1
• 当10 < x <=20时，返回(20 - x)/10
• 当x >20时，返回0

解决方案2：在特殊场景下复用三角函数逻辑 #

如果一定要用三角函数处理a == b的梯形场景，可以把该梯形拆解为平台段+右三角段，通过判断逻辑复用三角函数：

def trap_with_triangular_fallback(x, a, b, c, d):
    if a == b:
        # 左边界重合时，平台段[a,c]返回1，右段[c,d]用三角函数的下降段逻辑
        if x <= c:
            return 1.0 if x >= a else 0.0
        else:
            # 复用三角函数的下降部分
            return max((d - x)/(d - c), 0)
    elif c == d:
        # 右边界重合时，左段[a,b]用三角函数的上升段逻辑，平台段[b,d]返回1
        if x >= b:
            return 1.0 if x <= d else 0.0
        else:
            return max((x - a)/(b - a), 0)
    else:
        # 正常梯形逻辑
        return max(min((x - a)/(b - a), 1, (d - x)/(d - c)), 0)

额外优化：修复三角函数的除零问题 #

你的三角函数同样存在类似风险（比如a == b或b == c时），可以同步修改：

def triangular(x, a, b, c):
    if a == b:
        # 左边界重合，退化为右三角
        return max((c - x)/(c - b), 0) if x <= c else 0.0
    if b == c:
        # 右边界重合，退化为左三角
        return max((x - a)/(b - a), 0) if x >= a else 0.0
    # 正常三角逻辑
    return max(min((x - a)/(b - a), (c - x)/(c - b)), 0)

内容的提问来源于stack exchange，提问作者Ahmed Abd El-Samie