2015全国大学生数学建模竞赛B题

“互联网+”时代的出租车资源配置

摘要

随着“互联网+”时代的到来，针对当今社会“打车难”的问题，多家公司建立了打车软件服务平台，并推出了多种补贴方案，这无论是对乘客和司机自身需求还是对出租车行业发展都具有一定的现实意义。本文依靠ISM解释结构、AHP-模糊综合评价、价格需求理论、线性规划等模型依次较好的解决了三个问题。

对于问题一求解不同时空出租车资源“供求匹配”程度的问题，本文先将ISM模型里的层级隶属关系进行改进，将影响出租车供求匹配的12个子因素分为时间、空间、经济、其它共四类组合，然后使用经过改进的AHP-模糊综合评价方法建立模型，提出了出租车空载率这一指标作为评价因子的方案，来分析冬季某节假日市南岗区出租车资源“供求匹配”程度。通过代入由1-9标度法确定的各因素相互影响的系数，得出各个影响因素的权重大小，利用无量纲化处理各影响因素，得出最终评判因子为0.3062，根据“供求匹配”标准，得出市南岗区出租车资源“供求匹配”程度处于供需合理状态的结论。同理，也得到了市不同区县、不同时间的供求匹配程度，最后作出市出租车“供求匹配”程度图。

对于问题二我们运用价格需求理论建立模型，以补贴前后打车人数比值与空驶率变化分别对滴滴和快的两个公司的不同补贴方案进行求解，依次得到补贴后对应的打车人数及空驶率的变化，再和无补贴时的状态对比，最后得出结论：当各公司补贴金额大于5元时，打车容易，即补贴方案能够缓解“打车难”的状况；当补贴小于5元时，不能缓解“打车难”的状况。

对于问题三，在问题二的模型下，建立了一个寻找最优补贴金额的优化模型，利用lingo软件[1]进行求解算出最佳补贴金额为8元，然后将这个值带入问题二的模型进行验证，经论证合理后将补贴金额按照4种分配方案分配给司机乘客。关键词：ISM解释结构模型；AHP-模糊综合评价；价格需求理论；线性规划

一问题重述

交通是社会生活众多产业当中的一项基础产业，不但和社会的经济发展关系紧密，与人们的生活也是息息相关。而出租车作为交通工具中的重要组成部分，可以给人们的生活带来更便捷的服务。所以无论是国还是国外，对于与出租车相关的问题都有较深入的研究。作为居民日常出行的交通工具，出租车在人们生活中发挥着重要的作用。然而由于时间、空间等因素，导致出租车行业面临着巨大的挑战，与此同时，也会出现“打车难”的现象。但这也正促进了依托互联网建立的打车软件的发展以及多种出租车补贴方案的出现。当今，“互联网+”发展迅速，所以研究互联网与出租车有关的问题是很有意义，本文要求搜集相关数据建立数学模型求解下列问题：

问题一：建立合理的指标，分析在不同时间和空间条件下出租车资源的“供求匹配程度”。

问题二：分析各公司的出租车补贴方案是否对“缓解打车难”有所帮助？

问题三：若要创建一个新的打车软件服务平台，应该设计什么样的补贴方案，并论证方案的合理性。

二模型假设

假设1：城市中不出现大量的人口迁入与迁出。

假设2：城市中出租车的数量短期不会发生变化。

假设3：城区面积不发生大规模扩展。

假设4：城市道路发展程度不会发生大变化。

假设5：手机打车软件的使用者年龄主要集中在18～35周岁。

假设6：其它交通工具发展水平不变。

假设7：城市人均收入短期不变。

三符号说明

1,24)

1,24)

四问题分析

问题一的分析

