技術(shù)特征：

1.一種基于聲發(fā)射信號的金剛石壓機頂錘裂紋在線檢測方法，其特征在于，綜合利用BSA-SVM智能檢測方法和K-SPWVD時頻信號分析方法進(jìn)行頂錘裂紋檢測，具體步驟如下：

第一部分檢測依據(jù)的建立

步驟1：在已知是否有裂紋的頂錘上安裝聲發(fā)射傳感器，首先對已知是否有裂紋的頂錘進(jìn)行聲發(fā)射信號采集，通過電路對傳感器采集的信號進(jìn)行預(yù)處理，將預(yù)處理后的信號輸入計算機；在計算機中對頂錘正常狀態(tài)下的聲發(fā)射信號和裂紋狀態(tài)下的聲發(fā)射信號進(jìn)行收集；

步驟2：計算機對上述頂錘正常狀態(tài)下的聲發(fā)射信號和裂紋狀態(tài)下的聲發(fā)射信號，利用回溯搜索優(yōu)化算法(BSA)對支持向量機模型(SVM)進(jìn)行參數(shù)的優(yōu)化和訓(xùn)練，得到金剛石壓機頂錘裂紋智能檢測模型，該模型名詞定為：基于回溯搜索優(yōu)化算法的支持向量機模型(BSA-SVM)；

步驟3：

利用上述步驟1中采集到的已有的裂紋狀態(tài)下的聲發(fā)射信號，

(a)首先，取M個數(shù)據(jù)點的頂錘裂紋狀態(tài)聲發(fā)射信號作為模板信號；

(b)然后，利用平滑偽維格納-威利分布(SPWVD)處理分析該模板信號，根據(jù)得到的頻譜信息確定頂錘裂紋特征頻率；

第二部分實時檢測過程

步驟4：

將聲發(fā)射傳感器安裝到正在工作的頂錘上，針對采集到的正在工作中的頂錘聲發(fā)射實時信號，首先利用基于回溯搜索優(yōu)化算法的支持向量機模型(BSA-SVM)進(jìn)行識別，

如果檢測到是非裂紋信號，則執(zhí)行步驟6；

如果檢測到是裂紋信號，暫定為一級疑似裂紋信號，同時執(zhí)行步驟5；

步驟5：

將上述一級疑似裂紋信號利用柯爾莫諾夫-斯米爾諾夫檢驗(Kolmogorov-Smirnov test，簡稱K-S檢驗)與平滑偽維格納-威利分布(K-SPWVD)結(jié)合處理方法進(jìn)行處理，

該柯爾莫諾夫-斯米爾諾夫檢驗與平滑偽維格納-威利分布(K-SPWVD)結(jié)合處理方法的工作過程包括

(i)首先對采集的一級疑似裂紋信號利用進(jìn)行連續(xù)平移L點長度的加窗截取，窗口長度為M，利用柯爾莫諾夫-斯米爾諾夫檢驗(Kolmogorov-Smirnov test，簡稱K-S檢驗)方法與所述的模板信號進(jìn)行匹配；如果不匹配，說明是非裂紋信號，執(zhí)行步驟7；如果匹配，暫定為二級疑似裂紋信號，并轉(zhuǎn)到下一步(ii)；

(ii)再將上述二級疑似裂紋信號利用平滑偽維格納-威利分布(SPWVD)處理分析所截信號，根據(jù)得到的頻譜信息確定其特征頻率，如果能觀察到有步驟3(b)中的頂錘裂紋特征頻率，則判定為頂錘發(fā)生了裂紋，執(zhí)行步驟7，反之，則判定為非裂紋信號，執(zhí)行步驟6；

步驟6：

輸出頂錘沒有裂紋的判斷結(jié)果；

步驟7：

輸出頂錘有裂紋的結(jié)果。

2.根據(jù)權(quán)利要求1所述的一種基于聲發(fā)射信號的金剛石壓機頂錘裂紋在線檢測方法，其特征在于：所述步驟2中的基于回溯搜索優(yōu)化算法的支持向量機模型(BSA-SVM)的建立過程包括，

a.確定訓(xùn)練樣本

