CN109389281A - 一种基于高斯混合模型的采集终端生产厂商评价方法 - Google Patents

Abstract

本发明公开了一种基于高斯混合模型的采集终端生产厂商评价方法，包括步骤：获取采集终端的原始数据，筛选衡量采集终端的评估指标；采用高斯混合模型聚类算法提取典型评估指标数据，获取聚类中心，形成决策矩阵，并对其分别按效益型和成本型进行标准化处理；分别计算采集终端各项评估指标基于层次分析法、熵权法、皮尔逊相关系数法以及变异系数法的权重，用指标权重的组合优化模型获得组合权重，对决策矩阵进行加权处理从而获得采集终端综合评估的评估矩阵；用夹角度量法对采集终端的质量进行综合评估，按照从高至低的顺序将各个采集终端供应商的质量进行优劣排序，获得评价结果。其可以为用用电信息采集系统采集运维工作提供支撑。

Description

一种基于高斯混合模型的采集终端生产厂商评价方法

技术领域

本发明涉及电力系统领域，特别是涉及一种基于高斯混合模型的采集终端生产厂商评价 方法。

背景技术

自2011年国网开始大范围建设用电信息采集系统以来，用电信息采集系统越来越多得承 担了营销侧的计量业务工作，采集终端定时定期地采集电能表数据项，按要求将采集信息上 送系统，是用电信息采集系统的核心组成部分，其质量直接影响着用户用电信息的获取，针 对其故障导致的运维也浪费了大量的人力物力。

用电信息采集系统覆盖范围的逐步扩大，采集的对象和采集频率逐年提高，各类计量数 据和采集数据应用逐渐完善。以浙江的用电信息采集系统为例，用电信息采集每周数据增长 量近1TB，历年累积收集的数据十分丰厚，并且在2015年建立的采集运维模块收集了大量的 采集终端故障历史信息。但针对这些数据却缺少足够的数据挖掘手段，使数据金矿的价值难 以体现。目前浙电采集大数据分析平台已经具备初步的数据分析能力，能够通过大数据技术 分析计量装置运行数据，开展可靠性质量评估方法研究，进而提升公司计量精益化管理水平。

供电单位一般对采集终端的厂家运行质量评估停留在一些考核指标的基础上，以及针对 地市单位反映的批量性故障收集上，相对比较粗犷。

发明内容

本发明要解决的技术问题和提出的技术任务是对现有技术方案进行完善与改进，提供一 种基于高斯混合模型的采集终端生产厂商评价方法，以达到对采集终端质量作出科学、合理 的定量评价的目的。为此，本发明采取以下技术方案。

包括如下步骤：

1)获取每个批次和区域中采集终端的原始数据，确定衡量采集终端的8个评估指标：8 个评估指标分别为：负荷采集合格率α₁，电量数据采集率α₂，无故障工作时间α₃，反映采集 异常的严重异常次数α₄、一般异常次数α₅、轻微异常次数α₆，维修率α₇和在线率α₈：

式中：N_samp是采样时间段的总个数，T_samp,i是第i次采样的时间长度，即总采 样时间；M_i、m_i、C_i、c_i、b_i和B_i分别是第i个采样时间段内采集终端负荷数据总个数、负荷数据合格的个数、理论应采集的数据量、实际采集的数据量、在线的采集终端个数和采集终端总个数；F_ser,i、F_gen,i、F_un,i分别是第i个采样时间段内采集终端发生严重、一般和轻微故障的次数；N_samp,F是首次出现F_ser,i+F_gen,i+F_un,i≠0的采样时间段；T_rate和T_F分别是采集终端的额 定运行时间和故障停运时间；

2)采用高斯混合模型聚类算法提取采集终端海量数据中的评估指标数据，获取聚类中心， 形成采集终端综合评估的决策矩阵，并对其分别按效益型和成本型进行标准化处理；

3)分别计算采集终端各项典型评估指标基于层次分析法、熵权法、皮尔逊相关系数法以 及变异系数法的权重，然后用指标权重的组合优化模型获得组合权重，对决策矩阵进行加权 处理从而获得采集终端综合评估的评估矩阵；

