import cvxpy as cp
# 定义变量
x = cp.Variable()
y = cp.Variable()
# 定义参数（这里我们假设a和b的不确定性可以用区间表示）
# 例如，a在[1, 2]之间，b在[0.5, 1.5]之间
a_min, a_max = 1, 2
b_min, b_max = 0.5, 1.5
# 目标函数（最小化）
objective = cp.Minimize(a_max*x + b_max*y) # 使用最坏情况的参数
# 约束条件
constraints = [x + 2*y >= 1,
x + y >= 1,
x >= 0,
y >= 0]
# 定义问题
prob = cp.Problem(objective, constraints)
# 求解问题
result = prob.solve()
# 输出结果
print("Status:", prob.status)
print("Optimal value (x, y):", x.value, y.value)
print("Optimal objective:", objective.value)
