CN113807524A - 量子差分进化算法优化svm的滚动轴承故障诊断方法 - Google Patents
量子差分进化算法优化svm的滚动轴承故障诊断方法 Download PDFInfo
- Publication number
- CN113807524A CN113807524A CN202110923559.1A CN202110923559A CN113807524A CN 113807524 A CN113807524 A CN 113807524A CN 202110923559 A CN202110923559 A CN 202110923559A CN 113807524 A CN113807524 A CN 113807524A
- Authority
- CN
- China
- Prior art keywords
- quantum
- rolling bearing
- differential evolution
- evolution algorithm
- fault
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
- G01M13/045—Acoustic or vibration analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/10—Machine learning using kernel methods, e.g. support vector machines [SVM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明涉及一种量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,将待故障诊断的滚动轴承的输入特征矩阵P输入到所述故障诊断模型中,得到一维列向量Q,Q中出现1代表内圈故障,出现2代表滚动体故障,出现3代表外圈6点钟方向故障,出现4代表外圈3点钟方向故障,出现5代表外圈12点钟方向故障,出现6代表正常。本发明使用量子差分进化算法优化支持向量机模型,从而实现参数的自适应调整。相比传统的差分进化算法,量子的差分进化算法在优化支持向量机模型时,可以减少需要花费的时间,提高优化的精度,使得模型拥有更可靠的故障诊断性能。
Description
技术领域
本发明属于智能故障诊断技术领域,涉及量子差分进化算法优化支持向量机(SVM)的滚动轴承故障诊断方法。
背景技术
随着计算机性能的飞速发展,有越来越多的学者将机器学习中的各种算法应用在滚动轴承的故障诊断中。支持向量机作为机器学习中最为经典的算法之一,在故障诊断领域也有广泛的应用。支持向量机通过训练数据计算出超平面,再使用此超平面来划分数据的类别,从而实现故障诊断。超平面的设置不仅与训练数据相关,还与两个参数C和γ的取值相关,若这两个参数设置的不合适,将会极大的影响支持向量机的诊断精度。因此,使用优化算法对这两个参数进行自适应地调整有相当的必要性。但优化算法也有其局限性,如容易陷入局部最优、鲁棒性不高等问题。特别是在当今的大数据背景下,传统优化算法还需要面对一个不可回避的问题,即优化迭代时间过长,这将会极大的影响故障诊断的效率。
发明内容
为解决现有技术中存在的问题,本发明提供一种量子差分进化算法优化支持向量机(SVM)的滚动轴承故障诊断方法;
为达到上述目的,本发明采用的方案如下:
量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,将由振动传感器采集到的待故障诊断的滚动轴承的输入特征矩阵P输入到故障诊断模型中,得到一维列向量Q;一维列向量Q中的1代表内圈故障,2代表滚动体故障,3代表外圈6点钟方向故障,4代表外圈3点钟方向故障,5代表外圈12点钟方向故障,6代表正常;
故障诊断模型的构建过程为:将振动传感器采集的数据作为训练集对初始状态的支持向量机模型进行训练,通过不断的迭代寻找到精度最高的惩罚因子C和核函数参数γ的取值,得到所述故障诊断模型;
所述迭代的过程中,采用量子差分进化算法对支持向量机模型的惩罚因子C和核函数参数γ进行优化;所述量子差分进化算法是将差分进化算法与量子理论结合得到,具体为:首先,以量子比特的编码方式取代差分进化算法中的实数编码方式,每一个实数值使用m个量子比特来表示;然后,采用改进的量子旋转门来更新量子比特中的量子角度,取代现有的差分进化算法中的变异和交叉策略;
改进的量子旋转门更新后的量子角度与改进的量子旋转门的数学公式为:
其中,θ*代表改进的量子旋转门更新后的量子角度,θ代表改进的量子旋转门更新前的量子角度,Δ代表改进的量子旋转门旋转的角度,rand(0,1)为随机生成的0到1范围内的随机数,G代表量子差分进化算法的最大迭代次数,g代表当前迭代次数。
所述改进的量子旋转门的示意图如图2所示。
量子比特不同于经典比特。经典比特只能为1或者0,而量子比特在被观测之前,处于叠加态之中。处于叠加态的量子比特既是1也是0,即n个量子比特等于2n个经典比特。利用量子比特的特性可以实现运算的加速,提高算法的效率。
差分进化算法是基于种群的优化算法,种群的多样性极大的影响着差分进化算法的性能。理想的情况是:在差分进化算法迭代进化的前期,尽可能的拓展种群的多样性,在差分进化算法迭代进化的后期,尽可能的提高种群的收敛速度,量子差分进化算法同样也是如此。
现有的量子旋转门旋转的方向只有一个,不能根据迭代进化的次数自适应的调整旋转的方向。而本发明提出的改进的量子旋转门,可以在差分进化算法迭代进化的前期,更容易控制量子角度往π/4的方向旋转,对量子比特进行观测之后,量子比特坍塌向|0>的概率与坍塌向|1>的概率会尽可能的均等,保证了种群在前期的多样性,从而增加了量子差分进化算法搜索到最优解的可能性。在差分进化算法迭代进化的后期,更容易控制量子角度往0度或π/2的方向旋转,对量子比特进行观测之后,量子比特坍塌向|0>的概率或坍塌向|1>的概率会提高,保证了种群在后期的收敛性。因此,量子差分进化算法比传统差分进化算法拥有更丰富的种群多样性和更快的收敛速度,从而取得更高的故障诊断精度。
同时,由于量子具有并行性,在相同的迭代次数下,量子差分进化算法花费的时间比传统差分进化算法花费的时间要更短。因此,使用量子差分进化算法优化的支持向量机模型对滚动轴承进行故障诊断,能够取得更高的诊断精度和更快的诊断速度,从而提高模型故障诊断的可靠性以及工程实用性。
作为优选的技术方案:
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,故障诊断模型的构建的具体步骤为:
(1)使用振动传感器采集滚动轴承内圈故障数据、滚动轴承滚动体故障数据、滚动轴承外圈6点钟方向故障数据、滚动轴承外圈3点钟方向故障数据、滚动轴承外圈12点钟方向故障数据和滚动轴承正常状态运行数据;(采集的数据为振幅数据,振幅数据在数值上等于位移的大小,单位是厘米或毫米等长度单位,为标量)并划分训练集和测试集;
(2)构建量子差分进化算法优化支持向量机模型的初始状态的故障诊断模型,初始状态的故障诊断模型中,设置量子差分进化算法的初始种群内个体向量的个数为N,个体向量由支持向量机的两个参数,惩罚因子C和核函数参数γ组成,每一个C和γ都由m个量子比特所构成,量子比特中的量子角度为随机生成的(0,π/2)区间内的任意一个角度值;
(3)并以训练集的数据对初始状态的故障诊断模型进行训练,以最大迭代次数G为终止条件进行迭代,通过不断的优化,寻找到精度最高的惩罚因子C和核函数参数γ的取值,得到所述故障诊断模型,记为:
是关于量子差分进化算法迭代进化G次之后得到的对于支持向量机而言最优的参数组合ZG、输入特征矩阵为P和基于Ptrain和Qtrain建立的支持向量机滚动轴承故障诊断模型,其中,Q为故障诊断模型的输出,Ptrain为滚动轴承故障诊断模型的训练特征矩阵,P为待故障诊断的滚动轴承的输入特征矩阵,Qtrain为滚动轴承故障诊断模型的训练一维列向量;ZG为量子差分进化算法迭代G次之后得到的对于支持向量机模型来说最优的个体向量;
其中,W1为采集的滚动轴承内圈故障数据,W2为采集的滚动轴承滚动体故障数据,W3为采集的滚动轴承外圈6点钟方向故障数据,W4为采集的滚动轴承外圈3点钟方向故障数据,W5为采集的滚动轴承外圈12点钟方向故障数据,W6为采集的滚动轴承正常状态运行数据;
其中,S1为W1的标签向量,S2为W2的标签向量、S3为W3的标签向量、S4为W4的标签向量、S5为W5的标签向量、S6为W6的标签向量;S1、S2、S3、S4、S5、S6为:
其中,1为滚动轴承出现内圈故障的标签,2为滚动轴承出现滚动体故障的标签,3为滚动轴承出现外圈6点钟方向故障的标签,4为滚动轴承出现外圈3点钟方向故障的标签,5为滚动轴承出现外圈12点钟方向故障的标签,6为滚动轴承正常运行的标签。
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,所述量子差分进化算法包括使用量子比特对种群进行编码、使用改进的量子旋转门更新量子角度以及对量子比特进行观测。
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,步骤(3)中,50≤G≤100。
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,步骤(3)中,所述优化的具体步骤为:
(a)生成初始种群;支持向量机模型的惩罚因子C和核函数参数γ作为种群内的个体向量,且使用m个量子比特来表示C与γ的值;
(b)使用量子比特对初始种群进行编码,编码方式为:
sin2θ1,i,C+cos2θ1,i,C=1;
sin2θm,i,C+cos2θm,i,C=1;
sin2θ1,i,γ+cos2θ1,i,γ=1;
sin2θm,i,γ+cos2θm,i,γ=1;
其中,zi,g为量子差分进化算法在第g代种群中的第i个个体向量;g代表当前迭代次数,g∈(1,G);i∈(1,N),N代表种群内个体向量的数目;
Ci,g为第g代种群中第i个个体向量中惩罚因子的值;
γi,g为第g代种群中第i个个体向量中核函数参数的值;
sinθ1,i,C为第i个个体向量中惩罚因子C的第1个量子比特|0>上的概率幅值,cosθ1,i,C为第i个个体向量中惩罚因子C的第1个量子比特|1>上的概率幅值;
sinθm,i,C为第i个个体向量中惩罚因子C的第m个量子比特|0>上的概率幅值,cosθm,i,C第i个个体向量中惩罚因子C的第m个量子比特|1>上的概率幅值;
sinθ1,i,γ为第i个个体向量中核函数参数γ的第1个量子比特|0>上的概率幅值,cosθ1,i,γ为第i个个体向量中核函数参数γ的第1个量子比特|1>上的概率幅值;
sinθm,i,γ为第i个个体向量中核函数参数γ的第m个量子比特|0>上的概率幅值,cosθm,i,γ为第i个个体向量中核函数参数γ的第m个量子比特|1>上的概率幅值;
(c)使用改进的量子旋转门更新量子角度;
(d)每一次迭代之后都对每个量子比特进行观测;所述对每个量子比特进行观测是指:若sin2θ>rand(0,1),则此量子比特坍塌为经典比特1;若sin2θ≤rand(0,1),则此量子比特坍塌为经典比特0;其中,sinθ代表量子比特|0>上的概率幅值;
(e)将观测得到的经典二进制比特转化为十进制数并映射至一定的区间范围内,从而得到N组C和γ的取值。
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,步骤(e)中所述映射的公式如下:
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,所述精度最高的惩罚因子C和核函数参数γ的取值在进化迭代次数达到G次时,以测试集的数据从最后一代种群中的N组C和γ的取值中选出对支持向量机而言精度最高的个体向量,即为C和γ的最优值,记为:
其中,C为支持向量机中的惩罚因子,γ为支持向量机中的核函数参数,G为量子差分进化算法的最大迭代次数,CG,best、γG,best分别为量子差分进化算法迭代G次之后,从第G代种群中选出的对于支持向量机来说精度最高的参数。
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,m=8。
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,N=30。
如上所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,θ∈(0,π/2),Δ=0.01π。
本发明的原理是:
本发明主要由两个算法组成,故障诊断的基础模型为支持向量机,优化算法为量子差分进化算法,并使用此量子差分进化算法优化支持向量机中的两个参数,分别是惩罚因子C和核函数参数γ。由于现有的传统差分进化算法在精度和时间等方面存在一定的局限性,从而影响整个故障诊断模型的性能。因此,本发明将量子理论与差分进化算法相结合,并提出一种改进的量子旋转门,再使用此量子旋转门来更新量子角度。本发明提出的改进的量子旋转门可以控制量子角度在差分进化算法迭代进化的前期更容易往π/4的方向旋转,保证量子差分进化算法在前期的种群多样性,也可以控制量子角度在差分进化算法迭代进化的后期更容易往0度或π/2的方向旋转,保证量子差分进化算法在后期的收敛速度。同时,此量子旋转门也可以根据算法的不同阶段自适应的调整旋转方向,从而有效的提高算法的鲁棒性。
在每一次的迭代进化过程中,都使用此改进的量子旋转门更新种群内所有的量子角度,更新之后的量子角度作为下一代种群,再次使用此改进的量子旋转门更新量子角度,如此反复迭代。同时,在每一次迭代之后,都会对所有的量子比特进行一次观测。将观测之后得到的二进制数转化为十进制数,并映射至提前设置好的区间内,从而得到支持向量机中惩罚因子C和核函数参数γ的取值。最后,经过G次的迭代进化之后,得到最优的惩罚因子C和核函数参数γ的取值,从而完成量子差分进化算法优化支持向量机的滚动轴承故障诊断模型的训练。训练完成的模型便可以实现故障诊断的功能。
有益效果
(1)本发明提出一种量子差分进化算法优化支持向量机的滚动轴承故障诊断模型,相比现有的传统差分进化算法,量子差分进化算法在训练模型时花费的时间要更少,使得本发明提出的模型拥有更高的故障诊断效率;
(2)本发明提出一种改进的量子旋转门来更新量子角度,从而在量子差分进化算法迭代进化的前期增加种群多样性,在量子差分进化算法迭代进化的后期提高收敛速度,最终使得本发明提出的故障诊断模型拥有更高的故障诊断精度。
附图说明
图1为本发明量子差分进化算法的流程图;
图2为改进的量子旋转门的示意图。
具体实施方式
下面结合具体实施方式,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,如图1所示,包括如下过程:
构建故障诊断模型:
(1)使用振动传感器采集滚动轴承内圈故障数据、滚动轴承滚动体故障数据、滚动轴承外圈6点钟方向故障数据、滚动轴承外圈3点钟方向故障数据、滚动轴承外圈12点钟方向故障数据和滚动轴承正常状态运行数据;并划分训练集和测试集;
(2)构建量子差分进化算法优化支持向量机模型的初始状态的故障诊断模型;其中,所述量子差分进化算法是将差分进化算法与量子理论结合得到,具体为:首先,以量子比特的编码方式取代差分进化算法中的实数编码方式,每一个实数值使用m个量子比特来表示(m=8);然后,采用改进的量子旋转门来更新量子比特中的量子角度取代差分进化算法中的变异和交叉策略;其中,初始状态的故障诊断模型中,设置量子差分进化算法的初始种群内个体向量的个数为N,个体向量由支持向量机的两个参数,惩罚因子C和核函数参数γ组成,每一个C和γ都由m个量子比特所构成,量子比特中的量子角度为随机生成的(0,π/2)区间内的任意一个角度值;
(3)以训练集的数据对初始状态的故障诊断模型进行训练,以最大迭代次数G(50≤G≤100)(一般取值为50)为终止条件进行迭代,通过不断的优化,得到N组C和γ的取值;其中,迭代的过程中,采用量子差分进化算法对支持向量机模型的惩罚因子C和核函数参数γ进行优化,优化的具体步骤为:
(a)生成初始种群;支持向量机模型的惩罚因子C和核函数参数γ作为种群内的个体向量,且使用m个量子比特来表示C与γ的值;
(b)使用量子比特对初始种群进行编码,编码方式为:
sin2θ1,i,C+cos2θ1,i,C=1;
sin2θm,i,C+cos2θm,i,C=1;
sin2θ1,i,γ+cos2θ1,i,γ=1;
sin2θm,i,γ+cos2θm,i,γ=1;
其中,zi,g为量子差分进化算法在第g代种群中的第i个个体向量;g代表当前迭代次数,g∈(1,G);i∈(1,N),N代表种群内个体向量的数目;
Ci,g为第g代种群中第i个个体向量中惩罚因子的值;
γi,g为第g代种群中第i个个体向量中核函数参数的值;
sinθ1,i,C为第i个个体向量中惩罚因子C的第1个量子比特|0>上的概率幅值,cosθ1,i,C为第i个个体向量中惩罚因子C的第1个量子比特|1>上的概率幅值;
sinθm,i,C为第i个个体向量中惩罚因子C的第m个量子比特|0>上的概率幅值,cosθm,i,C第i个个体向量中惩罚因子C的第m个量子比特|1>上的概率幅值;
sinθ1,i,γ为第i个个体向量中核函数参数γ的第1个量子比特|0>上的概率幅值,cosθ1,i,γ为第i个个体向量中核函数参数γ的第1个量子比特|1>上的概率幅值;
sinθm,i,γ为第i个个体向量中核函数参数γ的第m个量子比特|0>上的概率幅值,cosθm,i,γ为第i个个体向量中核函数参数γ的第m个量子比特|1>上的概率幅值;
(c)使用改进的量子旋转门更新量子角度;
改进的量子旋转门更新后的量子角度与改进的量子旋转门的数学公式为:
其中,θ*代表改进的量子旋转门更新后的量子角度,θ代表改进的量子旋转门更新前的量子角度,θ∈(0,π/2);Δ代表改进的量子旋转门旋转的角度,Δ=0.01π,rand(0,1)为随机生成的0到1范围内的随机数,G代表量子差分进化算法的最大迭代次数,g代表当前迭代次数。
(d)每一次迭代之后都对每个量子比特进行观测;所述对每个量子比特进行观测是指:若sin2θ>rand(0,1),则此量子比特坍塌为经典比特1;若sin2θ≤rand(0,1),则此量子比特坍塌为经典比特0;其中,sinθ代表量子比特|0>上的概率幅值;
(e)将观测得到的经典二进制比特转化为十进制数并映射至一定的区间范围内,从而得到N组C和γ的取值;所述映射的公式如下:
(4)寻找到精度最高的惩罚因子C和核函数参数γ的取值:在进化迭代次数达到G次时,以测试集的数据从最后一代种群中的N组(N=30)C和γ的取值中选出对支持向量机而言精度最高的个体向量,即为C和γ的最优值,记为:
其中,C为支持向量机中的惩罚因子,γ为支持向量机中的核函数参数,G为量子差分进化算法的最大迭代次数,CG,best、γG,best分别为量子差分进化算法迭代G次之后,从第G代种群中选出的对于支持向量机来说精度最高的参数;
以C和γ的最优值得到所述故障诊断模型,记为:
其中,Q为故障诊断模型的输出,Ptrain为滚动轴承故障诊断模型的训练特征矩阵,P为待故障诊断的滚动轴承的输入特征矩阵,Qtrain为滚动轴承故障诊断模型的训练一维列向量;ZG为量子差分进化算法迭代G次之后得到的对于支持向量机模型来说最优的个体向量;
其中,W1为采集的滚动轴承内圈故障数据,W2为采集的滚动轴承滚动体故障数据,W3为采集的滚动轴承外圈6点钟方向故障数据,W4为采集的滚动轴承外圈3点钟方向故障数据,W5为采集的滚动轴承外圈12点钟方向故障数据,W6为采集的滚动轴承正常状态运行数据;
其中,S1为W1的标签向量,S2为W2的标签向量、S3为W3的标签向量、S4为W4的标签向量、S5为W5的标签向量、S6为W6的标签向量;S1、S2、S3、S4、S5、S6为:
其中,1为滚动轴承出现内圈故障的标签,2为滚动轴承出现滚动体故障的标签,3为滚动轴承出现外圈6点钟方向故障的标签,4为滚动轴承出现外圈3点钟方向故障的标签,5为滚动轴承出现外圈12点钟方向故障的标签,6为滚动轴承正常运行的标签。
依此模型进行故障诊断:
将由振动传感器采集到的待故障诊断的滚动轴承的输入特征矩阵P输入到故障诊断模型中,得到一维列向量Q;Q中出现1代表内圈故障,出现2代表滚动体故障,出现3代表外圈6点钟方向故障,出现4代表外圈3点钟方向故障,出现5代表外圈12点钟方向故障,出现6代表正常。
以凯斯西储大学滚动轴承数据中心的数据集为例。凯斯西储大学滚动轴承试验平台由一个1.5KW的电动机、一个扭矩传感器、一个功率测试计和电子控制器所组成。待测的轴承有两个,其中一个轴承位于驱动端,型号为SKF6025,采样频率为12KHz和48KHz;另一个轴承位于风扇端,型号为SKF6023,采样频率为12KHz。轴承的损伤是由电火花加工而来的单点损伤,损伤的直径有0.1778mm、0.3556mm等,损伤直径越大代表故障越严重。故障类型有三种,分别是轴承内圈、外圈和滚动体故障,其中外圈故障又分为6点钟方向、3点钟方向和12点钟方向的故障。本发明以12KHz采样频率下驱动端轴承的运行数据为例,数据的描述如表1所示:
表1
故障类型 | 标签 |
内圈 | 1 |
滚动体 | 2 |
外圈(6:00) | 3 |
外圈(3:00) | 4 |
外圈(12:00) | 5 |
正常 | 6 |
为了更好的表明本发明所提出的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法的优越性,本发明设置两组对比方法,一种为现有的传统差分进化算法优化支持向量机的滚动轴承故障诊断方法,另外一种为未经优化的支持向量机的滚动轴承故障诊断方法。
将凯斯西储大学滚动轴承数据集的百分之七十用作训练数据,即为Ptrain,并根据表1制作出Qtrain,剩余的百分之三十的数据作为测试数据。再使用上述量子差分进化算法优化的支持向量机模型(模型1)、现有的传统差分进化算法优化支持向量机模型(模型2)和未经优化的支持向量机模型(模型3)进行故障诊断,故障诊断的精度和花费的时间如表2所示:
表2
故障诊断模型 | 故障诊断精度 | 模型训练时间 |
模型1 | 100% | 5.65秒 |
模型2 | 98.87% | 21.77秒 |
模型3 | 97.19% | 0.09秒 |
从表2中可知,本发明提出的故障诊断模型取得的精度是最高的,同时训练所花费的时间相比现有的传统差分进化算法优化支持向量机模型有显著的降低。未经优化的支持向量机模型时间虽然花费的很短,但故障诊断精度是最低的,工程实用性不高。
为了进一步验证模型的鲁棒性,本发明再使用机器学习算法中的标准数据集UCIdataset作为验证数据集进行验证。使用的UCI dataset数据集的介绍如表3所示
表3
数据名称 | 样本个数 | 特征个数 | 标签个数 |
Iris | 150 | 4 | 3 |
Pima | 768 | 8 | 2 |
Glass | 214 | 13 | 6 |
Blood | 748 | 4 | 2 |
将Iris、Pima、Glass和Blood四个数据集分别使用上述三种模型进行分类,每一个数据集的百分之七十作为训练数据,百分之三十作为训练数据,分类精度和花费的时间如表4所示:
表4
从表4中可以看出,本发明提出的量子差分进化算法优化支持向量机的模型在上述四个数据集中取得的精度都是最高的,同时训练模型所花费的时间也要少于现有的传统差分进化优化支持向量机的模型。从而可以得出本发明提出的量子差分进化算法优化支持向量机模型拥有较高的鲁棒性。
Claims (10)
1.量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征是:将由振动传感器采集到的待故障诊断的滚动轴承的输入特征矩阵P输入到故障诊断模型中,得到一维列向量Q;
故障诊断模型的构建过程为:将振动传感器采集的数据作为训练集对初始状态的支持向量机模型进行训练,通过不断的迭代寻找到精度最高的惩罚因子C和核函数参数γ的取值,得到所述故障诊断模型;
所述迭代的过程中,采用量子差分进化算法对支持向量机模型的惩罚因子C和核函数参数γ进行优化;所述量子差分进化算法是将差分进化算法与量子理论结合得到,具体为:首先,以量子比特的编码方式取代差分进化算法中的实数编码方式,每一个实数值使用m个量子比特来表示;然后,采用改进的量子旋转门来更新量子比特中的量子角度取代差分进化算法中的变异和交叉策略;
改进的量子旋转门更新后的量子角度与改进的量子旋转门的数学公式为:
其中,θ*代表改进的量子旋转门更新后的量子角度,θ代表改进的量子旋转门更新前的量子角度,Δ代表改进的量子旋转门旋转的角度,rand(0,1)为随机生成的0到1范围内的随机数,G代表量子差分进化算法的最大迭代次数,g代表当前迭代次数。
2.根据权利要求1所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征在于,故障诊断模型的构建的具体步骤为:
(1)使用振动传感器采集滚动轴承内圈故障数据、滚动轴承滚动体故障数据、滚动轴承外圈6点钟方向故障数据、滚动轴承外圈3点钟方向故障数据、滚动轴承外圈12点钟方向故障数据和滚动轴承正常状态运行数据;并划分训练集和测试集;
(2)构建量子差分进化算法优化支持向量机模型的初始状态的故障诊断模型;
(3)以训练集的数据对初始状态的故障诊断模型进行训练,以最大迭代次数G为终止条件进行迭代,通过不断的优化,寻找到精度最高的惩罚因子C和核函数参数γ的取值,得到所述故障诊断模型,记为:
其中,Q为故障诊断模型的输出,Ptrain为滚动轴承故障诊断模型的训练特征矩阵,P为待故障诊断的滚动轴承的输入特征矩阵,Qtrain为滚动轴承故障诊断模型的训练一维列向量;ZG为量子差分进化算法迭代G次之后得到的对于支持向量机模型来说最优的个体向量;
其中,W1为采集的滚动轴承内圈故障数据,W2为采集的滚动轴承滚动体故障数据,W3为采集的滚动轴承外圈6点钟方向故障数据,W4为采集的滚动轴承外圈3点钟方向故障数据,W5为采集的滚动轴承外圈12点钟方向故障数据,W6为采集的滚动轴承正常状态运行数据;
其中,S1为W1的标签向量,S2为W2的标签向量、S3为W3的标签向量、S4为W4的标签向量、S5为W5的标签向量、S6为W6的标签向量;S1、S2、S3、S4、S5、S6为:
其中,1为滚动轴承出现内圈故障的标签,2为滚动轴承出现滚动体故障的标签,3为滚动轴承出现外圈6点钟方向故障的标签,4为滚动轴承出现外圈3点钟方向故障的标签,5为滚动轴承出现外圈12点钟方向故障的标签,6为滚动轴承正常运行的标签。
3.根据权利要求2所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征在于,所述量子差分进化算法包括使用量子比特对种群进行编码、使用改进的量子旋转门更新量子角度以及对量子比特进行观测。
4.根据权利要求3所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征在于,步骤(3)中,50≤G≤100。
5.根据权利要求4所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征在于,步骤(3)中,所述优化的具体步骤为:
(a)生成初始种群;支持向量机模型的惩罚因子C和核函数参数γ作为种群内的个体向量,且使用m个量子比特来表示C与γ的值;
(b)使用量子比特对初始种群进行编码,编码方式为:
sin2θ1,i,C+cos2θ1,i,C=1;
sin2θm,i,C+cos2θm,i,C=1;
sin2θ1,i,γ+cos2θ1,i,γ=1;
sin2θm,i,γ+cos2θm,i,γ=1;
其中,zi,g为量子差分进化算法在第g代种群中的第i个个体向量;g代表当前迭代次数,g∈(1,G);i∈(1,N),N代表种群内个体向量的数目;
Ci,g为第g代种群中第i个个体向量中惩罚因子的值;
γi,g为第g代种群中第i个个体向量中核函数参数的值;
sinθ1,i,C为第i个个体向量中惩罚因子C的第1个量子比特|0>上的概率幅值,cosθ1,i,C为第i个个体向量中惩罚因子C的第1个量子比特|1>上的概率幅值;
sinθm,i,C为第i个个体向量中惩罚因子C的第m个量子比特|0>上的概率幅值,cosθm,i,C第i个个体向量中惩罚因子C的第m个量子比特|1>上的概率幅值;
sinθ1,i,γ为第i个个体向量中核函数参数γ的第1个量子比特|0>上的概率幅值,cosθ1,i,γ为第i个个体向量中核函数参数γ的第1个量子比特|1>上的概率幅值;
sinθm,i,γ为第i个个体向量中核函数参数γ的第m个量子比特|0>上的概率幅值,cosθm,i,γ为第i个个体向量中核函数参数γ的第m个量子比特|1>上的概率幅值;
(c)使用改进的量子旋转门更新量子角度;
(d)每一次迭代之后都对每个量子比特进行观测;所述对每个量子比特进行观测是指:若sin2θ>rand(0,1),则此量子比特坍塌为经典比特1;若sin2θ≤rand(0,1),则此量子比特坍塌为经典比特0;其中,sinθ代表量子比特|0>上的概率幅值;
(e)将观测得到的经典二进制比特转化为十进制数并映射至一定的区间范围内,从而得到N组C和γ的取值。
8.根据权利要求5所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征在于,m=8。
9.根据权利要求5所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征在于,N=30。
10.根据权利要求1所述的量子差分进化算法优化支持向量机的滚动轴承故障诊断方法,其特征在于,θ∈(0,π/2),Δ=0.01π。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110923559.1A CN113807524B (zh) | 2021-08-12 | 2021-08-12 | 量子差分进化算法优化svm的滚动轴承故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110923559.1A CN113807524B (zh) | 2021-08-12 | 2021-08-12 | 量子差分进化算法优化svm的滚动轴承故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113807524A true CN113807524A (zh) | 2021-12-17 |
CN113807524B CN113807524B (zh) | 2023-07-25 |
Family
ID=78893540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110923559.1A Active CN113807524B (zh) | 2021-08-12 | 2021-08-12 | 量子差分进化算法优化svm的滚动轴承故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113807524B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4254083A1 (de) * | 2022-03-30 | 2023-10-04 | Siemens Aktiengesellschaft | Verfahren zur parametrisierung eines überwachungssystems, parametrisierungsvorrichtung und überwachungssystem |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2562945A1 (en) * | 2011-08-24 | 2013-02-27 | Ciena Corporation | Short-Term Optical Recovery Systems and Methods for Coherent Optical Receivers |
CN106053067A (zh) * | 2016-05-24 | 2016-10-26 | 广东石油化工学院 | 基于量子遗传算法优化支持向量机的轴承故障诊断方法 |
CN110261735A (zh) * | 2019-06-18 | 2019-09-20 | 西华大学 | 基于改进量子布谷鸟算法的配电网故障定位方法 |
US20190340316A1 (en) * | 2018-05-03 | 2019-11-07 | Lam Research Corporation | Predicting etch characteristics in thermal etching and atomic layer etching |
CN110837851A (zh) * | 2019-10-25 | 2020-02-25 | 西安交通大学 | 一种电静液作动器液压泵故障诊断方法 |
CN111323220A (zh) * | 2020-03-02 | 2020-06-23 | 武汉大学 | 风力发电机齿轮箱故障诊断方法及系统 |
CN111597651A (zh) * | 2020-04-30 | 2020-08-28 | 上海工程技术大学 | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 |
CN111598887A (zh) * | 2020-05-25 | 2020-08-28 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种基于lvq-gmm算法和多目标优化分割算法的航天器缺陷检测方法 |
CN112729826A (zh) * | 2020-12-21 | 2021-04-30 | 湘潭大学 | 一种人工鱼群-蛙跳优化极限学习机的轴承故障诊断方法 |
WO2021135630A1 (zh) * | 2019-12-31 | 2021-07-08 | 福州大学 | 基于grcmse与流形学习的滚动轴承故障诊断方法 |
CN113095355A (zh) * | 2021-03-03 | 2021-07-09 | 上海工程技术大学 | 改进差分进化算法优化随机森林的滚动轴承故障诊断方法 |
-
2021
- 2021-08-12 CN CN202110923559.1A patent/CN113807524B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2562945A1 (en) * | 2011-08-24 | 2013-02-27 | Ciena Corporation | Short-Term Optical Recovery Systems and Methods for Coherent Optical Receivers |
CN106053067A (zh) * | 2016-05-24 | 2016-10-26 | 广东石油化工学院 | 基于量子遗传算法优化支持向量机的轴承故障诊断方法 |
US20190340316A1 (en) * | 2018-05-03 | 2019-11-07 | Lam Research Corporation | Predicting etch characteristics in thermal etching and atomic layer etching |
CN110261735A (zh) * | 2019-06-18 | 2019-09-20 | 西华大学 | 基于改进量子布谷鸟算法的配电网故障定位方法 |
CN110837851A (zh) * | 2019-10-25 | 2020-02-25 | 西安交通大学 | 一种电静液作动器液压泵故障诊断方法 |
WO2021135630A1 (zh) * | 2019-12-31 | 2021-07-08 | 福州大学 | 基于grcmse与流形学习的滚动轴承故障诊断方法 |
CN111323220A (zh) * | 2020-03-02 | 2020-06-23 | 武汉大学 | 风力发电机齿轮箱故障诊断方法及系统 |
CN111597651A (zh) * | 2020-04-30 | 2020-08-28 | 上海工程技术大学 | 一种基于hwpso-svdd模型的滚动轴承性能退化评估方法 |
CN111598887A (zh) * | 2020-05-25 | 2020-08-28 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种基于lvq-gmm算法和多目标优化分割算法的航天器缺陷检测方法 |
CN112729826A (zh) * | 2020-12-21 | 2021-04-30 | 湘潭大学 | 一种人工鱼群-蛙跳优化极限学习机的轴承故障诊断方法 |
CN113095355A (zh) * | 2021-03-03 | 2021-07-09 | 上海工程技术大学 | 改进差分进化算法优化随机森林的滚动轴承故障诊断方法 |
Non-Patent Citations (4)
Title |
---|
刘海龙: "基于优化深度置信网络的电机轴承故障诊断方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, no. 6, pages 042 - 304 * |
刘祥 等: "基于EQPSO-GS优化SVM的入侵检测模型研究", 《通信技术》, vol. 53, no. 1, pages 180 - 185 * |
刘芳: "基于改进量子旋转门的量子进化数据聚类", 《电子学报》, vol. 39, no. 9, pages 1 - 6 * |
陆彦希 等: "基于改进CEEMDAN和TEO的轴承故障特征提取方法", 《噪声与振动控制》, vol. 40, no. 2, pages 109 - 114 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4254083A1 (de) * | 2022-03-30 | 2023-10-04 | Siemens Aktiengesellschaft | Verfahren zur parametrisierung eines überwachungssystems, parametrisierungsvorrichtung und überwachungssystem |
Also Published As
Publication number | Publication date |
---|---|
CN113807524B (zh) | 2023-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111323220B (zh) | 风力发电机齿轮箱故障诊断方法及系统 | |
CN108228716B (zh) | 基于加权极限学习机的SMOTE_Bagging集成污水处理故障诊断方法 | |
CN109902399A (zh) | 一种基于att-cnn的变工况下滚动轴承故障识别方法 | |
CN106124212B (zh) | 基于稀疏编码器和支持向量机的滚动轴承故障诊断方法 | |
US10824373B2 (en) | Effective quantum RAM architecture for quantum database | |
CN112883994B (zh) | 一种均衡分布适配的旋转机械变工况故障诊断方法 | |
CN110657984B (zh) | 一种基于强化胶囊网络的行星齿轮箱故障诊断方法 | |
CN112132102B (zh) | 一种深度神经网络结合人工蜂群优化的智能故障诊断方法 | |
CN112163640A (zh) | 基于深度卷积迁移学习行星齿轮箱故障诊断方法及系统 | |
Wang et al. | Intelligent fault diagnosis for planetary gearbox using transferable deep q network under variable conditions with small training data | |
CN114169110B (zh) | 基于特征优选和GWOA-XGBoost的电机轴承故障诊断方法 | |
CN111506862A (zh) | 一种基于多源加权集成迁移学习的滚动轴承故障诊断方法 | |
CN113807524B (zh) | 量子差分进化算法优化svm的滚动轴承故障诊断方法 | |
CN115859077A (zh) | 一种变工况下多特征融合的电机小样本故障诊断方法 | |
Liu et al. | A rotor fault diagnosis method based on BP-Adaboost weighted by non-fuzzy solution coefficients | |
Luo et al. | TRNet: A cross-component few-shot mechanical fault diagnosis | |
Xu et al. | Intelligent fault diagnosis of bearings under small samples: A mechanism-data fusion approach | |
Xing et al. | A novel periodic cyclic sparse network with entire domain adaptation for deep transfer fault diagnosis of rolling bearing | |
CN113869451A (zh) | 一种基于改进jgsa算法的变工况下滚动轴承故障诊断方法 | |
CN113269221A (zh) | 一种基于多尺度卷积神经网络的滚动轴承故障诊断方法 | |
CN117113198A (zh) | 一种基于半监督对比学习的旋转设备小样本故障诊断方法 | |
Yang et al. | Improved graph-regularized deep belief network with sparse features learning for fault diagnosis | |
CN114781448A (zh) | 一种轴承故障特征提取方法、系统、介质以及设备 | |
Ming et al. | rgfc-Forest: An enhanced deep forest method towards small-sample fault diagnosis of electromechanical system | |
CN112444395A (zh) | 基于CMWPE和SaE-ELM的机车轮对轴承故障诊断方法 |
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 |