權(quán)利要求書： 1.一種球磨機(jī)料位預(yù)測方法，其特征在于，包括：s v

向球磨機(jī)滾筒中填充物料，采集球磨機(jī)在不同料位的振動信號X及振聲信號X；

s v s v

根據(jù)采集到的振動信號X及振聲信號X計(jì)算振動信號X及振聲信號X的特征參數(shù)；

s v

將振動信號X及振聲信號X的特征參數(shù)作為TBM模型的輸入，訓(xùn)練TBM模型，通過多次訓(xùn)s v

練，在使TBM模型的優(yōu)化函數(shù)取最小值時(shí)，建立振動信號X及振聲信號X的特征參數(shù)與預(yù)測料位的TBM模型，并通過特征參數(shù)計(jì)算預(yù)測料位；具體包括：在訓(xùn)練集中尋找k個(gè)最鄰近樣本并得到每個(gè)樣本對應(yīng)的信度，第k個(gè)樣本對應(yīng)的信度為

2 2

其中， Φ(d)＝γexp(?d)，d＝g(x?xi) ，xi為鄰近值，x為訓(xùn)練集中的樣本，γ為常數(shù)；

對計(jì)算得到的k個(gè)信度 進(jìn)行DS合成，合成結(jié)果為：根據(jù)DS合成結(jié)果計(jì)算預(yù)測料位，計(jì)算結(jié)果為：重復(fù)前述步驟，當(dāng)優(yōu)化函數(shù)的值w達(dá)到最小值時(shí)，得到參數(shù)g，并計(jì)算得到預(yù)測料位。

2.根據(jù)權(quán)利要求1所述的球磨機(jī)料位預(yù)測方法，其特征在于，在根據(jù)采集到的振動信號s v s v

X及振聲信號X計(jì)算振動信號X及振聲信號X的特征參數(shù)的步驟中，包括步驟：s v s v

對輸入的振動信號X及振聲信號X 進(jìn)行預(yù)處理，得到振動信號X及振聲信號X的分幀信號；預(yù)處理方式包括：預(yù)加重、分幀和加窗處理；

s v s v

對預(yù)處理后的每幀信號進(jìn)行FFT變換，將時(shí)域信號x(n)、x(n)轉(zhuǎn)換為頻域信號x (k)、xs v

(k)，再計(jì)算其模的平方得到能量譜P(k)、P(k)；

s v

對所得能量譜P(k)、P(k)進(jìn)行Mel濾波，并取對數(shù)，得到：其中，其中L為濾波器個(gè)數(shù)，Hi(k)、Hj(k)三角濾波器的傳遞參數(shù)；

s

對梅爾濾波的結(jié)果進(jìn)行離散余弦變換，得到MFCC系數(shù)A(n)作為振動信號X 及振聲信號v

X的特征參數(shù)，計(jì)算公式為

3.根據(jù)權(quán)利要求1所述的球磨機(jī)料位預(yù)測方法，其特征在于，TBM模型訓(xùn)練的優(yōu)化函數(shù)為

其中 是預(yù)測料位，y是真實(shí)料位。

4.根據(jù)權(quán)利要求1所述的球磨機(jī)料位預(yù)測方法，其特征在于，在采集球磨機(jī)在相應(yīng)料位s v

的振動信號X及振聲信號X的步驟中，是通過傳感器采集球磨機(jī)筒壁或者軸承上的振動信號，及通過安裝在球磨機(jī)周圍麥克風(fēng)采集振聲信號。

5.根據(jù)權(quán)利要求4所述的球磨機(jī)料位預(yù)測方法，其特征在于，通過傳感器和麥克風(fēng)采集s v

的振動信號和振聲信號在計(jì)算振動信號X及振聲信號X的特征參數(shù)之前，將振動信號和振聲信號通過模數(shù)轉(zhuǎn)換轉(zhuǎn)換為數(shù)字信號。

