宮元娟，張建鵬，秦軍偉

(沈陽(yáng)農(nóng)業(yè)大學(xué)工程學(xué)院，遼寧沈陽(yáng)，110161)

胡蘿卜微粉碎后，胡蘿卜粉的孔隙率增大，比表面積也增大，具有良好的溶解性，易于人體消化吸收，可實(shí)現(xiàn)原料的全效利用［1-2］。胡蘿卜富含纖維素、糖分、蛋白質(zhì)及其他生物活性物質(zhì)，使得加工過(guò)程中抗熱性變?nèi)酢ざ仍龃?、韌性亦增強(qiáng)，而使胡蘿卜的粉碎工藝變得復(fù)雜［3-5］。本文通過(guò)采用Matlab數(shù)字圖像技術(shù)計(jì)算胡蘿卜表面積收縮率，研究了其干燥特性，進(jìn)而利用萬(wàn)能試驗(yàn)機(jī)對(duì)胡蘿卜片的物理機(jī)械特性進(jìn)行分析。以不同的切片厚度、切片方式、含水率、粉碎量和粉碎時(shí)間為試驗(yàn)因素，通過(guò)80目標(biāo)準(zhǔn)篩質(zhì)量比為試驗(yàn)指標(biāo)，研究胡蘿卜高速剪切粉碎工藝參數(shù)，通過(guò)單因素試驗(yàn)、正交試驗(yàn)和極差優(yōu)化分析，并用Matlab中stepwise逐步回歸建立數(shù)學(xué)模型與優(yōu)化，確定胡蘿卜采用高速剪切粉碎工藝參數(shù)。

1 材料與方法

1.1 材料

試驗(yàn)所用胡蘿卜品種為紅芯四號(hào)，生胡蘿卜的平均濕基含水率為90%。

1.2 儀器與設(shè)備

DHG-6023A型電熱鼓風(fēng)干燥箱，上海一恒科技有限公司;LD-Y400A型高速萬(wàn)能粉碎機(jī)，上海頂帥電器公司;JY2002型電子天平，上海雙旭電子有限公司;CM7-8102型微電子萬(wàn)能試驗(yàn)機(jī)，美斯特工業(yè)系統(tǒng)有限公司;8411型電動(dòng)振篩機(jī)，上虞市一峰儀器廠。

1.3 實(shí)驗(yàn)方法

1．3．1 胡蘿卜原料前處理

胡蘿卜片的制備工藝流程:

胡蘿卜→初選→清洗→去皮→切片→干燥［6］

胡蘿卜片的長(zhǎng)寬均為20 mm×20 mm，按照不同的切片方式和切片厚度分為6組，前3組:切片方式為軸向，厚度分別為為3、5、7 mm;第4～6組:切片方式為徑向，厚度分別為3、5、7 mm。

1．3．2 胡蘿卜片的干燥特性

以不同的切片方式和切片厚度為試驗(yàn)因素，以平均切片表面積收縮率最小為評(píng)價(jià)指標(biāo)。利用胡蘿卜干燥前后面積變化，用以計(jì)算胡蘿卜片干燥前后面積收縮率，從而分析其干燥特性。

1．3．3 胡蘿卜片的物理機(jī)械特性

采用CM7-8102型微電子萬(wàn)能試驗(yàn)機(jī)對(duì)胡蘿卜片進(jìn)行剪切和壓縮試驗(yàn)，微電子萬(wàn)能試驗(yàn)機(jī)對(duì)胡蘿卜片的預(yù)加載荷力均為0．1N，刀頭下降速度為5 mm/min，以不同的切片方式、切片厚度為因素，以切片平均剪切載荷力最小為評(píng)價(jià)指標(biāo)，分析胡蘿卜干片的物理機(jī)械特性［7-8］。

2 結(jié)果與分析

2.1 胡蘿卜干燥特性分析

測(cè)定胡蘿卜干燥過(guò)程中實(shí)際面積收縮率較復(fù)雜，本文采用與其相關(guān)的表面積收縮率作為指標(biāo)，分析其干燥特性。面積收縮率是干燥前物料面積與干燥后物料面積差與干燥前物料面積之比［9-12］。

式中:S1，干燥前物料的面積;S2，干燥后物料的面積。

已知鮮胡蘿卜片投影面積為4 mm2，只需計(jì)算干燥后胡蘿卜片的面積。利用Matlab圖像處理功能，可以分析圖像中各個(gè)不同特征(如顏色等)點(diǎn)的像素的個(gè)數(shù)。每張圖片都由一定點(diǎn)組成，而圖片中點(diǎn)的個(gè)數(shù)就像素。胡蘿卜干片面積的測(cè)量過(guò)程如圖1所示。

圖1 胡蘿卜面積測(cè)量圖像Fig．1 Image of carrot area measurement

