考虑属性关联的两阶段获胜者确定模型

杨 娜¹，廖貅武²，雷宏振¹

(1.陕西师范大学国际商学院，陕西 西安 710119；2.西安交通大学管理学院，过程控制与效率工程教育部重点实验室，陕西 西安 710049)

摘 要： 研究了多属性逆向拍卖的获胜者确定问题。考虑属性之间的两两关联，以定义在2-可加模糊测度上的Choquet积分表达拍卖人的偏好，其中模糊测度的值由拍卖人提供的偏好信息推测得出。由于一般情况下与偏好信息一致的模糊测度取值并不唯一，考虑所有一致的评分函数提出了两阶段获胜者确定方法。首先采用线性规划挑选出在任意评分函数下可能获胜的报价，再通过混合整数规划确定一个与所有一致的评分函数的评价结果最为接近的报价排序，以得分最高者为稳健获胜报价。仿真实验表明，大量的报价为不可能获胜报价，说明了在第一阶段进行筛选的必要性。与现有方法的比较表明了该方法的有效性，且在拍卖轮数较大、报价数目较多时，该方法在计算效率上更有优势。

关键词： 多属性逆向拍卖；获胜者确定问题；Choquet积分；偏好揭示

1 引言

逆向拍卖是一种由买方发起并主导，多个卖方参与竞价的拍卖机制。实践表明，企业采用这种方式采购产品和服务不仅可以提升采购效率，还能大幅度降低采购成本。近年来，逆向拍卖日益成为一种重要的采购解决方案。它与正向拍卖最主要的区别在于其多属性特征。例如，在运输服务采购中，托运商不仅根据运输价格来选择承运人，还需考虑及时性、设备的合适性以及业务的熟悉程度等其他因素。这种多属性特征引出了一个关键的研究问题——获胜者确定问题(winner determination problem, WDP)，即如何综合多个属性评价投标人的报价。

目前主要的研究思路是采用评分函数进行打分。Bichler和Kalagnanam^[1]以标准的加性价值函数为评分规则，研究了多单位采购情形下的获胜者确定问题，建立了混合整数规划模型筛选获胜者。Kameshwaran等^[2]在此基础上，将存在交易约束以及投标人成本函数具有批量折扣特点的单一产品大规模采购下的报价评价问题归结为一个混合整数优化问题，提出采用目标规划技术进行求解。Cheng Chibin^[3]将密封多属性逆向拍卖机制下买方的获胜者确定问题归结为一个多属性决策问题，并使用TOPSIS方法对投标人的报价进行评价。Ray等^[4,5]使用了拟线性评分函数，研究了中标者实际交付与承诺报价不一致情形下的获胜者确定问题。王明喜等^[6]基于加权和评分函数，提出了一种赢者决策方案，并给出了风险中性者和风险规避者的均衡投标策略。在这些以评分函数为基础的研究中，Pham等^[7]指出加性效用函数是目前最为普遍使用的评分函数。该函数必须以属性相互偏好独立为前提，这在现实中很难满足。

针对这一不足，已有一些学者进行了改进。Talluri等^[8]基于数据包络分析思想提出了一个考虑属性之间相互关联的价值函数。但该函数只能反映部分属性之间的交互关系，具有一定的局限性。呼大永等^[9]采用数据包络分析中的C²R模型确定最终获胜供应商，在解决传统指标体系评价方法中假设属性不相关以及人为设定权重等不足方面做出了有益尝试。Huang Min等^[10]考虑到拍卖属性之间的冲突性，引入了BOCR框架，并在其中嵌入前景理论刻画拍卖人的风险规避行为。所提出的PT-BOCR方法首先采用Choquet积分得到下层的评价值，再使用CCSD方法求解上层模型的权重，最后综合二者得到总评价值。刘卫锋等^[11]还提出了一种Heronian平均算子来刻画属性间的相互作用。在多准则决策领域，除了这些模型和方法，可反映属性间交互关系的偏好函数还有两类：一类是带惩罚和奖励的一般加性效用函数^{[12, 13]}；另一类是模糊积分，其中建立在k-可加模糊测度上的Choquet积分最为常用。Choquet被广泛应用于管理领域的现实问题^[14-16]。这一类型的偏好函数在应用上的最大困难是参数值的获得。虽然文献中提出了一些方法，如Marichal熵方法^[17]、菱形成对比较法^[18]以及序关系分析法与施密特正交马田系统相结合的方法^[19]，但这些方法过于强调理论推演，计算过程复杂，不适合应用于多属性逆向拍卖的现实问题。

文献中有一类研究涉及到拍卖人的偏好揭示。它们大都强调了揭示拍卖人偏好并构造恰当的评分函数的重要性^{[1, 20]}。最新的多属性在线逆向拍卖的综述性文章^[7]还专门将偏好揭示列为一项重要的研究内容。于红岩和刘仲英^[21]通过设定多属性的期望值和保留值作为投标参考点，提出了一种基于拍卖方偏好揭示的多属性网上拍卖模型。Karakaya和Köksalan^[22]在多轮次的多属性拍卖框架下提出了一种交互式决策方法，从拍卖人的历史偏好信息中间接地估计其偏好函数。拍卖过程中不断将新的偏好信息加入偏好约束集，从而进一步精确对偏好函数的估计。Yang Na等^[23]进一步扩展了这一研究，设计了交互式偏好引出框架，并提出了相应的偏好揭示方法。这些方法都为获得拍卖人的偏好函数进而确定获胜者提出了新的思路。