4)根据评估矩阵，用夹角度量法对采集终端的质量进行综合评估，按照从高至低的顺序 将各个采集终端供应商的质量进行优劣排序，获得评价结果。

作为优选技术手段：在步骤2)高斯混合模型聚类算法中：

假设GMM由K个高斯分布混合而成，则每个高斯分布称作一个“组分”，这些“组分”线 性相加即为GMM的概率密度函数：

式中：ω_k表示第k个多维单高斯分布的权值，N_k(x；μ_k；Σ_k)表示第k个多维单高斯分布的概 率密度函数，x表示数据样本列向量，μ_k表示第k个高斯模型的期望向量，Σ_k表示第k个高 斯模型的方差；

采集终端数据的GMM聚类算法流程为：

205)令l＝0，随机选择初始化的ω^(l),μ^(l),Σ^(l)，求取第i个样本点x_i属于第k类高斯模 型的初始化后验概率：

206)将第201)步得到的结果代入最大似然公式计算第l+1次迭代的GMM聚类参数：

式中：N为待聚类的数据点个数；

207)若成立(ε一般取10^-5)，则迭代结束，即为估计得到的GMM聚类参数，继续下一步；否则，回到第202)步；

208)采用贝叶斯概率公式计算第i个样本点x_i属于第k类高斯模型的概率：

根据贝叶斯概率最大准则，将第i个样本点x_i划分到使其概率p(μ_k,Σ_k|x_i)取得最大的那 类高斯模型中。

作为优选技术手段：在步骤2)，决策矩阵为：

式中：d_ij表示第i个采集终端供应商第j个指标的值，P为采集终端供应商个数，Q为衡量供 应商的采集终端的评价指标个数，其中Q等于8；

效益型指标标准化处理方法为：

成本型指标标准化处理方法为：

式中：和分别表示所有P个采集终端供应商中指标j的最小值和最大值；Ω₁表示效益型指标集合，其中Ω₁＝{1,2,3,8}；Ω₂表示成本型指标集合，其中Ω₂＝{4,5,6,7}。

作为优选技术手段：在步骤3)中，分别采用层次分析法、熵权法、相关系数法和变异 系数法计算采集终端指标权重，接着用指标的组合优化模型计算组合权重，其中：

层次分析法为：设Q个指标为U＝{u₁,u₂,…,u_Q}，每次取两个指标u_i和u_j，按1～9的比 例用a_ij表示指标u_i与u_j对采集终端的影响程度之比，数值越大表示u_i相对于u_j越重要；这 样即可获得一个采集终端评估问题的判断矩阵A＝(a_ij)_Q×Q，该判断矩阵具有如下性质：a_ij>0,a_ji＝1/a_ij,(i＝1,2,…,Q；j＝1,2,…,Q)；计算指标权重前首先对判断矩阵A进行一致性检验，如果一 致性检验未通过，则说明a_ij取值前后矛盾，需要重新赋值；如果通过，则说明a_ij取值前后一 致，判断矩阵A可用于求取指标权重；求取指标权重可以用列和求逆法，即：

将b_j归一化，即可求得归一化后的指标权重为

式中：且

采集终端指标的熵权为：

式中：κ表示常数κ＝1/ln P,并且假定f_ij＝0时,f_ijln f_ij＝0；且

计算皮尔逊相关系数权重包括：假设d′_x＝(d′_1x,d′_2x,...,d′_Px)^T和d′_y＝(d′_1y,d′_2y,...,d′_Py)^T表示标 准化决策矩阵D′的两个列向量，则第x个指标和第y个指标之间的皮尔逊相关系数定义为

式中：和分别是向量d′_x和d′_y的平均值；基于皮尔逊相关系数的权重定义为

式中：且

计算变异系数法权重包括：指标u_j的变异系数定义为：

式中：表示向量(d′_1j,d′_2j,...,d′_Pj)^T的平均值；基于变异系数法的指标权重定义为：

式中：且

指标权重的组合优化模型为：