6.根據(jù)權(quán)利要求1所述的球磨機(jī)料位預(yù)測方法，其特征在于，在得到預(yù)測料位后，采用均方根誤差(RMSE)作為模型性能評價(jià)指標(biāo)，判斷預(yù)測料位結(jié)果的準(zhǔn)確性；其中，均方根誤差(RMSE)的計(jì)算公式為：式中， 和yi分別表示第i個(gè)樣本的料位估計(jì)值和實(shí)際值，N為在線測試樣本個(gè)數(shù)。

說明書： 球磨機(jī)料位預(yù)測方法技術(shù)領(lǐng)域[0001] 本發(fā)明涉及數(shù)據(jù)檢測領(lǐng)域，尤其涉及一種球磨機(jī)料位預(yù)測方法。背景技術(shù)[0002] 滾筒式球磨機(jī)是一種基本的制粉設(shè)備。球磨機(jī)滾筒內(nèi)部的料位是影響球磨機(jī)能量消耗、生產(chǎn)效率和生產(chǎn)安全的關(guān)鍵因素，然而在實(shí)際工業(yè)生產(chǎn)過程中，當(dāng)球磨機(jī)系統(tǒng)運(yùn)轉(zhuǎn)

時(shí)，由于其筒內(nèi)環(huán)境惡劣并且情況比較復(fù)雜，球磨機(jī)料位難以準(zhǔn)確預(yù)測，導(dǎo)致球磨機(jī)系統(tǒng)難

以經(jīng)濟(jì)運(yùn)行并且容易發(fā)生安全事故。因此，有效預(yù)測球磨機(jī)的料位，對于球磨機(jī)系統(tǒng)高效和

安全穩(wěn)定運(yùn)行具有重大意義。

[0003] 目前，基于軟測量的球磨機(jī)料位測量方法主要是通過單一的信號源如壓電信號、差壓信號、振聲信號、振動信號或者功率信號來預(yù)測料位。從反復(fù)試驗(yàn)中我們發(fā)現(xiàn)單一信號

源具有一定的局限性例如：壓電信號使用時(shí)間的不斷增長，壓電應(yīng)變傳感器很容易出現(xiàn)零

點(diǎn)漂移現(xiàn)象，導(dǎo)致測量精度降低；差壓信號需要的輔助變量較多，各個(gè)變量值難以準(zhǔn)確獲

取，各變量之間存在冗余和耦合，并且系統(tǒng)復(fù)雜度較高，不方便進(jìn)行維護(hù)；振聲信號容易受

到來自毗鄰磨機(jī)噪音的干擾，以致于料位識別效率不高；振動信號容易受到電網(wǎng)、磨內(nèi)研磨

體損耗等因素的影響，從而使負(fù)荷控制點(diǎn)發(fā)生漂移現(xiàn)象；功率信號受電流影響較大，而電流

又受鋼球的磨損、鋼球總量的減少、煤質(zhì)的變化、風(fēng)量的變化的影響，因而推斷結(jié)果存在較

大誤差。

發(fā)明內(nèi)容[0004] 本發(fā)明的目的在于避免現(xiàn)有技術(shù)的不足之處而提供一種球磨機(jī)料位預(yù)測方法。[0005] 本發(fā)明的目的可以通過采用如下的技術(shù)措施來實(shí)現(xiàn)，設(shè)計(jì)一種球磨機(jī)料位預(yù)測方s

法，該方法的步驟包括：向球磨機(jī)滾筒中填充物料，采集球磨機(jī)在不同料位的振動信號X 及

v s v s v

振聲信號X；根據(jù)采集到的振動信號X及振聲信號X計(jì)算振動信號X 及振聲信號X的特征參

s v

數(shù)；將振動信號X及振聲信號X的特征參數(shù)作為TBM模型的輸入，訓(xùn)練TBM模型，通過多次訓(xùn)

s v

練，在使TBM模型的優(yōu)化函數(shù)取最小值時(shí)，建立振動信號X及振聲信號X的特征參數(shù)與預(yù)測

