馬衛(wèi)東，于亞萍※，張櫻飄，李留安，趙 喆

（1.天津農(nóng)學(xué)院工程技術(shù)學(xué)院，天津 300384；2.天津農(nóng)學(xué)院動(dòng)物醫(yī)學(xué)與動(dòng)物科學(xué)學(xué)院，天津 300384；3.金洼農(nóng)業(yè)科技發(fā)展有限公司，天津 300273）

0 引言

雞蛋是受大眾歡迎的食品之一，其具有人體需要的豐富營(yíng)養(yǎng)物質(zhì)。對(duì)雞蛋等級(jí)的評(píng)價(jià)一般有大小、新鮮度、完整和清潔度。雞蛋蛋殼具有保護(hù)蛋白、蛋黃的作用，但蛋殼容易在遭受碰撞后產(chǎn)生破損。蛋殼即使破損微小的裂縫，蛋殼上或周?chē)奈⑸镆矔?huì)迅速污染蛋殼內(nèi)，加速雞蛋內(nèi)部的腐壞，會(huì)嚴(yán)重影響消費(fèi)者的個(gè)人健康。雞蛋是否有裂縫，是判斷雞蛋的重要品質(zhì)指標(biāo)之一，因此，在雞蛋售賣(mài)前，進(jìn)行無(wú)損檢測(cè)，剔除有裂縫的雞蛋是非常重要的一步[1-2]。

國(guó)內(nèi)外對(duì)雞蛋破損檢測(cè)研究主要集中在3種方法。

（1）機(jī)器視覺(jué)技術(shù)。通過(guò)圖像采集蛋殼數(shù)據(jù)，對(duì)圖像進(jìn)行處理，然后通過(guò)建立數(shù)學(xué)模型進(jìn)行分析判斷[3]。但是，圖像檢測(cè)易受禽蛋蛋殼顏色影響，且不易快速、動(dòng)態(tài)進(jìn)行檢測(cè)。

（2）動(dòng)力學(xué)特性技術(shù)。利用動(dòng)力學(xué)特性來(lái)檢測(cè)雞蛋是否破損是根據(jù)敲擊雞蛋表面使其振動(dòng)，通過(guò)采集雞蛋的振動(dòng)信號(hào)來(lái)判別雞蛋是否破損[4-5]。

（3）聲學(xué)特性技術(shù)。蛋殼破損后，其內(nèi)部結(jié)構(gòu)便隨之遭到破壞。破損后蛋殼的聲學(xué)特性與破損前相比有十分明顯的差別[6-9]。敲擊完好的雞蛋，發(fā)出的聲音比較清脆；敲擊有裂紋的雞蛋，發(fā)出的聲音比較沉悶。

本設(shè)計(jì)主要是應(yīng)用聲學(xué)特性規(guī)律，設(shè)計(jì)了一套判斷雞蛋是否有裂縫的測(cè)試系統(tǒng)，通過(guò)敲擊雞蛋，將發(fā)出的聲音進(jìn)行采集，通過(guò)Matlab 對(duì)采集的數(shù)據(jù)進(jìn)行分析，實(shí)現(xiàn)對(duì)裂縫蛋和完好蛋的判別。

1 系統(tǒng)組成

整個(gè)系統(tǒng)主要由敲擊裝置、聲音采集裝置、數(shù)據(jù)采集卡、LabVIEW 人機(jī)接口界面和Matlab 軟件分析，系統(tǒng)組成如圖1所示。

圖1 系統(tǒng)組成框圖

敲擊裝置主要完成對(duì)每個(gè)雞蛋的敲擊，敲擊頭和敲擊棒是一體的，設(shè)計(jì)后通過(guò)3D打印機(jī)采用PLA（聚乳酸）材料將敲擊棒打印出來(lái)，該材料重量輕，敲擊發(fā)聲清晰，且不會(huì)對(duì)雞蛋造成損害，敲擊的角度固定30°。敲擊的位置在蛋小頭位置，蛋托上放置雞蛋，固定敲擊棒與雞蛋的距離。蛋托的另一端固定微音器及放大電路，將敲擊的聲音放大，然后通過(guò)數(shù)據(jù)采集傳送到上位機(jī)，上位機(jī)的LabVIEW軟件編寫(xiě)人機(jī)交互界面，采集數(shù)據(jù)和波形[10]，將采集到的數(shù)據(jù)通過(guò)Excel 存儲(chǔ)，截取有用數(shù)據(jù)后，利用Matlab軟件進(jìn)行多變量數(shù)據(jù)分析。

圖2 測(cè)試系統(tǒng)實(shí)物圖

