污水排放问题

污水排放问题

摘要：我国是一个干旱缺水严重的国家，淡水资源总量为28000亿立方米，占

全球水资源的6%，仅次于巴西、俄罗斯和加拿大，居世界第四位，但人均只有2300立方米，仅为世界平均水平的四分之一，美国的五分之一，在世界上排名121位，是全球13个人均水资源最贫乏的国家之一。而污水处理是我国走可持续开展路线亟待解决的问题。

本文研究了污水排放问题，从处理污水总开支最少考虑，首先对管道长与管道费用，排污量与建站费用数据处理，进展拟合得出相关方程，并在可承受的误差范围内，用得到的方程计算相关费用。根据实际问题，建立了五种方案。通过对五种方案分别计算并进展比拟得出最优方案。费用分担问题上综合各方面因素拟定了五种分担方案，多定方案有利于各工厂依据实际进展有针对的选择。

本模型利用了MATLAB软件模拟曲线、Lingo软件编程求解数据关系和最终解，费用分担使用了对实际的分析和Shapley合作对策公平分配的方法得出费用分担的五种方案，此模型能够解决污水的排放问题和费用的分担问题，具有较强的规律性，可以推广到业、商业、交通运输、工程技术、行政管理等领域。

在文章的最后，分析了在模型中还需进一步考虑的问题。 关键词：污水排放 开支最少 费用分担

- -可修编-

. -

1 问题重述

1.1问题的背景与简化

随着国民经济的快速开展和构造转型，企业在追求经济效益的同时，越来越重视环境保护问题。如何减少污染物的排放以保护环境，使经济得以稳健及可持续开展，是许多企业亟待解决的重要问题。

假设沿河有假设干工厂，每天都会排放一定量的污水，这些污水必须经过处理才能排入河中。通常的解决方法是建造污水处理站，将污水进展处理，使之到达排放标准后再予以排放。

污水处理站可以由每个工厂单独建造，也可以几个工厂联合建造。联合建造时，处理站必须建在下游位置，上游工厂将污水通过管道送往下游的处理站集中处理。处理站的建造费用与污水处理量及铺设的管道总长度有关，题目所给表中显示了不同污水处理量和不同管道铺设总长度的建造费用及管道铺设费用。

1.2问题提出

问题1：建立适当的数学模型，给出合理的污水处理站建造方案，从不同的方案中得出费用最低的方案。如果是联合建造，应给出建造费用的分担方法。

问题2：假设沿河从上游到下游有A，B，C三家工厂，各厂的排污量分别为4.5 t/s，2.5 t/s和6 t/s。AB之间的距离为20 km，BC之间的距离为40 km。用已经建立的模型给出具体的污水处理站建造方案和费用分担方法。

问题3：分析说明你所给方案的合理性。

2 问题分析

2.1问题1的分析

在对污水排放的问题上，由题可知污水处理站的建造费用与污水处理量及铺设的管道总长度有关且题中表格给出建站费和管道费两种费用。由实际情况可知建站费只可能与排污量有关，故求建站费的问题时只需要考虑排污量；管道费应该与排污量和管道长度都有关，故求解管道费问题时需要两个因素综合考虑。两项费用与管道的长度、排污量之间的关系可以由MATLAB中的cftool进展曲线拟合得出，将曲线拟合得到的关系式带入预设的方案那么可得出最优方案。根据实际和已有的经历拟定了五种费用分担方案，方案的最终确立一定要由工厂共同商量得出。

2.2问题2的分析

沿河从上游到下游有A，B，C三家工厂，三家工厂建立污水处理厂可用有五种不同的组合方案，通过问题1得出的方程式将五种方案的总费用lingo软件求解出来然后进展比拟从而得出最正确的方案。确定方案后拟定费用的分担方法，将公式输入lingo软件中最后求得各厂所需支付的总费用。

3 模型假设

1) 题目所给出的表格中可以反映全国污水处理站的费用与管道长度、排污量之

间的关系；

2) 实际中河边的工厂个数为有限个，可以使用列举法得出多种建厂方案； 3) 联合建站费按排污量或函数值之比进展分担；

4) 管道费按谁使用谁投资的原那么或直接按排污量、函数值的比例进展分担； 5) 由已有经历使用Shapley合作对策公平分配的方法进展合理分担费用。

