最小二乘法耐抗線回歸方程在制絲加工預(yù)測環(huán)節(jié)的融合分析

高淼 張曉峰 張春廣 楊光露 王迪

摘要：隨著我國技術(shù)水平的發(fā)展，人們對卷煙產(chǎn)品的質(zhì)量要求也越來越高，為了提高煙卷質(zhì)量和耐加工性，對于煙葉制絲過程各道工序質(zhì)量及數(shù)據(jù)進行了分析。因計算結(jié)果在進行回歸預(yù)測時容易受異常值影響，為避免各工序在加工過程中各項數(shù)據(jù)出現(xiàn)偏差，故設(shè)置相關(guān)實驗引入耐抗線方法。實驗結(jié)果顯示，煙葉在制絲加工過程中引用耐抗線方法，不僅提高了數(shù)據(jù)的準確性，還提高了預(yù)測結(jié)果的耐抗性，充分說明耐抗線方法在制絲加工領(lǐng)域的預(yù)測良好效果和應(yīng)用價值。

關(guān)鍵詞：耐抗線方法;制絲加工預(yù)測;最小二乘法回歸方程

中圖分類號：TS452

文獻標識碼：A

文章編號：1001-5922（2020）12-0118-03

0 引言

在制絲加工領(lǐng)域中，煙葉制絲與傳統(tǒng)機械制造領(lǐng)域的精加工不同，因設(shè)備情況與外部環(huán)境都會對煙葉或煙絲的物理特性產(chǎn)生影響，所以為避免異常數(shù)據(jù)值造成預(yù)測結(jié)果偏離情況，需利用最小二乘法計算線性回歸方程。全文為預(yù)測切后梗絲含水率，通過耐抗線回歸方程來探索煙梗水槽式回潮機和梗絲含水率兩者之間的相關(guān)聯(lián)性，并消除計算過程中異常值的影響。耐抗線法是通過繪制圖表、方程擬合和計劃某特征量等對數(shù)據(jù)的整體結(jié)構(gòu)和規(guī)律進行探索分析的一種方法，主要適用于線性歸擬合，其作用是消除異常值對回歸方程產(chǎn)生的影響，提高回歸方程的耐抗性，在舍棄異常值的同時迭代修正數(shù)據(jù)[1]。

1 實驗部分

1.1 實驗材料

將環(huán)境濕度控制在40% -60%，溫度控制在20℃-26℃，準備實驗所需材料及儀器進行試驗[2]。整個實驗過程中所需的原料選用配方煙梗，實驗所用儀器需準備密封袋、水槽式煙?；爻睓C、天平秤以及電熱恒溫干燥箱等，其中電子天平型號選擇為CP223C，干燥箱選擇型號為FED240。

1.2 實驗方法

根據(jù)2003版《卷煙工藝規(guī)范》詳細的操作步驟進行煙?；爻痹囼灐Ｔ诒ＷC溫度與濕度穩(wěn)定的環(huán)境中，水溫控制在55℃，實驗儀器及設(shè)備參數(shù)固定不變，對煙?；爻睓C進行分批固定工藝投料（ 1500kg/h），每批設(shè)置為2000kg，待貯存時間達到4h后，在出料口取煙梗取樣21，對其進行含水率檢測[3]。

1.3 耐抗線的計算方法

1.3.1 進行數(shù)據(jù)分組

將數(shù)據(jù)（Ai，Bi）按照數(shù)據(jù)自變量由小到大分為左、中、右3組的樣本量，分別是E、F、Q。對于樣本量無法準確判斷是不是3的倍數(shù)情況下，可對樣本數(shù)據(jù)進行分組處理。表1所示，當樣本量n取值不同時的分組方式[4]。

為了規(guī)避異常值對回歸方程造成的不良影響，利用以上試驗算法，同時也為了得到回歸系數(shù)，將回歸方程最大限度數(shù)據(jù)作為主體部分[6]。

2 結(jié)果與討論

2.1 建立耐抗線方程

將表2中21個樣本數(shù)據(jù)觀測值劃分為每7個一組，計算出各組數(shù)據(jù)的中位數(shù)作為總括點，詳情統(tǒng)計結(jié)過如表3所示。

為準確得到直線擬合方程的初始估計可利用總括點計算，該形式如下：

貯存后煙梗含水率=34.84- 0.0509（浸梗時間一51） =37.232-0.0509浸梗時間

為得到修正后的耐抗線方程，需利用殘差對該方程進行迭代修正：

貯存后煙梗含水率=30.1104+0.1163浸梗時間

由圖1可知當煙梗浸梗時間為73s時，有其中一個數(shù)據(jù)發(fā)生偏離，且偏離程度較大，其他數(shù)據(jù)未見異常偏離，既避免了傳統(tǒng)最小二乘法對異常觀測值的敏感性，又表現(xiàn)了線性數(shù)據(jù)關(guān)系，并且樣本的異常數(shù)據(jù)更具有耐抗性。

