2023 年数学建模国赛D 题 圈养湖羊的空间利用率Matlab代码

问题一:

% 定义线性规划问题
% 目标函数系数矩阵，此处为最大化年化出栏羊只数量
f = -1; 

% 不等式约束矩阵 A 和不等式约束右侧 b
% 这些约束基于问题描述中的条件
A = [
    1, -365/T, -2*365/T, 0, 0;
    0, 0, 0, -365/T, -P_l*365/T;
    0, 0, 0, 0, -P_k*365/T;
    -1, 365/T, 2*365/T, 365/T, P_l*365/T;
    0, -365/T, -2*365/T, -365/T, -P_l*365/T;
    0, 365/T, 0, 0, 0;
    0, 0, -365/T, 0, 0;
    0, 0, 0, 365/T, 0;
    0, 0, 0, 0, 365/T;
    P_p, 0, 0, 0, 0;
    -1, 0, 0, 0, 0;
    0, 0, 0, 0, -4;
    0, 0, 0, 0, 14;
    0, 0, 0, 0, 8;
    0, 0, 0, 0, 6;
];

b = [
    0;
    0;
    0;
    N_b;
    0;
    0;
    0;
    0;
    0;
    N_r;
    -N_r;
    4*N_s;
    -14*N_s;
    -8*N_s;
    -6*N_s;
];

% 不等式约束的上下界
lb = zeros(5, 1);  % 变量的下界，对应 Nb, Nr, Np, Nm, Nl
ub = Inf(5, 1);    % 变量的上界，可以是正无穷

% 线性规划求解
options = optimoptions('linprog', 'Displ