將胡蘿卜干片以白色為背景拍攝照片，如圖1-(a)所示;背景圖如圖1-(b)所示;利用差分法去除圖像的背景得到灰度圖如圖1-(c)所示;利用自動(dòng)閾值法對(duì)圖像進(jìn)行二值化處理得到圖1-(d)。利用Matlab程序計(jì)算出干燥后胡蘿卜片的面積，再根據(jù)面積收縮率的計(jì)算公式從而求得干燥前后的面積收縮率。計(jì)算結(jié)果如表1所示。

表1 胡蘿卜片的面積收縮率Table 1 The area shrinking rate of carrot

由表1所知，當(dāng)胡蘿卜片的切片方式軸向、切片厚度為7 mm時(shí)，表面積收縮率最小為30．2%，其干燥特性最好。

2.2 胡蘿卜物理機(jī)械特性分析

以不同的切片方式、切片厚度為因素，胡蘿卜干片在微機(jī)電子萬(wàn)能試驗(yàn)機(jī)上對(duì)進(jìn)行剪切和壓縮試驗(yàn)。剪切力與位移的變化如圖2所示，壓縮力與位移的變化如圖3所示。

從圖2中可以看出，當(dāng)胡蘿卜片為干燥徑向3 mm時(shí)，剪切力為1．17N。其中徑向片的平均剪切力比軸向片的小;從圖3中可以看出，當(dāng)胡蘿卜片為軸向3 mm時(shí)，此時(shí)壓縮力為4．8N。其中徑向片的平均壓縮載荷小于軸向片。表明胡蘿卜的切片方式為徑向、切片厚度為3 mm時(shí)，胡蘿卜片的剪切力和壓縮力最小。

圖2 剪切力與位移Fig．2 Shear stress and displacement

圖3 壓縮力與位移Fig．3 Compressions stress and displacement

2.3 胡蘿卜高速剪切粉碎工藝參數(shù)分析

2．3．1 單因素試驗(yàn)分析

分析各試驗(yàn)因素與指標(biāo)的關(guān)系，將復(fù)雜的多元問(wèn)題轉(zhuǎn)化為一元的問(wèn)題［13-14］。根據(jù)粉體粒度特性，以通過(guò)80目檢驗(yàn)篩的質(zhì)量百分比為評(píng)價(jià)粉碎效率指標(biāo)。為確定高速剪切粉碎工藝參數(shù)，進(jìn)行5項(xiàng)單因素試驗(yàn)。試驗(yàn)的技術(shù)參數(shù)如表2所示。

表2 單因素試驗(yàn)參數(shù)Table 2 Single factor test parameters

切片方式對(duì)粉碎效率的影響見(jiàn)圖4。不同的切片方式對(duì)粉碎效率影響顯著。胡蘿卜高速剪切粉碎時(shí)，切片方式為徑向時(shí)的粉碎效率明顯優(yōu)于軸向。

圖4 不同切片方式和粉碎時(shí)間對(duì)粉碎效率的影響Fig．4 Impact of different cutting method and crushing time on pulverizing efficiency

由圖4還可知，不同的粉碎時(shí)間對(duì)粉碎效率影響顯著。胡蘿卜高速剪切粉碎時(shí)，當(dāng)粉碎時(shí)間由3 min增加到4 min，通過(guò)80目標(biāo)準(zhǔn)篩的質(zhì)量百分比呈先增長(zhǎng)后下降的趨勢(shì)。當(dāng)粉碎時(shí)間在3．5 min左右時(shí)出現(xiàn)波峰，粉碎效果較好。

切片厚度對(duì)粉碎效率的影響如圖5所示。由圖5可知，胡蘿卜高速剪切粉碎時(shí)，隨著切片厚度的增加，通過(guò)80目標(biāo)準(zhǔn)篩的質(zhì)量百分比在迅速下降。當(dāng)切片厚度為3．5 mm左右時(shí)，粉碎效果較好。

圖5 切片厚度對(duì)粉碎效率的影響Fig．5 Impact of thickness of slice on pulverizing efficiency

含水率對(duì)粉碎效率的影響見(jiàn)圖6。由圖6可知，胡蘿卜高速剪切粉碎時(shí)，含水率從3%增加的7%，通過(guò)80目標(biāo)準(zhǔn)篩的質(zhì)量百分比呈先增長(zhǎng)后下降的趨勢(shì)。曲線的波峰出現(xiàn)在含水率為5%左右，粉碎效率較好。

粉碎量對(duì)粉碎效率的影響如圖7所示，由圖7可知，胡蘿卜高速剪切粉碎時(shí)，通過(guò)80目標(biāo)準(zhǔn)篩的質(zhì)量百分比呈先增長(zhǎng)后緩慢下降的趨勢(shì)。當(dāng)粉碎量為35g左右時(shí)，粉碎效果較好。