对于问题一，要想得到出租车资源的“供求匹配”程度，首要的问题是建立一个合理的评价指标。通过对影响出租车与乘客供求关系的广义因素进行分析，每种大的前提下又细分为其它的影响因子，也就是构造了两个层次，再将最底层的影响因子利用ISM解释结构模型[2]进行归类。利用层次分析-模糊综合评价模型对得到的归类进行分析和求解，得到综合评判集合，然后考虑结合一种出租车供需合理的标准，例如空载率这一指标对供求匹配程度进行分析。最后结合实际着重研究不同时间和空间前提下城市出租车资源的“供求匹配程度”。

问题二的分析

对于问题二，求各公司的补贴方案对“缓解打车难”是否有帮助，考虑到不同补贴方案归根到底是补贴金额的不同，因此考虑寻找一种补贴金额与打车难的关系，通过逆向思考，补贴金额可以等效看为出租车价格降低的金额，出租车价格变化与打车人数需求之间构成价格需求，于是可利用价格需求理论模型对此进行分析，但是单一的打车人数多少不足以表示打车是否困难，于是考虑增加一个

空载率指标与打车人数结合表示打车是否困难，最后评判打车困难时，由于打车难这是一个不可量化指标，因此要想得出打车难是否有缓解只能先建立一个标准，然后将价格需求理论模型的求解结果带入该标准。即可知道各公司的补贴整车对打车难是否有帮助。

问题三的分析

由于问题三是设计补贴方案，而问题二里我们建立了价格需求理论模型，求解了各个公司不同补贴方案对打车难缓解的影响，于是我们在问题三过建立一个补贴金额与乘客满意程度的双优化模型来设计一个补贴方案，然后利用问题二的求解结果对设计方案进行论证。

五、模型的建立与求解

5.1 问题一模型的建立与求解

5.1.1 利用ISM模型对影响因素分组

由于出租车资源供求匹配关系受到较多因素影响，其中很多因素又相互包含，必须全面考虑各个因素。所以我们采用ISM模型对相应影响因素进行分组、归类，使问题简化，方便求解。ISM模型是..

J N Warfield于1973年为了分析复杂的社会经济问题而提出的解释结构模型，是静态的定性模型,其特点是能够将复杂的系统逐级分解成若干个子系统。

为了分析出租车资源的供求匹配程度，我们考虑影响出租车与乘客供求关系的一些主要因素如表1。

表1 出租车与乘客供求关系主要影响因素

然后分析这些因素互相之间是否有关系，用0表示相互之间无影响，用1表示相互之间有影响，从而得到相互影响关系的邻接矩阵如下：

由于此矩阵中影响因素较多，所以运用Matlab软件对邻接矩阵进行求解，程序见附录一，得

1000000100101110110100101010110100101111110111110000100000000

0000100000000000010000000000001001000000000101000000000010000000000001000

1

1K ?? ? ? ? ? ? ? ? ?=

? ? ? ? ? ? ? ? ??

?

可达矩阵可说明两个因素之间是否存在路径，并能清楚说明两因素之间影响程度。定义集合()P I 为可达矩阵中要素I 一行中值为1的元素所在行元素集合，

()H I 为可达矩阵中I 这列值为1对应的列元素集合，当()()()H I P I H I =时，I 即为该层元素，然后剔除矩阵中的W 元素，进行下一层元素计算,可以得到最

终的分组：

{}1168,,U A A A =; {}22310,,U A A A =; {}371112,,U A A A =; {}4459,,U A A A =。 通过分析每组所包含的因素，我们发现分组1U 里面所包影响因素可理解为经济水平对出租车供求关系的影响，2U 里面所包影响因素可以看作为不同空间对出租车供求关系的影响，3U 里面所包影响因素可以看作不同时间对出租车供求关系的影响，结果如表2。

时间因素3u

是否是节假日5A 是否是高峰期4A

其它4u

乘客出行的紧急程度11A 旅游吸引能力12A 天气状况7A

5.1.2 问题一模型的建立

