CN108564136A - 一种基于模糊推理的空域运行态势评估分类方法 - Google Patents
一种基于模糊推理的空域运行态势评估分类方法 Download PDFInfo
- Publication number
- CN108564136A CN108564136A CN201810411484.7A CN201810411484A CN108564136A CN 108564136 A CN108564136 A CN 108564136A CN 201810411484 A CN201810411484 A CN 201810411484A CN 108564136 A CN108564136 A CN 108564136A
- Authority
- CN
- China
- Prior art keywords
- fuzzy
- antibody
- rule
- indicate
- affinity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24317—Piecewise classification, i.e. whereby each classification requires several discriminant rules
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
Abstract
本发明公开了一种基于模糊推理的空域运行态势评估分类方法,属于空域态势评估分类技术领域。包括如下步骤:步骤一:收集待处理扇区的空域运行态势样本;步骤二:基于待处理扇区的空域运行态势样本建立初步模糊推理系统;步骤三:基于多目标种群自适应免疫算法优化模糊推理系统的可解释性与准确性。通过利用本发明所提供的方法,可以针对大规模、高维度的扇区运行数据,围绕空域运行态势评估准确性和可解释性,使用多目标免疫优化算法,优化了空域态势评估的准确性,另外在实现免疫算法时避免了处理高维数据时模糊矩阵规模呈指数增长的情况,大大减少了算法所需的时间复杂度与空间复杂度,提高了收敛精度。
Description
技术领域
本发明属于空域态势评估分类技术领域,具体涉及一种基于模糊推理的空域运行态势评估分类方法。
背景技术
随着我国航空运输业的快速发展,航空业务量与日俱增,航班逐年增加,空域运行态势越发复杂。这些情况使得空中交通管制员工作负荷与航班运行风险不断增加,并由此成为航班延误、管制事故发生的重要原因。
在当前空中交通管理系统中,扇区是管制员对航空器进行指挥的空域基本单元。扇区的空域运行态势复杂程度的高低与空中交通管制员的工作负荷大小有着密切联系。过于复杂的空域态势将提高空中交通管制员错误操作的可能性,造成事故;而较低的复杂度则使得管理系统效率低下,资源浪费。为确保空域运行情况良好,保证空中交通管制员处于适当的工作负荷下,应当及时的对空域结构、飞行流量进行调整。为了实施有效的空域管理举措,空域态势评估成为空管领域中重要的研究课题和亟需解决的问题。
由于扇区的空域运行态势与扇区的数十种动、静态状态特征相关,因此,现有方法普遍使用机器学习、模糊推理系统针对具有多种态势特征的样本,建立分类模型得到总体态势指标。基于大批量数据样本的机器学习方法往往将精准性指标作为首要的评判标准,而对模型的可解释性较为忽略。另一方面,模糊推理系统通过建立知识表达形式和推理机制,使建立的模型具有较明显的物理意义,但现有的具有良好可解释性的模糊推理系统多建立在专家知识的基础上,而基于数据的模糊推理系统往往在精确性有所提升,但可解释性上有所欠缺。
发明内容
本发明的目的在于提供一种基于模糊推理的空域运行态势评估分类方法,在扇区的空域运行态势样本具有大规模、高维度性质的情况下,本发明建立了同时具备可解释性与精确性的分类模型,弥补现有扇区态势评估模型无法兼顾上述两者的不足。
本发明提供的基于模糊推理的空域运行态势评估分类方法,具体包括如下步骤:
步骤一:收集待处理扇区的空域运行态势样本;
获取待处理扇区的空域运行态势样本,形成空域运行态势样本集(简称为样本集),所述样本集中共k条样本,其中每条样本包含待处理扇区在某一单位时间内的n个空域运行态势特征值(简称为特征),并且每条样本被标定有一个态势分类标签(简称标签),表示不同的空域运行态势等级,共有m个不同类的标签。其中,k、n、m均为从1开始的正整数。
所述的特征,是指能够反映空域运行态势的航班航迹分布、空域航路结构、空管运行规则等方面的属性因素,一般用连续或离散数值表示。
步骤二:基于待处理扇区的空域运行态势样本建立初步模糊推理系统;
具体如下:
步骤2.1、哑编码:对m个不同类的标签进行哑编码,即建立m个m维单位正交向量:第j个m维单位正交向量的第j维为1,其余均为0(1≤j≤m)。将进行哑编码后的标签作为以后的标签表达方式;
步骤2.2、模糊化:使用聚类算法(如:Fuzzy C-Means算法)将待处理扇区的空域运行态势样本归并聚类,每一聚类都有一个聚类中心,每一聚类中心可以初始化该聚类每个特征的高斯函数作为每个特征的模糊隶属度函数,初始化该聚类每个标签的钟形函数作为每个标签的模糊隶属度函数,对所述样本的每个特征与标签进行模糊化,将所述样本每个特征与标签的精确值转化为模糊值;其中使用聚类算法可以归并数据以减少后续运算的运算量并将模糊隶属函数的参数初始化。
步骤2.3、模糊规则库:使用每个特征与每个标签的模糊值建立IF-THEN规则,特征的模糊值作为规则前件,标签的模糊值作为规则后件,经过哑编码后的m维单位正交向量转化为模糊值后的向量作为标签模糊向量,标签模糊向量的第j维作为该样本属于第j个标签的置信值;
步骤2.4、模糊推理:确定模糊规则库中的同一规则内部不同维度模糊值间和不同规则输出间的模糊运算符,生成模糊集合;
步骤2.5、解模糊:选用重心法作为解模糊方法,将模糊推理得到的模糊集合解模糊生成预测精确值,将不同维度间的预测精确值组成新的m维向量,若其中第j维分量值最大,则将样本预测为属于第j类;
步骤2.6、反向传播:对于经过上述步骤2.5生成的预测精确值,针对模糊推理系统的分类精确度进行优化,所述分类精确度是预测精确值组成的向量的精确度,包括:建立每个特征的模糊隶属度函数与约束条件,对样本集中的每条样本进行预测分类,计算预测分类与实际分类的误差(交叉熵或均方根误差),所述实际分类是指样本自带的真实分类,将所述误差作为精确度损失函数,判断误差是否达到设定误差阈值,当误差未达到设定误差阈值时,求解误差梯度,并沿误差梯度下降的方向使用反向传播算法,在参数合理范围内更新每个特征的相应模糊隶属度函数与每个标签的相应模糊隶属度函数之参数,以此提高空域态势评估分类的精确度,直到误差达到设定误差阈值。
步骤三:基于多目标种群自适应免疫算法优化模糊推理系统的可解释性与准确性;
对于经过反向传播算法调整过的模糊推理系统,使用多目标种群自适应免疫算法,针对模糊推理系统的可解释性与准确性进行多目标优化,包括:
步骤3.1、抗原识别:将待求解的多目标函数和约束条件作为多目标种群自适应免疫算法的抗原。所述多目标函数包括精确度损失函数与规则库复杂程度评估函数;所述的约束条件是指隶属度函数的参数范围为-1~1。
步骤3.2、抗体初始化:使用步骤二中生成的模糊规则库作为多目标种群自适应免疫算法的抗体,并在该模糊规则库周围随机生成多个模糊规则库作为抗体;为所有模糊规则库中的所有隶属度函数的参数使用实数编码成染色体结构;
步骤3.3、支配区分:对所有抗体进行多目标函数的比较,从中识别所有非支配抗体和支配抗体,所述非支配抗体是指在多目标函数中不存在其他抗体均优于该抗体的抗体,并从非支配抗体中随机取出一个作为标记抗体Abidentified;
步骤3.4、亲和度计算:分别计算标记抗体与非支配抗体的亲和度,标记抗体与支配抗体的亲和度,非支配抗体与支配抗体使用不同的亲和度计算方式;
步骤3.5、免疫选择:选择所有亲和度小于预设亲和度阈值δ的非支配抗体与支配抗体组成被选抗体集合,其余的非支配抗体与支配抗体组成未选抗体集合;
步骤3.6、抗体克隆:预设克隆规模最大值Ncmax,被选抗体集合中的抗体按亲和度高低进行排序克隆,越高的亲和度抗体克隆程度越高;未选抗体集合中的抗体不论亲和度高低则全部进行克隆;
步骤3.7、抗体变异(亲和度成熟):被选抗体集合中的抗体的一维产生变异,未选抗体集合中抗体的两维产生变异,变异程度同亲和度成比例。
步骤3.8、抗体简化:为了提高模糊推理系统的可解释性,在多目标种群自适应免疫算法中加入了简化抗体的步骤用以移除冗余的模糊规则与模糊集合,包括有:移除不重要规则、合并相似规则、移除近似通用模糊集合和合并相似模糊集合;
步骤3.9、抗体再选择:首先挑选非支配抗体;然后将支配抗体按照亲和度从小到大排序后,从亲和度最小的支配抗体开始取;直到新的被选抗体个数与初始化抗体个数相同;再选择完成后,计算两两抗体间距离,如距离大于预设距离阈值λ,则把两个抗体中亲和度大的那个删掉,亲和度小的那个抗体保留;
3.10、种群刷新:判断是否达到约束条件,如果否,重复步骤3.2-3.9直至达到约束条件。
本发明的优点以及带来的有益效果在于:
1、本发明针对大规模、高维度的扇区运行数据,围绕空域运行态势评估准确性和可解释性,使用多目标免疫优化算法,实现了对扇区态势评估的模糊推理系统,这在空域态势评估中是一种全新的方法;
2、本发明通过扇区运行数据建立空域运行态势评估模型,通过设立损失函数对模型中参数进行反向传播算法更新,并在优化算法中仍将其作为目标之一,使得所建立的模型充分利用扇区数据,使得空域态势评估的模糊系统预测精确性大幅度提高;
3、本发明针对初步建立的模糊推理系统模型可解释性,建立了多目标免疫优化算法,同时优化了空域态势评估的准确性,同时优化后的结果,其规则库可以给出具有物理意义、符合人们理解的评估规则;
4、本发明在实现免疫算法时实现了不定长染色体编码与新的距离定义方式、遗传过程种群自适应、将原始的模糊规则库作为抗体,避免了处理高维数据时模糊矩阵规模呈指数增长的情况,大大减少了算法所需的时间复杂度与空间复杂度,提高了收敛精度。
附图说明
图1是本发明基于模糊推理的空域运行态势评估分类方法示例的三步骤框架图;
图2是本发明基于模糊推理的空域运行态势评估分类方法示例步骤2的详细流程示意图;
图3是本发明基于模糊推理的空域运行态势评估分类方法示例步骤3的详细流程示意图。
具体实施方式
下面通过具体实施例并结合附图对本发明的技术方案做进一步的详细描述。
本发明提供的基于模糊推理的空域运行态势评估分类方法,如图1所示,具体包括如下步骤;
步骤一:收集待处理扇区的空域运行态势样本,具体包括:
将基于n个空域运行态势特征值对某一待处理扇区的空域运行态势等级进行计算。采集待处理扇区的空域运行态势样本形成样本集,每条样本中包括n个特征。
空域运行态势特征值是指能够影响或反映空域运行态势的航班航迹分布、空域航路结构、空管运行规则等方面的属性因素,一般用连续或离散数值表示。示例特征如表1所示:
表1 空域运行态势特征集
为找到待处理扇区对应的特征与空域运行态势之间的关联,本发明基于实际空管运行数据采集一定数量的待处理扇区的空域运行态势样本,形成空域运行态势样本集(简称为样本集),所述样本集中共k条样本,一条样本包含着待处理扇区在某一单位时间内的n个空域运行态势特征值,并让空中交通管制员对各条样本对应的空域运行态势等级进行标定,例如,共标定有m种空域运行态势等级,即有m个不同类的标签,这里假设m=3,即共有3种空域运行态势等级:低复杂度态势、中复杂度态势和高复杂度态势,分别简写为L、N和H,得到标定的空域运行态势样本集。其中,k、n、m均为从1开始的正整数。
步骤二:基于待处理扇区的空域运行态势样本建立初步模糊推理系统;
使用步骤一中的实际空管运行数据建立模糊推理系统,并使用反向传播算法进行初步精确度优化,如图2所示;
步骤2.1、哑编码:对上述假设的3个分类标签进行哑编码,即建立3个3维单位正交向量:第j个向量的第j维为1,其余均为0(1≤j≤3),即向量[1,0,0]、[0,1,0]、[0,0,1],并将其作为以后的标签表达方式;
步骤2.2、模糊化:使用聚类算法(如:Fuzzy C-Means算法)将待处理扇区的空域运行态势样本归并聚类成r类,r≥m,每一聚类都有一个聚类中心,每一聚类中心可以初始化该聚类n个特征相应的n个高斯函数作为该聚类相应n个特征的模糊隶属度函数,初始化该聚类m个标签相应的m个钟形函数作为该聚类相应m个标签的模糊隶属度函数,将第p(p=1,2,…,r)聚类所有的高斯函数记为Ap,将第p聚类所有的钟形函数记为Bp,对所述样本的每个特征与标签进行模糊化,将所述样本的每个特征与标签的精确值转化为模糊值,具体地,若现已知第i个样本的第j个标签为1,其余标签为0,则:
第i(i=1,2,…,k)个样本第s(s=1,2,…,n)个特征通过第p聚类第s个高斯函数转化为规则前件的隶属度函数:
第i个样本第j(j=1,2,…,m)个标签通过第p聚类第j个钟形函数转化为规则后件的隶属度函数:
其中,表示第p聚类中第s个高斯函数,表示第i个样本的第s个特征通过函数计算出的隶属度值,表示第i个样本的第s个特征,和分别为函数的参数(表示中心与宽度);表示第p聚类中第j个钟形函数,表示第i个样本的第j个标签通过计算出的隶属度值,表示第i个样本的第j个标签,和分别为函数的参数(表示中心与宽度)。
步骤2.3、模糊规则库:使用样本各维特征与标签的模糊值(隶属度值)建立IF-THEN规则,特征的模糊值作为规则前件,标签的模糊值作为规则后件,经过哑编码后的m维单位正交向量转化为模糊值后的向量作为标签模糊向量,标签模糊向量的第j维作为该样本属于第j个标签的置信值,每一个聚类建立一条规则,因此共有r条规则,建立的模糊规则具有如下形式:
其中,Rp表示第p条规则,xi表示样本,x1…xn表示样本第1到第n个特征,表示本规则下的模糊隶属函数,Cj表示第j类,x∈Cjwith CF=αj表示样本在本规则下属于第j个标签的置信值为αj。
步骤2.4、模糊推理:确定同一规则内部不同维度模糊值间和不同规则输出间的模糊运算符,以便对规则库中的模糊规则进行模糊推理,生成模糊集合,模糊运算符分别采用如下形式:(以样本在第p条规则下为例)
第i个样本规则前件的合并方式:
其中,表示第i个样本规则前件在模糊推理后的隶属度值,为第i个样本的第s个特征在第p条规则的第s个高斯函数下的隶属度值;表示第i个样本的第s个特征,与表示第p个规则中第s个高斯函数的中心与宽度;
第i个样本规则后件的合并方式:
其中,表示第i个样本规则后件在模糊推理后的隶属度值,为第i个样本的第j个标签在第p条规则的第j个钟形函数下的隶属度值,表示第i个样本的第j个标签;
第i个样本通过第p条规则的Mamdani推理生成的模糊集合(将规则前件与规则后件一起推理):
其中,μp(yi)表示最后推理生成的模糊集合。
步骤2.5、解模糊:选用重心法作为解模糊方法,将推理得到的模糊集合解模糊生成预测精确值,将不同维度间的预测精确值组成新的m维向量,若其中第j维分量值最大,则将测试样本的预测为属于第j类,具体表达式如下:
第i个样本生成的模糊集合解模糊生成判别第j个标签的置信值;
其中,gm(xi)表示第i个样本的解模糊后的向量,bf为对应规则钟形函数的中心,yU、yL∈[0,1]为预设常数,其余物理量与上相同。
最终选择的xi所属标签类别为:
步骤2.6、反向传播:对于经过上述步骤2.5生成的预测精确值,针对模糊推理系统的分类精确度进行优化,所述分类精确度是预测精确值组成的向量的精确度,包括:建立每个特征的模糊隶属度函数与约束条件,对样本集进行预测分类,计算预测分类与实际分类的误差(交叉熵或最小均方误差),所述实际分类是指样本自带的真实分类,将误差作为精确度损失函数,判断误差是否达到设定误差阈值,当误差未达到设定误差阈值时,求解误差梯度,并沿误差梯度下降的方向使用反向传播算法,在参数合理范围内更新每个特征的相应模糊隶属度函数与每个标签的相应隶属度函数之参数,以此提高空域态势评估分类的精确度,直到误差达到设定误差阈值。
步骤三:基于多目标种群自适应免疫算法优化模糊推理系统的可解释性与准确性;
为进一步改善初步模糊推理系统的精确性与可解释性,指导空域态势评估分类,使用多目标种群自适应免疫算法同时对初步模糊推理系统精确性与可解释性进行优化,获得其Pareto最优解。如图3所示,具体地,种群自适应免疫算法步骤如下:
步骤3.1、抗原识别:将待求解的多目标函数和约束条件作为免疫算法的抗原,所述多目标函数包括精确度损失函数与规则库复杂程度评估函数。
多目标函数如下(其中精确度损失函数选用交叉熵,均方根误差亦可):
Obj2:Complexity=Nrule+Nset+Rl
其中,h(i)表示样本i的原始所属类别,q(i)表示样本i在模型下的预测所属类别,Nrule表示模糊规则库中规则数目总和,Nset表示模糊集合数目总和,Rl表示每条模糊规则长度总和。
所述的约束条件是指隶属度函数的参数范围-1~1。
步骤3.2、抗体初始化:使用步骤二中生成的模糊规则库作为抗体,并在该模糊规则库周围随机生成多个模糊规则库作为抗体,为模糊规则库中的所有隶属度函数的参数使用实数编码成染色体结构;
步骤3.3、支配区分:从所有抗体中进行多目标函数的比较,从中识别所有非支配抗体和支配抗体,所述非支配抗体是指在多目标函数中不存在其他抗体均优于该抗体的抗体,并从非支配抗体中随机取出一个作为标记的抗体Abidentified;
步骤3.4、亲和度计算:分别计算标记的抗体Abidentified与非支配抗体的亲和度,标记抗体Abidentified与支配抗体的亲和度,非支配抗体与支配抗体将使用不同的亲和度计算方式;亲和度计算公式如下:
定义抗体间距离:
非支配抗体亲和度:
支配抗体亲和度:
Affinityd=dist(Abidentified,Abd)
其中,Abi,Abj表示两个不同抗体,各有k1,k2条规则,Rl表示每条模糊规则长度总和,l表示每条模糊规则内部的各个模糊隶属度函数,表示抗体Abi的第i1个规则,表示抗体Abj的第i2个规则,表示在抗体Abj中与抗体Abi的第i1个规则最接近的规则,为该规则在抗体Abj中的编号,表示在抗体Abi中与抗体Abj的第i2个规则最接近的规则,为该规则在抗体Abi中的编号,Abnd表示非支配抗体,N表示非支配抗体总数,Abd表示支配抗体;
步骤3.5、免疫选择:选择所有亲和度小于预设亲和度阈值δ的非支配抗体与支配抗体组成被选抗体集合,其余的非支配抗体与支配抗体组成未选抗体集合;
步骤3.6、抗体克隆:被选抗体集合中的抗体在预设的克隆规模最大值Ncmax下,按亲和度高低进行排序克隆,越高的亲和度抗体克隆程度越高;未选抗体集合中的抗体不论亲和度高低则全部进行克隆;
步骤3.7、抗体变异(亲和度成熟):被选抗体集合中的抗体的一维产生变异,未选抗体集合中抗体的两维产生变异,变异程度同亲和度成比例,具体如下:
Abnew(i)=Abold(i)+α·N(0,1),i=1,…,n;
其中,Abold(i)与Abnew(i)表示变异前后的抗体,N(0,1)为标准高斯分布,G表示当前代数,Gen表示预设的总代数,rand表示[0,1]间的随机数;Affinity表示抗体亲和度,r表示随代数演变的逐渐变小的比例系数,α表示抗体变异程度与亲和度相关性的比例系数。
步骤3.8、抗体简化:为了提高模糊推理系统的可解释性,在多目标种群自适应免疫算法中加入了简化抗体的步骤用以移除冗余的模糊规则与模糊集合(以下的操作都是在每一个抗体内部操作的,不同抗体之间相互没有影响),具体包括如下:
A、移除不重要规则:在不过度删除规则的情况下可将对模型精确性提高最少的规则删去以提高可解释性:
其中,HAR表示使用所有规则时预测结果的交叉熵,Hγ表示删除第γ条规则后的交叉熵,当如下不等式满足时,不重要的规则将被删除:
其中,cr表示当前模糊规则系统中的规则数,maxr表示模糊规则系统中可以达到的最大规则数,rand表示[0,1]间的随机数,随着迭代次数的改变而改变(下同),pm为第一预设阈值用以控制最少的规则数;
B、合并相似规则:若存在抗体中两个模糊规则相似度符合如下不等式,则可认为这两个模糊规则可用同一种方式进行表示:
其中,表示两个模糊规则系统(抗体)中对应模糊集合的相似度,表示抗体中第θ个规则的第β个高斯函数,表示抗体中第个规则的第β个高斯函数,c与σ表示各自隶属度函数的参数,和表示第θ个规则的第β个高斯函数的中心和宽度,和表示第个规则的第β个高斯函数的中心和宽度,Nrule表示模糊规则库中规则数目总和,Rl表示每条模糊规则长度总和,mr为第二预设阈值;
C、移除近似通用模糊集合:若存在模糊集合与通用模糊集合相似度符合如下不等式,则可认为该模糊集合为近似通用模糊集合且可以删去:
其中,U表示通用模糊集合(本例中为模糊集合的模糊隶属度函数宽度参数的值大于2),ufs为第三预设阈值;
D、合并相似模糊集合:若存在模糊前件的模糊集合或模糊后件的模糊集合符合如下等式,则可以认为两个模糊集合可以共同表示:
其中,和分别表示第θ、个规则的第z个钟形函数,表示模糊后件的模糊集合的相似度,sfs为第四预设阈值;
步骤3.9、抗体再选择:将新生成的抗体与原有抗体混合一起进行重新选择,即首先挑选非支配抗体,然后将支配抗体按照亲和度从小到大排序后,从亲和度最小的支配抗体开始取,直到新的被选抗体个数与初始化抗体个数相同;再选择完成后,计算两两抗体间距离,如距离大于预设距离阈值λ,则把两个抗体中亲和度大的那个删掉,亲和度小的那个抗体保留;
步骤3.10、种群刷新:判断是否达到约束条件?如果否,重复步骤3.2-3.9直至达到约束条件。
本发明提供了一种基于模糊推理的空域运行态势评估分类方法,通过空域态势样本特征的获取,建立反向传播算法调整过的初步模糊推理系统,使用多目标种群自适应免疫算法优化模糊推理系统,最终得到具有优异分类性能,同时兼具解释性规则的空域态势分类评估系统。
本发明实现了将模糊推理系统、反向传播算法、多目标种群自适应免疫优化算法与空域态势样本的有效结合,从而建立了分类准确并有解释含义的空域态势评估分类系统,对于保证空中交通管理系统运行安全性,提升空中交通管理系统运行效率和调控手段实施的精度具有较大意义。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施示例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施示例技术方案的精神和范围。
Claims (8)
1.一种基于模糊推理的空域运行态势评估分类方法,其特征在于,包括如下步骤:
步骤一:收集待处理扇区的空域运行态势样本;
获取待处理扇区的空域运行态势样本,形成空域运行态势样本集,所述样本集中共k条样本,其中每条样本包含待处理扇区在某一单位时间内的n个特征,并且每条样本被标定有一个标签,共有m个不同类的标签;
步骤二:基于待处理扇区的空域运行态势样本建立初步模糊推理系统;
步骤三:基于多目标种群自适应免疫算法优化模糊推理系统的可解释性与准确性。
2.如权利要求1所述的一种基于模糊推理的空域运行态势评估分类方法,其特征在于,所述步骤二具体包括:
步骤2.1、哑编码:对所述m个不同类的标签进行哑编码;
步骤2.2、模糊化:使用聚类算法将待处理扇区的空域运行态势样本归并聚类,每一聚类都有一个聚类中心,每一聚类中心初始化该聚类每个特征的高斯函数作为每个特征的模糊隶属度函数,初始化该聚类每个标签的钟形函数作为每个标签的模糊隶属度函数,对所述样本的每个特征与标签进行模糊化,将所述样本每个特征与标签的精确值转化为模糊值;
步骤2.3、模糊规则库:使用每个特征与每个标签的模糊值建立IF-THEN规则,特征的模糊值作为规则前件,标签的模糊值作为规则后件,经过哑编码后的m维单位正交向量转化为模糊值后的向量作为标签模糊向量;
步骤2.4、模糊推理:确定模糊规则库中的同一规则内部不同维度模糊值间和不同规则输出间的模糊运算符,生成模糊集合;
步骤2.5、解模糊:选用重心法作为解模糊方法,将模糊推理得到的模糊集合解模糊生成预测精确值,将不同维度间的预测精确值组成新的m维向量;
步骤2.6、反向传播:针对预测精确值组成的向量的精确度进行优化。
3.如权利要求2所述的一种基于模糊推理的空域运行态势评估分类方法,其特征在于,所述步骤2.6具体为:建立每个特征的模糊隶属度函数与每个标签的模糊隶属度函数参数取值范围限制,对样本集中的每条样本进行预测分类,计算预测分类与实际分类的误差,所述实际分类是指样本自带的真实分类,将所述误差作为精确度损失函数,判断误差是否达到设定误差阈值,当误差未达到设定误差阈值时,求解误差梯度,并沿误差梯度下降的方向使用反向传播算法,在参数合理范围内更新每个特征的相应模糊隶属度函数与每个标签的相应模糊隶属度函数之参数,直到误差达到设定误差阈值。
4.如权利要求1所述的一种基于模糊推理的空域运行态势评估分类方法,其特征在于,所述步骤三具体包括:
步骤3.1、抗原识别:将待求解的多目标函数和约束条件作为多目标种群自适应免疫算法的抗原;
步骤3.2、抗体初始化:使用步骤二中模糊推理系统的模糊规则库作为抗体,并在该模糊规则库周围随机生成多个模糊规则库作为抗体,为模糊规则库中的所有隶属度函数的参数使用实数编码成染色体结构;
步骤3.3、支配区分:对所有抗体进行多目标函数的比较,从中识别所有非支配抗体和支配抗体,并从非支配抗体中随机取出一个作为标记抗体;
步骤3.4、亲和度计算:分别计算标记抗体与非支配抗体的亲和度,标记抗体与支配抗体的亲和度,非支配抗体与支配抗体使用不同的亲和度计算方式;
步骤3.5、免疫选择:选择亲和度小于预设亲和度阈值的所有非支配抗体与支配抗体组成被选抗体集合,其余的非支配抗体与支配抗体组成未选抗体集合;
步骤3.6、抗体克隆:预设克隆规模最大值,将被选抗体集合中的抗体按亲和度高低排序并进行克隆;未选抗体集合中的抗体不论亲和度高低则全部进行克隆;
步骤3.7、抗体变异:被选抗体集合中的抗体的一维产生变异,未选抗体集合中抗体的两维产生变异,变异程度同亲和度成比例;
步骤3.8、抗体简化:移除冗余的模糊规则与模糊集合,包括有:移除不重要规则、合并相似规则、移除近似通用模糊集合和合并相似模糊集合;
步骤3.9、抗体再选择:首先挑选非支配抗体;然后将支配抗体按照亲和度从小到大排序后,从亲和度最小的支配抗体开始取;直到新的被选抗体个数与初始化抗体个数相同;选择完成后,计算两两抗体间距离,如距离大于预设距离阈值,则把两个抗体中亲和度大的那个删掉,亲和度小的那个抗体保留;
步骤3.10、种群刷新:判断是否达到约束条件,如果否,重复步骤3.2-3.9直至达到约束条件。
5.如权利要求4所述的一种基于模糊推理的空域运行态势评估分类方法,其特征在于,所述的多目标函数如下:
Obj1:
Obj2:Complexity=Nrule+Nset+Rl
其中,h(i)表示样本i的原始所属类别,q(i)表示样本i在模型下的预测所属类别,Nrule表示模糊规则库中规则数目总和,Nset表示模糊集合数目总和,Rl表示每条模糊规则长度总和;
约束条件是指隶属度函数的参数范围为-1~1。
6.如权利要求4所述的一种基于模糊推理的空域运行态势评估分类方法,其特征在于,所述非支配抗体与支配抗体使用不同的亲和度计算方式,亲和度计算公式如下:
定义抗体间距离:
非支配抗体亲和度:
支配抗体亲和度:
Affinityd=dist(Abidentified,Abd)
其中,Abi,Abj表示两个不同抗体,各有k1,k2条规则,Rl表示每条模糊规则长度总和,l表示每条模糊规则内部的各个模糊隶属度函数,表示抗体Abi的第i1个规则,表示抗体Abj的第i2个规则,表示在抗体Abj中与抗体Abi的第i1个规则最接近的规则,为该规则在抗体Abj中的编号,表示在抗体Abi中与抗体Abj的第i2个规则最接近的规则,为该规则在抗体Abi中的编号,nd表示非支配抗体,N表示非支配抗体总数,Abd表示支配抗体。
7.如权利要求4所述的一种基于模糊推理的空域运行态势评估分类方法,其特征在于,所述步骤3.7具体为:
Abnew(i)=Abold(i)+α·N(0,1),i=1,…,n;
其中,Abold(i)与Abnew(i)表示变异前后的抗体,N(0,1)为标准高斯分布,G表示当前代数,Gen表示预设的总代数,rand表示[0,1]间的随机数;Affinity表示抗体亲和度,r表示随代数演变的逐渐变小的比例系数,α表示抗体变异程度与亲和度相关性的比例系数。
8.如权利要求4所述的一种基于模糊推理的空域运行态势评估分类方法,其特征在于,所述步骤3.8具体包括:
A、移除不重要规则:在不过度删除规则的情况下将对模型精确性提高最少的规则删去以提高可解释性:
其中,HAR表示使用所有规则时预测结果的交叉熵,Hγ表示删除第γ条规则后的交叉熵,当如下不等式满足时,不重要的规则将被删除:
其中,cr表示当前模糊规则系统中的规则数,maxr表示模糊规则系统中达到的最大规则数,rand表示[0,1]间的随机数,pm为第一预设阈值;
B、合并相似规则:若存在抗体中两个模糊规则相似度符合如下不等式,则认为这两个模糊规则用同一种方式进行表示:
其中,表示两个模糊规则系统中对应模糊集合的相似度,表示抗体中第θ个规则的第β个高斯函数,表示抗体中第个规则的第β个高斯函数,与σ表示各自隶属度函数的参数,和表示第θ个规则的第β个高斯函数的中心和宽度,和表示第个规则的第β个高斯函数的中心和宽度,Nrule表示模糊规则库中规则数目总和,表示每条模糊规则长度总和,mr为第二预设阈值;
C、移除近似通用模糊集合:若存在模糊集合与通用模糊集合相似度符合如下不等式,则认为该模糊集合为近似通用模糊集合且删去:
其中,U表示通用模糊集合,ufs为第三预设阈值;
D、合并相似模糊集合:若存在模糊前件的模糊集合或模糊后件的模糊集合符合如下等式,则认为两个模糊集合共同表示:
其中,和分别表示第θ个规则、第个规则的第z个钟形函数,表示模糊后件的模糊集合的相似度,sfs为第四预设阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810411484.7A CN108564136B (zh) | 2018-05-02 | 2018-05-02 | 一种基于模糊推理的空域运行态势评估分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810411484.7A CN108564136B (zh) | 2018-05-02 | 2018-05-02 | 一种基于模糊推理的空域运行态势评估分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108564136A true CN108564136A (zh) | 2018-09-21 |
CN108564136B CN108564136B (zh) | 2019-05-03 |
Family
ID=63537744
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810411484.7A Active CN108564136B (zh) | 2018-05-02 | 2018-05-02 | 一种基于模糊推理的空域运行态势评估分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108564136B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109801681A (zh) * | 2018-12-11 | 2019-05-24 | 江苏大学 | 一种基于改进的模糊聚类算法的snp选择方法 |
CN109886352A (zh) * | 2019-03-04 | 2019-06-14 | 北京航空航天大学 | 一种空域复杂度无监督评估方法 |
CN109993225A (zh) * | 2019-03-29 | 2019-07-09 | 北京航空航天大学 | 一种基于无监督学习的空域复杂度分类方法及装置 |
CN110363382A (zh) * | 2019-06-03 | 2019-10-22 | 华东电力试验研究院有限公司 | 全能型乡镇供电所一体化业务融合技术 |
CN111951613A (zh) * | 2020-07-24 | 2020-11-17 | 北京航空航天大学 | 空地协同的广域空域安全态势评估方法 |
CN111985526A (zh) * | 2020-07-02 | 2020-11-24 | 华北理工大学 | 一种基于相似场景聚类的尾随间隔管理策略生成方法及其系统 |
CN113823125A (zh) * | 2021-10-11 | 2021-12-21 | 北京瀚文智远信息技术有限公司 | 一种扇区动态划分方法 |
CN114089632A (zh) * | 2021-11-15 | 2022-02-25 | 陕西师范大学 | 基于模糊逻辑的风湿免疫疾病特征识别方法及系统 |
CN115273564A (zh) * | 2022-09-29 | 2022-11-01 | 北京航空航天大学 | 一种基于多目标优化的空域复杂度调控方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6535860B1 (en) * | 1999-08-02 | 2003-03-18 | National Science Council | Design and hardware synthesis of adaptive weighted fuzzy mean image filter |
CN101106273A (zh) * | 2007-06-04 | 2008-01-16 | 天津大学 | 统一电能质量调节控制方法及装置 |
CN101887531A (zh) * | 2010-06-13 | 2010-11-17 | 北京航空航天大学 | 一种飞行数据知识获取系统及其获取方法 |
CN102663232A (zh) * | 2012-03-13 | 2012-09-12 | 江苏润和软件股份有限公司 | 一种面向用户能效评估的多维仿真分析系统及其方法 |
CN103942623A (zh) * | 2014-04-22 | 2014-07-23 | 中国民航大学 | 一种基于需求与容量不确定性的机场拥挤风险预测方法 |
CN104835178A (zh) * | 2015-02-02 | 2015-08-12 | 郑州轻工业学院 | 一种低信噪比运动小目标的跟踪与识别的方法 |
CN107818342A (zh) * | 2017-10-27 | 2018-03-20 | 重庆邮电大学 | 大数据环境下基于限定模糊规则的多分类系统及方法 |
-
2018
- 2018-05-02 CN CN201810411484.7A patent/CN108564136B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6535860B1 (en) * | 1999-08-02 | 2003-03-18 | National Science Council | Design and hardware synthesis of adaptive weighted fuzzy mean image filter |
CN101106273A (zh) * | 2007-06-04 | 2008-01-16 | 天津大学 | 统一电能质量调节控制方法及装置 |
CN101887531A (zh) * | 2010-06-13 | 2010-11-17 | 北京航空航天大学 | 一种飞行数据知识获取系统及其获取方法 |
CN102663232A (zh) * | 2012-03-13 | 2012-09-12 | 江苏润和软件股份有限公司 | 一种面向用户能效评估的多维仿真分析系统及其方法 |
CN103942623A (zh) * | 2014-04-22 | 2014-07-23 | 中国民航大学 | 一种基于需求与容量不确定性的机场拥挤风险预测方法 |
CN104835178A (zh) * | 2015-02-02 | 2015-08-12 | 郑州轻工业学院 | 一种低信噪比运动小目标的跟踪与识别的方法 |
CN107818342A (zh) * | 2017-10-27 | 2018-03-20 | 重庆邮电大学 | 大数据环境下基于限定模糊规则的多分类系统及方法 |
Non-Patent Citations (6)
Title |
---|
AYSE MERVE ACILAR, AHMET ARSLAN: "A novel approach for designing adaptive fuzzy classifiers based on the combination of an artificial immune network and a memetic algorithm", 《INFORMATION SCIENCES》 * |
HASSANALI MOSALMAN YAZDI 等: "Fuzzy Multi-Objective Genetic Algorithm in Determination of Optimum Mid Connection Location of Off-Centre Bracing System", 《ADVANCED COMPUTER THEORY AND ENGINEERING》 * |
T. Y. HWANG 等: "System Optimal Compensators Placement Via Immune Multi-Objective Algorithm", 《WSEAS TRANSACTIONS ON SYSTEMS》 * |
任杰: "终端区交通态势识别研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
林金星 等: "一种基于分层模糊控制的免疫遗传优化算法", 《东南大学学报(自然科学版)》 * |
钱锋: "10KV配电网规划研究", 《中国优秀博硕士学位论文全文数据库 (硕士) 工程科技Ⅱ辑》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109801681A (zh) * | 2018-12-11 | 2019-05-24 | 江苏大学 | 一种基于改进的模糊聚类算法的snp选择方法 |
CN109886352A (zh) * | 2019-03-04 | 2019-06-14 | 北京航空航天大学 | 一种空域复杂度无监督评估方法 |
CN109993225B (zh) * | 2019-03-29 | 2021-03-23 | 北京航空航天大学 | 一种基于无监督学习的空域复杂度分类方法及装置 |
CN109993225A (zh) * | 2019-03-29 | 2019-07-09 | 北京航空航天大学 | 一种基于无监督学习的空域复杂度分类方法及装置 |
CN110363382A (zh) * | 2019-06-03 | 2019-10-22 | 华东电力试验研究院有限公司 | 全能型乡镇供电所一体化业务融合技术 |
CN111985526A (zh) * | 2020-07-02 | 2020-11-24 | 华北理工大学 | 一种基于相似场景聚类的尾随间隔管理策略生成方法及其系统 |
CN111985526B (zh) * | 2020-07-02 | 2022-03-15 | 华北理工大学 | 一种基于相似场景聚类的尾随间隔管理策略生成方法及其系统 |
CN111951613A (zh) * | 2020-07-24 | 2020-11-17 | 北京航空航天大学 | 空地协同的广域空域安全态势评估方法 |
CN111951613B (zh) * | 2020-07-24 | 2021-08-31 | 北京航空航天大学 | 空地协同的广域空域安全态势评估方法 |
CN113823125A (zh) * | 2021-10-11 | 2021-12-21 | 北京瀚文智远信息技术有限公司 | 一种扇区动态划分方法 |
CN114089632A (zh) * | 2021-11-15 | 2022-02-25 | 陕西师范大学 | 基于模糊逻辑的风湿免疫疾病特征识别方法及系统 |
CN114089632B (zh) * | 2021-11-15 | 2023-08-18 | 陕西师范大学 | 基于模糊逻辑的风湿免疫疾病特征识别方法及系统 |
CN115273564A (zh) * | 2022-09-29 | 2022-11-01 | 北京航空航天大学 | 一种基于多目标优化的空域复杂度调控方法 |
CN115273564B (zh) * | 2022-09-29 | 2022-12-27 | 北京航空航天大学 | 一种基于多目标优化的空域复杂度调控方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108564136B (zh) | 2019-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108564136B (zh) | 一种基于模糊推理的空域运行态势评估分类方法 | |
CN111191732B (zh) | 一种基于全自动学习的目标检测方法 | |
CN108416394B (zh) | 基于卷积神经网络的多目标检测模型构建方法 | |
CN111178611B (zh) | 一种日电量预测的方法 | |
CN109143408B (zh) | 基于mlp的动态区域联合短时降水预报方法 | |
CN110059852A (zh) | 一种基于改进随机森林算法的股票收益率预测方法 | |
CN109902953A (zh) | 一种基于自适应粒子群聚类的电力用户分类方法 | |
CN108764540A (zh) | 基于并行lstm串联dnn的供水管网压力预测方法 | |
CN109993225A (zh) | 一种基于无监督学习的空域复杂度分类方法及装置 | |
CN111179592B (zh) | 基于时空数据流融合分析的城市交通预测方法和系统 | |
Liang et al. | A deep spatio-temporal fuzzy neural network for passenger demand prediction | |
CN110443448A (zh) | 一种基于双向lstm的飞机机位分类预测方法和系统 | |
CN111415049A (zh) | 一种基于神经网络及聚类的停电敏感性分析方法 | |
CN110490234A (zh) | 基于聚类分类联合机制的分类器的构建方法以及分类方法 | |
CN110837857A (zh) | 工业用电负荷预测方法、系统及其存储介质 | |
Chen | Estimating job cycle time in a wafer fabrication factory: A novel and effective approach based on post-classification | |
CN111985152A (zh) | 一种基于二分超球面原型网络的事件分类方法 | |
Eichhorn et al. | Soft computing for automated surface quality analysis of exterior car body panels | |
Gunawan et al. | Long Short-Term Memory Approach for Predicting Air Temperature In Indonesia | |
Fan | Data mining model for predicting the quality level and classification of construction projects | |
CN115393631A (zh) | 基于贝叶斯层图卷积神经网络的高光谱图像分类方法 | |
CN107894970A (zh) | 航站楼出港人数的预测方法和系统 | |
Li et al. | Scheduling strategy of semiconductor production lines with remaining cycle time prediction | |
CN102122349A (zh) | 基于巴氏距离和有向无环图构建多分类支持向量机分类器的方法 | |
CN112365120A (zh) | 一种基于三支决策的经营战略智能生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |