Processing math: 100%

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

信号交叉口饱和流率动态提取方法

杨东援 罗江邻 刘翀 段征宇

杨东援, 罗江邻, 刘翀, 段征宇. 信号交叉口饱和流率动态提取方法[J]. 交通运输工程学报, 2013, 13(1): 98-103. doi: 10.19818/j.cnki.1671-1637.2013.01.015
引用本文: 杨东援, 罗江邻, 刘翀, 段征宇. 信号交叉口饱和流率动态提取方法[J]. 交通运输工程学报, 2013, 13(1): 98-103. doi: 10.19818/j.cnki.1671-1637.2013.01.015
YANG Dong-yuan, LUO Jiang-lin, LIU Chong, DUAN Zheng-yu. Dynamic extraction method of saturation flow rate at signalized intersection[J]. Journal of Traffic and Transportation Engineering, 2013, 13(1): 98-103. doi: 10.19818/j.cnki.1671-1637.2013.01.015
Citation: YANG Dong-yuan, LUO Jiang-lin, LIU Chong, DUAN Zheng-yu. Dynamic extraction method of saturation flow rate at signalized intersection[J]. Journal of Traffic and Transportation Engineering, 2013, 13(1): 98-103. doi: 10.19818/j.cnki.1671-1637.2013.01.015

信号交叉口饱和流率动态提取方法

doi: 10.19818/j.cnki.1671-1637.2013.01.015
基金项目: 

国家863计划项目 2011AA110305

详细信息
    作者简介:

    杨东援(1953-), 男, 江苏南京人, 同济大学教授, 工学博士, 从事交通运输规划研究

  • 中图分类号: U491.51

Dynamic extraction method of saturation flow rate at signalized intersection