式中：表示待优化的第j个采集终端指标的组合权重，表示采用第k种方法求解得到的 第j个指标的权重，G表示权重确定方法的数目，表示的隶属度函数；显然，越大，权重与最终的组合权重越贴近。

作为优选技术手段：在步骤4)用夹角度量法对采集终端的质量进行综合评估时，包括 步骤：

405)形成P个待评估的供应商对应的Q个质量指标的评估矩阵

R＝(r_ij)_P×Q

式中：

406)计算采集终端评估中的理想点和负理想点，其分别为

式中：

407)分别计算每个供应商下的采集终端指标与理想点、负理想点的夹角距离，即

式中：r_i＝(r_i1,r_i2,...,r_iQ)表示评估矩阵R的第i个行向量；

408)计算每个供应商下的采集终端指标与理想点的夹角逼近程度，即

可以看出：第i个供应商生产的采集终端越好，γ_i愈接近于1；反之，质量越差，γ_i越接 近于0；因此，可以根据γ_i取值的大小得到采集终端的评估结果。

有益效果：本技术方案及时利用其用采多年存储的大数据，建立可靠的质量评估指标体 系，对采集终端质量作出科学、合理的定量评价，以提供物资招标的信息支持。

附图说明

图1为本发明流程图；

具体实施方式

为了更好地理解本发明的目的、技术方案以及技术效果，以下结合附图对本发明进行进 一步的讲解说明。

参考图1，图1所示为本实施例的基于计量大数据聚类模型的采集终端生产厂商评价方 法流程图，包括如下步骤：

S10，获取采集终端的原始数据，筛选8个衡量采集终端的评估指标：负荷采集合格率， 电量数据采集率，无故障工作时间，反映采集异常的严重异常次数、一般异常次数、轻微异 常次数，维修率和在线率；

在本实施例中：

采集终端分为供应商、型号、批次和单个设备四个层次，而计量自动化系统中一般以单 个设备为单元存储数据。目前，在计量自动化系统中已采集的采集终端数据包含有负荷采集 合格率α₁，电量数据采集率α₂，无故障工作时间α₃，反映采集异常的严重异常次数α₄、一 般异常次数α₅、轻微异常次数α₆，维修率α₇和在线率α₈，其含义分别可以为：

式中：N_samp是采样时间段的总个数，T_samp,i是第i次采样的时间长度，即总采样时 间；M_i、m_i、C_i、c_i、b_i和B_i分别是第i个采样时间段内采集终端负荷数据总个数、负荷数据 合格的个数、理论应采集的数据量、实际采集的数据量、在线的采集终端个数和采集终端总 个数；F_ser,i、F_gen,i、F_un,i分别是第i个采样时间段内采集终端发生严重、一般和轻微故障的次 数；N_samp,F是首次出现F_ser,i+F_gen,i+F_un,i≠0的采样时间段；T_rate和T_F分别是采集终端的额定运 行时间和故障停运时间。

S20，采用大数据分析技术中的高斯混合模型聚类算法提取采集终端海量数据中的典型指 标数据，获取聚类中心，形成采集终端综合评估的决策矩阵，并对其分别按效益型和成本型 进行标准化处理；

在本实施例中：

每个供应商生产的采集终端数量众多、型号复杂，因此从系统中导出的采集终端数据无 法直接应用于对每个供应商的质量评估。为此，可以用大数据分析技术中的数据聚类方法， 挖掘各个供应商采集终端的典型特征，从而实现对其质量的综合评估。

高斯混合模型(Gaussian Mixture Model,GMM)用高斯分布概率密度函数来描述数据，它 可以将一个大数据集分解为若干个符合高斯分布的模型。因此，高斯混合模型可用于数据聚 类，提取数据集中的典型特征。假设GMM由K个高斯分布混合而成，则每个高斯分布称作 一个“组分”，这些“组分”线性相加即为GMM的概率密度函数：

式中：ω_k表示第k个多维单高斯分布的权值，N_k(x；μ_k；Σ_k)表示第k个多维单高斯分布的概 率密度函数，x表示数据样本列向量，μ_k表示第k个高斯模型的期望向量，Σ_k表示第k个高 斯模型的方差。