- -可修编-

. -

4 符号说明

表1 模型符号

X Y 排污量单位〔立方米/秒〕 管道长度单位〔公里〕 管道费单位〔万元〕 建站费单位〔万元〕 总费用单位(万元) 工厂单独建站的费用单位(万元) 特征函数 单位(万元) Shapley值 单位(万元) W Z S P V 

5 模型的建立与求解

5.1模型的建立

5.1.1费用关系式的建立

通过题目所给的以下表格可以做出散点图：

表2 不同污水处理量和不同管道铺设长度的建造费用及管道铺设费用 排污量(t/s) 2 4 6.5 8 9 9.5 10 11 12.5 15 管道总长(km) 8 15 21 26 30 38 45 51 56 65 建站费(万元) 125.18 195.8 265.21 332.56 322.93 372.63 363.29 374.56 430.18 459.29 管道费(万元) 7.08 19.91 36.41 49.58 62.06 79.55 97.06 116.66 141.12 181.18 总费用(万元) 132.26 215.71 301.62 382.14 384.99 452.18 460.35 491.22 571.3 640.47 - -可修编-

. -

图 1 建站费与排污量关系散点图

图 2 管道费与管道长度关系散点图

- -可修编-

. -

图 3 管道费与排污量关系散点图

观察数据所得的散点图可以猜想方程模型为：

ZaXb〔1〕 WaXbY〔2〕

通过MATLAB中cftool可以模拟出一下曲线：

图 4 建站费与排污量的曲线拟合

- -可修编-

. -

图 5 管道费与排污量、管道长度的曲线拟合

从拟合曲线中可以得出方程：

Z=82.06X0.6451 〔3〕 W0.5452X0.6022Y 〔4〕

5.1.2费用分担的模型建立

因为实际情况中河边工厂的数量为有限个，所以建造污水处理站的方案可以用列举法一一列出，如果某工厂独立建造污水处理站那么费用自然由自己承当，如果几个工厂联合建造污水处理站那么由以下方案分担费用： 〔1〕联合建站费按排污量或函数值之比进展分担； 〔2〕管道费按谁使用谁投资的原那么或直接按排污量、函数值的比例进展分担； 〔3〕由已有经历使用Shapley合作对策公平分配的方法进展合理分担费用。

通过以上描述总共可以拟定五种费用分担方案，最终方案应由所有公司共同商议给出。

5.2模型的求解

5.2.1下面通过以上模型对问题1进展求解：

用Ni〔i=1,2,3…n〕表示n个工厂设每个工厂距离下一个工厂的距离分别为Yi〔i=1,2,3…n〕，每个工厂的排污量分别为Xi〔i=1,2,3…n〕对应的管道费用为Wi〔i=1,2,3…n〕对应的建厂费为Zi〔i=1,2,3…n〕，每个工厂单独建立污水处理厂的费用分别为Pi〔i=1,2,3…n〕，第i个和第j个工厂联合建立污水处理厂的费用为P。 ij〔i，j=1,2,3…n〕

〔1〕对n个厂的问题进展讨论可组合出以下几种方案： 方案1：各工厂分别建站,总费用：

S1P1P2...Pn 〔5〕

方案2：利用组合的方式对n个中某些工厂进展联合建立污水处理厂：

假设n个工厂中有k〔k=2,3,…,n〕个工厂采用联合建厂的方式，剩余(n-k)个单独建厂。

- -可修编-

. -

a．将n-k个单独建厂的工厂按从上游到下游的位置重新记为M(i)到M(j)，那么其建厂费用为:

PM(i)82.06Xi0.6451〔6〕 b．K个采用联合建厂的工厂重新按从上游到下游位置记为Q(i)〔i=1,2,…,k〕,用Yij表示第i个工厂与第j个工厂之间的距离，也就是管道长度。那么管道费用：

WQ(i)0.5452Xki0.6022YQ(i)Q(k) 〔7〕

总管道费为：

WWQ(i) 〔8〕

k2i1nk建站费用：

ZQi总联合建厂费为：

XQ(i)Xi1k82.06(XQ(i))0.6451 〔9〕

i1kQ(i)ZZQ(i) 〔10〕

k1i1nk总费用：