本文将从拍卖人偏好揭示的角度构造Choquet积分形式的评分函数。在获胜者确定方法的选择上，文献中有一些学者研究了两阶段拍卖机制^[24]，这种方式有利于降低拍卖的复杂性并易于被拍卖人所接受，受到这些研究的启发，本文采用一种两阶段的方法确定获胜者。这一方法首先考虑所有与拍卖人主观偏好一致的参数值筛选可能获胜报价，再以与整体评价结果最为接近的报价排序为依据确定稳健获胜报价。该方法的理论贡献主要反映在：以模糊测度和模糊积分表达拍卖人评价属性之间的关联性，避免了属性相互偏好独立的强假设，扩展了现有研究；以所有与主观偏好相一致的参数集表示拍卖人的评分函数，避免人为选择的随意性；侧重分析拍卖人更关注的潜在获胜报价，避免分析大量不可能获胜报价的偏好关系，提升了评价效率；提出了选择稳健获胜报价的规则。

2 问题的提出

多属性逆向拍卖中，拍卖人的评价属性集记为G ={g ₁,g ₂,…,g _n }，其中g _i 表示第i 个属性。考虑到属性之间的关联性，引入模糊测度和模糊积分表示拍卖人的偏好。相关的基本知识^[25]介绍如下：

定义1：P (G )表示G 的幂集，如果集函数μ :P (G )→[0,1]满足以下非负有界性条件和单调性条件：

(a -1)μ (Ø)=0，μ (G )=1；

(1)采用压缩空气吹扫皮带，将细粉硫磺吹起再次形成粉尘，增加了吹扫点的粉尘浓度。尤其是细粉硫磺在皮带停运期间形成干粉，在皮带启动瞬间，压缩空气将干粉全部吹起，形成高浓度的粉尘云。

(a -2)∀A ,B ∈P (G )，A ⊆B ，有μ (A )≤μ (B )，

智慧交通是随着现代计算机技术和网络技术的不断发展而产生和建立的，使用了多种具有创新性的技术，智慧交通可以使物联网技术在交通行业中得到大范围的运用。智慧交通系统在进行建立和发展的过程中，需要有高新技术的支持，未来的发展方向也比较偏向于智能化交通系统的建立，并且在此基础上与先进的计算机和网络技术结合，从而更好地解决在交通系统中出现的各种现实问题。

则称μ 是G 上的模糊测度。

对每个属性集A ⊆G ，μ (A )可理解为A 的权重或重要程度。

定义2：模糊测度μ 可由其默比乌斯变换唯一地表示为：

(1)

其中，集函数m (B )反映出忽略任意g _i ∉B 与g _j ∈B 之间的关联时，集合B 中属性之间的关联程度。任意集函数m ：p (G )→不一定对应一个模糊测度的默比乌斯变换形式。根据条件(a -1)、(a -2)以及式(1)，以下两个条件保证了m (·)能唯一地确定集合G 上的一个模糊测度：

在取得现场检测的采集视频资料后，通过管道检测软件对检测视频进行判读并作出量化分析，得出山前水库涵洞存在缺陷统计，见表1。

即这与是MIP 2的最优解矛盾，因此假设不成立，得证。

k -可加模糊测度可表示任意复杂度的模糊测度。k 值越大，模型表现力越强，但参数越多。由于现实拍卖问题中三个及三个以上的属性之间的关联并不多见，本文为简化建模，采用定义在2-可加模糊测度上的Choquet积分作为拍卖人的评分函数。该积分由默比乌斯变换表示为^[26]：

(2)

其中，b 表示投标人的报价，g _i (b )表示该报价在属性g _i 上的评价值。

使用这一函数评价投标人报价需要确定参数m ({g _i })和m ({g _i ,g _j })的值。这样的参数共有个。由于认知局限性，拍卖人直接提供这些参数值是相当困难的。相比之下，要求拍卖人提供一些能反映其潜在评分规则的偏好信息更为容易。本文将利用这些信息推测与拍卖人偏好一致的参数值，并根据推测值确定获胜者。一般情况下，满足拍卖人偏好信息的参数值很多，即参数取值不唯一。采用不同评分函数可能造成获胜者选择结果不同。为了避免参数选择对决策结果造成的影响，文中考虑所有可能的参数取值，提出一种两阶段获胜者确定方法。首先根据所有可能的模糊测度取值挑选出可能获胜报价，再对这些潜在获胜报价进行进一步分析确定稳健获胜者。

该方法在多轮次拍卖的背景下提出。t =0表示拍卖的初始时刻，此时由拍卖人定义评价属性集G 以及提供偏好信息。根据拍卖人的主观偏好将推测与其一致的评分函数每一轮截止时刻之前，每个投标人都有一次机会投标。投标人是否会提交新报价取决于他能否获利，理性的投标人在不能获利时会退出竞价。在每一轮的结束时刻，对本轮提交的所有报价进行分析，并将分析得到的备选获胜报价集W _t 以及推荐的稳健获胜报价输出给拍卖人，由其确定本轮的最终获胜者拍卖人的选择将生成一个反馈信息，该信息可在下一轮推测中进一步精确参数的估计结果，进而提高推荐的准确性。如果某一轮在事先确定的响应时间内没有收到任何新报价，整个拍卖过程结束。这样，当前的临时获胜者即为本次拍卖的获胜者，拍卖人与获胜投标人根据报价进行实际交付。这一多轮次拍卖流程如图1所示。

