CN105095516B - 基于谱聚类集成的广播电视用户分群系统及方法 - Google Patents

Abstract

本发明提供一种基于谱聚类集成的广播电视用户分群系统及方法，包括：输入单元，输入收视偏好的参数；节目数据库，存储节目播放信息；收视数据库，采集用户对节目的收看信息；收视偏好空间构建单元，根据输入单元输入的属性特征指标，从节目数据库和收视数据库调取数据源，得到各用户对各节目类型的属性特征指标数据，形成收视偏好矩阵；第一分群单元，基于收视偏好空间对用户进行多次分群；匹配单元，采用一致性函数对分群集合中的簇进行共识匹配，构建簇关系图；第二分群单元，将簇关系图转换为簇关系度矩阵，其作为相似度矩阵，采用谱聚类方法对簇进行分组；集成单元，将数据点在不同簇中出现次数最多的簇所在的组作为该数据点最终所在组。

Description

基于谱聚类集成的广播电视用户分群系统及方法

技术领域

本发明涉及广播电视领域，更为具体地，涉及一种基于谱聚类集成的广播电视用户分群系统及方法。

背景技术

随着广播电视行业的发展以及电视资源的不断丰富，用户在收视过程中形成了不同的节目类型偏好。通过对用户收视行为进行分析，对用户进行分群，找到具有不同收视习惯的用户群体，进行个性节目推荐，成为了广电行业研究焦点。

聚类分析，就是把对象按照性质上的亲疏程度分成多个类或簇，使得类或簇内的数据具有较高相似度，类或簇间的数据具有较高的相异程度，目前聚类分析已成功应用到信息检索、数据挖掘等多个领域。聚类集成利用集成学习技术，通过学习合并数据集的多个聚类结果,得到一个新的聚类结果。与单一聚类算法相比,聚类集成通过处理分布式数据源，不仅可以提高聚类结果的准确性和鲁棒性，而且还可以检测和处理孤立点以及噪音情况。

在聚类集成算法中，传统一致性函数是通过将聚类集体中的多个聚类结果进行结合(或称为集成)，以生成一个统一的聚类结果。传统聚类集成算法主要采用基于互联合矩阵、基于图形划分或基于最大似然估计的一致性函数对聚类成员进行匹配，具有很大的局限性。采用基于传统一致性函数的聚类集成方法对用户进行分群，并不能得到很好的效果。

发明内容

鉴于上述问题，本发明的目的是提供一种能够对广播电视用户进行准确分群的基于谱聚类集成的广播电视用户分群方法及系统。

根据本发明的一个方面，提供一种基于谱聚类的集成广播电视用户分群系统，包括：输入单元，用于输入确定广播电视用户收视偏好所需的参数，其中，所述参数至少包括：节目类型、属性特征指标、空间区域、时间区域、簇个数和组个数，其中，所述属性特征指标是与收视时长成正比的一个或多个广播电视参数；节目数据库，用于将各节目类型的节目播放信息存储为数据源或者从网站上采集各节目的节目标签，将各节目的节目标签与各节目的播放信息存储为数据源，其中，节目标签是对电视节目内容的多角度概括；收视数据库，通过用户机顶盒的收视信号采集用户对各节目类型的节目的收看信息，并存储为数据源；收视偏好空间构建单元，根据输入单元输入的属性特征指标，从节目数据库和收视数据库调取计算节目类型或节目标签的属性特征指标所需的数据源，并进行计算得到各用户对各节目类型或各节目标签的属性特征指标数据，构成收视偏好矩阵，即，收视偏好空间；第一分群单元，基于收视偏好空间对广播电视用户进行多次分群，包括第一相似度空间构建单元，利用高斯函数构建基于收视偏好矩阵的相似度矩阵，即相似度空间；第一特征向量空间构建单元，构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前k个最大的特征值对应的特征向量构成特征向量空间，其中，k为簇个数，为自然数，同时也是特征向量空间的列数；第一谱聚类分群单元，将特征向量空间中每一行作为一个数据点，每一次采用k个数据点作为k个簇的初始聚类中心对广播电视用户进行分群，直到特征向量空间中每一个数据点均作为初始聚类中心进行分群，多次分群结束，多次分群结果构成分群集合，其中，每一次分群结果称为一个谱聚类成员，每一次分群结果形成k个簇；匹配单元，采用基于谱聚类的一致性函数对第一分群单元形成的分群集合中的簇进行共识匹配，找到不同簇之间的对应关系，构建簇关系图；第二分群单元，将匹配单元形成的簇关系图转换为簇关系度矩阵，其作为相似度矩阵，采用谱聚类方法对簇进行分组；集成单元，集成第一分群单元和第二分群单元的分群结果，将数据点在不同簇中出现次数最多的簇所在的组作为该数据点最终所在组，从而确定每一组用户的收视偏好。

根据本发明的另一个方面，提供一种基于谱聚类集成的广播电视用户分群方法，包括：通过输入单元输入确定广播电视用户收视偏好所需的参数，其中，所述参数至少包括：节目类型、属性特征指标、空间区域、时间区域、簇个数和组个数；构建收视偏好空间，具体地，根据输入单元输入的属性特征指标，分别从节目数据库和收视数据库调取上述时间区域和空间区域范围内计算各用户对节目类型的属性特征指标数据所需的数据源，并计算得到各用户对各节目类型或各节目标签的属性特征指标数据，组合形成收视偏好矩阵；基于收视偏好空间对广播电视用户进行多次分群，具体地，采用高斯函数构建基于收视偏好矩阵的相似度矩阵，构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前k个最大的特征值对应的特征向量构成特征向量空间，其中，k为簇个数，为自然数，也是特征向量空间中的列数，将特征向量空间中每一行作为一个数据点，每一行的元素作为数据点，每一次采用k个数据点作为k个簇的初始聚类中心对广播电视用户进行分群，直到特征向量空间中每一个数据点均作为初始聚类中心进行分群，多次分群结束，多次分群结果组成分群集合，其中，每一次分群结果称为一个谱聚类成员，每一次分群结果形成k个簇；构建簇关系图，具体地，采用基于谱聚类的一致性函数对上述分群集合中的簇进行共识匹配，找到不同簇之间的对应关系，构建簇关系图；基于簇关系图采用谱聚类算法对不同分群结果中的簇进行分组，具体地，将簇关系图转换为簇关系度矩阵，其作为相似度矩阵，采用谱聚类方法对簇进行分组；集成基于收视偏好空间对广播电视用户进行多次分群和基于簇关系图对簇进行分组的分群结果，将数据点在不同簇中出现次数最多的簇所在的组作为该数据点最终所在组，从而确定每一组用户的收视偏好。

本发明所述基于谱聚类集成的广播电视用户分群系统及方法通过构建节目类型或节目标签的收视偏好空间，基于所述收视偏好空间采用谱聚类方法将广播电视用户分成多个簇，根据多个簇之间的关系构建簇关系图，基于簇关系图采用谱聚类方法再对簇进行分组，最后集成用户在不同簇中出现的次数作为其所属组，对广播电视用户进行了多次分群，提高了分群的准确性，并且将基于谱聚类的一致性函数应用到广播电视领域，提高了分群的集成性，可以产生更加准确的聚类集成效果。