我们从出租车空载率角度考虑出租车资源的“供求匹配”程度问题，当出租车空载率过低时，说明打车的人少，出租车量小于乘客的需求；当空载率过高时，表明打车的人较多，出租车量大于乘客需求，出租车空载率能很好地反映出租车与人之间“供求匹配” 程度。所以我们选取空载率这一指标作为模型最终评价因子分析，来分析不同时间和空间出租车资源的“供求匹配”程度问题。

我们利用()AHP —模糊综合评价方法建立模型，首先，利用()AHP 构造出一个层次分析模型，指标评价体系结构图如图1。

图1 层次分析结构图

1）第一、二层权重集的确定

第一层包括4个因素，即1234(,,,)U U U U U ，利用AHP 层次分析法比较几种指标间的关联度如图3。

表 3 第一层因素间关联度

然后确定第一层指标权重，利用1～9标度法求解判断矩阵，构造第一层的评判矩阵Y 具体形式如下：

11121314212223243132333441424344u u u u 11/41/41/4u u u u 4124=u u u u 41/2121/2

1/41/21u u u u Y ????

???

?????=

????

????

??

?? 其中：12u 表示1u 与2u 之间的关联度。

之后求解第一层的权重集W ={1W ,2W 3W ,4W }。方法如下：先计算判断矩阵Y 每行元素的乘积i Q ，再对i Q 求五次方根，得到

{}

{}

1234,,,0.232,0.821,1.736,0.427

W W W W W ==

归一化处理：利用公式 4

1

/

i

i i i W W W ==∑对()1234,,,W W W W W =做归一化处理，

得到1,234(,,)(0.100,0.462,0.305,0.133)W W W W W ==即为所求特征向量。 一致性检验：为了说明以上所求得特征向量是否能够合理的分配权重，需要进行一致性检验，方法如下：通过公式()4

max 14i

i

YW W λ=∑

求得判断矩阵的最大特征值，

得到最大特征值max

5.024λ=。然后利用公式CI CR RI

=和max

1n CI n λ-=-，通过代入

数据 1.12RI =,4n =，max 5.024λ= ，得到0.0350.1CR =<，这就表明评判矩阵Y 具有很好的一致性，所以()1234,,,W W W W W =中的各项均可以作为U 的权重系数。

同第一层权重的求解过程，对于第二层指标，由于第一层的每一个因素都包含()123,,i

i i i u u u u =3个因素，于是得到第二层级的各项权重集：

()()()()

'1'2'3'4=0.2860.0810.5670.0860.1770.764=0.0760.6810.267=0.1820.0510.727W W W W =

2）确定综合评判结果

根据模糊综合评价法可知，综合判别公式12

(,)m B W Y b b b =?=，其中

4

1

()(1,2

)j i ij i b w y j m ===∑，将数据代入计算公式4

'1

()i i i B W W ==∑得到所研究

问题的综合评判结果：()0.1160.3040.580B =。

3）综合评判矩阵B 的修正

影响出租车供求匹配的具体因素有些有具体实际数据，各种数据的单位并不统一，不容易量化，而另外一些因素如：旅游吸引能力、乘客出行的紧急程度等影响因素是模糊的量。为了使各项因素具有可比性，我们将这些有实际数据的影响因素进行标准化处理，消除量纲差异。利用如下公式 '(1,2,;1,2,,)ij j

ij j

i n j m x x x s ==-=

（1）

()

()211,1

111,2,,j j n

n

ij j

ij i i n n j m x x s x x =

-∑∑==-??=

=??

??

（2）

将市民人均输入，人口密度，出租车数量，私家车数量数据分别按照上式标准化，但是得到的数据并不全是处于[]0,1之间，所以应用平移-极差变换使得所有的原始数据都在[]01，之间。得出量化结果：人均收入()11A =，人口密度()50.87A =，出租车数量()60.96A =，私家车数量()81A =。

对于模糊的影响因素，我们用 []01，之间的数来描述，例如：0.2对于天气状况来说意味着“差”，而0.8意味着“好”。其他模糊因素同理，得具体情形如表4。

表4 模糊影响因素处理结果