料位的TBM模型，并通過特征參數(shù)計(jì)算預(yù)測料位。

[0006] 其中，在根據(jù)采集到的振動信號Xs及振聲信號Xv計(jì)算振動信號Xs及振聲信號Xv的s v

特征參數(shù)的步驟中，包括步驟：對輸入的振動信號X 及振聲信號X進(jìn)行預(yù)處理，得到振動信

s v

號X及振聲信號X的分幀信號；預(yù)處理方式包括：預(yù)加重、分幀和加窗處理；對預(yù)處理后的每

s v s v

幀信號進(jìn)行FFT變換，將時(shí)域信號x (n)、x (n)轉(zhuǎn)換為頻域信號x (k)、x (k)，再計(jì)算其模的

s v s v

平方得到能量譜P(k)、P(k)；對所得能量譜P(k)、P(k)進(jìn)行Mel濾波，并取對數(shù)，得到：

[0007][0008] 其中，其中L為濾波器個(gè)數(shù)，Hi(k)、Hj(k)三角濾波器的傳遞參數(shù)；對梅爾濾波的結(jié)s v

果進(jìn)行離散余弦變換，得到MFCC系數(shù)A(n)作為振動信號X及振聲信號X的特征參數(shù)，計(jì)算公

式為

[0009][0010] 其中，在將振動信號Xs及振聲信號Xv的特征參數(shù)作為TBM模型的輸入，訓(xùn)練TBM模s v

型，通過多次訓(xùn)練，在使TBM模型的優(yōu)化函數(shù)取最小值時(shí)，建立振動信號X 及振聲信號X的特

征參數(shù)與預(yù)測料位的TBM模型，并通過特征參數(shù)計(jì)算預(yù)測料位的步驟中，包括步驟：

[0011] 在訓(xùn)練集中尋找k個(gè)最鄰近樣本并得到相應(yīng)的信度 其中，2 2

Φ(d)＝γexp(?d)，d＝g(x?xi) ，xi為鄰近值，x為訓(xùn)練集中的樣本，γ為常數(shù)；對計(jì)算得到

的k個(gè)信度 進(jìn)行DS合成，合成結(jié)果為：

[0012][0013] 根據(jù)DS合成結(jié)果計(jì)算預(yù)測料位，計(jì)算結(jié)果為：[0014][0015] 重復(fù)前述步驟，當(dāng)優(yōu)化函數(shù)的值w達(dá)到最小值時(shí)，得到參數(shù)g，并計(jì)算得到預(yù)測料位。

[0016] 其中，模型訓(xùn)練的優(yōu)化函數(shù)為[0017] 其中 是預(yù)測料位，y是真實(shí)料位。[0018] 其中，在采集球磨機(jī)在相應(yīng)料位的振動信號Xs及振聲信號Xv的步驟中，是通過傳感器采集球磨機(jī)筒壁或者軸承上的振動信號，及通過安裝在球磨機(jī)周圍麥克風(fēng)采集振聲信

號。

[0019] 其中，通過傳感器和麥克風(fēng)采集的振動信號和振聲信號在計(jì)算梅爾頻率倒譜系數(shù)之前，將振動信號和振聲信號通過模數(shù)轉(zhuǎn)換轉(zhuǎn)換為數(shù)字信號。

[0020] 其中，在得到預(yù)測料位后，采用均方根誤差(RMSE)作為模型性能評價(jià)指標(biāo)，判斷預(yù)測料位結(jié)果的準(zhǔn)確性；其中，均方根誤差(RMSE)的計(jì)算公式為：

[0021][0022] 式中， 和yi分別表示第i個(gè)樣本的料位估計(jì)值和實(shí)際值，N為在線測試樣本個(gè)數(shù)。[0023] 本發(fā)明的目的可以通過采用如下的技術(shù)措施來實(shí)現(xiàn)，設(shè)計(jì)一種球磨機(jī)料位預(yù)測裝s

