for i in range(batch_size): # 更新需求U if i == 0: U = ui[0, :].copy() else: U += ui[i, :] # 模型求解（需要实现对应函数） Zjk, Wk, G = lastmile_model(phi_jk, tk, t_jk, U, beta1, beta2, beta3, K, J, c) # 路线排序 S = route_rank2(Wk, tk, Zjk, phi_jk, btw, t, headway) Temp_S = S.copy() # 车辆调度（需要实现对应函数） (btw, S, U, total_trip, total_traveltime, total_waitingtime, totoal_ridingtime, btw_vehicle, chengke) = vehicle_dispatch( btw, t_jk, S, U, Zjk, t, headway, total_trip, total_traveltime, total_waitingtime, totoal_ridingtime, btw_vehicle, chengke ) class VehicleRecord: def __init__(self): self.record = [] # 存储时间记录的二维列表 def vehicle_dispatch(btw, t_jk, S, U, Zjk, t, headway, total_trip, total_traveltime, total_waitingtime, totoal_ridingtime, btw_vehicle, chengke): U = U.astype(np.int32).copy() # 第一部分：更新车辆可用时间 for i in range(len(btw)): if btw[i] < t: vr = btw_vehicle[i] if len(vr.record) == 0: # 处理空记录情况 vr.record.append([t, t, 0, 0]) else: last_end = vr.record[-1][1] vr.record.append([last_end, t, 0, 0]) btw[i] = t # 第二部分：车辆调度循环 for current_time in range(t, t + headway + 1): # 找到可用车辆 available = [i for i, bt in enumerate(btw) if current_time > bt] sorted_available = sorted(available, key=lambda x: btw[x]) if sorted_available and np.sum(U) > 0: for bus_idx in sorted_available: if np.sum(U) <= 0: break # 更新统计量 total_trip[0] += 1 route_info = S[0] total_traveltime += route_info[2] # 乘客服务处理 served_pax = Zjk[:, route_info[0]-1] # 调整索引 totoal_ridingtime += np.sum(served_pax * t_jk[:, route_info[0]-1]) waiting_time = btw[bus_idx] - t total_waitingtime += np.sum(served_pax) * waiting_time # 更新乘客信息 stop_indices = np.where(served_pax > 0)[0] for stop in stop_indices: pax_mask = (chengke[:, 2] == stop+1) & (chengke[:, 9] == 0) pax_candidates = np.where(pax_mask)[0] if len(pax_candidates) > 0: num_pax = min(served_pax[stop], len(pax_candidates)) selected = pax_candidates[:num_pax] chengke[selected, 9] = 1 chengke[selected, 4] = btw[bus_idx] chengke[selected, 5] = chengke[selected, 4] - chengke[selected, 3] chengke[selected, 6] = t_jk[stop, route_info[0]-1] chengke[selected, 7] = route_info[0] chengke[selected, 8] = bus_idx + 1 # 车辆编号从1开始 # 更新车辆记录 vr = btw_vehicle[bus_idx] if not vr['record']: vr['record'].append([btw[bus_idx], btw[bus_idx] + route_info[2], route_info[0], route_info[2]]) else: last_end = vr['record'][-1][1] vr['record'].append([last_end, last_end + route_info[2], route_info[0], route_info[2]]) # 更新系统状态 btw[bus_idx] += route_info[2] U -= Zjk[:, route_info[0]-1].T del S[0] if np.sum(U) <= 0: break # 处理未满足的需求 if current_time == t + headway and np.sum(U) > 0: total_waitingtime += np.sum(U) * headway ValueError: cannot delete array elements