2018年中国研究生数学建模竞赛D题

第一篇：2018年中国研究生数学建模竞赛D题

2018年中国研究生数学建模竞赛D题

基于卫星高度计海面高度异常资料 获取潮汐调和常数方法及应用

1.潮汐潮流现象的研究意义

海洋潮汐是在天体引潮力作用下形成的长周期波动现象，在水平方向上表现为潮流的涨落，在铅直方向上则表现为潮位的升降。潮汐潮流运动是海洋中的基本运动之一，它是动力海洋学研究的重要组成部分，对它的研究直接影响着波浪、风暴潮、环流、水团等其他海洋现象的研究，在大陆架浅海海洋中，对潮汐潮流的研究更具重要性。

海岸附近和河口区域是人类进行生产活动十分频繁的地带，而这个地带的潮汐现象非常显著，它直接或间接地影响着人们的生产和生活。潮汐潮流工作的开展和研究，可为国防建设、交通航运、海洋资源开发、能源利用、环境保护、海港建设和海岸防护提供资料。例如，沿海地区的海滩围垦、农田排灌，水产的捕捞和养殖，制盐，海港的选址及建设，以至于潮能发电等活动，无不与潮汐潮流现象有着密切的关系。2.潮汐潮流数值模拟所面临的问题

区域海洋潮汐的数值模拟需要提供开边界的水位调和常数，而开边界的水位调和常数，或者来源于观测、或者来源于全球海洋潮汐的数值模拟；而全球海洋潮汐的数值模拟，相当耗费资源。虽然目前有国外学者或研究机构，能够提供区域海洋潮汐的调和常数，但实质上的评价结果难以令人满意。

从区域海洋潮汐的数值模拟的现状来讲，四个主要分潮（M2、S2、K1、O1）的单一分潮的数值模拟与同化可以得到令人满意的结果，但其它分潮（、、、等）的单一分潮的数值模拟与同化，结果却差强人意；这意味着其它分潮的数值模拟，只有与四个主要分潮同时进行数值模拟，才能得到可以接受的结果。从具体操作来讲，其它分潮由于相对较弱，导致模拟结果的精度难以提高。

长周期分潮（、、、）的获取，目前已有基于全球长周期分潮数值模拟手段的报道，但其面临的困境，与其它较弱分潮面临的困境没有差别。

从各分潮的调和常数获取的发展史来说，通过对已有观测结果进行插值曾经是首选，但发展过程中逐渐被数值模拟方法所取代。高度计资料的出现，引发部分学者开展了插值方法的研究，并取得了一些值得一提的结果，尽管被所谓的主流方式淹没，但也难掩其光芒所在。鉴于目前已有高度计资料作为支持，其它分潮及长周期分潮的调和常数获取的插 值方法研究大有可为。3.资料描述 3.1 地形数据

地形数据来自ETOP5，全球的分辨率为5'5'，图1的区域是 2~25N，99~122E。

图1 南海地形图

3.2 验潮站资料

中国近海及周边海域770个验潮点的资料，和56个验潮点的资料（是国际上公开的长期验潮站数据分析得到的调和常数），包括9个分潮（M2、S2、K1、O1、N2、K2、P1、Q1、Sa）的潮汐调和常数。

图2显示了上述资料点所在的位置，从图中可以看出上述验潮点主要分布在近岸或岛屿附近。

图2 验潮站资料的分布图

3.3 TOPEX/POSEIDON卫星高度计简介

卫星高度计是一种向卫星下方海洋发射脉冲的雷达，通过测量脉冲经海面反射之后的往返时间，获得卫星距海面的高度。主要用途：利用所得到的海面动力高度同化反演海洋重力场、流场、潮、大地水准面、海洋重力异常；根据回波强度获取风速资料；根据回波波形前沿斜率获取海面有效波高。

TOPEX/POSEIDON卫星是由美国国家航空航天局和法国空间局联合于1992年8月10日发射的，是世界上第一颗专门为研究世界大洋环流而设计的高度计卫星。其轨道高度达1336km，倾角为66°，覆盖面大，保证了资料的连续性。轨道的交点周期（绕地球一圈的时间）为6745.8s，轨道运行127圈以后精确重复，轨道重复周期为9.9156天。相邻最近的轨道之间在赤道上的间隔为360/1272.835。卫星在一个周期内的每一圈分为上行轨和下行轨两条轨道，一个完整的周期内共有254条轨道，沿轨道的两个相邻的星下观测点的距离5.75km。高度计系统的定规精度和测高精度较以前有显著提高，其测量精度约为5cm，是目前观测海面高度精度最高的卫星。

当然，本文只是涉及到TOPEX/POSEIDON卫星高度计资料与潮汐相关的研究，即海面高度异常产品。3.4 南海高度计资料

图3的给出了 2~25N，99~122E，TOPEX/POSEION卫星高度计星下观测点所在的轨道。一共有超过4000个数据点，每个点都对应一个海面高度异常的时间序列，从1992年到2017年，时间跨度为25年。

图3 南海TOPEX/POSEIDON高度计资料的星下轨迹

本题所附文件包括：

1.地形数据来自ETOP5，全球的分辨率为5'5'，此处数据度范围2~25N，99~122E（附数据说明）。

2.中国近海及周边海域验潮点的资料，包括9个分潮（M2、S2、K1、O1、N2、K2、P1、Q1、Sa）的潮汐调和常数。还包括56个验潮点的资料，同样给出了上述9个分潮的潮汐调和常数，是国际上公开的长期验潮站数据分析得到的调和常数（附数据说明）。

3.南海海面高度异常数据文件及说明。

4.调和分析方法简介（包括分潮的Doodson数、分潮角速度和交点因子与订正角m、mm、5列表）4本题所涉及的数据与结果，均应采用国际标准单位，即：时间单位为秒、距离单位为米、速度单位为米每秒等。根据以上介绍及提供的数据文件及表格，请你们团队研究下列问题：

1.根据沿轨道的星下观测点的海面高度异常值，提取所有星下观测点各主要分潮（M2、S2、K1、O1）的潮汐调和常数，注意能有效提取那些分潮的潮汐调和常数取决于相应的资料长度；对提取的潮汐调和常数，应利用潮汐验潮点的调和常数给予评价或检验，并给出评价结果的分析或评价。

2.得到所有星下观测点各主要分潮（M2、S2、K1、O1）的潮汐调和常数，沿轨道作图后，可发现潮汐调和常数在沿轨道方向，在空间有细结构，而此细结构是内潮对正压潮的调制；请设法对沿轨道的各分潮的潮汐调和常数进行正压潮和内潮的分离。

3.设计数据插值或拟合方法给出南海的各主要分潮的同潮图，并利用潮汐验潮点的调和常数给予评价或检验，并给出评价结果的分析或评价。如果你们还有时间和兴趣，还可考虑下列：

4.如果在对沿轨道的潮汐调和常数分离、插值或拟合的过程中，利用了特定的函数进行拟合，是否能够确定出需利用的特定函数的最佳（高）次数？上述结论是否对第3问有启示或帮助。

 本题要求提供可计算出所提交报告中答案的计算程序，所使用的语言和工具不限，但推荐使用CC++、Fortran、Matlab、...。