置，連接于待預(yù)測的球磨機(jī)，包括：采集模塊，用于采集球磨機(jī)在不同料位的振動信號X 及

v s v s

振聲信號X ；第一計(jì)算模塊，用于根據(jù)采集到的振動信號X及振聲信號X 計(jì)算振動信號X 及

v s v

振聲信號X 的特征參數(shù)；第二計(jì)算模塊，用于將振動信號X 及振聲信號X的特征參數(shù)作為

TBM模型的輸入，訓(xùn)練TBM模型，通過多次訓(xùn)練，在使TBM模型的優(yōu)化函數(shù)取最小值時(shí)，建立振

s v

動信號X及振聲信號X的特征參數(shù)與預(yù)測料位的TBM模型，并通過特征參數(shù)計(jì)算預(yù)測料位。

[0024] 其中，第一計(jì)算模塊包括：預(yù)處理單元，用于振動信號Xs及振聲信號Xv進(jìn)行預(yù)處理；其中，預(yù)處理的方式包括預(yù)加重、分幀和加窗處理；FFT變換單元，用于對預(yù)處理后的振動信

s v

號X及振聲信號X進(jìn)行快速傅里葉變換，將時(shí)域信號轉(zhuǎn)換為頻域信號，并計(jì)算頻域信號的模

的平方，得到能量譜；Mel濾波單元，用于對FFT變換單元處理后的能量譜進(jìn)行Mel濾波，并取

對數(shù)；DCT變換單元，用于對Mel濾波處理后的結(jié)果進(jìn)行離散余弦變換，得到梅爾頻率倒譜系

s v

數(shù)作為振動信號X及振聲信號X的特征參數(shù)。

[0025] 區(qū)別于現(xiàn)有技術(shù)，本發(fā)明的球磨機(jī)料位預(yù)測方法通過采集球磨機(jī)筒壁或者軸承上的振動信號及安裝在球磨機(jī)周圍麥克風(fēng)的振聲信號；振動信號、振聲信號經(jīng)調(diào)理電路、模數(shù)

轉(zhuǎn)換電路輸入數(shù)據(jù)采集卡；通過MFCC分別計(jì)算數(shù)據(jù)采集卡中振動、振聲信號的特征參數(shù)；重

新組合振動、振聲信號；數(shù)據(jù)采集卡收集重新組合后的信號；最后數(shù)據(jù)采集卡中的數(shù)據(jù)通過

歷史數(shù)據(jù)建立的TBM模型，預(yù)測相應(yīng)信號的料位。本發(fā)明能有效的預(yù)測球磨機(jī)的料位，其預(yù)

測料位精確度高、波動小，具有良好的實(shí)用價(jià)值和推廣前景。

附圖說明[0026] 圖1是本發(fā)明提供的一種球磨機(jī)料位預(yù)測方法的流程示意圖。[0027] 圖2是本發(fā)明提供的一種球磨機(jī)料位預(yù)測裝置的結(jié)構(gòu)示意圖。[0028] 圖3是本發(fā)明提供的一種球磨機(jī)料位預(yù)測裝置的硬件結(jié)構(gòu)示意圖。具體實(shí)施方式[0029] 下面結(jié)合具體實(shí)施方式對本發(fā)明的技術(shù)方案作進(jìn)一步更詳細(xì)的描述。顯然，所描述的實(shí)施例僅僅是本發(fā)明的一部分實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，

本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動的前提下所獲得的所有其他實(shí)施例，都應(yīng)屬于

本發(fā)明保護(hù)的范圍。

[0030] 參閱圖1，圖1是本發(fā)明提供的一種球磨機(jī)料位預(yù)測方法的流程示意圖。該方法的步驟包括：

[0031] S110：向球磨機(jī)滾筒中填充物料，采集球磨機(jī)在不同料位的振動信號Xs及振聲信v

號X。

[0032] 本發(fā)明的方法涉及MFCC和TBM兩種理論方法。MFCC是一種提取聲音特征的方法。其魯棒性強(qiáng)，可用于分析非線性信號。針對信號分析和特征提取困難的問題，將近年來廣泛研

