大飞机定价 数模

1、一、问题重述国产大飞机C919项目的上马和良好进展向世界展示了中国的大国实力。近 日，中国商用飞机公司销售经理陈进表示国产大飞机C919在燃油消耗方面将比 目前所有机型减少12-15%，且在上市时定价将会低于5000万美元/架，大大低 于波音和空客同类产品的定价。这将为国产大飞机开辟广阔的市场。在放手开发研制一种新飞机时，除了技术细节外，还有很多经济问题需要回 答。其中最重要的是飞机制造商需要知道他的原始投资是否能够收回来，多久才 能收回来。这就要预测一下飞机的上市价格和市场前景。毕竟对于飞机这样庞然 大物的高科技商品，高投入，也高风险，周期还长，其销售价格除了和制造商的 制造成本、航空公司的

2、运营成本有关外，还和同类型飞机的市场前景、竞争格局 等其它因素有关。而飞机的技术数据往往能够一定程度上反映其中的一些因素。 我们要解决的问题有：1、假如你是中国商用飞机公司销售经理，请你对比同类客机B737-800，以 航空公司的角度，结合飞机运营成本，评估一下国产大飞机C919未来 的市场潜力；2、假如你是飞机制造商，请你建立一合理的数学模型为将新开发的商用飞 机预测一个较合理的价格，并结合你的模型对C919的价格进行预测；3、以你所建立的模型，结合相关数据，计算当前A380-800和B737-800两 种商用飞机的价格，并说明所建价格预测模型的优缺点；二、C919市场潜力分析2.1基本假设

3、在问题（1）中，根据实际情况，可以做如下假设来简化问题：1、本题所进行的C919和B737-800的运营成本的对比，是在同一航空公司、同 样的航线、同一时间段下所进行的，也就是说，忽略各航空公司经营模式和 飞行天气的影响。而且，两者每年起降次数、飞行时间、里程相同。2、基于第一个假设，两种机型的运营成本中的机务人员工资、航空食品配送费、 飞机入库费用、起飞场地费等都可以认为相同。3、由于B737-800的座位数为经济舱158座，C919有三种类型，分别为156座、 158座、180座，可以认为两者座位数大致相等，因此在计算平均每座运营成 本时对二者取相同座数。4、不考虑汇率变动通货膨胀。5、C

4、919的价格按官方透露的比同款机型低10%来计算。2. 2问题分析要分析一种待上市产品的市场潜力，必须剖析它与同类产品相比所具有的竞 争力，并以现有产品的市场为基础来预估。考虑到航空运输业的特点，以下我们 将通过层次分析法按重要性大小从几方面对B737-800和C919进行对比。1、运营成本：航空公司生产的最终目的是在满足经济和社会需求的同时， 最大限度地获取利润。在交通技术日新月异的今天，各航空公司不断进行价格战， 机票已不再昂贵，各航空公司尤其是民营航空公司将不断降低成本作为实现航空 公司盈利目标最主要、最可靠的途径。购买飞机时，航空公司也就更看重飞机在 使用维护阶段的经济性，其评估标准是

