约束条件下LR公式最大分母对应x,y,z取值的Python代码求助

你的需求是找到x、y、z的取值组合，让公式LR = 500/(30x + 40y + 10z)的分母最大，约束条件为：

• x、y、z的取值范围都是0.8~1.2
• 分母不能超过500

你当前编写的代码存在几个明显问题：

• Python内置的range仅支持生成整数序列，无法处理0.8到1.2的浮点数区间
• 变量声明顺序错误，x、y、z还未定义就用来计算p
• 约束条件p < 500只是单独写了一句表达式，没有实际的判断逻辑
• 代码里的分母公式和你需求里的30x + 40y + 10z不匹配，返回值也不是优化后的最优结果

由于分母是线性表达式，三个变量的系数都是正数，理论上要让分母最大，x、y、z直接都取最大值1.2即可，我们可以先验证该结果是否满足分母≤500的约束，也可以通过遍历所有取值的方式做通用验证。

def find_optimal_xyz(step=0.01):
    # 变量取值边界
    min_limit = 0.8
    max_limit = 1.2
    max_denominator = 0
    best_x = best_y = best_z = None

    # 按指定步长遍历所有可能的取值组合
    x = min_limit
    while x <= max_limit + 1e-8:  # 加极小值避免浮点数精度问题漏掉最大值
        y = min_limit
        while y <= max_limit + 1e-8:
            z = min_limit
            while z <= max_limit + 1e-8:
                current_denominator = 30 * x + 40 * y + 10 * z
                # 筛选符合约束且分母更大的组合
                if current_denominator <= 500 and current_denominator > max_denominator:
                    max_denominator = current_denominator
                    best_x, best_y, best_z = x, y, z
                z += step
            y += step
        x += step
    
    if not best_x:
        return None
    lr_value = 500 / max_denominator
    return best_x, best_y, best_z, max_denominator, lr_value

# 调用测试
result = find_optimal_xyz(step=0.01)
if result:
    x, y, z, d, lr = result
    print(f"最优取值：x={x:.2f}, y={y:.2f}, z={z:.2f}")
    print(f"最大分母：{d:.2f}，对应LR值：{lr:.4f}")
else:
    print("无符合约束的解")

最优取值：x=1.20, y=1.20, z=1.20
最大分母：96.00，对应LR值：5.2083

可以看到三个变量都取最大值时的分母仅为96，远小于500的约束上限，因此该组合就是最优解。

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