全国大学生数学建模竞赛B题5335.pdf

《全国大学生数学建模竞赛B题5335.pdf》由会员分享，可在线阅读，更多相关《全国大学生数学建模竞赛B题5335.pdf（20页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 “互联网+”时代的出租车资源配置 摘 要 随着“互联网+”时代的到来，针对当今社会“打车难”的问题，多家公司建立了打车软件服务平台，并推出了多种补贴方案，这无论是对乘客和司机自身需求还是对出租车行业发展都具有一定的现实意义。本文依靠ISM解释结构、AHP-模糊综合评价、价格需求理论、线性规划等模型依次较好的解决了三个问题。对于问题一求解不同时空出租车资源“供求匹配”程度的问题，本文先将ISM模型里的层级隶属关系进行改进，将影响出租车供求匹配的 12 个子因素分为时间、空间、经济、其它共四类组合，然后使用经过改进的AHP-模糊综合评价方法建立模型，提出了出租车空载率这一指标作为评价因子的方案，

2、来分析冬季某节假日哈尔滨市南岗区出租车资源“供求匹配”程度。通过代入由 1-9 标度法确定的各因素相互影响的系数，得出各个影响因素的权重大小，利用无量纲化处理各影响因素，得出最终评判因子为 0.3062，根据“供求匹配”标准，得出哈尔滨市南岗区出租车资源“供求匹配”程度处于供需合理状态的结论。同理，也得到了哈尔滨市不同区县、不同时间的供求匹配程度，最后作出哈尔滨市出租车“供求匹配”程度图。对于问题二我们运用价格需求理论建立模型，以补贴前后打车人数比值与空驶率变化分别对滴滴和快的两个公司的不同补贴方案进行求解，依次得到补贴后对应的打车人数及空驶率的变化，再和无补贴时的状态对比，最后得出结论：当各

3、公司补贴金额大于 5 元时，打车容易，即补贴方案能够缓解“打车难”的状况；当补贴小于 5 元时，不能缓解“打车难”的状况。对于问题三，在问题二的模型下，建立了一个寻找最优补贴金额的优化模型，利用 lingo 软件1进行求解算出最佳补贴金额为 8 元，然后将这个值带入问题二的模型进行验证，经论证合理后将补贴金额按照 4 种分配方案分配给司机乘客。关键词：ISM解释结构模型；AHP-模糊综合评价；价格需求理论；线性规划 一 问题重述 交通是社会生活众多产业当中的一项基础产业，不但和社会的经济发展关系紧密，与人们的生活也是息息相关。而出租车作为交通工具中的重要组成部分，可以给人们的生活带来更便捷的服

4、务。所以无论是国内还是国外，对于与出租车相关的问题都有较深入的研究。作为居民日常出行的交通工具，出租车在人们生活中发挥着重要的作用。然而由于时间、空间等因素，导致出租车行业面临着巨大的挑战，与此同时，也会出现“打车难”的现象。但这也正促进了依托互联网建立的打车软件的发展以及多种出租车补贴方案的出现。当今，“互联网+”发展迅速，所以研究互联网与出租车有关的问题是很有意义，本文要求搜集相关数据建立数学模型求解下列问题：问题一：建立合理的指标，分析在不同时间和空间条件下出租车资源的“供求匹配程度”。问题二：分析各公司的出租车补贴方案是否对“缓解打车难”有所帮助？问题三：若要创建一个新的打车软件服务平

5、台，应该设计什么样的补贴方案，并论证方案的合理性。二 模型假设 假设 1：城市中不出现大量的人口迁入与迁出。假设 2：城市中出租车的数量短期内不会发生变化。假设 3：城区面积不发生大规模扩展。假设 4：城市道路发展程度不会发生大变化。假设 5：手机打车软件的使用者年龄主要集中在 1835 周岁。假设 6：其它交通工具发展水平不变。假设 7：城市人均收入短期内不变。三 符号说明 iW 影 响“供 求 匹 配”程 度 第 一 层 因 素 的 权 重(1,24)i iW 影响“供求匹配”程度第二层各因素的权重(1,24)i iW 修正后的影响“供求匹配”程度第二层各因素的权重(1,24)i M 用于