参考文献

[1] 陈宗镛，潮汐学，北京：科学出版社，1980．

[2] 方国洪，郑文振，陈宗镛，王骥，潮汐和潮流的分析与预报，北京：海洋出版社，1986． [3] 黄祖珂，黄磊，潮汐原理与计算，青岛，中国海洋大学出版社，2005 [4] 孙丽艳：渤黄东海潮汐底摩擦系数的优化研究［硕士学位论文］．青岛：中国海洋大学海洋环境学院，2006．

[5] 范丽丽：风暴潮数值同化研究和高度计资料拟合方法研究［硕士学位论文］．青岛：中国海洋大学海洋环境学院，2011．

[6] Mazzega, P., and M.Berge, 1994，Ocean tides in the Asian semi-enclosed seas from TOPEX/POSEIDON.J.Geophys.Res., 99: 24,867–24,881.[7] Yangi, T., A.Morimoto, and K.Ichikawa, 1997.Co-tidal and co-range charts for the East China Sea and the Yellow Sea derived from satellite altimetric data.J.Oceanography, 53: 303–309.[8] Fan, L.L.，Wang B.，and Lv X.Q.，2011，Cotidal Charts near Hawaii Derived from TOPEX/Poseidon Altimetry Data，J Atmos Ocean Technol，28，606-614． [9] Fang, G.H.，Wang Y.G., Wei Z.X., Choi B.H., Wang X.Y., and Wang J., 2004，Empirical cotidal charts of the Bohai, Yellow, and East China Seas from 10 years of TOPEX/Poseidon altimetry.J.Geophys.Res.，Vol.109:C11006,doi: 1029/2004JC002484.[10] Wang Y.H.，FANG G.H.，Wei Z.X.，Wang Y.G.，Wang X.Y., Xu X.Q.，Cotidal charts and tidal power input atlases of the global ocean from TOPEX/Poseidon and JASON-1 altimetry，Acta Oceanol.Sin.，2012，31，4，11-23．

第二篇：2018年中国研究生数学建模竞赛E题

2018年中国研究生数学建模竞赛E题

多无人机对组网雷达的协同干扰

组网雷达系统是应用两部或两部以上空间位置互相分离而覆盖范围互相重叠的雷达的观测或判断来实施搜索、跟踪和识别目标的系统，综合应用了多种抗干扰措施，具有较强的抗干扰能力，因而在军事中得到了广泛应用。如何对组网雷达实施行之有效的干扰，是当今电子对抗界面临的一个重大问题。

诸多干扰方式中较为有效的是欺骗干扰，包括距离欺骗、角度欺骗、速度欺骗以及多参数欺骗等。本赛题只考虑距离假目标欺骗，其基本原理如图1所示，干扰机基于侦察到的敌方雷达发射电磁波的信号特征，对其进行相应处理后，延迟（或导前）一定时间后再发射出去，使雷达接收到一个或多个比该目标真实距离靠后（或靠前）的回波信号。

图 1

对雷达实施距离多假目标欺骗干扰示意图

在组网雷达探测跟踪下，真目标和有源假目标在空间状态（如位置、速度等）上表现出显著的差异：对于真目标，其空间状态与雷达部署位置无关，在统一坐标系中，各雷达探测出的真目标空间状态是基本一致的，可以认为它们是源自于同一个目标（同源）；对于有源假目标，它们存在于雷达与干扰机连线以及延长线上，其空间状态由干扰机和雷达部署位置共同决定，不同雷达量测到的有源假目标的空间状态一般是不一致的，有理由认为其来自于不同目标（非同源），利用这种不一致性就可以在组网雷达信息融合中心将假目标有效剔除。这种利用真假目标在组网雷达观测下的空间状态差异来进行假目标鉴别的思想简称为“同源检验”，它是组网雷达对真假目标甄别的理论依据。

为了能对组网雷达实施有效干扰，现在可利用多架无人机对组网雷达协同干扰。如图2所示，无人机搭载的干扰设备对接收到的雷达信号进行相应处理后转发回对应的雷达，雷达接收到转发回的干扰信号形成目标航迹点信息，传输至组网雷达信息融合中心。由于多无人机的协同飞行，因此在融合中心就会出现多部雷达在统一坐标系的同一空间位置上检测到目标信号，基于一定的融合规则就会判断为一个合理的目标航迹点，多个连续的合理目标航迹点就形成了目标航迹，即实现了一条虚假航迹。通过协同控制无人机的飞行航迹，可在敌方的组网雷达系统中形成一条或多条欺骗干扰航迹，迫使敌方加强空情处置，达到欺骗目的。

图 2 多无人机协同干扰组网雷达系统示意图

某组网雷达系统由5部雷达组成，雷达最大作用距离均为150km，也就是只能对距雷达150 km范围内的目标进行有效检测。5部雷达的地理位置坐标分别为雷达1（80，0，0），雷达2（30，60，0），雷达3（55，110，0），雷达4（105，110，0），雷达5（130，60，0）（单位：统一为km）。雷达将检测到的回波信号经过处理后形成航迹点状态信息（本赛题主要关心目标的空间位置信息）传输到融合中心，融合中心对5部雷达获取的目标状态信息进行“同源检验”，只要有3部以上雷达的航迹点状态信息通过了同源检验，即至少有3部雷达同一时刻解算出的目标空间位置是相同的，融合中心就将其确定为一个合理的航迹点，20个连续的经融合中心确认的航迹点形成的合理航迹，将被组网雷达系统视为一条真实的目标航迹。所谓合理的航迹是要满足相应的目标运动规律，无论是运动速度还是转弯半径等均应在合理的范围内。

现考虑多架无人机对组网雷达系统的协同干扰问题。无人机的飞行速度控制在120km/h~180km/h，飞行高度控制在2000m~2500m，最大加速度不超过10m/s2。由于安全等因素的考虑，无人机间距需控制在100 m以上。鉴于无人机的RCS较小，也采用了若干隐身技术，在距雷达一定距离飞行时，真实目标产生的回波不能被雷达有效检测（本赛题可不考虑无人机产生的真实目标回波）；干扰设备产生的欺骗干扰信号经过了放大增强环节，能保证被雷达有效检测到。每架无人机均搭载有干扰设备，可独立工作。同一时刻一架无人机只能干扰一部雷达，但可在该部雷达接收机终端（雷达屏幕上）产生多个目标点，这些目标点均位于雷达与无人机连线以及延长线上，距雷达距离超过150 km的假目标信息直接被雷达系统删除；同一时刻多架无人机可以干扰同一部雷达。雷达同一时刻接收的多个目标点的状态信息均同时传送到信息融合中心。每架无人机不同时刻可干扰不同雷达。同一条航迹不同时刻的航迹点，可以由组网雷达系统中不同的三部雷达检测确定。

请建立相应的数学模型，研究下列问题：