图1 多轮次拍卖流程

3 两阶段获胜者确定模型

3.1 定义与偏好信息一致的评分函数

拍卖人的评分规则由一组偏好信息间接给出。这些偏好信息涉及到的报价称为参考报价，第t轮的参考报价集记为拍卖人可提供的偏好信息一般包括以下几种形式：(1)报价比较信息，对表示报价r _i 至少和r _j 一样好。拍卖人可明确区分≻和～关系时，r _i ≻r _j 表示报价r _i 优于r _j ，r _i ～r _j 表示报价r _i 与r _j 无差异；(2)报价偏好强度信息，对表示报价r _i 相对于r _j 的偏好强度属于第g 个偏好强度类。g 值越大表示偏好强度越弱。偏好强度类一般用语言形式表达，如“较强”，“一般”，“较弱”等；(3)属性权重信息，由线性不等式Aw ≥h 表示，其中属性权重向量w 满足w _i >0且属性关联类型信息，表示为(g _i ,g _j )^X ，其中X 为+，-和0时分别表示属性g _i 和g _j 之间存在互补性、冗余性和独立性；(5)属性关联程度比较信息，对g _i ,g _j ,g _k ,g _l ∈G ，Δ(g _i ,g _j )^X Δ(g _k ,g _l )^X 表示属性g _i 和g _j 之间的

关联程度大于g _k 和g _l 之间的关联程度。

由于y _ij 和v _ij 满足约束式y _ij -δ ≤Mv _ij 和y _ij ≥M (v _ij -1)+ε ₁，那么满足：

公共管理学以政府和其他公共组织的管理活动为研究对象, 其研究的范畴和领域非常广泛，从其现有的二级学科——社会保障、教育管理、公共卫生、土地资源分别着眼于不同的公共事务就能窥其一般。仅从社会保障来看，二级学科关注的话题几乎涉及到了国计民生的各个方面，但在课堂教学中，只能就养老保障、医疗保障、劳动就业、社会救助等方面的主要理论和基本政策措施进行传授，其他诸如教育保障、住房保障、军人保障、慈善失业等领域基本没有机会在课堂教学中涉及，因而导师的课外引领有助于学生将视野向更广更深的领域拓展。

(3)

其中，前三个不等式约束分别表示在拍卖人给出的三种形式的报价比较信息下报价得分所应满足的条件。δ 是拍卖人的无差异阈值，δ ≥0。若两个报价得分值之差低于该阈值，拍卖人不能区分它们的优劣。反之，若拍卖人明确给出优于关系，则得分值之差必然超过该阈值。因此，ε ₁取略大于δ 的数。第四个不等式反映的是属于较强偏好强度类的两个报价的得分差值也较大。属性权重满足第五个约束式，其中I _μ ({g _i })表示单个属性g _i 的相对重要程度。属性集A (A ⊆G )中各属性之间的关联程度用I (A )表示。由于文中只考虑属性两两关联，当|A |>2时，I (A )=0。第六至第八个约束式分别表示当属性之间存在互补性、替代性和相互独立时关联系数需满足的条件，其中ε ₂是较小的正数，目的在于保证严格不等式成立。第九个约束式对应属性关联程度比较信息。以上各式中的C _μ (·)和I (·)可由集函数m (·)表示。根据Grabisch对关联系数的定义及该系数与默比乌斯变换的关系^[27]，单个属性的相对重要程度以及属性两两关联系数分别表示为：

(4)

I ({g _i ,g _j })=m ({g _i ,g _j }),∀{g _i ,g _j }⊆G

(5)

最后两个约束式分别与条件(b -1)和(b -2)等价，保证了集函数m (·)能唯一地确定集合G 上的一个2-可加模糊测度。为方便所述，将这组约束集记为

3.2 获胜者选择分析

一般情况下，满足约束集的2-可加模糊测度有多组解，一组解对应一个评分函数。选择不同的评分函数可能得到截然不同的获胜者。因此，本文方法首先从报价集合B 中筛选出潜在的获胜报价。为便于表述，将所有一致的评分函数集合记为

3.2.1 可能获胜报价

定义4：如果使得对∀b _j ∈B\ {b _i }有C _μ (b _i )≥C _μ (b _j )-δ 成立，则称报价b _i 为可能获胜报价。

根据这一定义，建立如下的线性规划模型LP 1确定备选获胜报价集。

LP 1(b _i ):Max y _i

(6)

如果表明存在一个评分函数使得b _i 至少和其他报价一样好，因此b _i 有可能是获胜报价。对所有b _i ∈B 执行该线性规划，确定备选获胜报价集合，

定理1 如果b _i ∉W _t-1 ，必有b _i ∉W _t 。

例如，教师可以将学生分成不同小组，合作探究教材中的案例，自主搜集资料，从不同角度进行思考，无形中锻炼学生解决问题和分析整理的能力。此外，学生在小组学习中更容易积极表达自身的不同看法，在提高自主学习积极性和会计专业核心素养的同时，启迪智慧，将学生快速培养成社会所需人才。

