基于HMM和優(yōu)化的PF的數(shù)控轉(zhuǎn)臺(tái)精度衰退模型

王 剛， 陳 捷,2， 洪榮晶,2， 王 華,2

(1.南京工業(yè)大學(xué) 機(jī)械與動(dòng)力工程學(xué)院，南京 211800；2. 江蘇省工業(yè)裝備數(shù)字制造及控制技術(shù)重點(diǎn)實(shí)驗(yàn)室，南京 211800)

數(shù)控轉(zhuǎn)臺(tái)是高檔數(shù)控機(jī)床的重要功能部件，數(shù)控轉(zhuǎn)臺(tái)的分度精度和重復(fù)定位精度對(duì)數(shù)控機(jī)床的加工精度具有重要影響[1]。在數(shù)控機(jī)床中轉(zhuǎn)臺(tái)處于加工區(qū)域，工作環(huán)境惡劣，容易導(dǎo)致精度的過快衰減。已有的測(cè)量精度方法過于繁瑣，且在很多情況下限于環(huán)境因素難以測(cè)量。因此，找到一種以描述數(shù)控轉(zhuǎn)臺(tái)精度衰減為核心的建模方法不失為一種新的嘗試。目前，建模思路主要分為兩種[2]，一種是根據(jù)數(shù)控轉(zhuǎn)臺(tái)的物理結(jié)構(gòu)建立數(shù)學(xué)模型；Yu等[3]分析了一種大型數(shù)控轉(zhuǎn)臺(tái)的電液耦合的動(dòng)態(tài)特性，并且認(rèn)為數(shù)控轉(zhuǎn)臺(tái)蝸輪傳動(dòng)系統(tǒng)的綜合反向間隙是影響分度精度和動(dòng)態(tài)特性的關(guān)鍵因素之一；然而，建立一個(gè)能準(zhǔn)確且全面地反映數(shù)控轉(zhuǎn)臺(tái)的精度變化的數(shù)學(xué)模型十分困難；另一種是通過智能算法在振動(dòng)信號(hào)等數(shù)據(jù)驅(qū)動(dòng)的基礎(chǔ)上建立模型[4]，與數(shù)學(xué)模型相比智能算法不需要分解數(shù)控轉(zhuǎn)臺(tái)的結(jié)構(gòu)，而是通過振動(dòng)信號(hào)等數(shù)據(jù)作為觀察值，對(duì)蘊(yùn)含在數(shù)據(jù)中的反應(yīng)數(shù)控轉(zhuǎn)臺(tái)狀態(tài)變化的信息進(jìn)行發(fā)掘，從而轉(zhuǎn)化為對(duì)數(shù)控轉(zhuǎn)臺(tái)狀態(tài)的評(píng)價(jià)。Benkedjouh等[5]使用小波變換進(jìn)行降噪和提取特征向量，運(yùn)用支持向量機(jī)(Support Vector Machine， SVM)建立剩余壽命預(yù)測(cè)模型。Lu等[6]在大型回轉(zhuǎn)支承的研究中使用小波能量熵(Wavelet Energy Entropy， WEE)和檢驗(yàn)?zāi)B(tài)分解法(Empirical Mode Decomposition， EMD)進(jìn)行降噪和特征向量提取，使用主成分分析法(Principal Component Analysis， PCA)對(duì)多維特征向量進(jìn)行降維，最后運(yùn)用最小二乘支持向量機(jī)建立了衰退預(yù)測(cè)模型。然而，數(shù)控轉(zhuǎn)臺(tái)精度變化的趨勢(shì)信息在轉(zhuǎn)臺(tái)的振動(dòng)信號(hào)中往往非常微弱，很容易被周圍環(huán)境產(chǎn)生的白噪聲掩蓋。Matej等[7]在研究回轉(zhuǎn)支承時(shí)運(yùn)用EEMD-MSPCA(Ensemble Empirical Mode Decomposition Multi-scale Principal Component Analysis)進(jìn)行降噪并提取特征向量，為解決信號(hào)微弱問題提供了思路。在此基礎(chǔ)上Feng等[8]進(jìn)行了改進(jìn)，建立了EEMD-PCA算法，克服了其在降噪時(shí)經(jīng)驗(yàn)參數(shù)難以獲得的缺點(diǎn)。除此之外，數(shù)控轉(zhuǎn)臺(tái)的驅(qū)動(dòng)形式、摩擦特性、系統(tǒng)剛度、阻尼因素、消隙預(yù)載等眾多因素都會(huì)影響到其最終的動(dòng)靜態(tài)分度精度以及其精度保持性[9]，其中某些因素難以直接觀察，因此需要一些更具有針對(duì)性的算法。Wang等[10]在研究中將隱馬爾科夫(Hidden Markov Model, HMM)方法引入到齒輪箱的壽命預(yù)測(cè)中，比較有效地識(shí)別出某些難以觀察到的狀態(tài)。張西寧等[11]將主分量分析法(PCA)與隱馬爾科夫(HMM)結(jié)合，用到對(duì)軸承的診斷檢測(cè)中。

