cons3 = sum(sum (t_dock * x_ship_dock , 2), 3) <= c_dock; 这句代码还有错误，请通过其他方式进行修改

如果上述代码仍然出现错误，可以尝试使用 bsxfun 函数进行计算。bsxfun 函数可以对两个大小不同的矩阵进行逐元素计算，需要对其中一个矩阵进行扩展操作，扩展成与另一个矩阵相同的大小。具体而言，根据 t_dock 和 x_ship_dock 的维度，可以将 t_dock 扩展为一个 N_ship × N_dock × T_dock 的三维数组，然后再与 x_ship_dock 逐元素相乘，得到一个 N_ship × N_dock × T_dock 的三维数组，表示每种船只在每个码头的停靠时间。然后，使用 sum 函数对第二个和第三个维度进行求和，得到每个码头的总停靠时间。最后，使用 sum 函数对所有码头的总停靠时间进行求和，并与码头容量 c_dock 进行比较。修改后的代码如下：

cons3 = sum(sum(sum(bsxfun(@times, t_dock, x_ship_dock), 2), 3), 1) <= c_dock;

在这里，bsxfun 函数对 t_dock 和 x_ship_dock 逐元素相乘，得到一个 N_ship × N_dock × T_dock 的三维数组，然后使用 sum 函数对第二个和第三个维度进行求和，得到每个码头的总停靠时间，最后使用 sum 函数对所有码头的总停靠时间进行求和，并与码头容量 c_dock 进行比较。