附图说明

通过参考以下结合附图的说明及权利要求书的内容，并且随着对本发明的更全面理解，本发明的其它目的及结果将更加明白及易于理解。在附图中：

图1是本发明基于谱聚类集成的广播电视用户分群系统；

图2是现有技术电视节目分类体系的示意图；

图3是本发明节目标签的示意图；

图4是本发明基于谱聚类集成的广播电视用户分群方法的流程图；

图5是本发明构建收视偏好空间一个实施例的流程图；

图6是本发明构建收视偏好空间另一个实施例的流程图；

图7是本发明基于收视偏好空间对广播电视用户进行多次分群的一个实施例的流程图；

图8是本发明基于收视偏好空间对广播电视用户进行多次分群的另一个实施例的流程图；

图9是本发明采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群的一个实施例的流程图；

图10是本发明采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群的另一个实施例的流程图；

图11是本发明对基于收视偏好空间对广播电视用户分群形成分群集合构建簇关系图的流程图；

图12是本发明簇关系图的示意图；

图13是本发明基于簇关系图采用谱聚类算法对不同分群结果中的簇进行分组的流程图；

图14是本发明采用k-means算法基于特征向量归一化矩阵对簇进行分组的流程图；

图15是本发明采用聚类方法对簇进行分组的示意图。

在所有附图中相同的标号指示相似或相应的特征或功能。

具体实施方式

在下面的描述中，出于说明的目的，为了提供对一个或多个实施例的全面理解，阐述了许多具体细节。然而，很明显，也可以在没有这些具体细节的情况下实现这些实施例。以下将结合附图对本发明的具体实施例进行详细描述。

以下将结合附图对本发明的具体实施例进行详细描述。

图1是本发明基于谱聚类集成的广播电视用户分群系统，如图1所示，所述广播电视用户分群系统，包括：

输入单元110，用于输入确定广播电视用户收视偏好所需的各种参数，其中，所述参数至少包括：节目类型、属性特征指标、空间区域、时间区域、簇个数和组个数，其中，所述属性特征指标是与收视时长成正比的一个或多个广播电视参数，输入单元110可以是触摸屏、键盘、手写板、鼠标等，其中，所述节目类型包括节目分类体系(图2示出)中的一种或者多种类型，例如，可以为电视剧、电影、新闻、综艺、广告中的一种或多种；所述属性特征指标是与收视时长成正比的一个或多个广播电视参数，例如，收视时长、收视占比、忠诚指数等，其中，收视占比是指某类节目单位播出时间内，用户的收视时长，反映了用户对节目类型的偏爱程度，忠诚指数是指用户对某类节目单位视频次内的收视时长，反映了用户对该节目的忠诚程度；所述空间区域可以是省、市、区等，例如北京市；所述时间区域为选定各节目类型的计算时间段，例如，时间区域为2003年1月1日至2014年1月1日；

节目数据库120，用于将各节目类型的节目播放信息存储为数据源或者从网站上采集各节目的节目标签，将各节目的节目标签与各节目的播放信息存储为数据源，其中，节目标签是对电视节目内容的多角度概括，每一个节目标签体现电视节目一个方面的特性，如图3所示，节目标签包括节目类型标签、年代标签、参加人员标签等，其中节目类型标签可以是新闻类、电视剧类、动漫类等；年代标签可以是古代、现代、抗日年代等；参加人员标签可以是演员、导演、编剧等；

收视数据库130，通过用户机顶盒的收视信号采集用户对各节目类型的节目的收看信息，并存储为数据源；

收视偏好空间构建单元140，根据输入单元110输入的属性特征指标，从节目数据库120和收视数据库130调取计算节目类型或节目标签的属性特征指标所需的数据源，并进行计算得到各用户对各节目类型或各节目标签的属性特征指标数据，组合形成收视偏好矩阵，即，收视偏好空间，具体地，将在图5至6中进行详细描述；

第一分群单元150，基于收视偏好空间构建单元140构建的收视偏好空间对广播电视用户进行多次分群，包括：

第一相似度空间构建单元151，利用高斯函数构建基于收视偏好矩阵的相似度矩阵，即相似度空间；

第一特征向量空间构建单元152，构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前K个最大的特征值对应的特征向量构成特征向量空间，其中，k为簇个数，为自然数，同时也是特征向量空间的列数；

第一谱聚类单元153，将特征向量空间中每一行作为一个数据点，每一次采用k个数据点作为k个簇的初始聚类中心对广播电视用户进行分群，直到特征向量空间中每一个数据点均作为初始聚类中心进行分群，多次分群结束，多次分群结果组成分群集合，其中，每一次分群结果称为一个谱聚类成员，每一次分群结果形成k个簇，具体地，将在图7至10中进行详细描述；

匹配单元160，采用基于谱聚类的一致性函数对第一分群单元150形成的分群集合中的簇进行共识匹配，找到不同簇之间的对应关系，构建簇关系图，具体地，将在图11和12中进行描述；

第二分群单元170，将匹配单元160形成的簇关系图转换为簇关系度矩阵，其作为相似度矩阵，采用谱聚类方法对对簇进行分组，具体地，将在表1以及图13至15的描述中进行详细说明；

集成单元180，集成第一分群单元150和第二分群单元170的分群结果，将数据点在不同簇中出现次数最多的簇所在的组作为该数据点最终所在组，从而确定每一组用户的收视偏好。

优选地，第二分群单元170包括：第二相似度空间构建单元171，将簇关系图转换为簇关系度矩阵，其作为相似度矩阵；第二特征向量空间构造单元172，构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前a个最大的特征值对应的特征向量构成特征向量空间，其中，a为设定值，为输入单元110输入的组个数，为自然数；第二谱聚类单元173，将上述特征向量空间中每一行作为新数据点，以原点作为初始聚类中心，利用K-means方法对簇进行分组。

图4是本发明基于谱聚类的广播电视用户分群方法的流程图，如图3所示，所述广播电视用户分群方法包括：

首先，在步骤S410中，通过输入单元110输入确定广播电视用户收视偏好所需的各种参数，其中，所述参数至少包括：节目类型、属性特征指标、空间区域、时间区域、簇个数和组个数；

确定了收视偏好的参数以后，在步骤S420中，构建收视偏好空间，具体地，根据输入单元110输入的属性特征指标，分别从节目数据库和收视数据库调取上述时间区域和空间区域范围内计算各用户对节目类型或节目标签的属性特征指标数据所需的数据源，并计算得到各用户对各节目类型或各节目标签的属性特征指标数据，组合形成收视偏好矩阵，即，收视偏好空间，具体地，将在图5至6中进行详细描述；

构建了收视偏好空间以后，在步骤S430中，基于收视偏好空间采用谱聚类算法对广播电视用户进行多次分群，具体地，采用高斯函数构建基于收视偏好矩阵的相似度矩阵，将同一用户所有节目的相似度相加构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前k个最大的特征值对应的特征向量构成特征向量空间，其中，k为簇个数，为自然数，也是特征向量空间中的列数，将特征向量空间中每一行作为一个数据点，每一行的元素作为数据点，每一次采用k个数据点作为k个簇的初始聚类中心对广播电视用户进行分群，直到特征向量空间中每一个数据点均作为初始聚类中心进行分群，多次分群结束，多次分群结果组成分群集合，其中，每一次分群结果称为一个谱聚类成员，每一次分群结果形成k个簇，具体地，将在图7至10中进行详细描述。