证明：假设b _i ∉W _t-1 ，但b _i ∈W _t 。b _i ∈W _t 意味着使得对∀b _j ∈B\ {b _i }有C _μ (b _i )≥C _μ (b _j )成立。由于满足的评分函数C _μ 也必然满足即这样，当b _i ∈W _t 时必然使得对∀b _j ∈B\ {b _i }有C _μ (b _i )≥C _μ (b _j )，即b _i ∈W _t-1 ，这与假设相矛盾。因此假设不成立，得证。

推论1：如果b _i ∉W _t ，则对任意k >t ，b _i ∉W _k 。

+u ₂(i 74770)

3.2.2 稳健获胜报价

当存在多个可能获胜报价时，进一步对这些报价进行分析，选择其中的稳健获胜者推荐给拍卖人。稳健获胜报价是这样的一个报价：该报价为获胜者时其他的可能获胜报价也排在尽可能靠前的位置；当存在多个这样的报价使得其他可能获胜报价排序最靠前时，选择其中使得所有可能获胜报价与不可能获胜报价得分总差值最大的报价。图2以一个示例说明了挑选稳健获胜者的思路。该例中，报价集合中共有5个报价，其中b ₂、b ₃和b ₅为可能获胜报价。当b ₂为获胜者时，其他两个可能获胜报价b ₃和b ₅分别排在第4位和第2位；当b ₃为获胜者时，b ₂和b ₅分别排在第3位和第2位；当b ₅为获胜者时，b ₂和b ₃可能排在第2、3位，也可能排在第3、4位。可见，b ₃和b ₅为获胜者时可保证其他可能获胜报价排在最靠前的位置(即依次排在第2、3位)。接下来，再考虑二者分别为获胜者时所有可能获胜报价与不可能获胜报价之间的整体得分差异。在示例数据下，b ₅为获胜者时得分总差值更大，因此选择b ₅为稳健获胜者。在这一规则下，使稳健获胜报价得分最高的评分函数的评价结果与所有一致的评分函数的评价结果最为接近。在所有可能获胜报价都有最好的评价时，稳健获胜报价为最优中的最优。

图2 稳健获胜报价选择示意图

根据这一思路，建立如下的混合整数规划MIP 2挑选稳健获胜者。

b _j ∈B\W _t

(7)

其中，M 是一个特别大的正数。y _ij 表示可能获胜报价b _i 相比于不可能获胜报价b _j 的得分差值。v _ij 为对应的偏好关系标记，v _ij =1表示b _i ≻b _j ，v _ij =0表示约束y _ij -δ ≤Mv _ij 和y _ij ≥M (v _ij -1)+ε ₁可保证

(8)

目标函数为项加权重λ ，目的在于优先保证可能获胜报价优于所有不可能获胜报价，其次最大化得分差值。稳健获胜者

师资队伍是实施高校数据库课程体系改革的前提。目前，高等院校的数据库任课教师大多擅长传统的关系型数据库管理技术，但是缺乏对海量数据分析和数据系统实践的能力，这也使得师资队伍的建设成为实施数据库课堂体系改革的主要瓶颈。新的课程体系对于数据库任课教师的专业能力和水平提出了更高要求。因此，高等院校必须要根据海量数据分析和数据库系统实践的教学内容和要求，加强师资队伍建设。可以聘请企业精英作为讲师，也可以通过培训和自主学习提高现有师资的专业素养，从而提高高校数据库课程教学的团队水平。

定理2 当λ ≥|W _t |·|B\W _t |·(M +δ -ε ₁)-M 时，若是MIP 2的最优解，则也必然使最大。

证明：假设是MIP 2的最优解，而不是的最优解。那么，存在可行解使即

根据拍卖人所提供的偏好信息可推测得到2-可加模糊测度的参数值。如果由这些估计值所确定的偏好关系与拍卖人的主观判断相吻合，称这样的2-可加模糊测度与偏好信息相一致。这样，所有一致的2-可加模糊测度可由以下约束所定义：

这样，

当λ ≥|W _t |·|B\W _t |·(M +δ -ε ₁)-M 时，不等式两边同乘以

u ₁(3599)+u ₂(i 34170)+u ₃(4G )+u ₄(1000G )>u ₁(5499)+u ₂(i 74770)+u ₃(4G )+u ₄(1000G )

定义3: 如果一个模糊测度μ 的默比乌斯变换满足：对任意A ⊆G ，若|A |>k ，有m (A )=0，并存在A ⊆G ，|A |=k ，使得m (A )≠0，则该模糊测度被称为k-可加模糊测度。

稳健获胜报价作为推荐报价输出给拍卖人，由拍卖人最终决定哪个报价获胜。若拍卖人选择报价为获胜者，则其认为该报价优于其他可能获胜报价，则有如下反馈信息：

(9)

3.3 获胜者确定方法的基本步骤

通过上述分析，文中确定获胜者的基本步骤概括如下。

表2显示了在整个研究样本期间，上证综指收益率与各国(地区)指数收益率之间的相关系数。虽然传统相关系数测度具有一定的局限性，但是同样也能反映一些基本信息。无论是Pearson相关系数还是秩相关系数(Kendall相关系数和Spearman相关系数)，上证综指收益率与香港恒生指数收益率相关性最强，其次是日经225指数收益率，与标普500指数收益率相关性最弱。

