数学建模——公交调度优化
数学建模——公交调度优化
本文通过建立利润阈值模型鉴定高平峰期,综合考虑公交线路资源配置与乘客候车时间,建立多目标优化模型,通过人工免疫算法算法对公交调度方案进行优化,通过建立梯度提升树模型预测客流量,从而预测“高峰”和“平峰”时期。
考虑到市民出行并不是“均匀”的,“高峰”需要表征指的是乘客较多的时段,为了方便处理,我们将一段时间ΔT\varDelta TΔT内的运营总收入与运营总成本之差不小于某一个给定的“阈值”的时段。而问题的关键在于如何去建立时间与运营收入与成本的关系以及阈值的选择,本组通过建立净利润阈值模型,引进交通指数概念作为阈值选择的依据,从而合理定义了“高峰”与“平峰”的划分依据。\par
净利润阈值模型
考虑到运营公交车的成本有诸多因素,包括一段时间内可变成本C~\widetilde{C}C
——耗油量、磨损因度等,以及固定成本C‾\overline{C}C——司机工资以及购车费用,这里我们考虑为了判别高峰期或平峰期,本组首先将发车时间划分为若干紧密区间,对每个划分区间进行观测其净利润是否达到阈值:
Ω=[T0,T1]∪[T1,T2]∪⋯∪[TN−1,TN]\varOmega =\left[ T_0,T_1 \right] \cup \left[ T_1,T_2 \right] \cup \cdots \cup \left[ T_{N-1},T_N \right] Ω=[T0,T1]∪[T1,T2]∪⋯∪[TN−1,TN]
其中划分区间长度为:ΔT=Tn−Tn−1\varDelta T=T_n-T_{n-1}ΔT=Tn−Tn−1
设第iii班次行驶时间记为:[ti,ti+ttotal][t_i ,t_i +t_{total}][ti,ti+ttotal],因而每一班次在抽样划分区域[Tn,Tn+1]\left[ T_{n},T_{n+1} \right][Tn,Tn+1]上消耗的费用为cnic_{n}^{i}cni:
cni={[min(Tn+1,ti+1)−max(Tn,ti)]⋅v⋅f Tn⩽ti+1 and Tn+1⩾ti 0 Tn>ti+1 or Tn+1<tic_{n}^{i}=\begin{cases} \left[ \min \left( T_{n+1},t_{i+1} \right) -\max \left( T_n,t_i \right) \right] \cdot v\cdot f\,\, \qquad T_n\leqslant t_{i+1}\,\,and\,\,T_{n+1}\geqslant t_i\,\, \\ 0 \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad \quad \,\, T_n>t_{i+1}\,\,or\,\, T_{n+1}<t_i\\ \end{cases}cni={[min(Tn+1,ti+1)−max(Tn,ti)]⋅v⋅fTn⩽ti+1andTn+1⩾ti0Tn>ti+1orTn+1<ti
其中,第iii班次的公交起始时刻为tit_iti:
ti=T0+(i−1)Δti−1t_i=T_0+\left( i-1 \right) \varDelta t_{i-1}ti=T0+(i−1)Δti−1
则抽样划分区域[Tn,Tn+1]\left[ T_{n},T_{n+1} \right][Tn,Tn+1]上消耗的可变成本总费用为:
ΔC~N=∑i=1McNi\varDelta \widetilde{C}_N=\sum_{i=1}^M{c_{N}^{i}}ΔC
N=i=1∑McNi
因此,抽样划分区域 [Tn,Tn+1][ T_{n},T_{n+1}][Tn,Tn+1]上消耗的总费用 CCC 为:
ΔC=ΔC~+ΔC‾\varDelta C=\varDelta \widetilde{C}+\varDelta \overline{C}ΔC=ΔC
+ΔC
考虑到收入的影响因素为乘客的到达率,设第jjj个站点乘客随时间的变化的到达率为rjr_jrj,假设rjr_jrj服从均匀分布(人/分钟),则抽样划分区域$ \left[ T_{n},T_{n+1} \right]$上的总收入为:
ΔI=P⋅∑ΔtΔt⋅rj\varDelta I=P\cdot \sum_{\varDelta t}{\varDelta t \cdot r_j } ΔI=P⋅Δt∑Δt⋅rj
问题2的模型建立
假设乘客到达率服从均匀分布即:
rj∼U(tj−1,tj)r_j \sim U\left( t_{j-1},t_j \right) rj∼U(tj−1,tj)
乘客的平均等车时间为:
tˉ=ti−ti−12\bar{t}=\frac{t_i-t_{i-1}}{2} tˉ=2ti−ti−1
设第iii辆车发车时间为TiT_iTi,第iii辆车到达(或离开)第jjj站的时刻为:TijT_i^jTij则:
Tij={Ti+∑n=2jdn j>1Ti j=1T_{i}^j =\begin{cases} T_i+\sum_{n=2}^j{d_n}\,\,\quad \qquad j>1 \\ T_i\,\, \qquad \qquad \quad \qquad j=1\\ \end{cases}Tij={Ti+∑n=2jdnj>1Tij=1
考虑到顾客的时间利益、服务满意程度,引进候车超时率WTWTWT表现客户端的方便程度:
WT=候车时间超过15分钟的总人数调度期内上车总人数WT=\frac{\text{候车时间超过}15\text{分钟的总人数}}{\text{调度期内上车总人数}}WT=调度期内上车总人数候车时间超过15分钟的总人数
即为:
WT=∑Tij∑h=1hij{Wij(h)∣WTij(h)⩾15}∑TijPi(Tij)WT=\frac{\sum\limits_{T_{i}^j}{\sum\limits_{h=1}^{h_{ij}}{\left\{ W_{i}^j \left( h \right) |WT_{ij}\left( h \right) \geqslant 15 \right\}}}}{\sum\limits_{T_{i}^j}{P_i\left( T_{i}^j \right)}} WT=Tij∑Pi(Tij)Tij∑h=1∑hij{Wij(h)∣WTij(h)⩾15}
式中Wij(h)W_{i}^j (h)Wij(h)表现第iii辆车到达第jjj站时,该站上已经等待过hhh辆车,却仍未上车的人数,设这些人已经等待时间为:
WTij(h)=Tij−Ti−hj,Wij(0)=UijWT_{i}^j\left( h \right) =T_{i}^j-T_{i-h}^j,W_{i}^j \left( 0 \right) =U_{ij}WTij(h)=Tij−Ti−hj,Wij(0)=Uij
式中UijU_{i}^jUij表示在Ti−1jT_{i-1}^jTi−1j到TijT_{i}^jTij时间段内,新来到车站等车的人数:
WWij(0)=∫Ti−1jTijui(t)dtWW_{i}^{j}\left( 0 \right) =\int_{T_{i-1}^{j}}^{T_{i}^{j}}{u_i\left( t \right) dt}WWij(0)=∫Ti−1jTijui(t)dt
设hijh_{ij}hij表示第kkk辆车到达jjj站时,该站上等待最久乘客的候车辆数,则第iii辆车离开第jjj站时车上的乘客数可表示为:
Pi(Tij)={aij+∑r=0hijWij(r),hij∗=0,M,hij∗>0.P_i\left( T_{i}^{j} \right) =\begin{cases} a_{ij}+\sum_{r=0}^{h_{ij}}{W_{i}^{j}}\left( r \right) ,& \qquad h_{ij}^{*}=0,\\ M,& \qquad h_{ij}^{*}>0.\\ \end{cases}Pi(Tij)={aij+∑r=0hijWij(r),M,hij∗=0,hij∗>0.
其中aija_{ij}aij为第iii辆车到达第jjj站后,等乘客下车后,车上的剩余乘客数:
aij=max{(Pk(Tij−1)−Dj),0}a_{ij}=\max \left\{ \left( P_k\left( T_{i}^{j-1} \right) -D_j \right) \text{,}0 \right\}aij=max{(Pk(Tij−1)−Dj),0}
式中DjD_jDj为在第jjj站下车的乘客数。
设MMM为车辆最大承载人数,则此时可容纳的上车人数记为:
bij=M−aijb_{ij}=M-a_{ij}bij=M−aij
基于先到先上车的排队原则,在第k+1k+1k+1辆车到达第jjj站站时,仍然要等候车辆的最大等待次数hij∗h_{ij}^{*}hij∗:
hij∗=max{h∣∑r=0hijWij(r)⩽bij},h_{ij}^{*}=\max \left\{ h|\sum_{r=0}^{h_{ij}}{W_{ij}\left( r \right)}\leqslant b_{ij} \right\} , hij∗=max⎩
⎨
⎧h∣r=0∑hijWij(r)⩽bij⎭
⎬
⎫,
更多推荐
所有评论(0)