基于收视偏好空间对广播电视用户进行分群得到分群集合后，在步骤S440中，构建上述分群结果构成的分群集合的簇关系图，具体地，采用基于谱聚类的一致性函数对上述分群集合中的簇进行共识匹配，找到不同簇之间的对应关系，构建簇关系图，具体地，将在图11和12中进行描述；

得到表示分群集合中不同簇之间的对应关系的簇关系图后，在步骤S450中，基于谱关系图采用谱聚类算法对不同分群结果中的簇进行分组，具体地，将簇关系图转换为簇关系度矩阵，其作为相似度矩阵，采用谱聚类方法对簇进行分组，具体地，将在表1以及图13至15的描述中进行详细说明。

在步骤S460中，集成基于收视偏好空间对广播电视用户进行多次分群和基于簇关系图对簇进行分组的分群结果，将数据点在不同簇中出现次数最多的簇所在的组作为该数据点最终所在组，从而确定每一组用户的收视偏好。

图5是本发明构建收视偏好空间一个实施例的流程图，如图5所示，所述构建收视偏好空间的方法包括：

首先，在步骤S510中，确定节目类型的属性特征指标的计算方法，即选取节目类型中一个节目的属性特征指标作为节目类型的属性特征指标或者选取各用户对节目类型中多个甚至所有节目的属性特征指标的平均值作为各用户对该节目类型的属性特征指标，例如，输入单元110输入节目属性特征指标为收视占比，空间区域为北京市，时间区域为2003年1月1日至2014年1月1日，节目类型集合S＝{S₁,...,S_m}，m为节目类型总数，所述时间区域和空间区域内广播电视用户集合U＝{u₁,…,u_n}，n为用户总数，以节目类型S₁＝电视剧为例，电视剧类型的节目集合为S₁＝{sd₁,...,sd_t}，其中，t为电视剧节目总数，可以选取各用户对一个电视剧节目(例如节目sd₁)的收视占比作为各用户对电视剧节目类型的收视占比，也可以选取各用户对多个甚至所有电视剧节目的收视占比的平均值作为各用户电视机节目类型的收视占比。

确定各节目类型的节目属性特征指标的计算方法后，在步骤S520中，从节目数据库和收视数据库调用计算各节目类型的属性特征指标所需的数据源，如上例中，节目属性特征指标为收视占比，节目类型为电视剧，采用各用户对节目sd₁的收视占比作为各用户对电视剧类型的收视占比，则分别从节目数据库和收视数据库调用节目sd₁的播出总时长数据源和收视总时长数据源其中，表示节目sd₁在所述时间区域内的播出总时长，表示用户u_i在所述时间区域内对节目sd₁的收视总时长，1≤i≤n。

从节目数据库和收视数据库调用了各节目类型属性特征指标计算方法所需的数据源后，在步骤S530中，计算各节目类型属性特征指标，例如，采用各用户对节目sd₁的收视占比作为各用户对电视剧类型的收视占比，其中，表示用户u_i对节目sd₁的收视占比，表示用户u_i对节目类型S₁的收视占比，又如，采用各用户对所有电视剧节目的收视占比的平均值作为各用户对电视剧类型的收视占比，其中，表示用户u_i对节目sd_j的收视占比，表示用户u_i在时间区域内收看电视剧节目的收视占比总值；

得到各节目类型的属性特征指标以后，在步骤S540中，将各节目类型的属性特征组合形成收视偏好矩阵，即收视偏好空间P，其中，

，其中，表示用户u_i对节目类型S_j的偏好值，值越大表明用户对该节目类型的兴趣度越高，1≤i≤n，1≤j≤m，

优选地，采用各用户对各节目类型的收视指数作为各用户对各节目类型的的 属性特征指标，其中，收视指数是指用户u_i对 节目类型S_j的收视占比与用户u_i对所有节目类型的收视占比均值的比值做反正切运 算后，乘以圆周率倒数的2倍，数值在在0到1范围内，值越大表明用户对该节目兴趣越高，反 之则越小。

另外，优选地，收视偏好矩阵为各用户对各节目类型的收视指数作为各用户对各节目类型的的属性特征指标组成的矩阵，即

图6是本发明构建收视偏好空间另一个实施例的流程图，如图6所示，所述构建收视偏好空间方法包括：

首先，在步骤S610中，从节目数据库调出各节目类型中节目组成的节目集合，例如，节目类型集合S＝{S₁,...,S_m}，各节目类型中各节目组成的节目集合X＝[x₁...x_Z]，其中，x₁...x_Z为节目标号，z为节目总数。

得到各节目类型中各节目组成的节目集合之后，在步骤S620中，通过爬行技术从网络上采集各节目的特征词作为各节目的节目标签，其中，节目标签是对电视节目内容的多角度概括，每一个节目标签体现电视节目一个方面的特性，节目标签组成的集合为节目标签集合B＝[b₁...b_g]，其中，Tag表示节目标签集合,[tag_z,1,…tag_z,f]表示节目x_z的1至f个节目标签，f表示节目标签的个数，当节目标签不足f个时，为0；为节目标签集合Tag中所有非零和非重叠节目标签的集合，g表示节目x₁至x_z的节目标签的非零和非重叠个数，0<g≤zf。

得到各节目的节目标签组成的节目标签集合后，在步骤S630中，计算各节目标签对于各节目的标注频次，即词频，例如，词频表示标签b_j对于节目x_i标注的频次， 表示采用标签b_j对节目x_i标注的网站个数，表示节目x_i所有标签在不同网站出现的个数之和，1≤i≤z，1≤j≤g。

得到各节目标签对于各节目的标注频次后，在步骤S640中，计算各节目标签对于各节目的逆向节目频率 表示逆向节目频率，表征标签b_j的普遍重要程度，由节目总数z除以包含该标签的节目数量再将得到的商取对数得到。

得到各节目标签对于各节目的标注频次和逆向节目频率后，在步骤S650中，采用词频和逆向节目频率构建节目标签权重矩阵W，即采用词频-逆向节目频率算法(TF-IPF，term frequency–inverse program frequency)构建节目标签权重矩阵W

其中，W表示节目类型标签权重矩阵，为节目类型标签权重矩阵中元素，表示标签b_j对于节目x_i的重要程度。

得到各节目的节目标签权重矩阵W之后，在步骤S660中，通过计算各节目的收视占比或者收视指数构建各节目的收视占比矩阵R或者收视指数矩阵F，

其中，表示用户u_n对节目x_z的收视占比；表示用户u_n对节目x_z的收视指数。

得到各节目的收视占比矩阵R或者收视指数矩阵F后，在步骤S670中，根据各节目的节目标签权重矩阵W以及各节目的收视占比矩阵R或者收视指数矩阵F构建收视偏好矩阵P，

P＝RW或

P＝FW