S2PM(i)ZW〔11〕

i1nk(n1)!2种建厂方式。

(i1)!(ni)!i2方案3：所有工厂一起联合建立污水处理厂：

S3P1...n 〔12〕

通过组合共可得到n〔2〕费用分担方法由问题分析所给原那么可以列出五种方案：

方案1：管道费用根据谁使用谁投资的原那么，及由处于管道上游的工厂按排污量比例分担费用；建站费由每个工厂排污量的比例进展费用分担。 单独建厂费用：SiPi

联合建厂第i个工厂所需费用：

Xki表示k个联合建厂中的第i个厂的排污量，Zki表示k个联合建厂中的第i个厂的建站费，Wik表示k个联合建厂管道总费用，Wi表示k个联合建厂中的第i个厂的管道费用。

Xki 〔13〕 Zki=ZXk1Xk2...XkkWik0.5452(Xki)0.6022Yk(i1),kk 〔14〕

i1kWiWikXkiXi1k 〔15〕

kiSiZiWik 〔16〕

方案2：管道费用根据谁使用谁投资的原那么，及由处于管道上游的工厂按管道费方程W0.5452X0.6022Y中X0.6022的比例分担；建站费按建站费方程

- -可修编-

. -

Z=82.06X0.6451中X0.6451的比例分担费用。

a．单独建厂费用：

SiPi 〔17〕

b．联合建厂第i个工厂所需费用：

Xki表示k个联合建厂中的第i个厂的排污量，Zki表示k个联合建厂中的第i个厂的建站费，Wki表示k个联合建厂中的第i个厂的管道费用。

0.6451Xki 〔18〕 Zki=Z0.64510.64510.6451Xk1Xk2...XkikkXWki[0.5452(Xki,k(i1))0.6022Yki,k(i1)iki] 〔19〕

i1i1Xkjj1SiZiWik 〔20〕

方案3：联合建的建污水处理站总费用都按各个工厂的排污量之比分担。

XSiSni 〔21〕

Xii1方案4：联合建的污水处理站的建站费按各个工厂的X0.6451比例分担，管道费按各个工厂的X0.6022比例分担。 a．单独建厂费用：

SiPi 〔22〕

b．联合建厂第i个工厂所需费用：

Xki表示k个联合建厂中的第i个厂的排污量，Zki表示k个联合建厂中的第i个厂的建站费，Wik表示k个联合建厂管道总费用，Wi表示k个联合建厂中的第i个厂的管道费用。

0.6451Xki 〔23〕 Zki=Z0.64510.6451Xk1Xk0.6451...X2kiWik0.5452Xki0.6022Yki,k(i1)〔24〕

i1kWiWik0.6022XkiXi1k 〔25〕

0.6022kiSiZiWi 〔26〕

方案5：按Shapley方案。

定义1 设集合I={1,2,⋯,n}, 如果对于I 的任一子集s 都对应着一个实值函数V(s), 满足方程：V()0。

V(s1s2)V(s1)V(s2),s1s2称[I,V]为n 人合作对策, v 为对策的特征函数。

定义2 用Ni表示I的成员i从合作的最大效益V(I)中应得到的一份入。

N(N1,N2,,Nn)，叫作对策的分配( Imputation) ,满足NiV(I)，

i1nNiV(i),i1,2,n。

- -可修编-

. -

公式Shapley 值(V)(1(V),2(V),,n(V))为：

i(V)R(|s|)[V(s)V(s\\i)],i1,2,n 〔27〕

SSi(n|s|)!(|s|1)!其中Si是I 中包含i 的所有子集, |s|是子集s中的元

n!素数目( 人数) , R(|s|)是加权因子, s\\i表示s去掉i后的集合。 第i个工厂所需承当的费用为：

SiPii(V) 〔28〕 R(|s|)

5.2.2下面通过以上模型对问题2进展求解： 由

Z82.06X0.6451〔29〕

得

PA82.06X0.6451〔30〕

记Pi为一工厂单独建站的费用〔i=A,B,C〕那么： 可以算得：

PA=216.5304（万元） PB=148.1981（万元） PC=260.6855（万元）记Pij为i,j工厂合作,两家工厂需要的总费用那么：