（1）附件1给出了一条拟产生的虚假目标航迹数据，该虚假航迹数据包含20个时刻的虚假目标位置坐标信息，时间间隔为10秒。为实现较好的干扰效果，现限定每架无人机在该空域均做匀速直线运动，航向、航速和飞行高度可在允许范围内根据需要确定。请讨论如何以最少数量的无人机实现附件1要求的虚假目标航迹，具体分析每一架无人机的运动规律和相应的协同策略。

问题（1）的解算结果，需按附件4中规定格式具体给出每一架无人机对应时刻的空间位置坐标，存入文件“E队号_1.xls”中，作为竞赛附件单独上传竞赛平台，是竞赛论文评审的重要依据。

（2）对雷达实施有源假目标欺骗干扰时，干扰设备可同时转发多个假目标信息（本赛题限定每一架无人机同一时刻至多产生7个假目标信息），但它们均存在于雷达与无人机连线以及延长线上，延迟（或导前）的时间可根据实际需要确定。该组网雷达系统的每一部雷达的数据更新率为10秒（可直观理解为每间隔10秒获得一批目标的空间状态数据，无人机转发回对应雷达的假目标信息能及时获取）。协同无人机编队可产生出多条虚假航迹，以实现更好的干扰效果。实际中无人机可机动飞行，但为控制方便，无人机尽可能少做转弯、爬升、俯冲等机动动作，转弯半径不小于250m。请讨论由9架无人机组成的编队在5分钟内，完成附件1要求的虚假航迹的同时，至多还可产生出多少条虚假的航迹。给出每一架无人机的运动规律，并分析每一条虚假航迹的运动规律和合理性。

问题（2）的解算结果，需按附件4中规定的格式具体给出每一架无人机对应时刻的空间位置坐标和每一条虚假航迹的相关数据，存入文件“E队号_2.xls”中，作为竞赛附件单独上传竞赛平台。

（3）当组网雷达系统中的某部雷达受到压制干扰或其它因素的干扰时，可能在某些时刻无法正常获取回波信号，此时组网雷达系统信息融合中心可以采用下面的航迹维持策略：若之前与受干扰的雷达联合检测到目标的另2部雷达没有受到干扰，正常检测到回波信号，那么在融合中心就对这两部雷达检测的目标航迹点信息进行同源检验，若通过亦视为是合理的目标航迹点；若一条航迹中这类航迹点的个数不超过3个时（该航迹的其余航迹点仍需通过前面规定的“同源检验”），该航迹就被继续保留。针对上述航迹维持策略，协同无人机编队的飞行，有可能产生更多的虚假航迹。该组网雷达系统的每一部雷达的数据更新率仍为10秒。重新讨论由9架无人机组成的编队在5分钟内，完成附件1要求的虚假航迹的同时，至多还可产生出多少条虚假的航迹。给出每一架无人机的运动规律和协同策略，分析每一条虚假航迹的运动规律和合理性。

附件1 问题1虚假航迹点坐标数据.xls 附件2 问题1的结果.xls 附件3 问题2的结果.xls 附件4 对问题1、2提交结果的规定（请认真阅读，严格按照要求完成）

第三篇：2017年中国研究生数学建模竞赛F题（推荐）

地下物流系统网络

背景

交通拥堵是世界大城市都遇到的“困局”之一。2015年荷兰导航经营商TomTom发布了全球最拥堵城市排名，中国大陆有十个城市位列前三十名。据中国交通部2014年发布的数据，我国交通拥堵带来的经济损失占城市人口可支配收入的20%，相当于每年国内生产总值(GDP)损失5~8%。15座大城市的居民每天上班比欧洲发达国家多消耗28.8亿分钟。大量研究表明：“时走时停”的交通导致原油消耗占世界总消耗量的20%。高峰期，北京市主干线上300万辆机动车拥堵1小时所需燃油为240万~330万升。2015年城市交通规划年会发布数据显示：在石油消费方面，我国交通石油消费比重占到了消费总量的54%，交通能耗已占全社会总能耗10%以上，并逐年上升。高能耗也意味着高污染和高排放。

导致城市交通拥堵的主要原因是交通需求激增所带来的地面道路上车辆、车次数量巨增，其中部分是货物物流的需求增长。尽管货车占城市机动车总量的比例不大，但由于货运车辆一般体积较大、载重时行驶较慢，车流中如果混入重型车，会明显降低道路的通行能力，因此，其占用城市道路资源的比例较大。如北京，按常规的车辆换算系数（不同车辆在行驶时占用道路净空间的程度），货运车辆所占用的道路资源达40%。因此，世界各国都在为解决城市交通和环境问题进行积极探索，而处理好货运交通已成为共识。大量实践证明，仅通过增加地面交通设施来满足不断增长的交通需求，既不科学也不现实，地面道路不可能无限制地增加。因此“统筹规划地上地下空间开发”势在必行，“地下物流系统”正受到越来越多发达国家的重视。概念

地下物流系统（Underground Logistics System——ULS）是指城市内部及城市间通过类似地铁的地下管道或隧道运输货物的运输和供应系统。它不占用地面道路，减轻了地面道路的交通压力，从而缓解城市交通拥堵；它采用清洁动力，有效减轻城市污染；它不受外界条件干扰，运输更加可靠、高效。地面货车的减少同时带来巨大的外部效益，如路面损坏的修复费用，环境治理的费用，可以用于补偿地下物流系统建设的高投资。

图1 日本ULS概念图

图2 德国ULS概念图

图3 地下物流系统的技术形式

图1和图2描绘的分别是日本和德国对ULS的设计构想图。图1进一步描绘了地下节点的运作状况；图2侧重于ULS与其他地下设施的平行关系，在复杂的地下空间条件下实施ULS的空间并不富余。图3中展示的分别是运载车辆形式设计及运输通道的形式。

然而ULS的研究与实践还刚刚兴起，尚无成熟、成功的案例可供借鉴。尽管国内外在发展地下物流的必要性、可行性、技术系统和建造管理等方面取得了不少重要的研究成果，但地下物流系统是涉及地上与地下、物流与工程、技术与管理等多个学科领域的复杂巨系统，而且相比于其他物流方式，ULS复杂，造价高，风险大，不成功则损失巨大，而现有研究成果还处于“试验”阶段，相对单一，尚无系统的理论支撑其走向实际应用，尤其是迄今为止世界上还没有一个城市规模的成功案例，失败的情况反而有几个。因此，亟待进行城市规模的前瞻性研究，为发展地下物流系统的理论和实践进行探索。

我国人口众多、大城市密集、交通状况不佳已经到了迫切需要改善的程度，而且我国城市地铁网络的建设、高铁公路隧道的建造正大规模地进行，地下空间开发利用的规 模和速度已居世界前列，地下工程的技术水平也已基本满足需要，ULS应该提上议事日程。撇开可行性、工程技术问题，构建地下物流系统网络是建设“地下物流系统”必不可少的关键步骤。

地下物流系统的基本特征：

1.建造与运营成本高；网络形式多样（环型、树形、网状等等）且随区域特征变化。2.地下物流无法完全替代地面物流，“多式联运”是其重要运输组织方式。多式联运：由两种及两种以上的交通工具相互衔接、转运而共同完成的运输过程统称为复合运输，我国习惯上称之为多式联运。如地下物流-公路运输，公路-铁路。