其中，表示用户对u_n节目标签目b_g的偏好程度，值越大表明用户对该节目的兴趣度越高。

图7是本发明基于收视偏好空间对广播电视用户进行多次分群的一个实施例的流程图，如图7所示，所述基于收视偏好空间对广播电视用户进行多次分群的方法包括：

首先，在步骤S710中，采用高斯函数构建基于收视偏好矩阵的相似度空间A，其中，表示用户u_i 与用户u_j之间的基于节目类型收视偏好的欧氏距离，σ为尺度参数，优选地，采用收视偏好 空间中不同用户之间基于节目标签收视偏好的欧式距离，通过高斯函数构建基于收视偏好 矩阵的相似度，即，表示用户u_i与用户u_j之间的 基于节目标签收视偏好的欧式距离。

采用高斯函数构建了相似度空间以后，在步骤S720中，构建基于相似度矩阵的对角矩阵D，

其中：D表示基于相似度矩阵的对角矩阵,D_ij表示对角矩阵D的元素；

表示矩阵A中用户u_i与其他用户间高斯核相似度的和。

构建了基于相似度矩阵的对角矩阵D后，在步骤S730中，根据对角矩阵D构建拉普拉斯矩阵L。

对相似度矩阵进行了对角变换和拉普拉斯变化以后，在步骤S740中，找出拉普拉斯矩阵L的前k个最大的特征值对应的特征向量V_k构成特征向量空间。

L*V＝λ*V

其中：

L表示拉普拉斯矩阵；V表示矩阵L的特征向量；λ表示矩阵L的特征值；

V_k表示矩阵L的前k个最大特征值对应的特征向量，即特征向量空间，k为特征向量空间的列数也是簇个数，是自然数；

V_ij表示矩阵V_k的元素；

构建了特征向量空间以后，在步骤S750中，对特征向量空间进行归一化处理，得到特征向量归一化矩阵Y

得到特征向量归一化矩阵以后，在步骤S760中，采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群，具体地，将在图9中进行详细描述。

图8是本发明基于收视偏好空间对广播电视用户进行多次分群的另一个实施例的方法流程图，如图8所示，所述对收视偏好空间的广播电视用户进行分群的方法包括：

首先，在步骤S810中，采用高斯函数构建基于收视偏好矩阵的不同高斯核的相似度空间A，

或

其中：

d(u_i，u_j)表示用户u_i与用户u_j之间的基于节目类型收视偏好或者基于节目标签收视偏好的欧式距离；

A表示Ⅰ型、Ⅱ型、Ⅲ型、Ⅳ型高斯核相似度矩阵的统一标识；

A_ij表示矩阵A中的元素，可以采用Ⅰ型、Ⅱ型、Ⅲ型、Ⅳ型高斯核参数予以计算，即代表了A'_ij、A”_ij、A”'_ij、A””_ij；

A'_ij表示Ⅰ型高斯核相似度矩阵中元素，σ'_i表示Ⅰ型高斯核参数，dmax-dmin表示用户间基于节目类型收视偏好或者基于节目标签收视偏好的最大变化范围；

A”_ij表示Ⅱ型高斯核相似度矩阵中元素，σ”_i表示Ⅱ型高斯核参数，dmax表示用户间基于节目类型收视偏好或者基于节目标签收视偏好的最大欧氏距离。

A”'_ij表示Ⅲ型高斯核相似度矩阵中元素，σ”'_i、σ”'_j表示Ⅲ型高斯核参数，h表示设定的用户周围的用户数目，u_h表示用户u_i或u_j周围的第h个用户；

A””_ij表示Ⅳ型高斯核相似度矩阵，σ””_i、σ””_j表示Ⅳ型高斯核参数，u_v表示用户u_i周围的第v个用户，1≤v≤h，u_w表示用户u_j周围的第w个用户，1≤w≤h。

在步骤S820中，对具有不同高斯核的相似度空间进行对角变换和拉普拉斯变换。

在步骤S830中，找到经过对角变换和拉普拉斯变换后的不同高斯核的相似度空间内的前k个最大的特征值对应的特征向量构成具有不同高斯核的多个特征向量空间，并对所述特征向量空间进行归一化处理，形成多个特征向量归一化矩阵，其中，k为每一个特征向量归一化矩阵的列数，也是输入单元输入的簇个数，是自然数；

在步骤S840中，采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群，具体地，将在图10中进行详细描述。

图9是本发明采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群的一个实施例的流程图，如图9所示，所述采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群的方法包括：

首先，在步骤S910中，将特征向量归一化矩阵中的每一行作为一个数据点，随机选取k个数据点，作为k个簇的初始聚类中心。

在步骤S920中，计算数据点与初始聚类中心的欧式距离，并将各数据点指派到与其距离最近的初始聚类中心代表的簇；

在步骤S930中，计算用户分群均方误差E：

其中：

y_i表示数据点；

C_i表示簇；

m_i表示组C_i的聚类中心(数值为簇C_i中数据点的属性均值)；

|y_i-m_i|表示数据点y_i与簇y_i聚类中心m_i的欧式距离。

在步骤S940中，判断数据点所属簇是否发生变化或用户分群均方误差E收敛是否收敛；

如果数据点所属簇发生变化或用户分群均方误差E未收敛，在步骤S950中，取簇内所有数据点的属性均值作为初始聚类中心，返回步骤S920，直到数据点所属簇不发生变化且用户分群均方误差E收敛。

如果数据点所属簇没有发生变化且用户分群均方误差E收敛，则在步骤S960中，采用不同初始聚类中心，即，从特征向量归一化矩阵中选取不同的数据点作为初始聚类中心，重复上述过程完成广播电视用户的多次分群。

在步骤S970中，广播电视用户多次分群的分群结果组合构成分群集合C，

其中：C＝[C¹,…，C^M]表示对用户集体U＝{u₁,…,u_n}进行M次分群构成的分群集 合，C¹,…，C^M为M个谱聚类成员；表示对用户集体U＝{u₁,…,u_n} 的第1次分群结果，其中表示第1次分群结果中属于第k个簇的用户群体，例如，M次分 群中属于第k个簇的用户群体为u₁和u_n，则

图10采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群的另一个实施例的流程图，如图10所示，所述采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群的方法包括：

首先，在步骤S1010中，将基于不同高斯核相似度矩阵形成的四个特征向量归一化矩阵中的一个的每一行作为一个数据点，随机选取k个数据点，作为k个簇的初始聚类中心；

在步骤S1020中，计算该特征向量归一化矩阵中数据点与初始聚类中心的欧式距离，并将各数据点指派到与其距离最近的初始聚类中心代表的簇；

在步骤S1030中，计算用户分群均方误差E：

其中：

y_i表示数据点；

C_i表示簇；

m_i表示组C_i的聚类中心(数值为簇C_i中数据点的属性均值)；

|y_i-m_i|表示数据点y_i与簇y_i聚类中心m_i的欧式距离。

在步骤S1040中，判断数据点所属簇是否发生变化或用户分群均方误差E收敛；

如果数据点所属簇发生变化或用户分群均方误差E未收敛，则，在步骤S1050中，取簇内所有数据点的属性均值作为初始聚类中心，返回步骤S1020，直到数据点所属簇不发生变化且用户分群均方误差E收敛。