2 數(shù)據(jù)采集及數(shù)據(jù)分析

2.1 數(shù)據(jù)采集程序

數(shù)據(jù)采集卡選用NI USB-6211，12 位精度，輸入電壓范圍為-10～10V，結(jié)合NI-DAQmx 設(shè)備驅(qū)動(dòng)程序完成數(shù)據(jù)測(cè)量。首先創(chuàng)建一個(gè)模擬電壓輸入的虛擬通道，采用DAQmx Create Virtual Channel 模塊，接線端的配置設(shè)定為“單端接地”。通過(guò)定時(shí)模塊DAQmx Timing.vi對(duì)信號(hào)的采樣率、采樣數(shù)和采樣模式進(jìn)行設(shè)定。系統(tǒng)采樣率為2000 S/s，采樣模式為“有限采樣”；每個(gè)通道采樣數(shù)設(shè)定為1024。DAQmx Start Task.vi 來(lái)啟動(dòng)模擬輸入任務(wù)。然后從指定的虛擬通道和任務(wù)中讀取采樣值通過(guò)DAQmx Read.vi 來(lái)實(shí)現(xiàn)，測(cè)量完成后，通過(guò)DAQmx Stop Task.vi 模塊來(lái)停止測(cè)試，并由DAQmx Clear Task.vi 來(lái)消除任務(wù)，釋放任務(wù)占有的資源[11]。LabVIEW 上位機(jī)將測(cè)試圖形數(shù)據(jù)存儲(chǔ)到EXCEL文件中，為后續(xù)Matlab進(jìn)行多變量數(shù)據(jù)分析提供數(shù)據(jù)。數(shù)據(jù)采集程序如圖3所示。

圖3 數(shù)據(jù)采集程序框圖

2.2 數(shù)據(jù)測(cè)試

對(duì)60個(gè)雞蛋進(jìn)行測(cè)試，其中完好雞蛋30個(gè)，破損雞蛋30個(gè)。每一個(gè)樣品采集1024個(gè)數(shù)據(jù)。形成60 （樣品數(shù)）×1024（數(shù)據(jù)點(diǎn)）的矩陣。為了保證數(shù)據(jù)分析的準(zhǔn)確性，根據(jù)測(cè)試圖形的結(jié)果，去除部分無(wú)用數(shù)據(jù)。因?yàn)檠b置敲擊雞蛋時(shí)聲音幅度變化大，因此通過(guò)先求出聲音幅值的最大值，然后截取最大值的前50 個(gè)數(shù)據(jù)和后面的150 個(gè)數(shù)據(jù)，共得出60（樣品數(shù)）×201（數(shù)據(jù)點(diǎn)）進(jìn)行分析。

圖4 數(shù)據(jù)測(cè)試平臺(tái)

2.3 數(shù)據(jù)分析

圖5 主成分大小

對(duì)截取后的有用數(shù)據(jù)矩陣進(jìn)行分析，其主成分大小如圖5所示，前3個(gè)主成分分別占22.1%、20.8%和9.2%。因此，根據(jù)主成分1、主成分2、主成分3 畫(huà)出完好雞蛋與裂縫蛋的得分圖如圖6所示。圖中，五角星為裂縫蛋，星花“”為完好雞蛋，還可以看出，完好雞蛋主要集中在PC1＞-1 部分，而裂縫雞蛋集中在PC1＜-1 部分；在PC2 上，完好雞蛋主要集中在[-2，-4]之間；在PC3上，區(qū)分度不大。有3個(gè)裂縫雞蛋和5個(gè)完好雞蛋在得分圖上混合在一起，無(wú)法區(qū)分。

圖6 完好雞蛋和裂縫雞蛋得分圖

3 結(jié)束語(yǔ)

雖然有裂縫的雞蛋和完好雞蛋在雞蛋敲擊產(chǎn)生聲音上有所區(qū)別，根據(jù)縫隙的大小和縫隙所在位置的不同等因素，敲擊獲取聲音的結(jié)果不同。根據(jù)圖6所示的得分圖，裂縫蛋和完好雞蛋之間通過(guò)聲音敲擊有一定的識(shí)別度，通過(guò)單次敲擊獲取的信息識(shí)別度偏低。采用一個(gè)聲音信號(hào)進(jìn)行數(shù)據(jù)采集，系統(tǒng)辨識(shí)度低，接下來(lái)準(zhǔn)備對(duì)蛋中、雞蛋大頭、雞蛋小頭進(jìn)行敲擊檢測(cè)，采用信息融合技術(shù)進(jìn)行識(shí)別，提高其識(shí)別度。