由于三段氣流干燥過(guò)程的復雜性和不確定性,很難通過(guò)機理分析建立機理模型。用傳統的最小二乘法建立的多元線(xiàn)性回歸模型在自變量之間相關(guān)程度高,觀(guān)測值矩陣出現嚴重的多重共線(xiàn)時(shí),會(huì )使回歸系數的均方誤差偏大,從而影響模型擬合的優(yōu)度[6].本文采用主元分析的方法,對相關(guān)的過(guò)程變量進(jìn)行降維處理,挑選出能最大程度的反映原來(lái)多個(gè)變量信息的過(guò)程變量。為了提高模型的精度,采用改進(jìn)的BP神經(jīng)網(wǎng)絡(luò )對模型進(jìn)行修正。
 
    主元回歸建模主元分析主元分析是輸入數據降維處理的主要方法之一,又稱(chēng)為主成分分析或主元素分析,是多元統計過(guò)程控制(MSPC)方法最重要的數學(xué)工具。主元分析利用變量之間的線(xiàn)性相關(guān)關(guān)系對多維信息進(jìn)行統計壓縮,用少部分互不相關(guān)的主元變量描述多維空間的絕大部分的動(dòng)態(tài)信息,以減輕數據分析的復雜度。
 
    假設數據矩陣[X]pn,p代表測量采樣次數,n代表測量變量個(gè)數。主元分析法的基本思想為(1)將原始數據進(jìn)行標準化處理。
 
    xij=xij-MjSjj(1)式中,xij(i=0,1,,n;j=0,1,,m)為經(jīng)過(guò)自標準化的第i個(gè)樣本的第j個(gè)變量;xij為原始變量;Mj,Sj分別是第j個(gè)變量的算術(shù)平均值和標準偏差。
 
    (2)計算其協(xié)方差矩陣R.R=[rij]nn(2)式中,rij=1Ppk=1xkixkj(i,j=1,2,,n)。
 
    (3)計算R的特征值,特征向量。利用雅可比法求矩陣R的n個(gè)非負的特征值12n0,以及對應的特征向量:C(i)=[c1(i)c2(i)cn(i)]T(4)選擇主元。
 
    由特征向量組成n個(gè)新變量:z1=c1(1)x1 c2(1)x2  cn(1)xnz2=c1(2)x1 c2(2)x2  cn(2)xnMzn=c1(n)x1 c2(n)x2  cn(n)xn當前面m個(gè)變量z1,z2,,zm(m
 
    建立主元回歸模型回歸分析是處理變量間相關(guān)關(guān)系的有力工具,它不僅告訴人們怎樣建立變量間的數學(xué)表達式,而且還利用概率統計知識進(jìn)行分析討論,判斷出所建立的經(jīng)驗公式的有效性,從而可以進(jìn)行預測和估計[7].因此,在生產(chǎn)實(shí)際中得到廣泛的應用。
 
    根據現場(chǎng)工藝調查和對機理的定性分析,并考慮到變量的類(lèi)型、數目和測點(diǎn)位置,得到影響氣流干燥過(guò)程精礦含水率的因素有以下11個(gè):精礦量、濕礦含水率、煙氣量、煙氣溫度、燃油量、鼓風(fēng)量(是燃燒風(fēng)、稀釋風(fēng)和氮氣的總和)、熱風(fēng)溫度、機內負壓、混氣室出口溫度、回轉窯尾溫度及沉塵室溫度。對從現場(chǎng)采集的歷史數據進(jìn)行挑選并對它們進(jìn)行統計分析、數據濾波、3法則、目標范圍標準化和主元分析的處理后得到210組數據,其主元貢獻率如所以。
 
    從表中可以看出,前6個(gè)主成分的貢獻率達到了91,也就是說(shuō)前6個(gè)主成分可以描述91的原數據信息,因此我們選取前6個(gè)變量建立回歸模型。取150組數據作為訓練樣本,得到如下的主元回歸模型:y=-0304x1 14722x2 04551x3 03774x4-03445x5 00365x6-08607(3)利用另外60組數據進(jìn)行仿真分析,其預測結果如所示,干精礦含水率的實(shí)際值與預測值的最大相對差為8,預測值能很好地跟隨實(shí)際值在01的范圍內波動(dòng),反映實(shí)際值的變化趨勢,但是很少有實(shí)際值和預測值能夠完全吻合,預測精度有一定的起伏。
 
    主元回歸模型預測模型的校正由于上述模型在預測過(guò)程中部分結果誤差較大,必須對模型進(jìn)行校正。校正通常分為長(cháng)期學(xué)習和短期學(xué)習,如所示,用以克服模型結構算法的復雜性與過(guò)程實(shí)時(shí)性要求之間的矛盾。
 
    軟測量模型校正示意圖本文采用基于神經(jīng)網(wǎng)絡(luò )的自適應殘差補償方法[8]對主元回歸模型給出的預測值進(jìn)行補償。補償算法如下:yi=yi yi式中yi(i=1,2,,n)為干燥后精礦含水率的實(shí)際值,yi為主元回歸模型給出的預測值,yi為實(shí)際值與預測值之差,即殘差。yi=f(x1,x2,x3,x4,x5,x6)根據主元分析的結果,選取n組樣本,采用改進(jìn)的BP神經(jīng)網(wǎng)絡(luò )[9]進(jìn)行訓練,得到殘差的預測模型。利用神經(jīng)網(wǎng)絡(luò )進(jìn)行修正后的智能集成模型的預測曲線(xiàn)如所示。干精礦含水率的實(shí)際值與預測值的最大相對誤差為4,小于主元回歸模型的最大相對誤差,建立的集成模型可以實(shí)現水分的在線(xiàn)檢測。