步骤4：投标人进行投标；

步骤2：根据3.1节中的方法确定与拍卖人偏好一致的2-可加模糊测度集；

步骤3：检查拍卖结束条件是否已满足。若已满足，输出当前获胜者为最终获胜者，拍卖结束；若未满足，拍卖进入下一轮(t =t +1)；

步骤1：拍卖初始时刻(t =0)，拍卖人提出评价属性集，以及能反映自身偏好的参考报价信息和相应偏好信息；

这次拍摄能够如此顺利地结束，我们要特别感谢一位朋友，他就是Bristol道场的老板！向Simon O‘Brien先生致以诚挚的谢意，感谢他将整个道场借给我们使用，还为我们推荐了这两位如此出色的模特！如果你也喜欢武术格斗，请联系他吧！同时，也非常感谢两位模特精彩的配合，这给了我们的读者尽情创作的空间！总之，特殊的场地会带来新奇感，大家不妨也多多尝试！

步骤5：利用模型LP 1确定本轮的备选获胜报价集W _t ；

步骤6：利用模型MIP 2确定本轮的稳健获胜报价

步骤7：根据模型推荐结果，拍卖人选出本轮临时获胜者同时，相应的反馈信息生成，返回步骤2。

4 算例分析

4.1 仿真实验

某企业因业务需求要采购一批台式机用于日常办公。其采购部门拟通过一个网上拍卖项目采购这批商用台式机。根据企业的现实应用需求，该部门(即拍卖人)将以四个指标(即属性)评价供应商(即投标人)提供的台式机，分别是：价格(g ₁)、处理器性能(g ₂)、内存容量(g ₃)和硬盘容量(g ₄)。拍卖活动正式开始之前，他们考察了市场上部分商用台式机产品(即参考报价，如表1所示)，并向拍卖网站提交了一些关于这些台式机的评价信息(即偏好信息, 如表2所示)表达其采购偏好，期望拍卖网站根据这些信息筛选出符合其要求的报价。

表1 参考报价信息

注：产品信息来自太平洋电脑网产品报价频道

表2 采购部门提交的偏好信息

分析采购部门给出的报价比较信息，可以发现传统的加性效用函数不能刻画该采购部门的偏好。具体地，根据r ₅≻r ₂，有：

综上，必然有

(10)

同时，根据r ₆≻r ₁，有：

u ₁(10499)+u ₂(i 74770)+u ₃(2G )+u ₄(500G )>u ₁(3599)+u ₂(i 34170)+u ₃(2G )+u ₄(500G )

(11)

由(10)得到：

经济中心南移在两宋时期极为明显，其在一定程度上推动当时经济格局变化，也对后世经济态势变化产生一定影响。本次研究中对两宋时期的经济中心南移表现与南移的原因进行分析，同时指出经济中心南移所带来的影响，这些均反映出我国古代经济发展以及所做出的贡献不可估量。

一家地市级医院的发展空间有多大？在全国不少地市级医院切身感受到“夹心层”尴尬之时，沧州市中心医院在区域医疗版图中的引领角色已更加固化，并在新时代医疗势能与动能交汇下，打造出“大三环”医疗服务体系。

u ₁(3599)+u ₂(i 34170)>u ₁(5499)

该定理表明这种确定获胜者的方式保证了在某一轮不可能获胜的报价也不可能在以后各轮中获胜。因此，每一轮筛选可能获胜者时只需考虑上一轮的潜在获胜报价和本轮的新报价，即W _t ⊆(W _t-1 ∪B _t )。

(12)

而由(11)得到：

u ₁(3599)+u ₂(i 34170)<u ₁(10499)

+u ₂(i 74770)

(13)

根据边际价值函数u ₁(·)的单调非递增性质，即价格越高拍卖人的边际效用越低，有：

精密称取穿心莲内酯原料药、乳化剂及助乳化剂适量，置于锥形瓶中，于37℃水浴中溶解，振荡混匀；加入油相，于37℃水浴中振荡混匀，制得穿心莲内酯自微乳。另精密称取乳化剂及助乳化剂适量，置于锥形瓶中，于37℃水浴中溶解，振荡混匀；加入油相，于37℃水浴中振荡混匀，制得空白自微乳。

u ₁(10499)≤u ₁(5499)

(14)

因此不等式(12)和(13)明显地不一致，这表明拍卖人的偏好信息违反了偏好独立假设，传统的加性效用函数失效。

我国最早的道观藏书是简寂观道藏阁的藏书。南朝大明五年（公元461年），道教宗师陆修静在江西庐山建立简寂观，为庐山最大的道教宫观。陆修静在此修道、传教，整理道经与编撰道书，并编纂完成中国第一部道教经书总目《三洞经书目录》，共1090卷。简寂观是道教最大、最完备的藏经库。

本例采用Choquet积分表达采购部门的偏好。首先对初始偏好信息进行预处理。以权威IT产品评测机构Pconline发布的处理器横向测评结果为依据，采取3Dmake Vantage得分值作为处理器性能的评价值。再根据目前太平洋电脑网所有商用台式机产品在各属性上的取值范围，将参考报价各指标值转换到[0,100]统一标度上，0表示该指标上的最差得分，100表示最好得分。