2．3．2 胡蘿卜高速剪切粉碎正交試驗(yàn)方案

分析不同的切片厚度、切片方式、含水率、粉碎量和粉碎時(shí)間因素對(duì)通過(guò)80目標(biāo)準(zhǔn)檢驗(yàn)篩的粉體質(zhì)量百分?jǐn)?shù)的影響［15-19］。試驗(yàn)因素與水平如表3所示，正交試驗(yàn)方案及結(jié)果如表4所示。

圖6 含水率對(duì)粉碎效率的影響Fig．6 Impact of moisture content on pulverizing efficiency

圖7 粉碎量對(duì)粉碎效率的影響Fig．7 Impact of crushing amount on pulverizing efficiency

表3 試驗(yàn)因素與水平Table 3 Test factors and levels

表4 正交試驗(yàn)方案及結(jié)果Table 4 Orthogonal experiment methods and results

2．3．3 試驗(yàn)結(jié)果的極差優(yōu)化分析

試驗(yàn)結(jié)果分析如表5所示，由R的大小可以判斷各因素的主次順序。R的值越大，表明該項(xiàng)因素的水平變化對(duì)試驗(yàn)的指標(biāo)值影響就大，從而確定影響因的素最優(yōu)解［20-22］。研究結(jié)果表明，優(yōu)選方案為A2B1C2D3E2。

表5 極差分析Table 5 Range analysis

2．3．4 試驗(yàn)的Matlab回歸模型分析

Matlab的stepwise的逐步回歸分析可知，Intercept=β^0=-31．008 2，β^1=13．940 7，β^4=25．205 9，因B、C、E的t檢驗(yàn)不顯著，經(jīng)逐步回歸分析，剔除B、C、E 3個(gè)非顯著性影響因素［25］，剩下的2個(gè)因素回歸獲得的模型方程為:y^=-31．008 2+13．940 7x1+25．205 9x4，影響因素的大小順序?yàn)?x4＞x1。

模型的殘差分析及檢驗(yàn):由 Matlab中的 rcoplot命令可得出殘差分析如圖8所示。

圖8 殘差分析表Fig．8 Residual analysis table

運(yùn)算結(jié)果可得:r2=0．800，F(xiàn)=183．868 6，P=0．000 7．因 p=0．007 ＜ 0．05，可知回歸模型 y^=-31．008 2+13．940 7x1+25．205 9x4成立，顯著性較好。由Matlab中rstool命令找出極值，并結(jié)合實(shí)際情況得出優(yōu)化方案為:A2B1C2D3E2。

從殘差檢驗(yàn)分析圖8可以看出，除第12個(gè)數(shù)據(jù)外，剩下的17個(gè)數(shù)據(jù)的殘差離零點(diǎn)均較近，且零點(diǎn)包含在殘差的置信區(qū)間內(nèi)，表明回歸模型能良好的符合模擬原始數(shù)據(jù)，而第12個(gè)數(shù)據(jù)可視為異常點(diǎn)，應(yīng)剔除第12組數(shù)據(jù)，剔除后，運(yùn)行源程序得到修正的殘差分析如圖9所示，從而達(dá)到很好的回歸效果，表明回歸模型的顯著性較好。

圖9 修正的殘差分析表Fig．9 Correction of residual analysis table

3 結(jié)論

利用Matlab圖像處理技術(shù)分析胡蘿卜干燥前后面積收縮率，以面積收縮率為評(píng)價(jià)指標(biāo)研究其干燥特性。當(dāng)胡蘿卜片的切片方式軸向、切片為7 mm時(shí)，表面積收縮率最小為62．8%，其干燥特性最好。

采用微電子萬(wàn)能試驗(yàn)機(jī)對(duì)胡蘿卜片進(jìn)行剪切和壓縮試驗(yàn)，分析干燥后的胡蘿卜片的剪切和壓碎機(jī)械特性。表明胡蘿卜的切片方式為徑向、切片厚度為3 mm時(shí)，胡蘿卜片的剪切力和壓縮力最小。

胡蘿卜高速剪切粉碎的最佳工藝參數(shù):對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行極差優(yōu)化分析和Matlab的stepwise的逐步回歸分析，確定優(yōu)化方案為A2B1C2D3E2，即切片方式為徑向、切片厚度為3 mm、含水率為5%、粉碎量為40 g、粉碎時(shí)間為3．5 min時(shí)，胡蘿卜切片干燥粉碎后通過(guò)標(biāo)準(zhǔn)80目篩的質(zhì)量比最高。

［1］鄭瑤瑤，夏延斌．胡蘿卜營(yíng)養(yǎng)保健功能及其開(kāi)發(fā)前景［J］．包裝與食品機(jī)械，2006，24(5):35-37．