從步驟1中采集的頂錘聲發(fā)射信號集中選取n₁組已知正常信號和n₂組已知疑似裂紋信號作為無裂紋狀態(tài)，n₃組已知裂紋信號對應(yīng)裂紋狀態(tài)，此處已知疑似裂紋信號為頂錘沒有裂紋但信號波形與裂紋信號相似的信號；無裂紋狀態(tài)和裂紋狀態(tài)兩種類樣本共同組成訓(xùn)練樣本，利用小波包對信號進(jìn)行三層分解，提取各頻段的能量百分比作為樣本特征值，樣本特征值矩陣構(gòu)成支持向量機模型的輸入向量，模型的輸出為1或-1，分別代表無裂紋和裂紋狀態(tài)；

b.確定支持向量機的待優(yōu)化參數(shù)

支持向量機中有2個參數(shù)：核函數(shù)參數(shù)σ和懲罰函數(shù)參數(shù)C，標(biāo)記為p(σ，C)；

c.基于BSA的參數(shù)初始化

BSA的初始化是指初始化種群P和歷史種群oldP；種群包括N個個體，個體的長度設(shè)置為d＝2，代表參數(shù)的個數(shù)；種群P初始化表達(dá)式為：

P_i,j～U(low_j,up_j)

其中i＝1，2，3，…，N；j＝1，2，3，…，D；N和D分別代表種群的規(guī)模和問題維度，U(,)表示均勻分布函數(shù)，up_j和low_j分別表示第j維變量的上、下界；

歷史種群oldP初始化是防止首次運行篩選Ⅰ時，oldP數(shù)值為空集；歷史種群oldP由以下公式?jīng)Q定：

oldP_i,j～U(low_j,up_j)

d.篩選Ⅰ

根據(jù)概率運算重新定義oldP，其表達(dá)式為：

if a＜b then oldP＝P|a,b～U(0,1)

式中等號表示更新代替；新一代的歷史種群oldP是上一代或者是m代前的歷史種群，m＝1,2,…M-1，M為目前迭代次數(shù)；

e.突變

BSA通過下式形成試驗種群的初始形式Mutant

Mutant＝P+F·(oldP-P)

式中F作為變異尺度系數(shù)用于控制變異的幅度，F(xiàn)＝3·randn，其中randn～N(0,1)，N(,)表示標(biāo)準(zhǔn)正態(tài)分布函數(shù)；

f.交叉

BSA通過交叉過程產(chǎn)生試驗種群的最終形式T；該交叉方式是一種基于兩種交叉方式等概率調(diào)用的聯(lián)合交叉方式；

在交叉時，首先是交叉長度n的選擇，然后在歷史種群P中每個個體隨機篩選n個特征并與Mutant中等位置等維度的特征進(jìn)行交換，得到新的個體；n是(0,mixrate·D]中的整數(shù)，式中，mixrate是交叉概率，一般取為1，D是問題維度；

在每代的交叉中，n有兩種可能，由n的賦值方式?jīng)Q定，分為：

①交叉方式1：n＝1，即任意選取一維，互換歷史種群中個體與同位置試驗種群個體；

或者

②交叉方式2：給每個個體隨機生成一個randn～N(0,1)，把

n(i)＝[mixrate·randn·D]

定義成交叉長度；式中，mixrate是交叉概率，D是問題維度；

上述兩種交叉方式等概率隨機調(diào)用構(gòu)成BSA的交叉方式；

g.篩選II

將試驗種群帶入支持向量機，產(chǎn)生試驗種群的適應(yīng)度函數(shù)，其表達(dá)式為：

$F=\frac{1}{2}\underset{i=1}{\overset{n}{\Σ}}\left|abs(species-groupnames)\right|$

式中：F表示種群中每個個體的適應(yīng)度，species表示樣本分類結(jié)果，groupnames表示樣本對應(yīng)狀態(tài)；

比較當(dāng)前種群與試驗種群間的適應(yīng)度函數(shù)，用試驗種群中的優(yōu)秀個體替換當(dāng)前種群中相應(yīng)的個體；如果本次最優(yōu)個體Pbest具有比全局最優(yōu)個體更好的適應(yīng)度，則全局最優(yōu)個體就用Pbest代替，同樣其對應(yīng)的適應(yīng)度值也用本次最優(yōu)個體的適應(yīng)度值代替；

h.如果本次最優(yōu)個體滿足預(yù)先設(shè)定的分類精度或者迭代次數(shù)，不再進(jìn)行種群更新，退出訓(xùn)練；否則轉(zhuǎn)入步驟(d)，更新種群。