采集终端的质量受到多个独立因素的影响，且这些因素都不是决定性的。根据中心极限 定理，可认为同一批次的采集终端的指标数据符合高斯分布。对于某一采集终端供应商下的 K个批次的设备，可认为其质量的指标数据是K个高斯分布模型的叠加。因此，高斯混合模 型聚类算法适用于挖掘采集终端的指标数据。在进行GMM聚类之前，首先需要估计各个高 斯分布模型的参数，这里采用期望最大化算法进行估计。因此，采集终端数据的GMM聚类 算法流程可以为：

201)令l＝0，随机选择初始化的ω^(l),μ^(l),Σ^(l)，求取第i个样本点x_i属于第k类高斯模 型的初始化后验概率：

202)将第201)步得到的结果代入最大似然公式计算第l+1次迭代的GMM聚类参数：

式中：N为待聚类的数据点个数。

203)若成立(ε一般取10^-5)，则迭代结束，即为估计得到的GMM聚类参数，继续下一步；否则，回到第202)步；

204)采用贝叶斯概率公式计算第i个样本点x_i属于第k类高斯模型的概率：

根据贝叶斯概率最大准则，将第i个样本点x_i划分到使其概率p(μ_k,Σ_k|x_i)取得最大的那 类高斯模型中。

当所有的样本点均按照上述步骤划分到各自所属的高斯模型后，GMM聚类完成并获得 了相应的聚类中心。GMM聚类算法舍弃了采集终端评估原始指标数据中冗余和不重要的信 息，但保留了原始指标数据的重要数据和典型特征。因此，GMM聚类算法大大压缩了数据 的规模、减少了综合评估的计算量，却几乎不影响采集终端综合评估的准确性。

设D为采集终端评估问题的决策矩阵，D中的元素d_ij表示第i个采集终端供应商第j个 指标的值。因为不同指标间的量纲不同，所以不能直接对各个指标的重要程度进行比较，需 要进行归一化处理。此外，采集终端指标还分为效益型和成本型两类，效益型指标数值越大 代表质量越好，成本型指标数值越大代表质量越差。因此，在确定指标权重前首先要对D进 行归一化处理以得到归一化后的决策矩阵D′，其归一化处理公式可以为：

式中：和分别表示所有P个采集终端供应商中指标j的最小值和最大值；Ω₁表示效益型指标集合，其中Ω₁＝{1,2,3,8}；Ω₂表示成本型指标集合，其中Ω₂＝{4,5,6,7}。

S30，分别计算采集终端各项指标基于层次分析法、熵权法、皮尔逊相关系数法以及变异 系数法的权重，然后用指标权重的组合优化模型获得组合权重，对决策矩阵进行加权处理从 而获得采集终端综合评估的评估矩阵；

在本实施例中：

AHP法通过指标间两两比较重要程度，间接获得每个指标的权重，该方法属于一种主观 赋权法。设Q个指标为U＝{u₁,u₂,…,u_Q}。每次取两个指标u_i和u_j，按1～9的比例用a_ij表 示指标u_i与u_j对采集终端的影响程度之比，数值越大表示u_i相对于u_j越重要。这样即可获得 一个采集终端评估问题的判断矩阵A＝(a_ij)_Q×Q，该判断矩阵具有如下性质：a_ij>0,a_ji＝1/a_ij,(i ＝1,2,…,Q；j＝1,2,…,Q)。计算指标权重前首先对判断矩阵A进行一致性检验，如果一致性检 验未通过，则说明a_ij取值前后矛盾，需要重新赋值；如果通过，则说明a_ij取值前后一致，判 断矩阵A可用于求取指标权重。即可以为：

将b_j归一化，即可求得归一化后的指标权重为

式中：且

在信息学中，熵用来表征系统的无序程度和数据的离差程度。指标的信息熵越大则其在 综合评估中提供的信息就越少，相应的权重也应该越小；反之，其信息熵越小，则权重应当 越大。因此，采集终端指标u_j的熵H_j可以定义为