2.2 耐抗線法與最小二乘法回歸方程預(yù)估結(jié)果對比

利用傳統(tǒng)最小二乘法計算得出回歸方程：貯存后煙梗含水率=30.3368+0.1019浸梗時間。如圖1所示，虛線為最小二乘法回歸方程線，實線則為耐抗線，都是在同一散點處進行繪制的浸梗含水率，以此來進行對比。

據(jù)對比情況對最小二乘法回歸方程進行分析，該方法不具備耐抗性，主要原因是該方法在進行計算時，并未將納入計算的數(shù)據(jù)進行修正，在浸梗73s時出現(xiàn)的異常偏離在該方法計算中也產(chǎn)生了一定影響。

2.3 耐抗線法與最小二乘法回歸方程預(yù)測效果度量

用Q2，Q2 xdj和殘差標準值P計算對比，為進一步衡量兩個回歸方程的預(yù)測效果進行檢驗[7]。

2.3.1 耐抗線回歸方程檢驗與統(tǒng)計

回歸方程在進行顯著性檢驗過程中，J的觀測值落入了拒絕域，J= 23.939>J0.95（1，19） =4.46，結(jié)論顯示該回歸方程具有顯著性意義;且J2 Xdj≤J2數(shù)值相近，充分說明模型的擬合較好[10]。

2.3.3 對比預(yù)測效果

將耐抗線回歸方式?jīng)Q定系數(shù)的調(diào)整值與最小二乘法回歸方程決定系數(shù)的調(diào)整值相比較，耐抗線回歸方程值大于最小二回歸方程值數(shù)，并且耐抗線回歸方程的殘差標準值比最小二恢復(fù)方程的殘差標準值要小，以上這兩點充分表明了耐抗線回歸方程的預(yù)測效果要比最小二更好一些。

3 結(jié)語

根據(jù)耐抗線法對煙草制絲加工過程進行研究，既規(guī)避了最小二乘法對異常觀測值的敏感性，還使異常觀測值不會對耐抗線法建立的回歸方程產(chǎn)生影響。由于葉片、煙梗等原料等級的不同，潤透性和組分間的含水率也存在著一定的差異，從而導(dǎo)致在煙梗預(yù)處理和葉片處理加工過程中出現(xiàn)工藝要求和物理含水的偏差，為了避免出現(xiàn)數(shù)據(jù)偏差，可將耐抗線法運用至加工中，凸顯耐抗線方法的優(yōu)勢，提高檢驗數(shù)據(jù)的準確性，同時還減少了取樣含水率的異常觀測值。

參考文獻

[1]李達，李曉科，張濤，等，卷煙揮發(fā)性香氣成分測定及煙絲加香工藝優(yōu)化[J].食品與機械，2020，36（03）：84-90.

[2]張馨予，嚴莉紅，易小麗，等.不添加料液的混合型膨脹煙絲生產(chǎn)工藝及質(zhì)量分析[J].中國煙草學(xué)報，2020，26（ 01）.36-41+47.

[3]付亮，易斌，溫亞東，等.低溫真空干燥對煙絲及卷煙品質(zhì)的影響[J]。中國煙草學(xué)報，2019，25（03）：23-28.

[4]郭高飛，劉朝賢，李斌，等，煙粉和煙絲的熱解燃燒特性及反應(yīng)動力學(xué)分析[J].煙草科技，2019，52（06）：90-98.

[5]王超，黃飛，王艷麗，等.UPLC-IE法測定無煙氣煙草制品及卷煙煙絲中的羰基化合物[J].食品工業(yè)，2017，38（03）：256-261.

[6]周冰，唐軍，何邦華，等.卷煙成品煙絲中不同長度煙絲的化學(xué)特性[J].煙草科技，2015，48（09）：40-44.

[7]史近文，衛(wèi)青，付祺，等，造紙法再造煙葉著色劑的制備工藝研究[J].云南農(nóng)業(yè)大學(xué)學(xué)報（自然科學(xué)），2019，34（05）：802-807.

[8]葉建斌，呂麗文，王璐，等，造紙法再造煙葉工藝流程中中性香味物質(zhì)的變化[J].河南農(nóng)業(yè)大學(xué)學(xué)報，2019，53（03）：415-425.

[9]王充，龍章德，孟丹丹，等.造紙法再造煙葉濃縮液中微生物群落結(jié)構(gòu)分析及應(yīng)用[J].煙草科技，2019，52（03）：33-38.

[10]王小升，蘇丹丹，伏穎，等，打漿度對造紙法再造煙葉燃燒性能和煙氣指標的影響[J].河南農(nóng)業(yè)科學(xué)，2019，48（01）：146-152。