究的MFCC引入球磨機(jī)料位軟測量建模過程中，用于球磨機(jī)軸承振動信號以及振聲信號的有

效特征提取。

[0033] TBM是關(guān)于證據(jù)理論的一種解釋。它是一種雙層結(jié)構(gòu)?！癱redal層”位于底層，在該層中獲取信度并對其進(jìn)行量化、賦值和更新處理。“pignistic層”位于上層，它將credal層

上的信度轉(zhuǎn)換成pignistic概率，并由此做出決策。TBM比較適用于需要逐層進(jìn)行數(shù)據(jù)、特征

和決策層融合的數(shù)據(jù)融合系統(tǒng)。因此應(yīng)用于振聲信號、振動信號的融合。

[0034] 本發(fā)明所述的球磨機(jī)料位預(yù)測方法在具體實(shí)施時(shí)，首先在球磨機(jī)筒壁或者軸承上安裝振動傳感器收集振動信號，在球磨機(jī)周圍安裝麥克風(fēng)收集振聲信號。具體的，通過逐步

向球磨機(jī)料倉中填充物料，逐漸增加料位，直至最大料位，并通過振動傳感器、及作為振聲

s v

傳感器的麥克風(fēng)，采集振動信號X 、振聲信號X ，同時(shí)記錄相應(yīng)的料位值Y＝z，其中z＝1,

2,…,q。

[0035] S120：根據(jù)采集到的振動信號Xs及振聲信號Xv計(jì)算振動信號Xs及振聲信號Xv的特征參數(shù)。

[0036] 將采集到的振動信號Xs及振聲信號Xv作為輸入信號，并對輸入的信號進(jìn)行預(yù)加重、s v

分幀和加窗處理，得到若干分幀信號。x (n)、x (n)表示振動信號、振聲信號其中任意一幀。

預(yù)處理完成后，還包括對振動信號、振聲信號進(jìn)行模數(shù)轉(zhuǎn)換的步驟。

[0037] 對預(yù)處理后的每幀信號進(jìn)行FFT變換，將時(shí)域信號xs(n)、xv(n)轉(zhuǎn)換為頻域信號xsv s v

(k)、x(k)，再計(jì)算其模的平方得到能量譜P(k)、P(k)。

[0038] 對所得能量譜Ps(k)、Pv(k)進(jìn)行Mel濾波，并取對數(shù)，得到：[0039][0040] 其中，其中L為濾波器個(gè)數(shù)，Hi(k)、Hj(k)三角濾波器的傳遞參數(shù)。[0041] 對梅爾濾波的結(jié)果進(jìn)行離散余弦變換，得到MFCC系數(shù)A(n)，計(jì)算公式為[0042][0043] 式中A(n)為MFCC系數(shù)，J為MFCC階數(shù)，本發(fā)明中取為12。[0044] S130：將振動信號Xs及振聲信號Xv的特征參數(shù)作為TBM模型的輸入，訓(xùn)練TBM模型，s v

通過多次訓(xùn)練，在使TBM模型的優(yōu)化函數(shù)取最小值時(shí)，建立振動信號X及振聲信號X的特征

參數(shù)與預(yù)測料位的TBM模型，并通過特征參數(shù)計(jì)算預(yù)測料位。

[0045] 在訓(xùn)練集中尋找k個(gè)最鄰近樣本并得到相應(yīng)的信度 其中，2 2

Φ(d)＝γexp(?d)，d＝g(x?xi) ，xi為鄰近值，x為訓(xùn)練集中的樣本，可取步驟S120中計(jì)算

得到的MFCC系數(shù)A(n)，γ為常數(shù)。