5、飞机的直接运营成本(DOC)和间接运营成 本(LOC)。LOC整体呈现很不稳定状态，并且变化具有非线性的趋势，不好控制， 所以对于飞机使用维护经济性分析的依据选定为直接运营成本(D0C)DOC费用包 括拥有飞机费用(折旧费、利息、保险金、高价件摊销、经营性租赁等、空勤组 费用(驾驶员成本、客舱服务员成本)、维修成本(机体维修成本、发动机维修成 本)、运行费用、燃油成本和餐食费等。考察19952000五年直属航空公司的成 本结构可以发现，航油、维修、折旧、起降、旅客服务等费用占据全部DOC成本 的80%以上。又由于假设二的存在，应认为C919和B737-800的旅客服务费、起降 费用相同，因此我

6、们选择航行飞机每航行一公里的折旧费、燃油费、维修费和航 材费的总和来作为飞机的运营成本。2、安全系数：近年来由于飞机动力、操纵、起降系统故障和极端天气引发 的飞行事故频频出现，给航空公司和社会都造成巨大损失，安全问题已成为备受 关注的问题。因此，商用飞机的安全可靠性也成为航空公司选购时的一项重要指 标。由于C919尚未投入使用，没有可统计的事故率，这里我们采用两者发动机和 操纵系统的技术可靠系数来衡量飞机的安全可靠性。3、舒适度系数：由于各机型的设计方案不同，座椅排距、行李箱高度、机 身稳定性都影响着乘客舒适度，而乘客满意度是决定该航班满客率的重要因素， 也在很大程度上影响着航空公司的盈利水平

7、。所以，飞机的舒适度也是其竞争力 的一项重要指标。我们用客舱宽度、座椅排距来衡量舒适度。4、环保系数：随着全球变暖的趋势不断加剧，环保因素越来越多地被考虑 到工业生产和服务行业中来，环保力度也逐渐成为影响企业声誉的重要因素。这 里我们用航行每公里的二氧化碳排放量来衡量飞机的污染水平。基于以上考虑，我们将采用综合指数评价方法来比较B737-800和C919对购买 方的吸引力，并以B737-800为参考量预估C919的市场潜力。对于其中最关键的运 营成本，我们在已经成熟的兰德模型下用工程估算法来计算。2.3符号约定，：飞机使用时间（单位：年）； :飞机的年折旧率；k1 :飞机折旧N年后的成本核算；

8、七：飞机N年所消耗的航空燃油总费用；七:飞机N年所维修总费用；% :飞机购入时的价格（单位：百万）；以:飞机年耗油费用（单位：百万）；f （EC）:飞机的年维修费用；bi :飞机年维修费用中各组成部分的权重；C1 :航材费用；C2 :定检维修费用；C3 :机务人员工资；E综:飞机潜力综合指数；E,:第i个评价项目的值；q :每个评价项目在该机型潜力综合指数中的权重2.4模型建立与计算1、在计算飞机运营成本时，我们采用工程估算法。这种方法的思路是将飞 机的所需的维修费用项目细分，直到最小的本费用单元，然后自下而上累加而成。单架飞机运营成本y(a,t)= % (a,t)其中飞机第t年后由于折旧而付

9、出的成本k1= a0 - a0(1- a )t ；由于发动机在老化过程中燃料利用效率会不断下架，因此k2 = (1+ a)t ；飞机第t年后共用维修费用k = f (c) X t ；3在计算f (c)时，我们采用已经很成熟并被广泛使用的兰德DAPCA IV模型。该模型是飞机发展与采购费用(Development and Procurement Costs of AircraftDAPCA )模型的最终形式。由此得到f (c)= b + b C + b C + b C = 0.00348C + 0.000475C + 0.000671C - 181.437701 12 23 3123这个模型比较

10、简单，所以我们选择使用matlab分别作出C919和B737-800 的基于时间序列的运营成本曲线，如下：图1两者运营成本比较朱成营运从图中我们可以明显看出，两者的运营成本都随使用时间非线性上升，且斜 率越来越高。但C919的运营成本较B737-800低，且随着使用时间的增长二者差 距也随之增大，优势越来越明显。另外，鉴于全球航空汽油的价格呈上涨趋势，我们在两个增长后的航油价格 下得到以下两个结果：图2 航油价格增长5%时运营成本比较olm朱成营运图3航油价格增长10%运营成本比较nolm/成营运由图可知，当航油价格上涨时，二者运营成本都会有所升高，但C919的运 营成本增长缓慢，较之B737

11、-800的优势更明显。根据官方信息，将两种机型的使用寿命都取为25年，以目前的航油价格为标准，则在使用周期内两者运营成本之比为竺61 = 1.391。设B737-800的运营313.5成本效率为1，则C919运营成本相应为1.391。2、在比较安全系数时，鉴于两者制造工艺、材料、发动机、电子设备水平 大体相当，其安全系数皆取为1。3、在比较舒适度时，乘客的舒适感与空间拥挤度有直接的关系，B737-800 的两排座椅间距为30厘米，C919的间距为32厘米，设B737-800舒适度为1， 则C919舒适度为1X 32 =1.067。304、在比较环保系数时，由于二氧化碳排放量与单位航程燃油量成正

12、比，将 B737-800的环保系数设为1，则根据官方消息，C919排放量将是类似机型的一 半，即其环保系数为2。现在来求B737-800和C919的潜力综合评价指数：E = E qE综 i i其中按照各评价项目的重要程度取q1= 0.5，q2 = 0.3，q3 = 0.1，q4 = 0.1将两者参数带入计算，得到EC919 = 1.302 , Eb737颇=1销量与购买方对飞机的满意度成正比，所以可以认为 C919的销售量为 B737-800 的 1.302 倍。目前中国大陆共有商用飞机1383架，其中波音公司占有53%的市场份 额,2009年B737-800的销售量为128架。而未来五年中国

13、（包括香港、澳门） 的商用飞机市场将有3700架的需求量.按照我们计算的结果，根据比例可粗略预 计2016年C919将拿到446架的订单。2.5模型评价与改进我们认为综合指数更有利于客观全面反映水质的标准，它体现出了一种细分 的评测标准。当然，由于掌握资料广度所限，我们略去了一些因素，采取的假设 也比较理想化，选取的加权系数是否有一个通用的标准还有待商榷。基于相似机型B737-800的销售量的市场潜力预测简单易行，但是同时也略 显粗糙，忽视了未来几年市场可能的变化和其他竞争者的可能竞争。在工程估算法的应用中，将总系统费用分解为许多用费用方程联系起来的项 目细节，其反映的因果关系与实际情况更加接

14、近，因此具有如下优点：D结果准 确，能对竞争的各个方案研究其费用差异，允许进行详细地模拟和灵敏度 分析。但是它也具有对数据要求高、估算结果很难进行评价与鉴定等明显的缺 点，有待改进。对于本问中的模型，我们可以进一步结合民机发展的几十年中市场的变化， 建立基于时间序列的市场潜力预测模型对其进行改进。三、C919价格预测3.1基本假设1、由于商飞是国防性质的国有，为简化计算，不考虑税率问题。2、C919的研制不仅仅是一种商业行为，还关系到中国的经济安全和国力体 现，牵扯到众多的不易量化的政治因素。为了简化和明朗问题，这里我 们只考虑经济因素，不考虑政治利益。3、使用兰德模型计算成本时，只考虑中美人

15、工费用差异，不考虑技术差异。4、C919的研发情况由已经生产的国产运十、新舟60来类推。3.2问题分析对于民用飞机而言，其出售时要收回研究、发展、试验与鉴定费用和生产费 用两部分（包括合理的利润）。为了收回这两部分费用，需要合理确定民机售价。 国外各大飞机制造商，如波音公司，已经根据经验数据建立了非常合理的定价系 统，能够对新机型很快作出合理的成本和价格预测，如B737-800的定价。然而， 由于缺乏自主研发新机型的经验和有效数据，国内很少有人进行这方面的研究， 对飞机的采购价格预测也仅局限于军用飞机方面。因为国产商用飞机和军机同样 存在样本小、数据少、性能因素影响显著等特点，在解决本题时我们

16、可以借鉴军 机价格预测中广泛使用的偏最小二乘法、类比法和参数法，并使用具有学习训练 功能的神经网络来对其进行修正。价格，是价值的货币表现，价值是其决定、支配因素。商用飞机是一种高成 本、高科技、高价格的产品，对于制造商来说，出售飞机的目的是收回成本即在 研制期间花费的大量研究、试验经费以及鉴定费用，而经费的多少与研制难度即 飞机的各项性能指标有关；对于购买方（主要是航空公司）来说，飞机的运营成 本及其能够带来的利益是它们关心的主要问题，而这两方面也与飞机的燃油利用 效率、载重、巡航速度等技术参数息息相关。因此，从买卖双方来说，一架飞机 的价值都应该取决于它的功能，也就是说，飞机是一种基于功能的

17、产品。由此， 我们可以认为飞机的功能是其价格的主要影响因素。因为飞机的参数可以在很大 程度上反映其功能并且比较容易量化，我们选取一些性能参数作为变量来对价格 进行多元拟合。又由于对我们来说参数与价格的关系还是一个黑箱系统，所以选 用逐步回归来选择在价格函数中真正起作用的参数。步骤如下：1、采集产品数据。搜集目前仍在使用的所有商用飞机的各项性能参数，构 建性能参数及价格数据的初始样本集。2、性能参数数据的标准化。对所采集的样本数据进行整理，对于可以直接 量化且数据大小与产品价格明显相关的功能数据，可以不予标准化直接作为输入 进行计算，对于不可量化或数据大小与产品价格不明显相关的数据，则根据产 品

18、支持的同类功能的复杂程度采用01进行打分。3、对性能参数和价格进行逐步回归。根据偏相关系数筛选出对价格影响较 大的参数，最后确定价格的表达式。4、C919价格的预测。通过分析数据，我们发现不同型号的商用飞机之间功能、造价差异较大，四 发的大型客机样本的引入可能造成预测两发的支线飞机C919的价格时的误差，因 此通过估算法，参考和C919参数相近的A320、B737-800的功能、价格关系，比较 三者之间的差异，对C919价格进行预测，并与第一种方法比较、调整。另外，对于民用飞机而言，其出售时要收回研究、发展、试验与鉴定费用和 生产费用两部分（包括合理的利润）。为了收回这两部分费用，需要合理确定

19、民 机售价。因为研究、发展、试验与鉴定费用是确定的，所以必须根据市场需求， 对生产飞机的架数做出合理的预测，以确定每架飞机的售价中要包含多少成本。 先使用第一个问题中求出的C919的预估销量，计算该销量下的成本，与售价比较； 再求在保证收回成本的情况下的销售量，结合市场发展情况进行前景分析和对策 决策。3.3模型建立与计算1、基于功能驱动的产品价格预测模型根据问题分析，首先进行样本数据采集。商用飞机的性能参数有：几何参数：翼展、机长、高度；性能参数：经济/最大巡航速度、最大起飞重量、空重、业载、发动机台 数/型号、油箱容量、载客量。由于参数过多而样本太小，需要把内在相关的一些参数进行合理的整合

20、和选 择。翼展的大小主要由飞机的气动布局决定，不能较好地反映飞机的规格，所 以我们选择机身的剖面积来衡量飞机规格；商用飞机一般在经济巡航速度上使用，所以舍弃最大巡航速度，选择经济 巡航速度作为一个参数；航空公司最注重的是经济效益，所以飞机的最大载客量和最大业载也是影 响价格的重要因素；发动机是飞机的心脏，同时发动机成本和维修费用也是飞机成本的重要组 成部分。根据空气动力学，最大推力与最大起飞重量有关，为了减少参数交互， 舍弃最大推力和发动机台数。综上所述，使用Mat lab，我们先选取飞机规格X。、经济巡航速度X、最 大载客量Xp、飞机运营成本X。、最大业载X”作为参数对价格进行逐步回归。 部

21、分波音、空客商用飞机数据见附录。当然，这些数据之间还是存在一定的交互性，而且在这个黑箱模型下我们也 不能确定它们是否对价格函数作用。下面开始进行参数的筛选：第一步，将X。、X、Xp、X。、X”作为变量对价格P进行函数拟合P = E m X ( i = vo, ve, p,。, ” )并分别求每个回归变量系数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

22、.1974）。检查它们的置信区间发现，X,。和X,e的系数的置信区间包含了零点，表明 回归变量X,。和X,e对P的影响不是太显著，因此从模型中移出次变量X,。和X,e第二步，由于Xp、X、X”对P具有较明显的显著性，用更为精确的拟合曲线来拟合表示X、X、X和P的关系。做X、X、X分别对P的散点图，pewpew再分别进行拟合，结果如下：图4以Xp为变量的P的拟合曲线图5以X。为变量的P的拟合曲线图6 以X”为变量的P的拟合曲线具体表示式如下：P = 0.50059X + 93.07P = - 0.029436X2 + 10.372X - 636.85P3 = 0.0025453Xj - 0.45

23、209Xj + 24.218乂可-123.16对应的Norm of residuals分别为46.2、43.14、45.645，由此可以看出，对 应的拟合方程的拟合程度较高。所以，用拟合后的曲线来代替原来的Xp，Xc、X 得到如下回归模型：P = E wP（1）其中 w广 0.0705 , w2 = 0.0828 , w3 = 0.1260由于计算X时涉及到一些性能参数，所以X、X”之间存在相互制约，它 们的交互作用会对P有影响，不妨简单地用X、X”的乘积代表它们的交互作用， 于是，模型（1）增加一项，得到：P = w P + w P + w P + w X X（2）1 12 23 34 C

24、W再次进行回归分析，得到修正后的系数w广 0.0676, w2 = 0.0664, w3 = 0.1044, w4 = 0.0002相关系数为0.9773，说明该模型比较接近真实情况。至此，第一个模型已经建成。将C919的有关数据代入模型（2），得到第一 个预测价格5984.65万美元。2、类比法预测C919价格由第一个模型得到的C919价格为5984.65万美元，这与官方公布的低于5000 万美元有不小的差距。经过分析数据，我们发现样本中的四发大型商用飞机与两 发支线飞机相比，各项性能参数以及价格差距甚大，从而影响了模型（2）的通 用性。我们又将一些样本外的机型的数据代入进行检验，结果显示该

25、模型用于大 型飞机的价格预测时比较准确，而在对中小型飞机的预测中结果总是偏高。因此， 我们决定使用类比法，将样本划分为四发大型飞机样本和。两发中小型飞机样 本两部分，分别进行逐步回归。其中样本。选取与C919性能参数比较接近的B737 系列和A320系列的参数。具体步骤同上节，得到两个模型，形式与（2）相同， 系数不同：模型（3） P = wP + w P + w P + w X X（适用于样本）1 12 23 34 C W模型（4） P = wP + w P + w P + w XX（适用于样本）1 12 2334 C W运用模型（4）再一次计算C919的价格，结果为4710.53万美元，这

26、一结果比 较符合官方消息。3、成本导向定价模型飞机，尤其是商用飞机的生产和销售时是以盈利为目的的高投资的商业行 为，必须结合市场来制定合适的价格保证盈利。虽然C919项目有国家资金支持， 并且考虑到一定的政治利益因素，为了国家的利益，我们还是要保证成本的回收， 甚至一定的利润。因此，在问题一预测的销量下求解C919的利润以及在预测盈亏 平衡点的价格并与上一模型求解结果比较具有很大的实际意义。1、首先我们来计算C919的研制成本。兰德DAPCA IV模型是飞机发展与采购费用（Development and Procurement Costs of AircraftDAPCA）模型的最终形式，在飞

27、机寿命周期费用分析领域有相当的影响力。按照1986年定值美元，兰德DAPCA IV模型中工时、费用的计算公式如下:H = 5.18W 0.777V 0.894 Q 0.163Ht = 721W0777V 0.696 Q 0.263H = 10.46 W 0.82 V 0.484 Q 0.641Hq = 0.767 HmCd = 33.54W o.63oV 1.3C = 968.42Wo.325V 0.822 FTA121C = 15.54W0.921V 0.621 Q 0.799Ce = 1.5480.095T + 243.25M+ 0.54$ - 2228工程工时工艺装备工时制造工时质量控制

28、工时发展支援工时飞行试验费用制造材料费用发动机制造费用研究、发展、试验与鉴定费用+生产费用C = H R + H R + H R + H R + C + C + C + C N + C 成本 E E T T M M Q Q D F M Eng Eng av式中W :空重（kg）；V :最大飞行速度（km/h）Q :产量；FTA :飞行试验机架数（一般为26架）；N财:总产量乘以每架飞机的发动机台数；T :发动机最大推力（kg）M :发动机最大Ma数；七:涡轮进口温度（K）R、R、R、R :综合费率（数值见附录） E T M QCv :航空电子设备费用C919有关参数如见附录：代入计算得到C成本

29、=g（Q）鉴于C919应用了 20%的复合材料，在上述结果上应乘以相应的软糖系数 1? 0.8 1.5? 0.2 1.1 ；考虑到中国的人工费用与欧美的差异以及目前美元价值 与1986年的差异，在上述结果上再乘以一个换算系数1.522，最终得到成本C 关成本于产量Q的函数C 成本=1.647g （Q）（5）2、当使用平均成本加成定价法对C919价格进行预估时，根据千分之一法则 将利润率定为0.1%。在这种情况下得到P加 成法=C成本：Q X （1+0.001） =6078.65 （万美元）这表明若要获得千分之一的利润，定价需在六千万美元左右。3、鉴于国产大飞机C919的研制与B737-800相

30、比基础较差，国家前期投入较 高，因此收回成本是一项重要的任务。用盈亏平衡定价法来计算保本点的价格及产量：取模型（4）中计算的较为合理的单价4710.53万美元，代入模型（5）, 得到保本点销量Q保本=673 （架）。取问题一中的销售量446架，代入模型（5），得到P保本= 6072.58 （万 美元）从这两个结果我们可以看出，想要保证回收成本，必须上调价格或者提 高产量。但是，C919的最大优势在于低廉的价格，一旦价格上调，必然影响 销量，这样很难实现盈亏平衡。因此，只能提高它的销量。由于问题一的销 量是在中国2016年市场饱和情况下求得的，所以唯一的手段是大力开拓海外 市场，而我们认为其重点

31、应该放在发展中国家等价格和燃油对决策影响巨大 的民机市场上。四、A380-800和B737-800价格计算与模型的评价改进4.1A380-800和B737-800价格的计算这两种机型的有关数据见附录。其中，A380-800属于大型商用飞机，B737-800属于中小型商用飞机，所 以应该分别使用模型（3）和（4）进行计算，结果如下：匕380 颇=333.33 （million）（真实价格为327.4million，误差为 1.81%）CB737 800 = 67.77 （million）（真实价格为72.5million，误差为6.52%） 由此看出，我们的模型能够比较准确地计算各种机型的价格。

32、4.2模型评价与改进1、模型评价多元逐步回归针对本题中众多的原始参数和较为复杂的互相制约，有效地筛 选出对价格函数起作用的变量，能够较准确地反映各有效变量对价格函数的作 用，稳定性也比较强。另外，此种算法充分利用了Matlab工具箱的优势，简化了 运算。然而，由于各参数之间存在交互作用，参数之间的关系以及各参数对价格 函数的真实作用并不明朗，这个模型实际上是在黑箱下得到的结果，对不同参数 变化的灵敏度有很大差异。另外，大型飞机与中小型飞机两类分歧显著，对结果 的精确度造成了极大的影响。类比法正是针对这个问题而补充的，它将待求对象分为两类，分别通过逐步 回归得到不同的模型，提高的预测的准确度，更

33、为科学与合理。成本导向定价法的引入进一步解决了题干中的成本回收和盈利问题，为C919 的生产、上市提供了合理、可行的建议，提高了本文结论的实用价值。2.使用最小二乘支持向量机对模型进行优化逐步回归分析通常需要大样本容量空间的支持，并且模型泛化能力差，而本 题中能够使用的样本经过分类后数量过少，不足以保证模型的精准和稳定。为了 解决本模型在小样本下的失准问题，我们又选择了能够在Matlab工具箱中方便使 用的最小二乘支持向量机作为预测工具，对模型进行优化。这种算法基于以小样 本数据为主要研究对象的统计学习理论，采用结构风险最小化准则来代替传统的 经验风险最小化准则，根据有限的样本信息在模型的复杂

34、度和模型的学习能力之 间寻求最佳折衷，同时最小化VC维和经验风险，以获取期望风险的最小化,它不像 神经网络等其他预测方法一样需要大样本空间的支持，并最终转化为二次寻优问 题，可获得问题的全局最优解，较好地解决了易陷入局部最优的问题。最小二乘支 持向量机优化指标采用了平方项，只有等式约束，而没有标准支持向量机的不等 式约束，从而推出不同的一系列等式约束，将二次规划问题转变成线性方程组的 求解，简化了计算复杂性，因此具有较快的运算速度。模型建立如下：p(X)= E d K(x ,x.) + b式中，X为飞机的参数数值，K(x ,x.)为核函数，l .为待定参数，其值可由下 式计算：式中p =(p,

35、p,?,p)t; ei =(1,1,?,1)t; Q= f (x )f(x)= K(x,x )； x、x .为2iijiJi Ji J 飞机参数和飞机价格。具体代码和核函数求解见附录。经过多次学习训练后得到C919的预测价格为4967.56万美元，相比前述模 型更为可信。五、参考资料【1】数学模型(第三版)姜启源谢金星叶俊高等教育出版社2003 年【2】基于MATLAB和LINGO的数学实验肖华勇 西北工业大学出版社2008 年【3】人工神经网络原理及仿真实例(第二版)高隽 机械工业出版社 2007 年【4】飞行费用与效能分析【5】基于MATLAB的多元非线性回归模型董大校2009年【6】基于

36、功能驱动的产品价格预测 黄训江六、附录6.1飞机运营成本matlab代码%商用飞机的运营成本模型%前提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;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)/

