2010年认证杯SPSSPRO杯数学建模B题(第一阶段)交通拥堵问题全过程文档及程序

2010年认证杯SPSSPRO杯数学建模

交通拥堵问题

B题 Braess 悖论

原题再现：

  Dietrich Braess 在 1968 年的一篇文章中提出了道路交通体系当中的Braess 悖论。它的含义是：有时在一个交通网络上增加一条路段，或者提高某个路段的局部通行能力，反而使所有出行者的出行时间都增加了，这种为了改善通行能力的投入不但没有减少交通延误，反而降低了整个交通网络的服务水平。人们对这个问题做过许多研究，在城市建设当中也尽量避免这种现象的发生。但在复杂的城市道路当中，Braess 悖论仍然不时出现，造成实际交通效率的显著下降。在此，请你通过合理的模型来研究和解决城市交通中的 Braess 悖论。
  第一阶段问题：
  (1) 通过分析实际城市的道路交通情况（自行查询的数据需给出引用来 1源），建立合理的模型，判断在北京市二环路以内的路网中（包括二环路）出现的交通拥堵，是否来源于 Braess 悖论所描述的情况。
  (2) 请你建立模型以分析：如果司机广泛使用可以反映当前交通拥堵情况的GPS 导航系统，是否会缓解交通堵塞，并请估计其效果。

整体求解过程概述(摘要)

  本文首先通过对非合作网络中 Nash 平衡点与 Pareto 边界理论的分析，通过单起点单终点的简单路网入手, 分析 Braess 悖论的成因，Braess 悖论实质上是非合作网络中 Nash 平衡点不满足 Pareto 最优性时出现的现象。定量的给出出行时间与流量的关系，这里引入了路段延迟参数的概念，量化得到延迟参数与路段流量，出行时间，出行密度之间的关系。进而通过建立北京二环区交通路网模型，拟合出 flow ~ / L speed 关系曲线，得到各路段的延迟参数，进而得到单月内北京二环区 Pareto 最优解的判别条件。
  其次，根据 Pareto 边界理论，我们以年月为单位统计北京二环以内交通流量因素，得到城区交通年延迟参数，并计算了共 120 个月的延迟参数。接着，我们构建了 Nash 均衡原理模型，引入了路段 Nash 平衡比率的概念，即当某一时间范围内，路段 Nash 平衡比率低于路段延迟参数时，则发生了 Braess 悖论，当同时又满足 Nash平衡波动率超过 30%时，为严重 Braess 悖论，从而完善了 Nash 平衡比率理论方程。
  再次，根据已建立的 Nash 平衡比率理论方程，在此基础上构建小波神经网络来对北京二环区交通路网路段 Nash 平衡比率进行预测，这里我们将选择路段延迟参数为样本，根据路段延迟参数与路段 Nash平衡比率的动态关系式，对路段 Nash 平衡比率进行计算。得到结果:东城区，朝阳区路段常年持续存在 BRAESS悖论情况（平均每年超过 6 个月）。尽管我们发现坛路拓宽，普方路拓宽，永定门桥延修，天坛路改道会给北京二环区南部路段带来交通改善，但结果显示，整个北京二环路段的交通情况的平衡性变得更不稳定，BRAESS 悖论情况更加突出。譬如在 2006 年，2009 年的东城区，朝阳区路段堵塞问题更加尖锐，这也使得本来相对通畅的南段也出现 BRAESS 悖论情况。
  接着，我们建立了 GPS 动态导航的车辆运行模型，量化得到了时耗期望值和时耗方差两个目标函数并作为 ATM 路由选择的约束条件，以降低 Nash 平衡比率为目标，通过ATM 路由选择进行 GPS 导航对路段影响的预测。当我们控制业务源提高 GPS 使用率时，路段 NASH 平衡比率也持续上升，使用率升至 0.34 时，NASH 平衡比率上升幅度最大化。
  最后，我们通过灰色 Gompertz 模型对 ATM 路由选择算法进行误差分析，检验结果表明此抽检方案为合理方案，误差范围 5.84%。通过遗传算法对路段 Nash 平衡比率进行优化，得到最优化后的结果，证明了ATM路由选择算法对GPS动态导航研究的可靠性。

