免疫优化算法选址 1、免疫算法是一种具有生成+检测 (generate and test)的迭代过程的群智能搜索算法。 从理论上分析，迭代过程中，在保留上一代最佳个体的前提下，免疫算法是全局收敛的 2、结合选址理论方法，建立模型，得出最优物流中心选址方案。 3、程序内容丰富，输出选址方案、适应度函数图像，效果可见。 4、可修改坐标配送量等参数，方便学习和推广
物流中心的选址直接影响着企业的运营成本。传统的选址方法在应对多约束条件时常常力不从心，这时候咱们就需要搬出免疫优化算法这个黑科技了。
免疫算法的核心在于模拟生物免疫系统：通过抗体克隆、变异、选择等机制，在迭代中不断逼近最优解。咱们先看一个典型的适应度函数设计（以运输成本最小化为目标）：
def calculate_cost(location, facilities, demand_points): total_cost = 0 for dp in demand_points: min_dist = float('inf') for fc in facilities: # 计算欧式距离 dist = ((location[fc][0]-dp[0])**2 + (location[fc][1]-dp[1])**2)**0.5 if dist < min_dist: min_dist = dist total_cost += min_dist * dp[2] # 距离×配送量 return total_cost这个函数暗藏玄机——通过遍历所有需求点，找到最近的物流中心计算运输成本。代码中的dp[2]代表配送量权重，后续修改参数时这里就是关键入口。
接下来是免疫算法的核心架构。咱们用Python实现一个简化版：
class Antibody: def __init__(self, solution): self.solution = solution # 解向量，例如候选坐标集 self.affinity = 0 # 亲和度（适应度） self.concentration = 0 # 浓度控制 def immune_optimization(max_iter=100): population = initialize_population() # 随机生成初始抗体群 best_solution = None for epoch in range(max_iter): # 评估种群 for ab in population: ab.affinity = 1 / calculate_cost(ab.solution) # 成本越低适应度越高 # 克隆扩增（精英保留） sorted_pop = sorted(population, key=lambda x: x.affinity, reverse=True) elites = sorted_pop[:int(0.2*len(population))] clones = generate_clones(elites) # 按适应度比例克隆 # 超变异操作 mutated_clones = [] for clone in clones: if random.random() < 0.15: # 在坐标维度上施加随机扰动 mutated = clone.solution + np.random.normal(0, 0.1, clone.solution.shape) mutated_clones.append(Antibody(mutated)) # 更新种群 population = elites + mutated_clones[:len(population)-len(elites)] # 记录最优解 current_best = max(population, key=lambda x: x.affinity) if best_solution is None or current_best.affinity > best_solution.affinity: best_solution = deepcopy(current_best) return best_solution.solution这段代码有几个精妙之处：克隆扩增时采用适应度比例复制，既保证了优秀个体的保留，又避免早熟收敛。变异操作在坐标维度添加高斯噪声，这种局部扰动正是跳出局部最优的关键。
实际运行后，咱们可以通过Matplotlib绘制适应度变化曲线：
plt.plot(history['iteration'], history['best_affinity']) plt.xlabel('Iteration') plt.ylabel('Transport Cost (1/affinity)') plt.title('Optimization Process') plt.show()曲线呈现典型的收敛形态——前期快速下降，后期趋于平稳。如果曲线出现剧烈震荡，可能需要调整变异概率或种群规模。
参数修改入口设计示例：
# 需求点坐标及配送量（可在此处修改） demand_points = [ (12, 25, 150), # (x, y, 配送量) (30, 18, 200), (5, 40, 180) ] # 算法参数配置 config = { 'population_size': 50, 'mutation_rate': 0.15, 'max_iterations': 100 }通过修改demand_points的数据，可以快速验证不同配送场景下的选址方案。比如把某个点的配送量调高，观察最优选址是否向该点偏移——这正是权重参数在空间上的直观体现。
最终的选址方案输出建议包含帕累托前沿分析，用散点图展示多个候选解的运输成本与建设成本的权衡关系。这种可视化能让决策者直观看到不同预算约束下的最优选择。
