圈养湖羊空间利用率代码

圈养湖羊空间利用率的代码如下：

import numpy as np
from scipy.optimize import linprog

# 目标函数的系数
c = [-365 / T_育]

# 不等式约束矩阵 A_ub * x <= b_ub
A_ub = np.array([[-365 / T_配, 0, 0, 0, 0, 0],  # 配种周期
                 [0, -365 / T_孕, 0, 0, 0, 0],  # 孕期
                 [0, 0, -365 / T_分, 0, 0, 0],  # 分娩期
                 [0, 0, 0, -365 / T_哺, 0, 0],  # 哺乳期
                 [0, 0, 0, 0, -365 / T_育, 0],  # 育肥期
                 [0, 0, 0, 0, 0, -365 / T_休]])  # 休整期

# 不等式约束右侧 b_ub
b_ub = [-365 / T_配 * 2,  # 配种周期
        -365 / T_孕 * 2,  # 孕期
        -365 / T_分,  # 分娩期
        -365 / T_哺,  # 哺乳期
        -365 / T_育,  # 育肥期
        -365 / T_休  # 休整期

# 变量的边界
x_bounds = [(0, C_基础), (0, C_公), (0, C_基础), (0, C_基础), (0, C_基础), (0, C_基础)]

# 求解线性规划问题
res = linprog(c, A_ub=A_ub, b_ub=b_ub, bounds=x_bounds)

# 输出结果
print("最大年化出栏羊只数量：", -res.fun)
print("最优解：")
print("基础母羊数量：", res.x