如果数据点所属簇没有发生变化且用户分群均方误差E收敛，在步骤S1060中，采用上述方法对其他特征向量归一化矩阵对广播电视用户进行分群，得到四次分群结果；

在步骤S1070中，将四次分群的分群结果组合成分群集合C，

其中：C＝[C¹,…，C⁴]表示对用户集体U＝{u₁,…,u_n}进行4次分群构成的分群集合，C¹,…，C⁴为4个谱聚类成员；表示对用户集体U＝{u₁,…,u_n}的第1次分群结果，其中表示第1次分群结果中属于第k个簇的用户群体。

优选地，上述分群方法还包括：基于不同的高斯核相似度矩阵形成的多个特征向量归一化矩阵并采用不同初始聚类中心对广播电视用户进行多次分群，形成分群集合C，

图11是本发明对基于收视偏好空间对广播电视用户分群形成分群集合构建簇关系图的流程图，如图11所示，所述构建簇关系图的方法为采用基于谱聚类的一致性函数对第一分群单元形成的分群集合中的簇进行共识匹配，找到不同簇之间的对应关系，构建簇关系图，具体地，

首先，在步骤S1110中，将分群集合中的一列(谱聚类成员)作为簇关系图的一行，以每一个谱聚类成员中的簇为每一行中的节点，如图12所示，其中簇i_j代表第i次用户分群结果中第j个用户群体，是簇关系图中第i行第j个节点，例如簇2_3表示第2次分群结果中的第3个用户群体，为簇关系图中第二行的第三个节点，用户群体可以包括用户为u₁、u₃和u_n。

构建了簇关系图的行列和节点之后，在步骤S1120中，以簇之间拥有的相同用户占所有用户的比例作为节点间的边权重，通过与边权重对应的加权边连接所有节点，构建簇关系图，如图12中加权边表示各簇之间的关系，簇1_1与簇2_1、簇2_2、簇3_1之间的连线代表有用户同时归属于这些簇类，加权边的长度与边权重相对应，例如，簇1_3的用户群体为簇1_3的用户群体为u₁、u₂和u₄，簇2_3用户群体为u₁、u₃和u_n，则簇1_3和簇2_3的边权重为1/n，簇3_1的用户群体为u₃、u₅和u_n，则簇1_3和簇2_3的边权重为2/n，簇1_3和簇2_3加权边的长度是簇1_3和簇2_3加权边长度的2倍。

图13是本发明基于簇关系图采用谱聚类算法对不同分群结果中的簇进行分组的流程图，如图13所示，所述基于簇关系图采用谱聚类算法对不同分群结果中的簇进行分组的方法包括：

在步骤S1310中，将簇关系图转换为簇关系度矩阵Ac，其作为相似度矩阵，图10的簇关系图对应的簇关系度矩阵Ac的簇关系度矩阵表如表1所示，

表1

Ac

簇1-1

簇1-2

簇1-3

簇2-1

簇2-2

簇2-3

簇3-1

簇3-2

簇3-3

簇1-1

1

0

0

a<sub>1</sub>

a<sub>2</sub>

a<sub>3</sub>

a<sub>4</sub>

a<sub>5</sub>

a<sub>6</sub>

簇1-2

0

1

0

a<sub>7</sub>

a<sub>8</sub>

a<sub>9</sub>

a<sub>10</sub>

a<sub>11</sub>

a<sub>12</sub>

簇1-3

0

0

1

a<sub>13</sub>

a<sub>14</sub>

a<sub>15</sub>

a<sub>16</sub>

a<sub>17</sub>

a<sub>18</sub>

簇2-1

a<sub>1</sub>

a<sub>7</sub>

a<sub>13</sub>

1

0

0

a<sub>19</sub>

a<sub>20</sub>

a<sub>21</sub>

簇2-2

a<sub>2</sub>

a<sub>8</sub>

a<sub>14</sub>

0

1

0

a<sub>22</sub>

a<sub>23</sub>

a<sub>24</sub>

簇2-3

a<sub>3</sub>

a<sub>9</sub>

a<sub>15</sub>

0

0

1

a<sub>25</sub>

a<sub>29</sub>

a<sub>27</sub>

簇3-1

a<sub>4</sub>

a<sub>10</sub>

a<sub>16</sub>

a<sub>19</sub>

a<sub>22</sub>

a<sub>25</sub>

1

0

0

簇3-2

a<sub>5</sub>

a<sub>11</sub>

a<sub>17</sub>

a<sub>20</sub>

a<sub>23</sub>

a<sub>26</sub>

0

1

0

簇3-3

a<sub>6</sub>

a<sub>12</sub>

a<sub>18</sub>

a<sub>21</sub>

a<sub>24</sub>

a<sub>27</sub>

0

0

1

其中，a₁至a₂₇代表边权重，从表1可以看出簇关系度矩阵Ac为对称矩阵，包含M*k个元素。Ac矩阵对角线上元素为1，其余元素为簇之间的边权重；

在步骤S1320中，构建相似度矩阵的对角矩阵Dc。

其中：

Dc表示相似度矩阵的对角矩阵；

Dc_ii表示对角矩阵Dc的元素，1≤i≤n，1≤j≤n，n表示对角矩阵Dc的行数或列数(数值等于Ac矩阵中的用户群体数目)；

表示矩阵Ac中群体C_i与其他群体间簇关系度的和

在步骤S1330中，将对角矩阵Dc进行拉普拉斯变换，构建拉普拉斯矩阵Lc。

其中：

Lc表示拉普拉斯矩阵；

Ac表示关系度矩阵；

Dc表示对角矩阵；

Lc_ij表示拉普拉斯矩阵Lc的元素，1≤i≤n，1≤j≤n，n表示拉普拉斯矩阵L的行数或列数(数值等于簇对角矩阵Dc的行数或列数)。

在步骤S1340中，找出拉普拉斯矩阵Lc的前a个最大的特征值对应的特征向量Vc_a

Lc*Vc＝λc*Vc

其中：

Lc表示簇Laplacian矩阵；Vc表示矩阵Lc的特征向量；λc表示矩阵Lc的特征值；

Vca表示矩阵Lc的前a个最大特征值对应的特征向量构成的矩阵，即，特征向量空间，a值等于组个数；

Vc_ij表示矩阵Vca的元素，1≤i≤n，1≤j≤a，n、a分别表示矩阵Vca的行数和列数；

在步骤S1350中，对特征向量Vc_a进行行归一化处理，得到特征向量矩阵归一化Yc

其中：

Yc表示特征向量Vca的归一化处理结果；

Yc_ij表示矩阵Yc中的元素；1≤i≤n，1≤j≤a，n、a分别表示矩阵Yc的行数和列数；

表示对特征向量Vca中的每行元素求和后开平方。

在步骤S1360中，采用k-means算法基于特征向量归一化矩阵对簇进行分组，具体地，将在图14中进行详细描述。

图14是本发明采用k-means算法基于特征向量归一化矩阵随簇进行分组的流程图，如图14所示，所述采用k-means算法基于特征向量归一化矩阵随簇进行分组的方法包括：