[0046] 對計(jì)算得到的k個(gè)信度 進(jìn)行DS合成，合成結(jié)果為：[0047][0048] 根據(jù)DS合成結(jié)果計(jì)算預(yù)測料位，計(jì)算結(jié)果為：[0049][0050] 重復(fù)前述步驟，當(dāng)優(yōu)化函數(shù)的值w達(dá)到最小值時(shí)，得到參數(shù)g。[0051] 其中，TBM模型訓(xùn)練的優(yōu)化函數(shù)為[0052] 其中 是預(yù)測料位，y是真實(shí)料位。優(yōu)化函數(shù)取最小值時(shí)，得到預(yù)測料位 代入前s v

述公式，計(jì)算得到g，再次重復(fù)步驟，得到預(yù)測料位 與振動信號X及振聲信號X的特征參數(shù)

的TBM模型，并通過采集的特征參數(shù)進(jìn)行料位預(yù)測。

[0053] 在得到預(yù)測料位后，采用均方根誤差(RMSE)作為模型性能評價(jià)指標(biāo)，判斷預(yù)測料位結(jié)果的準(zhǔn)確性；其中，均方根誤差(RMSE)的計(jì)算公式為：

[0054][0055] 式中， 和yi分別表示第i個(gè)樣本的料位估計(jì)值和實(shí)際值，N為在線測試樣本個(gè)數(shù)。[0056] 區(qū)別于現(xiàn)有技術(shù)，本發(fā)明的球磨機(jī)料位預(yù)測方法通過采集球磨機(jī)筒壁或者軸承上的振動信號及安裝在球磨機(jī)周圍麥克風(fēng)的振聲信號；振動信號、振聲信號經(jīng)調(diào)理電路、模數(shù)

轉(zhuǎn)換電路輸入數(shù)據(jù)采集卡；通過MFCC分別計(jì)算數(shù)據(jù)采集卡中振動、振聲信號的特征參數(shù)；重

新組合振動、振聲信號；數(shù)據(jù)采集卡收集重新組合后的信號；最后數(shù)據(jù)采集卡中的數(shù)據(jù)通過

歷史數(shù)據(jù)建立的TBM模型，預(yù)測相應(yīng)信號的料位。本發(fā)明能有效的預(yù)測球磨機(jī)的料位，其預(yù)

測料位精確度高、波動小，具有良好的實(shí)用價(jià)值和推廣前景。

[0057] 參閱圖2，圖2是本發(fā)明提供的一種球磨機(jī)料位預(yù)測裝置的結(jié)構(gòu)示意圖。該裝置100連接于待預(yù)測的球磨機(jī)101，包括：

[0058] 采集模塊110，用于采集球磨機(jī)101在不同料位的振動信號Xs及振聲信號Xv；[0059] 第一計(jì)算模塊120，用于根據(jù)采集到的振動信號Xs及振聲信號Xv計(jì)算振動信號Xs及v

振聲信號X的特征參數(shù)；

[0060] 第二計(jì)算模塊130，用于將振動信號Xs及振聲信號Xv的特征參數(shù)作為TBM模型的輸s

入，訓(xùn)練TBM模型，通過多次訓(xùn)練，在使TBM模型的優(yōu)化函數(shù)取最小值時(shí)，建立振動信號X 及

v

振聲信號X的特征參數(shù)與預(yù)測料位的TBM模型，并通過特征參數(shù)計(jì)算預(yù)測料位。

[0061] 其中，第一計(jì)算模塊120包括：[0062] 預(yù)處理單元121，用于振動信號Xs及振聲信號Xv進(jìn)行預(yù)處理；其中，預(yù)處理的方式包括預(yù)加重、分幀和加窗處理；

[0063] FFT變換單元122，用于對預(yù)處理后的振動信號Xs及振聲信號Xv進(jìn)行快速傅里葉變換，將時(shí)域信號轉(zhuǎn)換為頻域信號，并計(jì)算頻域信號的模的平方，得到能量譜；

[0064] Mel濾波單元123，用于對FFT變換單元處理后的能量譜進(jìn)行Mel濾波，并取對數(shù)；[0065] DCT變換單元124，用于對Mel濾波處理后的結(jié)果進(jìn)行離散余弦變換，得到梅爾頻率s v