〔31〕 PAB82.06(X1X2)0.64510.5452X10.6022Y1314.9158（万元）〔32〕 PAC82.06(X1X3)0.64510.5452X10.6022(Y1Y2)454.9480（万元）〔33〕 PBC82.06(X2X3)0.64510.5452(X2X3)0.6022Y2396.7558（万元）对三个厂的问题进展讨论可组合出以下5种方案： 方案1：各工厂分别建站,总费用：

〔34〕 S1PAPBPC625.4140 （万元）方案2：工厂A,B联合建站,工厂C单独建站,总费用：

〔35〕 S2PABPC575.6013（万元）方案3：工厂A,C联合建站,工厂B单独建站,,总费用：

〔36〕 S3PACPB603.1461（万元）方案4:工厂B,C联合建站,工厂A单独建站,总费用：

〔37〕 S4PBCPA631.2862（万元）方案5：三厂联合在城3建站,总费用为：

S582.06(XAXBXC)0.64510.5452XA0.6022YAB0.5452(XAXB)0.6022YBC526.6438（万元）〔38〕

通过结果的比照可见方案5最优。

由以上所得结论可以制定以下费用分担方法：

方案1：管道费用根据谁使用谁投资的原那么，及由处于管道上游的工厂按排污量比例分担费用；建站费由每个工厂排污量的比例进展费用分担。 〔1〕建站的费用的分担方法： 三厂联合建站费用：

- -可修编-

. -

Z82.06(XAXBXC)=429.2758（万元） 〔39〕

按三城污水量比例9:5:12分担，厂A,B,C分别负担费用：

XA 〔40〕 ZA=Z=148.5955（万元）XAXBXCXB 〔41〕 ZB=Z=82.55303（万元）XAXBXCXC 〔42〕 ZC=Z198.1273（万元）XAXBXC〔2〕管道费用的分担方法：

工厂A到工厂B的管道费用由工厂A承当： 工厂A承当的费用：

WAB0.5452XA0.60220.6451YAB26.97431（万元） 〔43〕

工厂B到工厂C的管道费用由工厂A与工厂B按排污量共同分担：

因为

WBC0.5452(XAXB)0.6022YBC70.39371（万元）〔44〕

所以A工厂应承当费用：

WBCAWBCB工厂应承当费用：

XA 〔45〕 45.25310（万元）XAXBXB 〔46〕 25.14061（万元）XAXB根据费用分担的原那么可以得出每个工厂需要分担的总费用：

〔47〕 SAZAWABWBCA220.8229（万元）WBCBWBC 〔48〕 SBZAWBCB107.6936（万元） 〔49〕 SCZC198.1273（万元）

方案2：管道费用根据谁使用谁投资的原那么，及由处于管道上游的工厂按管道费方程W0.5452X0.6022Y中X0.6022的比例分担；建站费按建站费方程Z=82.06X0.6451中X0.6451的比例分担费用。〔例如：A，B两个工厂联合建污水处

0.60220.6451XAXA理站时费用A分担的方程为SAWAB0.6022〕 ZAB0.64510.60220.6451XAXBXAXB〔1〕建站的费用的分担方法： 三厂联合建站费用：

Z82.06(X1X2X3)0.6451=429.2758（万元） 〔50〕

按三个工厂X0.6451比例4.50.6451:2.50.6451:60.6451分担，厂A,B,C分别负担费用：

0.6451XA 〔51〕 ZA=Z0.6451=148.6236（万元）0.64510.6451XAXBXC- -可修编-

. -

0.6451XB 〔52〕 ZB=Z0.6451=101.7212（万元）0.64510.6451XAXBXC0.6451XC 〔53〕 ZC=Z0.6451=178.9310（万元）0.64510.6451XAXBXC〔2〕管道费用的分担方法：

工厂A到工厂B的管道费用由工厂A承当： 工厂A承当的费用：

WAB0.5452XA0.6022YAB26.97431（万元） 〔54〕

工厂B到工厂C的管道费用由工厂A与工厂B按排污量共同分担：

因为

WBC0.5452(XAXB)0.6022YBC70.39371（万元） 〔55〕

所以A工厂应承当费用：

WBCA0.6022XA 〔57〕 WBC0.602241.36185（万元）0.6022XAXB0.6022XB 〔58〕 WBCBWBC0.602229.03186（万元）0.6022XAXB根据费用分担的原那么可以得出每个工厂需要分担的总费用：

