发明名称 一种基于时空属性的交通网络节点重要度测算方法
摘要 本发明提供了一种基于时空属性的交通网络节点重要度测算方法,首先以道路交通路网的单元路段为节点,以组成构件上传感器采集的检测数据组成时间序列,基于时间序列的相关性对路网进行建模;随后,考虑节点之间的空间关系,对节点重要性进行测算,包括获取交通网络的邻接矩阵和转移矩阵;构造距离矩阵和概率矩阵,并标准化概率矩阵,构造Google矩阵G;用幂法求解矩阵G的主特征向量,即得到各节点重要度。本发明在模型建立和重要度测算时考虑了路网的时空属性,能更好地表征实际网络,测算结果也更符合实际情况。
申请公布号 CN106530687A 申请公布日期 2017.03.22
申请号 CN201610893842.3 申请日期 2016.10.13
申请人 北京交通大学 发明人 秦勇;杨艳芳;贾利民;董宏辉;张庆;张纪升;孙晓亮;李斌
分类号 G08G1/01(2006.01)I 主分类号 G08G1/01(2006.01)I
代理机构 北京德崇智捷知识产权代理有限公司 11467 代理人 王金双
主权项 一种节点重要度测算方法,特别是一种基于时空属性的交通网络节点重要度测算方法,其中,以单位路段为一个节点,其特征在于,包括如下步骤:1)建立基于时间序列的交通网络模型首先计算路网中任意两个路段的相关系数r(S<sub>i</sub>,S<sub>j</sub>):<maths num="0001"><math><![CDATA[<mrow><mi>r</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>,</mo><msub><mi>S</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>T</mi><mi>D</mi></mrow></msubsup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow></mrow><msqrt><mrow><msubsup><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>T</mi><mi>D</mi></mrow></msubsup><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msubsup><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>T</mi><mi>D</mi></mrow></msubsup><msup><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>-</mo><msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mfrac></mrow>]]></math><img file="FDA0001129808250000011.GIF" wi="1223" he="245" /></maths>其中,S<sub>i</sub>为第i条路段上的时间序列,S<sub>j</sub>为第j条路段上的时间序列,y<sub>ik</sub>为S<sub>i</sub>上时间序列第k个值,<img file="FDA0001129808250000012.GIF" wi="46" he="63" />为S<sub>i</sub>上时间序列的均值,y<sub>jk</sub>为S<sub>j</sub>上时间序列k个值,<img file="FDA0001129808250000013.GIF" wi="53" he="69" />为S<sub>j</sub>上时间序列的均值,T为观测的天数,D为一天内观测的间隔数;其次根据时间序列构造交通网络模型:TNW=(N,E,R)其中,N={n<sub>1</sub>,n<sub>2</sub>,....,n<sub>M</sub>}是交通网络中节点的有限集合,M为网络中节点的个数;E={e<sub>ij</sub>|i≠j,i,j∈{1,2,....,M}}是交通网络中的连接节点的边的有限集合;R={r<sub>ij</sub>|i≠j,i,j∈{1,2,....,M}}是交通路网中相关系数的有限集合。2)计算交通网络构件节点重要度21)构建交通网络的邻接矩阵A和转移矩阵A<sup>*</sup>,交通网络的邻接矩阵A和转移矩阵A<sup>*</sup>分别表示为:<img file="FDA0001129808250000014.GIF" wi="566" he="111" /><img file="FDA0001129808250000021.GIF" wi="590" he="102" />其中<img file="FDA0001129808250000022.GIF" wi="510" he="271" />22)构造距离矩阵D和概率矩阵K,并标准化矩阵K得到矩阵K<sub>N</sub>距离矩阵D表示为:<maths num="0002"><math><![CDATA[<mrow><mi>D</mi><mo>=</mo><msubsup><mrow><mo>(</mo><msub><mi>b</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><mo>,</mo><msub><mi>b</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mfrac><mn>1</mn><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mfrac></mtd><mtd><mrow><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mi>j</mi></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001129808250000023.GIF" wi="686" he="254" /></maths>其中,d<sub>ij</sub>为节点i和节点j在路网中最短路径的长度;概率矩阵K表示为:<maths num="0003"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mi>A</mi><mo>&CircleTimes;</mo><mi>D</mi><mo>+</mo><mi>&epsiv;</mi><mi>F</mi></mrow>]]></math><img file="FDA0001129808250000024.GIF" wi="389" he="55" /></maths>其中,<img file="FDA0001129808250000025.GIF" wi="150" he="55" />为邻接矩阵A和距离矩阵D的Hadamard乘积;<img file="FDA0001129808250000026.GIF" wi="165" he="135" />F为元素为1的M×M矩阵;设概率矩阵K=(k<sub>1</sub>,k<sub>2</sub>,...,k<sub>j</sub>,..,k<sub>M</sub>),其中k<sub>j</sub>为概率矩阵K的第j列向量,K通过如下公式进行标准化,得到矩阵K<sub>N</sub>:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>k</mi><mi>j</mi><mo>*</mo></msubsup><mo>=</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>k</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>k</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mfrac><msub><mi>k</mi><mrow><mn>2</mn><mi>j</mi></mrow></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>k</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mo>...</mo><mo>,</mo><mfrac><msub><mi>k</mi><mrow><mi>M</mi><mi>j</mi></mrow></msub><mrow><mo>|</mo><mo>|</mo><msub><mi>k</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mi>T</mi></msup><mo>;</mo></mrow>]]></math><img file="FDA0001129808250000027.GIF" wi="773" he="229" /></maths> 23)构造矩阵G=(1‑α)A<sup>*</sup>+αK<sub>N</sub>,α为阻尼因子;24)用幂法求解矩阵G的主特征根λ<sub>1</sub>=1的主特征向量X<sub>1</sub>:X<sub>1</sub>={x<sub>1</sub>(1),x<sub>1</sub>(2),....,x<sub>1</sub>(N)},X<sub>1</sub>则为交通网络中N个节点的重要度值。
地址 100044 北京市海淀区上园村3号