式中：κ表示常数κ＝1/ln P,并且假定f_ij＝0时,f_ijln f_ij＝0。据此，第j个 采集终端指标的熵权可以定义为

式中：且

在统计学中相关系数用于衡量两个变量之间的关联性，包括皮尔逊相关系数、斯皮尔曼 相关系数和肯德尔相关系数等，其中皮尔逊相关系数更适用于符合正态分布的数据，因此本 发明使用皮尔逊相关系数来计算指标间内在的相关程度。相关程度越大说明指标间信息的重 复性越大，相应指标的权重应越小；反之，相关程度越小，权重应越大。假设d′_x＝(d′_1x,d′_2x,...,d′_Px)^T和d_y′＝(d′_1y,d′_2y,...,d′_Py)^T表示标准化决策矩阵D′的两个列向量，则第x个 指标和第y个指标之间的皮尔逊相关系数可以定义为：

式中：和分别是向量d′_x和d′_y的平均值。因此，基于相关系数法的权重可以定义为

式中：且

变异系数法是一种根据指标数据间的对比强度来确定指标权重的客观赋权法。某个指标 的变异程度越大，说明其相对于其它指标的对比强度越大，则在综合评估中对评估对象的重 要性越高，从而其权重也应越大；反之，变异程度越小，权重应越小。指标u_j的变异系数可 以定义为：

式中：表示向量(d′_1j,d′_2j,...,d′_Pj)^T的平均值。因此，基于变异系数法的指标权重可以定义为

式中：且

虽然上述4种主、客观方法简单易行，但均存在一定的不足。主观赋权法仅依照专家的 经验，很容易导致主观偏好过强；客观赋权法仅凭借实际数据，但实际数据可能会出现一定 误差，从而可能导致评估结果不符合实际。为了充分考虑依据专家经验的主观打分，同时依 据数据本身的特征对权重进行修正，本发明基于权重隶属度最大，构建了采集终端指标权重 的组合优化模型，即可以为

式中：表示待优化的第j个采集终端指标的组合权重，表示采用第k种方法求解得到的 第j个指标的权重，G表示权重确定方法的数目，表示的隶属度函数。显然，越大，权重与最终的组合权重越贴近。

S40，根据评估矩阵，用夹角度量法对采集终端的质量进行综合评估，按照从高至低的顺 序将各个采集终端供应商的质量进行优劣排序，获得评价结果。

在本实施例中：

确定各个采集终端指标的组合权重后，采集终端的综合质量评估可以转化为一个多属性 决策问题。夹角度量法采用变量之间的夹角作为距离的测度，依据被评估对象与理想化目标 的逼近程度进行排序。因此，基于夹角度量法的采集终端综合质量评估的流程可以为：

401)形成P个待评估的供应商对应的Q个质量指标的评估矩阵

R＝(r_ij)_P×Q

式中：

402)计算采集终端评估中的理想点和负理想点，其分别为

式中：

403)分别计算每个供应商下的采集终端指标与理想点、负理想点的夹角距离，即

式中：r_i＝(r_i1,r_i2,...,r_iQ)表示评估矩阵R的第i个行向量。

404)计算每个供应商下的采集终端指标与理想点的夹角逼近程度，即

可以看出：第i个供应商生产的采集终端越好，γ_i愈接近于1；反之，质量越差，γ_i越接 近于0。因此，可以根据γ_i取值的大小得到采集终端的评估结果。

为了进一步理解本发明，以下采用国网浙江省电力公司宁波供电公司管辖的某地区的采 集终端数据进行算例仿真，该原始数据集共有11565条数据，经过数据清洗之后可用数据为 11312条，共有17个待评估的采集终端供应商，每个供应商下有若干设备批次，所有的设备 批次总数为58个。

将17个供应商及其58个批次的数据按顺序重新进行编号，以第1、3、12个供应商为例 分析GMM聚类算法的效果。第1个供应商下有2个批次，共142台采集终端；第3个供应 商下有4个批次，共496台采集终端，第12个供应商下有6个批次，共1514台采集终端。 每台采集终端都有8个指标。

