一、一、 问题重述问题重述国产大飞机 C919 项目的上马和良好进展向世界展示了中国的大国实力近日，中国商用飞机公司销售经理陈进表示国产大飞机 C919 在燃油消耗方面将比目前所有机型减少 12-15%，且在上市时定价将会低于 5000 万美元/架，大大低于波音和空客同类产品的定价这将为国产大飞机开辟广阔的市场在放手开发研制一种新飞机时，除了技术细节外，还有很多经济问题需要回答其中最重要的是飞机制造商需要知道他的原始投资是否能够收回来，多久才能收回来这就要预测一下飞机的上市价格和市场前景毕竟对于飞机这样庞然大物的高科技商品，高投入，也高风险，周期还长，其销售价格除了和制造商的制造成本、航空公司的运营成本有关外， 还和同类型飞机的市场前景、竞争格局等其它因素有关而飞机的技术数据往往能够一定程度上反映其中的一些因素我们要解决的问题有：1、 假如你是中国商用飞机公司销售经理，请你对比同类客机 B737-800，以航空公司的角度，结合飞机运营成本，评估一下国产大飞机 C919 未来的市场潜力；2、 假如你是飞机制造商，请你建立一合理的数学模型为将新开发的商用飞机预测一个较合理的价格，并结合你的模型对 C919 的价格进行预测；3、 以你所建立的模型，结合相关数据，计算当前A380-800 和 B737-800 两种商用飞机的价格，并说明所建价格预测模型的优缺点；二、二、 C919C919 市场潜力分析市场潜力分析2.12.1 基本假设基本假设在问题（1）中，根据实际情况，可以做如下假设来简化问题：1、本题所进行的 C919 和 B737-800 的运营成本的对比，是在同一航空公司、同第 1 页 共 28 页 样的航线、同一时间段下所进行的，也就是说，忽略各航空公司经营模式和飞行天气的影响。
而且，两者每年起降次数、飞行时间、里程相同2、基于第一个假设，两种机型的运营成本中的机务人员工资、航空食品配送费、飞机入库费用、起飞场地费等都可以认为相同3、由于 B737-800 的座位数为经济舱 158 座，C919 有三种类型，分别为156 座、158 座、180 座，可以认为两者座位数大致相等，因此在计算平均每座运营成本时对二者取相同座数4、不考虑汇率变动通货膨胀5、C919 的价格按官方透露的比同款机型低 10%来计算2 2．．2 2 问题分析问题分析要分析一种待上市产品的市场潜力， 必须剖析它与同类产品相比所具有的竞争力，并以现有产品的市场为基础来预估考虑到航空运输业的特点，以下我们将通过层次分析法按重要性大小从几方面对B737-800和C919进行对比1、运营成本：航空公司生产的最终目的是在满足经济和社会需求的同时，最大限度地获取利润 在交通技术日新月异的今天， 各航空公司不断进行价格战，机票已不再昂贵， 各航空公司尤其是民营航空公司将不断降低成本作为实现航空公司盈利目标最主要、最可靠的途径购买飞机时， 航空公司也就更看重飞机在使用维护阶段的经济性，其评估标准是飞机的直接运营成本(DOC)和间接运营成本(LOC)。
LOC整体呈现很不稳定状态，并且变化具有非线性的趋势，不好控制，所以对于飞机使用维护经济性分析的依据选定为直接运营成本(D0C) DOC费用包括拥有飞机费用(折旧费、利息、保险金、高价件摊销、经营性租赁等 )、空勤组费用(驾驶员成本、客舱服务员成本)、维修成本(机体维修成本、发动机维修成本)、运行费用、燃油成本和餐食费等考察1995～2000五年直属航空公司的成本结构可以发现，航油、维修、折旧、起降、旅客服务等费用占据全部DOC成本的80%以上又由于假设二的存在，应认为C919和B737-800的旅客服务费、起降费用相同，因此我们选择航行飞机每航行一公里的折旧费、燃油费、维修费和航材费的总和来作为飞机的运营成本第 2 页 共 28 页 2、安全系数：近年来由于飞机动力、操纵、起降系统故障和极端天气引发的飞行事故频频出现，给航空公司和社会都造成巨大损失，安全问题已成为备受关注的问题因此，商用飞机的安全可靠性也成为航空公司选购时的一项重要指标由于C919尚未投入使用，没有可统计的事故率，这里我们采用两者发动机和操纵系统的技术可靠系数来衡量飞机的安全可靠性3、舒适度系数：由于各机型的设计方案不同，座椅排距、行李箱高度、机身稳定性都影响着乘客舒适度，而乘客满意度是决定该航班满客率的重要因素，也在很大程度上影响着航空公司的盈利水平。
所以，飞机的舒适度也是其竞争力的一项重要指标我们用客舱宽度、座椅排距来衡量舒适度4、环保系数：随着全球变暖的趋势不断加剧，环保因素越来越多地被考虑到工业生产和服务行业中来，环保力度也逐渐成为影响企业声誉的重要因素这里我们用航行每公里的二氧化碳排放量来衡量飞机的污染水平基于以上考虑， 我们将采用综合指数评价方法来比较B737-800和C919对购买方的吸引力，并以B737-800为参考量预估C919的市场潜力对于其中最关键的运营成本，我们在已经成熟的兰德模型下用工程估算法来计算2.32.3 符号约定符号约定；t：飞机使用时间（单位：年）：飞机的年折旧率；k1：飞机折旧 N 年后的成本核算；k2：飞机 N 年所消耗的航空燃油总费用；k3：飞机 N 年所维修总费用；a0：飞机购入时的价格（单位：百万） ；b0：：飞机年耗油费用（单位：百万） ；f(∑C)：飞机的年维修费用；i：飞机年维修费用中各组成部分的权重；第 3 页 共 28 页 C1：航材费用；C2：定检维修费用；C3：机务人员工资；E综：飞机潜力综合指数；Ei：第 i 个评价项目的值；i：每个评价项目在该机型潜力综合指数中的权重2.42.4 模型建立与计算模型建立与计算1、在计算飞机运营成本时，我们采用工程估算法。
这种方法的思路是将飞机的所需的维修费用项目细分， 直到最小的本费用单元， 然后自下而上累加而成单架飞机运营成本y,t∑∑ki,ta0a0(1)t；其中飞机第t年后由于折旧而付出的成本k1由于发动机在老化过程中燃料利用效率会不断下架，因此k2飞机第t年后共用维修费用k3b0∑∑(1)t；f c×t；；在计算f c时，我们采用已经很成熟并被广泛使用的兰德 DAPCADAPCA IVIV 模型该模型是飞机发展与采购费用（Development and Procurement Costs ofAircraft――DAPCA）模型的最终形式由此得到f c011C2C23C30.00348C10.000475C20.000671C3181.4377这个模型比较简单，所以我们选择使用 matlab 分别作出 C919 和 B737-800的基于时间序列的运营成本曲线，如下：第 4 页 共 28 页 国 产 c919和 737-800运 营 成 本 比 较350国 产 C919运 营 成 本 曲 线737-800运 营 成 本 曲 线300250运营成本/million200150100500024681012运 营 时 间 t/年14161820图 1两者运营成本比较从图中我们可以明显看出，两者的运营成本都随使用时间非线性上升，且斜率越来越高。
但 C919 的运营成本较 B737-800 低，且随着使用时间的增长二者差距也随之增大，优势越来越明显另外，鉴于全球航空汽油的价格呈上涨趋势，我们在两个增长后的航油价格下得到以下两个结果：第 5 页 共 28 页 国 产 c919和 737-800运 营 成 本 比 较400350300运营成本/million国 产 C919运 营 成 本 曲 线737-800运 营 成 本 曲 线250200150100500024681012运 营 时 间 t/年14161820图 2航油价格增长 5%时运营成本比较国 产 c919和 737-800运 营 成 本 比 较400350300国 产 C919运 营 成 本 曲 线737-800运 营 成 本 曲 线运营成本/million250200150100500024681012运 营 时 间 t/年14161820图 3航油价格增长 10%运营成本比较由图可知，当航油价格上涨时，二者运营成本都会有所升高，但 C919 的运营成本增长缓慢，较之 B737-800 的优势更明显根据官方信息，将两种机型的使用寿命都取为 25 年，以目前的航油价格为第 6 页 共 28 页 标准，则在使用周期内两者运营成本之比为436.1313.51.391。
设 B737-800 的运营成本效率为 1，则 C919 运营成本相应为 1.3912、在比较安全系数时，鉴于两者制造工艺、材料、发动机、电子设备水平大体相当，其安全系数皆取为 13、在比较舒适度时，乘客的舒适感与空间拥挤度有直接的关系，B737-800的两排座椅间距为 30 厘米，C919 的间距为 32 厘米，设 B737-800 舒适度为 1，则 C919 舒适度为 1×32=1.067304、在比较环保系数时，由于二氧化碳排放量与单位航程燃油量成正比，将B737-800 的环保系数设为 1，则根据官方消息，C919 排放量将是类似机型的一半，即其环保系数为 2现在来求 B737-800 和 C919 的潜力综合评价指数：E综=∑iEi其中按照各评价项目的重要程度取将两者参数带入计算，得到EC91910.5，20.3，80030.1，40.11.302,EB7371销量与购买方对飞机的满意度成正比，所以可以认为C919 的销售量为B737-800 的 1.302 倍目前中国大陆共有商用飞机1383 架，其中波音公司占有53%的市场份额,2009 年 B737-800 的销售量为 128 架。
而未来五年中国（包括香港、澳门）的商用飞机市场将有 3700 架的需求量.按照我们计算的结果， 根据比例可粗略预计 2016 年 C919 将拿到 446 架的订单2.52.5 模型评价与改进模型评价与改进我们认为综合指数更有利于客观全面反映水质的标准， 它体现出了一种细分的评测标准当然，由于掌握资料广度所限，我们略去了一些因素，采取的假设也比较理想化，选取的加权系数是否有一个通用的标准还有待商榷基于相似机型 B737-800 的销售量的市场潜力预测简单易行，但是同时也略显粗糙，忽视了未来几年市场可能的变化和其他竞争者的可能竞争在工程估算法的应用中， 将总系统费用分解为许多用费用方程联系起来的项第 7 页 共 28 页 1结果准目细节，其反映的因果关系与实际情况更加接近，因此具有如下优点：○2能对竞争的各个方案研究其费用差异，○3允许进行详细地模拟和灵敏度确，○分析 但是它也具有对数据要求高、估算结果很难进行评价与鉴定等明显的缺点，有待改进对于本问中的模型，我们可以进一步结合民机发展的几十年中市场的变化，建立基于时间序列的市场潜力预测模型对其进行改进三、三、C919C919 价格预测价格预测3.13.1 基本假设基本假设1、由于商飞是国防性质的国有，为简化计算，不考虑税率问题。
2、C919的研制不仅仅是一种商业行为，还关系到中国的经济安全和国力体现，牵扯到众多的不易量化的政治因素为了简化和明朗问题，这里我们只考虑经济因素，不考虑政治利益3、使用兰德模型计算成本时，只考虑中美人工费用差异，不考虑技术差异4、C919的研发情况由已经生产的国产运十、新舟60来类推3.23.2 问题分析问题分析对于民用飞机而言，其出售时要收回研究、发展、试验与鉴定费用和生产费用两部分 （包括合理的利润） 为了收回这两部分费用， 需要合理确定民机售价国外各大飞机制造商，如波音公司，已经根据经验数据建立了非常合理的定价系统，能够对新机型很快作出合理的成本和价格预测，如B737-800的定价然而，由于缺乏自主研发新机型的经验和有效数据，国内很少有人进行这方面的研究，对飞机的采购价格预测也仅局限于军用飞机方面 因为国产商用飞机和军机同样存在样本小、数据少、 性能因素影响显著等特点，在解决本题时我们可以借鉴军机价格预测中广泛使用的偏最小二乘法、类比法和参数法，并使用具有学习训练功能的神经网络来对其进行修正价格，是价值的货币表现，价值是其决定、支配因素商用飞机是一种高成本、高科技、高价格的产品，对于制造商来说，出售飞机的目的是收回成本即在第 8 页 共 28 页 研制期间花费的大量研究、试验经费以及鉴定费用，而经费的多少与研制难度即飞机的各项性能指标有关；对于购买方（主要是航空公司）来说，飞机的运营成本及其能够带来的利益是它们关心的主要问题， 而这两方面也与飞机的燃油利用效率、载重、巡航速度等技术参数息息相关。
因此，从买卖双方来说，一架飞机的价值都应该取决于它的功能，也就是说，飞机是一种基于功能的产品由此，我们可以认为飞机的功能是其价格的主要影响因素 因为飞机的参数可以在很大程度上反映其功能并且比较容易量化， 我们选取一些性能参数作为变量来对价格进行多元拟合又由于对我们来说参数与价格的关系还是一个黑箱系统，所以选用逐步回归来选择在价格函数中真正起作用的参数步骤如下：1、采集产品数据搜集目前仍在使用的所有商用飞机的各项性能参数，构建性能参数及价格数据的初始样本集2、 性能参数数据的标准化对所采集的样本数据进行整理,对于可以直接量化且数据大小与产品价格明显相关的功能数据,可以不予标准化直接作为输入进行计算,对于不可量化或数据大小与产品价格不明显相关的数据,则根据产品支持的同类功能的复杂程度采用0～1 进行打分3、对性能参数和价格进行逐步回归根据偏相关系数筛选出对价格影响较大的参数，最后确定价格的表达式4、 C919价格的预测通过分析数据，我们发现不同型号的商用飞机之间功能、造价差异较大，四发的大型客机样本的引入可能造成预测两发的支线飞机C919的价格时的误差， 因此通过估算法，参考和C919参数相近的A320、B737-800的功能、价格关系，比较三者之间的差异，对C919价格进行预测，并与第一种方法比较、调整。
另外，对于民用飞机而言，其出售时要收回研究、发展、试验与鉴定费用和生产费用两部分（包括合理的利润）为了收回这两部分费用，需要合理确定民机售价因为研究、发展、试验与鉴定费用是确定的，所以必须根据市场需求，对生产飞机的架数做出合理的预测，以确定每架飞机的售价中要包含多少成本先使用第一个问题中求出的C919的预估销量， 计算该销量下的成本， 与售价比较；再求在保证收回成本的情况下的销售量， 结合市场发展情况进行前景分析和对策第 9 页 共 28 页 决策3.33.3 模型建立与计算模型建立与计算1 1、基于功能驱动的产品价格预测模型、基于功能驱动的产品价格预测模型根据问题分析，首先进行样本数据采集商用飞机的性能参数有：○ 1 几何参数：翼展、机长、高度；2性能参数：经济/最大巡航速度、最大起飞重量、空重、业载、发动机台○数/型号、油箱容量、载客量由于参数过多而样本太小， 需要把内在相关的一些参数进行合理的整合和选择1翼展的大小主要由飞机的气动布局决定，不能较好地反映飞机的规格，所○以我们选择机身的剖面积来衡量飞机规格；2商用飞机一般在经济巡航速度上使用，所以舍弃最大巡航速度，选择经济○巡航速度作为一个参数；3航空公司最注重的是经济效益，○所以飞机的最大载客量和最大业载也是影响价格的重要因素；4发动机是飞机的心脏，○同时发动机成本和维修费用也是飞机成本的重要组成部分。
根据空气动力学，最大推力与最大起飞重量有关，为了减少参数交互，舍弃最大推力和发动机台数综上所述，使用Matlab，我们先选取飞机规格Xvo、经济巡航速度Xve、最大载客量Xp、、飞机运营成本Xc、最大业载Xw作为参数对价格进行逐步回归部分波音、空客商用飞机数据见附录当然，这些数据之间还是存在一定的交互性，而且在这个黑箱模型下我们也不能确定它们是否对价格函数作用下面开始进行参数的筛选：第一步，将Xvo、Xve、Xp、Xc、Xw作为变量对价格P进行函数拟合P∑miXi（ivo,ve, p,c,w）第 10 页 共 28 页 并分别求每个回归变量系数m的置信区间得到六个回归系数的估计值及其置信区间如下：m1的置信区间为（-0.0025，0.0153）；m2的置信区间为（-0.0001，0.0512）；m3的置信区间为（0.0676，0.2004）；m4的置信区间为（0.2114，0.2634）；m5的置信区间为（0.0867，0.1974）检查它们的置信区间发现,Xvo和Xve的系数的置信区间包含了零点, 表明回归变量Xvo和Xve对P的影响不是太显著,因此从模型中移出次变量Xvo和Xve第二步，由于Xp、Xc、Xw对P具有较明显的显著性,用更为精确的拟合曲线来拟合表示Xp、Xc、Xw和P的关系。
做Xp、Xc、Xw分别对P的散点图，再分别进行拟合，结果如下：300280260240220200180160140120100data 1 linear150200250300350400图4以Xp为变量的P的拟合曲线第 11 页 共 28 页 300280260240220200180160140120100100data 1 quadratic120140160180200图5以Xc为变量的P的拟合曲线300data 1 cubic2502001501005010203040506070图6以Xw为变量的P的拟合曲线具体表示式如下：P1P2P30.50059XP30.0025453XW93.0710.372XC20.45209XW20.029436XC636.8524.218XW123.16对应的Norm of residuals分别为46.2、43.14、45.645，由此可以看出,对第 12 页 共 28 页 Xw,应的拟合方程的拟合程度较高 所以,用拟合后的曲线来代替原来的Xp、Xc、得到如下回归模型:P∑1iPi (1)0.0705，2其中0.0828，30.1260由于计算Xc时涉及到一些性能参数，所以Xc、Xw之间存在相互制约，它们的交互作用会对P有影响， 不妨简单地用Xc、Xw的乘积代表它们的交互作用，于是，模型（1）增加一项，得到：P1 1P22P33P4XCXW（2）再次进行回归分析，得到修正后的系数10.0676,20.0664,30.1044,40.0002相关系数为0.9773，说明该模型比较接近真实情况。
至此，第一个模型已经建成将C919的有关数据代入模型（2），得到第一个预测价格5984.65万美元2 2、类比法预测、类比法预测 C919C919 价格价格由第一个模型得到的C919价格为5984.65万美元，这与官方公布的低于5000万美元有不小的差距经过分析数据，我们发现样本中的四发大型商用飞机与两发支线飞机相比，各项性能参数以及价格差距甚大，从而影响了模型（2）的通用性我们又将一些样本外的机型的数据代入进行检验，结果显示该模型用于大型飞机的价格预测时比较准确， 而在对中小型飞机的预测中结果总是偏高 因此，1四发大型飞机样本和○2两发中小型飞机样我们决定使用类比法， 将样本划分为○2选取与C919性能参数比较接近的B737本两部分，分别进行逐步回归其中样本○系列和A320系列的参数具体步骤同上节，得到两个模型，形式与（2）相同，系数不同：模型（3）P'11P'2P2'3P3'41）XCXW（适用于样本○第 13 页 共 28 页 模型（4）P''11P''2P2''3P3''42）XCXW（适用于样本○运用模型（4）再一次计算C919的价格，结果为4710.53万美元，这一结果比较符合官方消息。
3 3、成本导向定价模型、成本导向定价模型飞机，尤其是商用飞机的生产和销售时是以盈利为目的的高投资的商业行为，必须结合市场来制定合适的价格保证盈利虽然C919项目有国家资金支持，并且考虑到一定的政治利益因素， 为了国家的利益， 我们还是要保证成本的回收，甚至一定的利润因此，在问题一预测的销量下求解C919的利润以及在预测盈亏平衡点的价格并与上一模型求解结果比较具有很大的实际意义1、首先我们来计算C919的研制成本兰德 DAPCA IV 模型是飞机发展与采购费用（Development and ProcurementCosts of Aircraft――DAPCA）模型的最终形式，在飞机寿命周期费用分析领域有相当的影响力按照 1986 年定值美元，兰德DAPCA IV 模型中工时、费用的计算公式如下：工程工时工艺装备工时制造工时质量控制工时发展支援工时飞行试验费用制造材料费用发动机制造费用HEv5.18We0.777V0.894Q0.1637.21We0.777V0.696Q0.263HTHM10.46We0.82V0.484Q0.641HQ0.767HMCDCFCM33.54We0.630V1.3968.42We0.325V0.822FTA1.2115.54We0.921V0.621Q0.799CEng1.548 0.095Tmax243.25Mamax0.54tn2228研究、发展、试验与鉴定费用＋生产费用C成本HEREHTRTHMRMHQRQCDCFCMCEngNEngCav第 14 页 共 28 页 式中W：空重（kg）；V：最大飞行速度（km/h）Q：产量；FTA：飞行试验机架数（一般为2~6 架）；NEng：总产量乘以每架飞机的发动机台数；Tmax：发动机最大推力（kg）Mmax：发动机最大Ma 数；tti：涡轮进口温度（K）RE、RT、RM、RQ：综合费率（数值见附录）Cav：航空电子设备费用C919有关参数如见附录：代入计算得到C成本g Q鉴于C919应用了 20%的复合材料，在上述结果上应乘以相应的软糖系数1 0.81.5 0.21.1；考虑到中国的人工费用与欧美的差异以及目前美元价值与1986年的差异， 在上述结果上再乘以一个换算系数1.522， 最终得到成本C成本关于产量Q的函数C成本1.647g Q（5）2、当使用平均成本加成定价法对C919价格进行预估时，根据千分之一法则将利润率定为0.1%。
在这种情况下得到P加成法C成本÷Q×1+0.001=6078.65（万美元）这表明若要获得千分之一的利润，定价需在六千万美元左右3、鉴于国产大飞机C919的研制与B737-800相比基础较差，国家前期投入较高，因此收回成本是一项重要的任务用盈亏平衡定价法来计算保本点的价第 15 页 共 28 页 格及产量：1取模型○（4） 中计算的较为合理的单价4710.53万美元， 代入模型 （5） ，得到保本点销量Q保本673（架）6072.58（万2取问题一中的销售量446架，代入模型（5），得到P保本○美元）从这两个结果我们可以看出，想要保证回收成本，必须上调价格或者提高产量但是，C919的最大优势在于低廉的价格，一旦价格上调，必然影响销量，这样很难实现盈亏平衡因此，只能提高它的销量由于问题一的销量是在中国2016年市场饱和情况下求得的， 所以唯一的手段是大力开拓海外市场， 而我们认为其重点应该放在发展中国家等价格和燃油对决策影响巨大的民机市场上四、四、A380-800A380-800 和和 B737-800B737-800 价格计算与模型的评价改进价格计算与模型的评价改进4.1A380-8004.1A380-800和和B737-800B737-800价格的计算价格的计算这两种机型的有关数据见附录。
其中，A380-800属于大型商用飞机，B737-800属于中小型商用飞机，所以应该分别使用模型（3）和（4）进行计算，结果如下：CA380CB737800333.33 (million)(真实价格为327.4million， 误差为1.81%)67.77(million) （真实价格为72.5million，误差为6.52%）800由此看出，我们的模型能够比较准确地计算各种机型的价格4.24.2模型评价与改进模型评价与改进1 1、模型评价、模型评价多元逐步回归针对本题中众多的原始参数和较为复杂的互相制约， 有效地筛选出对价格函数起作用的变量，能够较准确地反映各有效变量对价格函数的作用，稳定性也比较强另外，此种算法充分利用了 Matlab工具箱的优势，简化了运算然而，由于各参数之间存在交互作用，参数之间的关系以及各参数对价格第 16 页 共 28 页 函数的真实作用并不明朗，这个模型实际上是在黑箱下得到的结果，对不同参数变化的灵敏度有很大差异另外，大型飞机与中小型飞机两类分歧显著， 对结果的精确度造成了极大的影响类比法正是针对这个问题而补充的，它将待求对象分为两类，分别通过逐步回归得到不同的模型，提高的预测的准确度，更为科学与合理。
成本导向定价法的引入进一步解决了题干中的成本回收和盈利问题， 为C919的生产、上市提供了合理、可行的建议，提高了本文结论的实用价值2.2.使用最小二乘支持向量机对模型进行优化使用最小二乘支持向量机对模型进行优化逐步回归分析通常需要大样本容量空间的支持，并且模型泛化能力差，而本题中能够使用的样本经过分类后数量过少，不足以保证模型的精准和稳定为了解决本模型在小样本下的失准问题， 我们又选择了能够在Matlab工具箱中方便使用的最小二乘支持向量机作为预测工具，对模型进行优化这种算法基于以小样本数据为主要研究对象的统计学习理论,采用结构风险最小化准则来代替传统的经验风险最小化准则,根据有限的样本信息在模型的复杂度和模型的学习能力之间寻求最佳折衷,同时最小化VC维和经验风险,以获取期望风险的最小化,它不像神经网络等其他预测方法一样需要大样本空间的支持,并最终转化为二次寻优问题,可获得问题的全局最优解,较好地解决了易陷入局部最优的问题 最小二乘支持向量机优化指标采用了平方项 ,只有等式约束,而没有标准支持向量机的不等式约束,从而推出不同的一系列等式约束 ,将二次规划问题转变成线性方程组的求解,简化了计算复杂性,因此具有较快的运算速度。
模型建立如下：p X∑diK xi,xjbi式中,X为飞机的参数数值,K xi,xj为核函数,式计算：式中为待定参数,其值可由下pp1, p2,?··，piT；；e11,1,? ··,1；；Qij第 17 页 共 28 页TxixjK xi,xj；；xi、、xj为 飞机参数和飞机价格具体代码和核函数求解见附录经过多次学习训练后得到 C919 的预测价格为 4967.56 万美元，相比前述模型更为可信五、参考资料五、参考资料【1】 《数学模型》 （第三版）姜启源 谢金星 叶俊高等教育出版社 2003年【2】 《基于 MATLAB 和 LINGO 的数学实验》肖华勇西北工业大学出版社2008 年【3】 《人工神经网络原理及仿真实例（第二版） 高隽机械工业出版社2007 年【4】 《飞行费用与效能分析》【5】 《基于 MATLAB 的多元非线性回归模型》董大校 2009 年【6】 《基于功能驱动的产品价格预测》黄训江六、附录六、附录6.16.1 飞机运营成本飞机运营成本 matlabmatlab 代码代码 %商用飞机的运营成本模型%%前提1：作为运营成本的测评，我们假设参加测试的N种飞机在单位时间内航程是相同的%%变量说明 Xi:第i种飞机的运营成本%%alphai:第i种飞机的平均年折旧率%%ai:第i种飞机的价格区间%clcclear;alpha=[0.05,0.03];p=1+alpha(1);q=1-alpha(1);t=1:0.5:20;第 18 页 共 28 页 a0=72.5+8.5*rand(1);time=2500+1000*rand(1);b0=time*2650*8400/1000/6.8;a=a0-a0*q.^t;b1=b0*(1-p.^t)/(1-p);b=b1/1000000;c0=3;c=c0*t;x=a+b+c;plot(t,x,'r');h=1+alpha(2);i=1-alpha(2);t=1:0.5:20;u0=(72.5+8.5*rand(1))*0.9;time1=2500+1000*rand(1);v0=(time*2650*8400/1000/6.8)*0.87;u=u0-u0*i.^t;v1=v0*(1-h.^t)/(1-h);v=v1/1000000;w0=3;w=w0*t;y=w+u+v;plot(t,y,'-r.',t,x,'-b');xlabel('运营时间t/年');ylabel('运营成本/million');title('国产c919和737-800运营成本比较');legend('国产C919运营成本曲线','737-800运营成本曲线');grid on;%年维修费用的求法%%alpha=0.05;%x1=[3.5 5.3 5.1 5.8 4.2 6.0 6.8 5.5 3.1 7.2 4.5 4.9 8.0 6.5 6.6 3.7 6.27.0 4.0 4.5 5.9 5.6 4.8 3.9];第 19 页 共 28 页 %x2=[9 20 18 33 31 13 25 30 5 47 25 11 23 35 39 21 7 40 35 23 33 27 3415];%x3=[6.1 6.4 7.4 6.7 7.5 5.9 6.0 4.0 5.8 8.3 5.0 6.4 7.6 7.0 5.0 4.4 5.57.0 6.0 3.5 4.9 4.3 8.0 5.8];%y1=[33.2 40.3 38.7 46.8 41.4 37.5 39.0 40.7 30.1 52.9 38.2 31.8 43.3 44.142.5 33.6 34.2 48.0 38.0 35.9 40.4 36.8 45.2 35.1]';%数据元素的初始化%%a=ones(24,1);%c=[x1;x2;x3]';%X=[a,c];%[b,bint,r,rint,stats]=regress(y1,X,alpha);%fprintf('线性回归系数：\n');%fprintf('%6.4f\n',b);%fprintf('回归方程的统计量：\n');%5fprintf('%6.4f\n',stats);6.26.2 问题二中逐步回归问题二中逐步回归 matlabmatlab 代码代码clcclear;%以下是数据的初始化%%飞机长度%length=[70.570.570.756.370.763.744.448.554.947.330.533.4];%飞机高度%high=[19.319.319.319.919.318.412.915.915.913.511.311.1];for i=1:12 V(i)=length(i)*high(i);%飞机的横截面%end%经济巡航速度%velocity=[935935935935935935885850850850856856 ];%最大载客量%passenger=[291 291 375400300 200330110134375214150 ];第 20 页 共 28 页 %运营成本%cost=[149 154 168.6 172 140 123 156 106 110 197 133 118.4 ];%受发动机台数制约的最大业重%engine=[3065.239.5406315582017];%发动机台数的初始化，此处可根据需要修改为2%fadongji=2;%12个样本飞机的售价(单位：million dollars)%%Y矩阵包含着所有4发商用飞机的价格%y=ones(1,12);%用来储存每个样本的售价&%747-400com%y(1) = 234.0 +(266.5-234.0)*rand(1);%747-400%y(2) = 238+(268.0-238)*rand(1);%B747CMB%y(3)=273.0+(308.0-293.0)*rand(1);%B747SP%y(4)=281.5+(304.5-301.5)*rand(1);%B747F%y(5)=223.5+(276-223.5)*rand(1);%B777%y(6)=205.5+(231.0-205.5)*rand(1);%B707%y(7)=237.5+(263.5-237.5)*rand(1);%B767-200%y(8)=127.5+(139.0-127.5)*rand(1);%B767-300%y(9)=144.5 + (161.5-144.5)*rand(1);%B757-200%y(10)=252.5+(260.5-252.5)*rand(1);%B737-200%y(11)=194.0+(205.5-194.0)*rand(1);第 21 页 共 28 页65.52514 %B737-300%y(12)=150+(155.5-150)*rand(1);%Z矩阵代表的样本皆为2发样本， 由于样本值浮动较大所以此处我们皆取平均值%z=ones(1,12);%737-600%z(1)=73.1436;%737-700%z(2)=74.5521;%737-800%z(3)=78.4714;%737-900ER%z(4)=80.2720;%A316%z(5) = 72.0097;%A317%z(6) = 57.3962;%A318%z(7) = 68.4223;%A319%z(8) = 36.1367;%A320%z(9) = 43.2333;%A321%z(10)= 73.0657;%A322%z(11)= 53.7438;%A323%z(12) = 41.6015;%绘制5组数据的散点图%subplot(1,6,4),plot(V,y,'*');subplot(1,6,5),plot(velocity,y,'+');subplot(1,6,1),plot(passenger,y,'o');subplot(1,6,2),plot(cost,y,'+');subplot(1,6,3),plot(engine,y,'*');%1.首先对5个因素进行线性回归%第 22 页 共 28 页 alpha=0.05;%设置置信区间%X=[V;velocity;passenger;cost;engine]';%根据不同的发动机规格选择不同的样本值进行模拟%if(fadongji==4) y1=(y)';else y1=(z)';end%利用线性回归求解的回归系数%[b,bint,r,rint,stats]=regress(y1,X,alpha);fprintf('线性回归系数：\n');fprintf('%6.4f\n',b);fprintf('回归方程的统计量：\n');fprintf('%6.4f\n',stats);subplot(1,6,6),rcoplot(r,rint);%绘出时序残差图%%2.进行非线性回归%for j=1:12 finalpassenger(j)=passenger(j)*0.50059+93.07; finalcost(j)=-0.029463*(cost(j))^2+10.372*cost(j)-636.85;finalengine(j)=0.0025453*(engine(j))^3-0.45209*(engine(j))^2+24.218*engine(j)-123.16;end%非线性化后回归系数的求解过程%beta=0.05;e=[finalpassenger;finalcost;finalengine]';[b1,bint1,r1,rint1,stats1]=regress(y1,e,beta);fprintf('线性回归系数：\n');fprintf('%6.4f\n',b1);fprintf('回归方程的统计量：\n');fprintf('%6.4f\n',stats1);%进行系数交互%for k=1:12 new1(k)=finalengine(k)*finalcost(k);end%最终回归系数的求解%第 23 页 共 28 页 final=[finalpassenger;finalcost;finalengine;new1]';[b2,bint2,r2,rint2,stats2]=regress(y1,final,beta);fprintf('线性回归系数：\n');fprintf('%6.4f\n',b2);fprintf('回归方程的统计量：\n');fprintf('%6.4f\n',stats2);%模型检验%%根据模型预测C919的价格%passenger(12)=190;cost(12)=128.5;engine(12)=15;finalpassenger(12)=passenger(12)*0.50059+93.07;finalcost(12)=-0.029463*(cost(12))^2+10.372*cost(12)-636.85;finalengine(12)=0.0025453*(engine(12))^3-0.45209*(engine(12))^2+24.218*engine(12)-123.16;new1(12)=finalengine(12)*finalcost(12);C919=b2(1)*finalpassenger(12)+b2(2)*finalcost(12)+b2(3)*finalengine(12)+b2(4)*new1(12);fprintf('根据模型，我们预测C919的价格为=%6.4f\n',C919);%根据模型预测737-800的价格%passenger(12)=186;cost(12)=156.9;engine(12)=26;finalpassenger(12)=passenger(12)*0.50059+93.07;finalcost(12)=-0.029463*(cost(12))^2+10.372*cost(12)-636.85;finalengine(12)=0.0025453*(engine(12))^3-0.45209*(engine(12))^2+24.218*engine(12)-123.16;new1(12)=finalengine(12)*finalcost(12);P737=b2(1)*finalpassenger(12)+b2(2)*finalcost(12)+b2(3)*finalengine(12)+b2(4)*new1(12);fprintf('根据模型，我们预测737-800的价格为=%6.4f\n',P737);6.36.3 问题三中最小二乘向量机求核矩阵代码问题三中最小二乘向量机求核矩阵代码x =[70.5 70.5 70.7 56.3 70.763.744.448.554.947.330.533.431 46.746.454.159.4;第 24 页 共 28 页 19.3 19.3 19.3 19.919.318.412.915.915.913.511.311.111.1 15.8 15.816.516.7; 935 935 935 935935935885850850850856856850 850828 850850; 360 400 291 291375 375148214225200128132130 204228274375; 125 154 135 172 125 123 156 178 148197 179 165 176 154 146 143 153; 44 4 4 4 2 4 2 2 2 2 2 22]y=ones(1,17);%用来储存每个样本的售价&%747-400com%y(1) = 234.0 +(266.5-234.0)*rand(1);%747-400%y(2) = 238+(268.0-238)*rand(1);%B747CMB%y(3)=293.0+(308.0-293.0)*rand(1);%B747SP%y(4)=301.5+(304.5-301.5)*rand(1);%B747F%y(5)=223.5+(276-223.5)*rand(1);%B777%y(6)=205.5+(231.0-205.5)*rand(1);%B707%y(7)=237.5+(263.5-237.5)*rand(1);%B767-200%y(8)=127.5+(139.0-127.5)*rand(1);%B767-300%y(9)=144.5 + (161.5-144.5)*rand(1);%B757-200%y(10)=252.5+(260.5-252.5)*rand(1);%B737-200%y(11)=194.0+(205.5-194.0)*rand(1);%B737-300%y(12)=150+(155.5-150)*rand(1);%B737-500%y(13)=161.0+(171.5-161.0)*rand(1);第 25 页 共 28 页2 2 2 y(14)=167;y(15)=176;y(16)=187;y(17)=178;xmin1 =min(x(1,:));xmax1 =max(x(1,:));xmin2 =min(x(2,:));xmax2 =max(x(3,:));xmin3 =min(x(3,:));xmax3 =max(x(3,:));xmin4 =min(x(4,:));xmax4 =max(x(4,:));xmin5 =min(x(5,:));xmax5 =max(x(5,:));xmax6 =max(x(6,:));xmin6 =min(x(6,:));net.trainparam.epochs =2000;net =newff([xmin1,xmax1;xmin2,xmax2;xmin3,xmax3;xmin4,xmax4;xmin5,xmax5;xmin6,xmax6],[6,1],{'logsig','logsig'});net =train(net,x,y);x1 =[70.5 70.570.7; 19.3 19.3 19.3; 935 935 935; 360 400 291; 125 154 135 ; 4 44]y1 =sim(net,x1);6.46.4 相关机型参数数据表相关机型参数数据表机型B747-400B747-400B747CMBB747CMBB747SPB747SPB747FB747FB777B777B707B707B767-200B767-200长度高度经济巡航速度(km/h)935935935935935935935935935935970970885885850850最大座发动机台最大业载 （吨） 位(个)数（个）75.575.565.265.2636343.543.575.575.565.565.524243131360360400400291291291291300300375375148148214214生产国家及公司B747-400COM 70.5B747-400COM 70.519.319.370.570.519.319.370.770.719.319.356.356.319.919.970.770.719.319.363.763.718.418.444.444.412.912.948.548.515.915.94 4美波音公司美波音公司4 4美波音公司美波音公司4 4美波音公司美波音公司4 4美波音公司美波音公司4 4美波音公司美波音公司2 2美波音公司美波音公司4 4美波音公司美波音公司2 2美波音公司美波音公司第 26 页 共 28 页 B767-300B767-30054.954.915.915.98508503535225225B757-200B757-20047.347.313.513.585085022.922.9200200B737-200B737-20030.530.511.311.38568561414128128B737-300B737-30033.433.411.111.18568561515132132B737-500B737-500313111.111.185085015.515.5130130A310-300A310-30046.746.715.815.88508503434204204A310-200A310-20046.446.415.815.88288283232228228A300-600A300-60054.154.116.516.585085037.437.4274274A340-200A340-20059.459.416.716.789089050.850.8375375A320-200A320-20037.637.611.811.887687619.719.71791796.56.5 相关机型价格表相关机型价格表波音飞机目录价格表737-60051.5 -- 58.5737-70058.5 -- 69.5737-80072.5 -- 81.0737-900ER76.0 -- 87.0747-400/ -400ER234.0 -- 266.5747-400/ -400ER Freighter238.0 -- 268.0747-8293.0 -- 308.0747-8 Freighter301.5 -- 304.5767-200ER127.5 -- 139.0767-300ER144.5 -- 161.5767-300 Freighter155.0 -- 166.0767-400ER158.0 -- 173.0777-200ER205.5 -- 231.0777-200LR237.5 -- 263.5777-300ER257.0 -- 286.5777 Freighter252.5 -- 260.5787-3150.0 -- 155.5787-8161.0 -- 171.5空客飞机目录价格表（单位：百万美金）A31859.1A31970.3A32076.9A32190.3A330-200180.9A330-200F184.2第 27 页 共 28 页2 2美波音公司美波音公司2 2美波音公司美波音公司2 2美波音公司美波音公司2 2美波音公司美波音公司2 2美波音公司美波音公司2 2空中客车公司空中客车公司2 2空中客车公司空中客车公司2 2空中客车公司空中客车公司2 2空中客车公司空中客车公司2 2空中客车公司空中客车公司 A330-300A340-500A340-600A350-800A350-900A350-1000A380200.8215.5237.1249.4208.7240.6327.4第 28 页 共 28 页 。