在步骤S1410中，将特征向量矩阵归一化Yc中元素转化为非负值，具体地，如果矩阵中元素全部为非负值，则不做转化，执行步骤S1420,否则，将矩阵中元素减去数值最小的元素；

在步骤S1420中，将经过非负值转化的特征向量归一化矩阵Yc的每一行作为一个新数据点，计算每一个新数据点到原点的欧式距离，并按照所述欧式距离进行升序或者降序排列，其中，原点为经过非负值转化的归一化特征向量矩阵中元素数值为零的一行，优选地，按照降序排列；

在步骤S1430中，将排好序的新数据点平均分成a组，选择每组的中间那个新数据点作为该组的初始聚类中心点；

在步骤S1440中，计算每个新数据点到每个初始聚类中心的欧式距离，并将每个新数据点分配给离它最近的初始聚类中心所代表的组，形成a个新组；

在步骤S1450中，将每个新组中所有数据点的平均值作为每个新组的聚类中心，计算簇分组均方误差Ec

其中，G_i表示组，1≤i≤a；yc_i表示组G_i中的数据点；m_i表示组G_i的中心，其数值为组G_i中各数据点的均值；|yc_i-m_i|表示数据点yc_i与组G_i聚类中心的欧式距离

在步骤S1460中，判断簇分组均方误差是否收敛或新数据点所属组否变化；

若簇分组均方误差收敛且新数据点所属组无变化，则说明簇分组完成；

若簇分组均方误差不收敛或新数据点所属组变化，则在步骤S1470中，取簇内所有新数据点的属性均值作为初始聚类中心，返回步骤S1440，直到新数据点所属组不发生变化且簇分组均方误差收敛。

图15是本发明采用聚类方法对簇进行分组的示意图，如图15所示，基于偏好矩阵采用谱聚类方法对广播电视用户进行分群形成的分群集合经一致性函数共识匹配后，形成图12所示的簇关系图，采用聚类方法基于簇关系图对簇进行分组后分成3组，如图15所示，第一组中包括簇1_1、簇2_1和簇3_3，第二组包括簇1_2、簇2_3和簇3_2,第三组包括簇1_3、簇2_2和簇3_1，根据数据点(基于收视偏好的特征向量归一化矩阵Y中的一行)在不同簇中出现的系数来选择其所在的组，例如，基于节目类型收视偏好矩阵构建的特征向量归一化矩阵Y中的第一行元素作为一个数据点出现在簇1_1、簇2_1和和簇3_2中，由于簇1_1、簇2_1属于第一组，所以该数据点属于第一组，又如，Y中的第一行元素通过用户u₂对节目类型的偏好值得到，则用户u₂属于第一组，如果该节目类型为电视剧，则第一组的用户对节目类型的偏好为电视剧，同理可以判断用户对节目标签的偏好。

采用本发明所述的基于谱聚类广播电视用户分群系统和方法可以确定用户对节目类型的偏好，也可以确定用户对节目标签的偏好，也可以确定用户对同一节目类型的不同节目标签的偏好，确定用户对节目类型的偏好有利于不同厂商将不同类型的节目提供给对该节目类型偏好的用户，确定用户对节目标签的偏好，有利于同一厂商提供不同节目标签的节目给用户，例如，电影厂可以将电影提供给对电影有偏好的用户，电视台可以将新闻提供给对新闻有偏好的用户，电影厂可以将古装电影提供给对古装电影有偏好的用户，也可以经搞笑电影提供给对搞笑电影有偏好的用户。

综上所述，参照附图以示例的方式描述了根据本发明提出的基于谱聚类集成的广播电视用户分群系统及方法。但是，本领域技术人员应当理解，对于上述本发明所提出的系统及方法，还可以在不脱离本发明内容的基础上做出各种改进。因此，本发明的保护范围应当由所附的权利要求书的内容确定。

Claims (11)

1.一种基于谱聚类集成的广播电视用户分群系统，包括：

输入单元，用于输入确定广播电视用户收视偏好所需的参数，其中，所述参数至少包括：节目类型、属性特征指标、空间区域、时间区域、簇个数和组个数，其中，所述属性特征指标是与收视时长成正比的一个或多个广播电视参数；

节目数据库，用于将各节目类型的节目播放信息存储为数据源或者从网站上采集各节目的节目标签，将各节目的节目标签与各节目的播放信息存储为数据源，其中，节目标签是对电视节目内容的多角度概括；

收视数据库，通过用户机顶盒的收视信号采集用户对各节目类型的节目的收看信息，并存储为数据源；

收视偏好空间构建单元，根据输入单元输入的属性特征指标，从节目数据库和收视数据库调取计算节目类型或节目标签的属性特征指标所需的数据源，并进行计算得到各用户对各节目类型或各节目标签的属性特征指标数据，构成收视偏好矩阵，即，收视偏好空间；

第一分群单元，基于收视偏好空间对广播电视用户进行多次分群，包括第一相似度空间构建单元，利用高斯函数构建基于收视偏好矩阵的相似度矩阵，即相似度空间；第一特征向量空间构建单元，构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前k个最大的特征值对应的特征向量构成特征向量空间，其中，k为簇个数，为自然数，同时也是特征向量空间的列数；第一谱聚类分群单元，将特征向量空间中每一行作为一个数据点，每一次采用k个数据点作为k个簇的初始聚类中心对广播电视用户进行分群，直到特征向量空间中每一个数据点均作为初始聚类中心进行分群，多次分群结束，多次分群结果构成分群集合，其中，每一次分群结果称为一个谱聚类成员，每一次分群结果形成k个簇；

匹配单元，采用基于谱聚类的一致性函数对第一分群单元形成的分群集合中的簇进行共识匹配，找到不同簇之间的对应关系，构建簇关系图；

第二分群单元，将匹配单元形成的簇关系图转换为簇关系度矩阵，其作为相似度矩阵，采用谱聚类方法对簇进行分组；

集成单元，集成第一分群单元和第二分群单元的分群结果，将数据点在不同簇中出现次数最多的簇所在的组作为该数据点最终所在组，从而确定每一组用户的收视偏好。

2.根据权利要求1所述的基于谱聚类集成的广播电视用户分群系统，其中，所述属性特征指标包括收视总时长、收视占比、收视指数和忠诚指数，其中，收视总时长是指空间区域内广播电视用户在时间区域内收看节目类型内节目的总时间；收视占比是指收视总时长与节目在时间区域内播出总时长的比值；收视指数是指某用户对某一节目类型的收视占比与该用户对所有节目类型的收视占比均值的比值做反正切运算后，乘以圆周率倒数的2倍，数值在0到1范围内，值越大表明用户对该节目类型兴趣越高，反之则越小；忠诚指数是指用户对某节目单位视频次内的收视时长，反映了用户对该节目的忠诚程度。

3.根据权利要求1所述的基于谱聚类集成的广播电视用户分群系统，其中，所述第二分群单元包括：

第二相似度空间构建单元，将簇关系图转换为簇关系度矩阵，其作为相似度矩阵；

第二特征向量空间构造单元，构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前a个最大的特征值对应的特征向量构成特征向量空间，其中，a为组个数，是自然数，为设定值；

