MATLAB如何求解大量决策变量的多目标整数规划问题?

求问如何用MATLAB求解很多个决策变量的多目标整数规划问题?
如下图的目标函数

其中，集合I={1，…，9}，J={1，…，7}，决策变量为xj,yij
代入数值如下所示

决策变量为x1-7，y11-97，一共70个决策变量
全拆开写的话，约束条件至少有70条……
求问怎么能简化代码😭