［2］張永明．胡蘿卜超微粉制備技術(shù)的研究［J］．農(nóng)產(chǎn)品加工，2005，7(4):12-15．

［3］洪興華，周秀芳．我國(guó)胡蘿卜深加工技術(shù)取得突破性進(jìn)展［J］．中國(guó)科技成果，2002(2):18-19．

［4］向智男，寧正祥．超微粉碎技術(shù)及其在食品工業(yè)中的應(yīng)用［J］．食品研究與開(kāi)發(fā)，2006(2):88-90．

［5］傅蕙英．胡蘿卜應(yīng)用價(jià)值與開(kāi)發(fā)前景［J］．農(nóng)業(yè)產(chǎn)業(yè)化，2006(2):39-41．

［6］于海杰．超微粉碎技術(shù)生產(chǎn)紅棗、胡蘿卜復(fù)合保健營(yíng)養(yǎng)果茶的研究［J］．科技信息，2009(5):15-17．

［7］嚴(yán)怡紅．胡蘿卜食品的加工開(kāi)發(fā)［J］．中國(guó)食物與營(yíng)養(yǎng)，2004(12):43-44．

［8］吳茂玉，葛邦國(guó)．苦瓜超微膳食纖維粉的研究［J］．食品科技，2007，32(3):72-75．

［9］張憋，王亮．超微粉碎在食品加工中的研究進(jìn)展［J］．無(wú)錫輕工大學(xué)學(xué)報(bào)，2003，22(4):106-110．

［10］黃建蓉，李琳，李冰．超微粉碎對(duì)食品物料的影響［J］．糧食與飼料工業(yè)，2007(7):25-27．

［11］萬(wàn)燕君，吳紅霞，張?jiān)Ｖ校芸厍懈畛?xì)粉碎機(jī)理及相關(guān)參數(shù)的分析［J］．食品工業(yè)，2007(1):51-52．

［12］宮元娟，曾程．胡蘿卜微粉物理特性和營(yíng)養(yǎng)成分的影響因素［J］．農(nóng)業(yè)機(jī)械學(xué)報(bào)，2009，40(11):124-128．

［13］廖小軍，胡小松，劉裁．胡蘿卜粉的理化性質(zhì)及其應(yīng)用研究［J］．食品科學(xué)，2004(2):61-64．

［14］宮元娟．胡蘿卜微粉碎工藝及其相關(guān)參數(shù)試驗(yàn)研究［D］．沈陽(yáng):沈陽(yáng)農(nóng)業(yè)大學(xué)，2008:28-36．

［15］宮元娟，曾日新，田素博，等．胡蘿卜精細(xì)加工技術(shù)及其綜合應(yīng)用［J］．農(nóng)業(yè)工程學(xué)報(bào)，2006，22(4):199-203．

［16］宮元娟，于永偉，秦軍偉，等．以圖像計(jì)算收縮率優(yōu)化香菇真空緩蘇干燥工藝［J］．農(nóng)業(yè)工程學(xué)報(bào)，2010，26(12):352-358．

［17］袁媛，李淼，梁青，等．基于水平集的作物病葉圖像分割方法［J］．農(nóng)業(yè)工程學(xué)報(bào)，2011，27(2):208-212．

［18］宮元娟，孟繁壯，秦軍偉，等．北蟲(chóng)草干燥特性與粉碎工藝試驗(yàn)［J］．農(nóng)業(yè)機(jī)械學(xué)報(bào)，2012，43(5):142-146．

［19］李明，張長(zhǎng)利，房俊龍．基于圖像處理技術(shù)的小麥葉面積指數(shù)的提?。跩］．農(nóng)業(yè)工程學(xué)報(bào)，2010，26(1):205-209．

［20］馬莉．MATLAB數(shù)學(xué)實(shí)驗(yàn)與建模［M］．北京:清華大學(xué)出版社，2010:129-160．

［21］張強(qiáng)，王正林．精通MATLAB圖像處理［M］．北京:電子工業(yè)出版社，2012:97-115．

［22］涂望明，魏友國(guó)，施少敏．MATLAB在數(shù)字圖像處理中的應(yīng)用［J］．微計(jì)算機(jī)信息，2007，23(2):21-23．

［23］周元嬌．篩選逐步回歸方法的改進(jìn)研究［D］．揚(yáng)州，揚(yáng)州大學(xué)，2011:15-16．

［24］LIU Zhan-yu，HUANG Jing-feng，SHI Jing-jing，et al．Characterizing and estimating rice brown spot disease severity using stepwise regression，principal component regression and partial least-square regression［J］．Journal of Zhejiang University Science B，2007，8(10):738-744．

［25］JIN Long，HUANG Xiao-yan，SHI Xu-ming．A study on influence of predictor multicollinearity on performance of the stepwise regression prediction equation［J］．Acta Meteorologica Sinica，2010，24(5):593-601．