37、(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

38、，);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.2 7.0 4.0 4.5 5.9 5.6 4.8 3.9;%x2=9 20 18 33 31 13 25 30 5 47 25 11 23 35 39 21 7 40 35 23 33 27 34 15;%x3=6.1 6.4 7.4 6.7 7.5 5.9

39、 6.0 4.0 5.8 8.3 5.0 6.4 7.6 7.0 5.0 4.4 5.5 7.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);%f

40、printf(%6.4fn,b);%fprintf(回归方程的统计量：n);%5fprintf(%6.4fn,stats);6.2问题二中逐步回归matlab代码clcclear;%以下是数据的初始化%飞机长度length=70.570.544.448.570.756.354.947.370.730.563.733.4;%飞机高度high=19.319.312.915.919.319.915.913.519.311.318.411.1;for i=1:12V(i) =length(i)*high(i);% 飞机的横截面 end%经济巡航速度velocity=935935850850856935

41、935935856 ;935885850%最大载客量 passenger=291291110134375375400300214150;200330%运营成本cost=149 154 168.6 172 140 123 156 106 110 197 133 118.4 ;%受发动机台数制约的最大业重engine=3065.239.5406365.5251415582017;%发动机台数的初始化，此处可根据需要修改为2% fadongji=2;%12个样本飞机的售价(单位：million dollars)%丫矩阵包含着所有4发商用飞机的价格y=ones(1,12);%用来储存每个样本的售价&%

