大学生数学建模竞赛-航班延误问题研究.docx

航班延误问题研究摘要随着我国民航运输业的快速发展，我国民航业正在经历美国等西方发达国家在上个世纪七十年代放松航空管制后航空业务量急剧膨胀，导致机场和空域拥堵严重、航班延误快速增长的局面。引起的乘客和航空公司之间纠纷也逐渐增多，如果不能及时解决，会激发两者之间的矛盾，从而影响航空公司的声誉。面对目前的形势，我国民航局在2010年开展了航班正常性和大面积航班延误的专项整治工作，要求各民航各企事业单位采取措施努力做好航班延误的治理工作，降低航班延误对我国民航运输业的不利影响。本文基于收集得到的数据，分析国内航班延误的真实原因，并对航空公司及乘客如何应对航班延误提出合理的策略，构建模型并对其前景进行分析，最后，本文基于航班总数的时间序列数据，对未来十年民航市场的发展趋势做出适当预测。针对问题一，我们将通过各种渠道了解国内各航空公司航班记录，特别是延误的航班，以及后续措施的实施。记录实情，通过分析证明结论的正确性。针对问题二，我们首先对原始数据进行统计并处理，得到航班总数，正常航班数，不正常航班数的时间序列数据，而且，在次基础之上，对因各种因素导致的航班延误数进行统计分析，充分挖掘航班延误的几个主要原因是航空公司自身原因，流量原因和天气原因。同时，得到各个年份的导致航班延误影响因素的比例分布表，紧接着做出这个比例分布表的直方图，进而依据数据特征并结合现实具体情况来分析航班延误的四个主要影响因素，最后我们得出结论：日益增长的航空运输需求与有限的空域资源之间的矛盾是航班延误的主导原因。针对问题三， 我们从航班延误成本最小和航班延误时间最短两个点入手，构造动态规划模型，最后利用匈牙利算法，为航空公司在航班延误上提供了合理的管理措施，同时针对航班延误的变化规律也为乘客做出了合理的出行建议。为了对未来十年国内民航市场的发展趋势做适当的预测，我们考虑从航班总数的角度入手，做出散点图，分析其发展趋势，并基于MATLAB软件用多项式拟合的方法得到一条拟合曲线，经检验，一次拟合的效果比较好，所以，我们在此基础之上得到未来十年的航班总数的估计值，并得出结论，未来十年国内民航行业的发展具有非常广阔的潜力。关键字：统计 航班延误 延误成本 拟合一、问题重述1、中国的航班延误最严重，国际上航班延误最严重的10个机场中，中国占了7个的结论是否正确？2、我国航班延误的主要原因是什么？3、有什么改进措施？二、问题分析2.1问题一的分析问题一要求统计国内国际航班延误数据，进行合理处理，证明结论是否正确。首先，我们查阅国内外各大航空公司的网页和一些主要统计部门的相关信息，得到关于年度航班延误的一些统计指标，并在此基础之上，考虑利用MATLAB软件做出各种统计指标的散点图，对各航空公司的实际情况分析，证明结论的正确性。2.2问题二的分析问题二要求我们分析我国航班延误的主要原因是什么。显然，航班延误是当前国际民航业发展中的一大难题，也是顾客对航空服务质量不满意的主要内容。根据收集得到的数据，我们发现，导致航班延误有两大主要原因，一是航空公司自身的原因，涉及到航空公司自身的相关运行管理；另外一方面是非航空公司自身因素，即空管流量控制，恶劣天气，军事活动等非航空公司自身因素。为了问题分析的方便，考虑对数据进行更深层次的挖掘和处理，并且，有效结合实际情况，分析得出航班延误的真实原因。2.3问题三的分析问题三要求我们根据数据分析提出改进的措施（如航空公司的预定票策略，乘客购买航空延误保险或恰当选择出行方式等），我们通过分析历年我国航班延误率初步得出我国延误的大致水平，然后从航班延误成本和航班延误时长两个点入手，构造动态规划模型，最后为航空公司提供了一种合理的管理措施，即在延误时长一定的合理范围内，满足延误成本最小的建议。同时我们通过分析航班延误率和延误时长的发展规律，给乘坐飞机的乘客提出了几种合理的意见，如周六航班延误时间较长且延误的可能性更大，对于此种风险厌恶系数较大的乘客不建议在周六出行等。三、问题假设1、假设收集到的数据真实可靠；2、假设各航空公司能制定切实的政策和措施；四、符号定义与说明i飞机的指示j航班的下标af执行航班f的飞机af替换航班f的飞机Ti可用飞机的就绪时间集合TJ最早延误航班之后的航班按原计划到达时间集合F最早延误航班之后的航班集合A最早延误航班之后可用的飞机集合AIm能够在m机场维修的机型为I的飞机集合Z当天备用飞机和修复飞机的集合xijf时间对i到j的航班yf取消航班f的标志，1为取消，0为不取消Pf旅客的失望溢出成本v乘客数w该航班上的平均票价pijI时刻就绪的飞机执行j时刻的航班及后续航班的延误成本xbj1表示当天有可用飞机b指派给航班f，0表示没有cfb把航班f指派给备用飞机的成本nf航班f在时间对i和j之间经过的机场数五、模型的建立与求解5.1 问题一的分析下图是2015年5月2日北京南苑机场CA1835一天的航班飞行情况图1.1从今日CA1835航班起飞状况，共23个时间点，航班延误19个，占82.6%。从这可以看出国内航班延误的严重性。图1.22015年5月2号全国航班延误情况图示说明：表示机场正常表示机场部分航班延误表示机场大面积延误表示机场当前没有航班起降表示当前机场状态无法判断图示等级：大型机场 中型机场 小型机场从全国来看，航班延误的机场主要分布在北上广地区，银川和台北桃园机场出现大面积航班延误，因此全国航班延误情况依旧值得关注。5.2问题二的分析与处理航班延误是指航班降落时间比计划降落时间（航班时刻表上的时间）延迟30分钟以上或航班取消的情况。日常生活中航班延误不仅影响着乘客的心情，也影响着航空公司的运行效率和服务质量，所以我们一般使用准点率来衡量承运人运输效率和运输质量。准点率，又称正点率、航班正常率，是指航空旅客运输部门在执行运输计划时，航班实际出发时间与计划出发时间的较为一致的航班数量与全部航班数量的比率下表2-1是我国2006年至2011年的年度航班延误统计情况，可以看出航班的正常率普遍高于80%，而这一数据明显高于43家国际主要航空公司的航班平均准点率76.54%1：表2-1年度航班延误统计情况类别时间航班数正常航班数不正常航班数正常率航空公司原因流量 原因天气 原因其他20061530443 1254258 276185 81.95%117711 57570 75797 25107 20071613786 1331955 281831 82.54%126374 58741 79937 16778 20081528208 1274090 254140 83.37%116842 58516 59398 19384 20091759438 1437036 322601 81.68%135921 72544 75676 38460 20102010652 1617150 403511 80.43%163821 105611 78802 55278 20112204147 1861196 343050 84.44%128426 93911 61255 59458 图2.1是我国06-11年航班数的发展情况，以及不正常航班数的变化趋势，初步分析可以得出，随着我国经济的发展，飞机作为一种交通工具越来越普遍，而需求的增加势必引起供给的增加，但是航班数的增加，所带来的航班延误也同比小幅上升，而且居高不下，这确实需要航空公司的进一步合理规划。图2.1 06-11年我国航班变化情况一般来说，航班的延误主要有以下原因：1、 航空公司的运行管理2、 流量控制3、 恶劣天气影响 4、 军事活动影响5、 机场保障其中军事活动和机场保障所造成的航班延误概率较小，为方便分析，我们将这两类归为其他原因。下图2.2为四种原因的变化趋势图，为更好地观察变化，我们取半年为一个观测点，时间范围为2006-2007年。用airlines航空公司原因，用flow表示流量控制，用weather表示天气原因，用other表示其他原因，纵坐标表示四种原因的所造成的延误数。图2.2 各航班延误原因的变化趋势图观察上图可以看出，由于航空公司自身原因所造成的延误在过去几年一直都是维持在6000（件/半年）以上，且教稳定，而在2010年的时候波动较大。流量管制则在10年以前稳定在30000（件/半年）左右，且10年变化波动突然上升。天气原因则在4000（件/半年）波动，其他原因也一直维持在较少的次数。从上图2.2我们可以看出过去几年航班延误的各种原因的变化情况，为了进一步看出各中原因所占的比重，我们通过加总计算过去几年各种原因下航班延误发生次数的和，再计算其百分比，画出其饼状图，如下图2.3所示：图2.3 各航班延误原因占比图由上图可以看出在航班延误原因中由于航空公司自身原因所造成的原因占最大的比重，占比42.17%，而天气原因和流量管制所造成的航班延误则差不多，约为23%，其他原因所占的比重比较小，占比10.87%。航班延误是当前国际民航业发展中的一大难题，也是顾客对航空服务质量不满意的主要内容。从上面的分析中，我们可知航班延误的主要原因有：一、航空公司的运行管理；二、流量控制；三、恶劣天气影响；四、其他。其中军事活动和机场保障是比例比较小的，所以我们为了问题分析的方便所考虑将这两者归结为其他。经过处理后的数据如下表2-4所示。表2-4 航班延误影响因素比例结构表年份航空公司流量控制天气其他20060.480.220.230.0720070.470.280.150.1020080.430.190.270.1120090.390.230.190.1920100.410.240.230.1220110.370.280.200.1520120.360.220.210.21资料来源:2006-2013从统计看民航将表2-4中的数据以直方图的形式呈现，如下图2-1所示。图2-4 航班延误原因直方图由直方图我们可以清晰的看出，在航班延误影响的因素比例中，航空公司自身的影响是占比重最大的，但从2010年以来，这个比例在逐年下降；天气原因造成的航班延误基本保持在20%左右。从当前实际来看，导致航班延误的原因可以分成两大类，分别为航空公司自身因素，例如不合理的航班调配；另外一类为非航空公司因素，例如流量控制，天气原因，军事活动等等。在上述归类的四大原因中，除天气原因外，其他三方面原因只是航班延误的表层原因，并不是航班延误的深层次原因和实质性矛盾。表面看来，航空公司自身因素是航班延误的“罪魁祸首”，因为数据表明，其所占比重为40%左右，但由于航空运输的系统性，航班能否正常准点起飞，很大程度上取决于民航系统中其他相关单位的协调与配合，例如机场和空中交通管理部门，而且，目前的航班延误的统计也存在一定问题，致使一些不是航空公司自身原因导致的航班延误也计入航空公司自身因素里，例如空中交通管理部门实施的流量控制也会导致航班延误。由此可以得出导致航班延误的真正原因是：随着国家经济社会的发展和改革开放的深入，中国航空运输的需求量日益增加，而民航可使用的空域仅占中国全部空域的20%左右，大量空域被划为军航空域或者禁区，日益增加的需求量与优先使用的空域资源之间的矛盾是导致航班延误的真实原因。有数据显示：2011年中国人均乘机次数是0.2次，比2002年的0.07次增长了3倍，比1978年提高了100倍。然而改革开放以来，我国民用航空的空域资源一直被限制在20%左右，时至今天，两者之间的矛盾越来越恶化，这才是航班延误的真实原因。5.3模型的建立和求解航班延误问题的处理一直是航空公司的比较棘手的一件事，也是国际航空行业的一个痼疾，而目前我国针对航空延误的措施虽不断地在改进，如成立航班延误治理委员会，建立预警系统和取消航班时刻措施，在一定程度上减小了航班的延误率，但仍是收效甚微，其中1998至2008年的延误率如下表3-1：表3-1 98-08年我国航班延误率情况年份19981999200020012002200320042005200620072008延误率22.9%23.8%24.1%23.4%27%20.2%20.1%19.9%18.41%16.88%17.43%数据来源：中国民用航空局网可以看出我过航班的延误率大体在20%左右，波动较小。5.4乘客应对延误策略以及航空公司的改进措施：5.4.1乘客应对延误策略上述模型针对的是航空公司应对航班延误的策略模型，而乘客如何应对航班延误，同样仍是一个值得深究的问题，下面我们将通过分析航班的延误规律，为乘客提供一些参考的意见，下表3-1是我国15家航空公司一周内的日均航时，和平均延误时长： 表3-1 航空日均航时和延误时 单位：min类别周一周二周三周四周五周六周日航时14847129201884313712138593602014706平均延误时长45394538378041据上表做如下图3.1、图3.2：由图3.1可以看出日均航时在周六出现一个高峰，相比与其他工作日和周日，周六选择航班出行的乘客相对比较多，而在航班客座供给一定的情况下，势必会对航空公司的航行造成一定的压力，而这种压力恰好体现在了下面的图3.2.图3.1 日均航时 由图3.2可以看出日均延误时长在周六出现一个较大的向上波动，这也正是航流人数增多给航空公司造成压力的一个体现。 图3.2 日均延误时长综合上述可以初步得出，乘客若选择周六乘坐航班出行，遇上航班延误的可能性会增大，另外由于航班延误造成的延误时长也会偏长。故建议乘客尽量少选择航班时期在周六的。图3.3 月均延误率由上图3.3 2006年2月至2011年5月的月均延误率可以可以看出过去几年来在每年的5-8月，航班的延误率总会小幅走高，为20%以上。乘客可以根据上图合理安排自己的出行时间，尽量选择在航班延误概率较大的时段出行。5.4.2航空公司自身具体原因分析（一）跑道系统容量的限制枢纽机场运作的特点之一就是航班波峰时刻的形成，也就是航班在几个时段内集中抵达和集中离开。在航班的波峰时刻，由于受到机场系统处理能力的限制特别是跑道系统处理能力的限制，航班的地面排队现象往往比较严重。根据美国运输部统计局（Bureau of Transportation Statistics，简称“BTS”）的统计数据显示，2008年1月至6月美国部分机场的航班地面延误和地面滑行时间如下表：机场名称 （机场IATA代码） 到达 出发 架次 平均到达延误时间（分钟） 平均进港滑行时间（分钟） 架次 平均出发延误时间（分钟） 平均出港滑行时间（分钟） 亚特兰大(ALT) 171,5228.9610.37171,59011.2319.66芝加哥奥黑尔(ORD) 148,76820.078.96148,75520.5219.50达拉斯(DFW) 117,65811.0410.47117,60713.4916.73丹佛(DEN) 99,1827.288.2199,13910.2016.25休斯顿(IAH) 78,5029.018.5578,51711.4719.25旧金山(SFO) 57,84415.996.1757,85214.5517.67奥兰多(MCO) 57,1377.816.8257,14610.8214.32纽约拉瓜迪(LGA) 51,10918.468.3551,10812.3728.94纽约肯尼迪(JFK) 49,5019.779.7449,49810.8632.22华盛顿巴尔迪摩(BWI) 43,5605.905.1143,5599.2712.27迈阿密(MIA) 27,29613.977.8127,28015.7517.30从表中可以看这些机场的航班到达和出发的地面滑行时间都很长，如纽约的拉瓜迪和肯尼迪这两个机场出港航班的平均地面滑行时间分别达到了28.94分钟和32.22分钟。航班地面滑行时间大幅增加除了受到机场跑滑系统物理构形的影响外，主要是受跑道系统容量的限制，在波峰时刻出港航班的地面等待时间很长。（二）地面保障能力不足航班在抵达机场和离开机场的一个地面经停周期，有一系列的保障作业活动需要开展，如：上下客、加油、排污、客舱清洁、加水、配餐以及其他相关的保障作业活动。航班从着陆、滑行、进位、机位上保障作业、推出、滑行、起飞所经历的每一个环节都受到保障能力的限制，当航班的运行量超出任何一个节点的保障能力时，必然导致航班排队等待现象的形成。排队一旦形成，将使后续航班产生持续性的延误。由于航班的排队等待延误具有“波浪效应”（即传递效应），因此，当机场在一天之中连续几个时段内的航班运行量超出机场地面处理能力时，航班延误的现象将会变得非常严重。（三）机场设施设备故障机场的离港系统、行李处理系统、供油系统、供电系统、信息系统对于机场的正常运行保障非常关键，一旦这些关键系统发生故障，必然影响机场对航班和旅客的保障作业，导致航班延误的产生。（四）机场不满足最低运行标准受大雾、雷暴等恶劣天气或者跑道及其他关键硬件设施受损因素的影响，机场不满足最低运行标准，航班无法正常滑行和起降，因此，将导致机场大面积的航班延误的发生。航空公司的改进措施在导致航班延误的各种原因中，恶劣天气、突发事件等因素是很难通过人为改变的，但是，其他如跑道系统容量限制、地面保障能力不足、设施设备故障等因素可以通过提高机场的运行管理水平得到有效的改善。按照航班延误治理的目标，可将航班延误治理的对策分为预防性措施和应急处置措施。预防性措施可以有效地降低机场航班延误的程度，而应急处置措施可以有效地降低延误产生后对机场的不利影响。因此，机场想要做好航班延误的治理工作，就应当同时做好航班延误的预防和航班延误的应急处置这两个方面的工作。（一）航班延误的预防措施1、做好机场的总体规划，准确评估机场的当前有效容量，准确预测机场航空业务量发展趋势，根据需求和容量的评估结果选择适当的时间节点扩建机场，提高机场系统的有效容量。根据美国联邦航空局（简称FAA）对机场容量和延误的研究结果，机场的航班延误程度与机场的容量和需求有着密切的对应关系，具体见下图：需求指数（即：需求量与容量的比值）在0.85-0.9之间时，由容量引起的延误是可以接受的，当需求指数大于0.9时，延误将呈指数方式增长，按照FAA的研究结果，平均每架次飞机的延误时间一般不应超过4分钟。因此，机场当局应当及时掌握当前机场系统的使用状况，根据航空业务量发展的预测结果，及时扩建机场设施，提高机场系统容量以满足航空业务量的发展需求。2、做好航班延误状况和延误原因的统计分析工作，详细记录航班通过机场地面各节点的时间，准确分析航班延误的具体原因，通过对航班延误原因的分析，采取有针对性的治理措施。由于我国机场的信息化水平较低，很多机场对于航班保障细节信息还是采用人工填写方式，并且这些保障细节信息没有进行汇总录入到机场数据库中，信息的集成度较低，导致机场对于航班保障各节点的时间信息缺失。没有这些信息，航班延误的统计分析工作就无法开展，机场也很难分析影响航班延误的具体原因，也无法掌握问题到底发生在那个环节，使得机场在治理航班延误的工作中无从下手。因此，提高机场信息化水平，做好航班延误的统计分析工作也是航班延误治理很重要的一项措施。3、提高机场的地面保障能力，根据航班波峰时刻的保障需求，合理配置机场的保障设施、设备、人员等保障资源，提高机场的运行效率、运行管理水平和地面保障能力，避免因机场保障能力不足而导致航班延误的产生。4、航班时刻的协调管理，民航局在审批各航空公司航班时刻申请时，主要参考的是机场的容量评估的结果。当航空公司未按照民航局批准的时刻执行航班时，将导致：其一，机场的航班波峰时刻更加明显，高峰期航班更加集中，受机场容量和地面保障能力的影响，航班延误将迅速增加；其二，按照民航局航班正常性统计标准，这些航班将被视为不正常，导致机场航班的总体正常性降低（民航局的统计结果显示，在11种不正常原因中，公司原因所占的比率最大）；其三，影响航班时刻的管理秩序。虽然我国的航班时刻管理采用的是行政管制制度，机场在航班时刻协调管理中的话语权不大，但是，机场可以积极的监督和评估各航空公司的时刻执行率，协调局方对时刻执行率低的公司采取相应的措施，降低因航空公司航班时刻执行率低对机场运行造成的不利影响。（二）航班延误的应急处置措施机场应当提高航班延误的应急处置水平，在航班延误特别是大面积航班延误发生后，能够迅速地采取有效的措施，妥善处理受影响的航班和旅客，尽可能将航班延误所造成的影响降到最低。1、建立航班延误的快速反应机制同一架飞机所执行的航班在每一个机场的周转周期内产生的延误具有累积效应，如不采取及时的弥补措施，每个周转周期内的延误累积后，航班的延误现象会变得非常严重。因此，机场应当建立航班延误的快速反应机制，当进港航班在前站已经产生延误的情况下，及时发布航班延误的预警信息，各保障单位重点加强这些延误航班的保障，协调空中交通管制部门对延误航班优先放行，尽可能缩短航班的过站时间，以减少航班的累积延误。2、制定和完善机场的大面积航班延误应急处置预案制定和完善机场的大面积航班延误应急处置预案，提高机场的应急反应速度和应急管理水平，有助于机场在面临大面积航班延误的情况下，能够高效、有序的开展工作，妥善处理因航班延误引起的各种问题，减小航班延误的影响程度和影响范围。（三）航班延误成本与延误时间的分析目前我国国内对航班延误的研究有很多，王红、刘金兰、曹卫东、郇秀霞（2009）利用Markov链模型，对航班的延误进行预测，再利用定性加定量的AHP层次分析，对航班的延误进行了预警处理，得出一种可以帮助航空公司管理延误的措施。而李俊生、丁建立（2008），刘玉洁（2009）等则是从航班的延误的波及入手，利用贝叶斯网络的传播模型进行分析，结果同样是得出了一种可以帮助航空公司管理延误的措施。关于航班延误的管理，国内研究已经颇多，但都由于过于复杂比较难实现，且其中关于延误成本的概念，较少被提及，而本文正是从该概念入手，通过建立一种延误成本最小的航班调度模型，既在一定程度可以帮助航空公司减少航班延误的发生，也帮其在航班延误发生的情况下使得损失成本最小。为了更好地分析问题，下面给出一些符号的定义：i飞机的指示j航班的下标af执行航班f的飞机af替换航班f的飞机Ti可用飞机的就绪时间集合TJ最早延误航班之后的航班按原计划到达时间集合F最早延误航班之后的航班集合A最早延误航班之后可用的飞机集合AIm能够在m机场维修的机型为I的飞机集合Z当天备用飞机和修复飞机的集合xijf时间对i到j的航班yf取消航班f的标志，1为取消，0为不取消Pf旅客的失望溢出成本v乘客数w该航班上的平均票价pijI时刻就绪的飞机执行j时刻的航班及后续航班的延误成本xbj1表示当天有可用飞机b指派给航班f，0表示没有cfb把航班f指派给备用飞机的成本nf航班f在时间对i和j之间经过的机场数其中延误成本：pij=Pf+bZcfbxbj （3.1）其中： Pf=wv3t/60229 （3.2）构建如下目标函数：minfFiTi,jTjpij or minfFiTi,jTjnft （3.3）约束条件：iAZxif+yf=1 fF （3.4）保证了每个时间对上都有航班覆盖。fFxif1 fA （3.5）保证每个航班都有飞机执行，否则取消航班。fFxifbfI=1 iAIm，IA （3.6）保证用于替换的飞机型号满足替换要求。其中： yf0，1，bfI0，1，xif0，1 为求解上面3.1式最优解，利用匈牙利科学家柯尼格提出的匈牙利矩阵算法，该算法的思想是系数矩阵(Pij)的一行(列)各元素中分别减去该行(列)的最小元素,得到新的(Pij)矩阵,那么以新的(Pij)为系数矩阵求得的最优解和用原系数矩阵求得的最优解相同。所以首先构造延误时间置换矩阵Tij：Tij=t11t1ntm1tmn，i1，2，n；j1，2，m 其中tij表示i时刻航班的飞机执行地j时刻航班的任务所延误的时间，根据延误时间置换矩阵Tij，计算延误成本置换矩阵Pij，Pij=p11p1npm1pmn其中Pij表示i时刻飞机执行第j时刻航班的任务的延误成本，最终可由上述矩阵得出航班置换方案，当然航班的置换最终还是要权衡两者的大小，单纯使得延误成本最小，势必使得延误时间不是最优，而使延误时间最优，又可能造成延误成本偏大，故在延误时间一个合理的范围内求解出延误成本最低，才是航空公司的最终目标。5.5未来十年国内民航市场的发展趋势显然，我们可以从不同的角度来分析未来十年民航市场的发展趋势并以此为依据做出适当的预测，根据所收集到的数据，我们考虑从年度航班总数的变化趋势来预测未来十年的航班总数的变化，从而再从航班总数的变化来说明未来十年的航班总数的变化趋势。首先，将收集得到的数据进行整理，得到2010年2014年的航班总数的数据，列表如下表5-1所示。表5-1 2010年-2014年航班总时间航班总数20105,531,71620115,979,66420126,603,20720137,315,44020147,933,110为了对航班总数的变化趋势做出进一步的分析研究，考虑用MATLAB软件做出其散点图，具体代码参照附录二，结果如下图5-1所示。图5-1 2010年-2014年航班总数散点图由图5-1，可以发现，航班总数普遍呈现递增的趋势，所以，我们考虑先对原始数据进行标准化处理，以便降低误差，再对得到的数据用MATLAB软件进行多项式拟合，具体代码参照附录三，拟合曲线见图5-2。你从网上搜搜怎么画航班总数对数拟合图，这个我还没弄好，你比照着源文件图。图5-2 2006年-2011年航班总数对数拟合图另外，假设航班总数为Y，时间序列为t,其中t=1，2，,6，根据MATLAB运行结果，我们可知该条拟合直线的方程为y=0.577t+5.2143还有公式写下，我没弄对当然，为了检验该模型的拟合程度，有必要进行一定的检验，下面，我们给出该对数航班总数时间序列拟合之后的残差图，如下图5-3所示。图5-3 对数航班总数的残差图还有这个图，我没改，你搜搜可对图5-3显示，该一次拟合式子是通过我们的检验的，所以，进一步我们认为航班总数和时间之间的关系式为： y=0.577t+5.2143 （5.1） 所以，我们可以根据式子5.1将未来2012年至2021年的航班总数估计出来，将估计结果整理出来，如下表5-2所示。表5-1 2010年-2014年航班总时间航班总数20105,531,71620115,979,66420126,603,20720137,315,44020147,933,110表5-2 2012年-2021年航班总数的估计值时间航班总数估计值20126,603,20720137,315,44020147,933,11020158,676,30020169,253,30020179,830,300201810,407,300201910,984,300202011,561,300202112,138,300为了进行更加直观的描述未来航班总数的趋势，考虑将2010年至2021年航班总数的散点图做出来，具体代码参照附录四，MATLAB运行结果如下图5-4所示。图5-4 2010年-2021年航班总数散点图根据图5-4，可以发现，航班总数从2010年到2021年的趋势是逐年递增的，即可以说明，未来的航空运输市场是一个比较大的领域，但是，由于目前中国航空领域经常发生航班延误事件，严重侵害消费者的权益，甚至导致消费者付出更大的时间成本和机会成本。为了缓解乘机客户和航空公司之间的矛盾，保险公司积极介入并推出航班延误险，此举可以说在一定的程度上会促进民航行业的发展，同时，民航行业的发展也会给保险公司带来收益，因此，这对两者而言是双赢的。既然未来民航行业将有较大的发展空间，所以，也就意味着航班延误保险具有广阔的发展前景，保险公司服务的水平越高，航空公司的发展水平也应该越高，这样，双方的合作将会带来更多的利益，极大的促进我国民航行业的发展。六、模型评价与推广航班延误是当前国际民航业发展中的一大难题，也是顾客对航空服务质量不满意的主要内容，如果航空公司不能有效的提高内部管控能力，增强其对飞机起飞的运筹把握程度，很有可能给客户造成诸多不便，甚至会激发客户与航空公司之间的矛盾，对航空公司的声誉造成影响。本文针对目前国内民航行业的基本特征，构建模型对航空公司的预订票策略进行充分研究，并对航空公司的运营提供建议，然而，消费者作为航空公司的服务对象，有时因为航班的意外延误会对航空公司的声誉造成一定的影响，所以，本文从市场的角度出发，认定保险公司为消费者提供保险产品，并对保险公司的经营风险构建概率模型进行描述，因此，本文所提出的方法和模型对现实具有一定的意义，并且，可以在相对合理的条件推广到其它航空公司或者其它保险公司。参考文献1数据来源于美国航空数据网站Flight States公布的2014年5月全球航空公司的准点表现报告 1数学模型（第四版）M姜启源，谢金星，叶俊，2011，北京，高等教育出版社2概率论与数理统计(第四版)M 盛骤，2008.10，北京，高等教育出版社3赵秀丽，朱金福，郭梅，不正常航班延误调度模型及算法D 南京 2008.044王红，刘金兰，曹卫东，航空公司航班延误预警管理模型与分析D 西安 2009.045徐涛，荣耀，王建东，基于SOA的民航航班延误波及分析与预警系统D 天津 2009.076李俊生，丁建立，基于贝叶斯网络的航班延误传播分析D 天津 2008.117刘光才，刘雷，美国减少航班延误的有效途径及启示D 天津 2010.048董念清，中国航班延误的现状、原因及治理路径D 北京 2013.119丁建立，陈坦坦，刘玉洁，有色-时间Petri网航班延误模型及波及分析D 天津 2008.12附录附录一 2006-2011年我国航班变化情况MATLAB代码x=2006:2011;y=1530443 1613786 1528208 1759438 2010652 2204147;y1=276185 281831 254140 322601 403511 343050;plot(x,y,*r-,x,y1,ob-);title(06-11年我国航班变化情况);xlabel(年分);ylabel(航班总数与延误的航班数);legend(航班数,延误航班数);legend(boxoff);axis square;附录二 各航班延误原因的