針對(duì)以上的問題，結(jié)合數(shù)控轉(zhuǎn)臺(tái)串聯(lián)系統(tǒng)的結(jié)構(gòu)特點(diǎn)，利用隱馬爾科夫(HMM)算法既可以討論可觀察狀態(tài)又可以討論隱藏狀態(tài)的特點(diǎn)，提出用振動(dòng)信號(hào)訓(xùn)練過的隱馬爾科夫(HMM)模型進(jìn)行精度衰退診斷并建立健康狀態(tài)指標(biāo)，以粒子群算法(Particle Swarm Optimization, PSO)優(yōu)化過的粒子濾波(Particle Filtering, PF)算法作為預(yù)測(cè)模型計(jì)算數(shù)控轉(zhuǎn)臺(tái)精度衰退的剩余壽命，并將優(yōu)化前后的結(jié)果以及其他優(yōu)化方法的結(jié)果進(jìn)行對(duì)比分析。同時(shí)，在數(shù)控轉(zhuǎn)臺(tái)精度衰退加速壽命試驗(yàn)過程中使用激光干涉儀定期測(cè)量轉(zhuǎn)臺(tái)的重復(fù)定位精度，以驗(yàn)證該模型的有效性。

1 構(gòu)建數(shù)控轉(zhuǎn)臺(tái)精度衰退預(yù)測(cè)模型

1.1 隱馬爾科夫(HMM)算法

隱馬爾科夫算法是建立在馬爾科夫假設(shè)、不動(dòng)性假設(shè)和輸出獨(dú)立性假設(shè)的基礎(chǔ)上。其核心思想是通過可觀察到的狀態(tài)來推斷出隱藏狀態(tài)。對(duì)HMM有下面五個(gè)參數(shù)來表述[12]：馬爾科夫鏈的狀態(tài)數(shù)N；每種狀態(tài)對(duì)應(yīng)的觀察值的數(shù)目M；狀態(tài)轉(zhuǎn)移概率矩A={αij}；在狀態(tài)j時(shí)的觀察概率分布B={βj(k)}；初始狀態(tài)概率分布v={vi}。HMM的表達(dá)式可以簡(jiǎn)寫成λ=(A,B,v)。

在使用HMM建立觀察序列的模型時(shí)通常使用三種算法來解決三個(gè)基本的問題[13]：前向-后向算法解決評(píng)估問題，即已知觀察序列和概率模型，計(jì)算觀察序列的概率值P(O|λ)。Viterbi算法解決解碼問題，即已知觀察序列和概率模型，找到一個(gè)最可能的隱藏狀態(tài)序列S。Baum-Welch算法解決學(xué)習(xí)問題，即通過已知的觀察序列建立最有可能的HMM模型。