利用上表对12个因素的定量分析，我们可以对第二层级权重

()123,,(1,2,3,4)i i i i W W W W i ==进行修正，使之与这12个因素联系更加紧密，使结果更准确，第二层权重的修正如下：

''1123(0,2860.0810,567)W A A A = ''2456=(0.0860.0170.764W A A A ）

()''378900760.6810.267W A A A =， ()''41011120.1820.0510.727W A A A =

4)最终评判标准的确定

对于出租车的资源供求匹配是否合理，我们采用出租车空驶率来进行评判，空驶率与出租车供需关系如下表所示：

表5 城市出租空载与交通供求关系

[4]

阶P ，即

0.250.300.40P ??

??=??????

于是我们结合已经求出的评判矩阵B ，可得出最终的评价结果公式:

M B P =?

评价因子M 与空载率是同一量级，我们将M 与空载率标准（表5）进行比较，从而得出出租车资源“供求匹配”的程度。到此，问题一的模型建立完成。 5.1.3 问题一模型的求解

本文以市为例，根据相关数据进行求解哈市出租车资源的“供求匹配”程度问题。我们从统计局获取哈市市区包括下辖区各县相关数据见附件1。为了满足从空间角度考虑问题，首先，我们将哈市分为中心区与边缘区，中心区包括市辖区：南岗区、道里区、道外区、香坊区、平房区、松北区、呼兰区、区。边缘区包括下瞎县城：宾县、延寿县、通河县等。时间角度我们主要从大的方向考虑，比如：季节，是否是节假日，是否是高峰期等。

我们先以哈市冬季南岗区某节假日为例，求解出租车资源“供求匹配”程度问题。南岗区相关数据见附件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.8i A =

评判矩阵为：

''

0.2630.0710.1130.0690.1420.7330.6080.6270.2140.146

0.011

0.582i W ??

?

?= ?

?

??

利用Matlb 求解(程序见附录四)，得到()0.2436 , 0.2639 ,0.4152B =,最后计算出的评价因子0.3062M BP ==。带入表5中进行比较，比较结果：0.300.30620.40<<，得出匹配程度：供需合理。

从而得出结论：哈市冬季南岗区某节假日出租车资源“供求匹配”程度处于供需合理状态。

类似的，利用上述方法可以求解出其他区和所辖县的“供求匹配”程度状况。画出哈市出租车资源“供求匹配”程度图如图2。

对于问题二中“打车难”这个问题，经过分析，发现乘客打车的难易程度最终都可以反映在补贴之后打车人数和出租车空载率这两个指标上。因为城市出租车数量短时间是不会发生大的变化的，而打车的人数却是一个变化的数值，打车

的人多了，打车就变得相对困难。出租车的空载率变小，人们打车就变得相对困难，反之亦然。所以我们用补贴方案实施前后，打车的人数和出租车空驶率的变化两个指标来反映乘客打车的难易程度，由此间接评判出各公司的出租车补贴方案对缓解打车难的现象是否有帮助。

虽然补贴的方案不尽相同，但其本质都是补贴一定数额的钱给相应的出租车司机和乘客，我们将不同补贴方式的差异转化为补贴金额的差异，由于补贴金额是直接影响打车难的关键因素，所以我们假设在这个期间不会有其他因素造成打车人数和空驶率的改变，打车人数和空驶率的改变只受补贴金额的影响。 5.2.2问题二模型的建立

我们借鉴价格理论模型[5]进行建模，价格需求理论是用于研究价格与需求人数的关系的方法，所以在这里我们将其应用于出租车价格变化与打车人数的关系研究上，但是本题中并没有价格变动这一指标，不过，因为补贴方案的不同最终影响的是补贴金额的不同，所以我们这里用补贴金额代替价格变动这一指标。价格需求理论指出价格与需求成负指数关系，所以我们建立如下方程

1220.480733( 5.117)

(1)x P P P r e r --=?+?- （3）

