matlab求助帖！求大神告诉我这个用什么算法能最方便的解决出来

由于班轮为定期航班，临时改变船舶数量是不现实的，所以船舶数量也是既定的，所以在这里假定该航线只有五条4000TEU的船营运。
目标函数：(由于四个季节内部运输情况一致，所以我们仅取四个季节各一个航次求单舱日均利润)

变量及参数说明：
P:均天单舱利润
TC:船舶最大载箱量
G,g:分别表示按季节因素分的时间段集合和在这个集合中的顺序
R:总收益
Cp:靠港所用成本
Cs:航线运行成本
T总：表示船舶航行的总时间
rij:i港到j港的运价
xijg:g季节重箱从i港到j港的运量
Bij:航段(i,j)上是否有货物运输，以0-1变量表示，即可及航段（B是A的可达矩阵）
Aij:船舶是否在航段(i,j)上航行，以0-1变量表示
Cpe,Cpf:分别表示空箱、重箱的单箱装卸费用
Cpi:表示i港的港口使费
约束条件及说明：
航线数量约束
保证研究过程中只有一条航线，即令所有航段都为既定航线上的航段。

环绕运输
船舶应保证其运输轨迹的连贯性，即从哪个港驶入就要从哪个港驶出。

可及航段约束
在可及航段和航段之间建立起一定的关系，方便问题的求解。

服务水平约束

该约束要求按照航线方案下的可及航段运输时间不应大于航段两端港直达时间的一定倍数，以此来保证服务质量。

航段货物容量约束
由于船舶的舱容有限，所以在每个航段营运过程中，其载箱量都不应超过其最大载箱量。

相关变量及参数的说明：
N：表示港口所有节点的集合
S：表示始发港所有节点的集合
L：表示所有航段的集合
Xij：表示从i港到j港重箱的舱位分配量
Yij：表示从i港到j港空箱的舱位分配量
M：表示极大的正数值
λ：表示服务约束中的最大容忍倍数
Xm：表示经过m航段的某可及航段的舱位分配
Zi：表示在航线方案中是否挂靠i港口，以0-1表示
船舶实际作业量说明：
由于船舶舱容以及计划舱位分配量的约束，所以船舶在i、j港之间的实际重箱运输量应当遵循以下原则：

这里xijg表示实际运输量，Xij表示计划舱位分配量，Qijg表示不同季节货物运输的需求量。
航线运行成本及班期计算
本文航线设计中的船舶英语方式为租船营运，则航线运行成本分为租金成本、重油成本、轻油成本。

参数说明：
CS：表示一个周期内船舶租金的费用
cs：表示船舶的日租金
CF：表示一个周期内船舶的燃料
FC：表示重油单位时间的消耗量
FCP：表示重油价格
DC：表示轻油单位时间的消耗量
DCP：表示轻油价格
T：表示船舶绕挂靠港口航行时间
T总：表示船舶航行的总时间
Sl：表示船舶绕挂靠港口航行一周的总距离
Tpi：表示船舶在i港靠泊的平均时间
T班：班轮的班期