3.对于货物单

一、流量大且稳定、地面交通差的区域应优先建立线路。

4.地下物流系统可以适应不同货物的运输需求，管道或隧道直径可根据需求进行设计。运行速度可以达到20-60 公里/小时，同一线路上同向两班车运行间隔2分钟以上。

5.地下物流网络由一级、二级节点和节点间地下通道构成。各级节点均与地面衔接并实现多式联运。一级节点与物流园区相连且采用10吨的大型车辆地下运输，并可跨区域调运货物，从地面收发货物总量上限为4000吨，一级节点之间连通。二级节点从地面收发货物总量上限为3000吨，且与非本区域一级节点仅通过本区域一级节点连通。

车辆参数（供参考）

 一班车一般由四至八节车辆构成；

 两种型号的车辆设计载重分别为10吨、5吨；  运行速度13.5米/秒（约49公里/小时）；  加（减）速度：1 米/秒2；

 电压：三相380伏，频率：50赫兹，单个直线电机感应板电流：460安；  转弯半径70-80米；

 每个节点每12分钟可发车一班（含装卸货物、启动和等待时间），每小时最多发车5班，每天可以运营18小时（每天剩余为检修时间）。

 尽量采用双向双轨（双洞）的隧道形式，在货运量特别大的线路可设置双向四轨（双洞）。题目给出了南京市仙林地区的交通货运区域划分图和相应的货运OD（Origin Destination）流量矩阵（只考虑始发地和目的地的货运流量，不考虑在此之间的途径地）、各区域中心点及区域面积、各区域交通拥堵系数（为简化计算，部分数据经过处理）。其他相关数据可以自行查找，收集你们认为与建立该区域“地下物流系统”网络有关的数据资料。

请完成以下几项任务：

发展城市地下物流网络的两个直接目标：一是缓解交通拥堵直至交通畅通，至少基本畅通；二是降低物流成本。

1.地下物流节点选择；根据你们的观点和该区域的实际情况建立该区域节点选择模型，确定该区域地下物流网络节点群。

计算结果需要但不限于提交：

一、二级节点数及位置、各节点的服务范围（经该节点出、入地面货物的起或讫点形成区域）、各节点实际货运量、各一级节点的转运率。

① 交通拥堵指数取值范围为0至10，每2个数为一等级，分别对应“0-2畅通”、“2-4基本畅通”、“4-6轻度拥堵”、“6-8中度拥堵”、“8-10严重拥堵”五个级别，数值越高，表明交通拥堵状况越严重。由于该地区非人口高度密集区，可以近似认为区域交通拥堵指数与OD数据反映出来的区域总货运量（进、出之和）成正比，考虑到比例关系给出的指数最高值可以大于10。（为简化计算，拥堵指数计算仅根据货运量）。

② 一级节点的转运率（φ）：从物流园区经由最近的一级节点转运至其他所有一级节点的货物量占该物流园区总出货量的百分比，称为该一级节点的转运率。由于需要更换运输车辆，在满足运输要求前提下，转运率低可减少工作量。

③ 考虑到部分区域货运量与面积之比过小，若节点覆盖了某区域中心点即可视为对该区域进行了覆盖。

④ 所有节点的服务半径在 3 公里范围内自由选择，节点间距离不受限制。⑤ 货物从二级节点至地面后采用人力或小型车辆在节点服务区域内进行运输，可认为不影响交通。

⑥ 进出4个物流园区的货物尽最大可能放入地下运输，区域内部的货物根据拥堵情况自行考虑。

2.地下通道网络设计；请你们在地下物流网络节点群的基础上选择合适的地下路线以建立该区域的“地下物流系统”网络。在转运率变化不大的情况下，若考虑优化网 络，可适当调整一、二级节点位置。除园区至一级节点的地下通道外其他地下通道均采用5吨的地下运输车辆。

计算结果需要但不限于提交：网络构成（节点及通道位置）；各节点实际货运量；各级通道的位置和实际流量。

① 地下节点及通道内的货物每天要清仓。

② 每两个相邻节点间地下物流通道双向尺寸一致，以单向流量较大的为设计原则。③ 要求总成本最小。每天的总成本由货物的运输成本和地下物流隧道与节点的折旧构成。

假设每吨货每公里的平均运输成本始终相等，约为1元/吨·公里（已计入车辆和设备折旧），与所经过的隧道尺寸无关。

地下物流隧道与节点的建设成本为：双向四轨（双洞）（10吨）造价5亿元/公里，双向双轨（双洞）（10吨）造价4亿元/公里，双向四轨（双洞）（5吨）造价3.5亿元/公里，双向双轨（双洞）（5吨）造价3亿元/公里，一级、二级节点的建设成本分别约为1.5亿元/个、1亿元/个；通道与节点的设计年限100年，年综合折旧率均为1%。

④ 不考虑物流园区的建设及园区内的地下节点建设，但从园区出发的通道长度需要计入总成本。

3.网络改进；以上是分步设计网络，并未从全局出发，根据你们对运行情况的仿真，上面得到的网络有无修改的必要？能否通过增加、减少节点的个数，调整节点的位置或级别，增加、减少、改变路径的方法縮短货物运输总里程（同时节省运输时间），降低运输成本。

进一步从增强ULS的抗风险能力（如某通道中断，某方向货运量激增）考虑，需要对第二问的ULS作怎样的改动？

4.建设时序与动态优化；“地下物流系统”造价高，风险大，改建困难，所以应做好顶层设计，如果希望考虑满足该市近30年内的交通需求（可以认为需求量每年呈5%增长）并根据建设进度分八年完成“地下物流网络系统”的建设（每年可建设道路长度大致相等），请给出该市“地下物流系统”网络各线路的建设时序及演进过程，与你在第三题中设置的网络有什么差别？并比较优劣。

① 图5为地下物流整体网络可能的演进示意。

② 随着需求量逐年增长，你在第三题中设计的网络是否仍然可行？何时班次和每班车的车厢节数达到满载，如何进行扩容处理？例如：建设期考虑增加线网容量、增设 节点等。

整体

演进过程

图5 地下物流整体网络演进示意

附件包含：

南京市仙林区域的交通货运区域划分图和相应的货运OD流量矩阵（单位：吨，表中数据表示横轴对纵轴的发货量）、各区域面积及中心点坐标（单位：米）、各区域交通拥堵系数。

第四篇：2014研究生数学建模竞赛A题

2014年全国研究生数学建模竞赛A题

小鼠视觉感受区电位信号(LFP)与视觉刺激之间的关系研究