〔59〕 SAZAWABWBCA216.9597（万元） 〔60〕 SBZBWBCB130.7530（万元） 〔61〕 SCZC178.9310（万元）方案3：联合建的建污水处理站总费用都按各个工厂的排污量之比分担。 三厂联合建站费用：

Z82.06(XAXBXC)=429.2758（万元） 〔62〕

按三城污水量比例9:5:12分担，厂A,B,C分别负担费用：

XA 〔63〕 ZA=Z=148.5955（万元）XAXBXCXB 〔64〕 ZB=Z=82.55303（万元）XAXBXCXC 〔65〕 ZC=Z198.1273（万元）XAXBXC所需的总的管道费用：

0.6022WWABWBC0.5452XAYAB0.5452(XAXB)0.6022YAB97.36803（万元）0.6451〔66〕

按三城污水量比例9:5:12分担，厂A,B,C分别负担费用：

XAWA=W=33.70432（万元） 〔67〕

XAXBXCXBWB=W=18.72462（万元） 〔68〕

XAXBXC- -可修编-

. -

XC 〔69〕 =44.93909（万元）XAXBXC根据费用分担的原那么可以得出每个工厂需要分担的总费用：

〔70〕 SAZAWA182.2998（万元）WC=W 〔71〕 SBZBWB101.2777（万元） 〔72〕 SCZCWC243.0664（万元）方案4：联合建的建污水处理站的建站费按各个工厂的X0.6451比例分担，管道费

按各个工厂的X0.6022比例分担。 〔1〕建站的费用的分担方法： 三厂联合建站费用：

Z82.06(XAXBXC)=429.2758（万元） 〔73〕

按三个工厂X0.6451比例4.50.6451:2.50.6451:60.6451分担，厂A,B,C分别负担费用：

0.6451XA 〔74〕 ZA=Z0.6451=148.6236（万元）0.64510.6451XAXBXC0.6451XB 〔75〕 ZB=Z0.6451=101.7212（万元）0.64510.6451XAXBXC0.6451XC 〔76〕 ZC=Z0.6451=178.9310（万元）0.64510.6451XAXBXC〔2〕管道费用的分担方法： 所需的总的管道费用：

0.6022WWABWBC0.5452XAYAB0.5452(XAXB)0.6022YAB97.36803（万元）0.6451按三个工厂X0.6022比例4.50.6022:2.50.6022:60.6022分担，工厂A,B,C分别负担费用：

0.6022XA 〔77〕 WA=W0.6022=33.67908（万元）0.60220.6022XAXBXC0.6022XB 〔78〕 WB=W0.6022=23.63933（万元）0.60220.6022XAXBXC0.6022XC 〔79〕 WC=W0.6022=40.04962（万元）0.60220.6022XAXBXC根据费用分担的原那么可以得出每个工厂需要分担的总费用：

〔80〕 SAZAWA182.3026（万元） 〔81〕 SBZBWB125.3605（万元） 〔82〕 SCZCWC218.9806（万元）方案5：按Shapley方案。

定义联合建厂比单独建厂节约的投资为特征函数，于是Vi0，i=1,2,3；

〔83〕 VAB625.4140575.601349.8127（万元） 〔84〕 VAC625.4140603.146122.2679（万元） 〔85〕 VBC625.4140631.2862-5.8722（万元） 〔86〕 VABC625.4140526.643898.7702（万元）由MATLAB求解可得：

- -可修编-

. -

A46.8492（万元） B32.8242（万元） C19.0518（万元） SA216.530446.8492169.6362（万元） SB148.198132.8242115.3739（万元） SC260.685519.0518241.6337（万元）

6 模型的检验

检验方法是将所有费用分担方案得出的各厂所需费用与每个厂单独建立污

水处理厂的费用进展比拟。 〔1〕费用分担方案1的检验：

工厂A分担费用为SA220.8229（万元）>PA 工厂B分担费用为SB107.6936（万元）工厂C分担费用为SC198.1273（万元）工厂A分担费用为SA216.9597（万元）>PA

工厂B分担费用为SB130.7530（万元）工厂B分担费用为SB101.2777（万元）工厂C分担费用为SC218.9806（万元）-