问题分析

  首先分析单一路段在单一时间范围内的 BRAESS 悖论情况，然后以北京二环以内的综合路段为考虑对象，天坛路拓宽，普方路拓宽，永定门桥延修，天坛路改道的 NASH 平衡比率，并通过 MATLAB 工具 CHIP（）提取出月波动率超过 30%的比率数进行标定，同时，我们将近十年的北京二环以内主要路段的 BRAESS 悖论情况进行分析，并讨论 GPS 动态导航对路段BRAESS 悖论现象的影响。

模型的假设：

  （1）天气因素对 BRAESS 悖论的影响忽略不计
  （2）《北京交通发展年报》及《全市综合交通调查》提供了精确无误的数据
  （3）北京市区内地铁路线改造的影响忽略不计。

论文缩略图：

全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可

部分程序代码：(代码和文档not free)

% III 为预测点 
% X 为预测序列（一行） 
% X0X 为预测值（对有 III 而言） 
% Q0Q 为预测精度 
% PP 为小误差概率 
% C 为后验误差比 
% PP>0.95 或 C<0.35 为 1 级(好) 
% PP>0.80 或 C<0.50 为 2 级(合格) 
% PP>0.70 或 C<0.65 为 3 级(勉强合格) 
% PP<=0.70 或 C>=0.65 为 4 级(不合格) 
% 
% AUA=[A;U;UA;XUA]; 
% XEG=[XG;EG]; 
% PPC=[PP;C]; 
FUNCTION [X0X,Q0Q,PPC,PDDSTR,XII,AUA,XEG]=GM11(X,III); 
%%%%%%%%%%%%%%%%%%%%%%% 
 
%例子 
%X=[3.711 3.723 3.716 3.721 3.728] 
%X=[2.67 3.13 3.25 3.36 3.56 3.72] 
% X=[2.97 3.23 3.29 3.46 3.59 3.71] 
%X=[43.45 47.05 52.75 57.14 62.64 68.52] 
% X=[3.38 4.27 4.55 4.69 5.59] 
% X=[4.24 4.33 5.20 6.42 7.32 8.53 8.82 10.72] 
% 
% X=[3.38 4.27 4.55 4.69 5.59] 
% III=[ 5 6 7 8 9 10 11] 
% III=[] 
 
%%%%%%%%%%%%%%%%%%%%%%% 
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% %% 
%% 一、数据处理 %% 
%% %% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
 
%1 建立 X1 生成数列 
%FORMAT LONG 
A12=SIZE(X); 
N=A12(2); 
X1=ONES(1,N); 
X1(1)=X(1); 
FOR I=2:N 
 X1(I)=X1(I-1)+X(I); 
END 
X1; 
%2 构造数据距阵 B 和数据向量 YN: 
B=ONES(N-1,2); 
FOR I=1:N-1 
 B(I,1)=(-0.5)*(X1(I)+X1(I+1)); 
END 
YN=ONES(N-1,1); 
FOR I=2:N 
 YN(I-1,1)=X(I); 