倒譜系數(shù)作為振動信號X及振聲信號X的特征參數(shù)。

[0066] 如圖3所示，圖3是本發(fā)明提供的一種球磨機(jī)料位預(yù)測裝置的硬件結(jié)構(gòu)示意圖。如圖所示，通過傳感器和麥克風(fēng)采集球磨機(jī)滾筒的不同料位的振動信號和振聲信號，將振動

信號和振聲信號分別發(fā)送到調(diào)理電路以對振動信號和振聲信號進(jìn)行調(diào)理，調(diào)理完成后，將

振動信號和振聲信號傳輸?shù)侥?shù)轉(zhuǎn)換電路進(jìn)行模數(shù)轉(zhuǎn)換，將采集到的模擬信號轉(zhuǎn)換為數(shù)字

信號。轉(zhuǎn)換完成后，將數(shù)據(jù)發(fā)送到一數(shù)據(jù)采集卡1和數(shù)據(jù)采集卡2，為后續(xù)的MFCC處理做準(zhǔn)

備。進(jìn)一步，數(shù)據(jù)采集卡1和數(shù)據(jù)采集卡2將接收到的數(shù)字信號形式的振動信號和振聲信號

進(jìn)行預(yù)處理，預(yù)處理的方式包括預(yù)加重、分幀和加窗處理。預(yù)處理完成后，對數(shù)字信號形式

的振動信號和振聲信號進(jìn)行FFT變換，將時(shí)域信號轉(zhuǎn)換為頻域信號，并計(jì)算其模的平方得到

能量譜。對所得能量譜進(jìn)行Mel濾波，對Mel濾波結(jié)果取對數(shù)后進(jìn)行DCT變換，將計(jì)算結(jié)果作

s v

為振動信號X及振聲信號X的特征參數(shù)。將特征參數(shù)進(jìn)行參數(shù)重組后，輸入到數(shù)據(jù)采集卡3

和數(shù)據(jù)采集卡4，進(jìn)行TBM模型訓(xùn)練。訓(xùn)練的具體步驟為：在訓(xùn)練集中尋找k個(gè)最鄰近樣本并

得到相應(yīng)的信度，信度計(jì)算公式為：

[0067] Φ(d)＝γexp(?d2)，d＝g(x?xi)2，xi為鄰近值，x為訓(xùn)練集中的樣本，γ為常數(shù)；

[0068] 對計(jì)算得到的k個(gè)信度 進(jìn)行DS合成，合成結(jié)果為：[0069][0070] 根據(jù)DS合成結(jié)果計(jì)算預(yù)測料位，計(jì)算結(jié)果為：[0071][0072] 重復(fù)前述步驟，當(dāng)優(yōu)化函數(shù)的值w達(dá)到最小值時(shí)，得到參數(shù)g，將參數(shù)g代入前述的信度計(jì)算公式，通過計(jì)算，可得到預(yù)測料位

[0073]tr tr

[0074] 最終，得到預(yù)測料位 和信度m 的函數(shù)關(guān)系，而m 又與計(jì)算得到的振動信號和振聲信號的特征參數(shù)成函數(shù)關(guān)系，從而得到預(yù)測料位 與振動信號和振聲信號的特征參數(shù)，

通過采集球磨機(jī)的振動信號和振聲信號，通過MFCC方式計(jì)算特征參數(shù)，即可得到預(yù)測料位，

實(shí)現(xiàn)對球磨機(jī)料位的預(yù)測。

[0075] 以上僅為本發(fā)明的實(shí)施方式，并非因此限制本發(fā)明的專利范圍，凡是利用本發(fā)明說明書及附圖內(nèi)容所作的等效結(jié)構(gòu)或等效流程變換，或直接或間接運(yùn)用在其他相關(guān)的技術(shù)

領(lǐng)域，均同理包括在本發(fā)明的專利保護(hù)范圍內(nèi)。