本文使用對(duì)精度變化敏感的若干統(tǒng)計(jì)特征量組成觀察矩陣OB作為觀察值序列，人為設(shè)置了四個(gè)隱藏狀態(tài)，訓(xùn)練一個(gè)HMM模型。再使用訓(xùn)練過的模型進(jìn)行精度衰退的早期診斷。與以往使用HMM算法進(jìn)行狀態(tài)識(shí)別不同的是，本文在進(jìn)行早期診斷的過程中，使用HMM算法將數(shù)控轉(zhuǎn)臺(tái)的正常精度狀態(tài)與后面變化的精度狀態(tài)作比較，得出的是各個(gè)精度狀態(tài)下的精度與正常狀態(tài)下的精度的匹配程度，由此建立數(shù)控轉(zhuǎn)臺(tái)精度衰退的健康狀態(tài)指標(biāo)H。

1.2 粒子濾波算法(PF)

(1)

在建立粒子濾波算法預(yù)測(cè)數(shù)控轉(zhuǎn)臺(tái)精度變化的過程中，首先要建立其狀態(tài)空間模型方程。在已有的健康指標(biāo)H的基礎(chǔ)上，使用高斯函數(shù)對(duì)H進(jìn)行非線性回歸擬合，得到擬合的狀態(tài)空間模型[15]

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

1.3 粒子群(PSO)算法優(yōu)化粒子濾波算法

粒子濾波算法雖然在理論上可以表示任意形式的概率分布，但是依然存在一系列的問題[16]，上文中的三個(gè)狀態(tài)空間的參數(shù)ak,bk,ck在算法預(yù)測(cè)前就需要確定其初始值，現(xiàn)有的文獻(xiàn)中確定這個(gè)初始值時(shí)雖然可以利用粒子濾波的狀態(tài)跟蹤能力通過觀察數(shù)據(jù)跟蹤獲得，但獲得的最優(yōu)參數(shù)只是一個(gè)范圍，具體值的確定依然需要依靠經(jīng)驗(yàn)。本文提出使用粒子群(PSO)算法找出這三個(gè)參數(shù)的最優(yōu)值，達(dá)到優(yōu)化粒子濾波算法的預(yù)測(cè)能力的目的。

PSO可以描述為：對(duì)d維空間中隨機(jī)初始化的一個(gè)粒子群的速度和位置進(jìn)行迭代，每次迭代之后，計(jì)算其自身搜索到的歷史最優(yōu)解pbest和全部粒子的歷史最優(yōu)解gbest，并由適應(yīng)度函數(shù)對(duì)每個(gè)粒子進(jìn)行優(yōu)化評(píng)價(jià)。迭代過程中每個(gè)粒子的速度和位置公式為

v=wi×vp+c1×r1×(pbest-x)+c2×r2×(gbest-x)

(10)

x=x+vp

(11)

式中：wi為慣性權(quán)重，本文取為0.7[17]；c1為粒子自身路徑中的歷史最值權(quán)重，一般取2；c2為群體最優(yōu)值權(quán)重，一般取2；r1,r2為[0,1]內(nèi)隨機(jī)數(shù)。

本文中的適應(yīng)度函數(shù)為

(12)

式中：y為觀測(cè)值；yp為預(yù)測(cè)值；在經(jīng)過300次迭代之后，粒子群算法停止優(yōu)化。具體的優(yōu)化過程如表1所示。

表1 PSO優(yōu)化粒子濾波算法

2 轉(zhuǎn)臺(tái)精度衰退預(yù)測(cè)實(shí)現(xiàn)過程

轉(zhuǎn)臺(tái)精度衰退預(yù)測(cè)的具體實(shí)現(xiàn)過程如下：

(1) 信號(hào)降噪與重構(gòu)，使用EEMD-PCA方法對(duì)原始數(shù)據(jù)降噪并進(jìn)行信號(hào)重構(gòu)[18]，通過使用聚合經(jīng)驗(yàn)?zāi)B(tài)法(EEMD)將原始振動(dòng)信號(hào)分解為若干組IMF(Intrinsic Mode Function)分量，再對(duì)這些組的IMF分量進(jìn)行PCA，從而篩選出其中包含低頻有效信息的IMF，再對(duì)篩選過的IMF分量重構(gòu)成矩陣Ej，從而得到反映進(jìn)度衰退的特征信號(hào)，同時(shí)起到降噪的目的。