对每个供应商下的采集终端都进行GMM聚类分析，可以获得表征各个供应商的采集终 端的典型采集终端，然后将这些典型采集终端相应的指标取平均值即可得到如表1所示的决 策矩阵D′。

表1标准化后的采集终端决策矩阵

表2分别给出了基于AHP、基于熵权法、相关系数法、变异系数法和组合优化模型的指 标权重。从表2可以看出：客观权重对主观权重起到了一定的修正作用，按照专家经验，采 集终端的严重异常指标α₄应占有较大的权重，但另外三种客观赋权法都对α₄给出了较小的权 重，这是因为原始数据中各个供应商的α₄数值差异不大。因此，和其它指标相比，该指标难 以辨别各个供应商的采集终端优劣，故应当赋予较小的权重。

表2不同方法计算得到的指标权重

然后，基于得到的指标组合权重，形成采集终端的评估矩阵R，进而采用夹角度量法对 采集终端的质量进行综合评估，最后得到如表3所示的采集终端综合评估的结果。

表3采集终端综合评估结果

从表3可以看出：采集终端最好的前9家采集终端供应商分别为：10、9、6、5、7、13、4、11和16，其中第10个供应商在采集终端综合评估中的评估值最高，第8个供应商的评估值最低。

Claims (5)

1.一种基于高斯混合模型的采集终端生产厂商评价方法，其特征在于，包括如下步骤：

1)获取每个批次和区域中采集终端的原始数据，确定衡量采集终端的8个评估指标：8个评估指标分别为：负荷采集合格率α₁，电量数据采集率α₂，无故障工作时间α₃，反映采集异常的严重异常次数α₄、一般异常次数α₅、轻微异常次数α₆，维修率α₇和在线率α₈：

式中：N_samp是采样时间段的总个数，T_samp,i是第i次采样的时间长度，即总采样时间；M_i、m_i、C_i、c_i、b_i和B_i分别是第i个采样时间段内采集终端负荷数据总个数、负荷数据合格的个数、理论应采集的数据量、实际采集的数据量、在线的采集终端个数和采集终端总个数；F_ser,i、F_gen,i、F_un,i分别是第i个采样时间段内采集终端发生严重、一般和轻微故障的次数；N_samp,F是首次出现F_ser,i+F_gen,i+F_un,i≠0的采样时间段；T_rate和T_F分别是采集终端的额定运行时间和故障停运时间；

2)采用高斯混合模型聚类算法提取采集终端海量数据中的评估指标数据，获取聚类中心，形成采集终端综合评估的决策矩阵，并对其分别按效益型和成本型进行标准化处理；

3)分别计算采集终端各项典型评估指标基于层次分析法、熵权法、皮尔逊相关系数法以及变异系数法的权重，然后用指标权重的组合优化模型获得组合权重，对决策矩阵进行加权处理从而获得采集终端综合评估的评估矩阵；

4)根据评估矩阵，用夹角度量法对采集终端的质量进行综合评估，按照从高至低的顺序将各个采集终端供应商的质量进行优劣排序，获得评价结果。

2.根据权利要求1的一种基于高斯混合模型的采集终端生产厂商评价方法，其特征在于，在步骤2)高斯混合模型聚类算法中：

假设GMM由K个高斯分布混合而成，则每个高斯分布称作一个“组分”，这些“组分”线性相加即为GMM的概率密度函数：

式中：ω_k表示第k个多维单高斯分布的权值，N_k(x；μ_k；Σ_k)表示第k个多维单高斯分布的概率密度函数，x表示数据样本列向量，μ_k表示第k个高斯模型的期望向量，Σ_k表示第k个高斯模型的方差；

采集终端数据的GMM聚类算法流程为：

201)令l＝0，随机选择初始化的ω^(l),μ^(l),Σ^(l)，求取第i个样本点x_i属于第k类高斯模型的初始化后验概率：

202)将第201)步得到的结果代入最大似然公式计算第l+1次迭代的GMM聚类参数：

式中：N为待聚类的数据点个数；

203)若成立(ε一般取10^-5)，则迭代结束，即为估计得到的GMM聚类参数，继续下一步；否则，回到第202)步；