6、评价“供求匹配”程度的评价因子 max 判断矩阵的最大特征值 B 综合评判集 iU 层次分析结构中第一层指标(1,24)i 1P 价格补贴之前打车人数 2P 价格补贴之后打车人数 E 出租车日均载客里程 n 出租车辆数 出租车出勤率 l 每辆车平均的日行里程。S i 乘客不满意度 四 问题分析 问题一的分析 对于问题一，要想得到出租车资源的“供求匹配”程度，首要的问题是建立一个合理的评价指标。通过对影响出租车与乘客供求关系的广义因素进行分析，每种大的前提下又细分为其它的影响因子，也就是构造了两个层次，再将最底层的影响因子利用ISM解释结构模型2进行归类。利用层次分析-模糊综合评价模型对得到的归

7、类进行分析和求解，得到综合评判集合，然后考虑结合一种出租车供需合理的标准，例如空载率这一指标对供求匹配程度进行分析。最后结合实际着重研究不同时间和空间前提下城市出租车资源的“供求匹配程度”。问题二的分析 对于问题二，求各公司的补贴方案对“缓解打车难”是否有帮助，考虑到不同补贴方案归根到底是补贴金额的不同，因此考虑寻找一种补贴金额与打车难的关系，通过逆向思考，补贴金额可以等效看为出租车价格降低的金额，出租车价格变化与打车人数需求之间构成价格需求，于是可利用价格需求理论模型对此进行分析，但是单一的打车人数多少不足以表示打车是否困难，于是考虑增加一个空载率指标与打车人数结合表示打车是否困难，最后评判

8、打车困难时，由于打车难这是一个不可量化指标，因此要想得出打车难是否有缓解只能先建立一个标准，然后将价格需求理论模型的求解结果带入该标准。即可知道各公司的补贴整车对打车难是否有帮助。问题三的分析 由于问题三是设计补贴方案，而问题二里我们建立了价格需求理论模型，求解了各个公司不同补贴方案对打车难缓解的影响，于是我们在问题三中通过建立一个补贴金额与乘客满意程度的双优化模型来设计一个补贴方案，然后利用问题二的求解结果对设计方案进行论证。五、模型的建立与求解 5.1 问题一模型的建立与求解 5.1.1 利用ISM模型对影响因素分组 由于出租车资源供求匹配关系受到较多因素影响，其中很多因素又相互包含，必须

9、全面考虑各个因素。所以我们采用 ISM 模型对相应影响因素进行分组、归类，使问题简化，方便求解。ISM模型是.J NWarfield于 1973 年为了分析复杂的社会经济问题而提出的解释结构模型，是静态的定性模型,其特点是能够将复杂的系统逐级分解成若干个子系统。为了分析出租车资源的供求匹配程度，我们考虑影响出租车与乘客供求关系的一些主要因素如表1。表 1 出租车与乘客供求关系主要影响因素 1A 市民人均收入 2A 人口密度 3A 道路拥堵程度 4A 是否是高峰期 5A 是否节假日 6A 出租车数量 7A 天气情况 8A 私家车数量 9A 季节 10A 城市交通发展水平 11A 乘客出行紧急程度

10、 12A 城市旅游吸引能力 然后分析这些因素互相之间是否有关系，用 0 表示相互之间无影响，用 1表示相互之间有影响，从而得到相互影响关系的邻接矩阵如下：1A 2A 3A 4A 5A 6A 7A 8A 9A 10A 11A 12A 1A 0 0 0 0 0 0 0 1 0 0 0 0 2A 1 0 1 0 1 1 0 1 0 0 0 0 3A 1 0 0 0 1 1 0 0 0 0 0 0 4A 0 1 1 0 1 1 0 1 1 1 1 1 5A 0 0 0 0 0 0 0 0 0 0 0 0 6A 0 0 0 0 0 0 0 0 0 0 0 0 7A 0 0 0 0 0 0 0 0 0 0

11、 0 0 8A 0 0 0 0 0 0 0 0 0 1 1 0 9A 0 0 0 0 0 0 0 0 0 0 1 0 10A 0 0 0 0 0 0 0 0 0 0 0 0 11A 0 0 0 0 0 0 0 0 0 0 0 0 12A 0 0 0 0 1 0 0 0 0 0 0 0 由于此矩阵中影响因素较多，所以运用Matlab软件对邻接矩阵进行求解，程序见附录一，得 10000001001011101101001010101101001011111101111100001000000000000100000000000010000000000001001000000000101000000