(2) 選擇統(tǒng)計(jì)特征量，組成特征觀察值矩陣。通過Baum-Welch算法在已知的觀察序列的基礎(chǔ)上建立最有可能的HMM模型。本文通過觀察序列OB10建立一個(gè)HMM模型λ10，代入測(cè)試數(shù)據(jù)計(jì)算得到精度衰退的對(duì)數(shù)似然概率曲線并診斷是否發(fā)生衰退，在此基礎(chǔ)上建立健康狀態(tài)指標(biāo)H。本文選擇了峭度(KU)、偏度(SK)、裕度指標(biāo)(CL)、峰值指標(biāo)(CF)、脈沖指標(biāo)(IM)和波形指標(biāo)(SI)等六個(gè)統(tǒng)計(jì)特征量。

由這些統(tǒng)計(jì)特征量組成的觀察值矩陣為

(13)

式中：OBj=[OBj1,OBj2, …,OBjK]的j對(duì)應(yīng)下文的第j組數(shù)據(jù)。

定義描述轉(zhuǎn)臺(tái)精度衰退的健康狀態(tài)指標(biāo)為

H=P(OBj|λ10)

(14)

式中：λ10為用健康狀態(tài)下的觀察值矩陣OB10訓(xùn)練的HMM模型。

(3) 粒子濾波算法預(yù)測(cè)精度趨勢(shì)。在經(jīng)過粒子群算法優(yōu)化得到最優(yōu)的起始點(diǎn)參數(shù)ak,bk,ck的值后，使用粒子濾波算法對(duì)精度衰退趨勢(shì)和剩余壽命進(jìn)行預(yù)測(cè)，如表2所示。

表2 標(biāo)準(zhǔn)粒子濾波算法

其中，

(15)

訓(xùn)練的粒子數(shù)定為2 000，σ1，σ2，σ3，σ4分別為8，0.007, 0.005，9。粒子濾波算法根據(jù)狀態(tài)空間模型對(duì)精度狀態(tài)進(jìn)行迭代更新，輸出每一步的狀態(tài)估計(jì)值，在達(dá)到精度的閾值之后，根據(jù)精度的PDF分布以及預(yù)測(cè)失效點(diǎn)的位置計(jì)算數(shù)控轉(zhuǎn)臺(tái)的精度剩余壽命(Remaining Useful Life, RUL)。

轉(zhuǎn)臺(tái)精度衰退評(píng)估和剩余壽命計(jì)算的模型方案如圖1所示。

圖1 數(shù)控轉(zhuǎn)臺(tái)精度衰退模型方案Fig.1 Model of precision decline of NC rotary table

3 數(shù)控轉(zhuǎn)臺(tái)精度衰退加速壽命試驗(yàn)

3.1 試驗(yàn)設(shè)備簡(jiǎn)介

本次的試驗(yàn)對(duì)象為國內(nèi)某機(jī)床附件公司生產(chǎn)的TK13250E型數(shù)控立臥轉(zhuǎn)臺(tái)。工作臺(tái)最高轉(zhuǎn)速為11.1 r/min，定位精度為30"，重復(fù)定位精度為8"(調(diào)整前)，經(jīng)間隙補(bǔ)償后重復(fù)定位精度可達(dá)6"。如圖2(a)所示，該數(shù)控轉(zhuǎn)臺(tái)由西門子s7-200 224XP型PLC控制的步進(jìn)電機(jī)驅(qū)動(dòng)，該P(yáng)LC(Programmable Logic Controller)同時(shí)控制一個(gè)加載液壓缸，振動(dòng)加速度信號(hào)由PXI數(shù)據(jù)采集箱采集。加速度傳感器1安裝在蝸桿的軸線延長線位置上，加速度傳感器2安裝在電機(jī)傳動(dòng)齒輪的鏈接法蘭上。由于蝸輪副的狀態(tài)變化對(duì)轉(zhuǎn)臺(tái)的定位精度具有顯著影響，所以本文中利用加速度傳感器1采集的精度衰退全壽命信號(hào)進(jìn)行分析。本次試驗(yàn)為轉(zhuǎn)臺(tái)的變載荷加速壽命試驗(yàn)，在試驗(yàn)的前200 h以100%的載荷運(yùn)行轉(zhuǎn)臺(tái)，后520 h以200%載荷運(yùn)行轉(zhuǎn)臺(tái)，采樣頻率為2 048 Hz。將樣本分成了87組，其中以第10組樣本為訓(xùn)練樣本，其余為測(cè)試樣本，每組122 886個(gè)樣本。