第二谱聚类单元，将上述特征向量空间中每一行作为新数据点，以原点作为初始聚类中心，利用K-means方法对簇进行分组。

4.一种基于谱聚类集成的广播电视用户分群方法，包括：

通过输入单元输入确定广播电视用户收视偏好所需的参数，其中，所述参数至少包括：节目类型、属性特征指标、空间区域、时间区域、簇个数和组个数；

构建收视偏好空间，具体地，根据输入单元输入的属性特征指标，分别从节目数据库和收视数据库调取上述时间区域和空间区域范围内计算各用户对节目类型的属性特征指标数据所需的数据源，并计算得到各用户对各节目类型或各节目标签的属性特征指标数据，组合形成收视偏好矩阵；

基于收视偏好空间对广播电视用户进行多次分群，具体地，采用高斯函数构建基于收视偏好矩阵的相似度矩阵，构建相似度矩阵的对角矩阵，将对角矩阵进行拉普拉斯变换，得到前k个最大的特征值对应的特征向量构成特征向量空间，其中，k为簇个数，为自然数，也是特征向量空间中的列数，将特征向量空间中每一行作为一个数据点，每一行的元素作为数据点，每一次采用k个数据点作为k个簇的初始聚类中心对广播电视用户进行分群，直到特征向量空间中每一个数据点均作为初始聚类中心进行分群，多次分群结束，多次分群结果组成分群集合，其中，每一次分群结果称为一个谱聚类成员，每一次分群结果形成k个簇；

构建簇关系图，具体地，采用基于谱聚类的一致性函数对上述分群集合中的簇进行共识匹配，找到不同簇之间的对应关系，构建簇关系图；

基于簇关系图采用谱聚类算法对不同分群结果中的簇进行分组，具体地，将簇关系图转换为簇关系度矩阵，其作为相似度矩阵，采用谱聚类方法对簇进行分组；

集成基于收视偏好空间对广播电视用户进行多次分群和基于簇关系图对簇进行分组的分群结果，将数据点在不同簇中出现次数最多的簇所在的组作为该数据点最终所在组，从而确定每一组用户的收视偏好。

5.根据权利要求4所述的基于谱聚类集成的广播电视用户分群方法，其中，所述构建收视偏好空间的方法包括：

确定节目类型的属性特征指标的计算方法，即选取节目类型中一个节目的属性特征指标作为节目类型的属性特征指标或者选取各用户对节目类型中多个甚至所有节目的属性特征指标的平均值作为各用户对该节目类型的属性特征指标；

从节目数据库和收视数据库调用计算各用户对各节目类型的属性特征指标所需的数据源；

计算各用户对各节目类型的属性特征指标，构成收视偏好矩阵P，即收视偏好空间P，

其中，u_i表示广播电视用户，u_i∈U，U为广播电视用户集合，U＝{u₁,…,u_n}，n为用户总数，1≤i≤n；S_j表示节目类型，S_j∈S，S为节目类型集合，S＝{S₁,...,S_m}，m为节目类型总数，1≤j≤m；表示用户u_i对节目类型S_j的偏好值，值越大表明用户对该节目类型的兴趣度越高。

6.根据权利要求4所述的基于谱聚类集成的广播电视用户分群方法，其中，所述构建收视偏好空间的方法包括：

从节目数据库调出各节目类型中节目组成的节目集合X，其中，X＝[x₁...x_Z]，x₁...x_Z为节目标号，z为节目总数；

通过爬行技术从网络上采集各节目的特征词作为各节目的节目标签，构成节目标签集合Tag，

B＝[b₁...b_g]

其中，[tag_z,1,…tag_z,f]表示节目x_z的1至f个节目标签，f表示节目标签的个数，当节目标签不足f个时，为0；为节目标签集合Tag中所有非零和非重叠节目标签的集合，g表示节目x₁至x_z的节目标签的非零和非重叠个数，0<g≤z*f；

计算各节目标签对于各节目的标注频次，即词频其中，词频表示标签b_j对于节目x_i标注的频次， 表示采用标签b_j对节目x_i标注的网站个数，表示节目x_i所有标签在不同网站出现的个数之和，1≤i≤z，1≤j≤g；

计算各节目标签对于各节目的逆向节目频率其中，逆向节目频率表示逆向节目频率，表征标签b_j的普遍重要程度，即，由节目总数z除以包含标签b_j的节目数量再将得到的商取对数得到逆向节目频率

采用词频和逆向节目频率构建节目标签权重矩阵W，

其中，W表示节目标签权重矩阵，为节目标签权重矩阵中元素，表示标签b_j对于节目x_i的重要程度；

通过计算各节目的收视占比或者收视指数构建各节目的收视占比矩阵R或者收视指数矩阵F，

其中，表示用户u_n对节目x_z的收视占比；表示用户u_n对节目x_z的收视指数；

根据各节目的节目标签权重矩阵W以及各节目的收视占比矩阵R或者收视指数矩阵F构建收视偏好矩阵P，

P＝RW或

P＝FW

其中，或表示用户u_n对节目标签目b_g的偏好程度，值越大表明用户对该节目的兴趣度越高。

7.根据权利要求4所述的基于谱聚类集成的广播电视用户分群方法，其中，所述基于收视偏好空间对广播电视用户进行多次分群的方法包括：

步骤S710，采用高斯函数构建基于收视偏好矩阵的相似度空间；

步骤S720和S730，对相似度空间进行对角变换和拉普拉斯变换；

步骤S740和S750，找到经过对角变换和拉普拉斯变换后的相似度空间内的前k个最大的特征值对应的特征向量构成特征向量空间，并对所述特征向量空间进行归一化处理，得到特征向量归一化矩阵，其中，k为特征向量归一化矩阵的列数也是簇个数，是自然数；

步骤S760，采用k-means算法基于特征向量归一化矩阵对广播电视用户进行多次分群，包括：

步骤S910，将特征向量归一化矩阵中的每一行作为一个数据点，随机选取k个数据点，作为k个簇的初始聚类中心；

步骤S920，计算数据点与初始聚类中心的欧式距离，并将各数据点指派到与其距离最近的初始聚类中心代表的簇；

步骤S930，计算用户分群均方误差E：

其中：

y_i表示数据点；

C_i表示簇；

m_i表示组C_i的聚类中心，数值为簇C_i中数据点的属性均值；

|y_i-m_i|表示数据点y_i与簇y_i聚类中心m_i的欧式距离；

步骤S940，判断数据点所属簇是否发生变化或用户分群均方误差E收敛；

如果数据点所属簇发生变化或用户分群均方误差E未收敛，在步骤S950中，则取簇内数据点的属性均值作为初始聚类中心，计算数据点到初始聚类中心的欧式距离，将数据点指派到与其距离最近的初始聚类中心代表的簇，计算用户分群均方误差，直到数据点所属簇不发生变化且用户分群均方误差E收敛；

如果数据点所属簇没有发生变化且用户分群均方误差E收敛，本次广播电视用户分群结束，在步骤S960中，采用不同初始聚类中心，即，从特征向量归一化矩阵中选取不同的数据点作为初始聚类中心，重复步骤S910-S960完成广播电视用户的多次分群；