More Information
Article Text (Baidu Translation)
  • 摘要: 提出了利用感应线圈检测器动态提取饱和流率的方法, 前后车辆离开线圈的时间差为车头时距, 计算第4辆至最后一辆处于饱和状态车辆的平均饱和车头时距, 运用指数平滑法处理历史饱和车头时距与当前周期饱和车头时距。确定了车型及饱和车头时距判断阈值, 当线圈占用时间大于小型车平均占用线圈时间2倍时, 判断为大型车, 小型车的饱和车头时距判断阈值为历史平均值加1 s, 大型车的饱和车头时距判断阈值为历史平均值加5 s。用VISSIM软件进行仿真, 验证提取方法的有效性。仿真结果表明: 动态提取方法能减少饱和车头时距突变的影响, 当前周期车头时距骤减31.3%, 饱和流率仅增加5.6%, 5个周期的饱和流率分别为1 782、1 682、1 600、1 690、1 773 veh·h-1, 而HCM模型的计算结果为1 680 veh·h-1。与传统方法相比, 该方法能满足动态提取的需求, 实施成本低。

     

  • 城市道路交叉口是城市道路的咽喉, 是制约城市道路交通功能的瓶颈。影响交叉口通行能力至关重要的一个参数就是交叉口饱和流率。交叉口饱和流率是指在一次连续的绿灯信号时间内, 进口道上一列连续车队能通过进口道停止线的最大流量。饱和流率是计算交叉口通行能力, 评价交叉口服务水平以及进行信号配时设计等必不可少的输入参数。目前, 获取饱和流率的方法大致分为两类, 模型计算法与实地测量法。不少国家都建立了饱和流率计算模型, 如美国、英国、加拿大等, 其中以美国的道路通行能力手册(Highway Capacity Manual, HCM)模型运用最为广泛。HCM模型综合考虑车道宽度、大车比例、进口道坡度、停车情况、公交阻挡、区域类型、车道利用、车道功能以及行人和自行车等因素的影响, 在理想饱和流率的基础上乘以相应的校正系数得到饱和流率。实际运用中获取饱和流率最可靠的方法是直接通过现场测量来确定, 首先测量饱和车流的平均饱和车头时距, 然后根据其与流率的关系计算出饱和流率, 这种方法不需要根据研究地点的具体道路几何与交通环境等进行修正, 但无论是现场人工测量或视频录像的实地测量, 都需要耗费大量的人力物力, 如果要实现动态测量, 成本将更高。

    目前, 国内外对于饱和流率的研究多集中在饱和流率影响因素及计算模型方面。Khosla等利用5个交叉口的实际数据研究了长绿灯时间对饱和流率的影响, 结果表明没有显著影响[1]; Chen等对左转车道的饱与流率研究表明HCM与日本交通工程协会(Japan Society of Traffic Engineering, JSTE)模型均高估了左转车道饱和流率, 估计误差随着行人和自行车流量的增加而增大[2]; Shao等研究了交通组成、车道宽度、车道坡度以及左转弯半径对饱和流率的影响, 并给出了理想饱和流率推荐值与相关校正系数的计算方法[3]; 杨晓光等根据实际数据分析了饱和车头时距分布以及车道宽度对饱和车头时距及饱和流率的影响, 建立了车道宽度、饱和车头时距以及饱和流率的关系模型[4]; 刘意等建立了饱和流率、车道宽度、大车率以及车道数的一元回归关系模型及三元回归线性模型[5]; 王进等采用回归分析的方法提出许可相位下左转专用车道、直左共用车道和一般直行车道的饱和流率修正系数计算模型[6]; 董方等对饱和流率计算模型进行了研究[7-12]。上述研究主要在分析实际数据的基础上研究饱和流率的影响因素, 并以此为基础建立饱和流率计算模型。然而因为影响因素很多, 所以饱和流率并不是恒定不变的, 而是随交通条件不断变化。模型计算法的计算结果为恒定值, 不能反映饱和流率动态变化的特点, 而且通常只在不能通过实地测量取得时才考虑使用模型计算方法进行估计, 如新建交叉口。

    本文充分利用了目前城市道路交叉口广泛布设的线圈检测器所记录的数据来提取饱和流率, 该方法能够适应交通状态不断变化的特点, 动态采集交叉口饱和流率, 有助于实时掌握交叉口的通行能力, 为交通控制方案的制定提供支持。

    经车道划分并加以渠化以后, 交叉口进口车道饱和流率随车道数以及渠化方案的变化而变化, 所以实际获取交叉口某一进口车道的饱和流率时, 必须分别计算各条进口车道的饱和流率, 然后再累加成进口道的饱和流率。各车道饱和流率S可以通过饱和车头时距h得到

    S=3600/h(1)

    在实际研究中, 取队列中的第几辆车作为计算饱和车头时距的起点并没有一个确切的结论, 取第4~7辆车的都存在[1, 13-15]。一般认为队列中初始的3~5辆车存在启动延误, 饱和车头时距的计算起点取第4~6辆车比较符合实际情况, 因此, 本文取第4辆车作为饱和车头时距的计算起点。另外, 提取饱和车头时距时规定, 能够提取饱和车头时距的车道必须满足周期内到达的车辆数不少于7 veh, 而且处于饱和流状态的车辆数不少于4 veh。当前周期的饱和车头时距提取算法如下。

    Step 1:当该周期绿灯开启时, 记录绿灯开始时间t0, 用以计算第1辆车的车头时距。

    Step 2:记录各辆车依次经过线圈的到达时间和离开时间, 车辆i的到达与离开时间分别为tiTi

    Step 3:检验周期内该车道到达的车辆数是否大于等于7 veh, 如果少于7 veh, 则停止计算, 否则继续。

    Step 4:根据车辆到达与离开线圈时间, 计算各辆车占用线圈的时间, 第1辆车的占用时间O1取车辆离开线圈时间T1减去绿灯开启时间t0, 后续车辆i的线圈占用时间Oi

    Οi=Τi-ti(2)

    Step 5:根据车辆离开线圈时间计算饱和车头时距, 以车尾通过线圈为准。第1辆车的车头时距h1取该车辆离开线圈时间T1减去绿灯开启时间t0, 后续车辆i的饱和车头时距hi

    hi=Τi-Τi-1(3)

    Step 6:当车辆的线圈占用时间大于大型车的判断阈值时, 判断经过车辆为大型车, 否者为小型车。

    Step 7:在判断车辆类型的基础上, 根据大小车型的饱和车头时距判断阈值分别判断车辆是否处于饱和流状态, 判断从第4辆车开始的各车辆车头时距是否为饱和车头时距, 当某辆车的车头时距超过了设定的饱和车头时距阈值, 则认为从该辆车开始, 之后的车流为非饱和车流。

    Step 8:判断处于饱和流状态的车辆数量是否大于等于4 veh, 如果满足此条件, 则继续计算, 否则停止计算。

    Step 9:计算从第4辆车开始的各饱和车头时距的平均值, 作为当前周期的饱和车头时距ha

    ha=1n-3ni=4hi(4)

    式中: n为处于饱和流状态的最后一辆车的位置。

    至此, 当前周期的饱和车头时距已经得到。但是, 由于交通状况等条件的影响, 例如大型车的比例突然增加, 会造成当前周期的饱和车头时距产生突变, 如果直接将通过这样的饱和车头时距计算出来的饱和流率用于交叉口服务水平评价或信号配时设计, 会产生很大偏差。为了减少当前周期饱和车头时距突变带来的不良影响, 采用指数平滑法对历史平均饱和车头时距和当前周期饱和车头时距进行平滑处理

    ˉha=0.25ha+0.75Η(5)

    式中:ha为当前周期平滑处理后的平均饱和车头时距; H为历史平均饱和车头时距, 对于当前周期而言, H取前一周期平滑后的平均饱和车头时距。

    根据式(1)可以得到饱和流率为

    S=3600/ˉha(6)

    由饱和流率提取方法可知, 计算饱和车头时距的关键是判定车辆是否处于饱和流状态, 即判定饱和车流的最后一辆车的位置。本文通过相关文献以及建立仿真模型来确定判断饱和车流的阈值。

    大型车性能条件与小型车有很大的差异, 造成大型车与小型车的饱和车头时距有很大的差别, 因而准确判断出饱和车流中的大型车十分重要。当有大型车经过线圈检测器时, 饱和车头时距会有所增加, 相关文献发现大型车辆的饱和车头时距比小型车多2.0~5.0 s[1, 7, 9, 16], 当线圈检测到当前车辆的占用时间大于小型车平均占用时间的2倍时, 判断该车辆为大型车[14]

    使用VISSIM软件建立模型进行仿真, 确定饱和车头时距的判断阈值。建立一个双车道的仿真路段模型, 路段长300 m, 每车道宽3.5 m, 坡度为0, 车道功能为直行车道, 信号周期为60 s, 绿灯时间为30 s, 输入车辆全为小型车辆, 仿真10个周期, 在系统预热2个周期后开始记录数据, 计算其中一条车道每个周期处于饱和状态车辆的车头时距。仿真结果见图 12, 由结果可知小型车的饱和车头时距总体波动不大, 不过也有少数偏离均值的程度较大, 考虑将此设为饱和车头时距的上限。仿真得到的平均饱和车头时距为1.91 s, 饱和状态下, 偏离均值最大的饱和车头时距为2.98 s, 所以将小型车饱和车头时距的判断阈值定为平均饱和车头时距加1 s比较符合实际情况。

    图  1  各周期饱和车头时距
    Figure  1.  Saturation headways of each cycle
    图  2  平均饱和车头时距
    Figure  2.  Average saturation headways

    综合上述分析, 饱和流率提取算法中Step 6的大型车判断阈值取小型车平均占用时间的2倍, 即当线圈检测到当前车辆的占用时间大于小型车平均占用线圈时间的2倍时, 判断该车辆为大型车。小型车平均占用线圈时间取上一周期处于饱和状态的小型车线圈占用时间的平均值。算法中Step 7的饱和车头时距判断阈值取值如下: 小型车饱和车头时距的判断阈值定为历史平均饱和车头时距加1 s, 大型车饱和车头时距的判断阈值为历史平均饱和车头时距加5 s。根据饱和车头时距的阈值即可判断车流是否为饱和车流, 如果当前车辆的车头时距超过了相应车型的阈值, 则认为从该辆车起为非饱和车流, 计算平均饱和车头时距时从该辆车截断。

    为了验证算法的有效性, 仍然采取上文建立的仿真模型进行仿真, 输入车辆中包括大型车, 根据仿真得到的数据计算车头时距, 并用车头时距判断阈值来判断车流是否处于饱和状态, 最后计算出当前周期平均饱和车头时距, 并和历史平均值进行平滑处理。

    系统预热5个周期后记录接下来5个周期的车流数据。为了得到初始输入参数, 先观察判断第1个周期处于饱和状态的车流, 根据仿真数据计算出该周期的平均车头时距以及小型车平均占有时间, 以此作为后续计算的初始化参数。

    第1周期仿真结果见图 3, 其中第1~7辆车处于饱和状态, 因而计算第4~7辆车的平均饱和车头时距, 即当前周期平均饱和车头时距为2.02 s, 小型车平均占用线圈时间为0.66 s, 初始化参数Η1=ˉha1=ha=2.02sS为1 782 veh·h-1, O为0.66 s。

    图  3  第1周期仿真结果
    Figure  3.  Simulation result of the first cycle

    第2周期仿真结果见图 4, 第7辆车的占用线圈时间O7 > 2O, 因而判断为大型车辆, h7=5.19 < H1+5=7.02 s, 判断该车辆处于饱和流状态, 第10辆车的车头时距h10=7.15 > H1+1=3.02 s, 则认为从该车辆起的车流为非饱和车流, 故计算第4~9辆车的平均饱和车头时距, ha为2.50 s, 该周期O0.67sˉha2=0.25ha+0.75Η1=2.14s, 此时可以计算出饱和流率S为1 682 veh·h-1, 令Η2=ˉha2, 进行下一周期计算。

    图  4  第2周期仿真结果
    Figure  4.  Simulation result of the second cycle

    第3周期仿真结果见图 5, O8 > 2O, 因而判断第8辆车为大型车, h8=5.36 < H2+5=7.14 s, 判断该车辆处于饱和流状态, O10 > 2O, 第10辆车也为大型车, 车头时距h10=8.57 > H2+5=7.14 s, 则认为从该车辆起的车流为非饱和车流, 故计算第4~9辆车的平均饱和车头时距ha为2.56 s, O0.62sˉha3=0.25ha+0.75Η2=2.25sS为1 600 veh·h-1, 令Η3=ˉha3, 转入下一周期计算。

    图  5  第3周期仿真结果
    Figure  5.  Simulation result of the third cycle

    第4周期仿真结果见图 6, O3 > 2O, 因而判断第3辆车为大型车, O11 > 2O, 第11辆车也为大型车, 车头时距h11=7.58 > H3+5=7.25 s, 则认为从该车辆起的车流为非饱和车流, 故计算第4~10辆车的平均饱和车头时距ha为1.76 s, O0.66sˉha4=0.25ha+0.75Η3=2.13sS为1 690 veh·h-1, 令Η4=ˉha4, 转入下一周期的计算。

    图  6  第4周期仿真结果
    Figure  6.  Simulation result of the fourth cycle

    第5周期仿真结果见图 7, 该周期所有车辆都处于饱和状态, 所以计算第4~16辆车的平均饱和车头时距ha=1.72 s, O0.60sˉha5=0.25ha+0.75Η4=2.03s; 可以计算出饱和流率S为1 773 veh·h-1

    图  7  第5周期仿真结果
    Figure  7.  Simulation result of the fifth cycle

    由以上算例可见, 本文提出的车型判断阈值在计算的4个周期均能够判断出当前经过的车辆是否为大型车, 能够反映大型车对饱和流率的影响。采用的饱和车头时距判定阈值也能够判断出车辆是否处于饱和流状态, 从而找出处于饱和状态的最后一辆车, 以计算平均饱和车头时距。大型车的存在对平均饱和车头时距有较大影响, 本文采用的指数平滑法可以消除某一周期饱和车头时距骤变带来的影响, 能够阻止其放大大型车的效果。5个周期饱和车头时距和饱和流率变化见图 8, 可以看出, 在第2周期, 由于大型车辆的存在导致当前周期平均饱和车头时距增加了23.8%, 但是平滑处理后饱和车头时距变化没有那么剧烈, 仅增加了5.9%, 饱和流率的波动也是在较小的范围内, 减少了5.6%, 而第4周期平均饱和车头时距较第3周期骤减了31.3%, 但平滑后的饱和车头时距仅减少了5.3%, 饱和流率增加了5.6%。

    图  8  五个周期饱和车头时距和饱和流率变化
    Figure  8.  Saturation headways and saturation flow rates of five cycles

    将以上计算结果与HCM模型计算结果进行了比较。HCM计算法在理想饱和流率的基础上进行一系列修正, 对于该算例, 根据HCM中模型校正系数的取值方法, 理想饱和流率取1 900 pcu·h-1, 宽度修正系数为0.98, 大车修正系数为0.95, 停车修正系数为0.95, 其他修正系数取1, 可以求得饱和流率S为1 680 veh·h-1, 显然模型计算方法无法反应饱和流率随交通状态等因素变化而变化的特征。图 8中5个周期饱和流率分别为1 782、1 682、1 600、1 690、1 773 veh·h-1, 本文动态提取方法能够反应出饱和流率随交通状态变化而变化的特征。

    通过交叉口停车线前线圈采集的交通数据来获取饱和车头时距, 本文提出了获取饱和车头时距的详细方法, 为了避免当前周期饱和车头时距突变带来的不良影响, 采用指数平滑的方法对历史平均饱和车头时距和当前周期饱和车头时距进行平滑, 得到均值, 进而计算饱和流率。

    动态提取方法的关键环节为判断车流是否处于饱和流状态, 本文提出判断大型车的线圈占用时间阈值, 提出了判断车流是否处于饱和流状态的阈值。对饱和车头时距判断阈值的深入研究, 可通过调查不同城市的不同道路条件、不同转向、不同交通条件下各种车型的饱和车头时距分布情况进行综合分析。

    本文的动态提取方法利用的是目前广泛使用的感应线圈检测器, 基本无需额外成本即可实现, 同时可以满足动态采集饱和流率的需求, 而且能避免大量的人工调查成本, 具有一定的实用价值。

  • 图  1  各周期饱和车头时距

    Figure  1.  Saturation headways of each cycle

    图  2  平均饱和车头时距

    Figure  2.  Average saturation headways

    图  3  第1周期仿真结果

    Figure  3.  Simulation result of the first cycle

    图  4  第2周期仿真结果

    Figure  4.  Simulation result of the second cycle

    图  5  第3周期仿真结果

    Figure  5.  Simulation result of the third cycle

    图  6  第4周期仿真结果

    Figure  6.  Simulation result of the fourth cycle

    图  7  第5周期仿真结果

    Figure  7.  Simulation result of the fifth cycle

    图  8  五个周期饱和车头时距和饱和流率变化

    Figure  8.  Saturation headways and saturation flow rates of five cycles

  • [1] KHOSLA K, WILLIAMS J C. Saturation flow at signalized intersections during longer green time[J]. Transportation Research Record, 2006(1978): 61-67.
    [2] CHEN P, NAKAMURA H, ASANO M. Saturation flow rate analysis for shared left-turn lane at signalized intersections in Japan[J]. Procedia Social and Behavioral Sciences, 2011(16): 548-559.
    [3] SHAO Chang-qiao, RONG Jian, LIU Xiao-ming. Study on the saturation flow rate and its influence factors at signalized intersections in China[J]. Procedia Social and Behavioral Sciences, 2011(16): 504-514.
    [4] 杨晓光, 庄斌, 李克平. 信号交叉口饱和流率和启动延误的影响分析[J]. 同济大学学报: 自然科学版, 2006, 34(6): 738-743. doi: 10.3321/j.issn:0253-374X.2006.06.006

    YANG Xiao-guang, ZHUANG Bin, LI Ke-ping. Analysis of saturation flow rate and delay at signalized intersection[J]. Journal of Tongji University: Natural Science, 2006, 34(6): 738-743. (in Chinese). doi: 10.3321/j.issn:0253-374X.2006.06.006
    [5] 刘意, 张永, 任刚. 信号交叉口直行车道饱和流率研究[J]. 交通运输工程与信息学报, 2010, 8(3): 62-67. doi: 10.3969/j.issn.1672-4747.2010.03.012

    LIU Yi, ZHANG Yong, REN Gang. Analysis of saturation flow rate of the through lanes at signalized intersections[J]. Journal of Transportation Engineering and Information, 2010, 8(3): 62-67. (in Chinese). doi: 10.3969/j.issn.1672-4747.2010.03.012
    [6] 王进, 杨晓光. 许可相位下饱和流率修正系数研究[J]. 公路交通科技, 2010, 27(9): 94-99. doi: 10.3969/j.issn.1002-0268.2010.09.017

    WANG Jin, YANG Xiao-guang. Study on correction coefficient of saturation flow rate under permissive phase[J]. Journal of Highway and Transportation Research and Development, 2010, 27(9): 94-99. (in Chinese). doi: 10.3969/j.issn.1002-0268.2010.09.017
    [7] 董方. 北京城市交叉口饱和流量及通行能力研究[D]. 北京: 北京工业大学, 2010.

    DONG Fang. Research on saturated flow and capacity of signalized intersections in Beijing[D]. Beijing: Beijing University of Technology, 2010. (in Chinese).
    [8] 赵林, 邵长桥, 荣建. 交通组成对信号交叉口饱和流率影响研究[J]. 公路工程, 2008, 33(4): 88-91. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGL200804023.htm

    ZHAO Lin, SHAO Chang-qiao, RONG Jian. Analysis of traffic composition affecting saturation flow at signalized intersection[J]. Highway Engineering, 2008, 33(4): 88-91. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGL200804023.htm
    [9] 赵林, 邵长桥, 荣建, 等. 信号交叉口右转车道饱和流率部分影响因素研究[J]. 武汉理工大学学报: 交通科学与工程版, 2011, 35(5): 873-876. doi: 10.3963/j.issn.1006-2823.2011.05.001

    ZHAO Lin, SHAO Chang-qiao, RONG Jian, et al. Analysis of some factors affecting right-turning saturation flow at signalized intersection[J]. Journal of Wuhan University of Technology: Transportation Science and Engineering Edition, 2011, 35(5): 873-876. (in Chinese). doi: 10.3963/j.issn.1006-2823.2011.05.001
    [10] 赵林, 邵长桥, 王子雷, 等. 设计时速与直行车道数对饱和流率的影响分析[J]. 交通标准化, 2009(1): 89-92. https://www.cnki.com.cn/Article/CJFDTOTAL-JTBH200901027.htm

    ZHAO Lin, SHAO Chang-qiao, WANG Zi-lei, et al. Influence of design speed and number of through lanes on saturation flow rate[J]. Communications Standardization, 2009(1): 89-92. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JTBH200901027.htm
    [11] 赵星, 任刚. 信号交叉口车头时距特性分析[J]. 交通运输工程与信息学报, 2010, 8(1): 103-108. https://www.cnki.com.cn/Article/CJFDTOTAL-JTGC201001022.htm

    ZHAO Xing, REN Gang. Analysis of vehicle headway characters at a signalized intersection[J]. Journal of Transportation Engineering and Information, 2010, 8(1): 103-108. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JTGC201001022.htm
    [12] 吴震, 杨晓光. 信号交叉口右转车道饱和流量研究[J]. 同济大学学报: 自然科学版, 2008, 36(1): 37-41. https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ200801010.htm

    WU Zhen, YANG Xiao-guang. Research on saturation flow of right-turning lane at signalized intersection[J]. Journal of Tongji University: Natural Science, 2008, 36(1): 37-41. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-TJDZ200801010.htm
    [13] LI H, PREVEDOUROS P D. Detailed observations of saturation headways and start-up lost times[J]. Transportation Research Record, 2002(1802): 44-53.
    [14] HENDERSON I, WOOD K. Computerised real-time saturation flow measurement for signalised traffic junctions[R]. London: Transport Research Laboratory, 2005.
    [15] BESTER C J, MEYERS W L. Saturation flow rates[C]//SATC. Proceedings of the26th Southern African transport conference. Pretoria: SATC, 2007: 560-568.
    [16] 刘明君. 基于混合交通流的信号交叉口机动车车头时距研究[D]. 北京: 北京交通大学, 2010.

    LIU Ming-jun. Modeling on discharge headway with heterogeneous traffic at signalized intersections[D]. Beijing: Beijing Jiaotong University, 2010. (in Chinese).
  • 加载中
图(8)
计量
  • 文章访问数:  825
  • HTML全文浏览量:  256
  • PDF下载量:  1131
  • 被引次数: 0
出版历程
  • 收稿日期:  2012-08-19
  • 刊出日期:  2013-02-25

目录

/

返回文章
返回