為了驗(yàn)證該模型的有效性，在實(shí)驗(yàn)中使用激光干涉儀定期測(cè)量了轉(zhuǎn)臺(tái)的定位精度。實(shí)驗(yàn)設(shè)備為雷尼紹XL-80激光干涉儀系統(tǒng)，該設(shè)備回轉(zhuǎn)測(cè)量時(shí)分度精度為±5 μm/m，重復(fù)精度為1 μm/m。如圖2(b)所示，使用激光干涉儀每60 h測(cè)一次共測(cè)得12組數(shù)據(jù)。每次測(cè)量時(shí)需要將轉(zhuǎn)臺(tái)平放在操作臺(tái)上，安裝測(cè)量?jī)x器。試驗(yàn)時(shí)，數(shù)控轉(zhuǎn)臺(tái)以10°為步長，正反交替轉(zhuǎn)兩圈獲得一次測(cè)量值，每次試驗(yàn)測(cè)三次并取平均值。整個(gè)試驗(yàn)系統(tǒng)的原理圖如圖2(c)所示。

圖2 數(shù)控轉(zhuǎn)臺(tái)精度衰退加速壽命試驗(yàn)系統(tǒng)Fig.2 Accelerated life test system ofaccuracy decline of NC rotary table

3.2 試驗(yàn)結(jié)果分析

通過對(duì)比分析試驗(yàn)采集的87組原始數(shù)據(jù)，在第10組～第30組數(shù)據(jù)的振動(dòng)幅值十分平穩(wěn)，從第39組數(shù)據(jù)開始幅值出現(xiàn)較大的增加，并且出現(xiàn)規(guī)律性的變化，圖3(a)中對(duì)第10組、第39組和第50組數(shù)據(jù)進(jìn)行了對(duì)比，從中可以看出，振動(dòng)信號(hào)的幅值在全壽命周期中躍升了三個(gè)數(shù)量級(jí)。

圖3 第10組、第39組和第50組原始數(shù)據(jù)以及EEMD-PCA方法重構(gòu)前后的第50組數(shù)據(jù)Fig.3 Raw data at file number 10、39 and 50 and Fiftieth sets of data before and after reconstruction by EEMD-PCA

獲取了數(shù)據(jù)后使用EEMD-PCA方法對(duì)原始數(shù)據(jù)降噪并進(jìn)行信號(hào)重構(gòu)，圖3(b)中以第50組數(shù)據(jù)為例，給出了重構(gòu)前后的數(shù)據(jù)。從中可以看出，信號(hào)降噪和重構(gòu)之后出現(xiàn)了比較明顯的能反映精度變化的特征。同時(shí)，圖4給出了表2中的統(tǒng)計(jì)特征量在代入重構(gòu)信號(hào)時(shí)的計(jì)算結(jié)果，這些特征量對(duì)轉(zhuǎn)臺(tái)精度衰退的變化比較敏感，驗(yàn)證了表2的合理性。

將測(cè)試樣本輸入到經(jīng)過第10組數(shù)據(jù)訓(xùn)練的HMM模型中，計(jì)算出一組對(duì)數(shù)似然概率值，該值反映了測(cè)試樣本與訓(xùn)練樣本的匹配度，如圖5所示。選用第10組數(shù)據(jù)作為訓(xùn)練樣本主要是考慮到轉(zhuǎn)臺(tái)在運(yùn)轉(zhuǎn)初期有一個(gè)磨合的過程，其初期的精度不是轉(zhuǎn)臺(tái)運(yùn)轉(zhuǎn)的最佳狀態(tài)，經(jīng)過對(duì)比前后的數(shù)據(jù)，認(rèn)為第10組數(shù)據(jù)可以代表正常精度狀態(tài)下的振動(dòng)信號(hào)，并且由于本次試驗(yàn)的采樣頻率較高，所以一組數(shù)據(jù)的信息量足夠大。