步骤S970，广播电视用户多次分群的分群结果组合成分群集合C，

其中：C＝[C¹,…，C^M]表示对用户集体U＝{u₁,…,u_n}进行M次分群构成的分群集合，C¹,...，C^M为M个谱聚类成员；表示对用户集体U＝{u₁,...,u_n}的第1次分群结果，其中表示第1次分群结果中属于第k个簇的用户群体。

8.根据权利要求4所述的基于谱聚类集成的广播电视用户分群方法，其中，所述基于收视偏好空间对广播电视用户进行多次分群的方法包括：

采用高斯函数构建基于收视偏好矩阵的不同高斯核的相似度空间A，

或

其中：

d(u_i，u_j)表示用户u_i与用户u_j之间的基于节目类型收视偏好或者基于节目标签收视偏好的欧式距离；

A表示Ⅰ型、Ⅱ型、Ⅲ型、Ⅳ型高斯核相似度矩阵的统一标识；

A_ij表示矩阵A中的元素，采用Ⅰ型、Ⅱ型、Ⅲ型、Ⅳ型高斯核参数予以计算，即代表了A'_ij、A”_ij、A”'_ij、A””_ij；

A'_ij表示Ⅰ型高斯核相似度矩阵中元素，σ'_i表示Ⅰ型高斯核参数，dmax-dmin表示用户间基于节目类型收视偏好或者基于节目标签收视偏好的最大变化范围；

A”_ij表示Ⅱ型高斯核相似度矩阵中元素，σ”_i表示Ⅱ型高斯核参数，dmax表示用户间基于节目类型收视偏好或者基于节目标签收视偏好的最大欧氏距离；

A”'_ij表示Ⅲ型高斯核相似度矩阵中元素，σ”'_i、σ”'_j表示Ⅲ型高斯核参数，h表示设定的用户周围的用户数目，u_h表示用户u_i或u_j周围的第h个用户；

A””_ij表示Ⅳ型高斯核相似度矩阵，σ””_i、σ””_j表示Ⅳ型高斯核参数，u_v表示用户u_i周围的第v个用户，1≤v≤h，u_w表示用户u_j周围的第w个用户，1≤w≤h；

对具有不同高斯核的相似度空间进行对角变换和拉普拉斯变换；

找到经过对角变换和拉普拉斯变换后的不同高斯核的相似度空间内的前k个最大的特征值对应的特征向量构成具有不同高斯核的多个特征向量空间，并对所述特征向量空间进行归一化处理，形成多个特征向量归一化矩阵，其中，k为每一个特征向量归一化矩阵的列数，也是输入单元输入的簇个数，是自然数；

将一个特征向量归一化矩阵中的每一行作为一个数据点，随机选取k个数据点，作为k个簇的初始聚类中心；

计算该特征向量归一化矩阵中数据点与初始聚类中心的欧式距离，并将各数据点指派到与其距离最近的初始聚类中心代表的簇；

计算用户分群均方误差E：

其中：

y_i表示数据点；

C_i表示簇；

m_i表示组C_i的聚类中心，数值为簇C_i中数据点的属性均值；

|y_i-m_i|表示数据点y_i与簇y_i聚类中心m_i的欧式距离；

判断数据点所属簇是否发生变化或用户分群均方误差E收敛；

如果数据点所属簇发生变化或用户分群均方误差E未收敛，则取簇内数据点的属性均值作为初始聚类中心，计算数据点到初始聚类中心的欧式距离，将数据点指派到与其距离最近的初始聚类中心代表的簇，计算用户分群均方误差，直到数据点所属簇不发生变化且用户分群均方误差E收敛；

如果数据点所属簇没有发生变化或用户分群均方误差E收敛，则采用上述方法对其他特征向量归一化矩阵对广播电视用户进行分群，得到四次分群结果；

将四次分群的分群结果组合成分群集合C，

其中：C＝[C¹,...，C⁴]表示对用户集体U＝{u₁,...,u_n}进行4次分群构成的分群集合，C¹,...，C⁴为4个谱聚类成员；表示对用户集体U＝{u₁,...,u_n}的第1次分群结果，其中表示第1次分群结果中属于第k个簇的用户群体。

9.根据权利要求8所述的基于谱聚类集成的广播电视用户分群方法，其中，所述基于收视偏好空间对广播电视用户进行分群的方法还包括：基于不同的高斯核相似度矩阵形成的多个特征向量归一化矩阵并采用不同初始聚类中心对广播电视用户进行多次分群。

10.根据权利要求7或8所述的基于谱聚类集成的广播电视用户分群方法，其中，所述构建簇关系图的方法包括：

以谱聚类成员作为簇关系图的行；

以每一个谱聚类成员中的簇为簇关系图每一行中的节点；

以簇之间拥有的相同用户数占所有用户数的比例作为节点间的边权重；

以加权边连接所有节点，构建簇关系图，其中，加权边表示各簇之间的关系，其长度与边权重对应。

11.根据权利要求7或8所述的基于谱聚类集成的广播电视用户分群方法，其中，基于所述簇关系图对广播电视用户进行分群的方法包括：

将簇关系图转换为簇关系度矩阵，其中，簇关系度矩阵为对称矩阵，包括M*k个元素，M为分群次数，用户分矩阵对角线上元素为1，其余元素为簇之间的边权重；

将簇关系度矩阵作为相似度矩阵，对其进行对角变换和拉普拉斯变换，找到拉普拉斯变换后矩阵的前a个最大的特征值对应的特征向量构成特征向量矩阵，其中，a为特征向量空间的列个数，也是输入单元输入的组个数；

将特征向量矩阵进行归一化处理形成特征向量归一化矩阵；

将特征向量归一化矩阵中的元素转化为非负值，具体地，如果矩阵中元素全部为非负值，则不做转化，否则，将矩阵中元素减去数值最小的元素；

将经过非负值转化的特征向量归一化矩阵的每一行作为一个新数据点，计算每一个新数据点到原点的欧式距离，并按照所述欧式距离进行升序或者降序排列，其中，原点为经过非负值转化的归一化特征向量矩阵中元素数值为零的一行；

将排好序的新数据点平均分成a组，选择每组的中间那个新数据点作为该组的初始聚类中心点；

计算每个新数据点到每个初始聚类中心的欧式距离，并将每个新数据点分配与其距离最近的初始聚类中心所代表的组，形成a个新组；

将每个新组中所有数据点的平均值作为每个新组的聚类中心，计算簇分组均方误差Ec

其中，G_i表示组，1≤i≤a；yc_i表示组G_i中的数据点；m_i表示组G_i的中心，其数值为组G_i中各数据点的均值；|yc_i-m_i|表示数据点yc_i与组G_i聚类中心的欧式距离；

判断簇分组均方误差是否收敛或新数据点所属组否变化；

若簇分组均方误差收敛且新数据点所属组无变化，则说明簇分组完成；

若簇分组均方误差不收敛或新数据点所属组变化，则取簇内新数据点的属性均值作为初始聚类中心，计算新数据点到初始聚类中心的欧式距离，将新数据点指派到与其距离最近的初始聚类中心代表的组，计算簇分组均方误差误差，直到数据点所属组不发生变化且簇分组均方误差收敛。