以下采用仿真实验进行分析。可能对计算结果有影响的因素主要包括每次拍卖的轮数和每轮的报价数目，为方便表述，分别记为r和b。首先设定每一次采购进行3轮(即r=3)，每一轮参与竞价的报价数目10个(即b=10)。随机生成各报价在各属性上的评价值。重复进行100次拍卖(即重复计算100次)，可能获胜报价的比例如图3所示。单次运算中最多有36.67%的报价可能获胜，最少仅有13.33%，100次运算的平均值为21.93%。控制每一轮参与竞价的报价数目不变(b=10)，拍卖轮数由1轮增加到10轮，每种实验设置仍重复执行100次，平均可能获胜报价比例如图4所示。如图，可能获胜报价的比例基本上与拍卖轮数无关。在每轮报价数目为10的规模下，可能获胜报价所占的比例在22.14%左右。

风机基础作为风电机组的支撑体系，负责将基础过渡连接段传递的风机荷载和自身受到的波浪荷载、水流荷载以及靠泊力和撞击力等传递到地基土中，在保证海上风电系统整体安全方面起着十分重要的作用。在众多的海上风机基础形式中，单桩基础结构形式简单，受力明确，在具备打桩和起吊设备能力的前提下，施工速度快，在国内外风电场中得到广泛使用，占整个市场的65%以上。然而，为方便安装且同时保证电缆的安全，越来越多的安装基础选择将电缆管从桩基底部进入桩内，这就必须在钢管桩壁上开口，开孔导致桩体在开孔位置产生应大集中现象，此时钢管桩基础的疲劳强度问题将不容忽视。

图3 可能获胜报价比例

图4 不同拍卖轮数下平均可能获胜报价比例

控制每一次拍卖的轮数为3轮(r=3)，每一轮参与竞价的报价数目发生变化时，运算执行100次的平均可能获胜报价比例如图5所示。如图，可能获胜报价的比例将随着报价个数的增加而递减。当报价数目为3时，平均有52.44%的报价为可能获胜报价，但当报价数目增加到12个，可能获胜报价比例降低至18.58%。这与拍卖实践中的已有经验一致，当报价增多时，意味着供应商之间的竞争增强，将导致具有相对领先优势的报价减少。从图中还可以看到，这种递减趋势在报价数目持续增多时减缓，报价数目分别为10、11和12时，平均的可能获胜报价比例分别为21.93%、20.12%和18.58%。

图5 不同报价个数下平均可能获胜报价比例

综上，尽管可能获胜报价的比例会随着报价数目的变化而变化，但整体上，实验表明大量的报价为不可能获胜报价。即便在每一轮仅有3个报价的情形下，仍有超过45%的报价不可能获胜，且随着报价数目的增多，不可能获胜报价的比例持续增加，当一轮报价增加到10个时，不可能获胜报价高达80%。这说明本文提出的获胜者确定方法在第一阶段筛选出可能获胜报价非常必要，而且参与竞价的供应商越多，本文的方法越适用。

4.2 方法比较

Angilella 等^[28]提出了一种非可加稳健序数回归方法，该方法通过计算方案之间的必然和可能关系确定方案的优劣。根据必然偏好关系的定义，当对于都有C _μ (a )≥C _μ (b )+δ 时，认为a 必然比b 好，记为a ≻^N b 。如果对报价b _i ，存在报价b _j ∈B\ {b _i }，使得b _j ≻^N b _i 。那么在任意一致的评分函数下该报价b _i 都不可能获胜。因此，该方法通过两两比较也可确定出不可能获胜报价。本节将从计算效果和计算效率两方面对两种方法进行比较，为方便表述，将本文方法记为TSA，Angilella等提出的方法记为ROR。

设置一次拍卖进行3轮(r=3)，且保持相同的初始偏好信息、相同的随机报价信息和相同的反馈信息。报价数目在b=3到b=12之间变动时，两种方法在“可识别比例”和“未知报价比例”两个指标上的变化情况如表3所示。其中，“可识别比例”指的是在所有100次运算中，TSA可识别出由ROR方法确定的被占优报价的平均比例。从表中可以看到，由ROR两两比较所确定出的必然被占优报价100%地被TSA方法识别为不可能获胜报价，表明了TSA方法的有效性。此外，ROR方法的两两比较分析还可能产生一类“未知报价”，即既不必然优于其他报价、也不必然被其他报价占优的报价，这类报价无法确定其是否可能获胜。如表3所示，当一轮中参与竞价的报价数目较少时，未知报价的平均比例较高。因此，整体来说，在计算效果上TSA略优于ROR。

表3 不同报价数目下TSA的可识别比例和ROR的未知报价比例

为比较两种方法的计算效率，首先设置报价个数b=10，每次拍卖由1轮增加到10轮时，分别运算100次，两种方法平均每一轮的运行时间如图6所示。设置拍卖轮数r=3，报价数目由3个增加到12个时，同样执行运算100次，平均每一轮的运行时间如图7所示。可以看到，在r和b变化时，两种方法的运行时间变化呈现相似的规律。首先，TSA和ROR方法的平均运行时间都随着r和b的增加而增大。这是因为，当拍卖轮数增加时，各轮要处理的反馈信息增多；当报价数目增多时，每一轮的计算量增加。虽然两种方法都呈现出递增趋势，但在相同的问题规模下，TSA花费的时间代价总是比ROR小。尤其是当拍卖轮数和报价数目较大时，TSA在计算效率上的优势更为明显。