圖4 統(tǒng)計(jì)特征量計(jì)算結(jié)果Fig.4 Calculation results of statistical characteristic quantities

從圖中可以看出，在第39組數(shù)據(jù)時(shí)，匹配度有一個(gè)比較大的下降，可以判斷轉(zhuǎn)臺(tái)精度從此時(shí)開始逐漸下降，本文將從第39組開始計(jì)算出的概率值作為健康狀態(tài)指標(biāo)H。同時(shí)對(duì)比圖5和圖6，可以看出在用激光干涉儀測(cè)量結(jié)果中第8次測(cè)量與第9次測(cè)量之間，轉(zhuǎn)臺(tái)的精度超過了設(shè)計(jì)允許的誤差精度值8"，時(shí)間是在第480 h～第540 h，對(duì)應(yīng)HMM計(jì)算的似然函數(shù)是在第67組，因此可以認(rèn)為函數(shù)值在-100為精度失效的閾值。

這里需要說明的是，圖6是利用激光干涉儀在數(shù)控轉(zhuǎn)臺(tái)精度衰退全壽命試驗(yàn)中測(cè)量得到的，它反映了整個(gè)精度衰退的過程的全貌，而在轉(zhuǎn)臺(tái)實(shí)際運(yùn)轉(zhuǎn)工作過程中，利用激光干涉儀獲得的僅僅是某一個(gè)或幾個(gè)點(diǎn)的精度值，它難以獲得對(duì)整個(gè)運(yùn)轉(zhuǎn)周期完整的描述，只能判斷是否精度失效，而無法反映出轉(zhuǎn)臺(tái)當(dāng)前的精度狀態(tài)。本文提出的算法模型通過對(duì)轉(zhuǎn)臺(tái)實(shí)際運(yùn)轉(zhuǎn)過程中監(jiān)測(cè)獲得的歷史數(shù)據(jù)進(jìn)行分析，可以得到轉(zhuǎn)臺(tái)精度衰退的趨勢(shì)，描繪出轉(zhuǎn)臺(tái)所處的精度狀態(tài)，并且可以以某一個(gè)精度未失效的點(diǎn)為起始點(diǎn)，來預(yù)測(cè)轉(zhuǎn)臺(tái)精度的變化趨勢(shì)，判斷出精度的剩余壽命，由此判斷是否要對(duì)轉(zhuǎn)臺(tái)進(jìn)行必要的調(diào)整。

圖5 數(shù)控轉(zhuǎn)臺(tái)精度衰退健康指標(biāo)曲線Fig.5 Health indicators curve of accuracy decline of NC rotary table

圖6 激光干涉儀測(cè)量重復(fù)定位精度值Fig.6 Measurement accuracy of repeated positioning by laser interferometer

圖7 三種算法預(yù)測(cè)結(jié)果圖Fig.7 Prediction results of three algorithms

方法平均絕對(duì)誤差均方誤差計(jì)算時(shí)間/sPSO優(yōu)化初始參數(shù)0．5750．88943．67未優(yōu)化的粒子濾波1．4753．53110．33PSO優(yōu)化粒子集1．3173．28339．89

表4給出以不同信號(hào)分組作為預(yù)測(cè)起始點(diǎn)的情況下三種算法的預(yù)測(cè)結(jié)果。其中，從第50組數(shù)據(jù)開始預(yù)測(cè)時(shí)，結(jié)合圖7可以看出，三種方法的不確定度為50%時(shí)的預(yù)測(cè)結(jié)果分別是71、75和60，則可得到其剩余壽命(RUL)分別為21、25和10。通過激光干涉儀測(cè)量精度變化獲得的實(shí)際的精度剩余壽命為17，差值分別為4、8和7。綜合對(duì)比不同預(yù)測(cè)起始點(diǎn)的預(yù)測(cè)結(jié)果后認(rèn)為本文提出的方法更加接近實(shí)際的結(jié)果。圖8給出了本文提出的優(yōu)化粒子濾波初始參數(shù)的方法在不同預(yù)測(cè)起始點(diǎn)的情況下，預(yù)測(cè)結(jié)果的PDF的分布曲線，可以看出起始點(diǎn)越大其分布越集中，說明其不確定度在變小，與實(shí)際情況是一致的。

