表1木板的尺寸 木板 长度(mm) 宽度(mm) S1 3000 1500 表2 产品尺寸及生产任务 产品名称 长度(mm) 宽度(mm) 生产任务(件) 利润(元/件) P1 373 201 774 19.9 P2 477 282 2153 23.0 P3 406 229 1623 21.0 P4 311 225 1614 16.0 有关问题模型建立如下：（其中X1、X3分别为P1在一块木板上切割后的数量、P3在一块木板上切割后的数量，B为原木板的数量） min 3000*1500-373*201*X1-406*229*X3 约束条件： B*X1>=774 B*X3>=1623 3000*1500>=373*201*X1+406*229*X3 B>=35 其中B,X1，X3为非负整数 另一个约束条件： min z = B B*X1>=774 B*X3>=1623 3000*1500-373*201*X1-406*229*X3<=9000 B>=35 其中B,X1，X3为非负整数 写出相关解题过程及代码实现

这是一个整数规划问题，可以使用MATLAB中的intlinprog函数求解。具体过程如下：

首先定义目标函数系数向量和约束矩阵、约束右侧向量：

f = [0, 0, 0, -1];
Aeq = [373, 0, 406, 0; 0, 0, 0, 1; 0, 0, 0, 1; 0, 1, 0, 0];
beq = [0; 774; 1623; 35];
lb = [0; 0; 0; 0];
intcon = [1; 2; 3; 4];

其中，f为目标函数系数向量，Aeq为等式约束矩阵，beq为等式约束右侧向量，lb为变量下界，intcon为整数变量的索引。

然后调用intlinprog函数求解：

[x,fval] = intlinprog(f,intcon,[],[],Aeq,beq,lb);

解释一下输入参数：f为目标函数系数向量，intcon为整数变量的索引，Aeq和beq为等式约束矩阵和右侧向量，lb为变量下界，[]表示没有不等式约束。

最后输出结果：

B = ceil(x(1));
X1 = ceil(x(2));
X3 = ceil(x(3));
profit = -fval;

解释一下输出结果：B为最优解的原木板数量，X1和X3为P1和P3分别在一块木板上切割的数量，profit为最大利润。

完整代码如下：

f = [0, 0, 0, -1];
Aeq = [373, 0, 406, 0; 0, 0, 0, 1; 0, 0, 0, 1; 0, 1, 0, 0];
beq = [0; 774; 1623; 35];
lb = [0; 0; 0; 0];
intcon = [1; 2; 3; 4];

[x,fval] = intlinprog(f,intcon,[],[],Aeq,beq,lb);

B = ceil(x(1));
X1 = ceil(x(2));
X3 = ceil(x(3));
profit = -fval;

另一个约束条件的模型可以按照上述过程类似求解。