204)采用贝叶斯概率公式计算第i个样本点x_i属于第k类高斯模型的概率：

根据贝叶斯概率最大准则，将第i个样本点x_i划分到使其概率p(μ_k,Σ_k|x_i)取得最大的那类高斯模型中。

3.根据权利要求1的一种基于高斯混合模型的采集终端生产厂商评价方法，其特征在于，在步骤2)，决策矩阵为：

式中：d_ij表示第i个采集终端供应商第j个指标的值，P为采集终端供应商个数，Q为衡量供应商的采集终端的评价指标个数，其中Q等于8；

效益型指标标准化处理方法为：

成本型指标标准化处理方法为：

式中：和分别表示所有P个采集终端供应商中指标j的最小值和最大值；Ω₁表示效益型指标集合，其中Ω₁＝{1,2,3,8}；Ω₂表示成本型指标集合，其中Ω₂＝{4,5,6,7}。

4.根据权利要求1所述的一种基于高斯混合模型的采集终端生产厂商评价方法，其特征在于，在步骤3)中，分别采用层次分析法、熵权法、相关系数法和变异系数法计算采集终端指标权重，接着用指标的组合优化模型计算组合权重，其中：

层次分析法为：设Q个指标为U＝{u₁,u₂,…,u_Q}，每次取两个指标u_i和u_j，按1～9的比例用a_ij表示指标u_i与u_j对采集终端的影响程度之比，数值越大表示u_i相对于u_j越重要；这样即可获得一个采集终端评估问题的判断矩阵A＝(a_ij)_Q×Q，该判断矩阵具有如下性质：a_ij>0,a_ji＝1/a_ij,(i＝1,2,…,Q；j＝1,2,…,Q)；计算指标权重前首先对判断矩阵A进行一致性检验，如果一致性检验未通过，则说明a_ij取值前后矛盾，需要重新赋值；如果通过，则说明a_ij取值前后一致，判断矩阵A可用于求取指标权重；求取指标权重可以用列和求逆法，即：

将b_j归一化，即可求得归一化后的指标权重为

式中：且

采集终端指标的熵权为：

式中：κ表示常数κ＝1/lnP,并且假定f_ij＝0时,f_ijlnf_ij＝0；且

计算皮尔逊相关系数权重包括：假设d′_x＝(d′_1x,d′_2x,...,d′_Px)^T和d′_y＝(d′_1y,d′_2y,...,d′_Py)^T表示标准化决策矩阵D′的两个列向量，则第x个指标和第y个指标之间的皮尔逊相关系数定义为

式中：和分别是向量d′_x和d′_y的平均值；基于皮尔逊相关系数的权重定义为

式中：且

计算变异系数法权重包括：指标u_j的变异系数定义为：

式中：表示向量(d′_1j,d′_2j,...,d′_Pj)^T的平均值；基于变异系数法的指标权重定义为：

式中：且

指标权重的组合优化模型为：

式中：表示待优化的第j个采集终端指标的组合权重，表示采用第k种方法求解得到的第j个指标的权重，G表示权重确定方法的数目，表示的隶属度函数；显然，越大，权重与最终的组合权重越贴近。

5.根据权利要求1的一种基于高斯混合模型的采集终端生产厂商评价方法，其特征在于，在步骤4)用夹角度量法对采集终端的质量进行综合评估时，包括步骤：

401)形成P个待评估的供应商对应的Q个质量指标的评估矩阵

R＝(r_ij)_P×Q

式中：

402)计算采集终端评估中的理想点和负理想点，其分别为

式中：

403)分别计算每个供应商下的采集终端指标与理想点、负理想点的夹角距离，即

式中：r_i＝(r_i1,r_i2,...,r_iQ)表示评估矩阵R的第i个行向量；

404)计算每个供应商下的采集终端指标与理想点的夹角逼近程度，即

可以看出：第i个供应商生产的采集终端越好，γ_i愈接近于1；反之，质量越差，γ_i越接近于0；因此，可以根据γ_i取值的大小得到采集终端的评估结果。