表4 不同起始點(diǎn)剩余壽命預(yù)測(cè)結(jié)果

圖8 PSO優(yōu)化粒子濾波算法預(yù)測(cè)結(jié)果的PDF圖Fig.8 PDF results of PF initial parameters optimizedby PSO

4 結(jié) 論

為了解決數(shù)控轉(zhuǎn)臺(tái)的精度衰退缺乏有效的描述方法且現(xiàn)有測(cè)量?jī)x器只是測(cè)量工具而無法實(shí)現(xiàn)狀態(tài)評(píng)估等問題。本文引入數(shù)據(jù)驅(qū)動(dòng)的方法，提出數(shù)控轉(zhuǎn)臺(tái)重復(fù)定位精度衰退趨勢(shì)預(yù)測(cè)模型，并得出如下的結(jié)論：

(1)本文區(qū)別于以往的狀態(tài)識(shí)別，將HMM算法以新的角度運(yùn)用于數(shù)控轉(zhuǎn)臺(tái)的精度衰退的研究中，結(jié)合優(yōu)化后的粒子濾波算法，建立了一種預(yù)測(cè)轉(zhuǎn)臺(tái)精度衰退趨勢(shì)的模型。在以第50組為預(yù)測(cè)起始點(diǎn)時(shí)，預(yù)測(cè)剩余壽命為21，實(shí)際剩余壽命為17，相差4，說明預(yù)測(cè)結(jié)果準(zhǔn)確度較高。

(2)在使用粒子濾波算法的過程中，對(duì)粒子濾波算法的預(yù)測(cè)點(diǎn)初始參數(shù)使用PSO算法進(jìn)行了優(yōu)化，找出了最優(yōu)的結(jié)果，并與未優(yōu)化的粒子濾波算法以及用PSO優(yōu)化粒子濾波粒子集方法的結(jié)果進(jìn)行了對(duì)比分析，結(jié)果表明，雖然PSO優(yōu)化初始參數(shù)的程序計(jì)算過程需要更長的時(shí)間，但是以第50組為預(yù)測(cè)起始點(diǎn)時(shí)，其預(yù)測(cè)結(jié)果與實(shí)際結(jié)果最接近，相差4，其他兩種方法分別相差8和7。其他預(yù)測(cè)起始點(diǎn)時(shí)也是相同的情況。

[ 1 ] YU C J, HUANG X D, FANG C G, et al. Research on damping and vibration characteristic of the large and precision NC rotary table[J]. Advanced Materials Research, 2010, 97/98/99/100/101: 1216-1222.

[ 2 ] ?VOKELJ M, ZUPAN S, PREBIL I. Non-linear multivariate and multiscale monitoring and signal denoising strategy using kernel principal component analysis combined with ensemble empirical mode decomposition method [J]. Mechanical Systems & Signal Processing, 2011, 25(7): 2631-2653.

[ 3 ] YU C J, HUANG X D, FANG C G. Research on dynamic characteristics of NC rotary table considering leakage factors of its hydrostatic guideway[J]. Journal of Mechanical Engineering Science, 2012, 226(11): 2674-2685.

[ 4 ] HENG A, ZHANG S, TAN A C C, et al. Rotating machinery prognostics: state of the art, challenges and opportunities [J]. Mechanical Systems & Signal Processing, 2009, 23(3): 724-739.

[ 5 ] BENKEDJOUH T, MEDJAHER K, ZERHOUNI N. Remaining useful life estimation based on nonlinear feature reduction and support vector regression [J]. Engineering Applications of Artificial Intelligence, 2013, 26(7): 1751-1760.

