function fit=fitnessindividual计算个体适应度值individual input 个体fit output 适应度值城市坐标city_coordinate=1129806 282862;1131352 281783;1122973 280959;1126055 281648;1128584 280472;1134333 283

假设有 $n$ 个城市，$i$ 和 $j$ 分别表示第 $i$ 个和第 $j$ 个城市，$d_{i,j}$ 表示它们之间的距离，$c_i$ 表示第 $i$ 个城市的货物量，$x_i$ 表示第 $i$ 个城市是否被选中（$x_i=1$ 表示选中，$x_i=0$ 表示未选中），则适应度函数可以表示为：

$$ f(x)=\sum_{i=1}^{n}\sum_{j=1}^{n}d_{i,j}c_i x_i x_j + w\sum_{i=1}^{n}c_i x_i $$

其中第一项表示运输费用，第二项表示设施建设费用，$w$ 表示设施建设的固定费用，$x_i x_j$ 表示第 $i$ 个城市和第 $j$ 个城市都被选中时为 1，否则为 0。如果距离超过某个阈值，可以加上惩罚项。

在实际实现中，可以先计算出每个城市到已选中的城市的距离，然后根据货物量计算出每个城市的运输费用，并将所有城市的运输费用相加得到第一项。第二项可以直接计算，即将所有被选中的城市的货物量相加，再乘以一个固定费用 $w$。最后将两项相加即可得到适应度值。