图6 拍卖轮数变化时平均每轮运行时间

图7 报价数目变化时平均每轮运行时间

综上，无论是在计算效果还是计算效率上，本文方法均优于Angilella等的方法。通过仿真实验，表明了本文方法更适用于拍卖轮数较多、报价规模较大的拍卖项目。

5 结语

多属性逆向拍卖日益广泛地应用在采购活动中，其中的关键问题是获胜者确定问题。目前主要的方法是采用评分函数进行打分，但是现有研究假设属性不相关，具有局限性。本文根据模糊测度和模糊积分理论，提出了一种考虑属性两两关联的获胜者确定方法。由于参数值众多且拍卖人难以直接提供，文中从拍卖人提供的偏好信息中推测Choquet积分形式的评分函数。与偏好信息一致的评分函数可能不唯一，不同的评分函数可能导致不同的决策结果，因此文中考虑了所有一致的评分函数集确定获胜者。首先根据拍卖人的主观偏好定义了与其一致的2-可加模糊测度集，其次根据该评分函数集合挑选出备选获胜报价集，最后以与所有一致评分函数的评价结果最为接近的报价排序为依据，选出得分最高者作为稳健获胜报价推荐给拍卖人。拍卖人确认最终获胜者，产生的反馈信息进入下一轮的偏好信息集，进一步精确方法的推荐结果。

经过随机数据验证，拍卖中不可能获胜报价的比例远高于可能获胜报价比例，且随着供应商竞价增多呈增长趋势，因此在第一阶段筛选出少量的可能获胜报价非常重要。与Angilella等提出的方法进行比较发现，文中方法可将全部的必然被占优报价识别为不可能获胜报价，且Angilella等的方法总会存在一小部分的报价无法确知其是否可能获胜，因此在计算效果上本文方法略优。从运算效率上，在相同的问题规模下文中方法花费的时间更短，且随着拍卖轮数和报价数目的增多，该方法的优势更为明显。

参考文献：

[1] Bichler M, Kalagnanam J. Configurable offers and winner determination in multi-attribute auctions [J]. European Journal of Operational Research, 2005, 160(2): 380-394.

[2] Kameshwaran S, Narahari Y, Rosa CH, et al. Multiattribute electronic procurement using goal programming [J]. European Journal of Operational Research, 2007, 179(2): 518-536.

[3] Cheng Chibin. Solving a sealed-bid reverse auction problem by multiple-criterion decision-making methods [J]. Computers and Mathematics with Applications, 2008, 56(12): 3261-3274.

[4] Ray A K, Jenamani M, Mohapatra P K J. An efficient reverse auction mechanism for limited supplier base [J]. Electronic Commerce Research and Applications, 2011, 10(2): 170-182.

[5] Ray A K, Jenamani M, Mohapatra P K J. Supplier behavior modeling and winner determination using parallel MDP [J]. Expert Systems with Applications, 2011, 38(5): 4689-4697.

[6] 王明喜, 谢海滨, 胡毅. 基于简单加权法的多属性采购拍卖模型[J]. 系统工程理论与实践, 2014, 34(11): 2772-2782.

[7] Pham L, Teich J, Wallenius H, et al. Multi-attribute online reverse auctions: Recent research trends [J]. European Journal of Operational Research, 2015, 242(1): 1-9.

[8] Talluri S, Narasimhan R，Viswanathan S. Information technologies for procurement decisions: A decision support system for multi-attribute e-reverse auctions [J]. International Journal of Production Research, 2007, 45(11): 2615-2628.

[9] 呼大永, 冯玉强, 唐振宇, 等. 基于自组织神经网络和DEA的采购拍卖获胜者确定问题模型[J]. 系统工程理论与实践, 2012, 32(2): 398-404.

[10] Huang Min, Qian Xiaohu, Fang Shu-Cherng, et al. Winner determination for risk aversion buyers in multi-attribute reverse auction [J]. Omega, 2016, (59): 184-200.

[11] 刘卫锋, 常娟, 杜迎雪.(语言)Heronian平均算子及其决策应用[J]. 中国管理科学, 2017, 25(7): 174-183.

[12] Angilella S, Corrente S, Greco S, et al. MUSA-INT: Multicriteria customer satisfaction analysis with interacting criteria [J]. Omega, 2014, 42(1): 189-200.

[13] Greco S, Mousseau V, Sowiński R. Robust ordinal regression for value functions handling interacting criteria [J]. European Journal of Operational Research, 2014, 239(3): 711-730.

[14] 陈岩, 李庭, 鲍博. 基于Choquet积分的指标关联模糊多目标指派问题[J]. 系统工程理论与实践, 2017, 37(8): 2162-2170.

[15] 杨洁, 李登峰. 指标相关且信息模糊的合作创新伙伴选择方法[J]. 运筹与管理, 2016, 25(5): 53-58.

[16] 赵树平, 梁昌勇, 罗大伟. 基于VIKOR和诱导广义直觉梯形模糊Choquet积分算子的多属性群决策方法[J]. 中国管理科学, 2016, 24(6): 132-142.

[17] 章玲, 周德群. 基于k-可加模糊测度的多属性决策分析[J]. 管理科学学报, 2008, 11(6): 18-24.