其中：1P 为价格补贴之前打车人数，2P 为价格补贴之后打车人数，r 为受出租车价格影响打车的人数占总人数的比例。在这里受出租车价格影响，打车的人数考虑为使用手机软件打车的人数，而这部分人年龄分布大概在18-35之间，我们通过查找统计局的各年龄人口数据（见附件1），估计出地区26.7%r =。

为了更好地体现补贴前后打车人数变化情况，我们采用补贴前后打车人数比值12

P P 作为评价打车人数变化的指标，根据公式（3）得方程：

12

0.48073(50117)(1)x P r P re ---=+ （4）

利用Matlab 画出函数图象如图4。（程序见附录二）

图4 打车人数变化率与补贴金额的关系

空载率与出租车价格变化的关系式为：

1i

j

E

n l p p i η?

??=-

（5）

E 为出租车日均载客里程；n 为出租车辆数；η为出租车出勤率；l 为平均每辆车的日行里程。

将式（4）带入（5）得到：

空驶率与出租车价格的关系式为： 0.48073( 5.117)

1[(1)]x E r n l

r e

i η--??+-??=-

（6）

我国汽车协会统计结果得出，当空载率大于40%时，供过于求；当空载率小

于25%时，供大于求；当处于25%～40%之间时，处于供求平衡状态。供过于求，车多，可以认为打车容易。供小于求，车少，可以认为打车难。当处于30%～40%之间时，处于供需平衡状态，我们不考虑空载率对打车难的影响，只考虑打车人数变化关系。

最后我们综合考虑补贴前后打车人数比值12

P P 和空载率i 变化情况，建立综合

评价标准如表。

5.2.3 问题二模型的求解

根据参考资料(见附件1)，我们得到各公司补贴方案，考虑到我们采用的是补贴金额数来进行求解的，所以我们对各种补贴方案数据进行处理得到如下的表格。

将各个阶段的补贴金额数分别带进公式（4）（6）中，利用Matlab软件求解（程序见附录三）得：

表8滴滴打车补贴方案评价指标表

表9 快的打车补贴方案评价指标表

根据上述求解结果，对不同补贴方案，打车人数比值12

P P 与空载率i 都不相同，

且当补贴金额大于5元时，

1

2

P P 均小于1且i 均大于40%，而当补贴金额小于5元时

1

2

P P 大于1，且空驶率大幅下降。据此根据表6，我们得出结论：各公司的出租车补贴方案在一定程度上对“缓解打车难”是有帮助的，当补贴方案对应的补贴金额大于5元时，乘客打车变得容易，缓解了“打车难”的状况。当各个公司补贴方案对应的补贴金额小于5元时甚至不补贴时，打车较为困难，对缓解“打车难”没有帮助。

5.3问题三的模型建立及求解 5.3.1线性规划模型的建立

对于补贴方案的合理性，考虑问题二的求解方法，首先我们应用同样的方式，将补贴方案的差异，最终转化为补贴金额的差异。当我们创建一个打车平台并且自己设计补贴方案时，设计的方案要达到的目的是给司机和乘客补贴最少的钱，并使顾客达到最大满意度。于是我们可以将问题三中求解补贴方案合理性的目标转化为求解补贴金额最小以及不满意度最低的优化模型。

对于满意度，我们认为空驶率越高,打车越容易,相应的满意度越高.。在此我们假定乘客满意度就是由空驶率唯一确定, 并且与空驶率成正相关。但是由于我们要求解的优化模型是求解最小值，故而我们把求最大满意度问题转变为求解最低不满意度问题。

我们用0～1之间的数来描述乘客不满意度，最高为1，最低为0。当出租车空载率为30%左右时，打车容易，供需基本处于平衡状态；超过40%时，供过于求，低于25%时，供不应求，打车困难。不满意度与空驶率对应关系如下表：

表10 空驶率与不满意度对应表

通过对上述表中数据进行拟合，得到空驶率与乘客不满意度有如下函数关系：