-可修编-

. -

工厂A分担费用为SA169.6362（万元）工厂B分担费用为SB115.3739（万元）对于费用分担问题，拟定的五种方案通过以上检验可以看出方案1、方案2对工厂A太不公平所以工厂A不可能同意这两种方案。方案3、方案4、方案5没有明显的不公平性，最终方案可以按照各工厂的实际情况共同商议来确定。

7模型的评价与推广

〔1〕本模型利用了Lingo软件、MATLAB软件编程共同求解的方法；

〔2〕本模型着重考虑了实际问题，通过实际可能进展了有针对性的求解，具有一定的实际指导意义；

〔3〕本模型的假设具有一定合理性，为节省开支和费用的分担提供了便利； 〔4〕模型具有一般性，在一定条件下能使用于其他工业生产组织、经济方案、输气管道建立、联合配电、污水处理等一系列问题；

〔5〕模型也带有主观性，联合建立污水处理厂时费用分担问题给定的五种方案都具有一定的局限性，实际使用时可能会有一些问题，如Shapley方法的问题是它需要知道所有合作的获利，这在实际上常常做不到，因此还可以考虑使用，如协商解.均衡解。满意解。Raiffa解等等。

参考文献

[1]叶其孝.大学生数学建模竞赛辅导教材〔一〕.XXXX教育，1993。

[2]姜启源.数学建模.2版.：高等教育，1993。

[3]任福田等泽，道路通行能力手册，工业大学交通工程室：建筑工业, 1985。 [4]李进平,X乐文，交通流的流体动力学模型研究[J],XX理工大学学报(交通科学与工程版)，Vo1.26.No.6.,Dec．2002：748—751。 [5]甘应爱等，运筹学，：清华大学，2005。

[6]茆诗松，贺思辉，概率论与统计学，XX：XX大学，2021。 [7]叶其孝，姜启源译，数学建模，：机械工业，2005。

- -可修编-

. -

附录

散点图的绘制： clc,clear

load data2.txt

paiwu=data2([1],:);paiwu=paiwu';paiwu=paiwu(:);

changfei=data2([2],:);changfei=changfei';changfei=changfei(:); plot(paiwu(1:10),changfei(1:10),'*') 建厂费与排污量的散点图

clc,clear

load data1.txt

guanchang=data1([1],:);guanchang=guanchang';guanchang=guanchang(:); guanfei=data1([2],:);guanfei=guanfei';guanfei=guanfei(:); plot(guanchang(1:10),guanfei(1:10),'*') 管道费与管道长度的散点图

clc,clear

load data3.txt

paiwuliang=data3([1],:);paiwuliang=paiwuliang';paiwuliang=paiwuliang(:);

guanfei=data3([2],:);guanfei=guanfei';guanfei=guanfei(:); plot(paiwuliang(1:10),guanfei(1:10),'*') 管道费与排污量的散点图

总费用方案一求解： X1=4.5; X2=2.5; X3=6; Y1=20; Y2=40;

P1=82.06*X1^0.6451; P2=82.06*X2^0.6451; P3=82.06*X3^0.6451;

P12=82.06*(X1+X2)^0.6451+0.5452*X1^0.6022*Y1;

P13=82.06*(X1+X3)^0.6451+0.5452*X1^0.6022*(Y1+Y2); P23=82.06*(X2+X3)^0.6451+0.5452*(X1+X2)^0.6022*Y2; S1=P1+P2+P3; S2=P12+P3; S3=P13+P2; S4=P23+P1;

S5=82.06*(X1+X2+X3)^0.6451+0.5452*X1^0.6022*Y1+0.5452*(X1+X2)^0.6022*Y2;

- -可修编-

. -

Z=82.06*(X1+X2+X3)^0.6451; Z1=Z*9/26; Z2=Z*5/26; Z3=Z*12/26;

G12=0.5452*X1^0.6022*Y1;

G23=0.5452*(X1+X2)^0.6022*Y2; G231=G23*9/14; G232=G23*5/14; SA=Z1+G12+G231; SB=Z2+G232; SC=Z3;

总费用方案二求解： XA=4.5; XB=2.5; XC=6; YAB=20; YBC=40;