人类脑计划（Human Brain Project, HBP）是继人类基因组计划之后，又一国际性科研计划，其核心是神经信息学(Neuroinformatics）。该研究旨在努力探究数十亿个神经元的信息，以期对知觉、行动以及意识等有更进一步的了解。科学家们预期这是一条开发新技术的好途径，由此可能进一步认识像老年痴呆和帕金森综合症等疾病，有望为各种精神疾病研究出新的治疗方法。此外，该计划还可以更好地为人工智能服务。目前该计划已经取得一定进展，例如在2014年足球世界杯巴西开幕式上，脊髓损伤患者开球就是利用了美国杜克大学的神经生物学家、大脑-计算机界面研究的先驱、巴西人Miguel Nicolelis开发的一种由大脑脑电波控制的外骨骼系统。它既是人类脑计划的结晶，也是2014年Nature十大科学展望领域之一神经科学的代表性成就。1.脑电波介绍

脑是支配人和高级动物活动的司令部和信息中心，神经系统承担着感受外界刺激，产生、处理、传导和整合信号，实现各种认知活动（如知觉、学习、记忆、情绪、语言、意思和思维等），以及运动控制等众多功能。神经系统的基本结构单元是神经元，其放电活动涉及复杂的物理化学过程，表现出丰富的非线性动力学行为。神经系统整体可视为由数目众多的神经元组成的庞大而复杂的信息网络，通过对信息的处理、编码、整合，转变为传出冲动，从而联络和调节机体的各系统和器官的功能。神经元对信息的处理和加工是神经元集群共同完成的，而神经元集群的同步形成较强的电信号就是脑电波。

脑电波（Electroencephalogram,EEG）是大脑在活动时，脑皮质细胞群之间形成电位差，从而在大脑皮质的细胞外产生电流。它们是脑神经细胞的电生理活动在大脑皮层或头皮表面的总体反映，将大脑活动时这种电波的变化记录下来就得到脑电图。2.脑电波的应用

脑电波或脑电图是一种比较敏感的客观指标，不仅用于脑科学的基础理论研究，而且更重要的在于临床实践的应用。事实上，脑电波是诊断癫痫的重要依据，而且对于各种颅内病变，如脑中风、脑炎、脑瘤、代谢性脑病变等的诊断，亦有很大帮助。脑疾病的诊断主要是从脑电波的异常入手，结合临床，对颅内病变进行定位。

另外，脑电波也是人们思维活动的体现。人的大脑是由数以万计的神经交错构成的。神经相互作用时，脑电波模式就是思维状态。人的大脑平均每天产生7万个想法，而且每次神经活动时都会产生轻微的放电。单个神经产生的放电很难从头皮外测量到，但是许多神经共同放电产生的集体电波是可以通过脑电波技术测量到的。因此，我们测量得到的脑电波是由许多神经共同放电产生的集体神经活动决定的。经过近一个世纪的科学实验，神经系统科学领域的专家们已经发现了大脑中控制具体活动的部位。如控制四肢的区域位于大脑的顶部。而负责视力的区域位于大脑的后部。从进化论角度看，大多数动物的大脑也具备这些功能。所以在研究人脑思维时，可以借助于动物实验。

尽管我们知道大脑意识和神经冲动的基础都是电信号，这些电信号是如何精确表达一连串复杂动作一直为科学家所热衷。一旦解析出这种信息密码就可以实现人机互动，甚至用电脑控制动物的行为。Miguel Nicolelis[1]在2011年10月完成了一个猕猴意识控制机械臂的试验。他们将一特殊装置的电极放置到猕猴大脑的运动皮质区和躯体感受皮质区，前者是发出运动信号的区域，后者负责身体其他部分传来的信息。该试验的目标是为瘫痪病人设计可以由大脑控制义肢行为的装置。

该领域的研究目前处于起步阶段，特别是关于人脑的研究，远没有达到实际应用的阶段。所以进行脑电波分析的基础研究具有重要意义。由于人脑的复杂性，在研究脑电波形成机理时，科学工作者大多采用动物实验，比如小鼠等。

一般认为δ波是由大脑皮层和丘脑之间的内在网络所产生的。Ito等人今年4月在“Nature Communication”上发表论文“Whisker barrel cortex delta oscillations and gamma power in the awake mouse are linked to respiration”[2]，他们在对小鼠进行的研究中，发现δ波段的峰值震荡与清醒状态下的小鼠晶须桶状皮层局部场电位的活动被呼吸锁相（delta band oscillation in spike and local field potential activity in the whisker barrel cortex of awake mice is phase locked to respiration）。这也许表明δ波段可能与呼吸有关。（注：锁相是使被控振荡器的相位受标准信号或外来信号控制的一种技术。用来实现与外来信号相位同步，或跟踪外来信号的频率或相位）

此外，大量研究都表明脑电波信号与动物的认知等功能活动有关[3]。3.问题的难点

与脑电波对颅内病变诊断不同的是，利用脑电波分析人（动物）的行为与脑电波之间的关系，并反过来通过脑电波确定或引导人的行为，具有更大的挑战性。脑电波信号是无数神经放电的混合，我们不可能也没有必要将单个神经放电分离出来。宏观意义上，控制某个特定行为或想法的脑电波是一系列众多神经放电的迭加。而我们测量所得到的脑电波信号又是由许许多多构成不同想法和控制行为的脑信号合成的。研究表明这些信号的强弱差别很大。在实现人机交互时，我们也许只关注若干个行为或思想，而对应的脑电波可能很弱。这在信号处理领域，相当于弱信号检测。在数学领域，这可能属于不适定的反问题。显然只有将脑电波信号很好地分离才能从中确定某种脑电波与某种行为相对应。这也可以理解为盲源分离或半盲信号分离问题。但通常的盲源分离技术在这里很难奏效，或误差太大，因为脑电波这一混合信号是由尺度差异很大的信号所构成的。

4.具体问题描述

脑电波来自于大脑内部，一般认为大脑在活动时，脑皮质细胞群之间就会形成电位差，从而在大脑皮质的细胞外产生电流。它是脑神经细胞的电生理活动在大脑皮层或头皮表面的总体反映。而局部场电位（Local Field Potential, LFP）则反映来自神经元网络局部神经核团的活动状态，它也是一种神经集合的协同行为。所以LFP信号是脑内某局部大量神经元树突电位和的综合反映。LFP可能与大脑对行为的控制有关，如呼吸及视觉刺激等。对应于不同行为或思维的脑电波，我们称之为脑电波成份。事实上，当一个人面对一个物品或需要拿起一个物品时，我们希望知道对应脑电波的反应，即该脑电波成份。该工作具有深远意义，如果能分离出与行为相关联的脑电波，将有助于对大脑疾病的诊断及脑中风病人的生活自理。由于人的大脑非常复杂，研究人的思维也相对困难。为了容易建立脑电波与行为之间的关系模型，这里我们选取的研究对象为小白鼠。我们感兴趣的是：

(1)视觉感受区的局部电位是否有规律性的变化？与呼吸曲线的周期性之间是否有联系？(2)视觉感受区的局部电位是否与视觉刺激相关？具体的联系是什么？

为了建立模型，我们提供的数据包括两方面：（1）睡眠状态下；（2）清醒状态下：(a)没有视觉刺激；(b)有视觉刺激。

提供的数据如下：

(1)文件名称：LC01_20131204_Data10_V1_50s_70s_1kHz.mat 数据格式：Matlab(.mat;int16)采样频率：1000 Hz 数据长度：30秒

第1-5 道（Ch11, Ch12, Ch13, Ch14, Ch15）：小鼠大脑视觉皮层局部电位(睡眠状态)第6道（Ch17）：呼吸曲线[注：仅波峰和波谷点有生理意义，分别代表吸气和呼气末。呼 吸信号是通过热敏电极在鼻子附近记录的，主要表现是呼气时温度升高（曲线下降）。因为是 间接记录呼吸，因此只能用曲线的峰或 谷的时间点标记呼吸时程（吸气相或呼气相），曲线 具体幅度等没什么价值。](2)文件名称：V01_20131126_Data03_80s_110s_LFP_NoStim_1khz.mat 数据格式：Matlab(.mat;int16)采样频率：1000 Hz 数据长度：30秒 第1-5 道（Ch11, Ch12, Ch13, Ch14, Ch15）：小鼠大脑视觉皮层局部电位(清醒状态)第6道（Ch17）：呼吸曲线

(3)文件名称：V01_20131126_Data03_300s_330s_LFP_VisStim_1khz.mat 数据格式：Matlab(.mat;int16)采样频率：1000 Hz 数据长度：30秒

第1-5 道（Ch11, Ch12, Ch13, Ch14, Ch15）：小鼠大脑视觉皮层局部电位(清醒状态)第6道（Ch17）：呼吸曲线

第7道（Ch18）：两次CheckBoard视觉刺激（注：波峰处为黑白配对视觉刺激开始。具体实 现见附图说明：moveonoff.pdf）信号采集条件说明：

信号（LFP：局部电位）分别来自睡眠或清醒状态下小鼠大脑皮层的视觉感受区，5个电极（一排）同时记录，每个电极间距为0.25 毫米。（电极置于皮下5毫米处，这样可以减少肌电的影响）。

另外，为了研究视觉刺激对LFP的影响。在记录视觉感受区局部电位的同时，动物的眼睛给了视觉刺激。该刺激是通过配对出现的Checkboard，开关的相隔时间是随机的（在小鼠的眼前放置一个Checkboard,随机地打开或关闭）。

请研究如下问题（部分问题给出了猜想，但问题的解决不局限于这些猜想）：

（1）由于对呼吸的观测是间接的，能否通过分析呼吸的机理，建立数学模型反映小鼠在睡眠状态下与呼吸相关联的脑电波。（猜测：呼吸过程是由脑干部分发出“呼”和“吸”的命令，由神经元集群同步产生动作电位，该电位完成呼吸过程。）（2）一般认为：在睡眠状态下，小鼠脑电波的周期节律有可能与呼吸相关联[2]。该结论是否正确？通过对所给的视觉感受区的局部电位数据建立模型论证你的结论。（3）研究在清醒状态下，小鼠视觉感受区的局部电位信号是否有周期性的变化？该周期性的变化是否与小鼠呼吸所对应的脑电波的周期性的变化有关？是线性相关吗？如果不是线性相关，是否具有其他形式的相关性？（4）建立脑电波信号的分离模型，能否从前两种状态（睡眠状态、无视觉刺激的清醒状态）的局部脑电位信号LC01_20131204_Data10_V1_50s_70s_1kHz.mat、V01_20131126_Data03_80s_110s_LFP_NoStim_1khz.mat中分离出与小鼠呼吸相关联的脑电波信号？

（5）通过Checkboard随时间变化的曲线，分析小鼠视觉刺激的时间曲线及其功率谱，并与呼吸曲线的功率谱对比。利用问题(4)所建立的信号分离模型，从数据V01_20131126_Data03_300s_330s_LFP_VisStim_1khz.mat中分离出与Checkboard刺激相关的脑电波信号成份和可能与呼吸相关的脑电波信号成份？请验证所分离出来的刺激脑电波成份与视觉刺激之间的相关性，并说明该脑电波成份中是否包含图形形状因素。

参考文献：

[1] J.Wessberg, C.R.Stambaugh, J.D.Kralik, P.D.Beck, M.Laubach, J.K.Chapin, J.Kim, J.Biggs, M.A.Srinivasan, M.A.L.Nicolelis, Real-time prediction of hand trajectory by ensembles of cortical neurons in primates, Nature,408, 361-365, Nov.2000 [2].J.Ito, S.Roy, Y.Liu, M.Fletcher, L.Lu, J.D.Boughter,S.Grun, D.H.Heck, Whisker barrel cortex delta oscillations and gamma power in the awake mouse are linked to respiration, Nature Communication, 2014, April.[3]Mark E.Bear等，王建军等译，神经科学—探索脑，高等教育出版社，2004 [4]王青云，石霞，陆启韶，神经元耦合系统的同步动力学，科学出版社，2008

附录：“关于大脑”

关于大脑

1.脑电波

现代科学研究已经发现人脑工作时会产生自己的脑电波。经过研究证实大脑存在至少多个不同波段的脑电波。事实上，脑电波是一系列自发的有节律的神经电活动，其频率变动范围在每秒1－30次之间的，可划分为四个波段，即δ（1－3Hz）、θ（4－7Hz）、α（8－13Hz）、β（14－30Hz）。（这几种波的频率边界，学界还没有完全统一的标准。亦有学者认为δ波小于4Hz，θ波4～7Hz，α波8～12Hz，β波13～35Hz，并认为有大于35Hz的脑电波，并命名为γ波。长期处于该状态下的人会有生命危险。另外，这几种波的划分是针对人类的。对于其它动物频率边界会有变化）δ波，频率为每秒1－3次，当人在婴儿期或智力发育不成熟、成年人在极度疲劳和昏睡状态下，可出现这种波段。

θ波，频率为每秒4－7次，成年人在意愿受到挫折和抑郁时以及精神病患者这种波极为显著。但此波为少年（10－17岁）的脑电图中的主要成分。

α波，频率为每秒8－13次，平均数为10次左右，它是正常人脑电波的基本节律，如果没有外加的刺激，其频率是相当恒定的。人在清醒、安静并闭眼时该节律最为明显，睁开眼睛或接受其它刺激时，α波即刻消失。

β波，频率为每秒14－30次，当精神紧张和情绪激动或亢奋时出现此波，当人从睡梦中惊醒时，原来的慢波节律可立即被该节律所替代。

在人心情愉悦或静思冥想时，一直兴奋的β波、δ波或θ波此刻弱了下来，α波相 对来说得到了强化。因为这种波形最接近右脑的脑电生物节律，于是人的灵感状态就出现了。

2.人脑的构造

人脑的构造主要包括脑干、小脑与前脑三部分。

脑干（brainstem）上承大脑半球，下连脊髓，呈不规则的柱状形。经由脊髓传至脑的神经冲动，呈交叉方式进入：来自脊髓右边的冲动，先传至脑干的左边，然后再送入大脑；来自脊髓左边者，先送入脑干的右边，再传到大脑。脑干的功能主要是维持个体生命，包括心跳、呼吸、消化、体温、睡眠等重要生理功能，均与脑干的功能有关。脑干部位又包括以下四个重要构造：

1．延髓（medulla）延髓居于脑的最下部，与脊髓相连；其主要功能为控制呼吸、心跳、消化等。

2．脑桥（pons）脑桥位于中脑与延脑之间。脑桥的白质神经纤维，通到小脑皮质，可将神经冲动自小脑一半球传至另一半球，使之发挥协调身体两侧肌肉活动的功能。

3．中脑（midbrain）中脑位于脑桥之上，恰好是整个脑的中点。中脑是视觉与听觉的反射中枢，凡是瞳孔、眼球、肌肉等活动，均受中脑的控制。

4．网状系统（reticular system）网状系统居于脑干的中央，是由许多错综复杂的神经元集合而成的网状结构。网状系统的主要功能是控制觉醒、注意、睡眠等不同层次的意识状态。

小脑（cerebellum）位于大脑及枕叶的下方，恰在脑干的后面，是脑的第二大部分。小脑由左右两个半球所构成，且灰质在外部，白质在内部。在功能方面，小脑和大脑皮层运动去共同控制肌肉的运动，籍以调节姿势与身体的平衡。

前脑（forebrain）属于脑的最高层部分，是人脑中最复杂、最重要的神经中枢。前脑又分为视丘、下视丘、边缘系统、大脑皮质四部分。

1．视丘（thalamus）视丘呈卵圆形，由白质神经纤维构成，左右各一，位于骈胝体的下方。从脊髓、脑干、小脑传导来的神经冲动，都先终止于视丘，经视丘在传送至大脑皮质的相关区域。所以说视丘是感觉神经的重要传递站。此外，视丘还具有控制情绪的功能。

2．下视丘（hypothalamus）下视丘位于视丘之下，是自主神经系统的主要管制中枢，它直接与大脑中各区相连接，又与脑垂体及延髓相连。下视丘的主要功能是管制内分泌系统、维持新陈代谢正常、调节体温，并与生理活动中饥饿、渴、性等生理性动机有密切的关系。

3．边缘系统（limbic system）边缘系统一般认为包括视丘、下视丘以及中脑等在内的部分。边缘系统的主要功能为嗅觉、内脏、自主神经、内分泌、性、摄食、学习、记忆等。边缘系统有两个神经组织，即杏仁核与海马，前者关系情绪的表现，后者与记忆有关。

4．大脑皮质（cerebral cortex）是大脑的表层，由灰质构成，其厚度约为1到4mm，其下方大部分则由白质构成。大脑中间有一裂沟（大脑纵裂，longitudinal fissure），由前至后将大脑分为左右两个半球，称为大脑半球（cerebral hemisphere）。两个半球之间，由胼胝体（corpus collosum）连接在一起，使两半球的神经传导得以互通。

大脑皮质为中枢神经系统的最高级中枢，各皮质的功能复杂，不仅与躯体的各种感觉和运动有关，也与语言、文字等密切相关。根据大脑皮质的细胞成分、排列、构筑等特点，将皮质分为若干区。皮质运动区：位于中央前回（4区），是支配对侧躯体随意运动的中枢。它主要接受来自对侧骨骼肌、肌腱和关节的本体感觉冲动，以感受身体的位置、姿势和运动感觉，并发出纤维，即锥体束控制对侧骨骼肌的随意运动。返回

皮质运动前区：位于中央前回之前（6区），为锥体外系皮质区。它发出纤维至丘脑、基底神经节、红核、黑质等。与联合运动和姿势动作协调有关，也具有植物神经皮质中枢的部分功能。

皮质眼球运动区：位于额叶的8区和枕叶19区，为眼球运动同向凝视中枢，管理两眼球同时向对侧注视。

皮质一般感觉区：位于中央后回（1、2、3区），接受身体对侧的痛、温、触和本体感觉冲动，并形成相应的感觉。顶上小叶（5、7）为精细触觉和实体觉的皮质区。

额叶联合区：为额叶前部的9、10、11区，与智力和精神活动有密切关系。

视觉皮质区：在枕叶的距状裂上、下唇与楔叶、舌回的相邻区（17区）。每一侧的上述区域皮质都接受来自两眼对侧视野的视觉冲动，并形成视觉。返回

听觉皮区：位于颞横回中部（41、42区），又称Heschl氏回。每侧皮质均按来自双耳的听觉冲动产生听觉。

嗅觉皮质区：位于嗅区、钩回和海马回的前部（25、28、34）和35区的大部分）。每侧皮质均接受双侧嗅神经传入的冲动。

内脏皮质区：该区定位不太集中，主要分布在扣带回前部、颞叶前部、眶回后部、岛叶、海马及海马钩回等区域。

语言运用中枢：人类的语言及使用工具等特殊活动在一侧皮层上也有较集中的代表区（优势半球），也称为语言运用中枢。它们分别是：①运动语言中枢：位于额下回后部（44、45区，又称Broca区）。②听觉语言中枢：位于颞上回42、22区皮质，该区具有能够听到声音并将声音理解成语言的一系列过程的功能。③视觉语言中枢：位于顶下小叶的角回，即39区。该区具有理解看到的符号和文字意义的功能。④运用中枢：位于顶下小叶的缘上回，即40区。此区主管精细的协调功能。⑤书写中枢：位于额中回后部8、6区，即中央前回手区的前方。

3.脑的节律

地球是一个富有节律的环境，气温和白昼都随着季节的变化而变化，每天昼夜交替，潮涨潮落。为了生存，动物的行为必须和环境的节律协调。在进化过程中，脑内形成了各种各样的节律控制系统。睡眠和觉醒是最明显的节律行为。由脑控制的节律周期有些比较长，比如冬眠动物的一些节律活动；而另一些较短，如呼吸节律，行走步伐和大脑皮层的电节律等。

大脑皮层存在着大量的电节律。它们通常是大群神经元组同步产生的。比如，当你试图抓住一个篮球时，同时对形状、颜色、运动、距离甚至篮球的含义反应的不同神经元群组趋于同步的震荡。通过分散的细胞群组震荡的高度同步化，在某种程度上把它们 连成一个有意义的大组，以区别于附近的神经元，从而把“篮球问题”这个杂乱的神经碎片统一起来。

迄今为止，大脑皮层节律的功能在很大程度上还是一个谜。

第五篇：2013年全国研究生数学建模竞赛F题

可持续的中国城乡居民养老保险体系的数学模型研究

一、问题的提出

中国共产党第十八次全国代表大会政治报告中提出了“统筹推进城乡社会保障体系建设”的任务：“社会保障是保障人民生活、调节社会分配的一项基本制度。要坚持全覆盖、保基本、多层次、可持续方针，以增强公平性、适应流动性、保证可持续性为重点，全面建成覆盖城乡居民的社会保障体系”。

近年来，我国城镇职工养老保险金制度基本建立，参保覆盖面日益扩大，退休人员待遇逐年提高，养老保险基金规模不断上升，在制度层面上提前完成了城市和农村居民社会保障全覆盖，成为我们社会保障事业发展的重要里程碑。然而也面临巨大的挑战。我国的养老保险体系建设远远滞后于社会发展水平。一方面资金规模相对较少，保障能力有限，地区差距、城乡差距明显；另一方面资金管理分散，运营效率不高，保值增值的压力很大。进一步完善养老保险保障体系仍将是中国政府未来最重要的任务之一。

全国人大常委会法工委信春鹰副主任是《社会保险法》的制定者，她指出中国养老制度“最大的问题就是‘可持续’”，“我们看当下的欧洲，我们看现在出现财政危机的那些国家，他们的养老保险无一不遇到了极大的挑战和问题。和其他国家相比，中国有特殊的问题，像前面几位专家讲的，我们比人家穷，比人家人口多，‘未富先老’。从制度上来说，我们的制度本身就有很多的缺陷，本身就很青涩，它是对现实妥协的一个产物。但是同时，这个制度一经产生就遇到了人口结构快速变化，社会尚未准备好，这些大的问题都出现了，对我们这个制度是非常大的挑战。” 由于中国养老保险缴纳与支付同时开始，中国保监会副主席陈文辉在陆家嘴论坛上承认“中国养老金缺口确实非常大”。社会科学院世界社保研究中心主任郑秉文给出的数字是2011年城镇基本养老保险个人账户“空账”已经超过2.2万亿，较2010年增加约5000亿。世界银行的研究机构在研究我国的养老金问题研究报告中，认为到2013年，中国养老基金缺口达到18.3万亿元。在当前养老制度不变情况下，以后的年份缺口逐年放大。假设2013年后中国GDP的年增长率为6%，到2033年时养老金缺口将为惊人的68.2万亿元，占当年GDP的38.7%（见附件3、4，由于这是根据国家资产负债表估算，显然比较粗糙）。另据中国社科院财贸所学者高培勇和汪德华最近的一项研究指出，在现有社保制度框架下，要确保2020年中国每个退休者都能领取养老金，以替代率52.4%计算，养老金支付缺口都会持续增加。最早在2015年，中国城镇职工养老保险就会出现支付缺口。两位学者预计，如果现在不能未雨绸缪，及时调整政策，到2050年，养老金缺口累积将达到该年GDP的95%。这样巨大的资金缺口，必将使国家财政不堪重负。

而现有官方公布的数据却有明显差别。2012年6月4日，中国人力资源与社会保障部向社会公布了2011人力资源和社会保障事业发展统计公报。此公报是透视中国人力资源与社会保障工作成绩的重要窗口，也是解析社会保险最权威的数据来源。统计公报显示，截至2011年末，全国基本养老保险基金累计结余1.9497万亿元，全年城镇基本养老保险基金总收入1.6895万亿元，比上年增长25.9％，全年基金总支出1.2765万亿元，比上年增长20.9％。和2009年、2010年相比，城镇职工基本养老保险基金当期结余资金规模进一步扩大，已由2010年当期结余的2865亿元增加到4130亿元，增幅较大，这有利于缓解部分

省市养老金支付压力。但有分析认为：基金征缴收入增长较快是其中一个重要原因。上述观点分别反映中国养老金现状的不同侧面和未来趋势的不确定性。

二、现实的情况

中国在20世纪90年代确立了城镇职工基本养老保险“统账结合”制度，基本养老保险是由社会统筹和个人账户两部分组成。社会统筹部分由单位负担缴费，为单位职工工资总额的20%，个人账户则由职工个人缴费，为个人工资的8%。前者“现收现付”，用于支付已退休人员的养老金，后者实行长期封闭积累，产权个人所有的“完全积累”制，原则上不能调剂借用的方式。

由于在现行养老保险制度确立之前，企业员工基本无需缴纳养老保险费用，所以现有的养老保险基金中没有这部分职工的个人账户部分。在养老保险制度设立之后，这部分职工退休后却要从养老保险基金中领取养老金，这是我国社会保险制度必须承担的转轨成本的一部分。仅靠统筹账户不足以应对当期发放，加之各地财政实力不同，多数地区不得不在实际上采用了“现收现付制”的方法，即挪用个人账户的资金，用正在工作人合计缴纳的28%的月工资来支付现有退休人员的退休金，现在工作人账户仅仅用来记账，上述2.5万亿的记账额由此形成。为解决这一问题，中国从2000年开始了“做实”个人账户试点。截至到2011年底，参与试点的辽宁、江苏和山东等13个省份共积累个人账户基金达到2703亿元，但其与记账额之间的差额，仍达到2.23万亿元，此“空账”被舆论定义成了养老保险“缺口”概念而广为传播。

三、需要解决的问题

面临人口老龄化以及经济结构性减速、财政收入增长速度下降，社会保障体系可持续问题亟待解决，必须进行前瞻性研究，否则将形成社会问题。我们应该承认历史、立足现实、尊重规律、借鉴国外、留有余地，既量力而行，又尽力而为地完善我国的养老保险体系。本题要求运用数学模型方法研究、解决该问题。为此，请你们团队做好以下几项工作。（请研究生注意：可以引用有关参考文献的结论和数学模型，更鼓励创新，独立提出新模型。若引用他人的模型与结论，必须注明出处，否则视情节轻重，直至作为抄袭处理）

1，分别建立合乎国情、适应国力的中国城乡居民（含新农保）养老金收入、支出的宏观数学模型，至少包括替代率（基本养老保险人均养老金占城镇单位在岗职工平均工资比率）、缴费率（基本养老保险人均缴费占城镇单位在岗职工平均工资比率）、人口结构、分年龄段死亡率、经济增速、财政补贴、工资水平或物价指数、投资效益等主要因素，要做到模型结合现实，分多个层次（含企业基金等），体现“多缴多得，长缴多得”（不考虑分省、分地区模型）。如果使用国家统计局公布和题目提供数据以外的数据，必须指明出处并向竞赛网站申请张贴于“资料下载”处，首先公布者可适当加分。

2，根据你们的数学模型、对养老金缺口的理解和对未来有关情况的合理估计，估计从今年至2035年我国养老金缺口，并说明你们对养老金缺口分析的合理性。如果全部情况维持不变，按照你们的数学模型我国城乡居民养老保险收支矛盾最尖锐的情况发生在什么时间，严重程度如何？考虑到党的十八大提出的收入倍增计划，你们的数学模型哪些部分需要调整？

3，养老保险制度也是调节社会分配，请你们分析各国养老保险的不同模式，取其精华，去其糟粕，根据你们建立的数学模型和中国的实际情况，利用仿真手段寻找替代率和缴费率的合理区间以保证我国养老保险体系的可持续性（因为人口结构、分年龄段死亡率、经济增速、投资效益等主要因素几乎无法人为较大幅度改动）；在步入良性循环之前，在矛盾最尖锐到来前的过渡期内应该采取哪些政策措施实现平稳过渡并仿真预测相关政策的效果。

四．尝试建立第三问增加可调节变量的数学模型。

参考文献

1.《中国养老金发展报告2012》中国社会科学院世界社保研究中心(其中部分全国性数据扫描在附件中，并提供word等版本省去输入)2.中国新闻网

《中国养老金发展报告2012》在京发布式 3.马骏等《化解国家资产负债中长期风险》

4.曹远征《意见中国――网易经济学家访谈录》