()1

1.0165328.1908exp(23.973(0.25))

S i i =

++ （7）

综合分析以上数据，并结合式（4）、（6）我们建立优化模型如下：

()1

2

0.48073( 5.117)1221

2

;..

1;

1

1.0165328.1908exp(23.973(0.25))

(1);0.30.4;

0;0;x Min S x

E

s t i n l S i i P r P e r P i P P ds

di P P η--=-?

??=

++=??+-?≤≤>>

利用Lingo 软件对以上模型进行求解（程序见附录五），得到最佳解8x =元。 5.3.2设计方案

根据上面的求解结果，我们可知当补贴金额为8元时，乘客满意度最高，出租车供需较为合理，且打车较为容易。所以我们提供的补贴方案为每单返现8元。考虑到司机接单积极性因素的影响，可以将8元进行分配，由于乘客才是打车这个事件的主要影响因素，所以分配返现金额时，乘客所得金额应该大于出租车司机，具体分配方案如下

表11 出租车补贴方案（单位：元）

5.3.3论证其合理性

将我们得出的补贴方案代入问题二的模型求解，首先将补贴金额8元代入公

式（4）、（6）得出的人数变化率为0.7562，出租车空驶率0.8643，此时打车较为容易，且打车人数受价格变化影响较小。当然，我们代入更高的补贴金额也能达到这个目标，但却无法更进一步影响打车容易度，反而会增加打车软件公司的资产消耗。

最终我们得出结论，4种补贴方案较为合理。

六模型评价及推广

6.1 模型的评价

6.1.1 模型优点：

1）问题一中的模型合理地将ISM解释结构模型、层次分析法、模糊综合评价模型相结合，较好的解决了相应的问题，与单一模型相比较更具灵活性。

2）问题二中将补贴方案的差异具体化成补贴金额的差异，将打车难易的模糊说法具体为补贴前后打车人口数不同，这样做既可以说明问题的主要因素，又巧妙的简化了问题的分析与求解过程。

6.1.2模型的缺点：

1 ）问题一中各类元素相互之间的相关系数在确定过程中具有一定的主观性。

2）在确定影响出租车“供求匹配”程度因素的过程中，由于时间的原因不能将因素集考虑的更加全面。

6.2 模型的推广

问题一中建立的数学模型用于评判不同时空出租车资源的供求匹配程度，联系实际情况，有很多与之类似的情形。比如不同地点的饭店在不同时间段提供的席位与客人的供求匹配程度的相关问题，还有不同地点的城市图书馆提供的图书与阅览者的供求匹配程度等。这样类似的实际问题有很多，只要确定好衡量指标，再结合具体问题，都可以用问题一中的模型来解决，这也就说明了模型的普适性和推广的价值。

七、参考文献

[1]凤秋，善强，作宝，数学实验，：工业大学，2010.

[2]守建．基于ISM模型的标准信息化影响因素分析[J].工业大学学报，

2010,42(8):1 306-1 310.

[3]丁树良，淑群，汪文义.可达矩阵在认知诊断测验编制中的重要作用［J］．师

大学学报: 自然科学版，2010，34( 5) : 490-495．

[4]道勇，戴建军.基于 FCD 的市出租车汽车空驶率特点研究[A].第三届

中国智能交通年会学术委员会. 第三届中国智能交通年会论文集[C].:东南大学,2007.

[5]梁小民. 西方经济学导论. : 大学, 1993.

八、附录

附录一：

%根据邻接矩阵求解可达矩阵计算代码

A=xlsread('C:\Documents and Settings\Administrator\桌面\kedajuzhen.xls');

I=eye(12);

R=A+I;

K=0;

while 1

Rnew=R * (A+I)>0;

if isequal(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*10000*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.117))+1-0.267)/(336*0.9*1613 6)

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.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.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;

i>0.3;

i<0.4;

diff(s)>0;

H1/H2>0;

End