P1=82.06*XA^0.6451; P2=82.06*XB^0.6451; P3=82.06*XC^0.6451;

P12=82.06*(XA+XB)^0.6451+0.5452*XB^0.6022*YAB;

P13=82.06*(XA+XC)^0.6451+0.5452*XA^0.6022*(YAB+YBC); P23=82.06*(XB+XC)^0.6451+0.5452*(XA+XB)^0.6022*YBC; S1=P1+P2+P3; S2=P12+P3; S3=P13+P2; S4=P23+P1;

S5=82.06*(XA+XB+XC)^0.6451+0.5452*XA^0.6022*YAB+0.5452*(XA+XB)^0.6022*YBC;

WAB=0.5452*XA^0.6022*YAB;

WBC=0.5452*(XA+XB)^0.6022*YBC;

WBCA=WBC*XA^0.6022/(XA^0.6022+XB^0.6022); WBCB=WBC*XB^0.6022/(XA^0.6022+XB^0.6022); Z=82.06*(XA+XB+XC)^0.6451;

ZA=Z*XA^0.6451/(XA^0.6451+XB^0.6451+XC^0.6451); ZB=Z*XB^0.6451/(XA^0.6451+XB^0.6451+XC^0.6451); ZC=Z*XC^0.6451/(XA^0.6451+XB^0.6451+XC^0.6451); SA=ZA+WAB+WBCA; SB=ZB+WBCB; SC=ZC;

总费用方案三求解：

- -可修编-

. -

XA=4.5; XB=2.5; XC=6; YAB=20; YBC=40;

P1=82.06*XA^0.6451; P2=82.06*XB^0.6451; P3=82.06*XC^0.6451;

P12=82.06*(XA+XB)^0.6451+0.5452*XA^0.6022*YAB;

P13=82.06*(XA+XB)^0.6451+0.5452*XA^0.6022*(YAB+YBC); P23=82.06*(XB+XC)^0.6451+0.5452*(XA+XB)^0.6022*Y2; S1=P1+P2+P3; S2=P12+P3; S3=P13+P2; S4=P23+P1;

S5=82.06*(XA+XB+XC)^0.6451+0.5452*XA^0.6022*YAB+0.5452*(XA+XB)^0.6022*YBC;

SA=S5*9/26; SB=S5*5/26; SC=S5*12/26;

总费用方案四求解： XA=4.5; XB=2.5; XC=6; YAB=20; YBC=40;

P1=82.06*XA^0.6451; P2=82.06*XB^0.6451; P3=82.06*XC^0.6451;

P12=82.06*(XA+XB)^0.6451+0.5452*XB^0.6022*YAB;

P13=82.06*(XA+XC)^0.6451+0.5452*XA^0.6022*(YAB+YBC); P23=82.06*(XB+XC)^0.6451+0.5452*(XA+XB)^0.6022*YBC; S1=P1+P2+P3; S2=P12+P3; S3=P13+P2; S4=P23+P1;

S5=82.06*(XA+XB+XC)^0.6451+0.5452*XA^0.6022*YAB+0.5452*(XA+XB)^0.6022*YBC;

SA=S5*(9/26)^0.6451/((9/26)^0.6451+(5/26)^0.6451+(12/26)^0.6451); SB=S5*(5/26)^0.6451/((9/26)^0.6451+(5/26)^0.6451+(12/26)^0.6451); SC=S5*(12/26)^0.6451/((9/26)^0.6451+(5/26)^0.6451+(12/26)^0.6451);

- -可修编-

. -

总费用方案五求解： function v=eigen(x) s=sum(x);

if s==0 v=0; return; end if s==1 v=0; return; end

if x==[1,1,0] v=49.8127; return; end if x==[1,0,1] v=22.2679; return; end if x==[0,1,1] v=-5.8722; return; end if s==3 v=98.7702; end

function f=shapley(v,n)

x=zeros(1,n); f=x; e=ones(1,n); s=0; t=s; while 1 k=1;

while x(k)==1 p=feval(v,x); if s==1 w=1; else

if s>t w=w*t/(n- t); elseif sf=f+w*x*p;

if sx(k)=0; s=t- 1;

if k==n f=f/n; return; end k=k+1; end

x(k)=1; s=t+1; end

- -可修编-