42、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)

43、二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);%Z矩阵代表的样本皆为2发样本，由于样本值浮动较大所以此处我们皆取平均值% z=ones(1,12);%737-600%z(1)=73.1436;%737-700%z(2)=74.5521;%73

44、7-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),p

45、lot(passenger,y,o);subplot(1,6,2),plot(cost,y,+);subplot(1,6,3),plot(engine,y,*);%1.首先对5个因素进行线性回归alpha=0.05;%设置置信区间X=V;velocity;passenger;cost;engine;%根据不同的发动机规格选择不同的样本值进行模拟if(fadongji=4)y1=(y);elsey1=(z);end%利用线性回归求解的回归系数b,bint,r,rint,stats=regress(y1,X,alpha);fprintf(线性回归系数：n);fprintf(%6.4fn,b);fp

46、rintf(回归方程的统计量：n);fprintf(%6.4fn,stats);subplot(1,6,6),rcoplot(r,rint);% 绘出时序残差图%2.进行非线性回归for j=1:12finalpassenger(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%非线性化后回归系数的求解过

47、程beta=0.05;e=finalpassenger;finalcost;finalengine;b1,bint1,r1,rint1,stats1=regress(y1,e,beta);fprintf(线性回归系数：n);fprintf(%6.4fn,b1);fprintf(回归方程的统计量：n);fprintf(%6.4fn,stats1);%进行系数交互for k=1:12new1(k)=finalengine(k)*finalcost(k);end%最终回归系数的求解final=finalpassenger;finalcost;finalengine;new1;b2,bint2,r2,

48、rint2,stats2=regress(y1,final,beta);fprintf(线性回归系数：n);fprintf(%6.4fn,b2);fprintf(回归方程的统计量：n);fprintf(%6.4fn,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(1

49、2)=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.4fn,C919);%根据模型预测737-800的价格%passenger(12)=186;cost(12)=156.9;engine(12)=2

50、6;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(1

51、2)+b2(4)*new1(12);fprintf(根据模型，我们预测737-800的价格为二%6.4fn,P737);6.3问题三中最小二乘向量机求核矩阵代码x =70.570.570.756.3 70.763.744.448.554.947.330.533.431 46.7 46.454.159.4;19.319.319.319.9 19.318.412.915.915.913.511.311.111.1 15.8 15.816.516.7;935935935935935935885850850850856856850 850828850850;3604002912913753751482

52、14225200128132130 204228274375;125154135172125123156178148197179165176 154146143153;4444424222222 2222y=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-