[ 6 ] LU C, CHEN J, HONG R, et al. Degradation trend estimation of slewing bearing based on LSSVM model [J]. Mechanical Systems & Signal Processing, 2016, 76/77:

353-366.

[ 7 ] MATEJ Z, SAMO Z, IVAN P. Non-linear multivariate and multiscale monitoring and signal denoising strategy using kernel principal component analysis combined with ensemble empirical mode decomposition method[J]. Mechanical Systems and Signal Processing, 2011, 25(7): 2631-2653.

[ 8 ] FENG Y, HUANG X, HONG R, et al. Residual useful life prediction of large-size low-speed slewing bearings-a data driven method[J]. Journal of Vibroengineering, 2015, 17(8): 4164-4179.

[ 9 ] 于春健.多物理場(chǎng)耦合的大重型數(shù)控回轉(zhuǎn)工作臺(tái)動(dòng)靜態(tài)精度研究[D]. 南京：南京工業(yè)大學(xué)，2012.

[10] WANG D, MIAO Q, ZHOU Q, et al. An intelligent prognostic system for gear performance degradation assessment and remaining useful life estimation [J]. Journal of Vibration & Acoustics, 2015, 137(2): 12.

[11] 張西寧, 雷威, 李兵. 主分量分析和隱馬爾科夫模型結(jié)合的軸承監(jiān)測(cè)診斷方法[J]. 西安交通大學(xué)學(xué)報(bào),2017, 51(6): 1-7.

ZHANG Xining, LEI Wei, LI Bing. Bearing fault detection and diagnosis method based on principal component analysis and hidden markov model [J].Journal of Xi’an Jiaotong University, 2017, 51(6)： 1-7.

[12] 余龍華,王宏,鐘洪聲. 基于隱馬爾科夫模型的人臉識(shí)別[J]. 計(jì)算機(jī)技術(shù)與發(fā)展,2012, 22(2): 25-28.

YU Longhua，WANG Hong，ZHONG Hongsheng. Face recognition based on hidden markov model [J].Computer Technology and Development, 2012, 22(2): 25-28.

[13] RABINER L R. A tutorial on hidden Markov models and selected applications in speech recognition [J]. Readings in Speech Recognition, 1990, 77(2): 267-296.

[14] ARULAMPALAM M S, MASKELL S, GORDON N. A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking[J]. IEEE Transactions on Signal Processing, 2002, 50(2): 174-188.

[15] 王法勝, 魯明羽, 趙清杰, 等. 粒子濾波算法[J]. 計(jì)算機(jī)學(xué)報(bào),2014,37(8): 1679-1694.

WANG Fasheng, LU Mingyu, ZHAO Qingjie, et al. Particle filtering algorithm [J]. Chinese Journal of Computer, 2014, 37(8): 1679-1694.

[16] 方正, 佟國峰, 徐心和. 粒子群優(yōu)化粒子濾波方法[J]. 控制與決策,2007, 22(3): 273-277.

FANG Zheng, TONG Guofeng, XU Xinhe. Particle swarm optimized particle filter [J].Control and Decision, 2007, 22(3): 273-277.

[17] 王維博, 林川, 鄭永康. 粒子群算法中參數(shù)的實(shí)驗(yàn)與分析[J]. 西華大學(xué)學(xué)報(bào)(自然科學(xué)版),2008, 27(1): 76-80.

WANG Weibo, LIN Chuan, ZHENG Yongkang. Experiment and analysis of parameters in particle swarm optimization [J]. Journal of Xi Hua University (Natural Science), 2008, 27(1): 76-80.

[18] 封楊, 黃筱調(diào), 陳捷, 等. 大型回轉(zhuǎn)支承非平穩(wěn)振動(dòng)信號(hào)的EEMD-PCA降噪方法[J]. 南京工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 37(3): 61-66.

FENG Yang, HUANG Xiaodiao, CHEN Jie, et al. EEMD-PCA based denosing method for non-stationary vibration signals of large-size slewing bearings [J].Journal of Nanjing Tech University(Natural Science) ,2015, 37(3):61-66.