12、0000100000000000010000010000001K 可达矩阵可说明两个因素之间是否存在链接路径，并能清楚说明两因素之间影响程度。定义集合 P I为可达矩阵中要素I一行中值为 1 的元素所在行元素集合，H I为可达矩阵中I这列值为1对应的列元素集合，当()()()H IP IH I时，I即为该层元素，然后剔除矩阵中的W元素，进行下一层元素计算,可以得到最终的分组：1168,UA A A;22310,UA A A;371112,UA AA;4459,UA A A。通过分析每组所包含的因素，我们发现分组1U里面所包影响因素可理解为经济水平对出租车供求关系的影响，2U里面所包影响因素可以

13、看作为不同空间对出租车供求关系的影响，3U里面所包影响因素可以看作不同时间对出租车供求关系的影响，结果如表2。表 2 影响因素分组关系表 经济水平因素1u 市民人均收入1A 出租车数量6A 私家车数量8A 空间因素2u 道路拥堵程度3A 人口密度2A 城市交通发展水平10A 季节9A 时间因素3u 是否是节假日5A 是否是高峰期4A 其它4u 乘客出行的紧急程度11A 旅游吸引能力12A 天气状况7A 5.1.2 问题一模型的建立 我们从出租车空载率角度考虑出租车资源的“供求匹配”程度问题，当出租车空载率过低时，说明打车的人少，出租车量小于乘客的需求；当空载率过高时，表明打车的人较多，出租车量

14、大于乘客需求，出租车空载率能很好地反映出租车与人之间“供求匹配”程度。所以我们选取空载率这一指标作为模型最终评价因子分析，来分析不同时间和空间出租车资源的“供求匹配”程度问题。我们利用()AHP模糊综合评价方法建立模型，首先，利用()AHP构造出一个层次分析模型，指标评价体系结构图如图1。图 1 层次分析结构图 1）第一、二层权重集的确定 第一层包括 4 个因素，即1234(,)UU U U U，利用AHP层次分析法比较几种指标间的关联度如图 3。表 3 第一层因素间关联度 1U 2U 3U 4U 1U 1 1/4 1/4 1/4 2U 4 1 2 4 3U 4 1/2 1 2 4U 1/2

15、1/4 1/2 1 然后确定第一层指标权重，利用 19 标度法求解判断矩阵，构造第一层的评判矩阵Y具体形式如下：11121314212223243132333441424344u u u u11/41/41/4u u u u4124=u u u u41/2121/21/41/21u u u uY 其中：12u表示1u与2u之间的关联度。之后求解第一层的权重集W=1W,2W3W,4W。方法如下：先计算判断矩阵Y每行元素的乘积iQ，再对iQ求五次方根，得到 1234,0.232,0.821,1.736,0.427WW W W W 归一化处理：利用公式 41/iiiiWWW对1234,WW W W

16、W做归一化处理，得到1,234(,)(0.100,0.462,0.305,0.133)WW W W W即为所求特征向量。一致性检验：为了说明以上所求得特征向量是否能够合理的分配权重，需要进行一致性检验，方法如下：通过公式4max14iiYWW求得判断矩阵的最大特征值，得到最大特征值max5.024。然后利用公式CICRRI和max1nCIn，通过代入数据1.12RI,4n，max5.024，得到0.0350.1CR，这就表明评判矩阵Y具有很好的一致性，所以1234,WW W W W中的各项均可以作为U 的权重系数。同第一层权重的求解过程，对于第二层指标，由于第一层的每一个因素都包含123,ii

17、iiuuuu3 个因素，于是得到第二层级的各项权重集：1234=0.2860.0810.5670.0860.1770.764=0.0760.6810.267=0.1820.0510.727WWWW 2）确定综合评判结果 根据模糊综合评价法可知，综合判别公式12(,)mBW Yb bb，其中41()(1,2)jiijibw yjm，将数据代入计算公式41()iiiBW W得到所研究问题的综合评判结果：0.1160.3040.580B。3）综合评判矩阵B的修正 影响出租车供求匹配的具体因素有些有具体实际数据，各种数据的单位并不统一，不容易量化，而另外一些因素如：旅游吸引能力、乘客出行的紧急程度等影

18、响因素是模糊的量。为了使各项因素具有可比性，我们将这些有实际数据的影响因素进行标准化处理，消除量纲差异。利用如下公式 (1,2,;1,2,)ijjijjin jmxxxs （1）211,1111,2,jjnnijjijiinnjmxxsxx （2）将市民人均输入，人口密度，出租车数量，私家车数量数据分别按照上式标准化，但是得到的数据并不全是处于 0,1之间，所以应用平移-极差变换使得所有的原始数据都在 0 1，之间。得出量化结果：人均收入 11A，人口密度 50.87A，出租车数量 60.96A，私家车数量 81A。对于模糊的影响因素，我们用 0 1，之间的数来描述，例如：0.2 对于天气状况