[18] 武建章, 张强. 基于2-可加模糊测度的多准则决策方法[J]. 系统工程理论与实践, 2010, 30(7): 1229-1237.

[19] 常志朋, 程龙生. 灰模糊积分关联度决策模型[J]. 中国管理科学, 2015, 23(11): 105-111.

[20] Teich J E, Wallenius H, Wallenius J, et al. Emerging multiple issue e-auctions [J]. European Journal of Operational Research, 2004, 159(1): 1-16.

[21] 于红岩, 刘仲英. 基于拍卖方偏好揭示的多属性网上拍卖模型[J]. 系统工程, 2008, 26(3): 54-58.

[22] Karakaya G, Köksalan M. An interactive approach for multi-attribute auctions [J]. Decision Support Systems, 2011, 51(2): 299-306.

[23] Yang Na, Liao Xiuwu, Huang Wayne Wei. Decision support for preference elicitation in multi-attribute electronic procurement auctions through an agent-based intermediary [J]. Decision Support Systems, 2014, (57): 127-138.

[24] 周正龙, 马本江, 胡凤英. 随机需求条件下的P2P网络借贷拍卖机制[J]. 中国管理科学, 2018, 26(5): 21-30.

[25] Bonetti A, Bortot S, Fedrizzi M, et al. Modelling group processes and effort estimation in project management using the Choquet integral: An MCDM approach [J]. Expert Systems with Applications, 2012, 39(18): 13366-13375.

[26] Grabisch M, Labreuche C. A decade of application of the Choquet and Sugeno integrals in multi-criteria decision aid [J]. Annals of Operations Research, 2010, 6(1): 1-44.

[27] Grabisch M. K-order additive discrete fuzzy measures and their representation [J]. Fuzzy Sets and Systems, 1997, 92(2): 167-189.

[28] Angilella S, Greco S, Matarazzo B. Non-additive robust ordinal regression: A multiple criteria decision model based on the Choquet integral [J]. European Journal of Operational Research, 2010, 201(1): 277-288.

Two Stage Winner Determination Approach with Interacting Attributes

YANG Na ¹,LIAO Xiu -wu ²,LEI Hong -zhen ¹

(1. International Business School, Shaanxi Normal University, Xi’an 710119, China;2. School of Management, The Key Lab of the Ministry of Education for Process Control & Efficiency Engineering Xi’an Jiaotong University, Xi’an 710049, China)

Abstract ： Winner determination problem (WDP) is one of the most important research issues in multi-attribute reverse auctions. Presently, the additive value function is commonly used as the scoring rule to solve WDP, which ignores the interactions among the attributes in reality. To deal with the interacting effects, the Choquet integral defined on 2-order additive fuzzy measures is used in this paper to represent an auctioneer’s preference, where the parameters are inferred by a set of indirect preference information in accordance with the preference aggregation-disaggregation paradigm. Generally, the scoring functions which could restore the auctioneer’s preference are not unique. A two stage winner determination model is then proposed, which aimes at utilizing all the compatible scoring functions to determine the winner. In each round of an iterative auction, all the potential winning bids are first selected via a linear programming model. Then, a mixed integer linear programming model is developed to determine a ranking which is most close to the evaluations of all the compatible scoring functions. The bid with the highest score in this ranking will be the robust winning bid, which will be recommended to the auctioneer. The auctioneer’s final decision will generate a feedback information to further refine the elicitation. Finally, a simulation experiment is conducted to evince the practical applicability of the method. The simulation results show that in general a large amount of bids are impossible to win, and the percentage increases with the number of bids in each round, which illustrates the necessity of choosing the potential winning bids in the first stage. The comparison with an existing approach verifies the validity and efficiency of the method. According to the simulation, the proposed method can identify all the outranked bids determined by the existing approach as losing bids. Moreover, unlike the existing method which will generate a set of “unknown bids”, the proposed one can distinguish the potential winning bids from the offers. In computation efficiency, the running time of the proposed method is always less than the existing one in the same auction setting. Although both two methods will cost more time with larger rounds or more bids, the running time of the existing method increases more sharply, which proves the advantage of the proposed method in computation efficiency.

Key words ： multi-attribute reverse auctions; winner determination problem; Choquet integral; preference elicitation

文章编号： 1003-207( 2019) 10-0148-11

DOI: 10.16381/ j.cnki.issn1003-207x.2019.10.015

中图分类号： C934; F724.59

文献标识码： A

收稿日期： 2017-12-05； 修订日期： 2018-06-11

基金项目： 国家自然科学基金资助项目(91846110,71872144)；陕西省软科学研究计划资助项目(2018KRM051)；陕西省哲学社会科学基金资助项目(2018S25)；陕西师范大学中央高校基本科研业务费专项资金资助项目(16SZYB36)

通讯作者简介： 杨娜(1987-)，女(汉族)，陕西白水人，陕西师范大学国际商学院讲师，博士，研究方向：多准则决策、多属性电子拍卖，E-mail: yang.na@snnu.edu.cn.

标签：多属性逆向拍卖论文;  获胜者确定问题论文;  Choquet积分论文;  偏好揭示论文;  陕西师范大学国际商学院论文;  西安交通大学管理学院过程控制与效率工程教育部重点实验室论文;