✨ 本团队擅长数据搜集与处理、建模仿真、程序设计、仿真代码、EI、SCI写作与指导，毕业论文、期刊论文经验交流。
✅ 专业定制毕设、代码
✅如需沟通交流，点击《获取方式》

（1）雨天改进元胞传输模型参数标定与验证：

在传统CTM中引入降雨强度R（mm/h）与能见度V（m）作为影响因素，对自由流速度v_f进行衰减修正，定义v_f(R,V)=v_f0·(1-α·R^0.6)·(V/V_ref)^0.5，其中α=0.012，V_ref=300m。临界密度k_c和拥堵密度k_j也随路面附着力系数μ下降而线性压缩，μ与降雨量关系采用实测回归公式μ=0.9-0.004R。基于上海市某高速公路雨天交通流数据，抽取12场降雨事件的微波车辆检测器记录，采用半自动标定法，通过遗传算法最小化仿真流量与实际流量误差，标定出模型改参数β及波速w，最终仿真流量误差MARE为7.6%。改进后的雨天CTM能够较准确地再现降雨时通行能力折减现象，自由流速度下降幅度在暴雨（R>30mm/h）时可达28%，与实地观测的27.5%基本吻合，模型通过K-S检验，置信度95%，证明适用于可变限速策略仿真。

（2）基于模型预测控制的面向通行效率可变限速算法：

以改进的雨天CTM为预测模型，建立模型预测控制框架，目标函数最小化总行程时间TTS与限速值变化惩罚项之和，控制变量为各路段可变限速值VSL_i(k)，约束包括限速值在40km/h至120km/h之间，相邻控制周期内限速变化不超过20km/h，且下游限速不得高于上游限速。预测时域取10个时间步长（每步20s），控制时域为3步。采用遗传算法在约束条件下对VSL序列进行滚动优化，种群规模60，每步迭代40代。在模拟暴雨场景下，无控制时路段平均速度下降至52km/h，并出现走走停停的波传播，上游通过流量下降23.5%；启用MPC-VSL后，瓶颈路段速度提升至64km/h，总延误时间减少18.4%，平均停车次数降低31%。通过多次随机降雨场景仿真，MPC-VSL能将交通崩溃临界密度延迟30%，显现出良好的主动管控能力。

（3）面向安全的事故风险Logistic模型与协同限速算法：

利用同一高速公路三年交通事故数据和对应的降水、交通流数据，建立Logistic回归事故概率模型，自变量包括车道平均速度v_avg、速度标准差σ_v、降雨强度R及小时流量q。模型输出事故风险指数P_crash，在P_crash>0.35时视为高风险状态。基于此设计安全导向的可变限速策略Safety-VSL，当P_crash超过阈值时，主线下游限速值按风险指数梯度降低，同时配合匝道信号灯控制入口流量。用VISSIM和MATLAB联合仿真平台验证，暴雨R=32mm/h场景中，无控制时事故概率峰值0.62，安全VSL控制后可控制在0.38以下；同时，在安全干预时段内，速度标准差降低26%，事故仿真计数从每百万车公里1.8起降至0.6起。两种可变限速算法通过集成的C#调度程序在VISSIM COM接口中执行，实现了宏观与微观耦合仿真，结果证明面向安全与效率的可变限速策略均可显著改善雨天高速公路运行质量。

import numpy as np
from scipy.optimize import differential_evolution

# 雨天改进CTM模型
class RainyCTM:
    def __init__(self, cells, params):
        self.N = cells
        self.vf0 = params['vf0']  # 120 km/h
        self.kj = params['kj']
        self.w = params['w']
        self.dx = params['dx']
        self.dt = params['dt']
        self.density = np.zeros(self.N)

    def update_rain(self, R, V):
        alpha = 0.012
        vf = self.vf0 * (1 - alpha * R**0.6) * (V/300)**0.5
        mu = 0.9 - 0.004*R
        kj_rain = self.kj * mu
        return vf, kj_rain

    def step(self, inflow, R, V):
        vf, kj = self.update_rain(R, V)
        flow = np.zeros(self.N+1)
        for i in range(1, self.N):
            q_cap = vf * self.density[i]
            if self.density[i] <= kj/2:
                supply = q_cap
            else:
                supply = self.w * (kj - self.density[i])
            demand = min(vf * self.density[i-1], supply)
            flow[i] = demand
        self.density += (self.dt/self.dx) * (flow[:-1] - flow[1:])
        self.density[0] += (self.dt/self.dx) * inflow
        return self.density.copy()

# MPC可变限速优化
def mpc_vsl_optimization(ctm, R, V, horizon=10):
    def cost(vsl_seq):
        # vsl_seq每个元素为各路段时间步限速值
        total_tts = 0
        sim = RainyCTM(ctm.N, ctm.__dict__)
        for k in range(horizon):
            vf_adj = vsl_seq[k] / ctm.vf0 * ctm.vf0
            sim.vf0 = vf_adj
            dens = sim.step(500, R, V)
            total_tts += np.sum(dens * ctm.dx) * ctm.dt / 3600
        penalty = 0.1 * np.sum(np.diff(vsl_seq)**2)
        return total_tts + penalty
    bounds = [(40, 120)] * horizon
    result = differential_evolution(cost, bounds, maxiter=40, popsize=60)
    return result.x

# 事故概率Logistic模型
def crash_prob_logistic(speed_avg, speed_std, R, flow):
    z = -4.3 + 0.08*speed_avg + 0.12*speed_std + 0.033*R + 0.002*flow
    p = 1 / (1 + np.exp(-z))
    return p

# 安全限速生成
def safety_vsl(p_crash, current_limit):
    if p_crash > 0.35:
        reduction = 20 * (p_crash - 0.35)
        new_limit = max(40, current_limit - reduction)
    else:
        new_limit = current_limit
    return new_limit

# 示例
if __name__ == '__main__':
    params = dict(vf0=120, kj=100, w=25, dx=0.2, dt=10/3600.0)
    ctm = RainyCTM(50, params)
    R, V = 25, 250   # 中雨，能见度250m
    vsl_opt = mpc_vsl_optimization(ctm, R, V, horizon=5)
    print('优化的限速序列(km/h):', vsl_opt)
    p = crash_prob_logistic(60, 15, R, 1200)
    print('事故概率:', p)
    new_vsl = safety_vsl(p, 100)
    print('安全限速:', new_vsl)