19、来说意味着“差”，而 0.8 意味着“好”。其他模糊因素同理，得具体情形如表 4。表 4 模糊影响因素处理结果 利 用 上 表 对 12 个 因 素 的 定 量 分 析，我 们 可 以 对 第 二 层 级 权 重123,(1,2,3,4)iiiiWW WWi进行修正，使之与这12 个因素联系更加紧密，使结果更准确，第二层权重的修正如下：1123(0,2860.0810,567)WAAA 2456=(0.0860.0170.764WAAA ）37890 0760.6810.267WAAA，41011120.1820.0510.727WAAA 4)最终评判标准的确定 对于出租车的资源供求匹配是否合

20、理，我们采用出租车空驶率来进行评判，空驶率与出租车供需关系如下表所示：表 5 城市出租空载与交通供求关系4 空载率 小于25%30%左右 大于 40%匹配程度 明显供不应求 基本饱和，供需合理 明显供过于求 为了最终得出评判标准，我们将已划分好的空驶率写成一个3 1阶的矩阵P，即 道路拥堵程度4A 拥堵（0.2）顺畅(0.8）乘客出行的紧急程度7A 不着急（0.2）着急（0.8）天气状况9A 好（0.2）差（0.8）城市交通发展水平6A 差（0.2）好（0.8）是否是节假日10A 否（0.2）是（0.8）季节11A 春秋（0.2）冬夏（0.8）旅游吸引能力8A 弱（0.2）强（0.8）0.25

21、0.300.40P 于是我们结合已经求出的评判矩阵B，可得出最终的评价结果公式:MB P 评价因子M与空载率是同一量级，我们将M与空载率标准（表 5）进行比较，从而得出出租车资源“供求匹配”的程度。到此，问题一的模型建立完成。5.1.3 问题一模型的求解 本文以哈尔滨市为例，根据相关数据进行求解哈市出租车资源的“供求匹配”程度问题。我们从哈尔滨统计局获取哈市市区包括下辖区各县相关数据见附件1。为了满足从空间角度考虑问题，首先，我们将哈市分为中心区与边缘区，中心区包括市辖区：南岗区、道里区、道外区、香坊区、平房区、松北区、呼兰区、阿城区。边缘区包括下瞎县城：宾县、延寿县、通河县等。时间角度我们主

22、要从大的方向考虑，比如：季节，是否是节假日，是否是高峰期等。我们先以哈市冬季南岗区某节假日为例，求解出租车资源“供求匹配”程度问题。南岗区相关数据见附件 1。我们将评价所需的 12 个南岗区的相关数据代入所建建立的模型中，得到量化后的指标为：0.92,0.87,0.2,0.8,0.8,0.96,0.8,0.92,0.8,0.8,0.8,0.8iA 评判矩阵为：0.2630.0710.1130.0690.1420.7330.6080.6270.2140.1460.0110.582iW 利用 Matlb 求解(程序见附录四)，得到0.2436,0.2639,0.4152B,最后计算出的评价因子0.

23、3062MBP。带入表 5 中进行比较，比较结果：0.300.30620.40，得出匹配程度：供需合理。从而得出结论：哈市冬季南岗区某节假日出租车资源“供求匹配”程度处于供需合理状态。类似的，利用上述方法可以求解出其他区和所辖县的“供求匹配”程度状况。画出哈市出租车资源“供求匹配”程度图如图 2。对于问题二中“打车难”这个问题，经过分析，发现乘客打车的难易程度最终都可以反映在补贴之后打车人数和出租车空载率这两个指标上。因为城市出租车数量短时间内是不会发生大的变化的，而打车的人数却是一个变化的数值，打车的人多了，打车就变得相对困难。出租车的空载率变小，人们打车就变得相对困难，反之亦然。所以我们用

24、补贴方案实施前后，打车的人数和出租车空驶率的变化两个指标来反映乘客打车的难易程度，由此间接评判出各公司的出租车补贴方案对缓解打车难的现象是否有帮助。虽然补贴的方案不尽相同，但其本质都是补贴一定数额的钱给相应的出租车司机和乘客，我们将不同补贴方式的差异转化为补贴金额的差异，由于补贴金额是直接影响打车难的关键因素，所以我们假设在这个期间不会有其他因素造成打车人数和空驶率的改变，打车人数和空驶率的改变只受补贴金额的影响。5.2.2 问题二模型的建立 我们借鉴价格理论模型5进行建模，价格需求理论是用于研究价格与需求人数的关系的方法，所以在这里我们将其应用于出租车价格变化与打车人数的关系研究上，但是本题

25、中并没有价格变动这一指标，不过，因为补贴方案的不同最终影响的是补贴金额的不同，所以我们这里用补贴金额代替价格变动这一指标。价格需求理论指出价格与需求成负指数关系，所以我们建立如下方程 1220.480733(5.117)(1)xPPPr er （3）其中：1P为价格补贴之前打车人数，2P为价格补贴之后打车人数，r为受出租车价格影响打车的人数占总人数的比例。在这里受出租车价格影响，打车的人数考虑为使用手机软件打车的人数，而这部分人年龄分布大概在 18-35 之间，我们通过查找哈尔滨统计局的各年龄人口数据（见附件 1），估计出哈尔滨地区26.7%r。为了更好地体现补贴前后打车人数变化情况，我们采用

26、补贴前后打车人数比值12PP作为评价打车人数变化的指标，根据公式（3）得方程：120.48073(50117)(1)xPrPre （4）利用 Matlab 画出函数图象如图 4。（程序见附录二）图 4 打车人数变化率与补贴金额的关系 空载率与出租车价格变化的关系式为：1ijEnlppi （5）E为出租车日均载客里程；n为出租车辆数；为出租车出勤率；l为平均每辆车的日行里程。将式（4）带入（5）得到：空驶率与出租车价格的关系式为：0.48073(5.117)1(1)xErnlrei （6）我国汽车协会统计结果得出，当空载率大于 40%时，供过于求；当空载率小于 25%时，供大于求；当处于 25%

27、40%之间时，处于供求平衡状态。供过于求，车多，可以认为打车容易。供小于求，车少，可以认为打车难。当处于 30%40%之间时，处于供需平衡状态，我们不考虑空载率对打车难的影响，只考虑打车人数变化关系。最后我们综合考虑补贴前后打车人数比值12PP和空载率i变化情况，建立综合评价标准如表。表 6 综合评价标准 打车难 打车容易 打车人数前后比率12PP 1 1 空载率 i 30%40%5.2.3 问题二模型的求解 根据参考资料(见附件 1)，我们得到各公司补贴方案，考虑到我们采用的是补贴金额数来进行求解的，所以我们对各种补贴方案数据进行处理得到如下的表格。表 7 各公司补贴方案（单位：元）公司 补

28、贴阶段 滴滴打车 快的打车 一阶段补贴金额 10 10 二阶段补贴金额 10-15 11 三阶段补贴金额 12-20 13 四阶段补贴金额 6-15 10 五阶段补贴金额 3-5 3-5 六阶段补贴金额 0 0 将各个阶段的补贴金额数分别带进公式（4）（6）中，利用 Matlab 软件求解（程序见附录三）得：表 8 滴滴打车补贴方案评价指标表 评价标标 补贴阶段 打车人数变化率12PP 空载率i 一阶段补贴 0.7585 0.8604 二阶段补贴金额 0.7407 0.8717 三阶段补贴金额 0.7344 0.8756 四阶段补贴金额 0.7531 0.8638 五阶段补贴金额 1.1898

29、 0.5883 不补贴 3.8529-1.0953 表 9 快的打车补贴方案评价指标表 评价标标 补贴阶段 打车人数变化率12PP 空载率i 一阶段补贴金额 0.7585 0.8604 二阶段补贴金额 0.7488 0.8665 三阶段补贴金额 0.7390 0.8727 四阶段补贴金额 0.7585 0.8604 五阶段补贴金额 1.1898 0.5883 不补贴 3.8579-1.0953 5.2.4 结论 根据上述求解结果，对不同补贴方案，打车人数比值12PP与空载率i都不相同，且当补贴金额大于 5 元时，12PP均小于 1 且i均大于 40%，而当补贴金额小于 5元时12PP大于 1，

30、且空驶率大幅下降。据此根据表 6，我们得出结论：各公司的出租车补贴方案在一定程度上对“缓解打车难”是有帮助的，当补贴方案对应的补贴金额大于 5 元时，乘客打车变得容易，缓解了“打车难”的状况。当各个公司补贴方案对应的补贴金额小于 5 元时甚至不补贴时，打车较为困难，对缓解“打车难”没有帮助。5.3 问题三的模型建立及求解 5.3.1 线性规划模型的建立 对于补贴方案的合理性，考虑问题二的求解方法，首先我们应用同样的方式，将补贴方案的差异，最终转化为补贴金额的差异。当我们创建一个打车平台并且自己设计补贴方案时，设计的方案要达到的目的是给司机和乘客补贴最少的钱，并使顾客达到最大满意度。于是我们可以

31、将问题三中求解补贴方案合理性的目标转化为求解补贴金额最小以及不满意度最低的优化模型。对于满意度，我们认为空驶率越高,打车越容易,相应的满意度越高.。在此我们假定乘客满意度就是由空驶率唯一确定,并且与空驶率成正相关。但是由于我们要求解的优化模型是求解最小值，故而我们把求最大满意度问题转变为求解最低不满意度问题。我们用 01 之间的数来描述乘客不满意度，最高为 1，最低为 0。当出租车空载率为 30%左右时，打车容易，供需基本处于平衡状态；超过 40%时，供过于求，低于 25%时，供不应求，打车困难。不满意度与空驶率对应关系如下表：表 10 空驶率与不满意度对应表 空驶率 40%30%0;if i

32、sequal(R,Rnew)Rnew N=N+1 break end R=Rnew;N=N+1;end success=xlswrite(E:R.xls,Rnew)附录二：%画打车人数变化率与补贴金额函数图象程序代码 x=linspace(0,20,1000);N=0.267*exp(-0.48073*(x-5.117)+1-0.267;plot(x,N,r)title(打车人数变化率与补贴金额函数图象);xlabel(补贴金额/元);ylabel(打车人数变化率);hold on 附录三：%滴滴打车补贴求空载率程序代码 x=10,12.5,16,10.5,4,0;i=1-307.9*1000

33、0*0.267*(exp(-0.48073*(x-5.117)+1-0.267)/(336*0.9*16136)i=0.8604 0.8717 0.8756 0.8638 0.5883 -1.0953%滴滴打车补贴求人数变化率程序代码 x=10,12.5,16,10.5,4,0;G=0.267*exp(-0.48073*(x-5.117)+1-0.267 G=0.7585 0.7407 0.7344 0.7531 1.1898 3.8579%快滴打车补贴下求空载率程序代码 x=10,11,13,10,4,0;i=1-307.9*10000*0.267*(exp(-0.48073*(x-5.11

34、7)+1-0.267)/(336*0.9*16136)i=0.8604 0.8665 0.8727 0.8604 0.5883 -1.0953%快滴打车补贴下求人数变化程序代码 x=10,11,13,10,4,0;G=0.267*exp(-0.48073*(x-5.117)+1-0.267 G=0.7585 0.7488 0.7390 0.7585 1.1898 3.8579 附录四：%求解哈市冬季南岗区某节假日出租车资源“供求匹配”程度程序代码 W1=0.100,0.462,0.305,0.133;W2=0.263,0.071,0.113;0.069,0.142,0.733;0.608,0.

35、627,0.214;0.146,0.011,0.582;W2(1,:)*W1(1)ans=0.0263 0.0071 0.0113 W2(2,:)*W1(2)ans=0.0319 0.0656 0.3386 W2(3,:)*W1(3)ans=0.1854 0.1912 0.0653 U=0.0263,0.0071,0.0113;0.0319,0.0656,0.3386;0.1854,0.1912,0.0653;B=sum(U)W1=0.100,0.462,0.305,0.133;W2=0.263,0.071,0.113;0.069,0.142,0.733;0.608,0.627,0.W2=0.

36、2630 0.0710 0.1130 0.0690 0.1420 0.7330 0.6080 0.6270 0 W1=0.100,0.462,0.305,0.133;W2=0.263,0.071,0.113;0.069,0.142,0.733;0.608,0.627,0.214;0.146,0.011,0.582;W2(1,:)*W1(1)；W2(2,:)*W1(2)；W2(3,:)*W1(3)；U=0.0263,0.0071,0.0113;0.0319,0.0656,0.3386;0.1854,0.1912,0.0653;B=sum(U)B=0.2436 0.2639 0.4152 P=0.25;0.30;0.40;M=B*P M=0.3062 附录五：%求解补贴金额最优解的代码 Moudel:Min x i=1-307.9*10000/366*0.9*13636*（H1/H2）；S（i）=1/(1.0165+328.1908exp(23.973*i+0.25*23.973);H1=0.267*H1*exp(-0.48073(x-5.117)+0.733*H2;i0.3;i0;H1/H20;End