END 
B; 
YN; 
%3 计算 B'B,INV(B'B)和 B'YN；得到参数 A 和 U 
%BTB=B'*B;% 
%BTB=INV(BTB); 
%BTYN=B'*YN; 
%AU=BTB*BTYN; 
AU=INV(B'*B)*B'*YN; 
 
 
A=AU(1,1); 
U=AU(2,1); 
 
% FPRINTF('参数 A：%F\N',A); 
% FPRINTF('参数 U：%F\N',U); 
UA=U/A; 
XUA=X(1)-UA; 
%4 得出预测模型 
%预测第 I 的预测值 
%有问题 
 
% XI1=XUA*EXP((-A)*I)+UA;% 
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%% %% 
%% 二、模型检验 %% 
%% %% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%5:残差检验 
%51 计算 X1 
X1=ONES(1,N); 
FOR I=1:N 
 X1(I)=XUA*EXP((-A)*(I-1))+UA; 
END 
X1; 
%52 累减生成 X0 序列 
X0=ONES(1,N); 
X0(1)=X1(1); 
FOR I=2:N 
 X0(I)=X1(I)-X1(I-1); 
END 
X0; 
%53 计算绝对误差及相对误差序列 
DD=ABS(X0-X);%计算绝对误差序列 
DDD=(DD./X).*100;%计算相对误差序列 X% 
% FPRINTF('相对误差序列：'); 
% FPRINTF('%F%% ',DDD); 
% FPRINTF('\N'); 
%6 进行关联度检验 
%61 计算序列绝对误差 XX 
XX=ONES(1,N); 
FOR I=1:N 
 XX(I)=X(I)-X0(I); 
END 
% FPRINTF('序列绝对误差：'); 
% FPRINTF('%F ',XX); 
% FPRINTF('\N'); 
 
MINXX=MIN(DD); 
% FPRINTF('最小差：'); 
% FPRINTF('%F ',MINXX); 
% FPRINTF('\N'); 
 
MAXXX=MAX(DD); 
% FPRINTF('最大差：'); 
% FPRINTF('%F ',MAXXX); 
% FPRINTF('\N'); 
%62 计算关联系数 P=0.5 (则 R=0.6) 
P=0.5; 
NI=ONES(1,N);
FOR I=1:N 
 NI(I)=(MINXX+P*MAXXX)/(DD(I)+P*MAXXX); 
END 
%63 计算关联度,(P=0.5 时，则 R=0.6) 
R=(1/N)*SUM(NI); 
% FPRINTF('请查找 P=0.5 是的检验准则 R=0.6 是否大于%F\N',R); 
%7 后验差检验 
%71 
XM=MEAN(X); 
%72 求的均方差 
S1=(SUM((X-XM).^2)/(N-1))^(1/2); 
%73 计算残差的均值 
DDM=MEAN(DD); 
%74 计算残差的均方差 
S2=(SUM((DD-DDM).^2)/(N-1))^(1/2); 
%75 计算后验误差比 C: 
C=S2/S1; 
% FPRINTF('验误差比 C：%F\N',C); 
%76 计算小误差概率 
%%%%%%%%%%%%%% 
 
PR=ABS(XX-MEAN(XX))<0.6745*S1; %M 满足条件的样本 
PA=SIZE(FIND(PR==1)); 
PS=PA(1,2); %M 满足条件的样本个数 
PB=SIZE(PR); 
PS=PB(1,2); %M 总样本个数 
 
PP=PS/PS; %小误差概率 
% FPRINTF('小误差概率 PP：%F\N',PP); 
 
%%%%%%%%%%%%%%%% 
%检验 预测精度 YCD 
 
% FPRINTF(' PP>0.95 或 C<0.35 为 1 级(好)\N PP>0.80 或 C<0.50 
为 2 级(合格)\N PP>0.70 或 C<0.65 为 3 级(勉强合格)\N PP<=0.70 
或 C>=0.65 为 4 级(不合格)\N\N',C); 
IF PP>0.95 
 PD=1; 
% FPRINTF(' 因 PP>0.95 且 '); 
END 
IF PP<=0.95&PP>0.80 
 PD=2; 
% FPRINTF(' 因 PP<=0.95&PP>0.80 且 '); 
END
IF PP<=0.80&PP>0.70 
 PD=3; 
% FPRINTF(' 因 PP<=0.80&PP>0.70 且 '); 
END 
IF PP<=0.70 
 PD=4; 
% FPRINTF(' 因 PP<=0.70 且 '); 
END 
 
PD1=0; 
IF C<0.35 
 PD1=1; 
% FPRINTF('C<0.35'); 
END 
IF C>=0.35&C<0.50 
 PD1=2; 
% FPRINTF('C>=0.35&C<0.50'); 
END 
IF C>=0.50&C<0.65 
 PD1=3; 
% FPRINTF('C>=0.50&C<0.65'); 
END 
IF C>=0.65 
 PD1=4; 
% FPRINTF('C>=0.65'); 
END 
 
PDD=MAX(PD,PD1); 
 
IF PDD==1 
 PDDSTR='1 级（好）'; 
% FPRINTF(' 故根据经验，预测精度为 1 级（好）\N\N'); 
END 
IF PDD==2 
 PDDSTR='2 级（合格）' ; 
% FPRINTF(' 故根据经验，预测精度为 2 级（合格）\N\N'); 
END 
IF PDD==3 
 PDDSTR='3 级（勉强合格）'; 
% FPRINTF(' 故根据经验，预测精度为 3 级（勉强合格）\N\N'); 
END 
IF PDD==4 
 PDDSTR='4 级（不合格）'; 
% FPRINTF(' 故根据经验，预测精度为 4 级（不合格）\N\N');

全部论文及程序请见下方“ 只会建模 QQ名片” 点击QQ名片即可